Распределение Fe2+/Mg отношения в системе расплав - шпинель - оливин

(по опубликованным экспериментальным данным); неравновесность составов оливин – шпинелевых включений как показатель возможного твердофазного генезиса алмазов.

Г.П.Пономарев, М.Ю.Пузанков

В результате обработки большого количества опубликованных экспериментальных данных получены простые линейные зависимости равновесного распределения Fe²/Mg отношения между сосуществующими основным - ультраосновным расплавом и кристаллами шпинели и оливина для широкого диапазона составов, давлений до 1,5 Гпа, при варьирующей температуре. Численная оценка равновесности распределения Fe²/Mg отношения позволяет выделять природные оливин-шпинелевые парагенезисы магматического (расплавного) происхождения и отличать их от оливин-шпинелевых пар, изменивших свои составы при метаморфизме, или имеющих метаморфический генезис.

Такая оценка выявляет неравновесность с магматическим расплавом шпинелевых и оливиновых включений в алмазах Якутской кимберлитовой провинции, что указывает на возможный твердофазный рост вмещающих их кристаллов алмаза в мантийных(?) условиях. В коровых условиях твердофазным ростом микрокристаллов алмаза во время палеоземлетрясений можно объяснить их генезис на месторождении Кумды-Коль (Казахстан). Выделен ряд признаков для прогноза и поиска аналогичных месторождений в областях развития гранито-гнейсовых куполов. В частности, в Камчатском регионе подобные рудопроявления (месторождения?) могут существовать в Срединном Камчатском массиве (Хангарский гранито-гнейсовый купол).

В породах базит-гипербазитового ряда оливин является типичным породообразующим силикатом, а минералы группы шпинели - преобладающими акцессориями. В генетических моделях для пород этого ряда часто используются числовые зависимости, связывающие составы этих минералов с условиями их образования и последующими изменениями. Исследования ведутся по трем направлениям: эксперимент в контролируемых условиях; термодинамическое моделирование; анализ природных оливин - шпинелевых парагенезисов.

В области между ликвидусом и солидусом по экспериментальным данным получены распределения Fe²⁺/Mg отношения в системах расплав - оливин и расплав - шпинель; в субсолидусной - в системе оливин-шпинель путем моделирования и привлечения данных по природным парагенезисам. В системе расплав - оливин при атмосферном давлении в условиях контролируемых температуры и фугитивности кислорода калибровочная зависимость для Fe²⁺/Mg была получена П.Редером и Р.Эмсли (89). По опубликованным данным экспериментальных исследований с подобным же контролем условий Г.П.Пономаревым и его соавторами (53) установлена числовая зависимость, связывающая Fe²⁺/Mg отношение в основном - ультраосновном расплаве и шпинели. При этом, благодаря учету влияния содержаний Ti в шпинели, коэффициент корреляции составил 0,98. Распределение Fe²⁺/Mg между оливином и шпинелью в субсолидусной - солидусной области, с учетом влияния некоторых элементов, было откалибровано рядом исследователей как геотермометр (80,87,88,90,92) и предложено в качестве геоспидометра (88). Однако для многих базит-гипербазитовых пород применение указанных выше зависимостей ограничено влиянием неучтенного давления и ошибками метода. Кроме того, можно указать и на то, что при кристаллизации, гибридизме, остывании расплавов и последующих метаморфических преобразованиях пород происходит изменение составов стекол и расплавных включений в минералах, а также и самих минералов. Первоочередной задачей данной работы было частичное снятие этих ограничений и получение критерия для отличия магматических (равновесных и неравновесных) OL - SP парагенезисов от немагматических. Необходимость и важность такого критерия для корректировки петрологических выводов очевидна. Одно из приложений имеет следствия, которые, ввиду их большого прикладного значения, несомненно, заслуживают более подробного рассмотрения. Установленная на основании полученных в первой части этой работы зависимостей, неравновесность OL-SP включений в кристаллах алмаза Якутской кимберлитовой провинции, позволяет предположить их общий твердофазный генезис. Рост алмаза в твердой среде возможен не только в мантийных, но и в коровых условиях, о чем свидетельствует ряд экспериментальных результатов (12,24,82). Такая возможность, по-видимому, реализована в микрокристаллах алмазов месторождения Кумды-Коль (Казахстан). Обсуждение генетической модели роста таких микрокристаллов и критериев поиска месторождений, подобных месторождению Кумды-Коль, и завершает данную работу.

Использованные экспериментальные данные.

Из опубликованных работ были выбраны (77,83,85,91,93,97,98,99,100), удовлетворяющие следующим условиям: эксперименты проводились при атмосферном давлении с контролем температуры и фугитивности кислорода, приведены полные составы сосуществующих расплавов (стекол), шпинелей и оливинов. Был сформирован массив, состоящий из 92 точек многомерного пространства признаков. Каждая такая точка представляет собой как бы объединенный "анализ" - матрицу, состоящую из пересчитанных на атомные количества содержаний химических элементов в стеклах, шпинели и оливине, значений температуры и фугитивности.

Методика обработки данных. Включает пересчеты первичных составов стекол, шпинелей и оливина, полученных зондовым методом, и дальнейшую математическую обработку методом наименьших квадратов.

Составы стекол. В имеющихся составах стекол содержания FeO и Fe₂O₃ (если такое имелось, в виде весовых % окислов) пересчитывались в FeO (общее). Затем разделение на FeO и Fe₂O₃ проводилось по обобщенной формуле, предложенной в (10). После этого рассчитывались атомные количества элементов: сначала с учетом содержания кислорода отношение K/O, где:

K = Si + Ti + Al + Cr + Fe³⁺ + Fe²⁺ + Mn + Ca + Na + K; O - кислород.

K/O варьировало от 41/59 до 38/62. Затем количество атомов кислорода вычиталось, а содержания элементов вновь пересчитывались на 100%.

Составы шпинелей и оливина. В имеющихся составах шпинелей содержания FeO и Fe₂O₃ также пересчитывались на FeO (общее), которое вновь разделялось на окисную и закисную формы по стехиометрии, и с учетом ульвошпинелевого компонента. Затем, составы шпинелей и оливина пересчитывались на атомные проценты, содержание кислорода вычиталось и остаток катионов нормировался до 100%.

При дальнейшей обработке мы исходили из отсутствия структурного мотива в расплавах, статистически беспорядочного распределения катионов в кристаллах шпинели между октаэдрическими и тетраэдрическими позициями, и из того, что распределение Fe и Mg по позициям M₁ и M₂ в оливине близко к случайному. Исследовалась связь между совместными значениями Fe²⁺/Mg в разных фазах, аппроксимированная по методу наименьших квадратов (70) линейной зависимостью. По получаемому уравнению вида y = A + B x ; где: x - известное значение Fe²⁺/Mg в данной фазе; y - оцениваемое значение Fe²⁺/Mg в другой фазе; A и B - постоянные коэффициенты, вычислялось значение Fe²⁺/Mg в другой фазе. Сила связи оценивалась по величине коэффициента корреляции (R). По R оценивались также и влияние температуры (T), фугитивности кислорода (fO₂) и активности элементов в расплаве. Стандартная ошибка () определялась по уравнению:

( d_i ) / (n-2) ), где d_i = y_i (истинное) - y_i (расчетное)

Определялось также и значение абсолютной ошибки ( ). Распределение стандартной ошибки практически во всех случаях соответствовало гауссовскому, что позволило отбросить точки, для которых квадратичное отклонение было > 3. Их количество составляло ~ 4-7% от общего числа точек. Для оставшегося массива процедура расчетов вновь повторялась; эти результаты приведены в тексте. Все расчеты проводились с помощью программы "FASTVIEW" (автор - Ананьев В.В.).

Используемая форма представлений составов (атомные количества катионов) кажется более целесообразной в сравнении с традиционными. В таком виде данные точнее отражают как содержания элементов в каждой из фаз, так и отношения содержаний элементов между фазами, т.к. пропорциональны числу атомов, а не их весовой или молекулярной долям, что существенно для легких элементов типа Na и элементов с валентностями 2. Эта форма представления составов привычна для восприятия из-за небольших отличий цифровых значений от окисных весовых процентов, в отличие от атомных количеств с участием кислорода. Расчитываемое по валентностям количество кислорода является предельным для данной породы или расплава того же состава. В силикатных стеклах (расплавах) выделяются 3 разновидности кислорода: мостиковый (0⁰), немостиковый (0^1-) и свободный (0^2-) (9). Концентрации каждой из этих форм кислорода зависят от состава расплава, его структуры, состава флюида, температуры (T) и давления (P). Концентрация мостикового кислорода должна быть меньше предельной из-за существования наряду со связью Si-O-Si связи Si-Si "кислородная вакансия" (1) хотя не ясно, насколько велика их доля и как они зависят от вышеперечисленных параметров. Растворенная в расплаве Н₂О (на примере альбитовой системы), по данным М.Б.Эпельбаума (76), влияет на соотношение мостиковой и немостиковой форм кислорода. По расчетам (43) такие флюидные компоненты, как Н₂О и Н₂, в процессе дегазации базальтовых расплавов увеличивают фугитивность кислорода, что должно, вероятно, сказываться и на концентрации различных форм кислорода в расплаве. F и Cl, постоянные участники магматического процесса, в расплавах образуют ионные группировки с катионами металлов (2,33), т.е. выполняют роль "свободного" кислорода. Эти наблюдения позволяют рассматривать рассчитываемое по валентностям содержание кислорода в предполагаемом расплаве как максимально возможное и без ущерба исключать из данных по составам. Кроме того, операция вычитания рассчитанного количества кислорода (практически постоянная величина: 59-62) и нормирование до 100% остатка с дальнейшим поиском корреляции между процентными величинами численно соответствует (62) нормированию на постоянную величину, кратную содержанию кислорода, и поиску истинного коэффициента корреляции между процентными величинами. В перспективе желательно научиться рассчитывать истинные концентрации трех форм кислорода в расплавах и знать, как связана фугитивность кислорода с концентрациями их в расплаве, и как влияет концентрация каждой из форм кислорода в расплаве и фугитивность кислорода в целом на его вязкость и т.д.

Полученные результаты.

Для исследования поведения Fe²⁺/Mg в системе расплав (m) - шпинель (SP) - оливин (OL) был использован массив, сформированный из 92 точек. Выбраковано было 10 точек (по 5 для каждой из зависимостей: расплав - оливин и расплав - шпинель). Для дальнейших расчетов было использовано 82 точки, с диапазоном Fe²⁺/Mg отношений для расплава: 0,14-2,02. Были получены следующие результаты:

расплав-шпинель

(Fe²⁺/Mg)_m = 0,64 ((Fe²⁺- Ti ) / Mg )_sp + 0,13; R - 0,97; - 0,078; - 0,104.

Графически зависимость представлена на рис.1.

расплав-оливин

(Fe²⁺/Mg)_m = 2,71 ^. (Fe² / Mg)_ol + 0,05; R - 0,98; - 0,73; - 0,096.

Графически зависимость представлена на рис.1.

Обычно для системы расплав-оливин распределение Fe²⁺/Mg представляют в таком виде:

К_D = (Fe²⁺/Mg)_ol / (Fe²/Mg)_m K_D среднее - 0,3508; - 0,045; - 0,056;

интервал значений К_D: 0,26-0,5. K_D (Fe²/Mg) имеет одно и то же численное значение для молекулярных и атомных количеств содержаний элементов.

Численные значения коэффициента К_D пропорциональны содержаниям ряда элементов в расплаве: К_D-Si (R-0,8); K_D-Na (R-0,72); K_D-K(R-0,64). Связь между температурой и К_D имеет значение R=0,38. Невысокое (R-0,4) значение коэффициента корреляции имеет и связь отношения (Fe²⁺/Mg)_sp/(Fe²⁺/Mg)_m c температурой. Расчет (Fe²/Mg)_m по оливину может быть представлен и в такой форме:

(Fe²⁺/Mg)_m = {(Fe²⁺/Mg)_ol / 0,35}^. 0,95 + 0,05; R-0,98; -0,073; - 0,095.

Графически зависимость представлена на рис. 2.

Оливин - шпинель в расплаве.

(Fe²⁺/Mg)_ol = 0,23 ((Fe²⁺- Ti) / Mg)_sp+ 0,04; R-0,96; -0,036; d -0,048.

Графически зависимость представлена на рис.3

Шпинель - оливин в расплаве.

((Fe²⁺-Ti) / Mg)_sp= 4,06 ^. (Fe²⁺ /Mg)_ol-0,08; R-0,96; -0,147; -0,199.

Графически зависимость представлена на рис. 4.

Зависимость ((Fe²⁺- Ti) / Mg)_sp / (Fe²⁺/ Mg)_ol ( и наоборот) - Т (температура) имеет R-0,52 (интервал температур: 1100-1500^oС). Полученные результаты, на достаточно большой по объему (82 точки) выборке, демонстрируют хорошую корреляцию (R-0,96) Fe²/Mg отношений в OL-SP (или SP-OL) ассоциациях, находящихся в равновесии с расплавом при атмосферном давлении.

Каждая из зависимостей (m-SP; m-OL) может быть представлена в более точном виде, полученными по более объемным выборкам с учетом Р и Т. Это дает уточненные зависимости в распределении Fe²⁺/Mg для ассоциации OL-SP (или SP-OL) сосуществующих с расплавом.

Расплав (m) - шпинель (SP). Используются два варианта уравнения: "равновесное" - 126 точек (1); "неравновесное" - 208 точек (2) (55). Для разделения на "равновесные" и "неравновесные" составы фаз были использованы следующие характеристики условий опытов при атмосферном давлении: температура и продолжительность эксперимента (2).

(Fe²⁺/Mg)_m = 0,75 ((Fe²⁺- Ti) / Mg)_sp+ 0,04 R - 0,98; -0,034; - 0,044. (1)

(Fe²⁺/Mg)_m = 0,65 ((Fe²⁺- Ti) / Mg)_sp+ 0,1 R - 0,98; -0,05; - 0,065. (2)

Различия в коэффициентах уравнений (1) и (2) связаны с влиянием железистости расплава и "неравновесности" на вхождение Fe²⁺ и Ti в шпинель при ее кристаллизации (55). В интервале значений 1 атм. - 1,5 ГПа влияние давления незаметно (55). Для распределения ((Fe²⁺- Ti) / Mg)_sp / (Fe²⁺/Mg)_m ~ T ^oC в интервале 1100-1500^oC R= 0,2 , что свидетельствует о практическом отсутствии зависимости этого отношения от температуры (55).

Расплав-оливин. По результатам обработки имеющейся в нашем распоряжении выборки (82 точки, P=1 атм), можно утверждать следующее: (Fe²⁺/Mg)_ol / (Fe²/Mg)_m отношение (К_D) лежит в интервале значений: 0,26-0,5; среднее: 0,35; -0,045; -0,056. Зависимость этого К_D от температуры практически отсутствует (R-0,38). Вместе с тем существуют достаточно сильные корреляционные связи этого отношения с рядом элементов в расплаве: Si-К_D (0,8); Na-К_D (0,7); K- К_D (0,6).

Относительно истинного значения К_D и его зависимости от состава расплава, Т, Р пока единого мнения нет. В работах (78,89,96) приводятся значения К_D: 0,3-0,33 + 0,03. Значения К_D слабо зависят от температуры и давления. По данным (81,95), значения К_D лежат в интервале 0,25-0,38, и зависят от концентрации SiO₂ в расплаве и давления. Зависимость К_D от давления (до значений Р ~ З ГПа) можно практически не учитывать (95). По данным (101) в интервале давлений 10⁵ ПА до 1,5 ГПа значение К_D увеличивается на 0,03 относительно 0,3 (при 10⁵ ПА), т.е. находится в пределах ошибки.

На основании приведенных выше данных можно утверждать, что в системе основной (ультраосновной) расплав - шпинель - оливин распределение Fe²⁺/Mg между этими фазами практически не зависит от температуры и давления до 1,5 ГПа. Зная закон распределения Fe²⁺/Mg отношения для OL-SP (SP-OL) ассоциации, находящейся в равновесии с расплавом, можно различать эту ассоциацию магматического генезиса от метаморфического.

Более точные значения коэффициентов в уравнениях, связывающих Fe²⁺/Mg в OL-SP (SP-OL) парах, кристаллизовавшихся из магматического расплава, можно получить, объединив приведенные выше результаты. Для оливина использованы два значения К_D: 0,3 и 0,33, т.к. нормирование по SiO₂ пока не проведено. Для шпинелей используются уравнения (1) и (2), "неравновесное" уравнение (2) целесообразно использовать для случая предполагаемых высокожелезистых расплавов с Fe²⁺/Mg > 1,4. В итоге, для каждого из случаев и OL-SP (расчет состава оливина по шпинели) и SP-OL (расчет состава шпинели по оливину) получено по 4 представленных ниже уравнения:

OL-SP

SP-OL

Обсуждение результатов.

Полученные в виде 8 уравнений цифровые зависимости позволяют рассчитывать Fe²⁺/Mg отношение в оливине по "равновесной" с составом данного кристалла шпинели, или, наоборот, по известному составу оливина - (Fe²-Ti)/Mg в шпинели. Расчеты по уравнениям для случаев К_D =0,3 и К_D =0,33 в каждой из 4 пар дают различия меньше стандартного отклонения. Уравнения (5,6,9,10) - условно "неравновесные" целесообразно использовать для случая образования OL-SP пар из высокожелезистых расплавов с Fe²⁺/Mg>1,4, а также для кристаллов основной массы эффузивов, эволюционировавших по феннеровскому пути. Ко всем уравнениям приводятся численные значения абсолютной ошибки () и стандартного отклонения (), что позволяет с определенной долей вероятности ( - 70%, 2 - 95% (70)) судить о расплавном, или ином генезисе составов OL-SP пар в породе.

Влияние давления до 1,5 ГПа, температуры (1100-1500^oC) на распределение Fe²⁺/Mg отношения в системах расплав - шпинель и расплав - оливин не сказываются на численных значениях коэффициентов в уравнениях (55). Соответственно и в уравнениях, связывающих оливин - шпинелевые пары, "равновесные" с основным-ультраосновным расплавом. T и P, в указанных интервалах, не влияют на распределение Fe²⁺/Mg отношения.

Все приведенные выше уравнения могут быть использованы для кристаллов шпинели, содержащих менее 6% вес. TiO₂ и менее 50% вес. Al₂O₃ (55). В высокоглиноземистые шпинели (Al₂O₃ более 50% вес.) магний из расплава входит более охотно по сравнению с высокохромистыми разностями. Это сказывается на числовых значениях уравнений распределения Fe²⁺/Mg. Для случая расплав - высокоглиноземистая шпинель к настоящему времени удалось получить только оценочные значения коэффициентов (55).

Разбалансировка Fe²⁺/Mg отношения в системе OL-SP (SP-OL) относительно равновесного с расплавом (т.е. не отвечающая одному или нескольким уравнениям) происходит на всех этапах становления магматической породы: охлаждение, нагрев, метаморфизм. Причем это рассогласование может быть зафиксировано численно путем учета величины, и отчетливо проявляется из-за различного характера поведения Fe²⁺ и Mg в оливине и шпинели при нагреве, остывании. Это было выявлено экспериментально (84) и подтверждено на природных объектах (87). В работе (55) рассмотрены различные варианты этого явления и показано, что оно, вероятно, фиксируется и при метаморфических процессах.

Выявленные цифровые зависимости получены по экспериментальным данным, для которых понятия одновременность образования и равновесность составов OL-SP пар в расплаве однозначно связаны в подавляющем большинстве случаев, существование же отклонений может служить признаком некондиционности экспериментальных результатов (55). Иначе на природных объектах. В горной породе одновременность или разновременность образования кристаллов шпинели и оливина может быть установлена по их структурным взаимоотношениям. Если же характер срастания установить не удается, то вопрос остается открытым. Равновесность же или неравновесность, с учетом зональности, мозаичности и т.д. в составах кристаллов может быть уверенно зафиксирована с помощью полученных уравнений. Эта оценка, соотнесенная со структурными особенностями породы, может быть очень полезной, и дать генетическую информацию.

Приложение полученных результатов.

Для демонстрации возможностей применения полученных цифровых зависимостей, связывающих Fe²⁺/Mg в OL-SP (SP-OL) парах, равновесных с ультраосновным - основным расплавом, были использованы данные по составам OL-SP природных ассоциаций и содержащих их пород (56). Были рассмотрены достаточно известные по литературным источникам геологические объекты из двух фациальных групп, включающих эффузивные образования (лавовые потоки различного возраста, состава, из различных геолого-структурных обстановок) и интрузивные тела (гипербазиты современной океанической коры и офиолитовых комплексов, расслоенные интрузивные комплексы и т.д.). Для эффузивов, имеющих кайнотипный облик, равновесность вкрапленников оливина с включенными в них кристаллами шпинели явление характерное, тогда как в палеотипных разностях, например, в докембрийских коматиитах, составы разбалансированы. Неравновесность OL - SP ассоциации c расплавом может возникать и при гибридизме, примером, по-видимому, могут служить переуравновешенные шпинели в предположительно плейстоценовых пикритобазальтах Авачинского вулкана (53). Для интрузивных пород базит-гипербазитового ряда составы подобных OL-SP пар, как правило, перенормированы относительно равновесных с расплавом. Использование различных минеральных равновесий, отградуированных как геотермометры в солидусной области, для интрузивных пород этого ряда дает большой разброс численных значений в разных минеральных парах, и, вероятно, фиксирует температуры окончания обменных реакций (4,72). Это свидетельствует о метаморфических преобразованиях большей части минералов, слагающих интрузивные тела, хотя в некоторых случаях, благодаря мозаичному равновесию, могут сохраняться OL-SP пары, не лишенные первичных признаков кумулятивного или реститового происхождения (56). Fe²/Mg отношения в Ol - Sp парах, включенных в платиноиды из ряда ультраосновных массивов, свидетельствуют о возможности формирования вмещающих платиноидов как в ликвидусной, так и в солидусной температурных областях становления пород (56,57). Распределение Fe²/Mg отношения (или содержания магния) и хрома в шпинелях из платиносодержащих ультраосновных массивов может помочь как в выявлении наиболее перспективных из них, так и отдельных блоков в каждом массиве (53,54). Это подтверждается недавними ревизионными исследованиями платиноносных дунитов Нижне-Тагильского массива (15). Очень интересным объектом для приложения полученных зависимостей является оливин-шпинелевая ассоциация из включений в кристаллах алмазов. Вероятно наиболее глубинной (3-7 ГПа), соответствующей "законсервированным" условиям верхней мантии, можно считать перидотитовую ассоциацию. По данным (27) в алмазах из 4-х трубок Якутской кимберлитовой провинции эта ассоциация встречается в 4-10% кристаллов, тогда как самостоятельные оливин и шпинель встречены соответственно в 40-60% и 26-44% кристаллов. Мы сделали попытку оценить равновесность этой ассоциации с расплавом. Для этого использовались уравнения (3,4) в двух вариантах. В первом случае поправка на влияние давления (96,101) учитывалась только для оливина, во втором - вводилась поправка и для шпинели, равная поправке для оливина (т.к. данных по шпинели нет). В качестве первичных данных были использованы составы шпинелей и оливинов из алмазов Якутской кимберлитовой провинции. Мы нашли в литературе ~ 100 опубликованных анализов шпинелей из алмазов этой провинции. Обычно это хромиты, содержащие в среднем ~ 11-15 % MgO. Расчеты были выполнены для 2-х парных анализов SP-OL. В качестве первой пары было выбрано срастание наиболее магнезиальной шпинели (14,9% MgO) с оливином (MgO 52,2%) (68). В качестве второй пары были взяты наиболее магнезиальная шпинель (16,4% MgO) (66) и высокомагнезиальный оливин (52,8% вес MgO) (68) среди опубликованных составов шпинели и оливина из алмазов этой провинции. Оба варианта расчетов показали, что эти пары не являются равновесными с расплавом, причем с увеличением давления расхождение расчетного и реального Fe²⁺ /Mg для оливина увеличивается. Эти расчеты до некоторой степени условны (неизвестна реальная численная поправка для шпинели), но позволяют с определенной долей осторожности присоединиться к точке зрения (21), что алмаз не кристаллизуется в магматическом расплаве. В пользу этой гипотезы свидетельствуют данные о находках в образцах алмазоносных эклогитов линзовидных поликристаллических сростков алмазов, напоминающих друзы (57). По мнению (72), это свидетельствует о флюидном переносе углерода в мантии и переотложению в форме алмаза, что укладывается в рамки представлений о немагматическом генезисе.

Наиболее магнезиальными (MgO - 26,8% вес.) шпинелями из кимберлитов, вероятно, можно считать собственно шпинели (Al2O3 - 66,8% вес.) из алькремитов трубки "Удачная" (58). Судя по данным обзора (46), для высокоглиноземистых шпинелей мантийных пород типичны более высокие содержания MgO относительно хромистых разностей шпинели в мантийных породах. Это совпадает с выявленной нами тенденцией более охотного вхождения магния (относительно хромистых разностей) в состав высокоглиноземистой шпинели при росте ее или равновесии с расплавом. Но это не исключает образования подобного рода шпинелей в солидусных условиях.

Забегая вперед, обратим внимание еще на одно условие, возможно необходимое для роста кимберлитовых алмазов в твердой среде. Ниже приводятся данные по некимберлитовым алмазам, которые позволяют предположить, что, при определенных условиях причиной твердофазного перехода графит - алмаз могут быть сдвиговые деформации. Геологическая история площадей развития кимберлитового вулканизма на Сибирской платформе свидетельствует о его генетической связи с конвективными течениями в верхней мантии (47). Можно допустить, что собственно перемещения происходят вдоль ограниченных (10-100 км) латеральных и вертикальных зон, для которых характерны максимальные градиенты скоростей движения и наибольшие сдвиговые деформации. Многие глубинные (мантийные) ксенолиты перидотитовой и эклогитовой серий в кимберлитах как в целом по миру, так и на Сибирской платформе, несут отчетливые следы твердофазных деформаций. Они запечатлены не только в минералах, включая алмазы, но и выражены в структурных особенностях пород, слагающих ксенолиты (21,69). Можно предположить, оставляя в стороне происхождение графита, что преобразование его в алмаз происходит под действием движений в таких локальных зонах, можно допустить и то, что до определенной величины сдвиговых деформаций рост кристаллов алмаза угнетен, поэтому не все ксенолиты и кимберлиты алмазоносны. Кроме того, следы сдвиговых деформаций могут затушевываться последующими преобразованиями породы в мантийных условиях.

Возможность роста кристаллов алмаза в твердофазной среде и P,T условия такого роста представляют значительный интерес в связи с находками алмазов в некимберлитовых породах (30). Для алмазов из кимберлитов их высокобарический генезис (>3 ГПа) практически не вызывает споров, и разногласия в большинстве своем относятся к среде роста алмаза (расплав, твердофазная среда или их смесь) и к среде питания для роста кристаллов алмаза (расплав, флюид). Для алмазов из некимберлитовых пород (альпинотипные гипербазиты, высокобарические ультрамафиты и метаморфиды) основные разногласия касаются значений давления при образовании алмазов в этих породах. Для ряда массивов альпинотипных гипербазитов (28,30) приводятся значения давления их образования ~6,5 ГПа, оцененные по составам ортопироксенов из этих массивов. Накопилось достаточное количество данных, свидетельствующих о том, что составы минералов в этих массивах претерпели существенные изменения в процессе становления этих массивов. Поэтому условия, оцениваемые по составам этих минералов, могут быть неидентичными условиям раннего генезиса этих пород. Кроме того, коэсит, магнезиальный гранат (показатели высокобарического генезиса) не типичны для альпинотипных гипербазитов. Для высокобарических ультрамафитов, характерных для "корневых зон" складчатых поясов со сложной тектоникой и интенсивным метаморфизмом (39), вариации значений давления в процессе формирования массивов могут соответствовать как условиям для роста кристаллов алмаза, так и их графитизации. В гранатовых клинопироксенитах плутона Бени-Бушера (Марокко) найдены параморфозы графита по алмазу, составляющие местами до 15% объема породы (65). Причем внутри подобных высокобарических массивов, на примере массива Ronda (Испания) (86), можно на расстоянии 1,5 км наблюдать смену шпинель-гранатовых перидотитов на плагиоклазовые. Это дает основание считать перспективными на алмазоносность с глубиной выходы на дневную поверхность и плагиоклазовых разностей перидотитов из корневых зон складчатых областей.

Наиболее разноречивые генетические построения (мантийное или коровое образование) вызвали находки алмазов в метаморфических комплексах, примером которых могут служить месторождение Кумды-Коль и участок Барчинский, находящиеся на северном склоне Зерендинского гранитогнейсового купола в Кокчетавском срединном массиве (Казахстан) (17,20,23,34,42,45,51,64). Не вдаваясь во всю аргументацию в пользу возможности роста кристаллов алмаза в твердофазной среде в земной коре, мы хотим обратить внимание на некоторые факты в пользу такого генезиса, позволяющие предположить, что образование кристаллов алмаза могло происходить в моменты тектонических событий - палеоземлетрясений, в их гипоцентральной части (очаговой области). На этом месторождении алмазоносность приурочена к узкой, шириной от 45 до 250 м линейной зоне, внутри которой наибольшие концентрации алмазов наблюдаются в тектонических зонах мощностью от 5 до 20 м, представленных графитизированными высокотемпературными метасоматитами с бластокластическими структурами, образовавшимися по гнейсам, кальцифирам и эклогитам (23,52). Алмазы мелкие, представлены как отдельными кристаллами, так и их сростками, размером 1-400 мк, мелкие кристаллы (1-30 мк) в большинстве своем включены в кристаллы граната и, в меньшей степени, в зерна пироксена, амфибола, биотита, полевого шпата и в явно вторичные минералы (22). Описаны находки алмазов (15-100 мкм) в гранатах из разгнейсованной оболочки гранитного блока размером ~1,5м в двуслюдяных гнейсах (18). К числу особенностей кристаллов алмаза можно отнести их неоднородность, наличие или отсутствие графитовой рубашки вокруг кристаллов алмаза в одном образце или в одном кристалле граната (74). Р-Т условия многоэтапного формирования и существования породообразующей пироксен-гранат-плагиоклаз-кварцевой ассоциации лежат в следующих рамках: давление ~ 0,8-2,1 ГПа, температура ~ 600-1060^оС (51). Условия формирования минералов-включений (алмаз, коэсит, сфен, кианит, рутил и т.д.) в гранатах и цирконах из высокотемпературных метасоматитов, сложенных вышеуказанной минеральной ассоциацией, оцениваются как превышающие 4 ГПА и 850^оС (51).

Находки коэсита (67) не являются бесспорным аргументом в пользу мантийного генезиса минералов-включений. В экспериментах (82) в условиях одномерного сжатия при давлениях 0,5-2 ГПа, температурах 450-900^оС и скоростях деформаций 10^-4сек^-1 происходил метастабильный рост коэсита. Подобные скорости относятся к разряду очень быстрых геологичесикх деформаций (35) и характерны для локальных движений (участки < 10 км) непосредственно перед землетрясением в его эпицентральной зоне (48) в таких высокоподвижных областях, как Камчатка, Япония, Американские Кордильеры, относимых к зонам с 8-10-балльными землетрясениями (19). Разрывы сплошности для большинства пород земной коры происходят при деформациях, имеющих численную величину n^.(10^{-4 -} 10^-5) (61). Предельное значение напряжения, вызывающего различные виды разрушений в горных породах (скалывание, отрыв и т.д.), <0,3 Гпа при нормальных условиях; оно численно увеличивается с ростом всестороннего сжатия и уменьшается, приобретая характер пластично-хрупких деформаций, с увеличением длительности действия разрушающих напряжений и увеличением температуры (19). В настоящее время зависимость, связывающая воедино эти параметры, отсутствует. В работе (32) приводятся графические зависимости хрупко-пластичного перехода породы при одноосном сжатии, в зависимости от давления обжима, из которых следует, что такой переход происходит, когда напряжение сжатия в 2-3 раза превосходит давление обжима. Допуская понижающее влияние длительности воздействия тектонических напряжений и температуры в условиях нижних горизонтов земной коры при становлении гранито-гнейсовых формаций, можно предположить, что это отношение < 2. Минералогическими барометрами это влияние предположительно фиксируется. В работе (6) приводятся результаты исследований особенностей составов гранатов и оцениваются Р-Т условия формирования гранат - биотит-плагиоклазовой ассоциации из деформированных, в различной степени, слюдистых гнейсов, приуроченных к области сочленения Кольского и Беломорского геоблоков. "Зона главного разлома", имеющая мощность 100-150 м и прослеженная на расстоянии ~ 2 км, состоит из участков со слабо деформированными породами и зон с неоднократными, интенсивными деформациями, сложенными бластомилонитами. Авторами (6) выявлена прямая корреляционная зависимость между содержаниями кальция в гранатах и степенью дислоцированности вмещающих пород при отсутствии зависимости кальциевости гранатов от валового химического состава пород. В пластически слабо деформированных породах содержание гроссулярового минала 9-12%, и рассчитанные значения давления лежат в интервале 0,6-0,8 ГПа, что соответствует (6) литостатическому давлению (давлению обжима). Наиболее кальциевые гранаты (30-34% гроссуляра) и аномально высокие оцениваемые значения давления