Reduced CO2 fluid as an agent of ore-forming processes: a case study of dolomite-replacement skarns at the Yoko-Dovyren massif)

Full Text

Йоко-Довыренский расслоенный ультрабазит-базитовый интрузив приурочен к южной окраине Сибирского кратона. Этот интрузив является одним из элементов Сыннырской вулканно-плутонической структуры, образовавшейся 700–760 млн лет назад (Конников и др., 1999). Возраст интрузива, оцененный U-Pb методом по цирконам (Арискин и др., 2013), составляет 728 ± ± 3 млн лет. В Йоко-Довыренском интрузиве в 1950-ые годы были обнаружены рудопроявления металлов платиновой группы в приподошвенных залежах медно-никелевых руд. На поверхности обнажается полный поперечный срез моноклинально наклоненного интрузива мощностью около 3.5 км от нижней ультраосновной части до прикровельных гранофиров (рис. 1). В последние годы выявлено несколько дополнительных горизонтов малосульфидного оруденения (рифов) вблизи границы ультрамафитов с трактолитами. Несмотря на небольшие подтвержденные запасы руд, Йоко-Довыренская интрузия наиболее полно изучена среди расслоенных ультрабазит-базитовых интрузий России.

 

Рис. 1. Схема геологического строения Йоко-Довыренского расслоенного массива и положения платиноносного горизонта (Риф I). 1 – плагиоклазовые перидотиты; 2 – плагиоклазовые дуниты, верлиты; 3 – ритмическое чередование плагиоклазовых дунитов, троктолитов, оливиновых габбро; 4 – оливиновые габбро и габбронориты; 5 – силлы гранофировых габброноритов; 6 – положение платиноносного горизонта (Риф I); 7 – эффузивы (трахилипариты, трахидациты, андезиты); 8 — кварциты, сланцы, песчаники; 9 – карбонатные породы; 10 – апокарбонатные скарны; 11 – разрывные нарушения. Цифрами отмечены точки опробования: 1 – вмещающих доломитов и кварцита, 2 – диопсид-кальцит-волластонитового скарна.

 

Глубина становления Йоко-Довыренского массива (до наклона с образованием моноклинали на заключительной стадии тектонических деформаций) оценивалась всего лишь в 1.7 км по результатам изучения процессов скарнообразования и геологическим реконструкциям (Р = 500 бар) (Перцев, Шабынин, 1978). Поскольку Йоко-Довыренская интрузия имеет большой размер с поперечным сечением 3.5 × 26 км (мощность × длинная ось) (см. рис. 1), исходя из соображений механической устойчивости, глубина его становления должна быть больше, чем 1.7 км. Алтухов с соавторами (1990) на основании гравиметрических наблюдений установили наличие пород повышенной плотности (ультрабазитов) на глубинах 25, 18, 11 и 5 км под Йоко-Довыренским интрузивом. Эти наблюдения можно связать с образованием серии магматических камер, из которых Йоко-Довыренская – самая близповерхностная. В настоящей работе приводится новая уточненная оценка глубины становления по составу магматического амфибола.

На стадии магматической кристаллизации магмы Йоко-Довыренского интрузива ЭПГ концентрировались в сульфидах (Ariskin et al., 2016). По мнению ряда исследователей (Mathez, 1995; Ballhaus, Stumpfl, 1986; Boudreau, Meurer, 1999; Aird, Boudreau, 2013), в формировании платиновых рифов Бушвельда и Стилуотера, генетически подобных платиноносному рифу Йоко-Довырена, большую роль играют субсолидусные реакции кумулуса и флюида. Предполагается, что в раннюю постмагматическую стадию ЭПГ могут высвобождаться из сульфидов в результате реакции водного флюида с сульфидсодержащим кумулусом (Boudreau, Meurer, 1999) и переотлагаться в рифах. Чисто водный состав флюида противоречит высокому содержанию углеводородов в платиноносном рифе Йоко-Довыренского интрузива (Konnikov et al., 2000). Через весь интрузив трассируется горизонт с апокарбонатными скарнами (рис. 1). На состав флюида оказывало влияние выделение СО₂ в результате протекания реакций скарнообразования на контакте магмы с доломитами (Wenzel et al., 2002), реакции декарбонатизации происходили также и в термическом ореоле интрузива. Для понимания механизмов метасоматического переноса ЭПГ важно оценить валовый состав флюида, летучесть кислорода fO₂, соотношения различных компонентов в нем.

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

МЕТОДЫ АНАЛИЗА И РАСЧЕТОВ

Исследование состава пород и минералов методом ICP-MS было проделано В.К. Карандашевым в аналитическом сертификационно-испытательном центре ИПТМ РАН (г. Черноголовка). Для контроля точности ана- лиза использовался сертифицированный стан- дарт (габбро эссекситовое СГД-2А (ГСО 8670--2005)). Анализ минералов локальным рентгено-спектральным методом осуществлен в ИЭМ РАН на микроскопе TESCAN Vega TS5130MM с энергодисперсионной приставкой INCA Energy 450 (напряжение – 20 кВ, ток поглощения электронов на стандарте Со – –180 нА, размер зонда ~ 0.150 мкм).

Термодинамические расчеты в работе носят иллюстративный характер, демонстрируя возможные пути протекания природных процессов, и не претендуют на большую точность. Расчеты проводились с использованием базы термодинамических данных (Holland, Powell, 1998) и в отдельных оговоренных ниже случаях по (Robie et al., 1978) в среде символьного манипулятора Maple 9.5. Уравнения закона действующих масс решались в простых случаях аналитически, а для многокомпонентных флюидов численно. Термодинамические расчеты равновесий в флюидах проводились в приближении идеального смешения неидеальных компонентов с коэффициентами летучести СО, Н₂, СН₄, COS, H₂S по (Shi, Saxena, 1992). Коэффициенты фугитивности Н₂О и СО₂ вычислялись с учетом неидеальности смешения по (Duan, Zhang, 2006; онлайн вычисление http://models.kl-edi.ac.cn/models/h2o_co2/index.htm), малыми содержаниями других компонентов мы пренебрегли.

РЕЗУЛЬТАТЫ ИЗУЧЕНИЯ МИНЕРАЛОВ И ПОРОД

Нами были исследованы образцы диопсидового скарна из коллекции Э.Г. Коникова, отобранные из апокарбонатных скарнов в троктолит-плагиодунитовой зоне (ручей Белый). Для сравнения изучены карбонатные протолиты из нижнего контакта Йоко-Довыренского интрузива, отобранные во время полевых работ 2015– –2016 гг. Е.В. Кисловым. Точки опробования указаны на рис. 1.

Описание вмещающих пород. Рентгенофазовый анализ образцов, отобранных по профилю у нижнего контакта Йоко-Довыренского интрузива, показал, что доломиты частично разложились с образованием периклаза (позднее гидратированного во вторичный брусит с пластинчатой морфологией) и кальцита (аналогично в (Wenzel et al., 2002)):

CaMg(CO₃)₂ = MgO + CaCO₃ + CO₂. (1)

При давлении 2 кбар и ХСО₂ = 1 разложение происходит при 850^оС. Первичный брусит образуется при малой доле углекислоты (например, при 1 кбар и ХСО₂ < 0.1 (Ganino et al., 2013)). Степень разложения снижается в крайней точке опробования в 10 м от контакта. Изучение образцов карбонатов на микрозонде позволило обнаружить редкие включения лизардита (Mg₃Si₂O₉(OH)₂), ильменита, магнетита, пирротина. Кальцит содержит 2–3 мол.% MgCO₃, что отвечает уравновешиванию этих минералов при низкой температуре. Также исследован прослой кварцита из карбонатной толщи, отобранный в 100 м от контакта. Содержание кремнезема в кварците составляет около 92 мас.%. На микрозонде идентифицированы силикатные минералы хлорит, пренит, мусковит, а также акцессорные циркон, рутил и пирротин (табл. 1). Интерес представляют находки цинковых фаз – цинкового ильменита и силиката цинка. Наличие в парагенезисе пренита отвечает низкотемпературной (регрессивной) зеленосланцевой фазе становления: Т = 200–300^оС, давление до 2–3 кбар (Frey, Robinson, 1999). Обнаружение цинковых фаз говорит о том, что на определенной стадии концентрация цинка в циркулирующем флюиде была высокая.

 

Таблица. 1. Состав минералов Йоко-Довыренского массива

Компоненты	Карбонаты экзоконтакта
	Брусит (5)			Доломит (8)			Кальцит (5)			Серпентин (7)			Ильменит (3)			Магнетит (5)
	C, мас.%		1s	C, мас.%		1s	C, мас.%		1s	C, мас.%		1s	C, мас.%		1s	C, мас.%		1s
SiO2	-		-	-		-	0.	08	0.06	41.	31	1.43	-		-	-		-
TiO2	-		-	-		-	-		-	0.	20	0.16	63.	33	0.78	0.	11	0.11
Al2O3	0.	13	0.11	-		-	-		-	2.	84	1.09	0.	10	0.09	0.	08	0.05
Cr2O3	-		-	-		-	-		-	0.	12	0.14	0.	15	0.15	0.	09	0.08
FeO	0.	45	0.06	0.	03	0.06	0.	27	0.13	0.	67	0.12	16.	14	0.53	86.	06	3.41
MnO	-		-	0.	06	0.09	-		-	-		-	3.	72	0.21	-		-
MgO	59.	78	0.62	19.	57	1.51	1.	02	0.68	40.	88	1.01	20.	35	0.69	1.	78	0.23
CaO	0.	09	0.06	32.	13	1.05	54.	49	0.70	0.	19	0.14	0.	39	0.05	-		-
Na2O	-		-	-		-	-		-	-		-	-		-	-		-
K2O	-		-	-		-	-		-	-		-	-		-	-		-

 

Компоненты	Кварциты экзоконтакта
	Мусковит (3)*			Пренит (5)					Хлорит (5)				Циркон (4)				Zn-ильменит (3)				Zn-силикат (1)
SiO2	44.	91	0.34	38.	16		2.95		27.	57		1.39	30.	93	0.62		0.	47	0.	12	28.	23	-
TiO2	0.	10	0.09	-			-		0.	11		0.09	-		-		51.	76	54.	21	-		-
Al2O3	30.	89	0.12	27.	31		1.15		20.	91		1.09	0.	04	0.04		0.	33	0.	24	5.	68	-
Cr2O3	-		-	0.	04		0.06		0.	27		0.31	0.	05	0.08		0.	49	0.	00	0.	84	-
FeO	1.	05	0.17	-			-		20.	96		1.22	0.	21	0.11		24.	82	16.	52	0.	8	-
MnO	-		-	0.	15		0.10		0.	30		0.19	0.	17	0.06		8.	65	6.	44	0.	05	-
NiO	-		-	-			-		-			-	0.	12	0.10		-		-		-		-
ZnO	-		-	-			-		-			-	-		-		10.	84	20.	81	49.	3	-
MgO	1.	29	0.07	0.	06		0.06		10.	66		0.64	-		-		0.	00	0.	23	2.	31	-
CaO	0.	06	0.04	23.	83		0.95		0.	13		0.11	-		-		0.	00	0.	00	-		-
Na2O	0.	22	0.15	-			-		-			-	-		-		1.	66	2.	93	6.	15	-
K2O	11.	25	0.05	-			-		-			-	-		-		-		-		-		-
ZrO2	-		-	-			-		-			-	64.	12	0.99		-		-		-		-
HfO2	-		-	-			-		-			-	1.	59	0.26		0.	75	0.	74	-		-
Компоненты	Состав минералов из скарна
	Диопсид (10)			Монтичеллит (9)							Ранкинит (5)					Волластонит (1)			Ca5Si4O13 (1)			Апатит (1)			Куспидин (1)
SiO2	53.	04	0.85	36.		76		0.53			40.28		0.60			49.	55	-	46.	08	-	1.	75	-	31.	95	-
TiO2	0.	06	0.05	0.		03		0.05			-		-			-		-	-		-	0.	23	-	-		-
Al2O3	0.	15	0.09	0.		14		0.06			0.07		0.07			-		-	0.	14	-	-		-	0.	14	-
Cr2O3	-		-	-				-			0.17		0.23			-		-	-			-		-	-		-
FeO	0.	05	0.08	0.		10		0.07			-		-			-		-	-		-	0.	14	-	-		-
MnO	0.	06	0.07	-				-			-		-			-		-	-		-	-			-		-
MgO	17.	38	0.45	23.		91		0.39			0.25		0.39			0.	21	-	0.	09	-	0.	2		0.	18	-
CaO	27.	08	0.18	37.		86		0.39			53.95		0.85			50.	60	-	53.	19	-	56.	39	-	63.	7	-
Na2O	0.	10	0.06	0.		09		0.07			-		-			-		-	0.	16	-	0.	02	-	-		-
K2O	0		-	-				-			-		-			-		-	-		-	0.	06	-	-		-
F	0		-	0.		27		0.19			0.38		0.32			-		-	0		-	0.	45	-	10.	0	-
Cl	0		-	-				-			-		-			-		-	-		-	0.	26	-	-		-
P2O5	0		-	-				-			-		-			-		-	-		-	38.	65	-	-		-
S	0		-	-				-			-		-			-		-	-		-	0.	32	-	-		-

*В скобках указано количество усредненных анализов.

 

Голубой диопсид. Апокарбонатные скарны находятся внутри интрузива, максимальные температуры их образования приближались к магматическим. Эти скарны изучаются с 1960-х годов (см. в (Кислов, 1998)). Особенно подробно проанализирован процесс формирования магнезиальных скарнов, состоящих из брусита (псевдоморфозы по периклазу), оливина, шпинели (Wenzel et al., 2002). Эти скарны образовались при максимальном нагреве свыше 1100^оС, вызывающем плавление кальцита. Нами из- учены образцы диопсид-кальцит-волластонитового скарна из блоков бруситовых скарнов в троктолитах. Согласно результатам анализов методом ICP-MS (табл. 2), зона голубого диопсида характеризуется атомарным отношением Mg/Ca = 0.81, а секущая осветленная зона с высоким содержанием кальциевых силикатов – Mg/Ca = 0.72. Обе зоны обогащены кальцием относительно исходного доломита (Mg/Ca = 1.02 ± 0.02). Светлая зона носит признаки выноса кремнезема. На это указывает парагенезис монтичеллита (Mg_0.94Ca_1.06SiO₄) и переотложенного диопсида. Кристаллы голубого диопсида имеют изометричный габитус и размер зерен до 1 мм (рис. 2а). Внутри зерен диопсида встречаются включения безводного кальциевого силиката состава Ca₅Si₄O₁₃, также отмечены кристаллы волластонита. Кристаллы переотложенного диопсида имеют удлиненный габитус с поперечником до 2–5 мм и длиной до 150 мкм (рис. 2б). О высокой активности фтора и низкой воды во флюиде при формировании осветленной зоны свидетельствует находка безводного куспидина Ca₄(Si₂O₇)F₂. Согласно микрозондовым анализам, содержание фтора в монтичеллите (видимо, в форме микрофлюидных включений) также высокое и составляет 0.2–0.6 мас.%. Более поздний привнос кальция отражается в наложенной кристаллизации кальциевого силиката ранкинита (Ca₃Si₂O₇), который образует крупные кристаллы в осветленной зоне и мелкозернистую массу, заполняющую трещины между кристаллами голубого диопсида. Пониженные суммы оксидов по результатам микрозондового анализа ранкинита, видимо, свидетельствуют о его вторичной частичной гидратации с образованием Ca₃(SiO₃OH)₂·2H₂O (авфиллита). Кристаллы апатита в ранкините имеют повышенное содержание фтора и низкое гидроксила – Ca₅(PO₄)₃(OH_0.05-0.11,F_0.87-0.93,Cl_0.02). Апатит из монтичеллит-диопсидовой зоны содержит карбонат-ион, а также больше гидроксила, хлора и серы Ca₅(PO₄)_2.69(CO₃)_0.29(SO₄)_0.02(OH_1.15,F_0.12,Cl_0.04). Состав апатита отражает увеличение концентрации фтора в позднем флюиде, что отмечалось ранее (Кислов, 1998).

 

Таблица 2. Состав скарнов и контактовых пород

Компоненты	Карбонат 1		Карбонат 2		Кварцит		Голубой диопсид		Голубой диопсид*		Осветленная зона в голубом диопсиде		Осветленная зона в голубом диопсиде**
SiO2	-		-		91.	5	55.	20	54.	84	-		-
TiO2	0.	039	0.	026	0.	35	0.	1	0.	027	0.	0075	0.	0072
Al2O3	0.	48	0.	36	3.	6	0.	050	0.	079	0.	11	0.	059
Fe2O3	0.	26	0.	11	1.	7	1.	075	0.	11	0.	13	0.	12
MnO	0.	026	0.	014	0.	042	0.	0039	0.	0051	0.	011	0.	011
MgO	26.	6	25.	2	0.	71	16.	3	16.	4	18.	1	17.	9
CaO	35.	6	35.	3	1.	6	28.	1	28.	4	34.	7	34.	1
Na2O	0.	014	0.	016	0.	039	0.	053	0.	039	0.	014	0.	012
K2O	0.	0080	0.	0059	0.	12	0.014		0.	015	0.	011	0.	0067
P2O5	0.	15	0.	27	0.	015	0.041		0.	018	0.	051	0.	047
Sобщ	0.	004	< ПО		0.	0034	0.021		0.	0064	< ПО		< ПО
Н2О	7.	57	6.	27	-		-		-		-		-
Li	1.	6	1.	0	19.	8	0.33		0.	39	3.	7	3.	0
Be	1.	0	0.	25	0.	71	0.13		0.	19	2.	0	1.	6
Sc	0.	65	0.	44	3.	2	0.27		0.	32	0.	15	0.	14
V	5.	3	9.	2	19.	4	300		345		7.	3	6.	1
Cr	13.	8	32.	8	261		3.	8	80.	2	98.	6	88.	7
Co	1.	0	0.	92	2.	9	0.	79	1.	0	1.	2	1.	2
Ni	9.	1	23.	4	107		6.	3	55.	4	68.	8	59.	9
Cu	8.	6	1.	3	11.	7	3.	9	5.	9	7.	6	7.	2
Zn	26.	7	7.	2	34.	2	11.	0	6.	6	5.	8	4.	7
Ga	0.	43	0.	38	5.	6	0.	15	0.	10	0.	22	0.	21
As	< ПО		< ПО		0.	29	2.	1	0.	93	4.	2	2.	1
Rb	0.	17	0.	15	6.	9	0.	18	0.	38	0.	57	0.	20
Sr	139		114		41.	4	108		118		64.	2	65.	5
Y	2.	8	4.	2	10.	8	0.	25	0.	44	0.	89	0.	43
Zr	27.	3	12.	6	129		4.	0	20.	9	2.	6	2.	3
Nb	1.	2	0.	75	6		0.	12	0.	40	0.	36	0.	32
Mo	0.	13	0.	32	1.	6	0.	19	0.	68	1.	0	1.	0
Ag	0.	065	< ПО		0.	052	0.	28	0.	057	< ПО		< ПО
Cd	< ПО		< ПО		< ПО		0.	19	0.	26	0.	063	0.	062
Sn	0.	15	0.	082	0.	57	0.	75	0.	51	0.	72	0.	72
Sb	0.	17	< ПО		0.	15	0.	26	0.	27	0.	89	0.	22
Cs	0.	23	0.	077	0.	36	0.	040	0.	042	0.	017	0.	014
Ba	7.	9	8.	5	24.	9	3.	0	5.	8	5.	5	2.	3
La	1.	9	1.	5	12.	2	0.	16	0.	13	0.	72	0.	38
Ce	4.	0	4.	0	29.	7	0.	31	0.	32	1.	5	0.	64
Pr	0.	46	0.	41	3.	0	0.	042	0.	035	0.	15	0.	060
Nd	1.	9	1.	8	11.	4	< ПО		0.	088	0.	77	0.	32
Sm	0.	42	0.	45	2.	1	0.	035	0.	025	0.	14	0.	050
Eu	0.	076	0.	083	0.	54	0.	018	0.	010	0.	045	0.	026
Gd	0.	36	0.	44	2.	0	0.	011	0.	038	0.	17	0.	059
Tb	0.	059	0.	075	0.	32	< ПО		< ПО		0.	022	0.	0080
Dy	0.	41	0.	53	2.	2	0.	025	0.	023	0.	14	0.	051
Ho	0.	085	0.	12	0.	41	0.	0068	0.	0076	0.	025	0.	010
Er	0.	34	0.	43	1.	3	0.	022	0.	028	0.	080	0.	033
Tm	0.	040	0.	050	0.	21	0.	02	0.	0038	0.	0093	< ПО
Yb	0.	28	0.	33	1.	5	0.	021	0.	029	0.	067	0.	030
Lu	0.	039	0.	042	0.	23	< ПО		0.	075	0.	041	0.	043
Hf	0.	61	0.	27	3.	6	0.	08	0.	24	0.	050	0.	034
Ta	0.	075	0.	060	0.	39	< ПО		0.	13	0.	077	0.	15
W	0.	18	0.	34	1.	9	0.	84	0.	69	1.	5	1.	6
Re	< ПО		< ПО		< ПО		< ПО		< ПО		0.38		0.	40
Pt	< ПО		< ПО		< ПО		0.	18	0.	034	< ПО		< ПО
Pb	7.	7	3.	0	15.	8	1.	1	1.	4	1.	9	1.	9
Bi	< ПО		< ПО		0.	032	< ПО		0.	0085	< ПО		< ПО
Th	0.	29	0.	61	0.	49	0.	024	< ПО		0.	11	0.	030
U	0.	68	1.	7	1.	5	0.	16	2.	36	1.	2	1.	2

Примечание. Светлым тоном отмечены повышенные содержания, темно-серым — концентрации меньше предела обнаружения. Содержания оксидов даны в мас.%, элементов – в г/т. Содержания воды найдены методом KFT. *Голубой диопсид — — другой кусочек. **Осветленная зона — две навески из одного образца.

 

Голубой цвет диопсида обусловлен высоким содержанием ванадия в форме ванидил- иона (VO²⁺), замещающего кальций. Ванадий в природе присутствует в основном в степени валентности три и четыре (V³⁺, VО²⁺). Эти валентные формы как в растворах, так и в составе минералов имеют характерную окраску: V³⁺ – зеленую, ванидил-ион (VО²⁺) – синюю (например, (Satake, Mido, 2010)). В природе встречаются ванадиевые минералы кавансит и пентагонит Ca(VO)Si₄O₁₀·4(H₂O), которые имеют цвет, близкий к цвету голубого диопсида. Соотношение V³⁺/V⁴⁺ зависит от летучести кислорода, что позволяет использовать коэффициент распределения ванадия между кристаллическими фазами (оливином, шпинелью) и расплавом в качестве оксобарометра. Коэффициент распределения между оливином и расплавом (K_D,V) в зависимости от летучести кислорода близок к 0.25, что отвечает стехиометрии реакции окисления (V³⁺ + 1/4О₂ = V⁴⁺ + 1/2O^2-) (Canil, Fedortchouk, 2001). Таким образом, чтобы соотношение V³⁺/V⁴⁺ изменилось на порядок, летучесть кислорода должна измениться на 4 порядка (например, от NNO + 3 до QFM).

Среднее содержание ванадия в кристаллах диопсида (16 точек), найденное методом съемки на спектрометре, составило 345 г/т. Съемка по профилю свидетельствуют о зональном распределении ванадия в диопсиде с уменьшением содержания до величин меньше предела обнаружения на краях кристалла с осветленной зоной и с межзерновыми границами. Часть зон, видимых на BSE изображениях, коррелирует с содержанием ванадия, а часть – нет (рис. 2а). Можно отметить, что в некоторых точках соизмеримое и большое содержание в диопсиде имеет титан. Содержания железа, хрома, марганца намного меньше. Оба анализа диопсидового прослоя методом ICP-MS дали близкие содержания 300–350 г/т ванадия. Содержание ванадия в перекристаллизованном диопсиде в агрегате с монтичелитом находится на уровне предела обнаружения микрозонда.

 

 

Рис. 2. BSE-изображения диопсида из скарна. (а) – кристаллы зонального голубого диопсида, стрелкой отмечено положение профиля, по которому измерено содержание ванадия. Поздний ранкинит заполняет трещины между кристаллами диопсида и замещает мелкозернистую зону с монтичеллитом. (б) – профиль распределения содержания ванадия в кристалле голубого диопсида на профиле (см. стрелку на рис. 2а), концентрация ванадия снижается к краям кристалла. (в) – зона перекристаллизации диопсид-монтичеллитового состава с наложенным на регрессивной стадии ранкинитом. Безводный фтористый куспидин (Ca4(Si2O7)F2) образует включения в диопсиде и на границе кристалла ранкинита. Mts – монтичеллит, Rnk – ранкинит, Di – диопсид, Csp – куспидин (Ca4(Si2O7)(ОН,F)2).

 

Состав скарнов и приконтактных пород по микроэлементам. Анализ методом ICP-MS позволяет расширить число проанализированных элементов (табл. 2). Помимо отмеченной выше обогащенности диопсидовых скарнов ванадием, можно также отметить (см. табл. 1) повышенные содержания платины (0.18 и 0.034 г/т) и серебра (0.28 и 0.057 г/т). Причем, если содержание ванадия в среднем сохраняется на уровне около 300 г/т, то содержания платины и серебра испытывают значительные синхронные вариации. Сходно ведет себя мышьяк (2.1 и 0.93 г/т).

Содержание кадмия также максимальны в диопсидовой зоне, но не коррелированны с платиной (0.19 и 0.26 г/т). Осветленная зона характеризуется низким содержанием ванадия (7–6 г/т) и относительно высоким хрома (98–88 г/т). Она обогащена рением (0.38, 0.40 г/т), содержание которого в остальных изученных породах ниже предела обнаружения, при этом она обеднена серебром (ниже 0.02 г/т) и платиной (ниже 0.006 г/т). Также можно отметить ее относительное обогащение молибденом (1 против 0.6 г/т в голубом диопсиде) и вольфрамом (1.5 против 0.7–0.8 г/т). Другой особенностью светлой зоны являются низкое содержание ниобия и значение отношения Nb/Ta = 3 против нормальных Nb/Ta = 13–15 во всех остальных изученных породах (карбонаты, кварцит, голубой диопсид). Эту особенность можно связать с влиянием фтора, который образует более прочные соединения с ниобием и выносит его в большей степени. Карбонатная порода, наиболее удаленная от контакта, и кварцит имеют относительно большие концентрации цинка (26.7, 34.2 г/т) и серебра. В них также выше концентрация свинца (7.7, 15.8 г/т), а в кварците наблюдаются максимальные содержания лития (19.8 г/т) и хрома (261 г/т).

В общем виде изученные составы скарнов и вмещающих пород отражены на рис. 3. Примечательно, что редкоземельные спектры всех изученных образцов подобны и имеют плоские (с незначительным повышением нормированной концентрации Er у карбонатов) участки в области HREE и линейно растущий участок для LREE (рис. 3а). Общее содержание редкоземельных элементов наименьшее в диопсиде и наивысшее в кварците. Для сравнения на рис. 3а нанесены редкоземельные спектры магматических пород Йоко-Довыренского массива (Арискин и др., 2015). Редкоземельные спектры кварцита и ультраосновных пород Йоко-Довыренского массива практически идентичны. Также для сравнения нанесены данные по составу экспериментальных расплавов, сформировавшихся при участии углекислого восстановленного флюида при Р = 2 кбар иТ = 1000^оС (Симакин и др., 2018), их спектры практически идентичны спектрам габбро Йоко-Довыренского интрузива. Для более широкого обзора полученные составы пронормированы на средний состав коры (по https://earthref.org/GERM), а не примитивной мантии или хондрита и отображены на рис. 3б. При такой нормировке форма спектров по HFSE и REE стала близка к прямой линии и для магматических пород Йоко-Довыренского плутона (причем, кварцит и в этом представлении практически идентичен ультрабазитам Йоко-Довыренского массива), и для изученных пород. В то же время для LILE и Pb, Th, U характерны высокие вариации концентраций. Показано, что магматические породы Йоко-Довыренского интрузива обогащены коровыми архейскими изотопами (Арискин и др., 2015). Изложенные данные можно интерпретировать как проявление выравнивания состава осадочных и магматических пород по микроэлементам в потоке флюидов мантийного происхождения, насыщенных кратонным компонентом из-за обрыва плиты при коллизии. На важную роль глубинного флюида в рассматриваемом регионе указывают данные изотопного состава скарнов в эндоконтакте сиенита Тажеранского массива (Приольхонье, Прибайкалье). Установлено, что термическое разложение карбонатных пород происходило при участии стороннего, не связанного с разложением карбонатов (корового или мантийного), восстановленного углерода (Doroshkevich et al., 2017).

 

Рис. 3. Геохимические характеристики изученных образцов и магматических пород Йоко-Довыренского массива: (а) – нормированные на состав хондрита (Sun, McDonough, 1989) концентрации редких земель в магматических породах Йоко-Довыренского плутона (Арискин и др., 2015), экспериментальных стеклах (Симакин и др., 2018) и образцах из данной работы; (б) – спайдер-диаграмма с нормировкой на средний состав земной коры по (https://earthref.org/GERM/).

 

Оценка глубины становления Йоко-Довыренского плутона по минеральным барометрам. Для анализа условий образования скарнов необходима независимая (от предполагаемых условий скарноообразования) оценка давления становления интрузива. Опубликованные данные по составу сосуществующих орто- и клинопироксена в плагиоперидотитовой зоне интрузива (Кислов, 1998) позволяют оценить давление по двупироксеновому барометру (Putirka, 2005). Лишь одна пара составов из табл. 25 (Кислов, 1998) удовлетворяет условию равновесия и дает оценку давления в 2.8 кбар. Нами проанализированы амфиболы из плагиоклазовых лерцолитов плагиоперидотитовой зоны приподошвенной части Йоко-Довыренского массива. Амфиболы подвержены вторичной хлоритизации. Выбранные анализы первичного амфибола (см. табл. 2) по составу отвечают кальциевому амфиболу, иногда с несколько повышенным содержанием щелочей (Na₂O + K₂O до 3.65 мас.%). Отличительной чертой этих амфиболов являются необычно высокие содержания ванадия (до 0.25 мас.% V₂O₃) и хрома (до 2 мас.% Cr₂O₃). Для оценки давления мы использовали эмпирический барометр (Симакин, Шапошникова, 2017), который откалиброван для высокомагнезиальных андезитовых и базитовых расплавов. Совокупность из восьми анализов в координатах (Al^VI, Ti + Fe³⁺ + V + Cr) характеризуется ограничивающей ее сверху прямой (см. рис. 4), которая точно проходит через контрольную точку (Al^VI = 0, Ti + Fe³⁺ + V + Cr = 1.8) и отсекает на оси ординат Al^VIмакс = 0.6, что отвечает давлению около 2 кбар. Содержание многозарядных катионов в среднем составляет лишь 0.4, что отражает относительно низкую летучесть кислорода при кристаллизации амфибола (предположительно в районе QFM). Таким образом, глубина становления интрузива была не менее 7 км, а не 1.5 км, как считалось ранее (Перцев, Шабынин, 1978), что необходимо учитывать при термодинамических оценках условий формирования пород Йоко-Довыренского интрузива.

ТЕРМОДИНАМИЧЕСКАЯ МОДЕЛЬ

Представленные данные можно интерпретировать, исходя из результатов термодинамического моделирования состава флюида при формировании изученных образцов голубого диопсида. Мы рассчитали термодинамику процесса скарнообразования до температуры плавления кальцита в приближении инертных MgO, CaO. При значительном переносе кремнезема, как в рассматриваемом случае, этот подход является более адекватным для описания равновесия при вполне подвижном поведении SiO₂ и СО₂ (Ferry et al., 2011). Непосредственное превращение доломита в диопсид, минуя форстерит-кальцитовый парагенезис, возможно при относительно низкой температуре, относительно высокой мольной доле СО₂ во флюиде и высокой активности кремнезема (рис. 5а). Диаграмма на рис. 5а близка к таковой на fig. 8а из работы (Ferry et al., 2011). При высокой температуре (900^оС) образование монтичеллита из диопсида возможно при выносе кремнезема, обусловленном снижением активности кремнезема во флюиде при относительно высокой активности СО₂ (0.5–0.3) (рис. 5б). При небольшой температуре равновесие Di–Mts возможно лишь при очень низкой активности СО₂ (меньше 0.01). Наблюдения показывают, что температура не достигала значений свыше 1000–1100^оС, что привело бы к формированию чисто магнезиальных скарнов из-за плавления кальцита. Углекислота, образующаяся при реакции разложения доломита, является сильным окислителем, в связи с чем ванадий образует ванидил-ион VO²⁺, который занимает позицию кальция в диопсиде и придает ему голубой цвет. Со временем происходит восстановление водно-углекислого флюида с повышением его температуры за счет взаимодействия с окружающими кумулатами и ультраосновной магмой. Обнажение апокарбонатных пород с голубым диопсидом имеет размеры несколько метров в сечении, поэтому их прогрев при начальной и, возможно, последующих инъекциях магмы был постепенным. Данных для описания детальной термической истории и схемы перемещения флюида у нас нет. Высокая активность кремнезема на начальном этапе обеспечивается растворением кремневых стяжений в доломите (Кислов, 1998). На высокотемпературной стадии активность кремнезема буферируется растворимостью темноцветных минералов (прежде всего, оливина) из кумулуса ультраосновной магмы, основного интерстициального расплава. В скарнах на контакте с кислыми магмами (адамелитами) активность кремнезема во флюиде выше, и парагенезисов с монтичеллитом в них не наблюдается (Ferry et al., 2011).

 

Рис. 4. Состав амфиболов Йоко-Довыренского массива: звезды – из дунитов и плагиперидотитов по (Кислов, 1995), ромбы из плагиолерцолитов – на диаграмме октаэдрический алюминий – сумма высокозарядных катионов по (Симакин и Шапошникова, 2018). Звезда, отмеченная кружком, отвечает составу амфибола выросшего из раствора (автометасоматического или метаморфического происхождения).

 

Рис. 5. Фазовые диаграммы диопсид-монтечеллитовых скарнов при вполне подвижном поведении кремнезема,Р = 2 кбар: (а) при Т = 600оС, (б) при Т = = 950оС. Mts – монтичеллит, Per – периклаз, Cal – – кальцит, Dol –доломит, Fo – форстерит, Di – диопсид.

 

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

Запишем реакцию окисления оливина в присутствии водно-углекислого флюида в общем виде как:

Fe₂SiO_4SS +1/3O₂ = 2/3Fe₃O_4SS + SiO_2fl, (2)

где фаялитовый минал образует твердый раствор с форстеритом, а кремнезем растворен во флюиде. Растворимость кремнезема зависит от состава флюида и Р-Т параметров, и вместе с активностью фаялита в оливине однозначно задает (буферирует) в равновесии летучесть кислорода. Валовый состав флюида (например, в системе C-S-H-O) и летучесть кислорода, в свою очередь, определяют соотношение между различными его компонентами.

Концентрация трехвалентного железа (по данным монографии Кислова (1998)) в части шпинелей ультрамафитовой части разреза Йоко-Довыренского массива очень мала, отвечающее ей содержание магнетитового минала составляет 0.02 ± 0.008 (X_Mag, мол. доля). В то же время в основном встречаются более окисленные шпинели с X_Mag = 0.15 ± 0.06. Самые окисленные шпинели характеризуются X_Mag= 0.32 ± ± 0.06. Железистость оливинов в нижней ультрамафитовой плагиоперидотитовой и плагиодунитовой части разреза (по данным монографии Кислова (1998)) составляет f = 13.89 ± 1.66, в нижней части дунитового разреза f = 12.57 ± ± 1.80. Используем эти данные по составам фаз для оценки летучести кислорода в равновесии (2) при Т = 950^оС и Р = 2000 бар.

Для расчета активности магнетитового минала мы использовали термодинамическую модель шпинелевого твердого раствора (Sack, Ghiorso, 1991). Активности рассчитаны для средних по группам составов шпинелей с использованием on-line калькулятора, в котором реализована последняя версия модели (http://melts.ofm-research.org/CalcForms/index.html). Коэффициент активности при больших содержаниях магнетита равен 0.5–0.6, для самых восстановленных шпинелей он снижается до 0.08. Таким образом, при расчетах мы применяли три значения активности магнетита (а_Mag): 0.002, 0.07 и 0.19. Для расчета активности фаялитового компонента в магнезиальном оливине нами использовалась модель твердого раствора (Berman, Aranovich, 1996). Активности оливина, отвечающие средним составам, представленным выше, составляют (a_Fa): 0.038 и 0.044. Согласно этим оценкам летучесть кислорода при log(a_Qz) = –1.3 в равновесии с наиболее распространенными шпинелями (X_Mag= 0.15) составляет QFM-2.2 (QFM: fO₂ = (–25096.3/T)+ + 8.735 + 0.11(P – 1)/T, где температура по Кельвину, давление в барах (Frost, 1991). Эта оценка несколько ниже, чем летучесть кислорода, отвечающая буферу ССО (QFM-2.0 при Р = 200 МПа и Т = 950^оС (French, Eugster, 1965)). При наличии СО₂ во флюиде это значение является нижним пределом летучести кислорода. Наиболее окисленные шпинели отвечают fO₂ около QFM-1.3. Наиболее восстановленные шпинели, вероятно, являются первичными магматическими шпинелями, не вступавшими в реакцию с флюидом, формально рассчитанная по реакции (2) равновесная fO₂ составляет фантастически низкую QFM-5.3, что ниже летучести буфера IW.

Предсказываемая низкая летучесть кислорода качественно согласуется с составом амфибола (см. выше). Низкая летучесть кислорода в близком экзоконтакте интрузива оценена по паре ильменит—магнетит из карбонатов (составы см. в табл. 1), согласно расчетам по программе ILMAT (Lepage, 2003) минералы уравновешены при T = 398^oC и fO₂ около QFM-1.9.

Низкая летучесть кислорода в оливиновом кумулусе приводит к появлению восстановленных газов во флюиде. Пренебрегая углеводородами с большим молекулярным весом в системе С-О-Н можно ограничиться следующими компонентами: H₂O, CO₂, H₂, CH₄, CO. Летучесть кислорода будем считать заданной буферной реакцией (2). СО является активным компонентом флюида, образующим устойчивые p-комплексы с переходными металлами, включая ЭПГ. Метан в этом смысле является пассивным компонентом. С увеличением температуры происходит термическое разложение метана, а концентрация СО возрастает. На рис. 6 представлены термодинамические расчеты зависимости концентраций СО и СН₄ (мол. доли) от температуры для Р = 2 кбар и летучести кислорода fO₂ = QFM и QFM-1. Также для сравнения представлены данные для флюида при давлении 1 атм в равновесии с графитом (буфер ССО), в котором достигается максимальная концентрация СО и минимальная – метана. Отношения СО/СН₄ = 1 достигаются при более высокой температуре в более восстановительных условиях (Т » 750^оС при fO₂ = QFM-1). Таким образом, при Р = 2 кбар (как в магматической камере Йоко-Довырена) можно предполагать, что значительные концентрации СО ожидаются при высоких температурах, близких к температуре солидуса.

Сульфидное равновесие. Рассмотрим, как будет происходить реакция водно-углекислого флюида с пирротином, главным сульфидным магматическим минералом, который концентрирует ЭПГ практически полностью при отделении сульфидного расплава от силикатного. Вода и СО₂ являются сильными окислителями, способными окислять сульфиды. При реакции сульфидов с флюидом металлы платиновой группы извлекаются в раствор и могут перемещаться в толще кумулуса и накапливаться. Реакция пирротина с водно-углекислым флюидом промоделирована нами термодинамически при Р = 2000 бар и Т = 950^оС. Принято, что во флюиде в равновесии с FeS содержатся следующие компоненты системы С-О-S-Н: 1) H₂O, 2) H₂, 3) CO₂, 4) CO, 5) H₂S, 6) SO₂, 7) COS, 8) O₂, 9) S₂. Для простоты пренебрегаем образованием метана и тиолов. Для того чтобы найти концентрации девяти компонентов, мы используем уравнения действующих масс шести реакций:

1) 3FeS + 2O₂ = Fe₃O₄ + 3/2S₂,

2) H₂ + 1/2O₂ = H₂O,

3) CO + 1/2O₂ = CO₂,

4) H₂S = 1/2 S₂ + H₂,

5) CO + 1/2S₂ = COS,

6) 1/2S₂ + O₂ = SO₂.

 

Рис. 6. Рассчитанные в зависимости от температуры мольные доли метана и СО в восстановленном флюиде СО2-Н2О при начальном составе флюида XН2О = XCO2 = 0.5, P = 200 MПa, летучести кислорода указаны на графике. Для сравнения нанесены концентрации восстановленных форм углерода в водно-углекислом флюиде в равновесии с графитом при давлении 1 бар. Точки на графике маркируют равенство концентраций СО и СН4.

 

Помимо этого учитываются два уравнения баланса масс по углероду (уравнение 6) и водороду (уравнение 7). Количество кислорода и серы во флюиде переменно из-за реакции окисления пирротина. В пределах закрытой системы количество кислорода, фиксированное в магнетите, соответствует количеству серы, перешедшей во флюид, и составляет третье условие (уравнения 9).

7) XCOS + XCO₂ + XCO = 0.5,

8) XH₂O + XH₂S + XH₂ = 0.5,

9) 3/4(0.5 – SX_inO,i ) = SX_inS,i,

где Xi – мольная доля i-компонента, nO,i nS,i – число атомов кислорода и серы в компоненте i. Шесть уравнений закона действующих масс для реакций (1–6) и три условия баланса масс достаточно для того, чтобы определить концентрации компонентов флюида. При вычислениях кислород учитывался в уравнениях равновесия в для переменной летучести (fO₂). В балансе масс его количеством пренебрегали, поскольку его мольная концентрация ничтожно мала. Константы уравнений (1–6) вычислялись с учетом коэффициентов фугитивности газовых компонентов во флюиде. Система уравнений численно решена в среде MAPLE. Установлено, что при равновесии летучесть кислорода равна fO₂ = QFM + + 0.72. В реакции с избытком пирротина расходуется 4% воды и 1.8% СО₂. При этом образуется флюид, содержащий примерно 2 мол.% серосодержащих компонентов (H₂S — — 1.7, COS – 0.3, SO₂ – 0.1, S₂ – 0.03 мол.%). Содержание СО достигает 0.6 мол.%, а Н₂ – 0.24 мол.%.

Если вместо условия (9) принять, что летучесть кислорода равна некоторой фиксированной величине, то можно найти состав флюида в системе, открытой по кислороду. Примем fO₂ = QFM-2, т.е. равной величине, которая буферируется реакцией флюида с оливиновым кумулусом. Как уже упоминалось выше, эта летучесть также близка к значению буфера ССО при Р-Т условиях кумулусного процесса в Йоко- Довыренском интрузиве. В этих условиях метаном нельзя пренебрегать, и мы включили его в список компонентов. Суммарная концентрация серосодержащих газов в восстановленном флюиде на контакте с пирротином снижается до 0.56 мол.% (H₂S – 0.48, COS – 0.08, S₂ – 0.67·10^-7, SO₂ – 0.3·10^-7 мол.%). Концентрация СО во флюиде достигает 10.4 мол.%, метана – 3.1 мол.%, а водорода – 4.4 мол.%. Согласно нашим экспериментальным данным, во флюиде с такой высокой концентрацией монооксида углерода наблюдается высокая растворимость платины (Simakin et al., 2016).

При летучести кислорода несколько выше буфера QFM содержание серы во флюиде достаточно высоко, помимо сероводорода в нем могут присутствовать серо-органические соединения тиолы, являющиеся спиртами, в которых кислород замещен серой (таких как CH₃SH). Устойчивость соединений благородных металлов с тиолами при высоких Р-Т параметрах не изучена. При низких Р-Т параметрах тиолы образуют устойчивые комплексы с Cd, Ag, Zn (https://en.wikipedia.org/wiki/Thiol), т.е. с металлами, характерными для Холоднинского месторождения, расположенного вблизи северо-восточной границы Йоко-Довыренской интрузии (Дистанов и др., 1995).

МЕХАНИЗМ ВОССТАНОВЛЕНИЯ ФЛЮИДА В КУМУЛУСЕ УЛЬТРАОСНОВНЫХ ИНТРУЗИВОВ

Таким образом, полученные данные можно интерпретировать следующим образом. На начальной стадии скарнообразования происходил перенос флюидом кремнезема из стяжений кремня в доломитах с образованием диопсида. На второй стадии флюид выносил кремнезем с образованием диопсид-монтичеллитовой зоны. На заключительной регрессивной стадии флюид, обогащенный фтором и кальцием, вызвал кристаллизацию ранкинита. Окисленная форма ванадия в ванидил-ионе (VO²⁺) связана с изначально высокой летучестью кислорода в СО², выделяющегося при реакциях декарбонатизации (например, (Simakin et al., 2012)). Можно предположить, что выделяющийся углекислый флюид восстанавливается за счет реакции с кумулусным оливином и сульфидами и насыщается платиной, серебром, халькофильными элементами (As, Cd, Se). При фильтрации восстановленного углекислого флюида через диопсидовый скарн при 700 < Т < 950^oC происходило формирование осветленной зоны. Граница между зонами резкая. Из флюида, проникающего в зону голубого диопсида, происходило осаждение платины и серебра. В изученных образцах голубого диопсида концентрации Pt и Ag (а также мышьяка) коррелируют между собой при постоянном содержании ванадия. В осветленной зоне концентрация ванадия снижается до 7–6 г/т. Низкие концентрации платины, серебра и ванадия в осветленной зоне, видимо, свидетельствуют о том, что растворимость этих элементов во флюиде была велика. Необычный состав амфиболов из плагиолерцолитов (см. выше) подтверждает высокие концентрации ванадия и хрома во флюиде при температуре около 1000^оС (температура кристаллизации амфибола). Кристаллизация амфибола также свидетельствует о достаточно высокой доле воды в восстановленном флюиде. Особенности флюидно-магматического, флюидного раннего постмагматического взаимодействия важны для понимания механизма формирования оруденения платинометальных рифов малосульфидного типа Йоко-Довыренского массива (Кислов, 1998). Летучесть кислорода, оцененная по оливинам из ультраосновной части массива методом высокотемпературной электрохимии, существенно меньше QFM (до QFM-4) при субсолидусной температуре около 1000^оС (Арискин и др., 2017), что согласуется с представленными выше термодинамическими оценками. Можно предположить, что состояние минералов фиксирует субсолидусное переуравновешивание в восстановительных условиях.

Предложенный механизм восстановления флюида также подтверждается данными по составу шпинелей из плагиоперидотитов и плагиодунитов Йоко-Довыренского интрузива (Кислов, 1998). Как было отмечено Г.С. Николаевым (Николаев и др., 2017), при субсолидусных условиях шпинель окисляется и обогащается магнетитовым миналом (см. рис. 7а). Параллельно с перекристаллизацией шпинели обогащаются Zn, V, а также в некоторых случаях никелем, что можно отметить на рис. 7б.

 

Рис. 7. Составы шпинелей из базальной серии плагиоперидотитов Йоко-Довыренского массива (Кислов, 1995). При субсолидосной реакции оливинов с водно-углекислым флюидом: (а) – шпинели обогащаются магнетитом (с хромом), этот тренд отмечен штриховой линией; вариации вдоль основания Cr-Al отвечают направлению тренда фракционной магматической кристаллизации по (Николаев и др., 2017); (б) – взаимодействие с флюидом иногда приводит к обогащению шпинелей V, Zn, концентрация никеля также меняется.

 

Перечень элементов, концентрации которых испытывают значительные вариации в изученных скарнах и осадочных протолитах, а значит обладают хорошей растворимостью в восстановленном углеродсодержащем флюиде: Ni, V, Mo, Cr, Re, примечателен тем, что все эти элементы имеют устойчивые карбонилы при низких Р-Т параметрах. Их устойчивость при высоких Р-Т параметрах предстоит изучить. Растворимость платины в водном флюиде в виде хлоридных комплексов велика при высокой летучести кислорода (Mungal, Naldrett, 2008) и при низкой летучести кислорода в углеродсодержащем флюиде (карбонилы, возможно с хлором и серой). При формировании рифа Меренского, как и платиноносного рифа Йоко-Довыренского массива, летучесть кислорода оценивается на уровне QFM-2 (Ballhaus, Sylvester, 2000), что говорит о вероятном формировании карбонилов. Предложенный в данной работе механизм восстановления углекислого флюида при взаимодействии с оливиновым кумулатом при высокой температуре может быть актуален как для расслоенных ультрабазитовых интрузивов (Бушвельд и Йоко-Довыренский), так и для мантии.

ВЫВОДЫ

Летучесть кислорода во флюидах, образующихся при взаимодействии доломитов и ультраосновных магм Йоко-Довыренского массива, имеет значительный диапазон вариаций. 1) Голубой диопсид образуется на начальной стадии скарнообразования, которая характеризуется высокой летучестью кислорода выделяющегося СО₂, обеспечивающей образование VO²⁺ иона. 2) Согласно термодинамическим расчетам, восстановление водно-углекислого флюида до уровня около QFM-2 обеспечивается реакцией с оливиновым кумулатом. Диопсид образуется из доломита при высокой активности кремнезема во флюиде, обеспеченной стяжениями кремнезема в доломитах. Низкая активность кремнезема в водно-углекислом флюиде, буферируемая реакцией с ультраосновными породами, вызывает переход диопсида в монтичеллит. При реакции с оливиновым сульфидсодержащим кумулусом флюид восстанавливается с образованием СО и H₂S, а также, возможно, сероорганических соединений. Предполагается, что голубой диопсид обогащается платиной до 0.2 г/т в связи с проработкой водно-углекислым флюидом, выделившимся при реакциях декарбонатизации и экстрагировавшим металл при взаимодействии с магмой и магматическим кумулусом.

Благодарности. Авторы благодарны Ю.К. Карандашеву (ИПТМ РАН) за анализы природных образцов методом ICP-MS. Признательны Л.Я. Арановичу (ИГЕМ РАН) и А.С. Мехоношину (ИГХ СО РАН) за содержательные рецензии работы. Особенно хочется отметить квалифицированную критику Л.Я. Арановича, позволившую кардинально улучшить качество термодинамических расчетов.

Источники финансирования. На заключительном этапе работа была поддержана грантом РФФИ № 18-05-00597.

About the authors

A. G. Simakin

D.S. Korzhinskii Institute of Experimental Mineralogy of Russian Academy of Sciences; Institute of the Earth Physics of the Russian Academy of Sciences

Author for correspondence.
Email: simakin@iem.ac.ru
Russian Federation, 4, Academica Osypyana ul., Chernogolovka, Moscow region, 142432; Bolshaya Gruzinskaya str., 10-1, Moscow 123242, Russia

E. V. Kislov

Geological Institute, Siberian Branch of the Russian Academy of Sciences

Email: simakin@iem.ac.ru
Russian Federation, 6A, Sakhyanova street, Ulan-Ude, 670047

T. P. Salova

D.S. Korzhinskii Institute of Experimental Mineralogy of Russian Academy of Sciences

Email: simakin@iem.ac.ru
Russian Federation, 4, Academica Osypyana ul., Chernogolovka, Moscow region, 142432

O. Yu. Shaposhnikova

D.S. Korzhinskii Institute of Experimental Mineralogy of Russian Academy of Sciences

Email: simakin@iem.ac.ru
Russian Federation, Bolshaya Gruzinskaya str., 10-1, Moscow 123242

A. N. Nekrasov

D.S. Korzhinskii Institute of Experimental Mineralogy of Russian Academy of Sciences

Email: simakin@iem.ac.ru
Russian Federation, 4, Academica Osypyana ul., Chernogolovka, Moscow region, 142432

References

Supplementary files