Плюмовая и спрединговая ассоциации базальтов и их отражение в петрологических и геофизических параметрах северной части Срединно-Атлантического хребта

Л. В. Дмитриев

Институт геохимии и аналитической химии им. В. И. Вернадского РАН

С. Ю. Соколов

Геологический Институт РАН

В. Г. Мелсон, Т. О'Хирн

Смитсонианский Институт, США

Содержание

Аннотация

Для разработки количественных геодинамических моделей, прежде всего, необходимы данные, свидетельствующие о взаимосвязи между разными геологическими процессами, участвующими в формировании исследуемого объекта. Понятно, что их взаимосвязь, и, следовательно, их иерархия могут быть установлены и охарактеризованы на примере геологически активных, развивающихся объектов, в пределах которых результаты действия того или другого процесса или их комбинаций еще не замаскированы последующими событиями. Одним из таких объектов являются срединно-океанические хребты (СОХ), точнее их осевые зоны. Геологические структуры этих зон создаются в настоящее время, а их образование тесно связано с современным магматизмом и новейшими тектоническими движениями. Результаты взаимодействия этих процессов записываются в конфигурации геофизических полей.

Характеристика взаимосвязи процессов разного уровня, приводящих к образованию срединно-океанических хребтов, находит отражение в многочисленных работах последнего десятилетия, выполненных в рамках главной идеи проекта InterRIDGE, где срединно-океанические хребты рассматриваются как уникальные "Природные лаборатории". К числу таких работ относятся, например, исследования Е. Г. Мирлина и Н. М. Сущевской [Мирлин, Сущевская, 1990, 1992], Д. Форсайта [Forsyth, 1993], Я. Ниу и Р. Хекиньяна [Niu and Hekinian, 1997], Г. Дика [Dick, 1989], Дж. Кэси и др. [Casey et al., 1996] и др. В работах Ч. Лангмюра и др. [Langmuir et al., 1992] была убедительно показана взаимосвязь между петрологическими параметрами базальтового магматизма срединно-океанических хребтов и толщиной коры, записанной в рельефе рифтовой долины. В качестве петрологических параметров ими были предложены величины Na ₍₈₎ (степень частичного плавления мантии и мощность образующейся базальтовой коры) и Fe ₍₈₎ (глубина отделения расплава от мантийного источника). Эти же параметры были успешно использованы Э. Бонатти [Bonatti, 1996a, 1996b] при петрологическом районировании Срединно-Атлантического хребта (САХ). Связь между петрологическими параметрами магматизма и поверхностью геоида для всего Мирового океана недавно была охарактеризована в работе [Lecroart et al., 1997]. Индикаторная роль петрологических параметров была показана на примере взаимосвязи между геодинамической сегментацией САХ и масштабом адиабатического подъема вещества мантии до уровня его плавления [Дмитриев, 1998]. Перечисленные работы стимулируют исследование взаимосвязи между формированием геологических структур срединных хребтов и особенностями мантийного магматизма на количественном уровне.

В данной работе на примере наиболее изученного северного сектора САХ показан опыт сопоставления петрологических параметров базальтового магматизма, аномалий силы тяжести, сейсмичности и поверхности геоида осевой зоны хребта с использованием всей доступной информации.

Фактический материал и построение профилей вдоль САХ

Петрологические и геохимические параметры магматизма ТОР (Толеиты Океанических Рифтов - альтернативный термин БСОХ или MORB, предложенный нами [Дмитриев и др., 1978]), определены по составам их закалочных стекол из Каталога Смитсониановского Института [Melson et al., 1977], существенно дополненного за последнее десятилетие нашими оригинальными данными и опубликованными материалами [Сущевская, Волокитина, 1995; Сущевская и др., 1992, 1995, 1996; Niu and Batiza, 1993]. Для данного исследования были выбраны образцы современных базальтов, драгированных со дна и бортов осевой рифтовой долины САХ в интервале 0-80^o с. ш. Выборка представительных анализов включает 4481 состав.

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

В работе использованы данные по гравитационным аномалиям в свободном воздухе с разрешением 2' по [Sandwell and Smith, 1997]. Эти данные ограничены 72^o с. ш., поэтому они были дополнены данными по [McAdoo and Laxon, 1997], область которых ограничена 82^o с. ш. Указанные наборы данных позволили построить единый профиль от экватора до 80^o с. ш. Аномалии Буге изначально рассчитывались по сетке этих данных для трехмерной модели, и уже после этого было получено сечение аномалий по профилю. Методика вычислений аномалий Буге по данным альтиметрии изложена в [Мазарович, Соколов, 1999]. Отметим лишь, что данные по рельефу были взяты по сетке ЕТОРО5 с разрешением 5' [Global Relief CD, NOAA, 1993]; данные гравиметрии были пересчитаны на ту же сетку с понижением разрешения от 2' до 5'; была принята средняя плотность океанической коры 2,75 г/см³ ; поправка за рельеф вычислялась в области с радиусом 166 км. Информация о поверхности геоида использовалась по сетке с разрешением 15' [Lemoine et al., 1996] с получением соответствующего профиля значений вдоль оси САХ.

Информация о сейсмических событиях с магнитудой M_b>4 была получена по каталогу CNSS за период с 1960 г. по июнь 1997 г. (Магнитуда M_b , широта, долгота, время, категория глубин 0-13 км, 13-35 км и более 35 км). Сейсмический момент SM[J10¹⁷] (джоули) вычислялся для событий разных глубинных категорий по формуле:

Lg (SM) = 17,1 + 1,3 (M_b-5).

Значения сейсмических моментов каждого события и общее число событий суммировались для отрезков САХ длиной в один дуговой градус и отображались в виде профилей соответствующих значений (включая магнитуду) вдоль оси широт. Общая выборка событий для осевой зоны САХ от экватора до 80^o с. ш. составляет 4116 записей.

Магматизм и геодинамика САХ (краткий очерк)

Наиболее принципиальные результаты петрологического исследования ТОР, полученные к настоящему времени и используемые в настоящей работе, сводятся к следующему.

1. В конце 70-х была открыта бимодальность состава ТОР, устойчивая как в Атлантике, так и в Тихом океане. Тогда же возникла проблема геодинамического контроля магматизма, так как разрыв состава базальтов не объясним с позиций экспериментальной петрологии. В начале исследований, выделенные две группы базальтов получили название ТОР-1 и ТОР-2 [Дмитриев и др., 1978], а положение области разрыва составов было эмпирически установлено путем расчета восьмикомпонентной дискриминанты "D" [Сущевская и др., 1983]. Позже, по мере накопления нового материала были выделены и другие группы. Согласно предложенной нами классификации, основанной на кластерном анализе вариаций состава закалочных стекол ТОР (более 5000 микрозондовых анализов) среди них выделяется шесть групп [Дмитриев и др., 1990]. Четыре из них - ТОР-1, ТОР-2, ТОР-Fe и TOP-FeTi составляют вместе более 80% всех продуктов мантийного магматизма в Атлантическом и Тихом океанах. Меньше 20% приходится на группы ТОР-Na и TOP-K (обогащенные этими элементами соответственно).⁺⁺

Соотношения разных групп ТОР для северной части САХ в координатах Fe₍₈₎ Na₍₈₎ [Klein and Langmuir, 1987; Langmuir et al., 1992] приведено на рис. 1.

2. Распространение перечисленных групп ТОР в Атлантике и в Тихом океанах показано в табл. 1. (Для Индийского океана таких оценок не проводилось из за малого объема фактического материала).

Было установлено, также, что скорость спрединга контолирует характер пространственного распределения ТОР разных групп [Дмитриев и др., 1990]. Так, по простиранию СОХ с высокой скоростью спрединга (Восточно-Тихоокеанское Поднятие, Галапагосский Центр спрединга) и с переменной скоростью (САХ к югу от экватора) все перечисленные группы ТОР распределены в беспорядке, мозаично. В то же время вдоль САХ к северу от экватора (низкая скорость спрединга) разные группы ТОР образуют протяженные провинции, ограниченные крупными трансформными разломами, разделяющими хребет на тектоно-магматические сегменты разного порядка, детально обоснованные с учетом материалов по магматизму, тектонике и геофизике, и показанные на рис. 2 [Дмитриев, 1998; Дмитриев и др., 1994].

3. При новой классификации ТОР бимодальность всей их совокупности сохранилась и для Атлантики и для Тихого океана. Для ТОР осевой зоны САХ от экватора до 80^o с. ш. по данным 4481 составов было рассчитано новое значение дискриминанты D1 = 2,36 SiO₂ + 15,61 TiO₂ + 1,35 Al₂ O₃ - 0,80 FeO - 3,17 MnO + 3,19 MgO + 2,34 CaO + 21,35 Na₂ O - 29,15 K₂ O. На рис. 3 показана гистограмма D1 с хорошо выраженной узкой зоной минимума при значении D1₀ = 255,5. При дальнейшем использовании этой эмпирической дискриминанты выяснилось, что она является функцией важных петрологических параметров - таких как температура и давление котектической кристаллизации ТОР, степени их дифференциации и др. В качестве примера на рис. 4 приведен график, иллюстрирующий связь между D1 и Na₍₈₎ На этом рисунке вблизи D1₀ хорошо видна граница между двумя ассоциациями групп ТОР, которые показаны разными символами:

3.1. Ассоциация ТОР суперплюмов - Азорского и Исландского [Dosso et al., 1999; Schilling, 1973, 1975, 1986], представленных ТОР-1, ТОР-Fe и ТОР-К (TOP-FeTi в северной части САХ не установлен). В эту ассоциацию входят также сильно ограниченные в пространстве проявления ТОР-К (иногда с ТОР-1), локализованные вблизи границ между сегментами САХ разного порядка.

3.2. Ассоциация ТОР низкоскоростного спрединга, представленная ТОР-2 и ТОР-Na в виде обширной провинции в интервале 0-34^o с. ш. и в пределах небольшого блока между 48-52^o с. ш., разделяющего провинции Азорского и Исландского суперплюмов.

Сравнение разных групп ТОР и их плюмовых и спрединговых ассоциаций по степени их геохимического обогащения показано на диаграмме D1-K₂O/TiO₂ (рис. 5). Здесь видно, что подавляющее большинство базальтов (кроме ТОР-К) относятся к геохимически нормальным толеитам - N-MORB - по [Wilson, 1989]. Небольшая их часть принадлежит к переходному T-MORB. ТОР-1 незначительно обогащен по сравнению с ТОР-2. Наиболее обедненными являются ТОР-Na, а самыми обогащенными - ТОР-К, представленные преимущественно Е-MORB. Предполагается, что ТОР-К связан с "эмбриональной стадией" формирования внутриплитовых (или плюмовых) базальтов [Bougault et al., 1988].

Таким образом, рис. 5 подчеркивает, что степень обогащения ТОР возрастает с глубиной их мантийного источника. Это свидетельствует о том, что в мантии исследуемого объекта, как источника расплавов, отсутствуют какие либо значительные локальные геохимические неоднородности, связанные, например, с рециклингом вещества континентальной коры.

Распределение ТОР разных групп вдоль простирания САХ показано на рис. 6.

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

4. В табл. 2 приведены составы родоначальных расплавов ТОР плюмовой и спрединговой ассоциаций и Р-Т условия их отделения от мантийного источника. Для ТОР-1 и ТОР-2 они были восстановлены путем прямого эксперимента (гомогенизация расплавных включений) и численного моделирования [Дмитриев и др., 1978; Sobolev and Dmitriev, 1989]. Для ТОР-К прямые эксперименты не проводились. Здесь в качестве наиболее вероятного родоначального расплава ТОР-К используются экспериментальные и расчетные данные [Гуренко и др., 1991; Соболев, Никогосян, 1994] по субщелочным базальтам Исландии, сходными по составу с ТОР-К. Прямое определение родоначального расплава ТОР-Na также не проводилось. Здесь мы приводим данные о реконструкции его состава и условий его отделения от мантии полученные Л. В. Данюшевским по составу наиболее магнезиальных разностей закалочных стекол таких базальтов с использованием усовершенствованной программы "Petrolog" [Danyushevsky, 1998]. Подобные оценки, выполненные с применение программы ТОЛЕМАГ [Арискин и др., 1991], дают заниженные значения давления и температуры (Р - 5-6  кбар, Т - 1250^oС) [Сущевская и др., 1995].

Численным моделированием было показано, что ТОР-Fe и TOP-FeTi являются продуктами внутрикамерной дифференциации ТОР-1 в условиях пониженных давления и температуры (1-2 кбар, 1100-1200^o С) [Дмитриев и др., 1990].

5. Были восстановлены значения давления, температуры и содержания воды, соответствующие котектической стадии кристаллизации ТОР (MgO < 8,5%) по методу [Danyushevsky et al., 1996] с ошибкой Р2 кбар, Т5^oС и Н₂О0,2%.

На рис. 7 в координатах давление (Р кбар) - температура (Т^oС) показано положение продуктов котектической кристаллизации всех групп ТОР для северной части САХ (0-80^o с. ш.) относительно линии солидуса мантии [Herzberg and Zhang, 1996] и адиабаты 3^o С/кбар, принятой для декомпрессии первичных расплавов [Nisbet, 1982]. Даже с учетом ошибки расчетов график дает представление об условиях котектической кристаллизации расплавов после их отделения от мантийного источника. Хорошо видно, что наклон поля котектической кристаллизации заметно отличается от наклона адиабаты 3^oС/кбар. Это должно быть связано с значительной потерей тепла в окружающую среду на пути доставки расплавов к поверхности и их кристаллизации. Обращает внимание расширение поля котектической кристаллизации с уменьшением давления, благодаря чему его нижняя граница становится круче верхней относительно адиабаты, что также, вероятно, сопровождается дополнительной потерей тепла. В особенности это заметно в области самых низких давлений на примере связи ТОР-1 - ТОР-Fe, что мы объясняем внутрикамерной дифференциацией ТОР-1 (см. выше). Не исключено, что расширение всего поля котектической кристаллизации ТОР с уменьшением давления связано с более благоприятными условиями формирования промежуточных камер в верхних горизонтах литосферы.

Наиболее важным представляется резкое различие смены Р-Т условий эволюции расплавов при сравнении "плюмовой" ассоциации (ТОР-1 - ТОР-Fe - TOP-K) и "спрединговой" ассоциации (ТОР-2 - TOP-Na). Прежде всего, видно, что поле кристаллизации плюмовой ассоциации смещено в область более низких значений давления и температуры по сравнению со спрединговой ассоциацией. Более наглядно это показано на графике D1 - давление котектической кристаллизации Р, кбар (рис. 8).

Все упомянутые различия показаны на примере сравнения условий кристаллизации по средним составам всех групп ТОР и их родоначальных расплавов (табл. 2 и 3, рис. 9). Здесь видно, что в случае "плюмовой" ассоциации декомпрессия не менее 15 кбар (от 15-20 кбар до 1-3 кбар) сопровождается снижением температуры расплавов от адиабаты 3^oС/кбар до средних составов примерно на 200^o С, а для "спрединговой" ассоциации эволюция расплавов происходит в интервале 4-5 кбар и менее 100^oС соответственно. Для наглядности на рис. 9 схематически показаны поля эволюции магматизма обеих ассоциаций. Не исключено, что эти различия играют существенную роль при формировании разного теплового режима литосферы в условиях развития плюмов и при медленном спрединге, не осложненном другими процессами глобального масштаба. В первом случае, очевидно, литосфера должна быть более прогретой и менее вязкой по сравнению со вторым случаем, который по расчетам [Bonatti, 1996a, 1996b] отвечает холодной, жесткой литосфере ("холодный пояс литосферы Атлантики").

6. По данным о Р-Т условиях образования родоначальных расплавов ТОР (табл. 2) было установлено, что при двухслойной конвекции по [Richter and McKenzie, 1981], длина траекторий адиабатического подъема мантии от геотермы до пересечения солидуса с образованием родоначальных расплавов разных групп ТОР резко увеличивается при переходе от ТОР-2 (~200 км) к ТОР-1 (~400 км) и далее к ТОР-К (~700 км) [Дмитриев, 1998]. Расчет этих траекторий иллюстрируется рис. 10, который дополнен здесь данными об условиях образования родоначального расплава TOP-Na и траекториями подъема мантии при однослойной конвекции [Трубицын, Рыков, 1998]. Длина этих траекторий по существу отражает скорость подъема вещества мантии при условии постоянства хода геотермы. Если максимально возможная скорость этого подъема, например по оценкам [Reid and Jackson, 1981], достигает 10 см в год (здесь ТОР-К в плюмовой ассоциации при двухслойной конвекции), то ее минимальное значение (ТОР-2 в спрединговой ассоциации при однослойной конвекции) должно быть на порядок меньше - около 1 см в год. Неравномерность апвеллинга мантии, выраженная в длине траекторий ее подъема и показанная в плоской проекции разреза мантии вдоль оси САХ по данным о распределении разных групп ТОР хорошо согласуется с границами тектоно-магматических провинций и рассматривается здесь как одна из главных причин геодинамической сегментации [Дмитриев, 1998]. На рис. 11 этот разрез воспроизведен для северной части САХ и показан в двух вариантах - для двухслойной (а) и однослойной (б) конвекции. На этом же рисунке схематически показано положение зон разуплотненной мантии под двумя суперплюмами по [Zhang and Tanimoto, 1992].

Рис. 11 наглядно иллюстрирует корреляцию между неравномерностью мантийного апвеллинга, положением зон разуплотнения мантии и размещением "спрединговых" и "плюмовых" ассоциаций ТОР. Такая корреляция служит наглядным подтверждением взаимодействия глубинных мантийных процессов и магматизма в геодинамическом аспекте. Вместе с тем, совмещение в ограниченном пространстве и в современный период по времени столь контрастного проявления скоростей подъема мантии приводит к сомнению относительно конвекции в глобальном масштабе под осевой зоной хребта [Дмитриев, 1998]. Сравнение траекторий для двухслойной (а) и однослойной (б) конвекции (см. также рис. 10) показывает, что в первом случае длина траекторий примерно в два раза больше по сравнению со вторым случаем. Отсюда следует, что если конвекция существует, то она скорее однослойная чем двухслойная, так как при однослойной конвекции различие в длине траекторий (скоростях подъема) выражено менее контрастно по сравнению с двухслойной.

То же самое следует и из распределения изотерм, восстановленных по траекториям подъема мантии для обоих случаев (рис. 12). При двухслойной конвекции контрастное изменение теплового поля под срединным хребтом начинается на глубине порядка 700-800 км (а), а при однослойной всего на 200-250 км (б).

Сопоставление геофизических и петрологических параметров

На рис. 13 попарно показано сопоставление петрологических и геофизических параметров вдаль оси САХ. Вертикальными линиями обозначены разломы по границам между сегментами САХ разного порядка. Регионы плюмовой ассоциации ТОР заштрихованы. На верхней паре (a) сопоставлены давление котектической кристаллизации ТОР (Р, кбар) и поверхность геоида, рассчитанная по [Lemoine et al., 1996]. Это сопоставление показывает, что взаимодействие процессов, участвующих в аккреции литосферы САХ и записанное в петрологических параметрах магматизма влияет на формирование гравитационного поля планеты. Холодной литосфере спрединговой ассоциации соответствуют отрицательные значения геоида, а плюмовой ассоциации - положительные. На геоид влияет даже узкий блок холодной литосферы, разделяющий Азорский и Исландский суперплюмы.

На средней паре (б) сопоставлены дискриминаната D1 и кривая аномалий в свободном воздухе, по существу отражающая рельеф [Sandwell and Smith, 1997]. В этой аномалии эффект контрастной границы плотности кора-вода не снят и она в основном прямопропорциональна рельефу. Здесь также обе кривые зеркально противоположны, и показывают, что положительным формам рельефа отвечает плюмовая ассоциация, а понижению рельефа соответствует спрединговая ассоциация.

На нижней паре (в) мощность базальтовой коры, рассчитанная по параметру Na₍₈₎ сопоставлена с аномалией Буге [Мазарович, Соколов, 1999]. По существу мощность соответствует продуктивности магматизма (потеря базальтовой составляющей мантии в единицу времени). Видно, что продуктивность магматизма в районах проявления плюмов значимо выше, чем в условиях спрединга.

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

Рис. 13 показывает, что коэффициент корреляции между параметрами силы тяжести и петрологическими параметрами достаточно высок.

На рис. 14 приведено два примера более детального сопоставления аномалии силы тяжести в свободном воздухе и дискриминанты D1. Здесь видно, что соответствие этих параметров выполняется с разрешением не менее половины дугового градуса, т. е. порядка 40-50 км.

На рис. 15 на том же профиле сравнивается распределение магнитуды землетрясений с значениями M_b>4 и их частоты на двух категориях глубины (0-15 км и 15-35 км) за последние 37 лет. Хорошо видно, что в условиях спрединговой ассоциации (левая часть профиля) магнитуда землетрясений повышена (преобладают светлые кружки), а частота их на обоих категориях глубины понижена, по сравнению с условиями плюмовой ассоциации (правая часть профиля, где преобладают темные значки). Это находит объяснение в том, что в условиях холодной, жесткой литосферы землетрясения происходят сравнительно редко, но обладают повышенной магнитудой, а в условиях разогретой, пластичной и сравнительно легко проницаемой литосферы плюмов большие напряжения не накапливаются и землетрясения происходят чаще, но с меньшей магнитудой. Такое объяснение не противоречит существующим сейсмологическим моделям [Болдырев, 1998] и является одним из признаков корреляции петрологических и геофизических параметров. Также следует отметить, что в условиях плюмовой ассоциации наблюдается явное преобладание частоты слабых мелкофокусных событий над более глубинными, что говорит о большей подвижности верхней части литосферы САХ в зоне плюмов чем в зоне спрединговой ассоциации. Отмеченный выше общий тренд распределения сейсмичности в зоне плюмов нарушается узкими "всплесками" событий с большой магнитудой в районах Азорских островов, Исландии и разлома Чарли Гиббс, и сопутствующими им "всплесками" числа глубинных событий. Азорские острова в меньшей степени иллюстрируют указанные "всплески" как по магнитуде, так и по числу событий и находятся как бы в промежуточном положении по рассматриваемым характеристикам.

На рис. 16 представлены кривые распределения суммарного сейсмического момента на отрезок САХ длиной в один дуговой градус: 1 - Общий сейсмический момент по всем глубинным категориям. Хорошо видно, что за исключением указанных выше "всплесков", суммарный момент более интенсивен в области спрединговой ассоциации. По-видимому, это может быть связано с общей тенденцией к сдвигу реализации сейсмической энергии в область больших магнитуд по сравнению с плюмовой зоной, в которой преобладает реализация сейсмической энергии слабыми порциями. Кроме того, большое количество слабых событий не поддается удаленной регистрации и выпадает из общей картины, занижая оценки суммарного момента в плюмовой зоне с преобладанием слабых событий. 2 и 3 - сейсмические моменты для глубинных категорий 13-35 км и 0-13 км соответственно. Эти профили в общих чертах наследуют закономерности поведения общего сейсмического момента, но хорошо заметно, что "всплески" в области Исландии и разлома Чарли Гиббс резко выделены в более глубинной категории и имеют большую магнитуду событий.

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

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

Согласно петрологической модели [Klein and Langmuir, 1987; Langmuir et al., 1992] рельеф осевой зоны СОХ соответствует мощности образующейся базальтовой коры, которая в свою очередь зависит от продуктивности мантийного магматизма, т. е. от отношения объема базальтовой составляющей мантии, выведенной в единицу вpемени на поверхность или в промежуточные магматические камеры, к объему исходного вещества мантии. Согласно этой детально разработанной модели, продуктивность магматизма является функцией степени частичного плавления мантии и проницаемости литосферы, а выявленные связи между петрологическими параметрами и глубиной осевой зоны СОХ отражают "глобальный тренд" магматизма океанов. В этом исследовании вопрос о дискретности состава базальтов игнорируется. Вероятно, предполагается, что Р-Т условия плавления мантии и эволюции расплава, от которых зависят вариации состава базальтов, меняются непрерывно.

Дискретность состава ТОР (бимодальность, нарушение непрерывности во всем интервале составов, существование кластеров одного ряда) необъяснимая в рамках экспериментальной петрологии, обсуждается в наших работах [Дмитриев, 1998; Дмитриев и др., 1990, 1994]. Предполагается, что эта дискретность является функцией изменения внешних геодинамических условий океанского магматизма. В настоящем исследовании на большом фактическом материале по составу ТОР и геофизическим полям северной части САХ представляется возможность более детального обсуждения этого вопроса.

На рис. 17 приведены гистограммы ряда ключевых петрологических и геофизических параметров исследуемого отрезка осевой зоны САХ. Здесь, прежде всего, видно, что дискретность состава ТОР (гистограммы для D1, Na₍₈₎ , Р, кбар) подчеркивается дискретностью аномалий силы тяжести в свободном воздухе в мГал. На этой гистограмме просматривается небольшой минимум при значении мГал>0 и хорошо виден обособленный максимум для значений мГал>50. На профиле осевой зоны САХ (например, рис. 13) эта дискретность выражена в том, что все отрицательные значения мГал приходятся на те участки САХ, где размещены ТОР спрединговой ассоциации, а положительные - на участки с плюмовой ассоциацией. Значения мГал>50 локализованы в наиболее поднятой части Азорского сегмента, где распространены глубинные ТОР-К, а также на отрезках хребтов Рейкьянес и Кольбенсей, примыкающих к Исландии. Значения мГал>80 (небольшой максимум в правой части гистограммы) непосредственно совпадают с выходами ТОР-К. Здесь же фиксируются и признаки наибольшей продуктивности магматизма (минимальные значения D1, Na₍₈₎ и максимальная мощность коры). Это соответствие подчеркивается также и тем, что максимальные значения этих параметров и минимальная мощность коры (локализация ТОР-Na) совпадают в пространстве с положением минимальных значений мГал<-40 , хотя на гистограмме это и не выражено, вероятно, из за малого распространения TOP-Na.

Гистограмма для D1 отчетливо разделяет плюмовую и спрединговую ассоциации (см. выше). Кроме того, здесь виден небольшой максимум при самых низких значениях D1, соответствующий ТОР-К.

Гистограммы для Na₍₈₎ как функции продуктивности мантийного магматизма, и для давления котектической кристаллизации ТОР также содержат признаки их дискретности, хотя и с менее четко выраженными минимумами. Возможно это связано с тем, что Na₍₈₎ рассчитывается всего по двум элементам и на более ограниченном по объему и информативности материале, так как расчет ведется для составов с содержанием MgO<8,5%. То же относится и к давлению котектической кристаллизации, которое рассчитывается с ошибкой порядка 2 кбар.

Все упомянутое дает основание для вывода о том, что продуктивность мантийного магматизма и диапазон давления котектической кристаллизации ТОР в северной части САХ носит отчетливо дискретный характер. При формировании плюмовых ассоциаций ТОР продуктивность магматизма значимо выше, а давление (глубина) котектической кристаллизации значимо ниже, чем при формировании спрединговых ассоциациий. Это подчеркивается дискретностью аномалий силы тяжести в свободном воздухе. Другие гистограммы, приведенные на рис. 17, не показывают четкого разделения ТОР на плюмовую и спрединговую ассоциации. На гистограмме для K₂ O/TiO₂ просматривается лишь небольшой минимум (при значениях около 0,3), отделяющий ТОР-К от остальных ТОР (см. рис. 5). Слабо выраженный минимум для Fe₍₈₎ (около 10), вероятно, отделяет группу ТОР-Fe (рис. 1), а гистограмма для Т^o С показывает, что интервалы температуры котектической кристаллизации ТОР обеих ассоциаций перекрывают друг друга.

На гистограмме для аномалии Буге (распределение силы тяжести в литосфере со снятым влиянием рельефа, образованного базальтовой корой) виден небольшой минимум при ее значениях около 150 мГал. Этот минимум отделяет Исландию с примыкающими участками хребтов Рейкьянес и Кольбенсей (рис. 13), как и для аномалии силы тяжести в свободном воздухе при значениях мГал>50. Отсутствие минимума на гистограмме, разделяющего ппюмовую и спрединговую ассоциации находит объяснение в том, что разная продуктивность магматизма, создающего базальтовую кору разной мощности в спрединговой и плюмовой зонах, особо подчеркивается аномалией силы тяжести в свободном воздухе без снятого влияния рельефа, так как продуктивность магматизма является сильным рельефообразующим фактором. Расчет аномалии Буге снимает влияние рельефа, увеличивая вклад внутренних неоднородностей верхней части литосферы в аномальном поле, и, таким образом, маскирует вклад базальтовой коры, имеющей меньшую плотность и мощность, чем те зоны литосферы, которые подчеркиваются после расчета аномалии. Контраст плотности кора - верхняя мантия меньше, чем кора - водный слой, что также объясняет отсутствие на гистограмме модальности, разделяющей ассоциации базальтов коры.

Связь между петрологическими параметрами и значениями геофизических полей хорошо согласуется с данными о землетрясениях. Как было показано на рис. 12 в спрединговых участках САХ землетрясения отличаются сравнительно высокой энергией, но происходят относительно редко по сравнению с плюмовыми участками, где их энергия заметно меньше, но происходят они чаще. Это хорошо согласуется с представлением о разной продуктивности магматизма в тех и других участках хребта. Можно предположить, что в спрединговых участках хребта низкая продуктивность магматизма и сравнительно небольшое различие между Р-Т условиями отделения родоначальных расплавов от мантийного источника и их котектической кристаллизации не способствуют прогреву литосферы, которая здесь сохраняет жесткость и остается мало проницаемой ("холодный пояс литосферы Атлантики" по [Bonatti, 1996a, 1996b; Bonatti et al., 1993]). В этих условиях вывод на поверхность продуктов магматизма может быть связан с растрескиванием мощной литосферы при достижении критических значений усилий ее растяжения. Судя по близости Р-Т условий котектической кристаллизации магнезиальных разностей ТОР этой ассоциации и соответствующих им родоначальных расплавов, трещины в литосфере могут достигать здесь глубины генерации этих расплавов. Кроме того, все ТОР этой ассоциации кристаллизуются в сравнительно узком интервале температуры, но при переменном давлении (полибарическое фракционирование). Все это может быть признаком того, что вывод расплава на поверхность происходит быстро (для высокомагнезиальных разностей - почти мгновенно). Это не способствует его кристаллизационной дифференциации на пути подъема с отдачей тепла в литосферу. В результате слабо раскристаллизованные расплавы выливаются на поверхность дна при большой температуре и закаливаются морской водой с образованием афировых разностей базальтов с относительно толстыми корками стекла. Последнее подтверждается материалами многочисленных публикаций и экспедиционных отчетов. Вполне вероятно, что при такой геодинамической обстановке сейсмические события будут происходить относительно редко, но с высокой энергией.

Не исключено, что максимальные по энергии землетрясения будут проходить в верхних горизонтах литосферы, т. е. там, где она менее разогрета и обладает повышенной прочностью. Эти соображения подтверждаются при анализе рис. 11.

С этим же согласуются особенности проявления наиболее глубинного магматизма плюмовой ассоциации - ТОР-К - в пределах "холодного пояса литосферы". Эти проявления здесь редки, сильно локализованы и сфокусированы вблизи тектонически подвижных зон глубокого заложения по границам между сегментами САХ разного порядка. Весьма показательно, что по данным сейсмической томографии в участках спрединговой ассоциации ТОР отсутствуют признаки разуплотнения мантии [Zhang and Tanimoto, 1992]. Особенно эффектным здесь выглядит узкий блок холодной литосферы на границе между двумя суперплюмами.

Сейсмические события в районах суперплюмов происходят чаще, но с меньшей энергией по сравнению с участками спрединговых ассоциаций. Это, по-видимому, связано с тем, что благодаря высокой скорости подъема нагретой мантии и повышенной продуктивности магматизма, литосфера здесь более прогрета, обладает пониженной вязкостью и сравнительно легко проницаема. Это согласуется и с данными о повышенных значениях здесь теплового потока, например, [Попова, Смирнов, 1985] и с материалами сейсмической томографии о наличии здесь обширных областей разуплотненной мантии. Расплавы здесь чаще поднимаются в верхние горизонты литосферы, а их кристаллизация и остывание происходит в широком интервале давления и температуры. В самых верхних горизонтах литосферы создаются благоприятные условия для аккумуляции магмы и ее дальнейшей кристаллизации в промежуточных камерах при постоянном и низком давлении с образованием ТОР-Fe. Существование таких камер недавно подтвердилось по результатам сейсмического эксперимента [Constable et al., 1997]. Последнее более характерно для Исландского суперплюма с его железистыми базальтами на хребтах Рейкьянес и Кольбенсей. В этом отношении Исландский суперплюм представляется более развитой стадией формирования плюмов по сравнению с Азорским.

В целом же геодинамическая обстановка формирования суперплюмов в Северной Атлантике (при самой низкой скорости спрединга) весьма напоминает таковую в Тихом океане (при самой высокой скорости спрединга), где также распространены ТОР-Fe как продукты внутрикамерной дифференциации при низком давлении. С одной стороны, это может служить подтверждением того, что скорость спрединга не влияет на условия формирования мантийных расплавов и их эволюцию, а лишь контролирует характер их распределения в пространстве [Дмитриев и др., 1990], а также на "глобальный тренд" магматизма океанов [Langmuir et al., 1992]. С другой стороны, это ставит вопрос о механизме формирования плюмов на хребтах с низкой скоростью спрединга.

Различие продуктивности магматизма обеих главных ассоциаций ТОР согласуется также с фактическим материалом глубинного сейсмического зондирования и с реальными данными опробования. Так, например, по данным [Горячев, Павлов, 1985] мощность коры в районе Исландского плюма достигает 30 км, а в "холодном поясе" она может вообще отсутствовать, как, например к югу от разлома Зеленого мыса (участки "сухого спрединга", [Bonatti, 1996a, 1996b; Cannat, 1993; Karson, 1991]). Продуктивность магматизма отражается также на составе остаточного вещества мантии: составы базальтовых выплавок комплементарны составу остаточных перидотитов по степени обогащения [Силантьев и др., 1999; Bonatti, 1996a, 1996; Silantyev et al., 1995]. Гидротермальные источники в северной части САХ сосредоточены в районе развития плюмов. Редкие гидротермальные поля в спрединговой ассоциации расположены только вблизи сильно локализованных проявлений ТОР-К, а в самих блоках "холодной литосферы" они отсутствуют [Батуев и др, 1995; Rona and Scott, 1993]. Это можно объяснить тем, что для формирования гидротермальных реакторов необходимы стабильные источники тепла связанные, например, с магматическим камерами, расположенными на небольшой глубине. Очевидно такие камеры в спрединговой ассоциации отсутствуют и достаточно распространены в плюмовой асоциации. Это может служить независимым подтверждением геодинамических различий между двумя главными ассоциациями ТОР.

Корреляция петрологических, геофизических и геодинамических параметров при развитии современного магматизма САХ служит надежным подтверждением идеи о взаимосвязи всех процессов любого уровня, участвующих в формировании литосферы. В свою очередь это служит реальной основой для выявления причин резкого различия условий эволюции магматизма плюмовых и спрединговых ассоциаций и их сосуществования в пространстве и времени. Для решения этих задач, очевидно, следует прежде всего обратить внимание на наиболее контрастное проявление тех или других петрологических и геодинамических параметров при сравнении обеих ассоциаций. Как показывает это исследование, наиболее контрастным представляется различие в скоростях подъема вещества мантии до пересечения солидуса и образования расплавов ТОР разного состава. В тех сегментах САХ, где преобладает спрединговая ассоциация ТОР скорость подъема мантии минимальна (вероятно в пределах 1-3 см в год) и может быть согласована с принятой скоростью медленного спрединга. В обширных регионах развития суперплюмов скорость подъема мантии резко возрастает и может быть в несколько раз больше, чем в спрединговых сегментах. Соответственно, если исходить из моделей мантийной конвекции, скорость спрединга в плюмовых регионах должна быть близкой к максимальной, что не согласуется с общепринятыми оценками для северной Атлантики. Противоречие моделей конвекции и неравномерности мантийного апвеллинга обсуждалось в работе [Дмитриев, 1998] и здесь мы к этому возвращаться не будем. Добавим лишь, что если осевая часть САХ в интервале от экватора до 80^o с. ш. в современный период располагается над восходящим потоком мантии, связанным с конвекцией, то исходя из количественных данных петрологии базальтового магматизма согласующихся с геофизическими полями, придется признать, что этот поток имеет чрезвычайно сложный характер (рис. 11, 12). Это должно стимулировать уточнение существующих моделей мантийной конвекции и оценок геотермального режима (например, [Anderson, 1998; Cordery and Phipps Morgan, 1993; Dobretsov and Kirdyashkin, 1993; Sparks and Parmentier, 1993]). Благодаря недавним исследованиям В. П. Трубицына появляется возможность решения этой проблемы. Согласно этим исследованиям конфигурация восходящего конвективного потока сильно зависит от его температуры. При ее повышении возможно его расслоение и образование локализованных мантийных струй весьма малого сечения - порядка 200 км и менее (устное сообщение). Такие оценки уже приближаются к нашим оценкам минимального поперечного сечения индивидуального мантийного диапира, производящего только один тип родоначального расплава ТОР [Дмитриев, 1998]. В любом случае представляется, что "средние" геотермы мало применимы для конкретных геодинамических реконструкций. В этих случаях необходимы детальные геотермические исследования на отдельных сегментах хребтов и на профилях вкрест их простирания.

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

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

Контрастные проявления всех перечисленных процессов и их результатов могут быть причиной дискретности состава ТОР, петрологических параметров их эволюции, продуктивности магматизма, мощности образующейся базальтовой коры и, как следствие, характера геофизических полей. Обсуждение вопроса о том, какие глобальные геологические (или планетарные ?) процессы приводят к столь контрастному проявлению мантийного апвеллинга и сосуществованию спрединговых и плюмовых геодинамических обстановок при аккреции океанской литосферы становится весьма актуальным (например, [Douglass and Schilling, 1999]). Обсуждение этой проблемы, составляющей предмет самостоятельного углубленного исследования в условиях явного недостатка полезной информации, выходит за рамки настоящей работы.

Заключение

На данном конкретном примере было показано, что любые геологические процессы, участвующие в аккреции северной части САХ, тесно связаны между собой, а их корреляция может быть охарактеризована на количественном уровне и с большим разрешением в пространстве. Последнее оказывается возможным в силу того, что объект исследования формируется в современный геологический период, когда результаты взаимодействия процессов еще не замаскированы будущими геологическими событиями.

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

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

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

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

Несомненно, связь между петрологическими и геофизическими параметрами осуществляется и в ходе формирования современных структур всей мировой системы срединно-океанических хребтов. Исследование этой связи представляет интерес, прежде всего, для решения проблемы неоднородности мантии в глобальном масштабе, проблемы взаимодействия конвекции и спрединга и их роли в формировании разных геодинамических обстановок в ходе аккреции литосферы.

Благодарности

Литература

Арискин А. А., Цехоня Т. И., Френкель М. Я., ЭВМ-барометрия и генетическая интерпретация базальтовых стекол Центральной Атлантики, Геохимия, (7), 1038-1047, 1991.

Батуев Б. Н., Кротов А. Г., Марков В. Ф., Краснов С. Г., Лисицын Е. Д., Черкашев Г. А., Новые гидротермальные поля в осевой зоне Срединно-Атлантического хребта (14^o 45' с. ш.), Доклады РАН, 343, (1), 75-79, 1995.

Болдырев С. А., Сейсмогеодинамика Срединно Атлантического хребта, 124 с., МГК, Москва, 1998.

Горячев Ю. В., Павлов С. П., Мощность земной коры, Геофизические характеристики земной коры Атлантического океана, с. 135-138, Недра, Ленинград, 1985.

Гуренко А. А., Соболев А. В., Кононкова Н. Н., Щелочные рифтогенные базальты Исландии: новые данные по петрологии, Геохимия, (9), 1262-1274, 1991.

Дмитриев Л. В., Магматизм дна океана и вопросы формирования его литосферы, Геология и геофизика, (7), 81-85, 1986.

Дмитриев Л. В., Вариации состава базальтов срединно-океанических хребтов как функция геодинамической обстановки их формирования, Петрология, 6, (4), 340-362, 1998.

Дмитриев Л. В., Соболев А. В., Сущевская Н. М., Первичный расплав океанических толеитов и состав верхней мантии океана, ДАН СССР, 240, (1), 177-180, 1978.

Дмитриев Л. В., Соболев А. В., Рейснер М. Г., Мелсон В. Дж., Петрохимические группы закалочных стекол ТОР (Толеиты Океанических Рифтов) и их распределение в Атлантическом и Тихом океанах, Кн. Магматизм и тектоника океана, Ред. Ю. М. Пущаровский, c. 43-108, Наука, Москва, 1990.

Дмитриев Л. В., Соболев А. В., Данюшевский Л. В., Ахметов К. С., Сегментация Срединно-Атлантического хребта и ее связь с давлением кристаллизации и содержанием воды в базальтовых расплавах, кн. Нелинейная геодинамика, Ред. Ю. М. Пущаровский, с. 26-54, Наука, Москва, 1994.

Данюшевский Л. В., Соболев А. В., Дмитриев Л. В., Ортопироксенсодержащие низкотитанистые толеиты - новый тип толеитов океанических рифтов, ДАН СССР, 292, (6), 1449-1452, 1987.

Мазарович А. О., Соколов С. Ю., Разломы Ангольской котловины, Российский журнал наук о Земле, 1, (3), 1999, (Электронная версия: http://eos.wdcb.ru/rjes/rje98014/rje98014.htm).

Мирлин Е. Г., Сущевская Н. М., Пространственно-временная неравномерность океаногенеза и ее структурно-петрологические следствия (на примере Атлантического океана), Магматизм и тектоника океана, Ред. Ю. М. Пущаровский, с. 108-122, Наука, Москва, 1990.

Мирлин Е. Г., Сущевская Н. М., Тектономагматическая сегментация срединно-океанских хребтов, мантийный апвеллинг, асимметрия Земли, Геотектоника, (2), 3-21, 1992.

Попова А. К., Смирнов Я. Б., Геотермическая активность срединно-океанического хребта, Геофизические характеристики земной коры Атлантического океана, с. 145-151, Недра, Ленинград, 1985.

Силантьев С. А., Базылев Б. А., Дмитриев Л. В., Плечова А. А., Петрологическая оценка геодинамических параметров мантийного диапиризма в Срединно-Атлантическом хребте, "Тектоника, геодинамика и процессы магматизма и метаморфизма", материалы совещания, т. 2, с. 123-126, ГЕОС, Москва, 1999.

Соболев А. В., Никогосян И. К., Петрология магматизма долгоживущих мантийных струй: Гавайские о-ва (Тихий океан) и о-в Реюньон (Индийский океан), Петрология, (2), 131-168, 1994.

Соболев А. В., Шимизу Н., Сверхобедненные расплавы и проницаемость океанической мантии, Доклады РАН, 326, (2), 354-360, 1992.

Сущевская Н. М., Волокитина Л. П., Тектоника и магматизм приэкваториальной зоны Срединно-Атлантического хребта, Океанология, 35, (2), 252-265, 1995.

Сущевская Н. М., Дмитриев Л. В., Соболев А. В., Петрохимический критерий классификации закалочных стекол океанских толеитов, ДАН СССР, 268, (6), 1475-1477, 1983.

Сущевская Н. М., Пейве А. А., Цехоня Т. И., Сколотнев С. Г., Кононкова Н. Н., Петрология и геохимия магматизма активных частей разломов Романш, Сан-Паулу и сопряженных с ними частей Срединно-Атлантического хребта, Геохимия, (5), 697-719, 1995.

Сущевская Н. М., Удинцев Г. Б., Цехоня Т. И., Кононкова Н. Н., Генезис и геохимические особенности толеитового магматизма Срединно-Атлантического хребта (3-5^o с. ш.), ДАН СССР, 322, (1), 147-154, 1992.

Сущевская Н. М., Цехоня Т. Е., Дубинин Е. П., Мирлин Е. Г., Кононкова Н. Н., Формирование океанской коры в системе срединно-океанических хребтов Индийского океана, Геохимия, (10), 936-975, 1996.

Трубицын В. П., Рыков В. В., Мантийная конвекция и глобальная тектоника Земли, Вестник ОГГГГ РАН, (электронная версия) (1(3))б, 1998.

Anderson D. L., The scales of mantle convection, Tectonophysics, 284, (Iss. 1-2), 1-17, 1998.

Bonatti E., Long-lived oceanic transform boundaries formed above mantle thermal minima, Geology, 24, 803-806, 1996a.

Bonatti E., Anomalous openung of the Equatorial Atlantic due to an equatorial mantle thermal minimum, Earth Planet. Sci. Lett., 143, 147-160, 1996b.

Bonatti E., Seyer M., and Sushevskaya N. M., A cold suboceanic mantle belt at the Earth's equator, Science, 261, 315-320, 1993.

Bougault H., Dmitriev L. V., Shilling J. G. et al., Mantle heterogeneity from trace elements MAR tripple junction near 14^o N, Earth and Planet. Sci. Lett., 88, 27-36, 1988.

Cannat M., Emplacement of mantle rocks in the sea floor at Mid-ocean Ridges, Journ. Jeoph. Res., 98, (B3), 4163-4172, 1993.

Casey J. F., Zhou H., Pavlenko E., Xia C., Golod V., and Bryan W. B., Comparison of Gravity, Bathymetry, Mantle Tomography and Geochemistry of Basalts along the Northern Mid-Atlantic Ridge, FARA-IR Mid-Atlantic Ridge Symposium, 19th-22nd June 1996, Reykjavik, Iceland, Cambridge Publications Home Page, p. 767.

Constable S., Sinha M., MacGregor L., Navin D., Peirce C., White A., and Heinson G., RAMESSES finds a magma chamber beneath a slow spreading ridge, InterRidge News, 6, 18-22, 1997.

Cordery M. J., Phipps Morgan J., Convection and melting at mid-ocean ridges, Journ. Geoph. Res., 98, 19,477-19,503, 1993.

Danyushevsky L. V., The effect of small amount of H ₂ O on fractionation of mis-ocean ridge magmas, Abstr. AGU Spring Meeting. Eos 79, (17/suppl), p. 375, 1998.

Danyushevsky L. V., Sobolev A. V., and Dmitriev L. V., Estimation of pressure of crystallization and H ₂ O content of MORB and BABB glasses: calibration of an empirical technique, Mineralogy and Petrology, 57, 185-204, 1996.

Dick H. J. B., Abyssal peridotites, very slow spreading ridges and ocean ridge magmatism, in: Magmatism in the Ocean Basins, A. D. Sounders and M. J. Norry (eds.), pp. 71-105, Geol. Soc. Spec. Publ., London, 1989.

Dmitriev L. V., Sobolev A. V., Reisner M. G., and Melson W. G., Quenched glasses of TOR: Petrochemical classification and distribution in Atlantic and Pacific oceans, 28th IGC, pp. 1399-1401, Washington D.C. USA, Abstracts, 1989.

Dmitriev L. V., Melson W. G., and Reisner M. G., The distribution of petrochemical types of the tholeiitic abyssal glasses, International geological-geophysical Atlas of the Atlantic ocean, p. 110, IOC, UNESCO, 1990.

Dobretsov N. L., Kirdyashkin A. G., Experimental modelling of two layered mantle convection, Ofoliti, 18, (1), 61-81, 1993.

Dosso L., Bougault H., Langmuir C., Bollinger C., Bonnier O., and Etoubleau J., The age and distribution of mantle heterogeneity along the Mid-Atlantic Ridge (31-41^o N), Earth Planet. Sci. Lett., 170, 269-286, 1999.

Douglass J., and Schilling J.G., Plume-ridge interactions of the Discovery and Shona mantle plumes with the southern Mid-Atlantic Ridge (40-45^o S), Journ. Geoph. Res., 104, (B2), 2941-2962, February 1999.

Forsyth D. W., Crustal Thickness and the Average Depth and Degree of Melting in Fractional Melting Models of Passive Flow Beneath Mid-Ocean Ridges, Journ. Geoph. Res., 98, (B9), 16,073-16,079, 1993.

General Bathimetric Chart of the Oceans, Digital Atlas, British Oceanographic Data Center, IOC, Proudman Oceanologic Laboratory, Bidston Observatory Birkenchead, 1997 (CD-ROM).

Global Relief CD, NOAA, 1993.

Herzberg C., and Zhang J., Melting experiments on anhydrous peridotite KLB-1; composition of magmas in the upper mantle and transitional zone, Journ. Geoph. Res., B 101, (4), 8271-8295, 1996.

Karson J. A., Seaflor spreading on the Mid-Atlantic Ridge: Implication for the structure of ophiolites and oceanic lithosphere produced in slow-spreading environments, Proceedings of the Simposium "Troodos 1978", Edited by J. Malpas, E. Moores, A. Panayiotou, C. Xenofontos, pp. 547-555, Nicosia, Cyprus, 1991.

Klein E. M., and Langmuir C. H., Global correlation of ocean ridge basalt chemistry with axial depth and crustal thickness, Journ. Geoph. Res., 92, 8089-8115, 1987.

Langmuir C. H., Klein E. M., and Plank T., Petrological systematics of Mid-Ocean Ridge basalts: Constraints of melt generation beneath ocean ridges, Mantle Flow and Melt Generation at Mid-Ocean Ridges, J. Phipps Morgan et al. (eds.), pp. 183-280, AGU, Washington DC, 1992.

Lecroart P., Albarede F., and Cazenave A.,. Correlation of Mid-Ocean Ridge Basalt chemistry with the geoid, Earth Planet. Sci. Lett., 153, 37-55, 1997.

Lemoine F. G. et al., The Development of the NASA GSFC and DMA Joint Geopotential Model, International Symposium on Gravity, Geoid and Marine Geodesy (GraGeoMar96), Univ. of Tokyo, Tokyo, Japan, September 30-October 5 1996, Geoid Undulation Grid from EGM96, NASA-NIMA, 1996.

McAdoo D. C., and Laxon S. W., 1997, Antarctic Tectonics: Constraints from a new ERS-1 Satellite Marine Gravity Field, Science, 276, (5312), 556-561, 25 April, 1997.

Melson W.G., and O'Hearn T., Basaltic glass erupted along the Mid-Atlantic Ridge between 0-37^o N: relationship between composition and latitude, Deep Drilling Results in the Atlantic Ocean Crust, Edited by M. Talwani, C. G. Harrison and D. E. Hayes, pp. 249-260, AGU, Washington, D.C., 1979.

Melson W. G., O'Hearn T., "Zero-age" variations in the composition of abyssal volcanic rocks along the axial zone of the Mid-Atlantic Ridge, The Geology of North America, V. M, The Western North Atlantic Region, The Geological Society of America, pp. 117-136, 1986.

Melson W. G., Byerly G. R., Nelen J. A., et al., Catalog of major element chemistry of abyssal volcanic glasses, Mineral Sci. Investigation. Smithsonian Contrib. Earth Sci., 19, 31-60, 1977.

NEIC (National Earthquake Informational Center), USA, 1997.

Niu Ya., Batiza R., Chemical Variation Trends at Fast and Slow Spreading Mid-Ocean Ridges, Journ. Geoph. Res., 98, (B5), 7887-7902, 1993.

Niu Ya., and Hekinian R., Spreading-rate dependence of the extent of mantle melting beneath ocean ridges, Nature, 385, 326-329, 23 January 1997.

Nisbet E. G., Tectonic setting and petrogenesis of komatiite, Komatiites, edited by G.A.a., pp. 501-520, Unwin, London, 1982.

Plechova A. A., High alumina-water content tholeiite - the new type of mid-ocean ridge basalt, 6th Zonenshain Conference on Plate Tectonic, 20 p., Moscow, Abstracts, 1998.

Reid I., and Jackson H. R., Oceanic spreading rate and crustal thickness, Marine Geoph. Res., 5, 165-172, 1981.

Richter F. M., and McKenzie D. P., On some consequences and possible causes of layered mantle convection, Journ. Geoph. Res., 86, 6133-6142, 1981.

Rona P. A., and Scott S. D., A Special Issue on Sea-Floor Hydrothermal Mineralization: New Perspective. Preface, Economic Geology, (8), 1935-1975, 1993.

Sandwell D. T., and Smith W. H., Marine gravity anomaly from Geosat and ERS-1 satellite altimetry, Journ. Geoph. Res., 102, (B5), 10,039-10,054, 1997.

Schilling J. G., Iceland mantle plume, Nature, 246, 141-143, 1973.

Schiliing J. G., Azores mantle blob: rare earth evidence, Earth Planet. Sci. Lett., 25, 103-115, 1975.

Schilling J. G., Geochemical and isotopic variation along the Mid-Atlantic Ridge axis from 79^o N to 0^o N, The Geology of North America. V. M, The Western North Atlantic Region, Vogt P. R. and Tucholke B. E. (eds.), The Geological Society of America, pp. 137-157, 1986.

Silantyev S. A., Dmitriev L. V., Bazylev B. A., Casey  J. F., Bougault  H., Levsky L. K., Belyatsky B. V., and Ovchinikova G. V., An Examination of Genetic Conformity Between Co-existing Basalt, Gabbro and Residual Peridotites from 15^o20' N Fracture Zone, Centrasl Atlantic: Evidence from Isotope Composition of Sr, Nd and Pb, InterRidge News, 4, (2), 18-21, 1995.

Sobolev A. V., and Dmitriev L. V., Primary melts of tholeiites of oceanic rifts (TOR): evidence from studies of primitive glasses and melt inclusions in minerals, 28th IGC, pp. 147-148, Washington D.C. USA, Abstracts, 1989.

Sparks D. W., and Parmentier E. M., The structure of three-dimentional convection beneath oceanic spreading centres, Geophys. J. Iny., 112, 81-91, 1993.

Wilson M., Igneous Petrogenesis, London Unwin Hyman Dostal, 463 p., Boston, Sidney, Wellington, 1989.

Zhang Y. S., and Tanimoto T., Ridges, hotspots and their interaction, as observed in seismic velosity maps, Nature, 355, 45-49, 1992.