Природа геохимической зональности вкрест простирания Камчатской островной дуги

Чурикова Т., Дорендорф Ф., Вёрнер Г.,  Институт вулканической геологии и геохимии ДВО РАН, Россия, Геохимический институт Гёттингена, Германия

Макро- и микроэлементы, а так же изотопы Sr, Nd, Pb, U, Th и O в породах и S, Cl и F в расплавных включениях были проанализированы в лавах северного пересечения Камчатки, которое простирается на 220 км вкрест дуги от во фронтальной зоны (ВВФ) через Центральную Камчатскую депрессию (ЦКД) к Срединному хребту (СХ). Было опробовано 9 верхнеплейстоценовых и голоценовых стратовулканов и 2 больших лавовых поля расположенных в 110 до 400 км над поверхностью субдуцируемой плиты, что позволило охарактеризовать пространственные вариации пород, а так же относительное количество и состав субдукционного флюида, вовлеченного в магмогенезис.

Типичные Камчатские островодужные базальты, нормализованные к 6% MgO обнаруживают обогащение щелочами, LILE, LREE и HFSE от фронта к тылу дуги. При этом Ba/Zr- и Ce/Pb-отношения примерно постоянны вкрест дуги, предполагая близкую флюидную добавку в мантийные источники пород. La/Yb и Nb/Zr возрастают от ВВФ к ЦКД и далее к СХ. Породы ЦКД обнаруживают широкий спектр в 87Sr/86Sr отношении, которое варьирует от 0,70334 to 0,70366, а также наивысшие, но постоянные изотопные отношения Nd, что коррелирует с наивысшими по пересечению U/Th и 18О значениями. Изотопы Pb близки к источнику MORB, чуть уменьшаясь от ВВФ к СХ. Все это предполагает, что мантийный источник под ВВФ и ЦКД близок к N-MORB, однако обогащен (добавкой источника типа OIB) под СХ. Степень плавления, оцененная по отношению CaO6.0 и Na2O6.0, уменьшается от ВВФ к ЦКД и сохраняется постоянной в ЦКД и СХ.

S и Cl в Камчатских лавах в основном контролируются процессами дегазации. Обогащение ВВФ и ЦКД по Cl и S указывает на большую добавку флюида в этих зонах. Причиной высоких F и F/Cl в расплавных включениях СХ может быть мантийное плавление обогащенных фтором фаз или добавка глубинных флюидов.

Введение

Плавление мантийного вещества под островными дугами происходит при взаимодействии мантийного клина с водонасыщенными флюидами, отделяющимися от субдуцируемой океанической плиты [6,7 и др.]. Такие флюиды обогащены крупными литофильными элементами (LILE, в т.ч. Cs, Rb, K, Ba, Pb) и легкими редкоземельными элементами (LREE), но обеднены высокозарядными (HFSE, в т.ч. Nb, Ta, Zr, Hf) и тяжелыми редкоземельными элементами (HREE). Наиболее спорными вопросами сегодня являются состав мантийного клина, состав и количество субдукционного флюида, а так же природа его взаимодействия с мантией.

В ряде работ 19,22] c целью учесть вклад флюидной компоненты в мантийный источник, использовались содержания несовместимых редких элементов в примитивных островодужных вулканических породах. Однако, проблемы в этих расчетах связаны с наличием множества дополнительных факторов, влияющих на геохимический состав породы. Важнейшими из них являются: состав мантийного источника, различная степень его плавления, вклад осадочного материала субдуцируемой плиты, субконтинентальной литосферы или коровых пород. Изучение геохимических вариаций вкрест дуги дает возможность определить некоторые из этих факторов. Такие работы были осуществлены для Японской [26] и Курильской [1] островных дуг. Предварительные попытки изучения Камчатской дуги [3,15,29] принесли довольно противоречивые результаты, обусловленные в основном ограниченным количеством данных.

Камчатский полуостров, формирующий северную часть Курило-Камчатской дуги, является одним из наиболее вулканически активных регионов земли и включает более 200 четвертичных вулканов, 29 из которых активны по сей день. Исключительно высокая магмопродуктивность Камчатской дуги, наличие множества молодых и исторических извержений, а так же необычно высокий процент магм основного состава, практически не загрязненных осадочным материалом [5,17], позволяет исследовать относительно простые системы. Наиболее интенсивно голоценовая вулканическая активность Камчатки проявлена на Ключевской группе вулканов, расположенной в северной части Цетрально-Камчатской депрессии.

С целью изучения изменения геохимических характеристик пород вкрест простирания дуги и установления возможных причин этого разнообразия, нами была отобрана представительная коллекция образцов вулканических пород вдоль восточно-западного пересечения полуострова на широте Ключевской группы (рис.1), которая анализировалась на макро- и микро элементы, изотопы Sr, Nd, Pb, Th, U и O, а также содержания летучих в их расплавных включениях.

Геологическое положение изучаемых объектов и отбор образцов.

ris01.jpg (14164 bytes)

Рис. 1

 Камчатская дуга находится в северо-западной части конвергентного сочленения Евроазиатской и Тихоокеанской плит, последняя из которых в настоящее время субдуцирует со скоростью 9 см в год, увлекая за собой под Камчатку Императорский подводных хребет. Современная конфигурация плит под Камчаткой была сформирована только в позднем миоцене - раннем плиоцене. Наиболее активная фаза вулканизма отмечается в период с верхнего плейстоцена по голоцен. Четвертичный вулканизм на Камчатке (рис.1) проявлен в трех зонах, параллельных желобу дуги, а так же основному её простиранию: (1) Восточный Вулканический Фронт (ВВФ), (2) грабенообразная Центральная Камчатская Депрессия (ЦКД), включающая Ключевскую группу вулканов и (3) западная вулканическая зона Срединного хребта (СХ). В настоящее время северное окончание вулканической активности на Камчатке фиксируется на вулкане Шивелуч, что, вероятно, связано с изменением геометрии границы плит с ЮЗ-СВ конвергентного сочленения на трансформный разлом СЗ-ЮВ простирания [4,33]. Глубина сейсмофокальной зоны субдуцируемой плиты возрастает вкрест дуги от 100-140 км под ВВФ до 400 км под Ичинским вулканом в СХ [14].

Глубинное сейсмическое зондирование [2] показало, что мощность земной коры на Камчатке изменяется от 20 км до 42 км, возрастая с юга на север. Вкрест простирания дуги, на широте Ключевской группы вулканов, ее мощность меняется с запада на восток от 30 км под Срединным хребтом, увеличиваясь под ЦКД до 40-42 км.

9 верхнеплейстоценовых и голоценовых стратовулканов и 2 лавовых поля конусов опробовались вдоль 200-километрового траверса на севере Камчатки от ВВФ (Комарова, Гамчен, Шмидта, Кизимен) через вулканы Ключевской группы в ЦКД (Ключевской, Толбачик, Плоские Сопки, Камень) к СХ с моногенными центрами в районах Эссо и Ахтанг и изолированным Ичинским вулканом. Образцы были отобраны с особой тщательностью, чтобы избежать любых следов вторичных изменений. Глубина сейсмофокальной зоны и положение вулканов, на которых проводилось опробование, показаны на рисунке 1.

Аналитические методы

Все аналитические работы (кроме 18О) производились в геохимическом институте университета Гtттинген. Содержания в породах макро - элементов и некоторых микро элементов (Sc, V, Cr, Co, Ni, Zn, Ga, Sr, Zr, Ba) определялись рентгено-флюоресцентным анализом (РФА). Аналитические ошибки для макро элементов составили около 1% (за исключением Fe, Na: 2% и ППП: ~10%) и для микро элементов около 5%. Все остальные малые элементы определялись методом ICPMS. Ошибки, оцененные по стандартам JB3 и JA2, составляют для Nb и Ta около 15-20%, для других редких элементов менее 10%.

Изотопы Sr, Nd и Pb мерились на масс-спектрометре Finnigan MAT 262 RPQ II+ с использованием стандартов NBS987 (0.710245) для Sr, LaJolla (0.511847) для Nd и NBS981 (рекомендованные значения по [31]) для Pb. Общие ошибки (2) составили менее 0,004% для Sr и Nd и менее 0,1% для Pb. Изотопы U и Th измерялись на масс-спектрометре Finnigan MAT 262 с приставкой RPQ 2+. Измеренные изотопные отношения U и Th корректировались на фракционирование относительно U-стандарта "U-112" и Th-стандарта "Santa Cruz". Несмотря на то, что бланковые анализы были ниже, чем 0,3 ppb для U и 0,08 - 0,31 ppb для Th, общая ошибка несколько завышена - около 5%. Изотопные отношения кислорода в оливинах были измерены в университете Карнеги, с использованием серии Synrad 48 CO2 лазера. Изотопный состав кислорода определялся на массах 32, 33 и 34 на масс-спектрометре Finnigan МАТ-252. Стандартный газ был откалиброван по шкале SMOW, используя NBS-28 (deltasm.gif (63 bytes)18O=9,60 ). Внешняя погрешность метода менее чем +0,2 .

Макроэлементы, сера, хлор и фтор, а также валентность серы в расплавных негомогенизированных включениях в оливинах и клинопироксенах измерялись на микрозонде JEOL8900 WDS со стандартным набором синтетических и природных стандартов.

Более детальное описание всех методик можно найти в [9,10,11,12].

Результаты и обсуждения

Макро- и микроэлементы

ris02.jpg (10052 bytes)

Рис. 2

 Макро - и ряд микроэлементов определены в 152 представительных образцах пересечения, в 72 из них был измерен более широкий спектр редких элементов.

Породы ВВФ, включая вулкан Кизимен, относятся к средне-калиевым сериям (рис.2). Некоторые редкие низко-калиевые толеитовые породы встречаются на вулканах Гамчен и Шмидта. Наиболее высокие щелочи наблюдаются в породах СХ, лавы которого представлены средне-высоко-калиевыми известково-щелочными сериями. Вблизи основания стратовулкана Ичинский были опробованы шлако-лавовые базальтовые конуса, обогащенные HFSE элементами с внутриплитными геохимическими признаками (т.н. базальты внутриплитного типа - ВПТ). При этом сам стратовулкан сложен типично островодужной андезит - дацит - риолитовой серией. Большинство пород ЦКД среднекалиевые известково-щелочные. Некоторые лавы вулканов Плоские Сопки и Толбачик (в т.ч. Южный прорыв 1975-76 гг.) относятся к высоко-калиевым сериям, природа которых требует специального рассмотрения, и в предлагаемой статье излагаться не будет.

ris03sm.jpg (6898 bytes)

Рис. 3

 Распределения редких элементов на спайдерограммах для ВВФ, ЦКД и СХ показаны на рис.3 для пород с > 5% MgO. Все породы имеют типичные островодужные признаки с различным обогащением LILE и LREE и низкими HFSE. Исключением являются несколько моногенных конусов ВПТ пород (см. ниже). Концентрации LILE и HFSE возрастают от фронта к тылу дуги. Породы ВВФ и ЦКД обеднены Nb и Ta в сравнении с составом NMORB (north middle ocean ridge basalt). Интересной особенностью всех изученных пород Камчатки являются низкие концентрации HREE, которые значительно ниже, чем в NMORB и не меняются значимо во всех трех регионах. Базальты ВПТ Ичинского вулкана обогащены больше, чем островодужные породы стратовулкана по LILE и LREE с повышенными HFSE. Nb-Ta отрицательная аномалия в породах СХ выражена значительно слабее, чем в островодужных породах пересечения.

Ввиду ограниченности объема публикации, все первичные данные представлены на EPSL и J.Pet. Online Background Dataset.

Коррекция на фракционную кристаллизацию

ris04sm.jpg (14011 bytes)

Рис. 4

 Наиболее основные образцы в ВВФ, ЦКД и СХ содержат до 8.5%, 11.6%, и 9.2% MgO, соответственно. Однако очевидно, что большинство образцов претерпели процесс минерального фракционирования и прямое сравнение полученных концентраций микроэлементов в них невозможно. Чтобы уменьшить этот эффект, мы попытались скорректировать первичные данные к примитивным составам магм.

ris05sm.jpg (15199 bytes)

Рис. 5

 С этой целью мы использовали метод, описанный в [24]. Для каждого вулкана были построены зависимости содержания элементов от MgO и рассчитано содержание каждого элемента в точке пересечения трендов с MgO = 6%. Значение пересечения линии регрессии с 6% MgO принималось за скорректированное значение для каждого элемента.

По корреляциям макроэлементов с MgO (не показано), фракционирование плагиоклаза начинается, когда расплав обедняется примерно до 5% MgO. Поэтому в наших расчетах линий регрессии с последующей нормализацией концентраций элементов мы использовали только образцы, в которых содержание MgO превышало 5%. Нормализованные значения отмечены как 6,0.

Для островодужной серии вулкана Ичинский наши данные дополнены неопубликованными данными О.Волынца (ныне умершего; данные доступны по просьбе), так же привлечены литературные источники. Менее 5 образцов с > 5% MgO имеется для вулканов Камень (2), Гамчен (4), Ахтанг (3) и Шмидт (2). Для графиков зависимостей микроэлементов от MgO относительное отклонение от линий регрессий составило около 10-20%.

Вариации макро- и микроэлементов вкрест простирания Камчатской дуги

ris06sm.jpg (17330 bytes)

Рис. 6

 Планк и Лангмюр [24] предположили, что степень плавления мантийного клина зависит от мощности земной коры. Небольшие изменения в мощности коры под Камчаткой (30 - 40 км) не достаточны, чтобы объяснить наблюдаемую значительную разницу в распределении редких элементов в породах пересечения только этим эффектом.

Расстояние от глубоководного желоба до вулкана и глубина до поверхности субдуцируемой плиты связаны напрямую. Общеизвестное положение фронтальных вулканических зон не менее чем в 110-130 км над поверхностью субдуцируемой плиты предполагает, что плавление обусловлено флюидами, освобождающимися из этой плиты в результате реакции дегидратации при достижении определенных температур и давлений [25,29 и др.]. Поэтому в данной работе при региональном сравнении пород мы использовали глубину поверхности погружающейся плиты под вулканами [14]. На рисунках 4,5,6 показаны диаграммы зависимости некоторых макро- и микроэлементов, а так же их отношений от глубины погружающейся плиты.

ris07sm.jpg (13700 bytes)

Рис. 7

 Базальты ВПТ помечены на диаграммах специальными символами. Большинство нормализованых к MgO6.0 редких элементов, в т.ч. HFSE (Zr, Nb, Hf, Ta), LILE (Sr, Ba, Rb, Be, Pb, U, Th), LREE, некоторые макроэлементы (K, Na) и многие отношения элементов (K/Na, La/Yb, Sr/Y, Nb/Yb) обнаруживают положительную корреляцию с глубиной погружающейся плиты. Наилучшие корреляции получены для K2O, Ba, Sr и Rb, для которых значения возрастают более чем в два раза от фронта к тылу дуги. Na2O, LREE и HFSE стремительно возрастают от ВВФ к ЦКД, но заметно слабее далее к СХ. При этом Ti даже понижается в лавах СХ в сравнении с породами ЦКД. Для Y и HREE нет значимых корреляций, они остаются постоянными от фронта к тылу дуги. Эти результаты сравнимы с аналогичными исследованиями по Курильской [1] и Японской [26] вулканическим дугам.

Базальты ВПТ Срединного хребта отличаются по ряду элементов от островодужных трендов, выделенных затененным полем на рисунках 4,5,6. В отличие от Ичинского стратовулкана, они имеют высокие концентрации Na2O, TiO2, P2O5 и всех HFSE и REE, и обеднены по SiO2 и Pb. HFSE и LREE в этих породах значительно выше, чем во всех образцах пересечения. Для этих породы так же типичны высокие Ce/Pb, La/Yb и низкие U/Th и Ba/Nb отношения.

ris08sm.jpg (11325 bytes)

Рис. 8

 Нормированные к 6% MgO отношения U/Th во всех породах Камчатки выше, чем в NMORB и варьируют от 0,41 до 0,58 в ВВФ, от 0,57 до 0,71 в ЦКД и от 0,38 до 0,64 в породах СХ. Исключением среди пород ЦКД является образец 2310 с вулкана Камень, в котором первичное U/Th отношение достигает 0,79. Это наивысшее известное U/Th отношение в примитивных базальтах Камчатки. Наряду с этим, породы вулкана Камень характеризуется самыми низкими LILE и LREE и самыми высокими HREE в сравнении с другими вулканическими сериями ЦКД (рис.6Д, 9).

Sr-, Nd-, и Pb-изотопы

Sr-, Nd- и Pb-изотопные данные для пород пересечения приведены на рис.7 и 8. Фигуративные точки ложатся довольно близко к области MORB. Имея представительную коллекцию образцов, мы можем идентифицировать внутри Камчатского поля более мелкие структуры, характерные для каждого региона. В целом наблюдается возрастание 87Sr/86Sr и 143Nd/144Nd отношений от ВВФ к ЦКД и дальнейшее их убывание от ЦКД к СХ (рис.7 и 8). Поля точек ВВФ и СХ очень близки на Sr-Nd диаграмме, за исключением двух образцов с вулкана Комарова, имеющих повышенные значения 87Sr/86Sr. Наблюдается широкий диапазон значений по Nd изотопам для ВВФ и СХ, в то время, как Sr изотопные отношения близки. Наивысшее обогащение по 87Sr в пределах ЦКД найдено для лав Ключевского вулкана, где отношение 87Sr/86Sr достигает 0.70366.

По Pb изотопной систематике лавы ЦКД менее радиогенные, чем породы ВВФ, но близки к полю СХ (рис.7 В). Породы ВПТ идентичны с другими лавами СХ. Для сравнения, на рис.8 показаны данные для вулкана Бакенинг (ВВФ, 200 км к югу от пересечения, [12]). Эти породы еще менее радиогенны по Sr при сравнимых значениях Nd и Pb изотопов. Исследования клинопироксенов (предварительно очищенных) из мантийных ксенолитов Камчатки (неопубликованные данные) показало, что по Nd изотопам мантийные ксенолиты близки вулканическим породам, но по Sr отношениям разброс точек ещё более широкий.

Изотопы кислорода

Подробно изотопная систематика кислорода в породах Ключевского вулкана изложена в [11]. Поэтому здесь мы только подчеркнем наиболее важные результаты.

Изотопные отношения кислорода в одиночных зернах оливинов из различных пород Камчатского пересечения изменяются от 5,6 до 7,4 . Максимальные значения отмечены для вулканов Ключевской группы. Значения изотопов кислорода в расплавах Ключевского вулкана рассчитывались как среднее deltasm.gif (63 bytes)18O в оливине с учетом эффекта фракционирования оливин-расплав при 1100-1200o. Диапазон значений (6,2-7,5 ) оказался явно шире, а значения выше, чем для типичных мантийных расплавов [11].

ris09sm.jpg (10613 bytes)

Рис. 9

 Установлены положительные корреляции deltasm.gif (63 bytes)18O с отношением 87Sr/86Sr (рис.8), причем лавы исторических извержений имеют тенденцию к повышенным значениям обоих отношений в сравнении с более ранними голоценовыми извержениями. Установлены положительные корреляции между deltasm.gif (63 bytes)18O и Cs, Li, Sr, Ba, Rb, Pb, Th, U, LREE и K, то есть с подвижными во флюиде элементами, а так же с отношениями K2O/Na2O, Ba/Zr и La/Yb. Отрицательная корреляция deltasm.gif (63 bytes)18O была найдена с U/Th отношением, однако корреляции с подвижными в расплаве Al и HFSE отсутствуют.

U-Th систематика

Самые низкие значения U/Th элементных отношений (0,35 - 0,6) наблюдаются в породах СХ за исключением нескольких точек с более высокими отметками, а наивысшие - в породах ЦКД (0,5 - 0,9). В породах ВВФ это отношение меняется от 0,45 до 0,75. U-Th изотопная систематика представлена на диаграмме (230Th/232Th)-(238U/232Th) (рис.9). На графике очевидны те же закономерности: наименьшие значения типичны для пород СХ, породы ВВФ имеют средние значения и лавы ЦКД характеризуются наивысшими изотопными отношениями, достигая значений 2,29 для (238U/232Th) и 2,15 для (230Th/232Th) в породах вулкана Камень.

Значительное 238U-230Th изотопное неравновесие

Страницы: 1 2 3