CC BY
11
УДК 550.4: 548.4
ЭВОЛЮЦИЯ ФИЗИКО-ХИМИЧЕСКИХ УСЛОВИЙ КРИСТАЛЛИЗАЦИИ РАСПЛАВОВ ПРИ ФОРМИРОВАНИИ ДУНИТ-ПЕРИДОТИТ-ГАББРОВЫХ МАССИВОВ ВОСТОЧНОГО САЯНА
Ступаков Сергей Иванович1,
stupakov@igm.nsc.ru
Симонов Владимир Александрович1,2,3,
simonov@igm.nsc.ru
Мехоношин Алексей Сергеевич4,5,
mekhonos@igc.irk.ru
Колотилина Татьяна Борисовна4,5,
tak@igc.irk.ru
1 Институт геологии и минералогии им. В.С. Соболева Сибирского Отделения Российской Академии Наук, Россия, 630090, г. Новосибирск, пр. Академика Коптюга, 3.
2 Новосибирский национальный исследовательский государственный университет, Россия, 630090, г. Новосибирск, ул. Пирогова, 1.
3 Казанский федеральный университет, Россия, 420008, Казань, ул. Кремлевская, 18.
4 Институт геохимии им. А.П. Виноградова Сибирского Отделения Российской Академии Наук, Россия, 664033, г. Иркутск, ул. Фаворского, 1А.
5 Иркутский национальный исследовательский технический университет, Россия, 664074, г. Иркутск, ул. Лермонтова, 83.
Актуальность. Изучение расплавных включений в минералах является одним из наиболее достоверных путей получения прямой информации о характере магматических процессов, позволяет построить модель образования интрузий и проследить эволюцию физико-химических условий кристаллизации магматических комплексов.
Цель: получить достоверные сведения о составе расплава, что дает возможность проследить эволюцию физико-химических условий образования пород платиноносного дунит-верлитового массива Медек.
Объекты: расплавные включения в хромшпинелидах дунитов и верлитов Медекского базит-ультрабазитового массивов. Методы. Поиск расплавных включений проводился в монофракциях акцессорных хромшпинелидов из дунитов и верлитов массива. Полированные препараты просматривались под микроскопом, и отбирались зерна, содержащие многофазные микровключения определённой формы, размера (10-30 мкм) и положения внутри зерна. Далее проводились высокотемпературные эксперименты. Выбор температурного режима (1280-1300 °С) определялся задачей перевести включения в расплавленное состояние, затем его закалить, и полученное стекло проанализировать на рентгеновском микроанализаторе и сканирующем микроскопе.
Результаты. В результате анализа стекол гомогенизированных расплавных включений установлено, что состав родоначальной магмы, из которой кристаллизовались дуниты массива, отвечал пикритоидному расплаву (MgO 13,36 мас. %), который эволюционировал до низкокалиевого пикробазальтового расплава, наиболее близкого по петрохимическим характеристикам к данным по расплавным включениям в хромитах из океанических комплексов. Геотермобарометрические расчеты и моделирование с помощью программ COMAGMAT и PETROLOG свидетельствуют о становлении пород массива при давлении 3-5 кбар и температурном диапазоне 1400-1230 °С. Расчетные данные и результаты высокотемпературных экспериментов показали хорошую сходимость.
Ключевые слова:
Медекский массив, ультрабазиты, дуниты, верлиты, хромшпинелиды, расплавные включения.
Введение
При исследовании рудоносности базит-ультраба-зитовых комплексов большое внимание уделяется генетическим проблемам [1-4]. На оценку перспектив платиноносности, прежде всего, влияют представления об условиях формирования пород [5-8] и составе исходной магмы. Ранее было обосновано, что изучение расплавных включений в минералах является одним из наиболее достоверных путей получе-
ния прямой информации о характере магматических процессов [9-20]. На примере ряда массивов Сибирского кратона [13, 14, 16, 21], в том числе и ду-нит-перидотит-габбровых массивов Восточного Сая-на [22, 23], ранее была доказана эффективность применения разработанной нами методики [14] исследования расплавных включений в хромшпинелидах при изучении рудно-магматических систем, связанных с базит-ультрабазитовыми комплексами.
208
DOI 10.18799/24131830/2019/7/2196
В настоящей статье представлены новые данные по минералам и расплавным включениям в хромшпинелидах дунит-верлитового массива Ме-дек, результаты компьютерной обработки оригинальной информации по составам минералов и включений, с помощью минералогических термометров и барометров, а также программ COMAG-МАТ [24-26] и PETROLOG [27], которые позволили выяснить особенности эволюции физико-химических параметров магматических систем при его формировании.
Геолого-петрографические особенности пород
массива Медек
Массив Медек расположен в центральной части Восточного Саяна, в междуречье Мал. Бирюса -Уда (рис. 1, врезка). Форма массива в плане эллипсовидная, длинная ось вытянута в субмеридиональном направлении, площадь выхода составляет около 2 км2 (рис. 1). Хотя непосредственные границы с вмещающими породами не наблюдаются, о тектоническом характере взаимоотношений говорят отрицательные формы рельефа, развитые на контактах зоны милонитизации и катаклаза, сопровождающиеся метаморфогенным преобразованием пород.
Рис. 1. Местоположение (врезка) и схематическая геологическая карта массива Медек (рисунок составлен на основе оригинальных данных с использованием материалов из работ [22,28]). 1 - граниты, сиениты, диориты; 2 - уль-табазиты; 3 - гнейсы, сланцы, амфиболиты; 4 - положение скважин
Fig. 1. Location (inset) and schematic geological map of the Medek massif (figure is based on the original data using materials from the works [22,28]). 1 - granites, syenites, diorites; 2 -ultrabasites; 3 - gneisses, shales, amphibolites; 4 - the position of boreholes
Массив Медек сложен дунит-верлитовой дифференцированной серией пород с постепенными переходами между петрографическими разностями. Дуниты массива среднезернистые, микроструктура пород кумулятивная с чётко проявленным идиоморфизмом зёрен оливина. Интеркуму-
лусный пироксен по составу отвечает хромдиопси-ду. Главным акцессорным минералом является хромшпинелид. Для верлитов характерны кумулятивные и пойкилитовые структуры. В качестве акцессорного минерала в них кроме хромшпине-лида появляется ильменит, который иногда находится с ним в тесном срастании. Породы массива в различной степени преобразованы низкотемпературными постмагматическими процессами, вплоть до образования в краевых частях хлорит-серпентиновых и актинолитовых сланцев.
Оруденение в массиве представлено в основном вкрапленными сульфидными рудами, которые состоят главным образом из пентландита, в переменных количествах присутствует халькопирит, трои-лит и кубанит. Максимальное количество сульфидов достигает 10 об. %. Минералы платиновой группы представлены главным образом сперрили-том и различными висмуто-теллуридами палладия. Первичные минеральные ассоциации были преобразованы в условиях низких уровней фуги-тивности кислорода и серы в системе [23]. При этом в крайних случаях происходило образование ава-руита и самородной меди и замещение сперрилита орселитом, а также Pt-Fe-Ni и Pt-Fe-Cu-сплавами.
Образцы и методы исследования
Поиск расплавных включений проводился в монофракциях акцессорных хромшпинелидов из дунитов и верлитов массива. Химический состав хромшпинелидов определялся в Центре коллективного пользования «Изотопно-геохимические исследования» Института геохимии им. А.П. Виноградова Сибирского отделения Российской академии наук (ИГХ СО РАН).
Полированные препараты просматривались под микроскопом, и отбирались зерна хромшпинели-дов, содержащие многофазные микровключения определённой формы, размера (10-30 мкм) и положения внутри зерна. Далее проводились высокотемпературные эксперименты по методике [14, 29]. Выбор температурного режима (1280-1300 °С) определялся задачей перевести включения в расплавленное состояние, затем его закалить и полученное стекло проанализировать на рентгеновском микроанализаторе и сканирующем микроскопе.
Все экспериментальные исследования включений проводились в лаборатории геодинамики и магматизма Института геологии и минералогии им. В.С. Соболева СО РАН (ИГМ СО РАН, г. Новосибирск). Включения были изучены на электронном сканирующем микроскопе (СЭМ) LEO 1430 VP (ИГМ СО РАН). Пределы обнаружения компонентов: 0,1-0,25 мас. %. Составы включений и хромшпинелидов анализировались на рентгеновском микроанализаторе Camebax-Micro (ИГМ СО РАН). Пределы обнаружения (мас. %) компонентов этим методом следующие: SiO2 - 0,007, TiO2 -0,032, Al2O3 - 0,011, Cr2O3 - 0,019, FeO - 0,019, MnO - 0,034, MgO - 0,011, CaO - 0,008, Na2O - 0,017, K2O - 0,009, Cl - 0,017, P2O5 - 0,011.
Для оценки РТ параметров кристаллизации ультраосновных пород массива Медек были использованы геотермобарометры, основанные на составах клинонопироксенов [30-34] и амфиболов [35-40].
Определение физико-химических параметров формирования массива проводилось с помощью обработки оригинальной информации по составам расплавных включений программными комплексами COMAGMAT [24-26] и PETROLOG [27].
Составы минералов из верлитов массива Медек
Оливины. В слабо измененных верлитах оливин образует идиоморфные и гипидиоморфные зерна, размером 2-3 мм, в серпентинизированных он сохранился только в виде реликтов. Изучение состава оливина по разрезу массива показало, что содержание форстеритового компонента варьирует от 88-90 до 73 %, при этом содержания NiO уменьшаются от 0,3-0,4 до 0,1 мас. %, а MnO увеличиваются от 0,1 до 0,4 мас. % (табл. 1). Также сверху вниз по разрезу наблюдается уменьшение объемного содержания оливина с 98-96 до 65 %. Кроме того, вместо акцессорного хромшпинелида в виде включений в оливине появляется ильменит с низким содержанием MgO.
Таблица 1. Представительные анализы оливинов из верлитов массива Медек
Table 1. Representative analyses of olivines from wehrlite of the Medek massif
Образец Sample SiO2 FeO MnO MgO NiO Сумма Total Fo Глубина Depth
14-85 39,30 12,70 0,20 47,50 0,20 99,90 87,00 386
14-92 40,60 14,30 0,30 44,20 0,30 99,70 84,70 429
14-99 40,30 13,80 0,30 45,70 0,30 100,40 85,50 456
14-100 40,00 14,10 0,30 44,40 0,30 99,10 84,90 457
14-102 40,20 15,20 0,20 45,10 0,30 101,00 84,10 471,4
14-103 39,10 18,10 0,40 42,20 0,20 100,00 80,60 483
14-112 38,50 22,70 0,30 37,70 0,10 99,30 74,80 517
Примечание. Глубина - глубина отбора пробы из керна скважины, м; Fo - содержание форстеритового минала в оливине.
Note. Depth - the depth of sampling from the well core, m; Fo - the content of the forsterite component.
Клинопироксены чаще реликтовые, по составу отвечают сУбкальциевому авгиту (Wo42En52Fs3) и диопсиду (Wo49En48Fs3), с содержаниями Al2O3 -0,1-1,8 мас. %, Cr2O3 - до 0,7 мас. % (табл. 2).
Хромшпинелиды (в которых были найдены расплавные включения). Как показали проведенные нами ранее исследования ультраосновных комплексов [14], составы хромшпинелидов, содержащих расплавные включения и соответственно имеющих магматогенное происхождение, могут служить прямыми источниками сведений для петрологических построений. Поэтому особое внимание уделялось также и составу минералов-хозяи-нов. Проведенная оценка возможного влияния расплавных включений на состав хромшпинели-дов продемонстрировала его полное отсутствие (табл. 3).
Таблица 2. Представительные анализы клинопироксенов из верлитов массива Медек
Table 2. Representative analyses of clinopyroxenes from wehrlite of the Medek massif
Образец Sample SiO2 AI2O3 Cr2Oa FeO MgO CaO Na2O Сумма Total En
14-95 55,21 0,19 н.о. 2,02 17,09 25,28 н.о. 99,79 46,88
14-96 55,26 0,17 н.о. 2,03 17,41 25,43 н.о. 100,30 47,20
14-97 55,91 0,14 н.о. 1,93 17,29 25,37 н.о. 100,64 47,16
14-93 54,95 н.о. н.о. 1,62 17,09 25,51 н.о. 99,17 46,96
14-93 55,15 0,11 н.о. 1,75 17,61 25,51 0,11 100,24 47,41
14-93 54,83 0,20 0,11 1,82 17,44 25,12 0,13 99,64 47,47
14-93 54,94 0,16 н.о. 1,75 17,59 25,35 0,11 100,01 47,54
14-93 54,63 0,10 н.о. 1,83 17,50 25,26 0,13 99,44 47,42
14-99 55,09 н.о. н.о. 1,69 17,28 25,34 0,12 99,52 47,14
14-99 54,62 0,11 н.о. 1,82 17,25 25,29 н.о. 99,09 47,08
Примечание. En - содержание энстатитового компонента; н.о. -не определено.
Note. En - content of the enstatite component; н.о. - not determined.
Cr#
Рис. 2. Диаграмма Cr#-Mg# для хромшпинелидов из верлитов массива Медек (Восточный Саян). 1-7 - хромшпинелиды: из верлитов массива Медек (1); из дунитов Нижнетагильского массива, Урал (2); из дунитов Кондерского (3), Инаглинского (4) массивов; из дунитов Карашат-ского массива (5, офиолиты Южной Тувы); из троктоли-тов района Сьерра-Леоне (6) и из перидотитов зоны разлома 15'20' (7) в Центральной Атлантике. Mg#=Mg*100/(Mg+Fe2+). Cr#=Cr*100/(Cr+Al). Рисунок составлен на основе оригинальных данных с использованием материалов из работ [14,29,41]
Fig. 2. Diagram Cr#-Mg# for Cr-spinel from wehrlite of the Medek massif (East Sayan). 1-7 - Cr-spinel: from wehrite of the Medek massif (1); from dunites of the Nizhny Tagil massif, Ural (2); from dunites of the Konder (3), Inaglinsky (4) massifs; from dunites of the Karachat massif (5, ophiolites of the South Tuva); troctolites area Sierra Leone (6) and pe-ridotites from the Fracture Zone 15'20' (7) in the Central Atlantic. Mg#=Mg*100/(Mg+Fe2+). Cr#=Cr*100/(Cr+Al). The figure is based on the original data using materials from the works [14,29,41]
Проведенные исследования показали умеренную хромистость (Cr# до 50,7) рассмотренных хромшпинелидов, и по этому показателю они существенно отличаются от гораздо более хромистых минералов из ультрабазитов платиноносных массивов Сибирского кратона (Кондер и Инагли) и Урала
Таблица 3. Представительные анализы хромшпинелидов из верлитов массива Медек Table 3. Representative analyses of Cr-spinel from wehrlite of the Medek massif
Образец/Sample SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O NiO Сумма/Total Cr#
2-2A/2 0,06 0,16 27,93 38,74 22,89 0,16 10,67 0,02 0,05 0,03 100,70 48,20
2-2A/3 0,02 0,18 27,55 38,27 22,78 0,17 10,70 0,01 0,04 0,06 99,78 48,24
2-2A/4 0,02 0,24 27,70 38,44 22,66 0,18 10,67 0,02 0,04 0,06 100,02 48,21
2-1-3/1 0,05 0,27 25,88 39,66 22,95 0,20 10,95 0,10 0,05 0,08 100,17 50,69
2-1-3/2 0,04 0,32 26,98 38,43 22,69 0,14 11,17 0,03 0,06 0,07 99,92 48,86
2-1-3/3 0,03 0,22 27,36 38,78 22,78 0,13 11,06 0,02 0,04 0,08 100,52 48,74
2-13/4 0,02 0,26 27,33 38,58 22,90 0,19 11,09 0,03 0,04 0,08 100,54 48,64
2-1-2/1 0,05 0,33 26,21 39,58 22,95 0,15 10,64 0,03 0,03 0,08 100,04 50,32
2-1-2/2 0,01 0,23 26,83 38,83 23,18 0,15 10,99 0,02 0,03 0,08 100,36 49,26
2-1-2/3 0,03 0,23 26,83 38,99 22,94 0,19 10,93 0,02 0,04 0,04 100,22 49,36
2-1-2/4 0,05 0,35 26,30 38,87 23,59 0,20 10,59 0,08 0,03 0,13 100,18 49,79
Примечание/Мае. Ст#=Сг*100/(Сг+А1)
(Нижнетагильский массив). При этом они хорошо согласуются с данными по хромитам из базит-ги-пербазитовых комплексов современных и древних (офиолиты) океанов (рис. 2). По соотношению хрома и железа изученные хромиты из верлита массива Медек наиболее соответствуют минералам из ульт-рабазитов Срединно-Атлантического хребта (САХ) в районе Сьерра-Леоне. По соотношению ТЮ2 -М§# они наиболее близки к хромитам из дунитов Карашатского массива (офиолиты Южной Тувы). Необходимо отметить, что хромшпинелиды из рассмотренных базит-ультрабазитовых пород современных (Срединно-Атлантический хребет) и древних (Карашатский массив) океанических структур содержат расплавные включения и, соответственно, имеют магматогенное происхождение [29, 41].
Расплавные включения в хромшпинелидах
из верлитов массива Медек
Было установлено, что хромшпинелиды из всех изученных проб содержат микровключения, которым присущи некоторые признаки, характерные для расплавных включений. При тщательном просмотре выяснилось, что «включения» в хромитах из дунитов и части верлитов обладают очень невыдержанными формами и разнообразным содержимым, без характерных пузырьков газа, что не позволяет их использовать для дальнейшего физико-химического моделирования. В единственном образце, в хромшпинелидах из верлита, найденные силикатные включения по внешним признакам и химическому составу соответствовали включениям расплавов.
Исследования на СЭМ показали, что прогретые силикатные включения содержат флюидный пузырек, располагающийся в гомогенной среде и имеющий почти всегда фактически эталонную круглую форму (рис. 3, А, Б). Такого рода равновесные формы флюидных включений являются характерным признаком расплавных, которые при нагреве и быстром охлаждении превратились в стекло. Соотношение объемов флюидных пузырьков и включений свидетельствует о существенном содержании летучих компонентов в магматических системах.
В зависимости от количества внутренних фаз включения подразделяются на две группы: 1) со стеклом и флюидным пузырьком (рис. 3, А, Б) и
2) со стеклом, пузырьком и кристаллическими фазами (рис. 3, В, Г). Кристаллические фазы представляют собой реликты высокомагнезиального тугоплавкого оливина, вокруг которого в одном случае (на пример № 4 на рис. 3, В) образуется кайма, а в другом (№ 3 на рис. 3, Г) он обрастает со всех сторон новообразованными микрокристаллическими фазами. В обоих случаях светлые микрокристаллики образуют фактически структуру спинифекс (рис. 3, В, Г), сформировавшуюся практически мгновенно во время закалки микроконтейнера с пробой в воде. Кроме игольчататой, новообразованные кристаллы иногда приобретают изоме-тричную форму (№ 4 на рис. 3, Г), и по данным количественного анализа (табл. 4) находятся ближе всего к высоко титанистому амфиболу - керсутиту (субкальциевый субкремнистый ферро-керсутит) [40]. В то же время недоплавленные оливины во включениях имеют высокие значения форстерито-вого компонента, сравнимые с данными по оливинам из дунитов Медекского массива.
Таблица 4. Представительные анализы кристаллических фаз в прогретых включениях в хромшпинелидах из верлитов массива Медек
Table 4. Representative analyses of crystal phases in heated inclusions in the Cr-spinel from wehrlite of the Medek massif
Образец/Sample 1 2 3 4 5 6
SiO2 41,25 40,74 40,45 40,43 32,32 30,57
TiO2 н.о. н.о. н.о. н.о. 18,6 19,51
Al2O3 н.о. н.о. н.о. н.о. 8,23 7,43
ClA 1,01 1,02 1,02 1,23 1,4 1,17
FeO 7,87 8,99 10,45 10,09 16,49 17,99
MnO н.о. н.о. н.о. н.о. н.о. 0,33
MgO 49,46 48,79 47,46 47,99 7,91 11,06
CaO н.о. н.о. 0,19 0,26 10,84 8,59
Na2O н.о. н.о. н.о. н.о. 1,24 1,33
K2O н.о. н.о. н.о. н.о. 0,15 н.о.
NiO 0,42 0,47 0,42 н.о. н.о. н.о.
Сумма/Total 100,01 100,01 99,99 100 97,18 97,98
Fo 91,42 90,2 88,4 89,14 - -
Примечание. 1-4 - оливин; 5,6 - кристаллические фазы, близкие по химическому составу к амфиболам (керсутитам); н.о. - не определено.
Note. 1-4 are the olivine; 5,6 are the crystalline phases similar in chemical composition of the amphiboles (kaersutite); н.о. - not determined.
Рис. 3. Первичные силикатные микровключения в хромшпинелидах из верлитов массива Медек (после высокотемпературных опытов и закалки): А - первичное двухфазное расплавное включение первого типа после высокотемпературного эксперимента и закалки; Б - двухфазное расплавное включение второго типа после высокотемпературного эксперимента и закалки; В, Г - рас-плавные включения с недоплавленными и новообразованными кристаллическими фазами после высокотемпературных экспериментов и закалки. Сканирующий микроскоп. Отмечены участки анализов. Фл - флюидный пузырек; Ст - гомогенное силикатное стекло; Ол - недоплавленные реликты оливина; Ам - новообразованные микрокристаллиты амфибола
Fig. 3. Primary silicate micro-inclusions in Cr-spinel from wehrlite of the Medek massif (after high-temperature experiment and quenching): A - primary two-phase melt inclusion of the first type after high-temperature experiment and quenching; B - two-phase melt inclusion of the second type after high-temperature experiment and quenching; B, G - melt inclusions with under-melted and newly formed crystalline phases after high-temperature experiments and quenching. Scanning microscope. Marked areas of analysis. Фл - fluid bubble; Ст -homogeneous silicate glass; Ол - under-melted relicts of olivine; Ам - newly formed amphibole microcrystallites
Описанные выше особенности соотношений различных фаз, заполняющих расплавные включения, сформировались в ходе высокотемпературных экспериментов, в результате которых в зависимости от продолжительности происходило полное расплавление содержимого или (при недостатке времени, а также, возможно, и температуры) оставались реликты оливина, обраставшие при закалке новыми фазами.
В одном случае содержимое включений было полностью расплавлено и при закалке образовались только две фазы: стекло (закаленный силикатный расплав) и флюидный пузырек. В других случаях времени (и возможно температуры) экспериментов, по всей вероятности, не хватило для полного расплавления содержимого включений и внутри остались недоплавленные кристаллики оливина, на которые в ходе закалки нарастали новообразованные кристаллические фазы. Это позволяет сделать вывод, что для реконструкции физико-химических параметров реальной магматической системы необходимо использовать результаты изучения стекол двухфазовых включений.
Составы расплавов, участвовавших в формировании верлитов массива Медек
Количественный анализ на сканирующем электронном микроскопе и на микрозонде стекол из прогретых и закаленных двухфазных расплавных включений (стекло + круглый флюидный пузырек) в хромшпинелидах показал наличие двух типов силикатных расплавов (табл. 5). Расплавные включения первого типа обладают нормальной щелочностью, и по соотношению суммы №20+К20 (до 1.6 мас. %) и содержаний 8Ю2 относятся к толеитовой серии, что характерно также для включений в хромшпинелидах из базит-ультрабази-товых комплексов САХ и офиолитов Южной Тувы [29, 41]. Однако, сравнивая их состав с уровнем содержаний щелочей во включениях из хромитов массива Кондёр [14, 21], можно отметить, что последним присущи более высокие величины содержаний К20.
По уровню содержаний М§0 (10,5-13,4 мас. %) и 8Ю2 (44,8-46,5 мас. %) включения этого типа располагаются на границе полей оливиновых базальтов и пикробазальтов (рис. 4), наиболее близки к ним по соотношению этих элементов сте-
Таблица 5. Представительные анализы стекол прогретых включений в хромшпинелидах из верлитов массива Медек Table 5. Representative analyses of glasses of heated inclusions in Cr-spinel from wehrlite of the Medek massif
№ Образец/Sample SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O K2O P2O5 Сумма/Total
1 2A-2 46,14 0,87 13,36 1,18 9,47 0,11 13,36 14,96 1,18 0,27 0,06 100,96
2 2A-2/1 45,79 0,91 13,74 1,23 9,59 0,15 10,68 16,26 1,30 0,25 0,05 99,95
3 2A2/2 45,90 0,98 13,55 1,38 9,44 0,11 11,41 15,86 1,15 0,27 0,03 100,07
4 2A-2/4 45,87 0,96 13,37 1,26 9,47 0,17 11,08 15,99 1,27 0,27 0,06 99,77
5 2A-2/5 46,42 0,89 13,62 1,37 9,16 0,18 11,86 15,64 1,16 0,25 0,05 100,60
6 2-2A/7 45,58 0,88 13,65 1,18 9,00 0,12 11,19 15,39 1,27 0,28 0,07 98,60
7 2-2A/8 45,80 0,85 13,23 1,39 9,28 0,09 11,05 15,43 1,21 0,26 0,06 98,66
8 2-2A/9 46,54 0,96 13,75 1,22 9,07 0,10 10,45 15,86 1,19 0,27 0,07 99,48
9 2-2a-2/1 45,63 0,82 13,68 1,28 10,44 н.о. 10,54 16,06 1,26 0,30 н.о. 100,01
10 2-2a-2/2 44,79 0,85 13,02 1,82 10,41 н.о. 11,46 16,08 1,23 0,35 н.о. 100,01
11 2-1-2/1 32,74 3,86 6,69 1,73 21,16 0,39 9,44 19,68 0,36 0,02 2,24 98,31
12 2-1-2/2 32,43 3,78 6,75 1,98 20,87 0,36 9,29 19,87 0,34 0,02 2,27 97,96
13 2.-1-2/1 32,79 3,71 6,61 1,79 23,24 0,45 9,14 20,08 н.о. н.о. 2,20 100,01
14 2.-1-2/2 32,47 3,48 6,86 1,83 23,49 0,43 9,07 20,17 н.о. н.о. 2,21 100,01
15 2-1-3/1 31,95 4,01 6,67 1,97 21,67 0,36 9,38 19,55 0,37 0,03 2,17 98,13
16 2-3-1/2 32,60 3,87 6,85 1,69 21,57 0,37 9,32 19,42 0,35 0,02 2,17 98,23
17 2.-1-3/1 32,76 3,70 6,75 1,50 23,37 0,47 8,88 19,86 0,43 н.о. 2,29 100,01
18 2.-1-3/2 32,72 3,71 6,87 1,68 23,37 0,44 8,91 20,00 0,22 н.о. 2,09 100,01
Примечание. 1-10 - включения первого типа; 11-18 - включения второго типа; н.о. - не определено. Note. 1-10 are the inclusions of the first type; 11-18 are the inclusions of the second type; н.о. - not determined.
кла из расплавных включений Нижнетагильского и, в меньшей степени, Карашатского массивов.
15 -
о го :>
О го О
• Med
□ Коп
Д Cha
ONT
■ Kar
42 44 46 48
Si02, мас.%
Рис. 4. Диаграмма MgO-SiO2 для стекол двухфазных (после высокотемпературных экспериментов) включений первого типа в хромшпинелидах из верлитов массива Медек. Стекла расплавных включений в хромшпинелидах: из верли тов массива Медек (Med); из дунитов Кондерского (Kon) Чадского (Cha), Нижнетагильского (NT) и Карашатско го (Kar) массивов. Поля по [42]: 1 - пикриты, 2 - пикроба зальты, 3 - оливиновые базальты, 4 - базальты, 5 - сте кла расплавных включений в хромшпинелидах из трокто-литов района Сьерра-Леоне (Срединно-Атлантический хребет). Рисунок составлен на основе оригинальных данных с использованием материалов из работ [14,29,41]
Fig. 4. MgO-SiO2 diagram for glass of two-phase (after high-temperature experiments) inclusions of the first type in Cr-spinel from wehrlite of the Medek massif. Glass of melt inclusions in Cr-spinel: from wehrlite of the Medic massif (Med); from du-nite Conder (Kon), Chad (Cha), Nizhny Tagil (NT) and Ka-rashat (Kar) massifs. Fields: 1 - picrites, 2 - picrobasalts, 3 - olivine basalts, 4 - basalts, 5 - glass of melt inclusions in chrome spinel from troctolite of Sierra Leone (Mid-Atlantic ridge). The figure is based on the original data using materials from the works [14,29,41]
0 5 10 15 20 25 30 35
MgO, мас.%
Рис. 5. Соотношение CaO-MgO в составах стекол двухфазных (после высокотемпературных экспериментов) включений первого типа в хромшпинелидах из верлитов массива Медек
Fig. 5. CaO-MgO ratio in compositions of glass of two-phase (after high-temperature experiments) inclusions of the first type in Cr-spinel from wehrlite of the Medek massif
С помощью диаграмм, показывающих зависимость значений химических компонентов в расплавных включениях в хромитах от MgO, можно выяснить особенности эволюции магматических систем при кристаллизации расслоенных базит-ульт-рабазитовых плутонов. Подобная ситуация была нами рассмотрена на примере Карашатского базит-ультрабазитового массива [41]. Составы включений в хромитах целого ряда массивов трассируют путь последовательной внутрикамерной кристаллизации расплавов (рис. 5), что характерно и для расплавных включений этого типа в массиве Медек. При этом точки их составов находятся в тесной ассоциации с данными по породам верлит-клинопирок-сенит-габбровой серии Карашатского массива. Та-
FeO/MgO
• Med □ Коп ACha О NT ■ Kar
Рис. 6. Соотношение Al2O3-FeO/MgO в составах стекол двухфазных (после высокотемпературных экспериментов) включений первого типа в хромшпинелидах из верлитов массива Медек. Тренды: Ol - оливиновые кумуляты, Pl -плагиоклазовые кумуляты, RM - остаточные расплавы. 1 - поле стекол расплавных включений в хромшпинелидах из троктолитов района Сьерра-Леоне (Срединно-Атлан-тический хребет). Остальные условные обозначения см. на рис. 4. Рисунок составлен на основе оригинальных данных с использованием материалов из работ [14,29,41]
Fig. 6. Al2O3-FeO/MgO ratio in compositions of glass of two-phase (after high-temperature experiments) inclusions of the first type in chrome spinel from wehrlite of the Medek massif. Trends: Ol - olivine cumulates, Pl -plagioclase cumulate, RM - residual melts. 1 - field of glass of melt inclusions in chrome spinel from troctolite of Sierra Leone (Mid-Atlantic ridge). Other symbols are shown in Fig. 4. The figure is based on the original data using materials from the works [14,29,41]
Рис. 7. Вариационные диаграммы (мас. %) для стекол двухфазных (после высокотемпературных экспериментов) включений первого типа в хромшпинелидах из верлитов массива Медек. 1 - поле стекол расплавных включений в хромшпинелидах из троктоли-тов района Сьерра-Леоне (Срединно-Атлантический хребет). Остальные условные обозначения см. на рис. 4. Рисунок составлен на основе оригинальных данных с использованием материалов из работ [14,29,41]
Fig. 7. Variation diagrams (wt. %o) for glasses of two-phase (after high-temperature experiments) inclusions of the first type in Cr-spinel from wehrlite of the Medek massif. 1 - field of glass of melt inclusions in Cr-spinel from troctolite of Sierra Leone (Mid-Atlantic ridge). Other symbols are shown in Fig. 4. The figure is based on the original data using materials from the works [14,29,41]
кая ^ке картина наблюдается и на диаграмме А1203 и Ре0/М£0, где фигуративные точки составов включений первого типа лежат недалеко от линии тренда, характерного для оливиновых кумулатов, и тесно ассоциируют с полями составов включений из троктолитов района Сьерра-Леоне (рис. 6).
5 - i
4 - 1
О СО :> з- \ 1
О F 2 - \
1 -п •41 Д + виаЫ. d»
• Med OMed2 □ Коп A Cha О NT ■ Kar
MgO, мас.%
Рис. 8. Содержания титана в стеклах двухфазных (после высокотемпературных экспериментов) включений второго типа в хромшпинелидах из верлитов массива Медек. Стекла расплавных включений в хромшпинелидах: из верлитов массива Медек (первый тип - Med, второй тип - Med2); из дунитов Кондерского (Kon), Чадского (Cha), Нижнетагильского (NT) и Карашатского (Kar) массивов. 1 - поле стекол расплавных включений в хромшпинелидах из троктолитов района Сьерра-Леоне (Срединно-Атлантический хребет). Рисунок составлен на основе оригинальных данных с использованием материалов из работ [14,29,41]
Fig. 8. Content of titanium in glasses of two-phase (after high-temperature experiments) inclusions of the second type in Cr-spi-nel from wehrlite of the Medek massif. Glass of melt inclusions in Cr-spinel from wehrlite of the Medek massif (first type - Med, the second type - Med2); from dunites of the Kon-der (Kon), Chad (Cha), Nizhny Tagil (NT) and Karachat (Kar) massifes. 1 - the field of glass of melt inclusions in Cr-spinel from troctolite of Sierra Leone (Mid-Atlantic ridge). The figure is based on the original data using materials from the works [14,29,41]
На других вариационных диаграммах (рис. 7) стекла включений первого типа наиболее близки по основным петрохимическим компонентам (SiO2, Al2O3, FeO, Na2O, K2O) к данным по расплав-ным включениям в хромшпинелидах из пород океанических структур - Карашатский массив (офио-литы Южной Тувы) и район Сьерра-Леоне (Сре-динно-Атлантический хребет). Особенно это хорошо видно по соотношению K2O-SiO2, когда рас-плавные включения в хрошпинелидах массива Медек, офиолитов Южной Тувы и Срединно-Ат-лантического хребта образуют единый низкокалиевый (менее 1 %) тренд, в отличие от включений из платиноносных массивов, обладающих явно повышенными содержаниями К20 (рис. 7). Эти данные позволяют сделать вывод о том, что образование верлитов происходило при участии низкокалиевых пикробазальтовых расплавов нормальной
щелочности, наиболее близких по петрохимическим характеристикам к данным по расплавным включениям в хромитах из океанических комплексов.
В то же время наличие многофазных включений с недоплавленными кристалликами оливина свидетельствует в пользу предположения, что исходная магма для пород массива была более магнезиальной, чем пикробазальтовый расплав. Исходя из этого нами был рассчитан предполагаемый состав первичного расплава для пород массива Медек [22], который по уровню содержаний БЮ2 (45 мас. %) и М£0 (18 мас. %) отвечает пикритовому.
Расплавные включения второго типа. Фазовый состав и морфология (рис. 3, Б) расплавных включений второго типа идентичны таковым первого. Главные отличия заключаются в особенностях химического состава, а именно: высоких содержаниях ТЮ2 (макс. 3,7 мас. %), ЕеО (макс.
23.5 мас. %), Р205 (макс. 2,3 мас. %), при этом средние содержания БЮ2 составляют всего
32.6 мас. %. Полные аналоги магматических пород и расплавных включений такого состава нам не известны. Кроме того, наблюдаемые соотношения компонентов во включениях этого типа вряд ли могут соответствовать составу расплава, из которого осаждались изученные хромшпинелиды. Однако возможность существования таких, по сути, рудоносных расплавов подтверждается результатами изучения Ее-Т1 интрузивных комплексов Центральной Атлантики [29] (рис. 8).
Эволюция физико-химических условий
кристаллизации расплавов
при формировании массива Медек
Основой для выяснения особенностей эволюции физико-химических условий кристаллизации магматических систем массива Медек являются исследования расплавных включений в хромшпи-нелидах из верлитов. В то же время к решению генетических задач привлекались и другие методы, включая минералогические термометры и барометры, а также расчетное моделирование, так как только согласованные результаты различных методик могут послужить наиболее обоснованными выводами об условиях кристаллизации пород.
Прежде всего, наличие первичных силикатных расплавных включений в хромшпинелидах (и преобразование их содержимого после высокотемпературных опытов в равновесную систему гомогенное стекло + круглый флюидный пузырек) прямо свидетельствует о том, что рассмотренные верлиты массива Медек формировались из расплавов. Данные по стеклам прогретых включений (с содержанием магния до 13,4 мас. %), близкие к пикроба-зальтовым магмам, позволили рассчитать возможный состав исходного пикритоидного расплава (с М£0 - 17,85 мас. %) для массива Медек [22]. Именно эти составы расплавов послужили в дальнейшем основой для расчетов условий кристаллизации рассмотренных верлитов.
Как было отмечено выше, основой для решения генетических задач послужили данные по расплав-ным включениям. При этом дополнительную и весьма важную информацию удалось получить с помощью минералогических термометров и барометров. В частности, на основе составов клинопирок-сенов были рассчитаны с использованием минералогических термометров и барометров [31-33] РТ параметры кристаллизации верлитов массива Ме-дек: около 1255 °С при давлении 2,3-5,9 кбар. Очень полезную информацию несут находящиеся в расплавных включениях в хромшпинелидах микрокристаллики амфиболов (рис. 3, В, Г), опираясь на составы которых, по минералогическому барометру [36], удалось оценить давление: 3,2-4,2 кбар.
Температуры кристаллизации верлитов были получены в результате численного моделирования в программном комплексе PETROLOG [27]. Основой для расчетов послужили составы стёкол включений первого типа, граничные условия соответствовали давлению 4 кбар (установлено по усредненным оценкам по пироксеновым и амфиболо-вым геобарометрам), Н20 - 0,5 %, расчетная величина хромистости хромшпинелида корректировалась по реальным значениям и составила 48,2 %. Исходя из этого, было установлено, что температура кристаллизации оливинов снижалась от 1275 до 1240 °С, хромшпинелиды формировались в диапазоне 1260-1240 °С, а клинопироксены -1235-1220 °С, данные по которым хорошо согласуются с оценками по геотермометру.
Расчеты условий кристаллизации расплавов по программе COMAGMAт [24, 25] проведены на основе данных по составам расплавных включений первого типа в хромшпинелидах из верлита и с учетом результатов использования минералогических термометров и барометров.
Равновесная кристаллизация возможного исходного расплава (с М§0 - 17,85 мас. %) в ходе декомпрессии от 10 кбар (низы коры) до 2 кбар. Буфер QFM. Содержание воды - 0,5 мас. %. Расчеты показали широкие диапазоны возможных давлений кристаллизации минералов из верлитов массива Медек. В то же время данные по реальным составам оливинов ^о=80,6-85,5) и клинопироксе-нов (Еп=46,8-47,6) верлитов этого массива в ходе анализа результатов расчетного моделирования позволили уточнить давления при кристаллизации оливинов (2,6-4,9 кбар) и клинопироксенов (2,9-3,6 кбар). Обращают на себя внимание перекрывающиеся в значительной мере диапазоны давлений при формировании этих двух минералов, а также сходство с оценками по пироксеновым и ам-фиболовым барометрам.
Расчеты фракционной кристаллизации того же исходного расплава (с М§0 - 17,85 мас. %) при тех же условиях (в ходе декомпрессии от 10 кбар до 2 кбар; буфер QFM•, Н2О - 0,5 мас. %) показали реальные максимальные давления (с учетом составов минералов) для оливина - до 4,9 кбар, и для клинопироксена - до 4,2 кбар, что подтверждает
результаты равновесной кристаллизации и минералогических барометров.
Таким образом, проведенные расчеты по различным минералогическим барометрам и по различным моделям кристаллизации на основе составов расплавных включений показали, что давление при формировании верлитов массива Медек было около 4 кбар (в среднем 2,8-4,7 кбар).
Термометрия равновесий минерал-расплав. Изобарическая система - 4 кбар (согласно оценкам по минералогическим барометрам и равновесной декомпрессионной кристаллизации). Буфер QFM. Состав магмы - стекло расплавных включений пи-кробазальтового состава ^§0 - 13,36 мас. %) и рассчитанный состав исходного расплава ^§0 -17,85 мас. %). Содержание воды - 0,5 мас. %. В результате для исходного пикритоидного расплава определены следующие температуры равновесия: оливин - 1380 °С, клинопироксен - 1270 °С. В случае пикробазальтового расплава: оливин -1300 °С, клинопироксен - 1240 °С. Обращает на себя внимание соответствие рассчитанных в случае пикробазальтового расплава (данные по составу расплавных включений) температур кристаллизации клинопироксена с данными по пироксеновому термометру и по результатам расчетов с помощью программы PETR0L0G [27].
хххххххххххххххххх х х х::::::: :::::><
s"
о
oooooooooSSAAMtei?®®
10 -
,оО
оО^
о -
Рис. 9.
+44Ф+++++++++ +++ + + + + +
OOI
■ Срх XSi02 О AI203 ДМдО + Na4K
1380 1330 1280 1230 1180 Т°С
Фракционная кристаллизация исходного расплава (MgO -17,85 мас. %) для массива Медек. Компонент, % -изменения содержаний минералов (Ol - оливин, Cpx -клинопироксен) и химических компонентов (Na+K=Na2O+K2O, мас. %) в расплаве при снижении температуры. Расчетное моделирование с помощью программы COMAGMAT [25,26]
Fig. 9. Fractional crystallization of the initial melt (MgO -17,85 wt. %) for the Medek massif. Component, % - changes in mineral content (Ol - olivine, Cpx - clinopyroxene) and chemical components (Na+K=Na2O+K2O, wt. %) in the melt with decreasing temperature. Computational modeling using COMAGMAT [25, 26]
Расчеты фракционной кристаллизации исходного расплава (MgO - 17,85 мас. %) при тех же условиях (в изобарической системе - 4 кбар. Буфер QFM. Содержание воды - 0,5 мас. %) показали, что оливин начинает кристаллизоваться при 1375 °C (совпадает с данными по равновесию -1380 °C). Процессы снижения температуры и роста
количества оливина происходят относительно быстро - до 1235 °C, когда начинает кристаллизоваться клинопироксен (рис. 9). Далее интенсивно растет клинопироксен и скорость его накопления превышает кристаллизацию оливина практически в четыре раза - идет формирование верлита. Таким образом, мы видим, что дуниты массива Медек формируются начиная с 1380 °C, а верлиты образуются при более низких температурах -1235 C.
Одновременно с кристаллизацией оливина и клинопироксена закономерно изменяется состав остаточного расплава: падает MgO, растет количество щелочей. Интересно поведение SiO2 - с началом кристаллизации клинопироксена рост этого компонента сменяется падением при снижении температуры. Фактически проявляется реверсивный тренд, хорошо выделяющийся на диаграмме MgO-SiO2 (рис. 10, А) и отмеченный нами при изучении расплавных включений в хромшпинелидах из пород других ультраосновных массивов [21]. Очень важным является факт расположения составов расплавных включений вдоль расчетных трендов изменения состава расплава (рис. 10, А, Б), то есть здесь экспериментальные данные по включениям с одной стороны подтверждают достоверность проведенных расчетов, а с другой - такое соотношение показывает, что включения не являются случайностью, а фиксируют закономерную и теоретически обоснованную эволюцию расплава.
На диаграмме CaO-MgO (рис. 10, Б) точки составов стекол расплавных включений приурочены в основном к изгибу тренда изменения состава расплава, который практически совпадает с магнезиальной ветвью тренда расплавных включений в хромшпинелидах из дунитов Карашатского массива (Южная Тува).
Расчеты фракционной кристаллизации пи-кробазальтового расплава (состав стекла расплав-ного включения с максимальным М§0 -13,36 мас. %) в изобарической системе (4 кбар. Буфер QFM. Содержание воды - 0,5 мас. %) показывают, что оливин начинает кристаллизоваться около 1275 °С, но количество его (до 12 %) существенно меньше, чем клинопироксена (до 28 %), который образуется при более низких температурах, начиная с 1230 °С (рис. 11). Обращает на себя внимание сходство этих параметров со всеми полученными другими методиками температурными характеристиками для клинопироксена. В случае оливина температуры согласуются с данными моделирования по программе PETROLOG [27].
Расчеты условий кристаллизации пород по разрезу массива Медек были проведены для исходного расплава (М§0 - 17,85 мас. %) при условиях изобарической системы - 4 кбар. Буфер QFM. Содержание воды - 0,5 мас. %. Наличие данных по реальным составам оливинов в образцах, отобранных по разрезу скважины [28], позволяет оценить достоверность результатов расчетного моделирования.
В настоящий момент расположение массива Медек перевернуто по отношению к его первичному залеганию [28]. В связи с этим на рис. 12 с данными по расчетному моделированию процессов кристаллизации в исходной камере метраж по скважине идет в обратном направлении - снизу вверх. Соответственно максимальным температурам (около 1380-1350 °С) кристаллизации оливина (на нижних уровнях камеры) соответствует его состав с максимальными значениями форстерито-вого компонента. Далее при понижении температуры до 1300-1270 °С (в средних горизонтах камеры, отметки скважины 250-350 м) расчетное со-
о
m
uli
о ueit
• Ulf
мас.%
MgO, мас.%
Рис. 10. Эволюция состава исходного расплава (MgO - 17,85 мае. % ) в ходе его фракционной кристаллизации. Melt - состав расплава; Inc - состав стекла расплавных включений первого типа в хромшпинелидах из верлитов массива Медек. Пунктиром показан тренд составов расплавных включений в хромшпинелидах из дунитов Карашатского массива, Южная Тува [41]. Расчетное моделирование с помощью программы COMAGMAT [25,26]
Fig. 10. Evolution of the initial melt composition (MgO - 17,85 wt. %) during its fractional crystallization. Melt is the melt composition; Inc is the composition of glass of melt inclusions of the first type in Cr-spinel from wehrlite of Medek massif. The dashed line shows the trend of compositions of melt inclusions in Cr-spinel from dunites Karachat array, southern Tuva [41]. Computational modeling using CO-MAGMAT [25, 26]
держание Fo практически совпадает с данными по реальным анализам оливинов.
х х х х х х >::»>::.......;■ -....................„.
........................................
Рис
Fig.
тической камеры, формируя оливиновые и в меньшей степени оливин-пироксеновые кумуляты. Плагиоклаз, благодаря своей меньшей плотности, неизбежно оказывался на верхних уровнях. Так как массив Медек сейчас находится в перевернутом положении, изначально верхняя его часть (с плагиоклазом), в отличие от соседних массивов, имела мало шансов сохраниться.
1 290 1 270 1 250 1 230 1 210 1 190 1 170 Т°С
11. Фракционная кристаллизация пикробазальтового расплава состава стекла прогретого расплавного включения первого типа (с максимальным MgO - 13,36 мас. %) в хромшпинелиде из верлитов массива Медек
11. Fractional crystallization of the picrobasalt melt of the glass composition of the heated melt inclusion of the first type (with a maximum MgO - 13,36 wt. %) in Cr-spinel from wehrlite of the Medek massif
В последующем (с 1240-1230 °C, около отметки 380-400 м) начинается интенсивная кристаллизация клинопироксена, преобладающая над образованием оливина, состав которого (как расчетный, так и реальный) существенно меняется с падением форстеритового компонента - формируются верли-ты на самых верхних уровнях камеры (рис. 12). Таким образом, расчеты показывают возможность фракционной кристаллизации оливина (начиная с температур 1380 °C) из исходного пикритоидного расплава (MgO - 17,85 мас. %), приводящей к заполнению оливиновыми кумулятами (формирование дунитов) основного объема (до 400 из 500 м по вертикали) реальной камеры и сменяющейся при температурах 1240-1230 °C интенсивной кристаллизацией клинопироксенов, приводящей к образованию в верхних горизонтах камеры пироксен-оливиновых кумулятов - верлитов. Необходимо подчеркнуть, что данная модель подтверждается фактическим материалом - на рис. 12 отчетливо видна (особенно в верхней половине) хорошая согласованность расчетных и реальных данных по составам оливина.
В целом следует отметить, что при моделировании по программе COMAGMAT [25, 26] фракционной кристаллизации как пикробазальтового, так и исходного пикритоидного расплавов массива Ме-дек, расчеты показывают образование плагиоклаза, начиная примерно с 1200 °C. Возможность формирования плагиоклаз содержащих пород из подобных расплавов подтверждается данными по соседним (родственным) массивам: Медвежий Лог, Тартайском и Малая Шита, в которых присутствует оливиновое габбро [28]. В случае массива Медек, кристаллизовавшиеся первыми (при наиболее высоких температурах) оливины (и далее клинопи-роксены) осаждались в нижних горизонтах магма-
Рис. 12. Изменение состава оливина по разрезу массива Медек. Fo1 - значения реального форстеритового компонента, измеренного в образцах, отобранных по разрезу скважины; Fo2 - расчетные значения форстеритового компонента; T "C - расчетные температуры кристаллизации минералов из исходного расплава (MgO - 17,85 мас. %); Cpx+Ol и Ol - уровни кристаллизации минеральных па-рагенезисов соответственно верлита и дунита массива Медек; Lm - глубины пробоотбора керна скважины, м. Расчетное моделирование с помощью программы COMAGMAT [25,26]. Рисунок составлен на основе оригинальных данных с использованием материалов из работы [28]
Fig. 12. Changes in olivine composition in the section of the Medek massif. Fo1 are the values of the real forsterite component measured in samples taken from the well section; Fo2 are the calculated values of the forsterite component; T "C are the calculated crystallization temperatures of minerals from the initial melt (MgO - 17,85 wt. %); Cpx+Ol and Ol are the levels of crystallization of mineral assemblages, respectively wehrlite and dunite of the massif Medek; Lm are the well core sampling depths, m. Computational modeling using COMAGMAT program [25, 26]. The figure is based on the original data using materials from the work [28]
Заключение
1. В результате анализа стекол гомогенизированных расплавных включений установлено, что исходная магма для верлитов массива Медек отвечала по составу низкокалиевому пикроба-зальтовому расплаву, наиболее близкому по пе-
трохимическим характеристикам к данным по расплавным включениям в хромитах из океанических комплексов.
2. Геобарометрические расчеты и моделирование фазовых равновесий по программе COMAG-MAT показали, что давление при формировании пород массива было около 4 кбар (2,8-4,7 кбар).
3. Совместное использование минералогических термометров и расчетов по программам PETRO-LOG и C0MAGMAT, основанных на результатах экспериментальных исследований, свидетельствует о формировании дунитов массива из пикритоидного ^§0 - 17,85 мас. %) расплава, начиная с кристаллизации оливина при 1380 °С. При этом формирование верлитов происходило преимущественно из пикробазальто-вого ^§0 - 13,36 мас. %) расплава с кристаллизацией оливинаначиная с 1280-1275 °С, но в гораздо меньших количествах, чем клинопи-роксена, интенсивное образование которого
шло при более низких температурах (начиная с 1235-1230 С). 4. Сравнительный анализ результатов расчетного моделирования составов расплавов и минералов с данными по расплавным включениям и по реальным оливинам и клинопироксенам свидетельствует об их согласованности. В частности, важным является факт расположения составов рас-плавных включений вдоль расчетных трендов изменения расплава (рис. 10, А, Б), то есть здесь экспериментальные данные по включениям с одной стороны подтверждают достоверность проведенных расчетов, а с другой - такое соотношение показывает, что включения не являются случайностью, а фиксируют закономерную и теоретически обоснованную эволюцию расплава.
Работа выполнена в рамках государственного задания ИГМ СО РАН (проект № 0330-2016-0014) и ИГХ СО РАН (проект № 0350-2016-0030), при поддержке Министерства науки и высшего образования Российской Федерации, договора № 14.Y26.31.0029.
СПИСОК ЛИТЕРАТУРЫ
1. Barnes S.J., Lightfoot P.C. Formation of magmatic nickel-sulfide ore deposits and processes affecting their copper and platinum-group element contents // Economic Geology 100th Anniversary Volume. - 2005. - P. 179-214.
2. Platinum group element. Chromium and vanadium deposits in mafic and ultramafic rocks / R.G. Cawthorn, S.J. Barnes,
C. Ballhaus, K.N. Malich // Economic Geology 100th Anniversary Volume. - 2005. - P. 215-249.
3. The mineral system approach applied to magmatic Ni-Cu-PGE sulphide deposits / S.J. Barnes, A.R. Cruden, N.T. Arndt, B.M. Saumur // Ore Geology Reviews. - 2016. - V. 76. -P. 296-316.
4. Barnes S.J., Robertson J.C. Time scales and length scales in magma flow pathways and the origin of magmatic Ni-Cu-PGE ore deposits // Geoscience Frontiers. - 2019. - V. 10. - № 1. -P. 77-87.
5. Платинометальные месторождения мира. I. Платинометаль-ные малосульфидные месторождения в ритмично расслоённых комплексах / Д.А. Додин, Н.М. Чернышев, Д.В. Полферов, Л.Л. Тарновецкий. - М.: АО Геоинформарк, 1994. - 279 с.
6. Naldrett A.J. World-class Ni-Cu-PGE deposits: key factors in their genesis // Mineralium Deposita. - 1999. - V. 34. -P. 227-240.
7. Додин Д.А., Чернышев Н.М., Яцкевич Б.А. Платинометальные месторождения России. - СПб.: Наука, 2000. - 755 с.
8. Naldrett A.J. From the mantle to the bank: the life of a Ni-Cu-(PGE) sulfide deposit // South African Journal of Geology. - 2010. - V. 113.1. - P. 1-32.
9. Соболев А.В. Включения расплавов в минералах как источник принципиальной петрологической информации // Петрология. - 1996. - Т. 4. - № 3. - С. 228-239.
10. Primitive basaltic melts included in podiform chromites from the Oman ophiolite / P. Schiano, R. Clocchiatti, J.-P. Lorand,
D. Massare, E. Deloule, M. Chaussidon // Earth and Planetary Science Letter. - 1997. - V. 146. - № 3-4. - P. 489-497.
11. Состав магматогенных флюидов. факторы их геохимической специализации и металлоносности / Л.С. Борисенко, Л.А. Боровиков, Л.М. Житова, Г.Г. Павлова // Геология и геофизика. - 2006. - Т. 47. - № 12. - С. 1308-1325.
12. Рябчиков И.Д., Когарко JI.H., Соловова И.П. Физико-химические условия магмообразовання в основании Сибирского плю-
ма по данным исследования расплавных микровключений в меймечитах и щелочных пикритах Маймеча-Котуйской провинции // Петрология. - 2009. - Т. 17. - № 3. - С. 311-322.
13. Условия кристаллизации дунитов Кондерского платиноносно-го щелочно-ультраосновного массива, Юго-Восток Алданского щита / В.А. Симонов, В.С. Приходько, С.В. Ковязин,
A.В. Тарнавский // Тихоокеанская геология. - 2010. - Т. 29. -№5. - С. 44-55.
14. Симонов В.А., Приходько В.С., Ковязин С.В. Условия формирования платиноносных ультраосновных массивов Юго-Востока Сибирской платформы // Петрология. - 2011. - Т. 19. -№6. - С. 579-598.
15. Васильев Ю.Р., Гора М.П. Меймечит-пикритовые ассоциации Сибири, Приморья и Камчатки (сравнительный анализ, вопросы петрогенезиса) // Геология и геофизика. - 2014. - Т. 55. -№ 8. - С. 1211-1225.
16. Петрогенезис дунитов Гулинского ультраосновного массива (север Сибирской платформы) / В.А. Симонов, Ю.Р. Васильев, С.И. Ступаков, А.В. Котляров, Н.С. Карманов // Геология и геофизика. - 2016. - Т. 57. - № 12. - С. 2153-2177.
17. Chromian spinels from the Magura Unit (Western Carpathians. Eastern Slovakia) - their petrogenetic and palaeogeographic implications / K. Bonova, J. Spisiak, J. Bona, M. Kovacik // Geological Quarterly. - 2017. - V. 61. - № 1. - P. 3-18.
18. Renna M.R., Tribuzio R., Ottolini L. New perspectives on the origin of olivine-rich troctolites and associated harrisites from the Ligurian ophiolites (Italy) // Journal of the Geological Society. -2016. - V. 173. - №6. - P. 916-932.
19. Среднее содержание летучих компонентов. петрогенных и редких элементов в магматических расплавах главных геодинамических обстановок земли. I. Расплавы основного состава /
B.Б. Наумов, В.А. Дорофеева, А.В. Гирнис, В.В. Ярмолюк // Геохимия. - 2017. - № 7. - С. 618-643.
20. Миронов Н.Л., Портнягин М.В. Связь окислительно-восстановительных условий плавления мантии и содержаний меди и серы в первичных магмах на примере Толбачинского дола (Камчатка) и хребта Хуан де Фука (Тихий океан) // Петрология. -2018. - Т. 26. - № 2. - С. 140-162.
21. Физико-химические условия кристаллизации пород ультраосновных массивов Сибирской платформы / В.А. Симонов, В.С. Приходько, Ю.Р. Васильев, А.В. Котляров // Тихоокеанская геология. - 2017. - Т. 36. - № 6. - С. 70-93.
22. Условия формирования ультрабазитов Алхадырского террей-на (Восточный Саян. Сибирь) по результатам комплексного изучения состава хромшпинелидов / Ю.П. Бенедюк, В.А. Симонов, А.С. Мехоношин, Т.Б. Колотилина, С.И. Ступаков, А.А. Дорошков // Геология и геофизика. - 2015. - Т. 56. -№ 9. - С. 1664-1680.
23. Связь платиноносных ультрамафит-мафитовых интрузивов с крупными изверженными провинциями (на примере Сибирского кратона) / А.С. Мехоношин, Р. Эрнст, У. Седерлунд, М.А. Гамильтон, Т.Б. Колотилина, А.Э. Изох, Г.В. Поляков, Н.Д. Толстых // Геология и геофизика. - 2016. - Т. 57. -№5.- С. 1043-1057.
24. Арискин А.А., Бармина Г.С. Моделирование фазовых равновесий при кристаллизации базальтовых магм. - М.: Наука. «МАИК/Интерпериодика», 2000. - 363 с.
25. AriskinA.A., Barmina G.S. COMAGMAT: Development of a magma crystallization model and its petrologic applications // Geochemistry International. - 2004. - V. 42 (Supp. 1). -P. S1-S157.
26. The COMAGMAT-5: Modeling the Effect of Fe-Ni Sulfide Im-miscibility in Crystallizing Magmas and Cumulates / A.A. Ariskin, K.A. Bychkov, G.S. Nikolaev, G.S. Barmina // Journal of Petrology. - 2018. - V. 59. - № 2. - P. 283-297.
27. Danyushevsky L.V., Plechov P.Yu. Petrolog 3: Integrated software for modeling crystallization processes // Geochemistry. Geophysics. Geosystems. - 2011. - V. 12. - № 7. C. Q07021. DOI: 10.1029/2011GC003516.
28. Мехоношин А.С., Колотилина Т.Б., Дорошков А.А. Геохимическая модель формирования платиноносного дунит-верлито-вого массива Медек (Восточный Саян. Россия) // Геология и геофизика. - 2018. - Т. 59. - № 12. - С. 2011-2026.
29. Симонов В.А., Шарков Е.В., Ковязин С.В. Петрогенезис Fe-Ti интрузивных комплексов в районе Сьерра-Леоне. Центральная Атлантика // Петрология. - 2009. - Т. 17. - № 5. -С. 521-538.
30. Sobolev A.V., Danyushevsky L.V. Petrology and Geochemistry of Boninites from the North Termination of the Tonga Trench: Constraints on the Generation Conditions of Primary High-Ca Boninite Magmas // Journal of Petrology. - 1994. - V. 35. - P. 1183-1211.
31. Перчук Л.Л. Пироксеновый барометр и пироксеновые геотермометры // Доклады Академии наук. - 1980. - Т. 233. -№6.- C. 1196-2000.
32. Lindnsley D.H., Dixon S.A. Pyroxene thermometry // American Mineralogist. - 1983. - V. 68. - P. 477-493.
33. Ashchepkov I.V. Clinopyroxene Jd Barometer for mantle perido-tites and eclogites and thermal conditions of the lithospheric keels of cratons and their surroundings. A Geo Odyssey // GSA Annual Meeting. - Boston, 2001. - P. ID 11658.
34. Monomineral universal clinopyroxene and garnet barometers for peridotitic, eclogitic and basaltic systems / I.V. Ashchepkov, T. Ntaflos, A.M. Logvinova, Z.V. Spetsius, H. Downes, N.V. Vla-dykin // Geoscience Frontier. - 2017. - V. 8. - № 4. - P. 775-795.
35. Schmidt M.W. Amphibole composition in tonalite as a function of pressure: an experimental calibration of the Al-in-hornblende barometer // Contrib. Mineral. - Petrology. - 1992. - V. 110. - P. 304-310.
36. Hammarstrom J.M., Zen E.A. Aluminum in hornblende: an empirical igneous geobarometer // American Mineralogist. - 1986. -V. 71. - P. 1297-1313.
37. Anderson J.L., Smith D.R. The effects of temperature and fO2 on the Al-in-hornblende barometer // American Mineralogist. -1995. - V. 80. - P. 549-559.
38. Krawczynski M.J., Grove T.L., Behrens H. Amphibole stability in primitive arc magmas: effects of temperature, H2O content, and oxygen fugacity // Contribution to Mineralogy and Petrology. -2012. - V. 164. - P. 317-339.
39. Симакин А.Г., Шапошникова О.Ю. Новый амфиболовый геобарометр для высокомагнезиальных андезитовых и базальтовых магм // Петрология. - 2017. - Т. 25. - № 2. - C. 215-230.
40. Nomenclature of amphiboles: Report of the Subcommittee on Amphiboles of the International mineralogical Association, Commission on New Minerals and Mineral Names / B.E. Leake, A.R. Woolley, C.E.S. Arps, W.D. Birch, M.C. Gilbert, J.D. Gri-ce, F.C. Hawthorne, F. Kato, H.J. Kisch, V.G. Krivovichev, K. Linthout, J. Laird, J.A. Mandarino, W.V. Maresch, E.H. Nickel, N.M.S. Rock, J.C. Schumacher, D.C. Smith, N.C.N. Stephenson, L. Ungaretti, E.J.W. Whittaker, Y.Z. Guo // American Mineralogist. - 1997. - V. 82. - P. 1019-1037.
41. Симонов В.А., Шелепаев Р.А., Котляров А.В. Физико-химические параметры формирования расслоенного габбро-гиперба-зитового комплекса в офиолитах Южной Тувы // Ультрабазит-базитовые комплексы складчатых областей и связанные с ними месторождения: Материалы третьей международной конференции. - Екатеринбург: Институт геологии и геохимии УрО РАН, 2009. - Т. 2. - С. 195-198.
42. Магматические горные породы. - М.: Наука, 1983. - Т. 1. -766 с.
Поступила 04.06.2019 г.
Информация об авторах
Ступаков С.И., кандидат геолого-минералогических наук, старший научный сотрудник Института геологии и минералогии им. В.С. Соболева Сибирского Отделения Российской Академии Наук.
Симонов В.А., доктор геолого-минералогических наук, главный научный сотрудник Института геологии и минералогии им. В.С. Соболева Сибирского Отделения Российской Академии Наук; профессор Новосибирского национального исследовательского государственного университета; профессор Казанского федерального университета.
Мехоношин А.С., кандидат геолого-минералогических наук, старший научный сотрудник Института геохимии им. А.П. Виноградова Сибирского Отделения Российской Академии Наук; профессор Иркутского национального исследовательского технического университета.
Колотилина Т.Б., кандидат геолого-минералогических наук, старший научный сотрудник Института геохимии им. А.П. Виноградова Сибирского Отделения Российской Академии Наук; доцент Иркутского национального исследовательского технического университета.
UDC 550.4: 548.4
EVOLUTION OF PHYSICO-CHEMICAL CONDITIONS OF CRYSTALLIZATION OF MELTS DURING FORMATION OF DUNITE-PERIDOTITE-GABBROIC MASSIFS OF THE EASTERN SAYAN
Sergey I. Stupakov1,
stupakov@igm.nsc.ru
Vladimir A. Simonov1,2,3,
simonov@igm.nsc.ru
Sergey A. Mekhonoshin4,5,
mekhonos@igc.irk.ru
Tatyana B. Kolotilina4,5,
tak@igc.irk.ru
1 V.S. Sobolev Institute of Geology and Mineralogy SB RAS,
3, Academician Koptyug Avenue, Novosibirsk, 630090, Russia.
2 Novosibirsk National Research State University, 1, Pirogov street, Novosibirsk, 630090, Russia.
3 Kazan Federal University,
18, Kremlevskaya street, Kazan, 420008, Russia.
4 A.P. Vinogradov Institute of Geochemistry SB RAS, 1A, Favorsky street, Irkutsk, 664033, Russia.
5 Irkutsk National Research Technical University, 83, Lermontov street, Irkutsk, 664074, Russia.
The relevance of the research. The study of melt inclusions in minerals is one of the most reliable ways of obtaining direct information on the nature of magmatic processes, allows building a model of formation of intrusions and tracing the evolution of physico-chemical conditions of crystallization of magmatic complexes.
The main aim of the research is to prove the magmatic origin of rocks of this massif, which makes it possible to consider the distribution of noble metals in melts, based on the available experimental data.
Objects: melt inclusions in Cr-spinel from dunites and wehrlites Medek basite-ultrabasite massifs.
Methods. The search for melt inclusions was carried out in accessory Cr-spinel sample isolated from dunites and wehrlite, to determine the nature of the high temperature experiments. For this purpose, the samples were taken (about 50-100 grains of Cr-spinel fraction 0,2-0,5 mm) and placed in a graphite microcontainer with internal dimensions in the first millimeters. Micro containers were tightly closed with graphite covers. Thus, the Cr-spinel in these containers were at temperatures above 1000 °C under reducing conditions as a result of the reaction of air oxygen with graphite. The experiments were carried out at high temperatures on the basis of existing methods of investigation of melt inclusions. For maximum transformation of the melt in inclusions into homogeneous glass, quenching into water was carried out. The choice of the temperature regime (1280-1300 °C) was determined primarily by the fact that the main task of high-temperature experiments was to melt the contents of inclusions and to obtain glass during hardening, which was then analyzed on a microprobe and a scanning microscope.
Result. Conducted research of homogenised melt inclusions glass showed that crystallization of Cr-spinel and containing wehrlite of the ultrabasic Medek massif (East Sayan) were involved in low potassium picrobasalt of melts of normal alkalinity, of the most similar petrochemical characteristics to the data on melt inclusions in Cr-spinel from oceanic complexes. The calculations were made for various mi-neralogical barometers and different models of crystallization on the basis of compositions of glasses of heated melt inclusions in Cr-spinel in the COMAGMAT program, which showed that the pressure in formation of massive rocks was about 4 kbar (2,8-4,7kbar); the combined use of mineralogical thermometers and calculations on the basis of glass compositions of heated melt inclusions in Cr-spinel according to PETROLOG and COMAGMAT programs allows us to speak about the formation of dunites from picritoids (MgO 17,85 wt. %) of the melt, since the crystallization of olivine at 1380 °C, the formation of wehrlite occurred mainly from picrobasalt (MgO 13,36 wt. %) of the melt with the crystallization of olivine starting from 1280-1275 °C, but in much smaller quantities than clinopyroxene, the intensive formation of which was at lower temperatures (starting from 1235-1230 °C); a comparative analysis of the results of computational modeling of melt and mineral compositions with data on melt inclusions and real olivines and clinopyroxenes indicates their consistency, as evidenced by the fact of the location of melt inclusions along the calculated trends of melt change, that is, experimental data on inclusions on the one hand confirm the reliability of the calculations, and on the other - this ratio shows that the inclusions are not an accident, but record the natural and theoretically justified evolution of the melt.
Key words:
Medek intrusion, ultrabasic rocks, dunites, wehrlite, Cr-spinel, melt inclusions.
The research was carried out within the State Task IGM SB RAS (project no. 0330-2016-0014) and IGC SB RAS (project no. 0350-2016-0030), supported by the Ministry of Science and Higher Education of the Russian Federation, agreement no.14.Y26.31.0029.
REFERENCES
1. Barnes S.J., Lightfoot P.C. Formation of magmatic nickel-sulfide ore deposits and processes affecting their copper and platinum-group element contents. Economic Geology 100th Anniversary Volume, 2005, pp. 179-214.
2. Cawthorn R.G., Barnes S.J., Ballhaus C., Malich K.N. Platinum group element. chromium and vanadium deposits in mafic and ul-tramafic rocks. Economic Geology 100th Anniversary Volume, 2005, pp. 215-249.
3. Barnes S.J., Cruden A.R., Arndt N.T., Saumur B.M. The mineral system approach applied to magmatic Ni-Cu-PGE sulphide deposits. Ore Geology Reviews, 2016, vol. 76, pp. 296-316.
4. Barnes S.J., Robertson J.C. Time scales and length scales in magma flow pathways and the origin of magmatic Ni-Cu-PGE ore deposits. Geoscience Frontier, 2019, vol. 10, no. 1, pp. 77-87.
5. Dodin D.A., Chernyshov N.M., Polferov D.V., Tarnovetskii L.L. Platinometallnye mestorozhdeniya Mira. T. 1. Kniga 1. Malosul-fidnye mestorozhdeniya v ritmichno rassloennykh kompleksakh [PGM Deposits of the World. Vol. 1. B. 1. Low sulfide PGM deposits in rhythmically Layered Complexes]. Moscow, Geoinfrmmark Publ., 1994. 279 p.
6. Naldrett A.J. World-class Ni-Cu-PGE deposits: key factors in their genesis. Mineralium Deposita, 1999, vol. 34, pp. 227-240.
7. Dodin D.A., Chernyshov N.M., Yatskevich B.A. Platinometalnye mestorozhdeniya Rossii [PGM Deposits of Russia]. St. Petersburg, Nauka Publ., 2000. 755 p.
8. Naldrett A.J. From the mantle to the bank: the life of a Ni-Cu-(PGE) sulfide deposit. South African Journal of Geology, 2010, vol. 113.1, pp. 1-32.
9. Sobolev A.V. Melt inclusions in minerals as a source of principle petrological information. Petrology, 1996, vol. 4, no 3, pp. 209-220.
10. Schiano P., Clocchiatti R., Lorand J.-P., Massare D., Deloule E., Chaussidon M. Primitive basaltic melts included in podiform chromites from the Oman ophiolite. Earth Planetary Science Letter, 1997, vol. 146, no. 3-4, pp. 489-497.
11. Borisenko L.S., Borovikov L.A., Zhitova L.M., Pavlova G.G. Composition of magmatogene fluids and factors of their geoche-mical specialization and metal-bearing capacity. Russian Geology and Geophysics, 2006, vol. 47, no. 12, pp. 1282-1300.
12. Ryabchikov I.D., Kogsrko L.N., Solovova I.P. Physicochemical conditions of magma formation at the base of the Siberian plume: Insight from the investigation of melt inclusions in the mey-mechites and alkali picrites of the Maimecha-Kotui province. Petrology, 2009, vol. 17, no. 3, pp. 287-299.
13. Simonov V.A., Prikhodko V.S., Kovyazin S.V., Tarnavskiy V.A. Crystallization conditions of dunites in the Konder platiniferous alkaline-ultramafic massif of the southeastern Aldan Shield. Russian Journal of Pacific Geology, 2010, vol. 4, no. 5, pp. 429-440.
14. Simonov V.A., Prikhodko V.S., Kovyazin S.V. Genesis of Platini-ferous massifs in the Southeastern Siberian Platform. Petrology, 2011, vol. 19, no. 6, pp. 549-567.
15. Vasilev Y.R., Gora M.P. Meimechite-picrite associations in Siberia, Primorye, and Kamchatka (comparative analysis and petro-genesis). Russian Geology and Geophysics, 2014, vol. 55, no. 8, pp. 959-970.
16. Simonov V.A., Vasiliev Y.R., Stupakov S.I., Kotlyarov A.V., Karmanov N.S. Petrogenesis of dunites of the Guli ultrabasic massif (northern Siberian Platform). Russian Geology and Geophysics, 2016, vol. 57, no. 12, pp. 1696-1715.
17. Bonova K., Spisiak J., Bona J., Kovacik M. Chromian spinels from the Magura Unit (Western Carpathians. Eastern Slovakia) -their petrogenetic and palaeogeographic implications. Geological Quarterly, 2017, vol. 61, no. 1, pp. 3-18.
18. Renna M.R., Tribuzio R., Ottolini L. New perspectives on the origin of olivine-rich troctolites and associated harrisites from the
Ligurian ophiolites (Italy). Journal of the Geological Society, 2016, vol. 173, no. 6, pp. 916-932.
19. Naumov V.B., Dorofeeva V.A., Girnis A.V., Yarmoluk V.V. Mean concentrations of volatile components, major and trace elements in magmatic melts in major geodynamic environments on Earth. I. Mafic melts. Geochemistry International, 2017, vol. 55, no. 7, pp. 629-653.
20. Mironov N.L., Portnyagin M.V. Coupling of redox conditions of mantle melting and copper and sulfur contents in primary magmas of the Tolbachinsky Dol (Kamchatka) and Juan de Fuca Ridge (Pacific Ocean). Petrology, 2018, vol. 26, no. 2, pp. 145-166.
21. Simonov V.A., Prikhodko V.S., Vasilev V.S., Kotlyarov A.V. Physicochemical Conditions of Crystallization of Rocks from Ultrabasic Massifs of the Siberian Platform. Russian Journal of Pacific Geology, 2017, vol. 11, no. 6, pp. 447-468.
22. Benedyuk Yu.P., Simonov V.A., Mekhonoshin A.S., Kolotili-na T.B., Stupakov S.I., Doroshkov A.A. Genesis of ultramafic rocks of the Alkhadyr terrane (East Sayan, Siberia): implications from the data on Cr-spinel compositions. Russian Geology and Geophysics, 2015, vol. 56, no. 9, pp. 1308-1321.
23. Mekhonoshin A.S., Ernst R., Sdderlund U., Hamilton M.A., Ko-lotilina T.B., Izokh A.E., Polyakov G.V., Tolstykh N.D. Relationship between platinum-bearing ultramafic-mafic intrusions and large igneous provinces (examplified by the Siberian Craton). Russian Geology and Geophysics, 2016, vol. 57 (5), pp. 822-833.
24. Ariskin A.A., Barmina G.S. Modelirovanie fazovykh ravnovesiy pri kristallizatsii bazaltovykh magm [Modeling phase equilibrium of basalt magmas]. Moscow, Nauka, MAIK/Interperiodika Publ., 2000. 363 p.
25. Ariskin A.A., Barmina G.S. COMAGMAT: Development of a magma crystallization model and its petrologic applications. Geochemistry International, 2004, vol. 42 (Supp. 1), pp. S1-S157.
26. Ariskin A.A., Bychkov K.A., Nikolaev G.S., Barmina G.S. The COMAGMAT-5: Modeling the Effect of Fe-Ni Sulfide Immiscibi-lity in Crystallizing Magmas and Cumulates. Journal of Petrology, 2018, vol. 59, no. 2, pp. 283-297.
27. Danyushevsky L.V., Plechov P.Yu. Petrolog 3: Integrated software for modeling crystallization processes. Geochemistry, Geophysics, Geosystems, 29 July, 2011, vol. 12, no. 7, pp. Q07021. DOI: 10.1029/2011GC003516.
28. Mekhonoshin A.S., Kolotilina T.B., Doroshkov A.A. Geochemical model for the formation of the Medek platinum-bearing dunite-wehrlite massif (East Sayan, Russia). Russian Geology and Geophysics, 2018, vol. 59, no. 12, pp. 822-833.
29. Simonov V.A., Sharkov E.V., Kovyazin S.V. Petrogenesis of the Fe-Ti Intrusive Complexes in the Sierra Leone Region, Central Atlantic. Petrology, 2009, vol. 17, no. 5, pp. 488-502.
30. Sobolev A.V., Danyushevsky L.V. Petrology and Geochemistry of Boninites from the North Termination of the Tonga Trench: Constraints on the Generation Conditions of Primary High-Ca Boni-nite Magmas. J. Petrol., 1994, vol. 35, pp. 1183-1211.
31. Perchuk L.L. Pyroxene barometer and pyroxene geotermometers. Doklady Ras, 1980, vol. 233, no. 6, pp.1196-2000. In Rus.
32. Lindnsley D.H., Dixon S.A. Pyroxene thermometry. Amer. Mineralogy, 1983, vol. 68. pp. 477-493.
33. Ashchepkov I.V. Clinopyroxene Jd Barometer for Mantle Perido-tites and Eclogites and Thermal Conditions of the Lithospheric Keels of Cratons and their Surroundings. A Geo Odyssey. GSA Annual Meeting, November 1-10, Boston, 2001. P. ID 11658.
34. Ashchepkov I.V., Ntaflos T., Logvinova A.M., Spetsius Z.V., Downes H., Vladykin N.V. Monomineral universal clinopyroxene and garnet barometers for peridotitic, eclogitic and basaltic systems. Geoscience Frontier, 2017, vol. 8, no. 4, pp. 775-795.
35. Schmidt M.W. Amphibole composition in tonalite as a function of pressure: an experimental calibration of the Al-in-hornblende barometer. Contrib. Mineral. Petrology, 1992, vol. 110, pp. 304-310.
36. Hammarstrom J.M., Zen E.A. Aluminum in hornblende: an empirical igneous geobarometer. American Mineralogist, 1986, vol. 71, pp. 1297-1313.
37. Anderson J.L., Smith D.R. The effects of temperature and fO2 on the Al-in-hornblende barometer. American Mineralogist, 1995, vol. 80, pp. 1297-1313.
38. Krawczynski M.J., Grove T.L., Behrens H. Amphibole stability in primitive arc magmas: effects of temperature, H2O content, and oxygen fugacity. Contribution to Mineralogy and Petrology, 2012, vol. 164, pp. 317-339.
39. Simakin A.G., Shaposhnikova O.Yu. New amphibole geobarome-ter for high-magnesium andesite and basalt magmas. Petrology, 2017, vol. 25, no. 2, pp. 215-230.
40. Leake B.E., Woolley A.R., Arps C.E.S., Birch W.D., Gilbert M.C., Grice J.D., Hawthorne F.C., Kato F., Kisch H.J., Kri-vovichev V.G., Linthout K., Laird J., Mandarino J.A., Ma-resch W.V., Nickel E.H., Rock N.M.S., Schumacher J.C.,
Smith D.C., Stephenson N.C.N., Ungaretti L., Whitta-ker E.J.W., Guo Y.Z., Nomenclature of amphiboles: Report of the Subcommittee on Amphiboles of the International mineralogi-cal Association, Commission on New Minerals and Mineral Names. American Mineralogist, 1997, vol. 82. pp. 1019-1037.
41. Simonov V.A., Shelepaev R.A., Kotlyarov V.A. Fiziko-khi-micheskie paramety formirovaniya rassloennogo gabbro-giperba-sitovogo kompleksa v ofiolitakh Yuzhnoy Tuvy [Physico-chemical parameters of formation of layered gabbro-ultrabasic complex in the ophiolites of the southern Tuva]. Ultrabazit-bazitovye komplek-sy skladchatykh oblastey b svyazannye s nimi mestorozhdeniya [Ultrabasic-basic complexes of folded regions and related deposits]. Ekaterinburg, IGG UB RAS Publ., 2009. Vol. 2, pp. 195-198.
42. Izverzhennye gornye porody [Igneous rocks]. Moscow, Nauka Publ., 1983. Vol. 1, 766 p.
Received: 4 June 2019.
Information about the authors
Sergey I. Stupakov, Cand. Sc., senior researcher, V.S. Sobolev Institute of Geology and Mineralogy SB RAS.
Vladimir A. Simonov, Dr. Sc., chief researcher, V.S. Sobolev Institute of Geology and Mineralogy SB RAS; professor, Novosibirsk National Research State University; professor, Kazan Federal University.
Aleksey S. Mekhonoshin, Cand. Sc., senior researcher, A.P. Vinogradov Institute of Geochemistry SB RAS; professor, Irkutsk National Research Technical University.
Tatyana B. Kolotilina, Cand. Sc., senior researcher, A.P. Vinogradov Institute of Geochemistry SB RAS; associate professor, Irkutsk National Research Technical University.
