Экспериментальное исследование плавления в форстеритнормативной части системы СаО-MgO-Al2O3-SiO2 в связи с существованием "эклогитового барьера"

Размещено на http://www.allbest.ru/

Институт геологии и минералогии СО РАН, г. Новосибирск

Экспериментальное исследование плавления в форстеритнормативной части системы СаО-MgO-Al2O3-SiO2 в связи с существованием "эклогитового барьера"

Банушкина С.В., Гартвич Ю.Г.,

Голицына З.Ф., Сурков Н.В.

Аннотация

Проведено экспериментальное исследование ликвидусных фазовых взаимоотношений в форстеритнормативной области фазовой диаграммы системы СаО-MgO-Al₂O₃-SiO₂ в области давления 2,5-3,0 ГПа. Установлено существование двух перитектических реакций плавления Opx+Fo+L=Cpx+Ga и Opx+L=Cpx+Ga и одной эвтектической реакции L=Cpx+Opx+Ga+Fo. Показано, что при давлении 2,5-2,6 ГПа и температуре 1500 ^OC существует нонвариантное равновесие (L,Cpx,Opx,Ga,Fo), в точке существования которого состав жидкости находится в плоскости Mg₂Si₂O₆-Al₂O₃-Ca₂Si₂O₆. Из точки этого нонвариантного равновесия в сторону высоких давлений выходит луч реакции L=Cpx+Opx+Ga+Fo, в сторону высоких температур и высоких давлений луч реакции Opx+L=Cpx+Ga, а в сторону низких давлений луч реакции Opx+Fo+L=Cpx+Ga. Показано, что сингулярное нонвариантное равновесие (L,Cpx,Opx,Ga,Fo) является моделью “эклогитового барьера”, которое обсуждалось в работах В.А. Жарикова и В.С. Соболева в связи с проблемой эволюции состава первичных ультраосновных магм до гранитных расплавов.

Ключевые слова: магматизм, фазовая диаграмма, нонвариантное равновесие, экспериментальное исследование, высокие давления.

Abstract

We conducted experimental study of the liquidus phase relationships in the forsterite-normative region of the phase diagram of the СаО-MgO-Al₂O₃-SiO₂ system in the pressure range of 2.5-3.0 GPa. The existence of two peritectic melting reactions Opx+Fo+L=Cpx+Ga и Opx+L=Cpx+Ga and one eutectic reaction L=Cpx+Opx+Ga+Fo is established. It is shown that at the pressure of 2.5-2.6 GPa and temperature of 1500 ^OC, there is a non-variant equilibrium (L, Cpx, Opx, Ga, Fo). At the point of its existence, the liquid composition is in the Mg₂Si₂O₆-Al₂O₃-Ca₂Si₂O₆ plane. From the point of view of this non-variant equilibrium, the reaction beam is L = Cpx + Opx + Ga + Fo in the direction of high temperatures and high pressures, while the reaction beam Opx + L = Cpx + Ga leaves the point of this non-invariant equilibrium towards high pressures. The reaction beam in the direction of low pressures is Opx + Fo + L = Cpx + Ga. It is shown that the singular non-variant equilibrium (L, Cpx, Opx, Ga, Fo) is a model of the “eclogite barrier,” which was discussed in the works of V.А.Zharikov and V.S.Sobolev in connection with the problem of the composition of primary ultra-basic magmas evolution to granite melts.

Keywords: magmatism, phase diagram, non-invariant equilibrium, experimental investigation, high pressures.

Одним из нерешённых вопросов современной петрологии является проблема, каким путём из материала ультраосновной мантии формируется Земная кора, состоящая преимущественно из пород гранитного состава. К настоящему времени проведено множество экспериментальных исследований и построен ряд умозрительных моделей на эту тему. В части работ этот вопрос обсуждался под названием “эклогитовый барьер” [1], [2], но физико-химического объяснения этой проблемы до сих пор не было предложено.

В связи с этим возникла задача разработать на основе физико-химических исследований модель этого явления, связанную с фазовыми диаграммами модельных алюмосиликатных систем при высоких давлениях. Для этого в качестве объекта для экспериментальных исследований выбраны фазовые взаимоотношения в области плавления в форстеритнормативной части фазовой диаграммы системы СаО-MgO-Al₂O₃-SiO₂.

Фазовая диаграмма системы СаО-MgO-Al₂O₃-SiO₂ (рис. 1) является наиболее популярной моделью для таких исследований. Оксиды Si, Mg, Al и Са составляют 80-90 мас. % от содержания в таких породах, как лерцолиты, дуниты, гранатовые лерцолиты, эклогиты, гарцбургиты и т. п.; кроме того в системе СаО-MgO-Al₂O₃-SiO₂ образуются практически все фазы, соответствующие минералам глубинных пород. Физико-химической моделью перехода состава магматического расплава от ультраосновных к основным и средним может служить форстеритнормативная область составов фазовой диаграммы системы СаО-MgO-Al₂O₃-SiO₂. Из результатов экспериментальных исследований при высоких давлениях известно, что ассоциация форстерит+клинопироксен+ортопироксен+гранат при давлениях 3,0-4,0 ГПа испытывает плавление эвтектического типа по реакции L=Fo +Срх+Орх+Ga [3], [4], [5]. При более низких давлениях эта эвтектика исчезает, и плавление носит более сложный характер [6], [7].

Рис. 1 _ Фазовые соотношения в системе CaO-MgO-Al₂O₃-SiO₂ при давлении 3,0 ГПа и температуре 1200 ^OC

В связи с этим возникла задача исследования НОН- и моновариантных равновесий возникающих при разрушении реакции эвтектического плавления L=Fo+Срх+Орх+Ga на фазовой диаграмме системы СаО-MgO-Al₂O₃-SiO₂ при понижении давления.

Методика экспериментальных исследований

Исследование проводили на аппарате типа “поршень-цилиндр” [8], [9] по методике, мало отличающейся от общепринятой для этого типа аппаратов. В качестве нагревательного устройства использована ячейка на основе хлорида натрия, подобная той, которую предложил Мирвальд и др. [10]. Выбор именно этой ячейки связан с тем, что она позволяет, не только наиболее точно контролировать давление, но и не содержит водосодержащих материалов, что даёт возможность проводить эксперименты по исследованию ликвидуса системы.

Давление определяли по усилению нагрузки за вычетом поправки на трение, которую вычисляли по разнице усилия нагрузки и расчётного усилия, необходимого для фазового перехода Bi_I_Bi_III(2,54 ГПа) при комнатной температуре. Для калибровки висмутовую проволочку длиной 3 мм в таблетке из хлористого серебра помещали непосредственно на поршень. Точность поддержания и измерения давления нами не определялась, но по оценкам разных исследователей [10] на “ячейке низкого трения” составляет около ±0,03 ГПа. Для определения и регулировки температуры использовали платинородиевую термопару PtRd6_PtRd30, управление температурой производили при помощи высокоточного регулятора ВРТ-3 (точность поддержания ±1 ^OС). Влияние давления на Э.Д.С. термопары не учитывали. Точность измерения температуры по этой методике оценивается не хуже ±10 ^OC [9].

Исходные вещества приготовлены весовым методом из прокаленных (1100_1400 ^OC; 5_6 часов) оксидов марки Ос.Ч. Кроме того, использовали оксид алюминия, который получен прокаливанием гидрооксида алюминия марки Ос.Ч.-18-3. Из них были приготовлены смеси, отвечающие по стехиометрии CaMgSi₂O₆, Mg₃Al₃Si₃O₁₂ и Mg₂SiO₄, которые были раскристаллизованы при 1450--1550 ^OC и 1 атм в ассоциации фаз субсолидуса. Прокаливание смесей через 8_10 ч чередовали с растиранием в ступке из твёрдого сплава до полного исчезновения исходных оксидов.

Диопсид и форстерит легко синтезируются при атмосферном давлении. Смесь, отвечающая по составу пиропу, была расплавлена в прозрачное стекло при атмосферном давлении и температурах несколько выше 1650 ^OС. Поскольку для исследования требовалась смесь со стехиометрией, смещённой в сторону форстерита, то эта смесь, обозначенная как “Di40Pyr605Fo”, приготовлена добавлением к смеси состава 40 масс.% Di 60 масс.% Pyr - ещё 5 масс.% Fo, дополнительно прибавленных к её изначальному весу. Характеристики использованных смесей приведены в таблице 1. Для проведения эксперимента образец помещали в платиновую ампулу, просушивали при 500_600 ^OC в течение 7_8 ч., затем герметизировали электросваркой.

Для исследования образцов изготавливали двусторонний полированный шлиф, который изучали обычными петрографическими методами с помощью поляризационного микроскопа. При изучении ликвидуса оптическое исследование шлифа является единственным методом, показывающим однозначные признаки плавления состава. Шлиф всегда изготавливали ориентировано, из поперечного скола образца. Дополнительно анализ фазового состава продуктов экспериментов проводили рентгенофазовым методом, а анализ состава фаз проводили на электронном микроанализаторе “Comebax-Micro”.

Зерна фаз в продуктах опытов даже в области плавления имели размер 5_20 микрон. Это затрудняет изучение состава фаз с помощью электронного микроанализатора. Увеличение продолжительности экспериментов не даёт заметного возрастания размера зёрен фаз. Кроме того, при плавлении зёрна фаз отседают на дно ампулы и часто не попадают в область среза петрографического шлифа. По этим причинам в части экспериментов не удаётся получить удовлетворительные анализы всех обнаруженных фаз (к рассмотрению приняты анализы с суммой более 97-98 %).

Жидкость закаливается в виде стекла не во всех экспериментах, при температурах вблизи солидуса происходит раскристаллизация жидкости в закалочные фазы, представленные преимущественно клинопироксеном, реже форстеритом. Закалочные фазы имеют характерную структуру типа “петушиного хвоста” либо “мороза на окнах”, по которой их легко отличить от других фаз. Анализ состава закалочных фаз с помощью электронного микроанализатора не даёт информации о составе исходной жидкости, поскольку состав заметно изменяется от точки к точке, и в некоторых случаях общая сумма оксидов оказывается значительно меньше 100 % (до 95 %). Последнее явление, скорее всего, связано с невозможностью хорошо отполировать поверхность закалочных фаз. На их поверхности остаётся заметная шагрень. В более высокотемпературных опытах жидкость удаётся частично закалить в стекло. Однако с течением времени (1_2 месяца) стекло постепенно переходит в агрегат закалочных фаз.

Таблица 1 _ Характеристика исходных веществ

Результаты экспериментов

Реакция L=Орх+Срх+Ga+Fo изучена нами при давлении 2,8 и 3,0 ГПа. Результаты экспериментов представлены в таблице 2 и на рис. 2, и свидетельствуют о крайне пологом наклоне луча изучаемой реакции в этом интервале давлений. Наши результаты не противоречат данным Девиса [3], [4]. Более высокие температуры плавления в его исследовании вполне объяснимы за счёт более высоких давлений и изменения состава всех участвующих в равновесии фаз.

Моновариантные реакций и реакция Орх+L=Ga+Срх изучены при помощи смесей, имеющих состав 40 масс.% Di 60 масс.% Pyr и 30 масс.% Di 70 масс.% Pyr. Состав стекла в экспериментах, в которых зафиксировано плавление, не совпадает с исходным. Всегда наблюдается обогащение оксидом кальция, что свидетельствует о присутствии магнезиальных твёрдых фаз типа ортопироксена или оливина. Условия плавления ассоциации Ga+Срх при давлении 3,0 ГПа хорошо совпадают с данными О'Хары и Йодера [5].

Рис. 2 _ Результаты экспериментов по исследованию лучей моновариантных реакций Opx+Fo+L=Cpx+Ga, Opx+L=Cpx+Ga и L=Cpx+Opx+Ga+Fo. 1 - конода Cpx+Ga для состава S-67; 2 - плавление для состава S-67; 3 - конода Cpx+Ga для состава S-69; 4 - плавление для состава S-69; 5 - ассоциация Cpx+Opx+Ga+Fo; 6 - плавление для ассоциации Cpx+Opx+Ga+Fo

Таблица 2 _ Условия проведения экспериментов

При давлении выше 2,6 ГПа в качестве ликвидусной фазы продиагностирован только ортопироксен, а при более низких давлениях к нему добавляется форстерит. Ликвидусные фазы видны в шлифах в придонных частях образца в виде одиночных, хорошо огранённых мелких кристаллитов, плавающих в стекле. Размер этих кристаллов меньше толщины шлифа и не всегда видимый кристалл фазы выходит на поверхность шлифа, что не позволило определить состав фазы с помощью электронного микроанализатора во всех опытах. Однако форстерит и ортоэнстатит значительно различаются по двупреломлению, что позволяет легко их различать при петрографическом изучением шлифа. Следует отметить, что в экспериментах ниже 2,6 ГПа меняется и облик области шлифа с ликвидусными фазами: помимо появления изометричных с высоким двупреломлением кристаллов форстерита, кристаллиты ортопироксена укрупняются и становятся более изометричными. По этим визуально наблюдаемым признакам при давлении 2,6 ГПа происходит смена характера плавления в обеих смесях состава 40 масс.% Di 60 масс.% Pyr и 30 масс.% Di 70 масс.% Pyr. Это даёт основания предположить, что эксперименты при давлении выше 2,6 ГПа соответствуют продуктам реакции типа Орх+L=Ga+Срх, а при более низком давлении реакции Fo+Орх+L=Ga+Срх.

Таблица 3 _ Составы фаз в продуктах экспериментов, проанализированные с помощью электронного микроанализатора

Топологический анализ фазовых соотношений и обсуждение результатов

Полученные экспериментальные результаты свидетельствуют, что ассоциация Cpx+Opx+Ga+Fo плавится эвтектически согласно реакции L=Cpx+Opx+Ga+Fo. Эти данные совпадают с выводами из экспериментальных исследований Дависа и Шайрера [3], [4].

Плавление ассоциации клинопироксена и граната при давлении 3,0 ГПа хорошо совпадает с результатами О'Хары и Йодера [5]. При давлениях выше 2,6 ГПа в области ликвидусной фазой является ортопироксен, при более низких давлениях появляется форстерит.

Из этих фазовых взаимоотношений следует возможное существование двух перитектических реакций плавления Opx+Fo+L=Cpx+Ga и Opx+L=Cpx+Ga. Если учесть существование эвтектической реакции L=Cpx+Opx+Ga+Fo, то в районе давления 2,5_2,6 ГПа однозначно должно существовать нонвариантное равновесие (L,Cpx,Opx,Ga,Fo), в точке существования которого состав жидкости находится в плоскости Mg₂Si₂O₆-Al₂O₃-Ca₂Si₂O₆. Так как составы клинопироксена, ортопироксена и граната находятся в этой же плоскости, то это равновесие имеет вырожденный характер, система становится условно трёхкомпонентной. Форстерит становится сингулярной фазой.

Анализ топологии этого участка фазовой диаграммы указывает, что при давлении около 2,6 ГПа и температуре 1500 ^OC существует точка нонвариантного равновесия (L,Cpx,Opx,Ga,Fo), из которой в сторону высоких давлений выходит луч реакции L=Cpx+Opx+Ga+Fo, в сторону высоких температур и высоких давлений луч реакции Opx+L=Cpx+Ga, а в сторону низких давлений луч реакции Opx+Fo+L=Cpx+Ga (рис. 2).

Строение исследованного участка фазовой диаграммы системы СаО-MgO-Al₂O₃-SiO₂ полностью подтверждает данные К. Хейзенберг, а также Т. Като и М. Кумазава [11, 12, 13, 14] о смешении состава расплава на эвтектике L=Cpx+Opx+Ga+Fo при понижении давления практически от чистого оливина до плоскости Mg₂Si₂O₆-Al₂O₃-Ca₂Si₂O₆ при 2,6 ГПа и далее в области низких давлений перехода состава расплава в кварцнормативную область (рис. 3).

кора мантия эклогитовый барьер

Рис. 3 _ Тренд изменения состава ультраосновного расплава при понижении давления. Точка е _ точка пересечения составом расплава плоскости Mg₂Si₂O₆-Al₂O₃-Ca₂Si₂O₆

Сингулярное нонвариантное равновесие (L,Cpx,Opx,Ga,Fo) является моделью “эклогитового барьера”, обсуждаемого в работах В.А. Жарикова и В.С. Соболева [1, 2], и позволяет на физико-химической основе объяснить эволюцию состава первичных ультраосновных магм в сторону составов с более высоким содержанием кремнезёма, по-видимому, вплоть до составов, соответствующим гранитным магмам.

Основные результаты

Экспериментально исследованы моновариантные реакции плавления в форстеритнормативной части системы СаО-MgO-Аl₂O₃-SiO₂ при высоких давлениях. Изучены лучи моновариантных реакций: Fo+Opx+L=Cpx+Ga, Opx+L=Cpx+Ga и L=Fo+Cpx+Opx+Ga Установлены составы сосуществующих фаз в экспериментах на ликвидусе и солидусе.

1. На основе экспериментальных данных при помощи топологического анализа проведена увязка лучей моновариантных реакций и установлено положение нонвариантного равновесия (L,Cpx,Opx,Ga,Fo).

2. Показан путь эволюции состава магматического расплава при понижении давления от ультраосновного до кварцнормативного, позволяющий преодолеть “эклогитовый барьер”.

Работа выполнена в рамках проекта НИР (0330-2016-0016) “Моделирование тепловой и гидродинамической структуры плюмов для определения условий формирования магматических расплавов и их влияния на вещественный состав и структуру литосферы Северной Евразии”.

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

1. Жариков В. А. Эклогитовый барьер и клинопироксены высоких давлений / В. А. Жариков, Р. А. Ишбулатов, Л. Т. Чудиновских // Геология и геофизика. _ 1984. _ № 12. _ С. 54_63.

2. Соболев В. С. Состав глубинных пироксенов и проблема эклогитового барьера / В. С. Соболев, А. В. Соболев // Геология и геофизика. _ 1977. _ № 12. _ С. 46_59.

3. Davis В. Т. C. System diopside-forsterite-pyrope at 40 kilobars / В. Т. C. Davis // Carnegie Institution Washington Yearbook. _ 1964. _ V. 63. _ P. 165_171.

4. Davis В. Т. C. Melting relations in join diopside-forsterite- pyrope at 40 kilobars and one atmosphere / В. Т. C. Davis, J. F. Schairer // Carnegie Institution Washington Yearbook. _ 1965. _ V. 64. _ P. 123_126.

5. O'Hara M. J. Formation and fractionation of basic magma at high pressures / M. J. O'Hara, H. S. Jr. Yoder // Scottish Journal of Geology. _ 1967. _ V. 3. _ №1. _ P. 67_117.

6. Presnall D. C. Liquidus phase relations on join diopside-forsterite-anorthite from 1 atm to 20 kbar: their beaming on the generation and crystallization of basaltic magma / D.C. Presnall, S.A. Dixon, J.R. Dixon et al. // Contributions to Mineralogy and Petrology. _ 1978. _ V. 66. _ №2. _ P. 203_220.

7. Presnall D. C. Generation of mid- ocean-ridge tholeiites / D. C. Presnall, J. R. Dixon, Т. H. O'Donnell et al. // Journal of Petrology. -- 1979. -- V. 20. _ № 1. -- P. 3--35.

8. Годовиков А. А. Аппарат для создания давления до 40 кбар при температурах до 1700 ^ОС / А. А. Годовиков, С. А. Смирнов, И. Ю. Малиновский и др. // Приборы и техника эксперимента. _ 1971. _ №6. _ С. 159_160.

9. Boyd F. R. Apparatus for phase-equilibrium measurements at pressures up to 50 kilobars and temperatures up to 1750 °C / F. R. Boyd, J. L. England // Journal of The Geophysical Research. _ 1960. _ V. 65, _ №2. _ P. 741_748.

10. Mirvald P. W. Low-friction cell for piston- cylinder high-pressure apparatus / P. W. Mirvald, I. C. Getting, G. C. Kennedy // Journal of The Geophysical Research. _ 1975. _ V. 86, _ № 11. _ P. 1519_1525.

11. Herzberg C. Origin of Mantle Peridotite: Constraints From Melting Experiments to 16.5 GPa / C. Herzberg, T. Gasparik, H. Sawamoto // Journal of the Geophysical Research. _ 1990. _ V. 95. _ №10. _ P. 15779-15803.

12. Herzberg C. T. Solidus and liquidus temperatures and mineralogies for anhydrous garnet_lherzolite to 15 GPa / C. T. Herzberg // Physics of the Earth and Planetary Interiors. _ 1983. _ V. 32. _ №2. _ P. 193_202.

13. Herzberg C. T. Origin of mantle peridotite and komatiite by partial melting / C. T. Herzberg, M. J. O'Hara // Geophysical Research Letters. _ 1985. _ V. 12. _ № 9. _ P. 541-544.

14. Kato T. Effect of high pressure on the melting relations in the system Mg₂SiO₄-MgSiO₃ Part I. Eutectic relations up to 7 GPa / T. Kato, M. Kumazawa // Journal of Physics of the Earth. _ 1985. _ V. 33. _ №6. _ P. 513_524.

Список литературы на английском языке / References in English

1. Zharikov V. A. Jeklogitovyj bar'er i klinopirokseny vysokih davlenij [An eclogite barrier and clinopyroxenes at high pressures] / V. A. Zharikov, R. A. Ishbulatov, L. T. Chudinovsky // Geologija i geofizika [Geology and geophysics.]. _ 1984. _ №12. _ P. 54-63. [in Russian]

2. Sobolev V. S. Sostav glubinnyh piroksenov i problema jeklogitovogo bar'era. [Composition of deep pyroxenes and the problem og eclogite barrier] / V. S. Sobolev, A. V. Sobolev // Geologija i geofizika [Geology and geophysics]. _ 1977. _ №12. _ P. 46-59. [in Russian]

3. Davis В. Т. C. System diopside-forsterite-pyrope at 40 kilobars/ В. Т. C. Davis // Carnegie Institution Washington Yearbook. _ 1964. _ V. 63. _ P. 165_171.

4. Davis В. Т. C. Melting relations in join diopside-forsterite- pyrope at 40 kilobars and one atmosphere / В. Т. C. Davis, J. F. Schairer // Carnegie Institution Washington Yearbook. _ 1965. _ V. 64. _ P. 123_126.

5. O'Hara M. J. Formation and fractionation of basic magma at high pressures / M. J. O'Hara, H. S. Jr. Yoder // Scottish Journal of Geology. _ 1967. _ V. 3. _ №1. _ P. 67_117.

6. Presnall D. C. Liquidus phase relations on join diopside-forsterite-anorthite from 1 atm to 20 kbar: their beaming on the generation and crystallization of basaltic magma / D.C. Presnall, S. A. Dixon, J. R. Dixon et al. // Contributions to Mineralogy and Petrology. _ 1978. _ V. 66. _ №2. _ P. 203_220.

7. Presnall D. C. Generation of mid- ocean-ridge tholeiites / D. C. Presnall, J. R. Dixon, Т. H. O'Donnell et al. // Journal of Petrology. -- 1979. -- V. 20. _ № 1. -- P. 3--35.

8. Godovikov A. A. Apparat dlja sozdanija davlenija do 40 kbar pri temperaturah do 1700 ^OC [Apparat for creation of pressure to 40 кбар at temperatures to 1700 ^OC] / A .A. Godovikov, S. A. Smirnov, I. Ju. Malinovsky and others // Pribory i tehnika jeksperimenta [Devices and technics of experiment]. _ 1971. _ №6. _ P. 159-160. [in Russian]

9. Boyd F. R. Apparatus for phase-equilibrium measurements at pressures up to 50 kilobars and temperatures up to 1750 °C / F. R. Boyd, J. L. England // Journal of The Geophysical Research. _ 1960. _ V. 65, _ №2. _ P. 741_748.

10. Mirvald P. W. Low-friction cell for piston- cylinder high-pressure apparatus / P. W. Mirvald, I. C. Getting, G. C. Kennedy // Journal of The Geophysical Research. _ 1975. _ V. 86. _ № 11. _ P. 1519_1525.

11. Herzberg C. Origin of mantle peridotite: constraints from melting experiments to 16.5 GPa / C. Herzberg, T. Gasparik, H. Sawamoto // Journal of the Geophysical Research. _ 1990. _ V. 95. _ №10. _ P. 15779-15803.

12. Herzberg C. T. Solidus and liquidus temperatures and mineralogies for anhydrous garnet_lherzolite to 15 GPa / C. T. Herzberg // Physics of the Earth and Planetary Interiors. _ 1983. _ V. 32. _ №2. _ P. 193_20213.

13. Herzberg C. T. Origin of mantle peridotite and komatiite by partial melting / C.T. Herzberg, M.J. O'Hara // Geophysical Research Letters. _ 1985. _ V. 12. _ № 9. _ P. 541-544.

14. Kato T. Effect of high pressure on the melting relations in the system Mg₂SiO₄-MgSiO₃ Part I. Eutectic relations up to 7 GPa. / T. Kato, M. Kumazawa // Journal of Physics of the Earth. _ 1985. _ V. 33. _ №6. _ P. 513_524.

Размещено на Allbest.ru