Бесплатный автореферат и диссертация по наукам о земле на тему

Плотностные и геодинамические модели литосферы в зоне сочленения плиты Филиппинского моря и Азиатского континента в районе о-ва Тайвань

ВАК РФ 25.00.28, Океанология

Автореферат диссертации по теме "Плотностные и геодинамические модели литосферы в зоне сочленения плиты Филиппинского моря и Азиатского континента в районе о-ва Тайвань"

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

ГИЛЬМАНОВА Гульшат Забировна

ПЛОТНОСТНЫЕ И ГЕОДИНАМИЧЕСКИЕ МОДЕЛИ ЛИТОСФЕРЫ В ЗОНЕ СОЧЛЕНЕНИЯ ПЛИТЫ ФИЛИППИНСКОГО МОРЯ И АЗИАТСКОГО КОНТИНЕНТА В РАЙОНЕ о-ва ТАЙВАНЬ

Специальность 25.00.28 - океанология

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

Владивосток - 2004

Работа выполнена в Тихоокеанском океанологическом институте им. В.И. Ильичева и Вычислительном Центре Д В ОРАН Научные руководители:

доктор геолого-минералогических наук Р.Г. Кулинич. (гравитационное моделирование)

доктор геолого-минералогических наук, профессор Л.А. Маслов (геодинамическое моделирование) Официальные оппоненты:

доктор геолого-минералогических наук, профессор Ю.Ф. Малышев кандидат геолого-минералогических наук В.М. Никифоров

Ведущая организация: кафедра геофизики и геоэкологии Института инженерной и социальной экологии ДВГТУ

Защита диссертации состоится "24 " декабря 2004 г. на заседании диссертационного совета Д 005.017.02 при Тихоокеанском океанологическом институте им. В.И. Ильичева ДВО РАН.

С диссертацией можно ознакомиться в библиотеке ТОИ ДВО РАН

Отзывы в двух экземплярах, заверенные подписями и печатью учреждения, просим направлять по адресу: 690041, г. Владивосток, ул. Балтийская, 43, ТОИ ДВО РАН, ученому секретарю диссертационного совета.

Автореферат разослан "19 "ноября 2004 г.

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

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

к.г.н.

Ф.Ф. Храпченков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Окраинные моря и островные дуги составляют главные элементы переходной зоны от Тихого океана к континентальному обрамлению и характеризуются активным геодинамическим режимом Изучение причин и механизмов активного взаимодействия океанической и континентальной тектоносфер имеет длительную историю, но до сих пор остается одной из ключевых проблем в геолого-геофизических и геодинамических исследованиях В рамках этой проблемы зона сочленения северо-западного угла плиты Филиппинского моря с континентальной окраиной Азии в районе о-ва Тайвань, где сходятся островодужные системы Рюкю и Лусон, является одним из наиболее сложных звеньев всей Западно-Тихоокеанской переходной зоны. Здесь активизирована тектоническая и магматогенная деятельность, носящая порой экстремальный характер В силу этого указанный район является одним из важных и привлекательных объектов для изучения процессов конвергентного сочленения океанической и континентальной тектоносфер Его исследованию посвящено большое количество работ, однако многие аспекты его глубинного строения и природы геодинамических событий еще во многом не ясны (рис 1)

20е

22°

24°

26'

Рис. 1. Схема расположения района исследования. Линии АВ, А1В1.А2В2-профили структурно-плотностного и геодинамического моделирования

11в°

120°

122°

124°

4-24* _____

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

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

Задачи исследований. Для реализации поставленной цели сформулированы следующие задачи:

1. Построение структурно-плотностных моделей земной коры и подкоровой литосферы, на основе гравитационного моделирования, по профилям, пересекающим зону сочленения Филиппинской плиты с окраиной континента в области островодужных систем Рюкю и Лусон; исследование особенностей структурно-шютностной дифференциации изучаемой зоны.

2. Расчет полей напряжений и оценка вклада в современное геодинамическое состояние района плотностных неоднородностей, выделенных по результатам структурно-плотностного моделирования.

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

Защищаемые положения:

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

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

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

4. По результатам гравитационного и геодинамического моделирования Рюкю-Тайваньская зона сочленения Филиппинской плиты с окраиной континента, к западу и востоку от подводного хребта Гагуа, имеет различные структурно-плотностные и геодинамические характеристики. Предполагается, что хребет Гагуа является звеном трансструктурного тектоно-магматического линеамента, протягивающегося вдоль 123° в.д., сыгравшего существенную роль в развитии местных геодинамических процессов и формировании структурного плана изучаемого района.

Научная новизна полученных результатов состоит в следующем:

1. Впервые построены детальные структурно-плотностные модели литосферы по серии профилей, пересекающих зону сочленения плиты Филиппинского моря с окраиной континента в пределах островодужных систем Рюкю и Лусон-Тайвань, которые являются основой для различных геологических, геодинамических и тектонических построений.

2. Установлена плотностная латеральная дифференциация литосферы. По различию структурно-плотностных и геодинамических характеристик выделены три типа литосферных блоков: океанический, островодужный и окраинно-континентальный.

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

4. По результатам гравитационного и геодинамического моделирования Рюкю-Тайваньская зона сочленения Филиппинской плиты с окраиной континента, к западу и востоку от подводного хребта Гагуа, имеет различные структурно-плотностные и геодинамические характеристики. Предполагается, что хребет Гагуа является звеном трансструктурного тектоно-магматического линеамента, протягивающегося вдоль 123° в.д., сыгравшего существенную роль в развитии местных геодинамических процессов и формировании структурного плана изучаемого района.

5. Выполнены расчеты внутренних напряжений и движений вещества литосферы, обусловленных установленными плотностными неоднородностями. Показано, что объемные силы, создаваемые собственными плотностными неоднородностями литосферы, вносят значимый вклад в ее современное геодинамическое состояние.

5. Выполнены расчеты внутренних напряжений и движений вещества литосферы, обусловленных установленными плотностными неоднородностями. Показано, что объемные силы, создаваемые собственными плотностными неоднородностями литосферы, вносят значимый вклад в ее современное геодинамическое состояние.

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

Исходные материалы и личный вклад автора. В работе использованы геофизические данные, полученные ТОИ ДВО РАН в притайваньском районе Филиппинского, Восточно - Китайского и Южно-Китайского морей в экспедициях на НИС «Академик А. Несмеянов» (1989 г.) и «Профессор Гагаринский» (1993-1994 г.г.). Комплекс исследований включил непрерывное сейсмическое профилирование (НСП), гравиметрию и магнитометрию. Автор настоящей работы принимала участие в измерениях, а также в обработке и интерпретации полученных гравиметрических и магнитометрических данных. В работе использованы многочисленные опубликованные данные российских и зарубежных авторов, в частности, результаты сейсмических работ, представленных в работе (Wang K.T. et al., 2001). Обработка, анализ и интерпретация результатов моделирования выполнены автором самостоятельно. Структурно-плотностное моделирование коры и подкоровой литосферы выполнено с использованием программ, разработанных в лаборатории гравиметрии ТОИ ДВО РАН (Т.Н. Колпащикова) и лаборатории региональной геологии и геофизики ИТиГ ДВО РАН (к.г.-м.н. В.Я. Подгорный). Расчет напряженно-деформированного состояния коры и верхней мантии выполнен д. г.-м.н. Л.А. Масловым (ВЦ ДВО РАН) в системе MAPLE по программам, разработанным им лично и Д.А. Сырниковым.

Апробация работы и публикации. Результаты выполненной работы докладывались на двух международных конференциях: "New Concepts in Global Tectonics", Otero Junior College, La Junta, Colorado, May 5-11, 2002, США и IUGG-2003 (International Union of Geodesy and Geophysics), Саппоро, Япония; на IV Косыгинских чтениях: "Тектоника,

ВЦ ДВО РАН г. Хабаровск, на заседании ученого совета ИТИГ ДВО РАН г. Хабаровск и отдела геологии и геофизики ТОИ ДВО РАН.

По теме исследований опубликовано 12 работ.

Структура и объем работы. Диссертационная работа состоит из Введения, трех глав и Заключения, содержит 107 страниц текста, 38 рисунков, список литературы из 182 наименований. Работа выполнена в лаборатории гравиметрии ТОИ ДВО РАН под руководством д.г.-м.н. Кулинича Р.Г. и лаборатории математического моделирования в геологии и геофизике ВЦ ДВО РАН под руководством д.г.-м.н. Маслова Л.А. Автор благодарен сотрудникам лаборатории гравиметрии, магнитометрии и сейсмометрии ТОИ ДВО РАН: Т.Н. Колпащиковой, СМ. Николаеву, Б.Я. Карпу, М. Валитову, П. Зимину, вместе с которыми прошла школу морской геофизики, и которые оказывали большую помощь в процессе сбора, обработки и интерпретации геолого-геофизической информации. Положительному завершению работы автор обязана творческим и деловым контактам и просто дружескому отношению большого числа геофизиков, геологов и других специалистов. Большую поддержку оказал чл.-корр. СИ. Смагин. Выражаю искреннюю благодарность к. г.-м. н. В. Я. Подгорному, оказавшему существенную помощь в процессе плотностного моделирования и анализа полученных результатов. Большой вклад в осознание геологических и тектонических проблем исследуемого региона, а также истолкование полученных результатов внес д.г.-м. н. В.Г. Варнавский.

СОДЕРЖАНИЕ РАБОТЫ

Глава 1 носит обзорный характер, в ней приводится краткий обзор истории геолого-геофизических исследований Филиппинского моря и общая геолого - геофизическая характеристика района исследований. Основной вклад в развитие теоретических моделей формирования и эволюции Филиппинского моря внесли Д. Кариг, С. Уеда, Ц. Бен-Аврахам, К. Мрозовский, Д. Хейес, Т. Хилде, Т.Сено и С. Маруяма, X. Токуяма, А.Г. Родников, Е.А. Константиновская и многие другие. Регулярные геолого-геофизические исследования Филиппинского моря до конца 90-х годов прошлого столетия проводились и учеными Дальневосточного Отделения АН СССР (в настоящее время ДВО РАН): Васильевым Б.И., Кулиничем Р.Г., Пущиным И.К., Евлановым Ю.Б., Карпом Б.Я., Леликовым Е.П., Точилиной СВ., Сьединым В.Т., Маляренко А.Н., Щека С.А., Сваричевским А.С, Мельниченко Ю. И., (ТОИ). Исследуемый район в целом рассматривается как область конвергентного сочленения океанической и континентальной литосферных плит и характеризуется сложной тектонической и

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

Глава 2 посвящена описанию гравитационного (структурно-плотностного) моделирования земной коры и подкоровой литосферы. Реализация этой задачи была выполнена вдоль трех профилей, два из которых пересекают островодужную систему Рюкю, а третий - островную дугу Лусон (см. рис. 1, профили АВ, А1В1 и А2В2). В тексте диссертации указанные профили получили наименования: АВ - Субмеридиональный, -А1В1 - Меридиональный и А2В2 - Субширотный.

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

Моделирование было выполнено в четырех структурных вариантах: а) литосфера, как однослойная толща; б) литосфера, разделенная на земную кору и подкоровую часть; в) литосфера с земной корой, разделенной на консолидированную ее часть и осадочно-вулканогенную толщу; г) литосфера с земной корой, расчлененной на сейсмические слои. Значения плотности вычислялись по известным зависимостям скорости упругих продольных волн от плотности среды их прохождения (Nafe J. Е. & Drake С. С, 1963 и др.). Далее в главе приводится подробное описание и анализ результатов структурно-плотностного моделирования, выполненного по каждому профилю.

Структурно-плотностная модель Субмеридионального профиля (АВ). Начинаясь в Восточно-Китайском море профиль пересекает трог

Окинава, островную дугу Рюкю, впадину Нанао и хребет Яэяма, расположенные на южном склоне островной дуги, желоб Рюкю и выходит в Западно-Филиппинскую котловину (см. рис. 1). Крайняя северная часть профиля,

расположенная в пределах Восточно-Китайского моря, проходит над Синцзи-Тайваньским складчатым поясом. Профиль проходит в той части островодужной системы Рюкю, где ее простирание начинает меняться с юго-западного на широтное. Это происходит в зоне меридиана 123° в.д. Разрез вдоль профиля охватывает литосферу и верхнюю часть астеносферного слоя (рис.2). Выявлен ряд структурно-плотностных особенностей, имеющих как общий, так и локальный характер. Среди них наиболее существенными являются следующие:

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

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

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

Рис.2 Плотностная модель литосферы по Субмеридиональному профилю АВ с элементами интерпретации.

Условные обозначення:1-гравитационная аномалия от вычисленного распределения плотности; 2- измеренные аномалии в свободном воздухе; 3* значения вычисленной плотности; 4-границы вычисленных плотностных неоднородностей; 5-11- элементы интерпретации: 5-границы между выделенными блоками литосферы: ВКБ-Восгочно-Китайским, ОДБ • Островодужным, ЗФБ - Западно-Филиппинским; б- то же, выделенные неуверенно; 7- неоднородности пониженной плотности первого порядка; 8- неоднородности пониженной плотности второго порядка; 9-неоднорсшности повышенной плотности первого порядка; Ю- неоднородности повышенной плотности второго порядка.11- границы выделенных неоднородностей повышенной и пониженной плотности.

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

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

Структурно-плотностная модель Меридионального профиля (А1В1). Северный конец профиля располагается в районе северного борта трога Окинава 30' В.Д.). Пересекая указанный трог, дугу Рюкю, впадину Нанао, хребет Яэяма и желоб Ркжю (см. рис. 1), он заканчивается в пределах котловины Гагуа (22е 10* С.Ш. И 122е 30* в.д.). Профиль располагается западнее уже упоминавшегося меридиана в

существенно иной, более активной структурно - и сейсмотектонической обстановке.

На уровне литосферных и коровых масс здесь удалось выделить ряд существенных особенностей плотностной дифференциации (рис.3).

1. По средним значениям и характеру латерального распределения плотности литосфера разделена на два блока - Западно-Филиппинский и Островодужный. Ранее выделенный Восточно-Китайский блок здесь отсутствует из-за недостаточной протяженности на север рассматриваемого профиля. Западно-Филиппинский блок включает в себя котловину Гагуа, Островодужный - трог Окинава, дугу Рюкю, бассейн Нанао и хребет Яэяма. Четкой границы между выделенными блоками не существует. Переход от Островодужного блока к Западно-Филиппинскому происходит в широкой зоне, включающей бассейн Нанао и хребет Яэяма. Латеральное изменение плотности литосферы здесь повторяет тренд, зафиксированный в пределах Субмеридионального профиля: ее понижение от Западно-Филиппинского блока в сторону континентальной окраины.

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

В пределах коры выделяется ряд зон с переменной плотностью, имеющих один общий признак: их общее погружение в северном направлении. В результате создается впечатление о поддвигании структур плиты Филиппинского моря под островодужную систему и континентальную окраину. По этому признаку коровый разрез данного профиля кардинально отличается от разреза предшествующего профиля, где выявлено налегание островодужных структур на окраину континента. Этот фактор рассматривается в качестве одного из признаков существенного различия структурно-тектонического и динамического состояния участков, располагающихся под Меридиональным и Субмеридиональным профилями. Граница между этими разнородными участками проходит примерно вдоль 123° в.д. В зоне этого меридиана протягивается и подводный хребет Гагуа. Помимо этого здесь же обнаружена пространственная корреляция зон высокой концентрации землетрясений с высокоградиентными зонами латерального изменения плотности и мощности структур в коре и мантии.

Мантийное основание Западно-Филиппинского блока неоднородно. Этот блок подстилается астеносферой с минимальной для всего профиля плотностью. Под котловиной Гагуа из астеносферы в литосферу восстает вертикальная зона разуплотнения вещества, имеющая свое продолжение и в коре этого бассейна. Над ее выходом к поверхности морского дна зафиксирована локальная положительная аномалия магнитного поля. На стыке котловины Гагуа и желоба Рюкю, там, где основание океанической коры погружается под островодужную систему, в ее консолидированной части хорошо выделилась субвертикальная зона разуплотнения, уходящая своими «корнями» в мантию. Вероятно, эта аномальная зона маркирует участок начальной фазы разрушения сплошности литосферы на стыке Филиппинской плиты с окраиной континента. Наличие здесь активного разломообразования неоднократно предполагалось многими предшественниками

Структурно-плотностная модель Субширотного профиля (А2В2). Профиль располагается в южной части исследуемого района и протягивается практически в широтном направлении от Западно-Филиппинской котловины (21° 26' С.Ш., 123" 31* в.Д.)

через Лусон-Тайваньскую островодужную систему до северо-восточной окраины Южно-Китайского моря (21° 15' С.Ш., 119° 53' В.Д.). Профилем пересекаются: подводный хребет Гагуа, котловина Гагуа, Северо-Лусонский хребет, Северо-Лусонский трог, хребет Хенчун и северное окончание Манильского желоба (см. рис. 1).

Результаты гравитационного моделирования вдоль профиля показаны на рис.4. Анализ полученных данных позволил выделить в литосфере и ее коровой части ряд существенных структурно-плотностных особенностей.

1. По средним значениям и характеру распределения плотности литосфера профиля делится на три блока - Южно-Китайский (Приконтинентальный), Островодужный и Западно-Филиппинский. По среднему уровню плотностей Островодужный и Южно-Китайский блоки различаются незначительно. Их большее различие проявляется в характере внутреннего распределения этого параметра.

2. Внутренние структурно-плотностные характеристики выделенных блоков имеют индивидуальные особенности. В мантии и коре Южно-Китайского блока на стыке с островодужной Лусон-Тайваньской системой фиксируется сквозная субвертикальная зона разуплотнения, уходящая «корнями» в астеносферу. Она представляется своеобразным апофизом более обширной области понижения плотности подкоровой части литосферы Островодужного. Это свидетельствует о преимущественно деструктивном характере тектонических процессов, формировавших это звено изучаемой конвергентной зоны.

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

Северо-Лусонский хребет на всех видах плотностных разрезов коры и подкоровой литосферы также выделяется как объект, сложенный массами с относительно пониженной плотностью. Плотность раннемиоценовых - плейстоценовых вулканогенных пород (Huang С. Y., et al., 1992), слагающих его вершинную часть, находится в полном соответствии с расчетными данными. Несмотря на указанную плотностную характеристику и достаточно большую мощность земной коры (около 19 км) этот хребет

Юж н о • Китайское

мор е Филиппинское море

•аиъ ХаЯкаяьский ЗапаДНО-ФиЛНППННСКа* КОТЛОВИНА

нг* Семро* Северо-

Маннльский хребет Лусонский Лусонский бассейн хребет желоб Ненчун трог хребет Гжгуа Гагуа

П2

и

Рис 4 Плотностная модель литосферы по Субширотному профилю А2В2 с элементами интерпретации.

Условные обозначения I-гравитационная аномалия от вычисленного распределения плотности, 2- измеренные аномал ни.в свободном воздухе, 3- значения вычисленной плотности, 4-фанииы вычисленных плотнострых неоднородностей, 5-11- элементы интерпретации 5-границы между выделенными блоками литосферы ЮКБ-Южно-Китайским, ОДБ • Островодужным, ЗФБ - Западно-Филиппинским, б- то же, выделенные неуверенно, 7- неоднородности пониженной плотности первого порядка, 8- неоднородности пониженной плотности второго порядка, 9-неоднородности повышенной плотности первого порядка,Ю- неоднородности повышенной плотности второго порядка 11- границы выделенных неоднородностей повышенной и пониженной плотности

вызывает значительные по амплитуде возмущения в аномальном гравитационном поле (около 120 мГал). Это дает основание предполагать его изостатическую неуравновешенность. Ранее А. Уоттс (Watts et all., 1980), а затем Р.Г. Кулинич (Karp et al., 1997) уже указывали, что этот хребет является структурой, некомпенсированной изостатически, а его топографические массы поддерживаются внешними сжимающими усилиями. Результаты выполненного моделирования подтверждают данное заключение.

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

Подводный хребет Гагуа является одной из самых примечательных структур Западно-Филиппинского блока. Мощность консолидированной коры под ним составляет около 10 км, что сопоставимо с мощностью коры северо-западной части Тихого океана (Karp et al., 1997). С другой стороны, указанная величина существенно отличает хребет Гагуа от других асейсмичных океанических хребтов, где мощность коры значительно превышает указанную. Кора конформно деформирована в виде двух асимметричных, падающих на запад складок. Это сопровождается локальным подъемом границы Мохо под западным склоном хребта до глубин 7-8 км от уровня моря. В целом внутренняя структура этого сооружения значительно шире ее морфологического проявления на поверхности. В вещественном составе коры хребта Гагуа преобладают мафиты, основную часть ее объема составляет третий океанический слой, мощность которого здесь больше, чем в пределах Западно-Филиппинской котловины. Ранее перечисленные структурно-вещественные особенности хребта послужили основой для заключения о том, что эта структура сформирована поднятием мафических и ультрамафических масс низов коры и верхней мантии в результате тектонического скучивания (Karp et al., 1997). Кроме того, хребет Гагуа является границей между районом Западно-Филиппинской котловины с тонкой корой на востоке и котловиной Гагуа с более мощной и деформированной корой на западе. Помимо этого, на северном продолжении хребта, примерно вдоль меридиана 123° в.д., проходит граница между участками с существенно различной сейсмичностью в притайваньском районе конвергентной зоны. Обнаруженная под хребтом сквозная зона уплотнения вещества коры и мантии свидетельствует, что хребет Гагуа есть внешнее

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

Глав» 3 содержит описание процедур геодинамического моделирования, выполненного на основе вышеописанных структурно-плотностных моделей, и анализ полученных результатов.

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

C.К. Ранкорн, И.К. Туезов, В.П. Болдырев, А.В. Ефимов, Л.И. Лобковский, W.M. Elsasser,

D.L. Turcotte, J.M. Vanpe и др.

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

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

Геодинамическая модель Субмеридионального профиля (АВ). Модель представлена векторным полем скоростей смещения вещества и изолиниями касательных деформаций в коре и подкоровой мантии (рис. 5а). По характеру распределения рассчитанных параметров профиль делится на три блока, аналогичных тем, что были выделены в результате структурно-плотностного моделирования: Западно-Филиппинский, Островодужный и Восточно-Китайский. Основной геодинамической характеристикой Западно-Филиппинского блока является незначительный подъем коровых масс и опускание всего мантийного субстрата, скорость которого с глубиной увеличивается.

Западная граница этого режима образует наклонную криволинейную поверхность, погружающуюся от верхних горизонтов коры Западно-Филиппинской котловины вблизи желоба Рюкю под островодужную мантию. Вероятно, указанный нисходящий процесс связан с релаксацией мантийного выступа под Западно-Филиппинской котловиной и латеральным перетеканием глубинного вещества. Касательные деформации в этом блоке фиксируются в зоне контакта плиты Филиппинского моря с желобом Рюкю и, видимо, маркируют восточную границу зоны ее конвергенции с окраиной континента. Геодинамический режим Восточно-Китайского блока прямо противоположен предшествующему. Здесь верхняя часть коры (до глубин 6-7 км), находится в режиме опускания. При этом горизонтальная составляющая векторов скорости смещений направлена в сторону трога Окинава. Большая часть нижележащей коры и мантия находятся в режиме подъема. Максимальные скорости этого процесса наблюдаются севернее трога Окинава. В целом режим подъема коромантийных масс прекращается на границе с островным хребтом Рюкю, уступая место режиму погружения. Верхняя граница области восходящих движений погружается под указанную островодужную систему. Пространство, занимаемое этой областью, включает в себя зону разуплотнения, вещества коры и мантии выявленной по результатам структурно-плотностного моделирования на стыке Островодужного и Восточно-Китайского блоков. В пределах выделенного блока участки максимальных касательных деформаций зафиксированы в зоне разломов Тунгхай ограничивающих с севера трог Окинава. Островодужная система (блок), располагающаяся между вышеописанными блоками, резко отличается от них, прежде всего, наличием плотной сети касательных напряжений и отсутствием единого режима скоростей и направлений смещения вещества коры и мантии. Верхняя часть коры (до глубин 6-7 км), подстилающая трог Окинава и его непосредственное обрамление, находится в режиме опускания, при этом максимальные скорости погружения приурочены к его центральной части, где располагаются наиболее молодые структуры растяжения - грабены. Касательные напряжения сосредоточены вдоль обоих бортов указанного трога и совпадают не только с известными системами разломов, ограничивающих эту морфоструктуру, но и с зонами высоких градиентов мощности коры, определяя тем самым области наиболее вероятного разломообразования. Все это, помимо других известных геолого-геофизических данных, свидетельствует о существовании здесь активных деструктивных процессов и продолжающемся формировании трога Окинава. Островной хребет Рюкю находится в режиме опускания, а хребет Яэяма и внутридуговая впадина Нанао - в режиме подъема. При этом указанный хребет подвергнут более активному процессу указанного направления. Коровые массы

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

Результаты выполненных расчетов в пределах данного профиля показали, что собственные плотностные неоднородности коры и мантии, их релаксация могут создавать напряжения до 1010 Па (при ее вязкости 102<> Па-с). Как известно, этого достаточно, для образования вертикально ориентированных разломов или ослабленных зон. На этой основе делается вывод, что рассмотренные плотностные неоднородности вносят значимый вклад в геодинамику района и могут быть одним из структурообразующих факторов. Распределение областей максимальных касательных напряжений подчиняется схеме, соответствующей теоретическим расчетам, а именно -максимальные касательные напряжения концентрируются в областях максимального градиента мощности слоя.

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

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

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

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

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

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

По характеру пространственного распределения указанных характеристик в пределах профиля выделяются три блока, аналогичные вышеописанным: Западно-Филиппинский, Островодужный и Южно-Китайский (аналог Восточно-Китайского). Западно-Филиппинский блок выделяется по общему погружению мантийного субстрата и незначительному подъему коровых масс. Границей между этими режимами служит поверхность Мохо. Указанное геодинамическое состояние наблюдается под всей глубоководной частью рассматриваемого района. Западная граница данного блока имеет сложный характер. В коровой части она проходит вдоль зоны сочленения котловины Гагуа с Северо-Лусонским хребтом, а затем, погружаясь, уходит с переменным углом наклона под указанный хребет. Внутреннее геодинамическое состояние разных участков рассматриваемого блока неоднородно. Хребет Гагуа, как и остальная часть коры котловины, находится в режиме подъема. Здесь наблюдается высокая концентрация касательных напряжений, свидетельствующих о существовании в

пределах данной морфоструктуры сколовых деформаций и соответствующих разломов. В поле изолиний горизонтальных напряжений хребет Гагуа проявляется как структура сжатия, что соответствует предшествующим косвенным оценкам. Наконец, по разные стороны от указанного хребта скорости подъема коровых, а частично и подкоровых масс, различны: к востоку от него этот режим проявлен более интенсивно, чем западнее. Таким образом, обнаруживается еще один признак барьерной роли хребта Гагуа в геодинамике и тектонике изучаемого района. Южно-Китайский блок выделяется, прежде всего, по режиму восходящих движений вещества в мантии и коре. Его условную границу с островодужной системой можно провести от западной периферии Манильского желоба в пределах верхней части коры вниз под хребет Хенчун. Геодинамический режим внутри блока однороден. При общем режиме подъема коромантийных масс здесь можно отметить лишь интенсификацию этого процесса в направлении от глубинных мантийных горизонтов в верхние слои коры под Южно-Китайским морем. В пределах блока отсутствуют также участки значительных касательных и горизонтальных напряжений и деформаций. Лишь в мантии под западным склоном хр. Хенчун и соответствующим бортом Манильского желоба фиксируется локальный участок касательных напряжений, возникновение которых вполне объяснимо на границе указанных радикально разнородных морфоструктур и соответствует теоретическим расчетам (Маслов, 1997). Между выделенными Западно-Филиппинским и Южно-Китайским блоками располагается Островодужный блок. Основной отличительной особенностью его геодинамического режима являются: интенсивные нисходящие движения в верхней части коры и практически статичное состояние подстилающей мантии и низов коры. Блок характеризуется широким развитием касательных и горизонтальных напряжений, сосредоточенных в коре на стыке разнородных морфоструктур. На этом фоне наиболее интенсивными нисходящими движениями коровых масс характеризуются хребты Северо-Лусонский и Хенчун. Разделяющий их Северо-Лусонский трог характеризуется менее интенсивными нисходящими движениями. Это «отставание», на наш взгляд, может быть одной из причин возникновения касательных напряжений в зонах соприкосновения указанного трога с пограничными хребтами. Вдоль западного и восточного бортов Северо-Лусонского хребта фиксируются зоны геодинамического сжатия, за счет которых, возможно, этот хребет удерживается в неуравновешенном состоянии. Однако этих сил, видимо, недостаточно для сохранения такого состояния и указанный хребет, опускаясь, постепенно приближается к своему изостатическому равновесию. Наибольшая концентрация сжимающих усилий приходится на Северо-Лусонский трог. Причиной этого, видимо, является боковое давление со стороны смежных хребтов. Численные

значения возникающих здесь напряжений также достигают что достаточно

для формирования разломных или ослабленных зон.

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

1. Выделенные первоначально по признакам структурно-плотностной дифференциации три литосферных блока: Западно-Филиппинский, Островодужный и Окраинно-континентальный (на профиле АВ - Восточно-Китайский, на профиле А2В2 -Южно-Китайский) выделены также и по различию геодинамического режима. Из этого следует важное заключение: коромантийные неоднородности являются одним из источников геодинамических напряжений, деформаций и перемещения масс в указанной среде.

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

3. Конфигурация глубинных границ литосферных блоков, выделенных на обоих профилях, свидетельствует о некотором поддвигании мантийного субстрата Западно-Филиппинского блока под островные дуги Рюкю и Лусон-Тайвань.

4. К существенным различиям рассматриваемых профилей можно отнести неодинаковый объем глубинных «геодинамических корней» островодужных систем. В частности, область островодужных геодинамических характеристик (концентрация касательных напряжений и минимум смещений мантийных масс) под дугой Рюкю гораздо больше, чем под дугой Лусон-Тайвань. Кроме этого, максимальной концентрацией касательных напряжений также характеризуется островодужная система Рюкю. Из этого следует, что активность геодинамических процессов в пределах этой системы существенно выше, нежели в зоне дуги Лусон-Тайвань.

ЗАКЛЮЧЕНИЕ

В результате выполненных исследований построены структурно-плотностные и геодинамические модели литосферы по трем профилям, пересекающим зону сочленения плиты Филиппинского моря и Азиатского континента в районе острова Тайвань. Основные результаты исследований сводятся к следующему:

1. Установлена структурно-плотностная дифференциация литосферы. По различию средних плотностей и особенностям латерального распределения плотностных неоднородностей в коре и мантии выделены три типа литосферных блоков: океанический, островодужный и окраинно-континентальный.

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

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

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

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

источников геодинамических напряжений, деформаций и перемещения масс в указанной среде.

5. На всех профилях зафиксирован режим нисходящего движения мантии Западно-Филиппинского блока; напротив, окраинно-континентальная литосфера испытывает подъем. Островодужные системы Рюкю и Лусон-Тайвань различаются высокой концентрацией касательных напряжений и деформаций на стыках разнородных морфоструктур. Основной причиной этого, на наш взгляд, является пограничное расположение островодужных систем между литосферными плитами с разнонаправленным движением слагающих их масс. Основной объем мантийного вещества под островодужными системами во всех случаях находится в нейтральном состоянии, без значимых перемещений.

6. Конфигурация глубинных структурно-плотностных и «геодинамических» границ литосферных блоков свидетельствует о некотором поддвигании мантийного субстрата Западно-Филиппинского блока под островные дуги Рюкю и Лусон-Тайвань.

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

По теме диссертации опубликованы следующие работы:

1. Кулинич Р.Г., Маслов Л.А., Гильманова Г.З., Комова О.С. Плотностная модель и напряжения в земной коре северной части Японского моря // Тихоокеанская геология. -1998. - Том 17, № 2. С. 108-114 (обработкаматериала, моделирование, участие в интерпретации результатов).

2. Гильманова Г.З., Маслов Л.А., Кулинич Р.Г., Комова О.С. Геодинамика земной коры зоны коллизии Тайваня // Тихоокеанская геология. - 2002. - Том 21, № 3. С. 40-49 (обработка материала, моделирование, участие в интерпретации результатов.)

3. Гильманова Г.З., Подгорный В.Я. Плотностная модель литосферы зоны перехода от Филиппинской плиты к Азиатской, на примере островной дуги Рюкю // IV Косыгинские чтения: Сб.докладов. - Хабаровск, 2003. С. 238-258 (обработка материала, построение разрезов, интерпретация результатов, доклад).

4. Подгорный В.Я., Гильманова Г.З. Плотностная модель литосферы зоны перехода от Филиппинской плиты к Азиатской на примере островной дуги Рюкю // Материалы 30-ой сессии Международного семинара им. Д.Г.Успенского под редакцией академика

В.Н.Страхова: Тез. докл. - Москва: ОИФЗ РАН, 2003. С. 95-96 (обработка материала, построение разрезов, интерпретация результатов).

5. Gilmanova G., Maslov L., Kulinich R. Density model and crustal stresses in the Collision zone ofTaiwan // New concepts in global tectonics: Proceedings

book, 2002. P. 213-224 (обработка материала, моделирование, участие в интерпретации результатов).

6. Кагр В., Gilmanova G., Kulinich R. Complexities in the tectonic evolution of the eastern Taiwan offshore area // 30th International Geological

Congress: Abstract - Vol. 2. Beijing, China, 1996. P.224 (обработка материала, построение разрезов, интерпретация результатов).

7. Gilmanova G., V. Podgorny. Results of density modeling ofthe lithosphere in

the Philippine - to Asian plate transition zone, as exemplified from the Ryukyu Island arc. // Proc. of IUGG. Vol. 2. - Japan, 2003. P. ВЗО (обработка материала, построение разрезов, интерпретация результатов, доклад).

8. Gilmanova G., R. Kulinich, L. Maslov. Structure and geodynamics ofthe

Earth's crust and upper mantle of Japan and Philippine marginal seas according to the complex ofgeological and geophysical data // Proc. oflUGG. Vol. 2. - Japan, 2003. ВЗО (обработка материала, построение разрезов, интерпретация, доклад).

9. Podgorny V., Gilmanova G., Gravitational model ofthe structure ofthe lithosphere in the zone of the Asian and Philippine plates interaction, as exemplified from the Luzon island arc // Proc. of 32nd International Geological Congress IGC. - Florence - Italy, 2004. Abstracts (part 1). P. 607 (обработка материала, построение разрезов, интерпретация результатов).

10. Gilmanova G., Maslov L., Podgorny V. Geodynamic and density modeling ofthe lithosphere ofthe collision zone ofTaiwan // Proc. of 32nd International Geological Congress IGC. -Florence - Italy, 2004. Abstracts (part 2). P. 1461. (обработка материала, построение разрезов, интерпретация результатов).

11. Maslov L., Gilmanova G. Mechanical Stresses in the Lithosphere ofthe Pacific Tectonic Belt // Proc. ofWestern Pacific Geophysics Meeting. - Honolulu, 2004. P. 194. (обработка материала, интерпретация результатов).

12. Maslov L.A., Gilmanova G.Z., Simekov D.A., The stress field and structure-forming processes caused by density and topography variation in the Pacific continent-ocean transition zone// Proc. of AGU Fall Meeting. - San Francisco, 2004. T43C-1343 1340h. (обработка материала, интерпретация результатов)

ГИЛЬМАНОВА Гульшат Забировна

ПЛОТНОСТНЫЕ И ГЕОДИНАМИЧЕСКИЕ МОДЕЛИ ЛИТОСФЕРЫ В ЗОНЕ СОЧЛЕНЕНИЯ ПЛИТЫ ФИЛИППИНСКОГО МОРЯ И АЗИАТСКОГО КОНТИНЕНТА В РАЙОНЕ о-ва ТАЙВАНЬ

(специальность 25.00.28 - океанология)

Автореферат диссертации на соискание ученой степени кандидата геолого-минералогических наук

Подписано к печати 29.10.04 г. Формат 60x84/16

Печать офсетная. Усл. изд. л. 0,8. Тираж 100 экз. Заказ № 46.

Отпечатано в ОНТИ ТОИ ДВО РАН. 690041, Владивосток, Балтийская, 43

»2430»

290

Содержание диссертации, кандидата геолого-минералогических наук, Гильманова, Гульшат Забировна

Введение.

Глава 1. Геолого-геофизическая характеристика района исследований.

1.1. Общие сведения.

1.2. Геологическая характеристика основных морфоструктур района исследований.

1.2.1. Западно-Филиппинская котловина.

1.2.1.1. Подводный хребет Гагуа.

1.1.2.2. Котловина Гагуа.

1.2.2. Островодужная система Рюкю.

1.2.2.1. Желоб Рюкю.

1.2.2.2. Островодужная гряда (хребет) Рюкю.

1.2.2.3. Трог Окинава.

1.2.3. Остров Тайвань.

1.2.4. Островодужная система Лусон.

1.2.5. Южно-Китайское море.

1.3. Геофизическая характеристика основных морфоструктур района исследований.

1.3.1. Гравитационные аномалии в свободном воздухе.

1.3.2. Аномалии геоида.

1.3.3. Аномалии магнитного поля.

1.3.4. Тепловой поток.

1.3.5. Электропроводность.

1.3.6. Сейсмичность.

1.3.7. Современные движения.

1.4. Нефтегазоносность района исследований.

1.5. Выводы.

Глава 2. Гравитационное моделирование.

2.1. Состояние проблемы моделирования.

2.2. Построение модели первого приближения, использованные методики моделирования.

2.3. Структурно-плотностная модель Субмеридионального профиля (АВ).

2.4. Структурно-плотностная модель Меридионального профиля (А1В1).

2.5. Структурно-плотностная модель Субширотного профиля (А2В2).

2.6. Выводы.

Глава 3. Геодинамическое моделирование.

3.1. Современные методы расчетов напряжений в литосфере.

3.2. Методика расчетов смещений и деформаций.

3.3. Оценка влияния вариаций плотности в разрезе на магнитуды компонент тензора напряжений.

3.4. Геодинамическая модель Субмеридионального профиля (АВ).

3.5. Геодинамическая модель Меридионального профиля (А1В).

3.6 Геодинамическая модель Субширотного профиля (А2В2).

3.7 Выводы.

Введение Диссертация по наукам о земле, на тему "Плотностные и геодинамические модели литосферы в зоне сочленения плиты Филиппинского моря и Азиатского континента в районе о-ва Тайвань"

Актуальность темы.

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

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

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

Настоящая работа посвящена структурно-плотностному моделированию и изучению напряженно-деформированного состояния коры и литосферы в зоне конвергенции плиты Филиппинского моря с окраиной Азиатского континента в районе острова Тайвань. В морфоструктурном плане этот район включает западную часть Западно-Филиппинской котловины, юго-западную оконечность островодужной системы Рюкю и северную оконечность Лусон-Тайваньской островной дуги (рис. 1,2).

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

Выполненная работа базировалась на следующих теоретических и методологических предпосылках. Зона сочленения Филиппинской плиты с окраиной континента, независимо от какой-либо концептуальной позиции, объективно является устойчивой системой конвергентного сочленения океанской и континентальной литосферы. Формирование и эволюция данной системы должна приводить к ее существенной плотностной неоднородности. Гравитационные силы, которые возникают при неоднородном распределении масс, помимо других геотектонических факторов, могут создавать значительный вклад в напряженное состояние среды и, как следствие, существенно влиять на формирование геолого-структурной обстановки. В условиях активных геодинамических процессов возможно нарушение изостатического равновесия и усиление гравитационной неустойчивости в тектоносфере. В результате этого должны возникать дополнительные источники напряжений и их соответствующая реализация. С этой позиции исследования должны разделяться на два этапа, на первом из них необходимо выполнить структурно-плотностное моделирование на базе инструментальных гравиметрических данных и другой геолого-геофизической информации, а затем - рассчитать напряженно-деформированное состояние геологической среды. В последнем случае решается прямая задача математической физики с объемными силами, вызванными плотностными неоднородностями.

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

Задачи исследований. Для реализации поставленной цели сформулированы следующие задачи:

1. Построение структурно-плотностных моделей земной коры и подкоровой литосферы, на основе гравитационного моделирования, по профилям, пересекающим зону сочленения Филиппинской плиты с окраиной континента в области островодужных систем Рюкю и Лусон; исследование особенностей структурно-плотностной дифференциации изучаемой зоны.

2. Расчет полей напряжений и оценка вклада в современное геодинамическое состояние района плотностных неоднородностей, выделенных по результатам структурно-плотностного моделирования.

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

Защищаемые положения:

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

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

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

4. По результатам гравитационного и геодинамического моделирования Рюкю-Тайваньская зона сочленения Филиппинской плиты с окраиной континента, к западу и востоку от подводного хребта Гагуа, имеет различные структурно-плотностные и геодинамические характеристики. Предполагается, что хребет Гагуа является звеном трансструктурного тектоно-магматического линеамента, протягивающегося вдоль 123° в.д., сыгравшего существенную роль в развитии местных геодинамических процессов и формировании структурного плана изучаемого района.

Научная новизна полученных результатов состоит в следующем:

1. Впервые построены детальные структурно-плотностные модели литосферы по серии профилей, пересекающих зону сочленения плиты Филиппинского моря с окраиной континента в пределах островодужных систем Рюкю и Лусон-Тайвань, которые являются основой для различных геологических, геодинамических и тектонических построений.

2. Установлена плотностная латеральная дифференциация литосферы. По различию структурно-плотностных и геодинамических характеристик выделены три типа литосферных блоков: океанический, островодужный и окраинно-континентальный.

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

4. По результатам гравитационного и геодинамического моделирования Рюкю-Тайваньская зона сочленения Филиппинской плиты с окраиной континента, к западу и востоку от подводного хребта Гагуа, имеет различные структурно-плотностные и геодинамические характеристики. Предполагается, что хребет Гагуа является звеном трансструктурного тектоно-магматического линеамента, протягивающегося вдоль 123° в.д., сыгравшего существенную роль в развитии местных геодинамических процессов и формировании структурного плана изучаемого района.

5. Выполнены расчеты внутренних напряжений и движений вещества литосферы, обусловленных установленными плотностными неоднородностями. Показано, что объемные силы, создаваемые собственными плотностными неоднородностями литосферы, вносят значимый вклад в ее современное геодинамическое состояние.

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

Исходные материалы и личный вклад автора. В работе использованы геофизические данные, полученные ТОЙ ДВО РАН в притайваньском районе Филиппинского, Восточно - Китайского и Южно-Китайского морей в экспедициях на НИС «Академик А. Несмеянов» (1989 г.) и «Профессор Гагаринский» (1993-1994 г.г.). Комплекс исследований включил непрерывное сейсмическое профилирование (НСП), гравиметрию и магнитометрию. Автор настоящей работы принимала участие в измерениях, а также в обработке и интерпретации полученных гравиметрических и магнитометрических данных. В работе использованы многочисленные опубликованные данные российских и зарубежных авторов, в частности, результаты сейсмических работ, представленных в работе (Wang К.Т. et al., 2001). Обработка, анализ и интерпретация результатов моделирования выполнены автором самостоятельно. Структурно-плотностное моделирование коры и подкоровой литосферы выполнено с использованием программ, разработанных в лаборатории гравиметрии ТОЙ ДВО РАН (Т.Н. Колпащикова) и лаборатории региональной геологии и геофизики ИТиГ ДВО РАН (к.г.-м.н. В.Я. Подгорный). Расчет напряженно-деформированного состояния коры и верхней мантии выполнен д. г.-м.н. JI.A. Масловым (ВЦ ДВО РАН) в системе MAPLE по программам, разработанным им лично и Д.А. Сырниковым.

Апробация работы и публикации. Результаты выполненной работы докладывались на двух международных конференциях: "New Concepts in Global Tectonics", Otero Junior College, La Junta, Colorado, May 5-11, 2002, США и IUGG-2003 (International Union of Geodesy and Geophysics), Саппоро, Япония; на IV Косыгинских чтениях: "Тектоника, глубинное строение и геодинамика востока Азии", г. Хабаровск, 2003; на конференции "Закономерности строения и эволюции геосфер", г. Хабаровск 2003; на семинарах лаборатории математического моделирования в геологии и геофизике и ученого совета ВЦ ДВО РАН г. Хабаровск, на заседании ученого совета ИТИГ ДВО РАН г. Хабаровск и отдела геологии и геофизики ТОЙ ДВО РАН.

По теме исследований опубликовано 12 работ.

Структура и объем работы. Диссертационная работа состоит из Введения, трех глав и Заключения, содержит 107 страниц текста, 38 рисунков, список литературы из 182 наименований. Работа выполнена в лаборатории гравиметрии ТОЙ ДВО РАН под руководством д.г.-м.н. Кулинича Р.Г. и лаборатории математического моделирования в геологии и геофизике ВЦ ДВО РАН под руководством д.г.-м.н. Маслова JI.A. Автор благодарен сотрудникам лаборатории гравиметрии, магнитометрии и сейсмометрии ТОЙ ДВО РАН: Т.Н. Колпащиковой, С.М. Николаеву, Б.Я. Карпу, М. Валитову, П. Зимину, вместе с которыми прошла школу морской геофизики, и которые оказывали большую

Заключение Диссертация по теме "Океанология", Гильманова, Гульшат Забировна

3.7. Выводы.

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

1. Выделенные первоначально по признакам структурно-плотностной дифференциации три литосферных блока: Западно-Филиппинский, Островодужный и Окраинно-континентальный (на профиле АВ - Восточно-Китайский, на профиле А2В2 -Южно-Китайский) выделены также и по различию геодинамического режима. Из этого следует важное заключение: коромантийные неоднородности являются одним из источников геодинамических напряжений, деформаций и перемещения масс в указанной среде.

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

3. Конфигурация глубинных границ литосферных блоков, выделенных на обоих профилях, свидетельствует о некотором поддвигании мантийного субстрата Западно-Филиппинского блока под островные дуги Рююо и Лусон-Тайвань и налегании островодужных коровых масс на окраину смежного континента.

4. К существенным различиям рассматриваемых профилей можно отнести неодинаковый объем глубинных «геодинамических корней» островодужных систем. В частности, область островодужных геодинамических характеристик (концентрация касательных напряжений и минимум смещений мантийных масс) под дугой Рююо гораздо больше, чем под дугой Лусон-Тайвань. Кроме этого, максимальной концентрацией касательных напряжений также характеризуется островодужная система Рююо. Из этого следует, что активность геодинамических процессов в пределах этой системы существенно выше, нежели в зоне дуги Лусон-Тайвань.

131

ЗАКЛЮЧЕНИЕ

В результате выполненных исследований построены структурно-плотностные и геодинамические модели литосферы по трем профилям, пересекающим зону сочленения плиты Филиппинского моря и Азиатского континента в районе острова Тайвань. Основные результаты исследований сводятся к следующему:

1. Установлена структурно-плотностная дифференциация литосферы. По различию средних плотностей и особенностям латерального распределения плотностных неоднородностей в коре и мантии выделены три типа литосферных блоков: океанический, островодужный и окраинно-континентальный.

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

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

4. Показано, что выделенные первоначально по признакам структурно-плотностной дифференциации три литосферных блока выделяются также и по различию геодинамического режима. Численные значения напряжений, возникающих в литосфере под воздействием ее плотностной дифференциации, достаточны для формирования разломных или тектонически ослабленных зон (109 - Ю10 Па). На этой основе сделано заключение, что, коромантийные плотностные неоднородности являются одним из источников геодинамических напряжений, деформаций и перемещения масс в указанной среде.

5. На всех профилях зафиксирован режим нисходящего движения мантии Западно-Филиппинского блока; напротив, окраинно-континентальная литосфера испытывает подъем. Островодужные системы Рюкю и Лусон-Тайвань различаются высокой концентрацией касательных напряжений и деформаций на стыках разнородных морфоструктур. Основной причиной этого, на наш взгляд, является пограничное расположение островодужных систем между литосферными плитами с разнонаправленным движением слагающих их масс. Основной объем мантийного вещества под островодужными системами во всех случаях находится в нейтральном состоянии, без значимых перемещений.

6. Конфигурация глубинных структурно-плотностных и «геодинамических» границ литосферных блоков свидетельствует о некотором поддвигании мантийного субстрата Западно-Филиппинского блока под островные дуги Рюкю и Лусон-Тайвань.

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

133

Библиография Диссертация по наукам о земле, кандидата геолого-минералогических наук, Гильманова, Гульшат Забировна, Владивосток

1. Абрамов В.А. Структура тектоносферы Алданского щита и его мезозойских рудно-магматических узлов / (автореферат докторской диссертации) Иркутск: 1992. 35 с.

2. Абрамов В.А. Глубинное строение Центрально Алданского района. Владивосток: Дальнаука, 1995. 179 с.

3. Алакшин A.M. Строение земной коры западной части региона БАМ по гравиметрическим данным / ( автореф. канд. дисс.) Иркутск: 1986. 19с.

4. Алексидзе М.А., Картвелишвили К.М. Плотностные модели Земли. Нормальное гравитационное поле, обусловленное моделью Земли концентрических элипсоидальных слоев // Гравитационная модель коры и верхней мантии Земли. Киев: Наук. Думка, 1979. С. 15-20.

5. Артемьев М.Е. Структура гравитационного поля // Гравитационная модель коры и верхней мантии Земли.- Киев: Наук. Думка, 1979. С. 20-28.

6. Артюшков Е.В. Геодинамика//М.: Наука, 1979. 327с.

7. Балакина JI.M. Механизмы очагов промежуточных землетрясений Курило-Камчатской фокальной зоны // Изв. АН СССР. Физика Земли. 1981. № 8 С. 3-25.

8. Болдырев С.А., Ефимов А. В. Сейсмотектонический эффект мантийных неоднородностей активных окраин океана // Докл. АН СССР. 1985. Т. 285 N4

9. Брянский Л.И. Плотностаая структура земной коры и верхов мантии Восточной окраины Азиатского континента. Владивосток: Дальнаука, 1995. 142 с

10. Брянский Л.И. Стандартные гравитационные модели земной коры (на примере юга Дальнего Востока). Владивосток: ДВО АН СССР, 1988. 144 с.

11. Васильев Б.И., Чой Д.Р. Геология глубоководных желобов и островных дуг Тихого океана // Моногр. Владивосток. Дальнаука. 2001. 182с.

12. Ващилов Ю.Я. Глубинная структура Северо-Востока СССР по гравиметрическим данным // Гравитационная модель коры и верхней мантии Земли. Киев: Наукова Думка. 1979. с. 209-214.

13. Гайнанов А.Г. Аномалии силы тяжести и плотностные модели литосферы Тихого океана// Тихоокеанская геология, 1982, № 5 . С. 21-30.

14. Гатинский Ю. Г. Кайнозой юго-востока Азиатского континента и некоторые вопросы процесса рифтогенеза//Изв. Вузов. Геология и разведка 1980 N3 С. 19-31 (ст. 1); N 7. С. 28-36 (ст.2).

15. Гатинский Ю.Г. Связь геодинамики континентальных окраин ЮВ Азии с развитием океанов . Геология Мирового океана:. Докл. 27 МГК; Наука, 1984. Т. 6., ч. 1.

16. Геологический словарь т.1 // М.: Недра, 1973. 486с.

17. Геология дна Филиппинского моря // М: Наука, 1980.

18. Геотраверс Китай Филиппинское море - Тихий океан // Тихоокеанская геология. 1991. №4. С. 3-18.

19. Гильманова Г.З., Маслов JI.A., Кулинич Р.Г., Комова О.С. Геодинамика земной коры зоны коллизии Тайваня // Тихоокеанская геология. Том 21,N 3,2002 С.40-49.

20. Гильманова Г.З., Подгорный В.Я. Плотностная модель литосферы зоны перехода от Филиппинской плиты к Азиатской, на примере островной дуги Рюкю // Сборник докладов на IV Косыгинских чтениях. Хабаровск 2003 С. 238-258.

21. Гнибиденко Г. С. Тектоника дна окраинных морей Дальнего Востока. М: Наука, 1979. 164 с.

22. Евланов Ю.Б. Филиппинское море. Геологическое строение. Эволюция магматизма и осадконакопления. Владивосток. Дальнаука. 2000. 91с.

23. Исаев В.И. Принципиальная схема комплексной интерпретации гравиметрических съемок на примере профиля Средние Лангары Мухто (Северный Сахалин) // Тихоокеанская геология, 1983, № 1 , С. 107-111.

24. Кадастр зарубежных стран, обладающих природными ресурсами нефти и газа // Т.2 . Ленинград "Недра" Л.О.1983 с. 260-264 (364с.)

25. Кириллова Г.Л. Глубоководные желоба Филиппинского моря: строение и эволюция // Тихоокеанская геология. N 4, С. 19-31.1991.

26. Кириллова ГЛ. Кайнозойские осадочные бассейны Восточно-Китайского моря: строение и эволюция.// Тихоокеанская геология. N 5, С. 64-83.1991.

27. Комова О.С. Моделирование глубинных геодинамических процессов в зоне перехода от континента к океану с помощью численных методов // Тихоокеанская геология. N 3 1984. С. 110-115.

28. Комова О.С., Маслов JI. А., Туезов И. К. Решение прямой трехмерной задачи геомеханики для аномальных плотностных тел // Тихоокеанская геология. 1984. № 3. С. 110-115.

29. Константиновская Е.А. Тектоника восточных окраин Азии: Структурное развитие и геодинамическое моделирование // Автореферат на соискание ученой степени доктора геолого-минералогических наук. М.: ГИН РАН. 2002. 46с.

30. Косыгин В.Ю., Исаев В.И. Плотностная модель тектоносферы вдоль геотраверса Япония Гавайи // Тихоокеанская геология, 1986, № 5.

31. Кочергин Е.В., Павлов Ю.А., Сергеев К.Ф. Геомагнитные аномалии Курильской и Рюкю островных систем . "Наука" 1980. 115с.

32. Красовский С.С. Отражение динамики земной коры континентального типа в гравитационном поле. Киев, Наук. Думка, 1981.

33. Красовский С.С. Гравитационное моделирование глубинных структур земной коры и изостазии. Киев, Наук. Думка, 1989. 248с.

34. Кулинич Р.Г. Роль сквозных тектонических линеаментов в строении и эволюции Филиппинской плиты // Геодинамические исследования. N11. Тектоника Восточно-Азиатских окраинных морей. Москва 1988. С. 95-101.

35. Кулинич Р.Г., Маслов Л.А., Гильманова Г.З., Комова О.С. Плотностная модель и напряжения в земной коре северной части Японского моря // Тихоокеанская геология. Том 17, №2, 1998. С. 108-114.

36. Лобковекий Л.И. Геодинамика зон спрединга, субдукции и двухъярусная тектоника плит. М.: Наука, 1988. 252 с.

37. Маслов Л.А. Простая модель расчетов напряжений в земной коре и литосфере и результаты ее интерпретации // Тихоокеанская геология, 1997, том 16, №2, С. 38-45.

38. Маслов Л.А. Геодинамика литосферы Тихоокеанского подвижного пояса // Хабаровск-Владивосток, 1996 С. 199.

39. Маслов Л.А., Косыгин В.Ю., Комова О.С. Изучение глубинных геодинамических процессов в Курило-Камчатской зоне перехода на основе ее плотностных моделей // Тихоокеанская геология.- 1988. N 2 С. 3-10.

40. Маслов Л.А., Строев П.А., Комова О.С. Геодинамика Япономорской переходной зоны // Гравиметрические исследования на море. М.: Наука, 1988. С. 92-97.

41. Маслов Л.А., Туезов И.К. Расчет механических напряжений в литосфере востока Азиатского континента и зоны перехода на основе гравитационных данных // Тихоокеанская геология.-1982. N6. С. 41-47.

42. Муравьев А.В., Смирнов Я.Б., Сугробов В.М. Тепловой поток вдоль международного геотраверса через Филиппинское море по 18° с.ш. // Докл. АН СССР. 1988 Т. 299,N1. С. 189-193.

43. Новацкий В. Теория упругости. М: Мир 1975.

44. Петрищевский А.М. Статистические гравитационные модели литосферы Дальнего Востока// Владивосток: ДВ ун-т, 1988. 168 с.

45. Подгорный В.Я. Послойное плотностное моделирование литосферы (на примере юга российского Дальнего Востока и Северо-Востока Китая) // Автореферат кандидатской диссертации. Хабаровск, 1999. 28 с.

46. Ранкорн С.К. Динамические процессы в нижней мантии // Верхняя мантия. М.: Мир, 1975. С. 286-298.

47. Строев П.А., Гайнанов А,Г. 1984. Глубинное строение литосферы Япономорской переходной зоны по гравиметрическим данным // Морские гравиметрические исследования,М„ МГК АН СССР, С. 102-111.

48. Сурков B.C., Мореин П.И. Определение мощности земной коры на основе гипотезы изостазии//Геология и геофизика. 1971,№ 10. С. 63-71.

49. Туезов И.К. Геоэлектрический разрез литосферы и астеносферы Северо-Восточной Азии и прилегающих частей Тихого океана// Владивосток: Дальнаука, 1994. 300с.

50. Хэйес Д.Е. Окраинные моря юго-восточной Азии: их геофизические характеристики и структура // История и происхождение окраинных и внутренних морей / Докл. На 27-м Междун. Геол. Конгр. М.Наука. 1984. Т.6.С.30-44.

51. Щека С.А., Кулинич Р. Г., Высоцкий С.В. и др. Новые данные по геологии разломов Яп и Центральный В Филиппинском море // Докл. АН СССР. 1986. Т.286.№ 2. С.417-421.

52. Щукин Ю.К., Солодилов Л,Н. В.В. Федынский и изучение глубинного строения земной коры и верхней мантии // Геофизика Спец. Вып., 1998. С. 20-30.

53. Шараськин А.Я. строение и тектоно-магматическая эволюция дна Филиппинского моря // 27-й Международный геологический конгресс. Доклады. Том 6. часть II. "Наука". Москва 1984.С. 44-58.

54. Aiba J., Sekiya Е. Distribution and characteristics of the Neogene Sedimetary Basins aroundthe Nansei-Shoto (Ryukyu Islands) // J. Jpn. Assoc. Pet. Technol. N 44(5). 1979. pp. 97-108.

55. Ben-Avraham Z., Bowin C. An extinctspreading center in the Philippine sea // Nature. Vol. 240. 1972. pp.453-455.

56. Ben-Avraham Z. The evolution of marginal basins and adjacent shelves in East and Southeast Asia // Tectonophysics. Vol. 45. 1978. pp. 269-288.

57. Ben-Avraham Z., Emery К. O. Structural framework of the Sunda Shef // Bull. Amer. Assoc. Petrol. Geol. Vol.57, N 20. 1973. pp. 2323-2366.

58. Ben-Avraham Z., Uyeda S. The evolution of the China basin and the Mesozoic paleogeography of Borneo // Earth and Planet Sci.Lett. Vol. 18. 1973. pp. 365-376.

59. Big, C., Collision. Taiwan-style // Met. Geol. Soc. China, 4,1981. pp.91-102.

60. Bonafede M., Dragoni M., Quareni F. Displacement and stress field produced by a center of dilanion and by a pressure source in visco-elastic half-space // Jeophys. J. R. Astron. Soc. 1986. V. 87. P. 455-486.

61. Bowin C., Lu R.S., Lee C.S., Shouten H. Plate convergence and accretion in Taiwan -Lusonregion// Bull. Amer. Assoc. Petrol. Geol. 1978. Vol. 62, N 13. pp. 1645-1672.

62. Chang C.-P., J.Angelier, T.-Q. Lee, C.-Y. Huang. From continental margin extension to collision orogen: structural development and tectonic rotation of Henghun peninsula, southern Taiwan//Tectonophysics. 361 . 2002. pp. 61-82.

63. Coblentz, Richardson, On the gravitational potential of the Earth's lithosphere// Tectonics, vol. 13, N4,1994. pp. 929-945.

64. De Bremaecker J.C. Is the oceanic lithosphere elastic or viscous? // J.Geophys. Res. 1977. V. 82. N. 14. P. 2001-2004.

65. Elsasser W.M. Sea-floor spreading as thermal convection // J.Geophys. Res. V. 76. 1971. pp. 427-428.

66. Fabbri O., Fournier M. Extension in the southern Ryukyu arc (Japan): Link with oblique subduction and back arc rifting // Tectonics. V. 18. N 3. 1999. pp. 486-497.

67. Faure M., Monie P., Fabbri O. Microtectonics and Ar/Ar datings of the high pressure metamorphic rocks in the south Ryukyu Ark and their bearings on the pre-Eocene geodynamic evolution of eastern Asia // Tectonophysics, 156. 1988. pp. 133-143.

68. Font, Y., Liu, C. S., Schurle, P., and Lallemand, S., Constraints on backstop geometry of the southwest Ryukyu subduction based on reflection seismic data // Tectonophysics 333,2001. pp.135-158.

69. Geology of the Northern Philippine Sea. Geological Results of the GDP Cruises of Japan // Ed. T. Shiki. Tokay Univ. Press. Japan. 1985. 288 p.

70. Gilmanova G., Maslov L., Kulinich R Density model and crustal stresses in the Collision zone of Taiwan // New concepts in global tectonics. Proceedings book. 2002 pp. 213-224.

71. Hamilton W. Tectonics of the Indonesian region // Geol. Surve. Prof. Pap. Usgpo, Wash., 1979. 245 p.

72. Hayes D. Margin oceans of South-West Asia: Their geophysical Characteristics and structure // History and origin of margin and midland seas. Moscow: Nauka, pp. 30-44. 1984.

73. Hayes D.E. and Lewis, S. D., A geophysical study of the Manila Trench, Luzon, Philippines, Crustal structure, gravity, and regional tectonic evolution // J. Geophys. Res., 89 (Bll): 1984. pp. 9171-9195.

74. Herman В. M., Anderson R., Truchan N. Extensional tectonics in the Okinawa Trough // Mem. Am. Assoc. Pet. Geol. Bull. N 29. 1978. pp. 199-208.

75. Hess H.H. Major structural features of the western North Pacific, an interpretation of H.O. 5485, bathymetric chart, Korea to New Guinea // Bull. Geol. Soc.Amer., vol.59, 1948. pp.417-466.

76. Hilde T. Isezaki N., Wageman J. Mezozoic sea floor spreading in the North Pacific// Amer. Geophys. Union. Monogr. 1976.Vol. 19.pp. 205-226.

77. Hilde T. and Lee C. 1984, Origin and evolution of the west Philippine Basin: a new interpretation//Tetonophysics, 102, pp. 85-104.

78. Ho, C.Y., A synthesis of the geologic evolution of Taiwan // Tectonophysics 125, 1986. pp. 1-16.

79. Holoway N.H. North Palawan Block. Philippines its Relation to Asian Mailand and Role in Evolution of South China Sea // Bull. Amer. Assoc. Petrol. Geol. 1982. Vol. 66 N.9.

80. Honkura Y., Isezaki N., Yashawa K. Electrical conductivity structure beneath the Northwestern Philippine sea as inferred from the island effect on Minami-Daito island // J. Gepmagn. And Geoelec. 1982. Vol. 33, N 6. pp. 365-377.

81. Hotta H. A crustal section across Izu-Ogasawara arc and trench // J. Phys. Earth, 1970, vol.18, pp. 123-141.

82. Huang C. Y., Shyu C.T., Lin S.B., LeeT.Q. & Sheu D. D. 1992. Marine geology in the arc-continent collision zone off southeastern Taiwan: Implications for late Neogene evolution of the Coastal Range // Marine Geology 107, pp. 183-212.

83. Huang C., Peter В., Lin C., Wang Т., Chang C., Geodynamic processes of Taiwan arc-continent collisionand comparision with analogs in timor, Papua New Guinea, Urals and Corsica //Tectonophysics 325. 2000. pp.1-21.

84. Huang Pei-Hua and Fu Rong Shan. The mantle convection pattern and force source mechanism of recent tectonic movement in China // Physics of the Earth and Planetary Interiors. 1982. V.28. pp. 261-268.

85. Hutchison C. S. Tectonic Evolution of Sudaland: A Phanerozoic Synthesis// Bull. Geol. Soc. Malaysia. 1973. N 6. pp. 61-86.

86. Karig D.E. Origin and development of marginal bassins in the western Pacific // J. Geophys. Res. 1971. Vol. 76. pp. 2543-2461.

87. Karig D Structural history of the Mariana islands arc system // Bull. Geol. Soc. Amer. 1971. Vol 82. pp. 323-344.

88. Karig D. Plate cjnvergence between the Philippines and the Ryukyu Islands // Marine geology. 1973. 14. pp. 153-168

89. Karig D. Basin genesis in the Philippine sea // Initial reports of the Deep Sea Drilling Progect. Wash. (D.C.): US Gov. print. Off., 1975. Vol. 31. pp.857-879.

90. Karp B.Ya., Karnaukh and Madvedev S.N. Crustal Structure of the Urdaneta Plateau and Its Tectonic Implication. // Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , pp.125 -136.

91. Karp Boris Ya, Ruslan Kulinich Chuen-Tien Shyu and Chengsung Wang. Some features of the arc-continent collision zone in the Ryukyu subduction system, Taiwan Junction area // The Island Arc. 6, 1997, pp. 303-315.

92. Katsumata M., Sykes L. Seismicity and tectonics о the western Pacific: Izu-Mariana-Caroline and Ryukyu-Taiwan regions // J.Geophys. Res., 1969, V. 74, pp. 5923-5948.

93. Kimura M, Tectonic evolution of the East-China sea and the Okinawa trough // Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , pp.145 -153.

94. Kimura M., Honza E., Miymoto M. Submarine geology around the southern Ryukyu Islands, with special reference to the Okinawa Trough // Geol. Stud. Ryukyus. 1979. N 4. pp. 7993.

95. Kimura M„ Kaneoka I., Kato Y. et al. U.C. Report on DELP 1984 Cruises in the middle Okinawa Trough // Bui. Earthquake Research Institute.- 1986. N 61, pp. 269-310.

96. Kimura M., Ueda S., Kato Y. et al. Aktive hydrothermal mounds in the Okinawa Trough back-arc basin, Japan // Tectonophys. 1988 V. 145. pp. 319-324.

97. Kizaki K., Geology and tectonics of the Ryukyu Islands // Tectonophysics,125, 1986, pp. 193-207.

98. Kizaki K., Tectonics of the Ryukyu Island Arc // J. Phys. Earth, 1978, 26 suppl., pp. 301307.

99. Kohlstedt D. L., Evans В., Mackwell S.J. Strength of the lithosphere: Constrains imposed by laboratory experiments // Journal of Geopthys. Res. 1995. V.100. NB9. pp. 17587 17602.

100. Konishi K. Geotectonic framework of the Ryukyu islands ( Nansei-Shoto) // Geol. Mag. 1965. V.71 pp.437-458.

101. Kuramoto S., and Konishi, The southwest Ryukyu Arc is a migrating microplate ( forearc sliver)// Tectonophysics, 1989, V.163 , pp. 75-91.

102. Magnetic anomaly map of East Asia 1:4 000 000 CD-ROM Version by Geological Survey of Japan and Coordinating Committee for Coastal and offshore Geoscience Programmes in East and Southeast Asia (ССОР). Published by Geological Survey of Japan. 1996.

103. Malyarenko A.N. and Lelikov E.P. Granites and Associated Rocks in the Philippine Sea and The East China Sea // Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , pp.311-328.

104. Matsumoto Т., Tomoda Y. Numerical simulationof the mutual interaction between a trench and a seamount // J. Phys. Earth.- 1983. V. 31. pp. 281-297.

105. Melnichenko Y., I.Popova, V. S'edin and S. Tochilina. Geomorphologic characteristics of the Kyushu-Palau Ridge, Philippine Sea // Z. Geomorphology N.F. Suppl.-Bd. Berlin . Stuttgart. November 1999. pp. 183-192 .

106. Mrozowski C.L., Hayes D.E. The evolution of the Parese Vela basin, Eastern Philippine sea //Earth and Planet. Sci. Lett. 1979. Vol.46.pp. 49-67.

107. Mrozowski C.L., Hayes D.E., Taylor B. Multichannel seismic reflection surveys of Leg 60 sites, DSDP // Initial reports of the Deep sea Drilling Progect. Wash. ( D.C.): US Gov. print. Off., 1982a. Vol. 59. pp. 57-69.

108. Mrozowski C.L., Lewis S.D., Hayes D.E. Complexities in the tectonic evolution of the west Philippine basin // Tectonophysics. 1982 b. Vol. 82. pp. 1-24.

109. Murauchi S., Den N., Asano S.et al. Crustal structure of the Philippine sea.- Earth and Planet // Sci. Lett., 1979, vol 46, pp. 49-67.

110. Murphy R.W. The Manila trench west Taiwan fold belt: a flipped subduction zone // Bull. Geol. Soc. Malasia, 1973, N 6.

111. Nakamura M., Yoshida Y., Zhao D,, Katao H., Nishimura S. Three-dimensional P-and S-wave velocity structures beneath the Ryukyu arc // Tectonophysics 369 (2003) pp. 121-143.

112. Nafe J. E. & Drake С. C. 1963. Physical properties of marine sediments // In Hill M. N. Ed . The Sea, vol. 3, P. 794-815. Wiley Interscience, New York.

113. Nishimura K. A schematic model concerning the formation of the Benioff zones as inferred from examination of global-scale geoid undulation // Дзисин "Землетрясения". 1982 .V. 35.

114. Obzhirov A.I. Gas Distribution in the Philippine Sea. // Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , pp.379-383.

115. Oldenburg, D. W., Conductivity structure of oceanic upper mantle beneath the Pacific plate, Geophys.J.Roy.astr. Soc., 65,1981, pp.359-394.

116. Oshima S., Takanashi M. Geological and geophysical survey in the Okinawa trough and the Nansei-shoto areas // Bull. Hydrogr. Dep. 1988. Vol.24, pp. 26-27.

117. Richard H. Rapp A Global atlas of surface heights based on the adjusted seas at altimeter data // The Ohio State University Department of Geodetic Science and Surveying 1982.

118. Runcorn S.K. Flow in the mantle inferred from the low-degree harmonics of the geopotential // Geophys. J. Roy. Astron. Soc.-1967. V.14.- N.3. pp. 375-385.

119. SaTdwell D.T., Smith W.H.F. Marine gravity anomaly from Geosat and ERS-1 satellite altimetry // J. Geophys. Res. 1997. V.102, pp. 10039-10054

120. Scholl D. W., Von Huene R., Vallier T. L., Howell D. G. Sedimentary masses and concepts about tectonic processes at under thrust ocean margins // Geology. 1980 V 8, N 12. pp. 564-568.

121. Sclater J.G., Karig D.E., Lawer L.A., Louden K. Heat floor, depths and crustal thickness of the marginal basins of the South Phylippine sea // J.Geophys. Res. 1976a. Vol 81. pp. 309318.

122. Sclater J.G., Karig D.E., Lawer L.A., Louden K. Heat flor, depths and crustal thickness of the marginal basins of the South Phylippine sea // J.Geophys. Res. 1976b. Vol 81. pp. 309-318.

123. Seno Т., Kurita K. Focal mechanisms and tectonics in the Taiwan-Philippine region // J.Phys. Earth, 26.pp.249-263.

124. Seno Т., Stein, and Gripp, A.E.,. A model for the motion of the Philippine Sea plate with NUVEL-1 and geological data//J. Geophys. Res., 1993, 98, pp. 17941-17948.

125. Seno Т., The instantaneous rotation vector of the Philippine Sea plate relative to the Eurasian plate // Tectonophysics, 1977,42, pp.209-226.

126. Seno T. and Maruyama S., 1984, Paleogeographical reconstruction and origin of the Philippine Sea // Tetonophysics, 102, pp.53-84.

127. Shu-Kun Hsu. Subduction collision Complexities in the Taiwan-Ryukyu Junction Area: Tectonics of the Northwestern Corner in the Philippine Sea Plate // TAO, Supplementaty Issue, May 2001, pp. 209-230.

128. Sibuet J.-C , Letouzey J., Barbier F. et al. Back-arc extension in the Okinawa Trough // J. Geophys. Res. 1987 V. 92, N B13. pp. 14041-14063.

129. Svininnikov A.I. Physical Properties of Rocks and Sediments from the Philippine Sea// Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , pp.137-144.

130. Taylor В., Hayes D.D. The tectonic evolution of the South China Basin // The tectonic and Islands. Edited by Hayes D.E. Geophys. Monogr. Amer. Geophys. Union. Wash. D.C. 1980. Vol. 23. pp. 89-104

131. Taylor В., D.E. Hayes . Origin and history of the South China Basin, The Tectonic and Geologic Evolution of South-east Asian Seas and Islands // Part 2, Geophys. Monogr. Ser., 1983. vol 27,edited by D.E.Hayas, pp. 23-56, AGU, Washington, D.C.

132. Toksoz M.N. Hsui A.T. Numerical studies of back-arc convection and the formation of marginal basins // Tectonophysics. 1978. V. 50. N. 2/3. pp. 177-196.

133. Tokuyama H. et al., 1986, Marine geology and subcrustal structure of the Shikoku Basin and the Daito Ridges Region in the northern Philippine Sea // Bull. Ocean Res. Inst., Univ. Tokio, 14, pp.1-158.

134. Tokuyama H. Origin and development of the Philippine sea // Geology and geophysics of the Philippine sea 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405, pp. 155-163.

135. Utada H., Y. Hamano, Segawa J. Conductivity anomaly around the Japanese islands // Geology and geophysics of the Japan sea. Japan-USSR Monograph Series, Vol.1, pp. 103-149. 1996.

136. Uyeda S. Formation of back- arc basins// DELP News. 1986. N 12. pp. 8-12. Uyeda S., Ben-Avraham Z. Origin and development of the Philippine sea // Nature. 1972. Vol. 240.pp. 176-178.

137. Vanpe J.M. Thermo-mechanical convection in a subduction zone and initiation of back-arc spreading // Annals Geophysical. 1984. V. 2. N. 3. pp. 343-352.

138. Wageman J. M., HildeT.W., Emery K. Struktural framework of the East China Sea And Yelow Sea // Bull. Am. Assoc. Pet. Geol. 1970. V. 54(9). N 4 pp. 1-10.

139. Wang С., S.W. Chuang, M.L. Li, W.B. Cheng. Lithospheric Structure of the Philippine Sea plate near the western end of the Ryukyu subduction zone and some of its tectonic effects // TAO, Supplementary Issue, 287-304, May 2001. pp. 287-304.

140. Wang 1CT., ICMcIntosh, Y.Nakamura, C.Liu and H.W.Chen. Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data // Marine Geophysical Researches 2001. 22: pp. 265-287.

141. Watanabe Т., Langseth M.G., Anderson R.N. Heat floor in back-arc basins of the Western Pacific // Island arcs, deep sea trenches and back-arc basins wash. (D.C.): Amer. Geophys. Union, 1977. Vol.1 pp. 137-161.

142. Watts, A.B., J.K. Weissel and R.L. Larson. Sea-floor spreading in the marginal basins of the western Pacific // Tectonophysics 1977., 37, pp.167-182.

143. Wu F. T. 1970. Focal mechanisms and tectonics in the vicinity of Taiwan // Bulletin of Seismological Society of America 60, pp. 2045-56.

144. Wu F. T. 1978 Recent tectonics of Taiwan // Journal of Physics of the Earth 26 (Suppl.), pp.256-99.

145. Xi Jiebo, Shen Yi, Wu Jiansheng, Wang Jialing. On Evolution of the NNE and NW Trending Faults in the East China Sea // J. Tongji Univ. 1988. V. 16, № pp. 1-10.

146. Yamano M., Kinishita M. Heat flor in the Philippine sea. // Geology and geophysics of the Philippine sea. 1995. Japan-Russia-China Monograph (JRCM). TERRAPUB, Tokio . Vol. 405 , P.59-75.

147. Yu , S. В., H.Y. Chen, and L. C. Kuo, Velocity field of GPS stations in the Taiwan area // Tectonophysics. 1997.274, pp.41-59.