Бесплатный автореферат и диссертация по наукам о земле на тему
Поведение изотопно-кислородной системы в процессах контаминации МАГМ и взаимодействия флюид-порода
ВАК РФ 25.00.09, Геохимия, геохимические методы поисков полезных ископаемых

Автореферат диссертации по теме "Поведение изотопно-кислородной системы в процессах контаминации МАГМ и взаимодействия флюид-порода"

[ I Г)

На правах рукописи

АВДЕЕНКО АННА СЕРГЕЕВНА

ПОВЕДЕНИЕ ИЗОТОПНО-КИСЛОРОДНОЙ СИСТЕМЫ В ПРОЦЕССАХ КОНТАМИНАЦИИ МАГМ И ВЗАИМОДЕЙСТВИЯ ФЛЮИД-ПОРОДА (НА ПРИМЕРЕ НЕОГЕНОВЫХ ГРАНИТОИДОВ И МИНЕРАЛИЗОВАННЫХ ВОД РАЙОНА КАВКАЗСКИХ МИНЕРАЛЬНЫХ ВОД)

Специальность 25.00.09 - геохимия, геохимические методы поисков полезных ископаемых

АВТОРЕФЕРАТ диссертации на соискание ученой сте кандидата геолого-минералогических

Москва-2009

003481985

Работа выполнена в Учреждение Российской Академии Наук Институте геологии рудных месторождений, петрографии, минералогии и геохимии (ИГЕМ) РАН.

Научные руководители: кандидат геолого-минералогических наук

Дубинина Елена Олеговна доктор геолого-минералогических наук Носова Анна Андреевна

Официальные оппоненты: доктор геолого-минералогических наук

старший научный сотрудник Покровский Борис Глебович доктор геолого-минералогических наук старший научный сотрудник Курчавов Анатолий Михайлович

Ведущая организация: Институт геохимии и аналитической химии

им. В.И. Вернадского РАН

Защита состоится 26 ноября 2009 г. в 1400 на заседании диссертационного совета Д 002.122.01 ИГЕМ РАН по адресу 119017, Москва, Старомонетный пер., д.35

С диссертацией можно ознакомиться в библиотеке ИГЕМ РАН Автореферат разослан 21 октября

Ученый секретарь

диссертационного совета

кандидат геолого-минералогических наук

Юдовская М. А.

Актуальность темы. Гранитоидный магматизм проявляется на всех стадиях формирования континентальной коры, но наибольший масштаб гранитообразования характерен для пред-, син- и постколлизионных обстановок. С этим типом магматизма генетически связаны разнообразные (Cu-порфировые, Sn, W-Mo, редкометальные и др.) месторождения. В последние годы выявлена большая металлогеническая роль гранитов высоко-Ba-Sr (латитового) типа, высокодифференцированные разновидности которых (Li-F гранитоиды) обладают высокой рудоносностью с широким спектром месторождений - молибден-меднопорфировых, редкометальных, урановых.

Образование гранитов высоко-Ba-Sr типа на сегодняшний день связывают с плюмовым магматизмом [Тагпеу, 1994]. Для этих гранитов характерны высокие концентрации Ва, Sr и Р, низкие концентрации Rb, Th, U, Y, Nb и тяжелых редкоземельных элементов (HREE), а так же высокое отношение K/Rb. Источником таких пород могут быть «обогащенные» мантийные расплавы. Однако одной из проблем, возникающей при исследовании гранитов высоко-Ba-Sr типа, являются довольно радиогенные изотопные отношения Sr, умеренно радиогенные Nd и высокие значения величины 5|80, которые не характерны для магматических пород, в образовании которых существенная роль отводится мантийным расплавам. Формирование изотопных параметров этого типа гранитов в основном объясняют процессами фракционной кристаллизации и ассимиляции осадочных пород «обогащенными» расплавами [Fowler, Henney 2001; Qian, 2003; Reiners, Nelson 1996]. Однако, этими процессами не всегда удается полностью описать поведение изотопных систем в высоко-Ba-Sr гранитах.

Данная работа посвящена исследованию проблемы формирования изотопных (Sr, Nd, О) параметров высоко-Ba-Sr гранитоидов. Для этого в качестве объекта исследования были выбраны высоко-Ba-Sr гранитоиды района Кавказских Минеральных Вод (КМВ), с которыми связано редкометальное и урановое оруденение. Изучение гранитоидов КМВ проводилось на основе исследования изотопной геохимии Sr, Nd и О с привлечением необходимых петрологических, геохимических и геологических данных. Комплексное изучение поведения изотопных систем (О, Sr, Nd,) позволяет решать широкий круг задач: исследование механизмов и установление степени контаминации магм, определение роли флюидной фазы в постмагматических процессах, решать ряд других вопросов.

Район Кавказских Минеральных Вод, в котором расположены высоко-Ba-Sr гранитоиды, отличается современной термальной активностью, выражающейся в многочисленных проявлениях деятельности минерализованных водных

источников. Параллельно с изучением самих гранитоидов КМВ было проведено изотопное (О, Н) изучение минеральных вод данного района с целью установления их типа питания и механизма формирования изотопных параметров. Полученные данные представляют интерес как для прогноза дальнейшего геохимического и гидрогеологического состояния уникального природного объекта - минеральных вод района КМВ, так и для решения общих проблем генезиса минерализованных вод.

Настоящее исследование имело следующие цели: (1) Определить роль процесса контаминации расплавов при формировании изотопно-кислородного состава высоко-Ba-Sr гранитоидных пород на примере гранитоидов КМВ; (2) Выявить роль взаимодействия водный флюид-порода при становлении гранитных интрузивных тел по данным изучения изотопно-кислородной системы гранитоидов КМВ; (3) Установить степень взаимодействия вода-порода при низких температурах на примере изучения изотопного состава водной фазы минеральных источников района КМВ.

Исходя из намеченных целей, были поставлены следующие задачи исследования: (1) Усовершенствование методики фторирования для повышения точности измерения изотопного состава кислорода в валовых пробах и мономинеральных фракциях силикатных пород, модернизация метода изотопного уравновешивания для определения б180 в природных водах; (2) Анализ изотопного состава кислорода в породах и минералах-вкрапленниках гранитоидных тел КМВ; (3) Анализ изотопного состава кислорода и водорода минерализованных и пресных вод района КМВ; (4) Интерпретация данных по изотопному составу кислорода и радиогенным изотопам (Sr, Nd) с использованием моделей, описывающих их поведение в процессах взаимодействия расплав-ассимилянт и вода-порода.

Фактический материал. В основу диссертации легли результаты исследований, проведенных автором в лаборатории изотопной геохимии и геохронологии ИГЕМ РАН. Образцы пород для изотопных исследований были отобраны при участии автора в течение двух полевых сезонов (2003-2004 гг.) в районе Кавказских Минеральных вод (КМВ), часть образцов любезно предоставлены А.Я. Докучаевым и A.A. Носовой (ИГЕМ РАН). Изотопные исследования сопровождались петрографическим и геохимическим исследованием пород.

Главные методы исследования. При исследовании гранитных и вмещающих пород района КМВ было проанализировано 23 образца пород Rb-Sr методом, 18 образцов пород Sm-Nd методом, в 39 образцах пород определен изотопный состав кислорода и в 14 образцах пород определен изотопный состав

углерода. Для 18 минеральных фракций из гранитных пород изучался изотопный состав кислорода. Было выполнено 20 парных определений изотопного состава водорода и кислорода в образцах минерализованных вод района КМВ. Изотопные исследования пород сопровождались микроскопическим изучением образцов гранитных пород и микрозондовым исследованием составов минералов (геологический факультет МГУ, H.H. Коротаева), породообразующие оксиды в породах определялись методами химического анализа в ИГЕМ РАН (О.Г. Унанова). Редкие элементы в породах определялись методом ICP-MS в ИМГРЭ (Д.З. Журавлев) и методом инструментального нейтронно-активационного анализа (INAA) в ИГЕМ РАН (А.Л. Керзин).

Научная новизна работы. Впервые для гранитоидов КМВ проведено комплексное изотопно-геохимическое (О, Sr и Nd) и геохимическое изучение. Изотопный состав кислорода в породах КМВ изучен впервые. Показана принципиальная возможность преимущественно корового происхождения высоко-Ba-Sr гранитоидов с формированием их геохимической специфики в результате химического взаимодействия кислого расплава с высоко-Mg карбонатной породой. Проведены расчеты формирования изотопных параметров гранитоидов КМВ в рамках одностадийного процесса AFC(+X) без привлечения дополнительного источника мантийных расплавов. Предложен механизм флюидно-магматического взаимодействия, приводящего к формированию высоких величин 5|80 в лейкогранитах. Подтверждено инфильтрационное происхождение минерализованных углекислых вод КМВ, впервые получена количественная оценка гипсометрической высоты области их питания.

Практическая значимость работы. С гранитами высоко-Ba-Sr типа часто связаны рудные месторождения: молибден-порфировые (Клаймакс), полиметаллические и золоторудные (Забайкалье (Приаргунье)), урановые (район КМВ), месторождения бора (Дальний Восток) и ряд других. Проведенные исследования условий формирования гранитов высоко-Ba-Sr типа дают информацию, которая может быть использована при построении генетических моделей рудообразования, связанного с позднеколлизионным магматизмом. Выявление роли взаимодействия флюид-порода и установление механизма контаминации магм, которым посвящена значительная часть работы, имеют большое значение для расшифровки источника и процессов транспорта рудного вещества.

Для рационального использования и обоснования мероприятий по защите основного природного ресурса района - месторождений минерализованных вод, результаты мониторинга их изотопных параметров во времени могут дать

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

Основные защищаемые положения.

1. Впервые определен изотопный состав кислорода в валовых пробах и сквозных минералах-вкрапленниках гранитоидов КМВ. Установлены закономерности поведения изотопно-кислородной системы пород, коррелирующиеся с их минералого-петрологическим типом: для группы граносиенитов и сиенитов характерен узкий диапазон вариаций величины 5180 (от 9.1 до 10.4 %о), для группы амфиболовых гранитов диапазон вариаций величины 5180 шире и характерны в целом ее более низкие значения (от 7 до 9.1 %о), наибольший размах вариаций и наиболее высокие значения величины 5180 характерны для группы лейкогранитов (от 8.9 до 13.8 %о). Изотопно-кислородное равновесие в сквозных минералах-вкрапленниках всех исследованных типах гранитоидов отсутствует.

2. На основе комплексного изотопно-геохимического (О, N(1, Бг) и геохимического изучения гранитоидных пород структуры КМВ установлено, что в истории их формирования имел место процесс контаминации исходных расплавов осадочными породами на глубине. Поведение изученных изотопных систем может быть описано в рамках реакции кислого силикатного расплава с породами карбонатного состава.

3. На основе расчета по линейной модели изотопного обмена в процессе кристаллизации расплава показано, что высокие значения и значительный интервал вариаций величин б|80 в лейкогранитах КМВ могут быть следствием взаимодействия флюид-порода на стадии становления интрузивных тел в близповерхностных условиях.

4. Усовершенствованы методики изотопного анализа кислорода и водорода в образцах природных вод. С применением этих методик изучен изотопный состав минерализованных и пресных вод района КМВ. Показано, что основным компонентом минеральных источников КМВ являются инфильтрационные метеорные воды, испытавшие процесс изотопно-кислородного обмена с силикатными породами при низких температурах.

Структура и объем работы. Диссертация состоит из 5-ти глав, введения, заключения и дополнительных материалов, включающих результаты химических анализов пород (Приложение 1). Основной материал работы изложен на 133 страницах, которые содержат 15 таблиц и 50 рисунков. Список цитируемой литературы включает 125 наименований.

Апробация работы. По теме диссертации в рецензированных журналах опубликовано 3 статьи и одна сдана в печать. Результаты исследований докладывались на 3-х российских и одной международной конференции: XVI Симпозиуме по геохимии изотопов им. акад. А.П. Виноградова, (Москва, 2001); 14th V.M. Goldschmidt conference (Copenhagen, 2004); XVIII Симпозиуме по геохимии изотопов им. акад. А.П. Виноградова (Москва, 2007) и IV Российской конференции по изотопной геохронологии (Санкт - Петербург, 2009).

Благодарности. Автор выражает искреннюю благодарность своим научным руководителям Е.О. Дубининой и А.А. Носовой за советы и консультации в процессе работы над диссертацией. Автор благодарен заведующему лабораторией изотопной геохимии и геохронологии ИГЕМ член-корреспонденту РАН И.В. Чернышеву за постоянную помощь и ценные советы при написании работы. Автор благодарен Э.Д. Баировой, Т.Н. Олейниковой и Ю.В. Гольцману за помощь в проведении изотопного анализа стронция и неодима. К.Н. Шатагину, А.В. Самсонову, Ю.О. Ларионовой и М.М. Аракелянц за ценные комментарии и полезные замечания, полученные при обсуждении результатов исследования. Особая благодарность А.В. Перкову, Н.И. Сердюку, Г.И. Сердюку и А.В. Леляеву, за помощь в усовершенствовании и поддержании работоспособности масс-спектрометрического оборудования лаборатории изотопной геохимии и геохронологии ИГЕМ РАН.

Автор признателен А.Я. Докучаеву за помощь в проведении и организации полевых работ, Л.В. Сазоновой за помощь в минеральных и петрологических исследованиях каменного материала, Ю.А. Костицыну (ГЕОХИ РАН) за ценные замечания при обсуждении отдельных положений диссертации.

Глава 1. Происхождение и изотопно-геохимические характеристики высоко-

Ba-Sr гранитоидов

Данная глава представляет собой литературный обзор, отражающий современные представления о происхождении высоко-Ba-Sr типа гранитов и гипотезы образования гранитоидов района КМВ, как типичных представителей этого типа гранитов.

Геохимическими критериями выделения высоко-Ba-Sr гранитов в отдельный тип являются высокие содержания Ва и Sr в этих породах по сравнению с традиционными S-, I- и А- типом гранитов. Происхождение этого типа гранитов чаще всего связывают с плюмовым магматизмом [Ташеу, 1994], считая, что источником высоко-Ba-Sr гранитов могут быть «обогащенные» мантийные расплавы, которые образовались в результате: 1) плавления вещества

субдуцированных океанских островов или плато, обогащенных Ва и Sr (и Р);

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

3) плавления вещества низов литосферы при внедрении карбонатных расплавов из астеносферы. Для указанного типа гранитов характерны высокие значения

87 86 18 16

изотопных отношений Sr/ Sr и О/ О, которые многие исследователи [Reiners, Nelson 1996; Qian, 2003; Fowler, Henney, 1996 и др.] считают результатом процессов фракционной кристаллизации и ассимиляции мантийных расплавов. Однако данными процессами не всегда удается объяснить наблюдаемые изотопные параметры высоко-Ba-Sr гранитов.

Исследование высоко-Ba-Sr гранитов района КМВ имеет богатую историю, однако почти все предшествующие исследования, достаточно разносторонние, обошли вниманием методы изотопной геохимии, весьма информативные в вопросах происхождения гранитоидов. Исключением является работа [Поль и др., 1993], в которой приведены результаты изучения Sr и РЬ изотопных систем для отдельных интрузивов КМВ.

Анализ литературных данных показывает, что на сегодняшний день существуют три гипотезы образования гранитоидов КМВ: 1) ассимиляция гранитной магмой вмещающих карбонатных пород в зоне становления интрузивных тел [Соболев и др., 1959; Короновский и др., 1986]; 2) смешение кислой трахириолитовой магмы с мантийными щелочными базальтами и их дериватами [Поль и др., 1993]; 3) магматическая дифференциация расплавов, которая сопровождалась незначительной контаминацией вмещающими породами и участием флюидной фазы на поздних стадиях дифференциации [Носова и др., 2005]. Тем не менее, выдвинутые гипотезы не дают исчерпывающих объяснений полученным в настоящей работе новым данным, в частности, по поведению изотопно-кислородной системы в гранитоидах КМВ.

Глава 2. История изучения, геологическая, геохимическая и минералогическая характеристика гранитоидов района Кавказских Минеральных вод

Район КМВ является составной частью Кавказского сегмента Альпийского складчатого пояса и располагается на северном фланге Эльбрусской вулканической области. Гранитоиды КМВ являются одними из самых ранних проявлений постколлизионного магматизма в регионе: их внедрение произошло в позднем миоцене (8.3 млн. лет. назад [Лебедев и др., 2005]). Все гранитоиды района КМВ приурочены к крупному Транскавказскому субмеридиональному поперечному

Л(мцюры

/ЗмёйкЗ.

Грахитсвые Холмы

т. Лысзя

¡Ыти.-щкь X-

'ссенмук/

г. Золотой^Ку^н

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

Выходы интрузивов района КМВ образуют структуру, близкую к кольцевой (рис.1). Эта структура хорошо выражена в гравитационном и магнитном полях [Смирнова, 1999]. Структура имеет зональное строение: во внешней части распространены интрузивы сиенитового и граносиенитового состава (гг. Верблюд, Лысая, Золотой Курган, Змейка, Развалка, Железная, Юца, Острая, Тупая, Медовая, Кинжал). Ядро структуры слагают интрузивы гранитного и лейкогранитного составов (гг. Бештау, Козьи Скалы, Бык, Шелудивая).

1 ; I 1 1 Жемезно/о<)ск ./ч?,

V ''' V ГЛ^Йв^ЛнЙЛмяй

' * 7 ч '-л'Л-,

V г. Острая +,

фл г ¡Бештау

Рг/с. 7. Геологическая схема района Кавказских Минеральных Вод

Условные обозначения: 1-4 - выходящие на поверхность сиениты (1), субщелочные граниты (2), граносиениты (3) и лейкограниты (4); 5-8 - невскрытые эрозией (прерывистый контур - по геоморфологическим данным, состав - предполагаемый по аналогии со вскрытой частью) сиениты (5), субщелочные граниты (6), граносиениты (7), лейкограниты (8); 9 - установленные в пределах интрузивов разломы; 10-13 - осадочные образования: четвертичные (10), неогеновые глины (11), палеогеновые глины и мергели (12), меловые песчаники, аргиллиты, алевролиты, мергели и известняки (13); периферическая и центральная часть кольцевой структуры КМВ (14). Схема составлена с использованием [Носова и др., 2005]

Во многих интрузивах района встречается большое количество ксенолитов (иногда до 20% объема породы), имеющих осадочное, метаморфическое и магматическое происхождение. Среди ксенолитов кристаллических пород доминируют две группы: габброиды и включения апатит-карбонат-клинопироксен-флогопит-плагиоклазовой (Ap-Carb-CPx-Phlog-Pl) породы.

Все гранитоиды КМВ имеют порфировидную структуру и содержат от 10-15 до 20-30% минералов-вкрапленников, среди которых наиболее распространены плагиоклаз, K-Na полевой шпат, кварц, флогопит, амфибол и клинопироксен. По минеральному составу гранитоидные породы КМВ можно разделить на три группы. Первая группа -амфиболовые граниты, в которых из вышеперечисленных минералов-вкрапленников практически отсутствуют флогопит и клинопироксен. Вторая группа - сиениты и граносиениты, в которых нет кварца и роговой обманки. Третья группа - лейкограниты, для которых характерно отсутствие такого минерала-вкрапленника как клинопироксен. Массивы, сложенные первой и второй группами пород, отличаются друг от друга положением в структуре и составом ксенолитов. Если группа амфиболовых гранитов располагается в центре структуры КМВ и в ней нет ксеногенных включений Ap-Carb-CPx-Phlog-Pl породы, то группа сиенитов и граносиенитов располагается на периферии структуры КМВ и характеризуется большим количеством ксеногенных включений Ap-Carb-CPx-Phlog-Pl состава.

По геохимическим особенностям гранитоиды КМВ являются типичными представителями гранитов высоко-Ba-Sr типа. Для гранитоидных пород КМВ характерны высокие концентрации Ва (до 2360 мкг/г), Sr (до 1380 мкг/г), РЬ (до 400 мкг/г) и умеренные - Rb (до 290 мкг/г). Спектры РЗЭ гранитоидов КМВ имеют сильно фракционированный характер ((La/Yb)N > 20).

Глава 3. Методы изотопных исследований

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

Основное внимание в этой главе уделено методам определения изотопного состава кислорода и водорода в изучаемых образцах. Все измерения изотопного состава водорода и кислорода выполнялись на масс-спектрометре DELTAplus фирмы ThermoQuest (Finigan). Подготовка образцов для определения величин 5180 в породах и минералах проводилась с использованием метода фторирования [Vennemann and Smith., 1990], каждый образец подвергался разложению в присутствии BrF5 и анализировался не менее 4-6 раз. Правильность полученных значений 5|80 контролировалась регулярными измерениями международных

стандартов NBS-28 (кварц) и NBS-30 (биотит). Погрешность полученных значений величин 5180 находилась на уровне ± 0.3 %о для мономинеральных фракций и не хуже ± 0.3 - 0.4 %о для образцов валовых проб. Определение изотопного состава кислорода в водных образцах проводилось методом изотопного уравновешивания кислорода воды с С02 в одноразовых шприцах [Авдеенко, Дубинина, 2001]. Правильность определения величины 6180 в водных образцах контролировались повторными измерениями международных образцов сравнения МАГАТЭ, погрешность использованного метода составила ± 0.2 %о. Изотопный состав водорода в минеральных водах КМВ анализировался путем разложения водных проб на горячем (800 °С) хроме с помощью периферийного устройства H/Device фирмы ThermoQuest (Finigan). Точность определения величин D в водных пробах составила ± 0.3 %о.

Измерения изотопного состава неодима и стронция валовых проб пород КМВ проводились из тех же порций, что и анализ изотопного состава кислорода. Измерения изотопных отношений 87Sr/86Sr и 143Nd /144Nd образцов проводились на многоколлекторном термоионизационном масс-спектрометре Micromass Sector 54 (Великобритания). Погрешность анализа (±2а) не превышала 0.003 % для 87Sr/86Sr и 0.002 % для 143Nd/'44Nd. Правильность измерений изотопного состава стронция контролировалась по международному стандартному образцу SRM-987. При изотопном анализе Nd использовался внутрилабораторный стандартный образец с изотопным отношением l43Nd/,44Nd = 0.512406, нормализованным к международному стандартному образцу La Jolla. Среднее значение, полученное по внутрилабораторному стандартному образцу в период измерения проб гранитоидов КМВ, составило 143Nd/144Nd = 0.512401 ± 0.000003.

Глава 4. Формирование изотопных (Sr, Nd, О) параметров высоко-Ba-Sr гранитоидов позднемиоценовых интрузивов района КМВ

Для гранитоидов КМВ характерен узкий интервал начальных изотопных отношении (87Sr/86Sr)0 = 0.7083-0.7086 и небольшие отрицательные значения величины eNd(T) - от -4.2 до -2.1. Близость начальных изотопных отношений Sr в гранитоидах КМВ (в совокупности с другими геологическими и петрологическими данными) позволяет сделать вывод о том, что расплав, из которого формировались данные породы, по-видимому, был когда-то одним и тем же. В общей последовательности поведения изотопных отношений стронция и неодима эффузивных и интрузивных пород в Кавказском регионе (рис. 2) гранитоиды КМВ занимают крайнее положение - в области наиболее радиогенных изотопных

составов Sr и наименее радиогенных Nd, что указывает на минимальное участие в их формировании расплавов с мантийными изотопными характеристиками (Рис. 2).

Практически все вулканические породы Большого Кавказа характеризуются значительными концентрациями стронция (в среднем около 550 ррт). Это обогащение Sr связывают с ассимиляцией корового вещества мантийными расплавами [Бубнов, 2003, Лебедев и др., 2006], которая могла протекать как до, так и после внедрения расплавов.

12 10

8 6

Р 4

■О 7

Z 2 ы о

-2 -4 -6

-10

- DM^ \ \чр о " \ □ 10 203 □ ° • 4 + 5 IS!*, □

Oofl о (87Sr/86Sr)o

0.702 0.703 0.704 0.705 0.706 0.707 0.708 0.709 0.71

Рис. 2. Изотопный состав Ш и 5У гранитоидов КМВ и неогеновых магматических пород Кавказа

1- вулканические породы Восточной Анатолии и Большого Кавказа [Реагсе е1 а1., 1990; Бубнов, 2003; Лебедев и др., 2003 и др.];

2-вулканиты Эльбрусской области [Бубнов, 2003, Лебедев и др., 2006]; 3 - Эльджуртинские граниты (Носова, неопубликованные данные]; 4, 5- гранитоиды КМВ и ксеногенные включения в них [данная работа]

По сравнению с вулканитами Большого Кавказа, концентрации Бг в гранитоидах КМВ еще более высокие (до 1380 ррт), что заставляет предполагать ассимиляцию пород, богатых 8г, например, осадочных карбонатов. Как показало сопоставление изотопно-геохимических параметров гранитоидов КМВ и вмещающих интрузивы пород, гипотеза об ассимиляции боковых пород расплавами гранитоидов КМВ в зоне становления интрузивов оказалась несостоятельной. Карбонатные породы, находящиеся в непосредственном контакте с интрузивами, характеризуются отношениями 875г/865г (0.7078-0.7085) близкими и частично перекрывающимися с диапазоном начальных отношений Бг в гранитоидах КМВ, но при этом обладают низкими содержаниями стронция (в среднем около 700 ррт) по сравнению с гранитными породами (1000-1300 ррт). Следовательно, боковые породы не могут рассматриваться как источник повышенных концентраций Бг, характерных для гранитоидов КМВ.

Вариации величин 3|80 в гранитоидах КМВ являются существенными и зависят от минералого-петрологического типа исследованных пород (рис. 3). Для амфиболовых гранитов установлен интервал от 7.0 до 9.1 %о; для сиенитов и граносиенитов - от 9.1 до 10.4 %о; для лейкогранитов - от 8.9 до 13.8 %о. Наличие существенных вариаций и закономерное поведение величин 5180 в пределах структуры КМВ (рис. 3) может также указывать на ассимиляцию осадочных пород

расплавами гранитоидов. Изотопный состав кислорода минералов-вкрапленников (кварца, полевого шпата, флогопита, амфибола и диопсида) амфиболовых

15

и

13 -12 -II 10 9 -8

7 О 6 -

§180,!

Алсйкограниты О амфиболовые граниты О сиениты и граносиениты

0 А

. о

Расстояние от Большого Бештау, км

Рис. 3. Изотопный состав кислорода валовых проб гранитоидов КМВ в зависимости от расстояния от центральной части структуры КМВ (г. Большой Бештау)

гранитов, сиенитов, граносиенитов и лейкогранитов КМВ указывает на отсутствие изотопно-кислородного равновесия между минералами-вкрапленниками всех изученных типов пород. На рис. 4 нанесены рассчитанные температуры изотопного уравновешивания «сквозных» минералов-вкрапленников - флогопита и полевого шпата для четырех интрузивов КМВ. Полученные температуры ( 900°С) нереальны для условий кристаллизации гранитного расплава.

Результаты изучения изотопного состава кислорода и стронция явно указывают на определенную роль в истории формирования гранитоидов КМВ

ю.о

Рис. 4. Изотопный состав кислорода «сквозных» минералов-вкрапленников флогопита и полевого шпата в гранитоидах КМВ

процессов ассимиляции, которые, однако, не могли протекать в зоне становления интрузивных тел, вследствие низких концентраций стронция и неконтрастных величин 875г/865г во вмещающих интрузивные тела породах.

Вероятно, что ассимиляция могла иметь место еще до внедрения интрузивных тел, на более значительных глубинах. Информацию о подобных процессах могут предоставить обнаруженные в сиенитах и граносиенитах КМВ ксеногенные включения пород разного состава. Среди ксенолитов, присутствующих в интрузивных телах КМВ, доминируют включения габброидов и пород Ар-СагЬ-СРх-Р11^-Р1 состава. Последние заслуживают особого внимания, поскольку по своим геохимическим параметрам (высокие концентрации Бг (около 4000 ррт), Ва (около 1700ррт) и Р (около 8000 ррт)) они близки к «обогащенным» мантийным расплавам, с которыми обычно связывают образование высоко-Ва-8г пород. Ксеногенные включения этого типа встречаются практически во всех интрузивных массивах КМВ (исключением являются только амфиболовые граниты группы Бештау).

Проведенное изотопное исследование включения Ар-СагЬ-СРх-РЫо§-Р1 породы (образец 3-03-17) показало, что данный тип ксеногенных включений в гранитодах КМВ имеет скорее осадочное, а не магматическое происхождение. Макро- и микроскопически установлено, что включение 3-03-17 имеет зональное строение - в центре его преобладает кальцит, а по краям - преимущественно флогопит. В области контакта включения с вмещающими граносиенитами наблюдаются микропустоты, к которым приурочены выделения барита. Измеренный в образце 3-03-17 изотопный состав кислорода карбонатной фазы находится в интервале, характерном для измененных осадочных карбонатных пород ( |80=18.8%о), а изотопный состав серы барита из этого же образца (5345=13.5 %о) характерен для сульфатов древнего морского происхождения. Все изложенное указывает на то, что ксеногенные включения пород Ар-СагЬ-СРх-РЬ^-Р1 состава могут являться реститами захваченных расплавом выcoкo-Mg карбонатных пород, которые взаимодействовали с кислым расплавом по реакции типа:

Во! + ЦЛр, <2ь + Г + Н20) РМоц + £>/ + Р1 + Ар + СОг Т (1)

где Оо/ - ассимилянт, представленный высоко-М§ карбонатной породой, /. - кислый расплав. Высокие концентрации Бг, Ва, величина 875г/865г отношения (0.70867), близкая к интервалу начальных отношений Бг в гранитоидах КМВ, наличие карбонатной фазы и барита с «осадочными» изотопными метками, указывают на породы Ар-СагЬ-СРх-РЫо£-Р1 состава как на остаточный материал ассимилированных расплавами карбонатных пород.

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

происхождении включений Ар-СагЬ-СРх-РЫо§-Р1 состава. Изотопный состав серы барита в образце 3-03-17, предоставляет возможность оценить возрастные рамки, используя хемостратиграфическую кривую для сульфатов морского происхождения, а изотопный состав стронция, полученный для этого же образца, может быть сопоставлен с аналогичной кривой для стронция морской воды. Сопоставление полученных данных с возрастными кривыми, представляющими изменение величин :,43 сульфата морской воды и 875г/8б5г морской воды в течение фанерозоя, приведено на рис.5, пунктирная линия соответствует изотопным характеристикам образца 3-03-17.

СишПв ы | [лловй

j Списали | Tali*,

25

600 W

200

Age (Ma)

Рис. 5.а. Изотопный состав серы сульфата морской воды в фанерозое [Kump, 2003]

fSiKt

C«obri«n jOtdovicienj SiL j Der. jCeboraf. | Peim. |Тпм j Jurttnc jCreUceou» j Terüry "j'i

Ш m Ase(Ma) 200

Puc. 5.6. Распределение отношений 87Srf6Sr в морских карбонатах в фанерозое [Kump, 2003]

Из двух пересечений с возрастной кривой на рис. 5а, более реальной является граница триаса и юры, поскольку в геологическом разрезе осадочных отложений, известных в районе КМВ, более древних пород не наблюдается. Изотопный состав стронция морских осадков в это же время мог оказаться близким к измеренному в образце 3-03-17 (рис. 56). Именно в это время для изучаемого района было широко проявлено образование эвапоритовых толщ [Шевченко, 1992; Богашова, 2007 и др.], что могло обеспечить

геохимическую специфику ассимилянта (высокие

концентрации Бг и Ва).

Моделирование поведения изотопных систем в процессе контаминации Количественная оценка вклада ассимилянта была проведена в рамках модифицированной модели AFC(+X) [Ray, 1998], которая, в отличие от классической модели AFC [DePaolo, 1981], подразумевает, что на фоне ассимиляции происходит не только кристаллизация, но и обособление

дополнительной фазы X. В приложении к гранитоидам КМВ, под фазой X понималось образование «остаточных» фаз после взаимодействия высоко-К^ карбонатной породы с кислым (рБр+СНг) расплавом по реакции (1), аналогичных включениям Ар-СагЬ-СРх-РЬ^-Р1 состава. Под процессом кристаллизации понималось формирование ассоциаций минералов-вкрапленников (флогопита и полевого шпата). Все модельные расчеты описывают измеренные величины в области 1 > Г > 0.95 (где {- доля оставшегося расплава), т.е. когда количество вкрапленников в соответствии с реакцией (1) должно составлять от 0 до 8 %, что не превышает наблюдаемое содержание вкрапленников в гранитоидах КМВ.

Результаты расчетов по модели АРС(+Х) для Яг и N<1 систем показали (рис. 6), что формирование изотопных параметров и содержаний Бг и N(1 в гранитоидах КМВ может быть описано в рамках одностадийного процесса АБС(+Х) при исключительно коровых параметрах исходного кислого расплава без привлечения дополнительного мантийного компонента.

0.7095 0.7093 0.7091 - — 1 — 2

1 о грнигоидыКМВ • 3-03-17

0.7089 ■ 1

0.7087 - \ [Паль и др., 1993]

0.7085 0.7083 0.7081 ■ <Ä>° О 3-,ррш

1

0 2000 4000 6000 8000

ЕЩТ)

о/

I о — 1 .......... 2

О г{ВшпоишК1*В

Nd.ppm • 3-03-17

-1-1-

0 50 100 150 200 250 300 350

Рис. 6. Сопоставление расчета поведения изотопного состава и концентраций Sr и Nd по модели AFC(+X) [Ray, 1998] с измеренными значениями в гранитоидах КМВ и Ap-Carb-CPx-Phlog-Pl включении (образец 3-03-17)

Линии 1 и 2 - расчет изменения составов гранитоидного расплава и Ap-Carb-CPx-Phlog-Pl включений соответственно. Параметры, использованные при расчете - для исходного расплава: 87Sr/ 6Sr = 0.71;

Nd = -5; C(Sr)=320, C(Nd) = 27 ppm при DSf = 2.5 и Djvij = 0.1. Для высоко-Mg карбонатной породы, принятой в качестве ассимилянта: 87Sr/86Sr = 0.7083;

Nd = 0; C(Sr) = 3500 ppm и C(Nd) = 14 ppm. Коэффициент распределения «фаза X - расплав» для Sr принят равным 5.3, для Nd = 5.2

Количественные оценки поведения изотопно-кислородной системы были основаны на сравнительном анализе величин 5|80 «сквозных» минералов-вкрапленников - полевого шпата и флогопита в зависимости от изменения изотопного состава валовых проб, которые в случае порфировых пород с

небольшим содержанием вкрапленников могут быть отождествлены с составом расплава. Расчеты проводились таким образом, чтобы описать полученные для гранитоидов КМВ зависимости в координатах «6180 флогопит - 6180 вал» и «5 0 полевой шпат - 5180 вал» (рис. 7). При расчете 6|80 исходного расплава был принят равным составу амфиболовых гранитов Большого Бештау (+7%о), величина 5180 ассимилирующихся карбонатных пород принята равной 24.5 %о. Коэффициенты фракционирования «риолит-полевой шпат» и «риолит-флогопит» при Т = 750800 °С рассчитаны по уравнениям Жао и Ченга [Zhao, Zheng, 2003]. При расчетах допускалось, что формирование минералов-вкрапленников начинается одновременно с началом ассимиляции расплавом карбонатной породы, при этом кристаллизующийся вкрапленник находится в локальном равновесии с

расплавом. По мере протекания реакции (1) величина 5|80 расплава возрастает. Если допустить, что диффузионный обмен между растущим

минералом-вкрапленником и расплавом отсутствует, становится очевидным, что суммарный изотопный состав минерала будет постепенно отклоняться от равновесия кристалл-

расплав. Как показывают расчеты, только такой вариант поведения изотопно-кислородной системы способен объяснить наблюдаемые для гранитоидов КМВ линейные зависимости в координатах «5180 минерал - 5|80 вал», наклон которых отличается от 45 "(рис. 7). Расчеты показывают, что в процессе ассимиляции кислым расплавом карбонатных пород по реакции (1) могли сформироваться минералы-вкрапленники (ПШ и флогопит) с изотопным составом кислорода, неравновесным с изотопным составом расплава.

Формирование изотопно-геохимических параметров лейкогранитов КМВ

Среди гранитоидов КМВ группа лейкогранитов занимает особое положение. В истории образования лейкогранитов ведущая роль отводится процессам фракционной кристаллизации [Носова и др., 2005], а не процессам ассимиляции, как было предложено выше для остальных гранитоидов КМВ. Микроэлементный

Рис. 7. Расчет изменения ô180 вкрапленников полевого шпата и флогопита по мере возрастания ô'80 расплава (линии 1 и 2)

На диаграмму нанесены измеренные значения для четырех гранитоидных тел КМВ

состав лейкогранитов хорошо моделируется процессами фракционной кристаллизации из расплава, отвечающего составу амфиболовым гранитам группы Б. Бештау. Величины |80 амфиболовых гранитов КМВ находятся в интервале от 7 до 9.2 %о, в то время как лейкограниты имеют значения величины 5|80 от 8.9 до 13.8 %о. Объяснить такой существенный изотопный сдвиг процессами фракционной кристаллизации невозможно, в силу высоких температур и малых величин коэффициентов фракционирования кристалл-расплав. Процесс заметной ассимиляции вмещающих (боковых) пород в зоне остывания интрузивов, как уже было показано, не мог иметь место. Рассмотренный выше процесс контаминации до момента внедрения расплавов, хорошо описывающий поведение изотопных параметров сиенитов и граносиенитов, также не способен объяснить высокие значения величин 6180 в лейкогранитах. Остается предполагать, что наблюдаемый изотопно-кислородный сдвиг в лейкогранитах возникает в процессе взаимодействия флюид-порода на магматической стадии или на стадии остывания интрузивного тела.

5180, %о

$ Бык 9 Шелудивая

тт

гх

Рис. 8. Изменение изотопного состава кислорода по мере кристаллизации интрузивного тела

I- 6180 кристаллизующегося расплава, 2-6180 остывающей породы на периферии интрузива в момент закрытия изотопно-кислородной системы минералов ее основной массы

Для интрузивов лейкогранитов КМВ был проведен расчет по линейной модели

взаимодействия флюид-порода, изложенного в работе [Дубинина и др., 2009 в печати]. Схема взаимодействия подразумевает, что

остывание внедрившегося расплава приводит к движению фронта

кристаллизации от

периферийных частей интрузива к его центру, а

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

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

При расчетах было принято исходное содержание воды в расплаве, равное 5%, б|80 расплава = 7%о, Т=700 °С кристаллизации расплава и Т=300 °С - закрытия изотопно-кислородной системы полевых шпатов основной массы лейкогранитов. На рис. 8 показана динамика изменения величин 6180 кристаллизующегося расплава (линия 1), и минералов основной массы породы на периферии интрузива в момент закрытия их изотопно-кислородной системы (линия 2). Величина X, отложенная по оси абсцисс - массовая доля закристаллизованной породы в системе, т.е. доля оставшегося расплава в системе равна (1-Х). Проведенный расчет показывает, что измеренные величины 6180 в лейкогранитах массивов Бык, Шелудивая и Козьи Скалы могут являться следствием расположения отобранных образцов в периферийных зонах этих плохо вскрытых эрозией интрузивных тел (рис. 8). Из расчета следует, что изотопный сдвиг в процессе взаимодействия остывающей породы с «собственным» флюидом может быть весьма значителен на ранних стадиях кристаллизации интрузивного тела (х < 0.2). Далее изотопный сдвиг становится менее 1 %о, и ко времени кристаллизации половины расплава почти исчезает.

Глава 5. Происхождение вод минеральных источников КМВ

Процесс взаимодействия флюид-порода, отражается не только на изотопном составе пород, но и на изотопном составе флюида или природных вод, и его исследования можно проводить, изучая изотопные сдвиги как в породах, так и в природных водах, которые участвовали во взаимодействии с горными породами. Район КМВ, для которого известно большое количество теплых минеральных источников [Углекислые..., 1963], является уникальным объектом для исследований. В литературе имеется ограниченное количество данных по изотопному составу вод района КМВ, в основном, эти данные получены для проб, отобранных в 80-х годах прошлого века. Наши исследования проведены на пробах минерализованных и пресных вод (полевые работы 2003-2004 гг) спустя 25 лет, в течение которых шла интенсивная эксплуатация этих минеральных источников.

Измеренные нами величины Б и 180 в углекислых водах КМВ варьируют в довольно узком интервале: от (-79.2) до (-63.8) %о и от (-12.2) до (-10.4) %о соответственно. Изотопный состав водорода соответствует интервалу среднегодовых осадков, выпадающих в этом районе на высотах 1500 м и ниже (рис. 9), что подтверждает метеорное происхождение вод минеральных источников КМВ [Ляшенко, Потапов, 1984; Федоров, 1999; Яковлев, Поляк, 1997].

На рис. 9 видно, что полученные величины 5D и |80 в минерализованных и пресных водах КМВ, расположены немного левее как глобальной линии метеорных вод, предложенной Крейгом [Craig, 1961], так и локальной линии метеорных вод (JIJIMB), рассчитанной для района КМВ, с использованием базы данных МАГАТЭ. Составы вод в районе КМВ с таким же изотопно-кислородным сдвигом обнаруживались и ранее [Fedorov, 1989], это объяснялось протеканием процессов подземного испарения. Полярной точки зрения придерживается группа исследователей из ВСЕГИНГЕО [Поляков, Соколовский, 2005], которые предполагают сугубо инфильтрационный характер вод КМВ. Новые данные по геохимии микроэлементов и РЗЭ в водах района КМВ и Приэльбрусья, также свидетельствуют об активности процессов взаимодействия в системе «вода-порода» [Дубинина и др., 2005]. В этом случае положение точек слева от JIJIMB может быть следствием изотопно-кислородного обмена вод с силикатными и карбонатными породами при низкой температуре в условиях доминирования кислородного резервуара пород. Расчеты показывают, что изотопный состав кислорода воды, находящейся в равновесии с породами карбонатного (+20 %0) и силикатного (+10 %о) составов при температурах 70°С и 40°С (оценки температуры минерализованных вод КМВ на глубине 1000 м [Поляков, Соколовский, 2005]), может показать наблюдаемый сдвиг влево от JIJIMB только при взаимодействии с силикатными породами.

Рис. 9. Изотопный состав минерализованных и пресных вод района КМВ

Линия локальных метеорных вод и среднегодовые составы осадков рассчитаны с использованием базы данных МАГАТЭ [Во\уеп, Яеуепа^Ь, 2003; Во\уеп й а1., 2002; Вои^л е( а!., 2005]. Цифры около точек среднегодовых атмосферных осадков -абсолютные высоты, для которой произведен расчет.

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

5Э %о ® Среднегодовые атмосферные осадки

-Глобальная ЛМВ

----Локальная ЛМВ

О Минерализованные воды О Пресные источники О

А Речные воды ® ®0 ® °

Проведенные расчеты показали, что если вода, прошедшая стадию изотопного обмена с силикатными породами при низких температурах, в дальнейшем будет разбавляться локальными атмосферными водами, то можно получить составы, близкие к составам вод КМВ. Разбавляемый компонент при этом должен обладать пониженными величинами О. Сопоставление полученных нами данных по изотопному составу минерализованных вод района КМВ и с имеющимися в литературе [Поляков, Соколовский, 2005] показало, что таким компонентом могут являться талые ледниковые или атмосферные воды. На рис. 10 видно, что пересечение двух трендов, построенных по литературным и нашим данным, дает общую точку (А2) с координатами (-14.5; -93). Предполагаемый механизм смешения единого для всех компонента с разными типами поверхностных или захороненных вод позволяет рассматривать координаты этой точки как изотопные характеристики общего компонента. Его природа, несомненно, метеорная, поскольку при экстраполяции на ЛЛМВ (точка А1) получаются составы, соответствующие ледникам Большого Кавказа. Интерполяция с использованием базы данных МАГАТЭ показывает, что атмосферные осадки, имеющие такой изотопный состав водорода, должны выпадать на более высоких

о -ю -20 -30 -40 -50 -60 -70 -80 -90 -100 -110 -120

-20 -15 -10 -5 0 5

ЛЛМВ - линия локальных метеорных вод, поле «Н» - хлоридные воды Нагута, «Ч», «К» -усредненный состав вод Черного [Иконникова и др., 2007] и Каспийского [Поляков, Соколовский, 2005] морей, «СВК» - седиментационные воды флишевых толщ Северо-Западного Кавказа [Поляков, Соколовский, 2005]. 1- минеральные и пресные воды КМВ, данная работа, 2 -группа вод, в формировании которых принимает участие седиментогенный водный компонент [Поляков, Соколовский, 2005], 3 - ледниковые воды Большого Кавказа [Поляков, Соколовский, 2005], 4 - фигуративные точки инфильтрационного компонента в области питания (Al) и в зоне формирования вод КМВ (А2)

SD, %о

А2 -—Л*?

й'Ч),«.

Рис.10 Изотопный состав вод КМВ по данным настоящей работы и опубликованным данным [Поляков, Соколовский, 2005]

гипсометрических отметках, чем район КМВ - около 2400 м над уровнем моря.

Таким образом, областью питания минерализованных вод КМВ могут являться воды атмосферных осадков или тающих ледников, которые попадают в моноклинальные структуры осадочного чехла на склонах горных сооружений Большого Кавказа. Дальнейшая их фильтрация приводит к изотопно-кислородному сдвигу (А1—>А2 на рис. 10) за счет обмена с терригенными толщами, слагающими моноклинальные структуры осадочного чехла. Разбавление этих вод локальными атмосферными водами в зоне разгрузки источников КМВ приводит к уменьшению изотопно-кислородного сдвига.

Разбавление компонента А происходит не только локальными атмосферными водами, о чем свидетельствует группа вод Нагутского месторождения, точки которых расположены вдоль пунктира на рис. 10.

Объяснить изотопные параметры этих вод невозможно без привлечения вод седиментогенной природы - захороненных морских вод или рассолов древних эвапоритов. Наиболее контрастны в этом ряду хлоридные воды Нагута (поле «Н» на рис.10), для которых отмечается аномально высокая минерализация и специфический солевой состав [Поляков, Соколовский, 2005].

На сегодняшний день можно лишь сказать, что эти воды не соответстуют изотопному составу вод Черного и Каспийского морей, и не соответствуют полю седиментационных вод Северо-Запада Кавказа, захороненных во флишевых отложениях (поле «СВК» на рис.10).

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

Заключение

Проведенные изотопно-геохимические исследования гранитоидов КМВ показали, что все породы изученных интрузивных тел могли сформироваться из единого исходного расплава. Изотопно-кислородные данные коррелируют с геолого-петрологическим разделением пород интрузивных тел КМВ на три группы - амфиболовых гранитов, группу сиенитов и граносиенитов и группу лейкогранитов.

Процесс ассимиляции карбонатных пород кислыми расплавами, который обусловил формирование геохимических и изотопных особенностей сиенитов и

граносиенитов КМВ, по-видимому, протекал по механизму химического взаимодействия гранитного расплава с высоко-№^ карбонатными породами, содержащими эвапоритовый компонент. Реститами ассимилянта, вероятно, являются обнаруженные в сиенитах и граносиентах КМВ ксеногенные включения Ар-СагЬ-СРх-РИ^-Р1 пород. Проведенные изотопные (873/863г и 8 Б) исследования карбонатной фазы и барита во включении этого типа позволяют предполагать возраст ассимилянта, приходящийся на границу триаса и юры. В это время для Предкавказья характерно образование эвапоритовых толщ, что могло обеспечить геохимическую специфику ассимилянта - высокие концентрации стронция и бария.

На основе расчетов, проведенных по изотопным (Бг, N(1 и О) данным, показано, что формирование изотопных параметров и содержаний Бг и N<1 в гранитоидах КМВ может быть описано в рамках одностадийного процесса АРС(+Х) при исключительно коровых параметрах исходного кислого расплава без привлечения дополнительного мантийного компонента. Расчеты показывают, что в данном процессе могли сформироваться минералы-вкрапленники (ПШ и флогопит) с наблюдаемым изотопным составом кислорода, неравновесным с изотопным составом кислорода основной массы пород.

Формирование изотопно-кислородной системы лейкогранитов не подчиняется описанию процессами кристаллизационной дифференциации и ассимиляции боковых пород, и взаимодействием гранитного расплава с выcoкo-Mg карбонатными породами на глубине. Расчетным путем показано, что при формировании расплавов лейкогранитов КМВ, по-видимому, ведущую роль играло флюидно-магматическое взаимодействие, имевшее место на стадии остывания интрузивных тел.

Проведенные изотопные (80 и о180) исследования минеральных и пресных вод района КМВ показывают, что основным источником их водного компонента являются метеорные воды. Область питания удалена от выходов источников КМВ и находится на гипсометрических отметках около 2400 м над уровнем моря. Наблюдаемый дефицит 180 в водах КМВ может быть объяснен процессом взаимодействия инфильтрационного водного компонента с силикатными породами при низких температурах. Разнообразие в солевом и изотопном составе вод КМВ определяется процессами смешения метеогенных инфильтрационных вод либо с атмосферными водами в зоне выхода источников, либо с водами седиментационного генезиса.

Список работ, опубликованных по теме диссертации:

1. Авдеенко A.C., Дубинина Е.О. Изотопный анализ кислорода воды с применением метода уравновешивания с С02 в одноразовых шприцах // Тезисы докладов XVI Симпозиума по геохимии изотопов им. акад. А.П. Виноградова. Москва, 2001. С. 152.

2. Лаврушин В.Ю., Дубинина Е.О., Авдеенко A.C.. Костенко O.E. Углекислые воды Северного Кавказа: происхождение и условия формирования // Проблемы гидрогеологии XXI века: наука и образование. М.:РУДН , 2003, С. 431-443

3. Dubinina Е.О., Avdeenko A.S.. Lavrushin V. Yu. The genesis of the North Caucasian mineral waters: an isotope and geochemistry study // 14th V.M Goldschmidt conference: Copenhagen: Denmark, 2004, V. 68, Suppl. 1, A251

4. Лебедев B.A., Чернышев И.В., Авдеенко A.C.. Носова A.A., Докучаев А.Я, Олейникова Т.И., Гольцман Ю.В. Неуравновешенность начального изотопного состава Ar и Sr в сосуществующих минералах миоценовых гипабиссальных гранитоидов района Кавказских Минеральных Вод // Докл. РАН, 2006, Т.410, №1, С. 95-100

5. Авдеенко A.C.. Дубинина Е.О. Оценка эффекта "памяти" при измерении изотопного состава водорода в водных пробах с помощью приставки H/Device // Тезисы докладов XVIII Симпозиума по геохимии изотопов им. акад. А.П. Виноградова. Москва, 2007, С. 23

6. Авдеенко A.C.. Дубинина Е.О., Носова A.A., Гольцман Ю.В., Олейникова Т.И. Источник и эволюция расплавов позднемиоценовых гранитоидов района Кавказских Минеральных Вод // Тезисы докладов XVIII Симпозиума по геохимии изотопов им. акад. А.П. Виноградова. 2007, С. 24-25

7. Авдеенко A.C.. Дубинина Е.О., Носова A.A., Гольцман Ю.В., Олейникова Т.И. Изотопные (Nd, Sr, 5|80) характеристики как отражение природы протолита позднемиоценовых высоко-Sr-Ba гранитоидов района Кавказских Минеральных Вод (Северный Кавказ) // Докл. РАН, 2008, Т. 422, №. 1, С. 1 -5.

8. Авдеенко A.C.. Дубинина Е.О., Носова A.A. Неравновесность изотопного состава кислорода минералов-вкрапленников порфировидных гранитоидов КМВ (Северный Кавказ) // Тезисы докладов IV Российской конференции по изотопной геохронологии. С-Петербург, 2009, T.I, С. 17-20.

9. Дубинина Е.О., Носова A.A., Авдеенко A.C. О механизме контаминации расплавов гранитоидов карбонатными породами на примере высоко-Sr-Ba гранитоидов района КМВ (Северный Кавказ) // Тезисы докладов IV Российской конференции по изотопной геохронологии. С-Петербург, 2009, T.I, С. 169-171

Дубинина Е.О., Носова A.A., Авдеенко A.C., Аранович Л.Я. Формирование изотопных (Sr, Nd,0) параметров высоко-Sr-Ba гранитоидов позднемиоценовых интрузивов района КМВ (Северный Кавказ) // Петрология (в печати)

Подписано к печати 13.10.2009 г. Печать цифровая. Формат 60x90 1/16 Усл. печ. л. 1,5. Тираж 150 экз. Отпечатано в типографии «ИП Скороходов В.А.» Москва, Старомонетный пер., д.31 тел.(495) 950-30-39 Заказ N2 7902.

Содержание диссертации, кандидата геолого-минералогических наук, Авдеенко, Анна Сергеевна

Содержание

Введение

Глава 1. Происхождение и изотопно-геохимические характеристики высоко

Ba-Sr гранитоидов

1.1. Ba-Sr гранитоиды: положение в общей классификации гранитоидов и взгляды на происхождение

1.2. Существующие представления о происхождении высоко-Ba-Sr гранитоидов

КМВ

Глава 2. История изучения, геологическая, минералогическая и геохимическая характеристика гранитоидов района Кавказских Минеральных вод

2.1. История изучения района КМВ

2.2. Региональная геолого-текгоническая позиция КМВ

2.3. Основные этапы развития новейшего магматизма на Северном Кавказе

2.4. Геологическое строение района КМВ

2.5. Геологическая характеристика интрузивов КМВ

2.6. Особенности минерального состава гранитоидов КМВ

2.7. Геохимические особенности гранитоидов КМВ

2.7.1. Общая геохимическая характеристика гранитоидов КМВ

2.7.2. Вариации главных и малых элементов в гранитоидах района КМВ

2.8. Геодинамическая типизация гранитоидов КМВ

2.9. Характеристика ксенолитов гранитоидов района КМВ

2.10. Характеристика гидротермально-метасоматических образований

Глава 3. Методика изотопных исследований

3.1. Введение

3.2. Методика изотопных измерений и подготовки проб

3.2.1. Подготовка проб силикатного состава к изотопному анализу кислорода

3.2.1.1. Химическая основа метода

3.2.1.2. Описание экспериментальной установки и методических процедур

3.2.1.3. Методы контроля правильности анализа

3.2.2. Подготовка проб природных вод к изотопному анализу кислорода и водорода

3.2.2.1. Изотопное уравновешивание кислорода воды с СО2 в одноразовых шприцах 58 3.2.2.2. Изотопный анализ водорода воды методом разложения на металлическом хроме с использованием автоматической приставки

H/Device (ThermoQuest, Finnigan)

3.2.3. Методика определения изотопного состава Sr и Nd

3.2.3.1. Химическая подготовка проб

3.2.3.2. Измерение изотопного состава Rb -Sr и Sm - Nd отношений

Глава 4. Формирование изотопных (Sr, Nd, О) параметров высоко- Sr -Ва гранптоидов позднемиоценовых интрузивов района КМВ

4.1. Результаты изотопных исследований пород и минералов гранитоидов и вмещающих пород

4.1.1. Изотопные характеристики валовых проб гранитоидов КМВ

4.1.2. Изотопный состав кислорода минералов-вкрапленников

4.1.3. Изотопные параметры ксеногенных включений в граносиенитах

4.1.4. Изотопные параметры вмещающих пород

4.1.5. Сравнение полученных изотопных данных с неогеновыми магматическими породами Кавказского региона

4.2. Интерпретация полученных результатов

4.2.1. Изотопно-геохимические особенности формирования сиенитов, граносиенитов и лейкогранитов структуры КМВ

4.2.2. Формирование изотопно-геохимических параметров группы сиенитов и граносиенитов структуры КМВ

4.2.2.1. Вещество ксеногенных включений как индикатор процесса контаминации

4.2.2.2. Предполагаемый возраст контаминанта

4.2.2.3. Моделирование поведения изотопных систем в процессе контаминации

4.2.2.3.1. Поведение изотопов стронция и неодима

4.2.2.3.2. Поведение изотопно-кислородной системы

4.2.3. Формирование изотопно-геохимических параметров лейкогранитов 98 КМВ

4.2.4.1. Роль процессов дифференциации и ассимиляции

4.2.4.2. Флюидно-магматическое взаимодействие при формировании 102 изотопно-кислородных параметров лейкогранитов

Глава 5. Происхождение вод минеральных источников КМВ

5.1. Общие сведения по гидрогеологии района КМВ

5.2. Материалы исследований и результаты

5.3. Особенности изотопных параметров вод КМВ 112 Заключение 119 Библиографический список использованной литературы 120 Приложение№

Введение Диссертация по наукам о земле, на тему "Поведение изотопно-кислородной системы в процессах контаминации МАГМ и взаимодействия флюид-порода"

Актуальность темы. Гранитоидный магматизм проявляется на всех стадиях формирования континентальной коры, но наибольший масштаб гранитообразования характерен для пред-, син- и постколлизионных обстановок. С этим типом магматизма генетически связаны разнообразные (Cu-порфировые, Sn, W-Mo, редкометальные и др.) месторождения. В последние годы выявлена большая металлогеническая роль гранитов высоко-Ba-Sr (латитового) типа, высокодифференцированные разновидности которых (Li-F гранитоиды) обладают высокой рудоносностью с широким спектром месторождений -молибден-меднопорфировых, редкометальных, урановых.

Образование гранитов высоко-Ba-Sr типа на сегодняшний день связывают с пшомовым магматизмом [Тагпеу, 1994]. Для этих гранитов характерны высокие концентрации Ва, Sr и Р, низкие концентрации Rb, Th, U, Y, Nb и тяжелых редкоземельных элементов (HREE), а так же высокое отношение K/Rb. Источником таких пород могут быть «обогащенные» мантийные расплавы. Однако одной из проблем, возникающей при исследовании гранитов высоко-Ba-Sr типа, являются довольно радиогенные изотопные отношения Sr, умеренно радиогенные Nd и высокие значения величины 6180, которые не характерны для магматических пород, в образовании которых существенная роль отводится мантийным расплавам. Формирование изотопных параметров этого типа гранитов в основном объясняют процессами фракционной кристаллизации и ассимиляции осадочных пород «обогащенными» расплавами [Fowler, Henney 2001; Qian, 2003; Reiners, Nelson 1996]. Однако, этими процессами не всегда удается полностью описать поведение изотопных систем в высоко-Ba-Sr гранитах.

Данная работа посвящена исследованию проблемы формирования изотопных (Sr, Nd, О) параметров высоко-Ba-Sr гранитоидов. Для этого в качестве объекта исследования были выбраны высоко-Ba-Sr гранитоиды района Кавказских Минеральных Вод (КМВ), с которыми связано редкометальное и урановое оруденение. Изучение гранитоидов КМВ проводилось на основе исследования изотопной геохимии Sr, Nd и О с привлечением необходимых петрологических, геохимических и геологических данных. Комплексное изучение поведения изотопных систем (О, Sr, Nd,) позволяет решать широкий круг задач: исследование механизмов и установление степени контаминации магм, определение роли флюидной фазы в постмагматических процессах, решать ряд других вопросов.

Район Кавказских Минеральных Вод, в котором расположены высоко-Ba-Sr гранитоиды, отличается современной термальной активностью, выражающейся в многочисленных проявлениях деятельности минерализованных водных источников.

Параллельно с изучением самих гранитоидов КМВ было проведено изотопное (О, Н) изучение минеральных вод данного района с целью установления их типа питания и механизма формирования изотопных параметров. Полученные данные представляют интерес как для прогноза дальнейшего геохимического и гидрогеологического состояния уникального природного объекта - минеральных вод района КМВ, так и для решения общих проблем генезиса минерализованных вод.

Настоящее исследование имело следующие цели: (1) Определить роль процесса контаминации расплавов при формировании изотопно-кислородного состава высоко-Ва-Sr гранитоидных пород на примере гранитоидов КМВ; (2) Выявить роль взаимодействия водный флюид-порода при становлении гранитных интрузивных тел по данным изучения изотопно-кислородной системы гранитоидов КМВ; (3) Установить степень взаимодействия вода-порода при низких температурах на примере изучения изотопного состава водной фазы минеральных источников района КМВ.

Исходя из намеченных целей, были поставлены следующие задачи исследования: (1) Усоверш енствование методики фторирования для повышения точности измерения изотопного состава кислорода в валовых пробах и мономинеральных фракциях силикатных пород, модернизация метода изотопного уравновешивания для определения

1 Я

5 0 в природных водах; (2) Анализ изотопного состава кислорода в породах и минералах-вкрапленниках гранитоидных тел КМВ; (3) Анализ изотопного состава кислорода и водорода минерализованных и пресных вод района КМВ; (4) Интерпретация данных по изотопному составу кислорода и радиогенным изотопам (Sr, Nd) с использованием моделей, описывающих их поведение в процессах взаимодействия расплав-ассимилянт и вода-порода.

Фактический материал. В основу диссертации легли результаты исследований, проведенных автором в лаборатории изотопной геохимии и геохронологии ИГЕМ РАН. Образцы пород для изотопных исследований были отобраны при участии автора в течение двух полевых сезонов (2003-2004 гг.) в районе Кавказских Минеральных вод (КМВ), часть образцов любезно предоставлены А .Я. Докучаевым и А.А. Носовой (ИГЕМ РАН). Изотопные исследования сопровождались петрографическим и геохимическим исследованием пород.

Главные методы исследования. При исследовании гранитных и вмещающих пород района КМВ было проанализировано 23 образца пород Rb-Sr методом, 18 образцов пород Sm-Nd методом, в 39 образцах пород определен изотопный состав кислорода и в 14 образцах пород определен изотопный состав углерода. Для 18 минеральных фракций из гранитных пород изучался изотопный состав кислорода. Было выполнено 20 парных определений изотопного состава водорода и кислорода в образцах минерализованных вод района КМВ. Изотопные исследования пород сопровождались микроскопическим изучением образцов гранитных пород и микрозондовым исследованием составов минералов (геологический факультет МГУ, Н.Н. Коротаева), породообразующие оксиды в породах определялись методами химического анализа в ИГЕМ РАН (О.Г. Унанова). Редкие элементы в породах определялись методом ICP-MS в ИМГРЭ (Д.З. Журавлев) и методом инструментального нейтронно-активационного анализа (INAA) в ИГЕМ РАН (А.Л. Керзин).

Научная новизна работы. Впервые для гранитоидов КМВ проведено комплексное изотопно-геохимическое (О, Sr и Nd) и геохимическое изучение. Изотопный состав кислорода в породах КМВ изучен впервые. Показана принципиальная возможность преимущественно корового происхождения высоко-Ba-Sr гранитоидов с формированием их геохимической специфики в результате химического взаимодействия кислого расплава с высоко-Mg карбонатной породой. Проведены расчеты формирования изотопных параметров гранитоидов КМВ в рамках одностадийного процесса AFC(+X) без привлечения дополнительного источника мантийных расплавов. Предложен механизм флюидно-магматического взаимодействия, приводящего к формированию высоких величин 5180 в лейкогранитах. Подтверждено инфильтрационное происхождение минерализованных углекислых вод КМВ, впервые получена количественная оценка гипсометрической высоты области их питания.

Практическая значимость работы. С гранитами высоко-Ba-Sr типа часто связаны рудные месторождения: молибден-порфировые (Клаймакс), полиметаллические и золоторудные (Забайкалье (Приаргунье)), урановые (район КМВ), месторождения бора (Дальний Восток) и ряд других. Проведенные исследования условий формирования гранитов высоко-Ba-Sr типа дают информацию, которая может быть использована при построении генетических моделей рудообразования, связанного с позднеколлизионным магматизмом. Выявление роли взаимодействия флюид-порода и установление механизма контаминации магм, которым посвящена значительная часть работы, имеют большое значение для расшифровки источника и процессов транспорта рудного вещества.

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

Основные защищаемые положения

1. Впервые определен изотопный состав кислорода в валовых пробах и сквозных минералах-вкрапленниках гранитоидов КМВ. Установлены закономерности поведения изотопно-кислородной системы пород, коррелируклциеся с их минералого-петрологическим типом: для группы граносиенитов и сиенитов

1Я характерен узкий диапазон вариаций величины 5 О (от 9.1 до 10.4 %о), для группы

1 R амфиболовых гранитов диапазон вариаций величины 5 О шире и характерны в целом ее более низкие значения (от 7 до 9.1 %о), наибольший размах вариаций и наиболее высокие значения величины 6180 характерны для группы лейкогранитов (от 8.9 до 13.8 %о). Изотопно-кислородное равновесие в сквозных минералах-вкрапленниках всех исследованных типах гранитоидов отсутствует.

2. На основе комплексного изотопно-геохимического (О, Nd, Sr) и геохимического изучения гранитоидных пород структуры КМВ установлено, что в истории их формирования имел место процесс контаминации исходных расплавов осадочными породами на глубине. Поведение изученных изотопных систем может быть описано в рамках реакции кислого силикатного расплава с породами карбонатного состава.

3. На основе расчета по линейной модели изотопного обмена в процессе кристаллизации расплава показано, что высокие значения и значительный интервал

1 о вариаций величин 5 О в лейкогранитах КМВ могут быть следствием взаимодействия флюид-порода на стадии становления интрузивных тел в близповерхностных условиях.

4. Усовершенствованы методики изотопного анализа кислорода и водорода в образцах природных вод. С применением этих методик изучен изотопный состав минерализованных и пресных вод района КМВ. Показано, что основным компонентом минеральных источников КМВ являются инфильтрационные метеорные воды, испытавшие процесс изотопно-кислородного обмена с силикатными породами при низких температурах.

Структура и объем работы. Диссертация состоит из 5-ти глав, введения, заключения и дополнительных материалов, включающих результаты химических анализов пород (Приложение 1). Основной материал работы изложен на 133 страницах, которые содержат 15 таблиц и 50 рисунков. Список цитируемой литературы включает 125 наименований.

Апробация работы. По теме диссертации в рецензированных журналах опубликовано 3 статьи и одна сдана в печать. Результаты исследований докладывались на 3-х российских и одной международной конференции: XVI Симпозиуме по геохимии изотопов им. акад. А.П. Виноградова, (Москва, 2001); 14th V.M. Goldschmidt conference (Copenhagen, 2004); XVIII Симпозиуме по геохимии изотопов им. акад. А.П. Виноградова (Москва, 2007) и IV Российской конференции по изотопной геохронологии (Санкт -Петербург, 2009).

Благодарности. Автор выражает искреннюю благодарность своим научным руководителям Е.О. Дубининой и А.А. Носовой за советы и консультации в процессе работы над диссертацией. Автор благодарен заведующему лабораторией изотопной геохимии и геохронологии ИГЕМ член-корреспонденту РАН И.В. Чернышеву за постоянную помощь и ценные советы при написании работы. Автор благодарен Э.Д. Баировой, Т.И. Олейниковой и Ю.В. Гольцману за помощь в проведении изотопного анализа стронция и неодима. К.Н. Шатагину, А.В. Самсонову, Ю.О. Ларионовой и М.М. Аракелянц за ценные комментарии и полезные замечания, полученные при обсуждении результатов исследования. Особая благодарность А.В. Перкову, Н.И. Сердюку, Г.И. Сердюку и А.В. Леляеву, за помощь в усовершенствовании и поддержании работоспособности масс-спектрометрического оборудования лаборатории изотопной геохимии и геохронологии ИГЕМ РАН.

Автор признателен А.Я. Докучаеву за помощь в проведении и организации полевых работ, Л.В. Сазоновой за помощь в минеральных и петрологических исследованиях каменного материала, Ю.А. Костицыну (ГЕОХИ РАН) за ценные замечания при обсуждении отдельных положений диссертации.

Заключение Диссертация по теме "Геохимия, геохимические методы поисков полезных ископаемых", Авдеенко, Анна Сергеевна

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

Прежде всего, следует отметить, что методы исследований, применявшиеся в работе, были значительно усовершенствованы. В результате, изотопный анализ кислорода и водорода в образцах природных вод и кислорода в силикатах (валовых пробах и мономинеральных фракциях гранитоидов) был проведен с необходимой точностью. Новые методические разработки [Авдеенко, Дубинина, 2001] позволили проводить анализ изотопного состава кислорода природных вод с высокой экспрессностью, сохраняя при этом необходимую воспроизводимость и точность анализа. Данные методики могут быть рекомендованы к использованию в других лабораториях. К достоинствам методической части работы можно отнести детальный анализ алгоритма работы и источников ошибок [Авдеенко, Дубинина, 2007], возникающих при работе с новым устройством - H/Device для анализа изотопного состава водорода в водных пробах.

Наверное, самым значимым итогом данной работы является то, что формирование изотопных и геохимических параметров высоко-Ba-Sr типа гранитов можно описать без привлечения дополнительного мантийного компонента. Как показано на примере гранитоидов района КМВ, специфические особенности этого типа гранитов могли сформироваться в результате химического взаимодействия корового кислого расплава с высоко-Mg карбонатными породами, содержащими эвапоритовый компонент. Проведенное изотопно-кислородное изучение лейкогранитов (гранитов Li-F типа) района КМВ, с которыми связано урановое оруденение, позволило выявить в их формировании наличие процесса флюидно-магматического взаимодействия в процессе кристаллизации интрузивных тел.

Процесс взаимодействия флюид-порода исследовался в данной работе не только на основе изучения изотопного состава пород, но и на основе изучения изотопных параметров природных вод. Проведенное изучение изотопного состава минерализованных вод района КМВ показало, что изотопные характеристики вод КМВ можно объяснить процессом взаимодействия инфильтрационных метеорных вод с терригенными породами при низких температурах. Дальнейшее разбавление такого компонента локальными атмосферными водами описывает изотопные характеристики минерализованных вод района КМВ.

Библиография Диссертация по наукам о земле, кандидата геолого-минералогических наук, Авдеенко, Анна Сергеевна, Москва

1. Аракелянц М.М., Борсук A.M., Шанин JI.JI. Новейшая гранитоидная вулкано-плутоническая формация Большого Кавказа по данным калий-аргонового датирования//Докл. АН СССР.- 1968,-Т. 182. -№ 5. С. 1157-1160.

2. Богатиков О.А., Гурбанов А.Г., Коваленко В.И. и др. Верхнечегемский кальдерный комплекс на Северном Кавказе // Изв. АН СССР. Сер. геол.- 1992,- № 1.- С. 5-21.

3. Богашова Л.Г. Роль галогенновых вод в формировании месторождений полезных ископаемых// М.: ГЕОС, 2007. 168с.

4. Борсук A.M. Мезозойские и кайнозойские магматические формации Большого Кавказа / М.: Наука, 1979.- 300с.

5. Бубнов С.Н. Хронология извержений и источники расплавов новейших вулканических центров Большого Кавказа: Автореф. дис. . канд. геол.-минерал, наук: 25.00.04 / С.Н. Бубнов. М., 2003. - 27с.

6. Ветштейн В.Е. Изотопы кислорода и водорода природных вод СССР / Л.: Недра, 1982.-216с.

7. Волков Г.А. Некоторые вопросы гидрогеохимии процессов гипергенеза урановых месторождений Бештау и Бык: Дисс. . канд. геол.-минерал. наук / Г.А. Волков. М.: Экспедиция № 1 ИГЕМ АН СССР, 1964. - 270с.

8. Геологический Атлас Кавказских Минеральных Вод /Министерство природных ресурсов Российской Федерации федеральное государственное унитарное геологическое предприятие (ФГУГП) «Кавказгеолсъемка» // Главный редактор Н.И. Пруцкий, Ессентуки, 2003г.

9. Греков И.И., Компанией М.А., Литовко Г.В. и др. Геолого-геофизическая модель Эльбрусской вулканической области (Центр льный Кавказ). // Эволюциятектонических процессов в истории Земли. Новосибирск: СО РАН «Гео», 2004. -Т. 1.-С. 124-127.

10. Давыдова В. В. Петрография и условия образования субщелочных гранитоидов массива Бештау (Кавказские Минеральные Воды): Курсовая работа / В. В. Давыдова. МГУ: М., 2004. 72с.

11. Дейтерий и кислород-18 в подземных водах (масс-спектрометрические исследования) / Селецкий Ю.Б., Поляков В.А., Якубовский А.В., Исаев В.М. М.: Недра, 1973. - 144с.

12. Дубинина Е.О. Поведение стабильных изотопов при редкометальном рудообразовании на месторождении Акчатау: Автореф. дис. . канд. геол.-минерал. наук: 04.00.02 / М., 1993. - 24с.

13. Дубинина Е.О., Лаврушин В.Ю., Коваленкер В.А. и др. Изотопы кислорода и водорода в минеральных источниках Приэльбрусья // Геохимия. — 2005. №10. - С. 1078-1089.

14. П.Дубинина Е.О., Носова А.А., Авдеенко А.С., Аранович Л.Я. Формирование изотопных (Sr, Nd,0) параметров высоко- Sr -Ва гранитоидов позднемиоценовых интрузивов района КМВ (Северный Кавказ) // Петрология (в печати).

15. Елманова Н.М., Шолпо В.Н. Анализ взаимосвязи проявления углекислых вод с глубинными процессами в горно-складчатых регионах (на примере Кавказа) // Бюл. Моск. О-ва испытателей природы. Отд. Геол., 1981. -Т. 56. - вып. 5. - С. 118-130.

16. Есиков А.Д. Изотопная гидрология геотермальных систем. М.: Наука, 1989 — 208с.

17. Иванов Д.А., Бубнов С.Н., Волкова В.М. и др. Изотопный состав Sr и Nd в четвертичных лавах Большого Кавказа в связи с проблемой их петрогенезиса / Геохимия. 1993. - №3. - С.343-353.

18. Катастрофические процессы и их влияние на природную среду /Под ред. акад. Н.П. Лаверова. М.: ОИФЗ РАН, 2002. - Т. 1: Вулканизм. - 434с.

19. Короновский Н.В., Демина Л.И. Коллизионный этап развития Кавказского сектора Альпийского складчатого пояса: геодинамика и магматизм //Геотектоника.- 1999.-№2.- С. 17-35.

20. Короновский Н.В., Молявко В.Г., Остафийчук И.А. Петрохимические особенности и условия формирования неогеновых интрузивов района Кавказских Минеральных Вод // Изв. АН СССР. Сер. геол. 1986. - № 6. - С. 39-51.

21. Косарев B.C. Перспективы нефтеносности верхнеюрских отложений западной части Терско-Каспийского прогиба // Геология нефти и газа. 1982. - №2

22. Костицын Ю.А. Условия становления Эльджуртинского гранита по изотопным данным (кислород и стронций) в вертикальном разрезе // Геохимия. 1995. - №6. -С. 780-797.

23. Кучер М.И., Фридман А.И., Островский А.Б., Динабург В.Н. Изучение происхождения природных газов и минеральных вод изотопными методами (на примере Большого района КМВ) // Геол. и разведка. 1989. - № 2. - С. 53-59.

24. Лебедев В. А., Бубнов С. Н., Чернышев И. В., Гольцман Ю. В. Основной магматизм в геологической истории Эльбрусской неовулканической области (Большой Кавказ): К—Аг- и Sr-Nd-изотопные данные / Докл. АН. 2006. -Т. 406. - №1. - С. 78-82.

25. Лебедев, В.А., Чернышев И.В., Авдеенко А.С. и др. Неуравновешенность начального изотопного состава Аг и Sr в сосуществующих минералах миоценовых гипабиссальных гранитоидов района Кавказских Минеральных Вод // Докл. АН.-2006.-Т. 410.-№ 1.-С. 95-100.

26. Лебедев В.А., Чернышев И.В., Аракелянц М.М. и др. Геохронология неоген-четвертичного дацитового вулканизма северо-западной части Малого Кавказа (Грузия) // Стратиграфия. Геологическая корреляция. 2004.-Т.12.-№ 1,- С. 96-115.

27. Лебедев В.А., Чернышев И.В., Дудаури О.З.и др. Самсарский вулканический центр как очаг новейшего вулканизма на Малом Кавказе: К-Ar геохронологические и Sr-Nd изотопные данные // Докл. АН. 2003. - Т. 393. - № 6. - С. 802-808.

28. Ляхович В. В., Гурбанов А. Г. Геохимия и условия становления Эльджуртинского массива (Северный Кавказ) // Геохимия. 1992. - №6. - С.

29. Ляшенко С.И., Потапов Е.Г. Дейтерий и кислород-18 в подземных и поверхностных водах Большого района Кавказских Минеральных вод // Курортные ресурсы Северного Кавказа. Пятигорск, 1984. - С. 76-90.

30. Милановский Е.Е., Короновский Н.В. Орогенный вулканизм и тектоника Альпийского пояса Евразии.- М.: Недра, 1973. 280с.

31. Ниф Г. Масс-спектрометрический анализ простых соединений водорода // Успехи масс-спектрометрии. М., 1963. - С. 504-513.

32. Пантелеев И.Я., Масуренков Ю.П., Пахомов С.И. О происхождении СОг в подземных водах // Изв. АН СССР, Сер. Геол. 1962. - № 6. - С.95-98.

33. Поль И.Р., Хесс Ю.С., Кобер Б. и др. Происхождение и петрогенезис миоценовых трахириолитов (А-тип) из северной части Большого Кавказа // Магматизм рифтов и складчатых поясов. М.: Наука, 1993. - С. 108-125.

34. Поляков В.А., Соколовский Л.Г. Генезис и динамика минеральных вод Кавказа по результатам изотопно-геохимических исследований // Гидрогеология, инж. геология. Обзор / ООО «Геоинформмарк»: М., 2005. 65с.

35. Поляков В.А., Якубовский А.В., Селецкий Ю.Б., Исаев Н.В. Масс-спектрометрический метод измерений дейтерия и его использование в качестве природного индикатора для изучения генезиса подземных вод // Ядерная геофизика. М., 1972 - С. 229-237.

36. Природные изотопы гидросферы /Под редакцией В.И. Ферронского М.: Недра, 1975.-280с.

37. Розен О.М., Федоровский B.C. Коллизионные гранитоиды и расслоение земной коры.- М.: Научный мир, 2001.- 188с.

38. Сазонова JI.B., Носова А.А., Докучаев А.Я., Гурбанов А.Г. Латитовый тип позднеколлизионных гранитоидов (Северный Кавказ): геохимические и минералогические особенности // Докл. РАН.- 2003.- Т. 393.- № 2.- С. 1-5.

39. Смирнова М.Н. Минераловодская тектоно-магматическая кольцевая структура // Тектоника, геодинамика и процессы магматизма и метаморфизма: Материалы совещ. 2-5 февр. 1999 г. М.: ГЕОС, 1999. Т. II. - С. 153-155.

40. Соболев Н.Д., Лебедев-Зиновьев А.А., Назаров А.С. Неогеновые интрузии и домезозойский фундамент района Кавказских Минеральных Вод // Тр. ВИМС. Нов. сер. 1959. - Вып. 3. -211с.

41. Сомин М.Л. Структурная позиция и геодинамические условия формирования метаморфических комплексов Большого Кавказа и Кубы: Автореф. дис. .д-ра геол.-минерал. наук : 25.00.01 / М. Л. Сомин. М., 2007. - 56 с.

42. Станкевич Е.К. Новейший магматизм Большого Кавказа. Л.: Наука, 1976. - 231 с.

43. Суворова В.А., Дубинина Е.О. Новое устройство для конверсии в СОг кислорода с целью его изотопного анализа // Геохимия. 1994. - № 2. - С. 286-289

44. Сысоев А.А., Артаев В.Б., Кащеев В.В. Изотопная масс — спектрометрия. М.: Энергоатомиздат, 1993.-288с.

45. Таусон Л.В. Геохимические типы и потенциальная рудоносность гранитоидов. -М.: Недра, 1977.- 279с.

46. Тейлор С.Р., Мак-Леннан С.М. Континентальная кора: её состав и эволюция: пер. с анг. М.: Мир, 1988.-384с.

47. Углекислые минеральные воды Северного Кавказа. / Под ред. И.Я. Пантелеева. -М.: изд-во АН СССР, 1963. 190с.

48. Устинов В.И. Принципы и методы изучения изотопной неравновестности минеральных ассоциаций: Автореф. дис. . д-ра хим.наук: 04.00.02 / В. И. Устинов. -М., 1991.-50 с.

49. Федоров Ю.А. Стабильные изотопы и эволюция гидросферы. Гидрохимический институт Росгидромета российский фонд фундаментальных исследований: Истина, 1999. - 367с.

50. Ферронский В.И., Поляков В.А. Изотопия гидросферы. М.: Наука, 1983. - 280с.

51. Фор Г. Основы изотопной геологии: пер. с анг. М.: Мир, 1989. — 590с.

52. Ханель М., Липпольт X. Й., Кобер Б. и др. О раннепалеозойском возрасте метагранодиоритов в зоне Главного хребта Большого Кавказа // Петрология. Т. 1.- №5. - С. 487-.

53. Шевченко В.И. Условия образования верхнеюрских хемогенных отложений Предкавказья // Бюллетень МОИП. Отд. геологии. 1992. - Т. 67. - Вып. 2. - С. 104119.

54. Шемпелев А.Г. Неудавшиеся вулканы Пятигорья // Вулканизм и геодинамика. Материалы II Всероссийского симпозиума по вулканологии и палеовулканологии. 9-18 сентября 2003 г. Екатеринбург: Институт геологии и геохимии УрО РАН, 2003.-С. 331 -335.

55. Юра Кавказа/ К.О. Ростовцев, В.Б. Агаев и др.; под ред. К.О. Ростовцев.- СПб.: Наука, 1992. 192с.

56. Яковлев JI.E., Поляк Б.Г. Природа изотопно-гелевой аномалии в Северном Приэльбрусье // Вулканология и сейсмология. 1997. - №6. - С. 3-14.

57. Andersen D.J., Lindsley D.H. New (and final!) models for the Ti-magnetite/ilmenite geothermometer and oxygen barometer // Abstr. AGU 1985 Spring Meeting Eos Transactions. American Geophysical Union. 1985. - V. 66 (18). - P. 416.

58. Barnes C.G., Petersen S.W., Kistler R.W., Murray R., Kays M.A. Source and tectonic implication of tonalite-trondhjemite magmatism in the Klamath Mountains // Contrib. Mineral. Petrol. 1996. - Y.123. - P.40-60.

59. Bindeman I.N., Eiler J.M., Yogodzinski G.M. et al. Oxygen isotope evidence for slab melting in modern and ancient subduction zones // Earth. Plan. Sci. Lett. 2005. - V. 235 - P. 480—496.

60. Bottinga Y., Javoy M. Comments on stable isotope geothermometry: the system quartz-water // Earth. Plan. Sci. Lett. 1987. - V. 84. - P. 406-414.

61. Bowen G. J., Revenaugh J. Interpolating the isotopic composition of modern meteoric precipitation // Water Resources Research. 2003. - V. 39. -№10. P. 1299.

62. Bowen G. J., Wassenaar, L. I., Hobson K. A. Global application of stable hydrogen and oxygen isotopes to wildlife forensics // Oecologia, 2005.

63. Bowen G. J., Wilkinson B. Spatial distribution of 8180 in meteoric precipitation // Geology.-2002.- V. 30. №4. - P. 315-318.

64. Buket E., Temel A. Major-element, trace-element, and Sr-Nd isotopic geochemistry and genesis of Varto (Mus) volcanic rocks, Eastern Turkey / J. of volcanology and geothermal research. 1998. - V. 85. - P. 405-422.

65. Chappell B.W. and White A.J.R. Two contrasting granite types // Pacific Geol. 1974. -V. 8.-P. 173-174.

66. Clayton R.N., Mayeda Т.К. The use of bromine pentafluoride in the extraction of oxygen from oxides and silicates for the isotopic analysis // Geohim. Cosm. Acta. 1963. V. 27. - P. 43-52.

67. Cole D.R., Chakraborty S. Rates and mechanisms of isotopic exchange // Stable Isotope Geochemistry. Reviews in mineralogy and geochemistry. / Eds. J.W. Valley, D.R. Cole. -2001.-V. 43.-P. 83-224.

68. Coleman M.L., Thomas J.S., John J.D., Moore G.R. Reduction of water with zinc for hydrogen isotope analysis // Analytical Chemistry. 1982. - V. 54. - P. 993 - 995.

69. Collins W.J., Beams D.S., Whitea .J.R., Chappelbl W. Nature and origin of A-type granites with particular reference to southeastern Australia // Contrib. Mineral. Petrol. -1982.-V.80.-P. 189-200.

70. Craig H. Isotopic variations in meteoric waters // Science.-1961.-V.133. P. 1702-1703.

71. DePaolo D.J. Trace element and isotopic effects of combined wallrock assimilation and fractional crystallization // Earth. Planet. Sci. Lett. -1981. V. 53. - P. 189-202.

72. Dodson M.H. Closure temperature in cooling geochronological and petrological systems // Contrib. Miner. Petrol. 1973. - V. 40. - P. 259-274.

73. Drummond M.S., Defant M.J. A model for trondhjemite-tonalite-dacite genesis and crustal growth via slab melting: Archean to modern comparisons // J. Geophys. Res. -1990. V.95. - P. 21503-21521.

74. Dubinina E.O., Lakshtanov L.Z. A kinetic model of exchange in dissolution-precipitation processes// Geochim. Cosmochim. Acta. 1997. - V. 61. - P. 2265-2273.

75. Eby G.N. The A-type granitoids: a review of their Occurrence and chemical characteristics and speculations on their pedogenesis // Lithos.-1990.-V.26 P. 115-134.

76. Eiler J.M., Schiano P., Kitchen N., Stolper E.M. Oxygen isotope evidence for recycled crust in the sources of mid oceanic ridge basalts // Nature. 2000. - V.403. - P. 530-534.

77. Fedorov Yu. A. Isotopic-chemical description of natural waters at the large district of CMV // In abstract of 5th warking meeting «Isotopes in nature». Leipzig.- Academy of Sciences of the GDR, 1989. P.23.

78. Fowler M.B., Henney P.J. Mixed Caledonian appinite magmas: implications for lamprophyre fractionation and high Ba-Sr granite genesis // Contrib. Miner. Petrol. -1996.-V.126.-P. 199-215.

79. Fowler M.B., Henney P. J., Darbyshire D.P.F., Greenwood P.B. Petrogenesis of high Ba-Sr granites: the Rogart pluton, Sutherland // J. Geol. Soc. London. 2001. - V.158. -P. 521-534.

80. Gehre M., Hoefling R., Kowski P., Strauch G. Sample preparation device for quantitative hydrogen isotope analysis using chromium metal // Analytical Chemistry. 1996. - V. 68. - №. 24.-P. 4414-4417.

81. Gerdes A., Montero P., Bea F. et al. Peraluminous granites frequently with mantle-like isotope compositions: the continental-type Murzinka and Dzhabyk batholiths of the eastern Urals // J. Earth Sciences. 2002. - V. 91. - P. 3-19.

82. Godfrey J.D. The deuterium content of hydrous minerals from the East-Central Sierra Nevada and Yosemite National Park // Geohim. Cosm. Acta. -1962. V. 26. - P. 12151245.

83. Goldstein S.J., Jacobsen S.B. Nd and Sr isotopic systematics of rivers water suspended material: implications for crustal evolution // Earth. Planet. Sci. Lett. 1988. - V.87. - P. 249-265.

84. Jacobsen S.B., Wasserburg G.J. Sm-Nd evolution of chondrites and achondrites // Earth Planet. Sci. Lett. 1984. - V.67. - P.137-150.

85. Kump L. R. The Geochemistry of Mass Extinction // Treatise on Geochemistry/ Eds. H.D. Holland and K.K. Turekian. Elsevier-Pergamon, Oxford, 2003. - V.7. - P. 351363.

86. Leake В., Woolley A., Arps C. et al. Nomenclature of amphiboles: report of subcommittee on amphiboles of the International Mineralogical Association, commission on new minerals and mineral names // Canad. Mineralogist. 1997.-V.35 (l).-P. 219-246.

87. Le Bas M.J., Le Maitre R.W., Streckeisen A., Zanettin B. A chemical classification of volcanic rocks based on the total alkali-silica diagram // J. Petrol. 1986. - V.27. - P. 745-750.

88. MacCrea J.M. On the isotopic chemistry of carbonates and a paleotemperature scale / J.Chem.Phys. 1950. - V.18. - P.849-857.

89. Martin H. Adakite magmas: modern analogues of Archean granitoids // Lithos. 1999. -V.46.-P. 411—429.

90. Nelson D. R. Isotopic characteristics of potassic rocks: evidence for the involvement of subducted sediments in magma genesis // Lithos. 1992. - V. 28. - P. 403-420.

91. O'Neil J.R., Adami L.H., Epstein S. Revised value for the 180 fractionation between C02 and H20 at 25°C / J. Res. U.S. Geol. Surv. 1975. - V. 3. - P. 623-624.

92. Pearce J. A. Sources and settings of granitic rocks // Episodes.-1996.-V. 19,- P. 120-125.

93. Pearce J.A., Bender J.F., De Long S.E. et al. Genesis of collision volcanism in Eastern Anatolia, Turkey // J. Volcanol. Geotherm. Res. 1990. - V. - 44. - P. 189-229.

94. Peccerillo A., Taylor S.R. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey // Contrib. Mineral. Petrol. 1976. - V. 58. - № 1. - P. 63-81.

95. Philip Y., Cisternas A., Gvishiani A,, Gorshkov A. The Caucasus: An actual example of the initial stages of continental collision // Tectonophysics. 1989. - V.161. - P. 1-21.

96. Polyak B.G., Tolstikhin I.N., Yakovlev L.E. et all. Helium isotopes, tectonics and heat flow in the Norten Caucasus // Geohim. Cosm. Acta. 2000. -V. 64. - № 11. - P. 19251944.

97. Qian Q., Chung S.L., Lee T.Y., Wen D.J. Mesozoic high-Ba-Sr granitoids from North China: geochemical characteristics and geological implications // Terra Nova. 2003. -V.15. -№4. - P. 272-278.

98. Rare Earth element geochemistry / Developments in geochemistry. // Ed. P. Henderson. — Elsevier: Amsterdam, Oxford, 1984. -V.2 510pp.

99. Ray J.S. Trace element and isotope evolution during concurrent assimilation, fractional crystallization, and liquid immiscibility of a carbonated silicate magma // Geochim.Cosm.Acta. -1998. -V.62. -№19/20. -P.3301-3306.

100. Reiners P.W., Nelson B.K., Nelson S. W. Grustal Contamination of Magma from Compositionally Zoned Plutons and Associated Ultramafic Intrusions of the Alaska Range // J. Petrol. 1996. - V. 37. - №2. - P. 261-292.

101. Rudnick R.L., Gao S. Composition of the Continental Crust // Treatise on Geochemistry/ Eds. H.D. Holland and K.K. Turekian. Elsevier-Pergamon, Oxford, 2003.-V.3.-P.1-64

102. Rye R.O., Ohmoto H. Sulfur and carbon isotopes and ore genesis: A Review // Econ. Geol. 1974 - V. 69 - №.6. -P. 826-842.

103. Sharp Z. Principles of Stable Isotope Geochemistry //Pearson Prentice Hall, 2007 344 P.

104. Steni C.R., Kilian R. Role of the subducted slab, mantle wedge and continental crust in the generation of adakites from the Andean Austral Volcanic Zone // Contrib. Mineral. Petrol. 1996. - V.123. - P. 263-281.

105. Sun S.-s., McDonough W.F. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes // Magmatism in the Ocean Basins / Eds. A.D. Saunders, M.J. Norry Geol. Soc.: London, 1989. - № 42. - P. 313-345.

106. Tanweer A., Han L.F. Reduction of microliter amounts of water with manganese for D/H isotope ratio measurement by mass spectrometry // Isot. Environ. Health Stud. 1996. -V.32.-P. 97-103.

107. Tarney J., Jones C.E. Trace element geochemistry of orogenic igneous rocks and crustal growth models // J. Geol. Soc. London. 1994. - V.151. -P. 855-868.

108. Taylor H.P., Jr. Water/rock interactions and the origin of H2O in granitic batholits // J. Geol. Soc. London. 1977. - V. 133. - P. 509-558.

109. Taylor H.P., Jr. Oxygen and hydrogen isotope studies of plutonic granitic rocks // EPSL. 1978.-V. 38. -P. 177-210

110. Tobias H.J., Goodman K.J., Blacken C.E., Brenna J.T. High-precision D/H measurement from hydrogen gas and water by continuous flow isotope ratio mass spectrometry // Anal. Chem. 1995. - V. 67. P. 2486-2492.

111. Tsunogae Т., Osanai Y., Owada M. et al. High fluorine pargasite in ultrahigh temperature granulites from Tonagh Island in the Archean Napier Complex, East Antarctica // Lithos. -2003.-V. 70.-P. 21-38.

112. Valley J.W. Stable isotope thermometry at high temperatures //Stable Isotope Geochemistry. Reviews in mineralogy and geochemistry. / Eds. J.W.Valley, D.R. Cole. -2001.-V. 43.-P. 365-414.

113. Veizer J., Mackenzie F. T. Evolution of Sedimentary Rocks // Treatise on Geochemistry / Eds. H.D. Holland and K.K. Turekian. Elsevier-Pergamon, Oxford, 2003. - V.7. - P. 369-407.

114. Zhao Z.F., Zheng Y.F. Calculation of oxygen isotope fractionation in magmatic rocks // Chem. Geol. 2003. - V. 193. - P. 59-80.

115. Zheng Y.F. Calculation of oxygen isotope fractionation in anhydrous silicate minerals // Geochim. Cosmochim. Acta. 1993. - V. 57. - P. 1079-1091.

116. Zheng Y.F. Calculation of oxygen isotope fractionation in hydroxyl-bearing silicates // Earth. Plan. Sci. Lett. 1993. -V. 120. - P. 247-263.