Скоростное строение Юго-Восточной Европы, Малой Азии и Восточного Средиземноморья по сейсмологическим данным тема автореферата и диссертации по физике, 01.04.12 ВАК РФ

1. Введение

2. Глава I. Строение Юго-Восточной Европы, Малой Азии и

Восточного Средиземноморья по данным геофизических исследований

§ 1.1 Геологическое строение

§ 1.2 Гравитационное и магнитное поля

§ 1.3 Тепловой поток

§ I.4 Сейсмические и сейсмологические исследования . 24 Выводы к главе I

3. Глава 2. Методы решения обратной кинематической задачи сейсмологии в двумерном и трехмерном случаях

§ 2.1 Методы решения обратной кинематической задачи для двумерно и трехмерно неоднородных сред, использующих параметризацию среды

§ 2.2 Метод Бэйкуса-Гильберта для решения одномерной обратной кинематической задачи

§ 2.3 Обощение метода Бэйкуса-Гильберта на двумерный случай (обратная задача для поверхностных волн)

§ 2.4 Обобщение метода Бэйкуса-Гильберта на трехмерный случай (обратная задача для объемных волн)

Выводы к главе

4. Глава 3. Распределение групповых скоростей поверхностных волн Релея и Лява в Юго-Восточной Европе и Малой Азии

§ ЗЛ Исходный материал наблюдений и его предварительный анализ

§ 3.2 Выбор ориентации координатной системы

§ 3.3 Распределение групповых скоростей волн Релея

§ 3.4 Распределение ррупповых скоростей вон Лява

Выводы к главе

5. Глава 4. Распределение скоростей объемных Р-волн в верхней мантии Юго-Восточной Европы, Малой Азии и

Восточного Средиземноморья

§ 4.1 Использованный материал сейсмологических наблюдений

§ 4.2 Особенности отклонений времен пробега Р волн в Юго-Восточной Европе, Малой Азии и Восточном

Средиземноморье

§ 4.3 Выбор координатной системы

§ 4.4 Скоростное строение района I

§ 4.5 Скоростное строение района П . П

§ 4.6 Скоростное строение района Ш

§ 4.7 Скоростное строение района Юго-Восточной Европы, Малой Азии, Восточного Средиземноморья и Закавказья

6. Глава 5. Связь поученных распределений скорости с особеностями строения района, полученных по другим геофизическим данным

§ 5.1 Сравнение распределений групповой скорости с особенностями строения земной коры и верхней мантии

§ 5.2 Сравнение распределений скорости объемных Р волн с аномалиями геофизических полей

Выводы к главе

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

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

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

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

ЦЕЛЬ РАБОТЫ. Целью работы являлось исследование горизонтальных неоднородностей в распределении упругих характеристик коры и верхней мантии в районе Юго-Восточной Европы, Восточного Средиземноморья и Малой Азии, которое включало в себя решение следующих задач:

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

2. Построение трехмерного распределения скоростей Р волн в верхней мантии по данным о невязках времен пробега этих волн.

3. Сравнительный анализ полученных скоростных распределений и результатов других геофизических исследований в данном районе (гравиметрических исследований, измерений теплового потока, данных ГСЗ).

ИСХОДШЙ МАТЕРИМ. Данные о групповых скоростях волн Релея и Лява в интервале периодов 10-30 с были получены автором диссертации путем обработки сейсмограмм станций Афины, Белград, Будапешт, Кишинев, Носов, Львов, Симферополь, Сопроы, Ужгород, Триест. Данные о невязках времен Р волн были взяты из Международного Сейсмологического Бюллетеня [79] за период с 1970 по 1981 г. включительно и подвергнуты статистической обработке.

НАУЧНАЯ НОВИЗНА. Впервые получены распределения групповых скоростей поверхностных волн для периодов 10-30 с в районе Юго-Восточной Европы и Малой Азии, на основании которых делаются выводы о горизонтальных неоднородностях строения земной коры в этом районе. Построены поправки к скоростному разрезу Джеффриса для Р волн в верхней мантии до глубины 300 км. Рас пределение скоростей в верхней мантии обнаруживает корреляцию с распределением аномалий силы тяжести и тепловым потоком.

ПРАКТИЧЕСКАЯ ЗНАЧИМОСТЬ. Выявление аномалии скоростей поверхностных волн и объемных Р волн могут быть положены в основу геодинамических концепций и гипотез об эволюции земной коры и верхней мантии в районе Юго-Восточной Европы и Малой Азии. Полученное распределение скоростей Р волн может быть использовано для построения регионального годографа, что важно для уточнения координат гипоцентров землетрвсений в высокосейсмичных областях Румынии, Греции, Турции, Болгарии и др.

АПРОБАЩЯ РАБОТЫ. Основные положения диссертационной работы докладывались на заседании Рабочей группы МСССС по поверхностным волнам (Ленинград, 1963), на научных семинарах лаборатории сейсмологии кафедры физики Земли ЛГУ (1985 г.), на заседаниях секции сейсмологии Болгарской Академии Наук (1963, 1984 гг., София), на XIX Генеральной Ассамблее Европейской Сейсмологической Комиссии (Москва, 1984).

ПУЕЛИКАЩИ. По теме работы опубликована статья в Геофизическом журнале Болгарской АН, 1984.

СТРУКТУРА РАБОТЫ И ОБЪЕМ. Диссертация состоит из введения, пяти глав, заключения, списка литературы, состоящего из 109 названий, и приложения. Работа содержит страниц машино

ВЫВОДЫ К ШВЕ 4

Показано, что имеют место

I. Аномалии времен пробега Р-волн в районе Юго-Восточной

Европы, Малой Азии и Восточного Средиземноморья. На некоторых трассах они достигают до -5 сек, +3-Г+4 сек. Эти данные получены на основе статистического анализа исходного материала, что позволяет считать их надежными.

2) Решение обратной трехмерной задачи позволило выявить основные; скоростные неоднородности в районе, которые сводятся к следующему: а) пониженными скоростями распространения объемных Р-волн характеризуются области Паннонского бассейна, Эгейского моря и Закавказья. Поправки к скоростному разрезу Джеффриса достигают в этих областях соответственно -0.13,-0.14,-0.25 км/сек. б) повышенными скоростями характеризуются Восточное Средиземноморье, южная часть Малой Азии - Тавриды, Черное море и

Восточные Карпаты, Динариды и Пелагониды. Поправки достигают максимальных значений на глубине 100 км и соответственно равны +0,25, +0.18, +0.18 км/сек. в) зоны близких к нулю поправок к скорости расположены над центральной частью Балканского полуострова и над Северо-Анато-лийским разломом. Они соединяют между собой области пониженных значений скорости. г) на глубинах 200—250 км отмечается увеличение областей, в которых наблюдаются отрицательные поправки к скорости. Они охватывают почти целиком Юго-Восточную Европу, кроме Карпат и Пе-лагонид. Это явление, вероятно, связано с существованием зоны пониженных скоростей на этих глубинах. Такимх особенностей также не наблюдается под Восточным Средиземноморьем и Малой Азией.

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

4) Предложенный в [201] метод переориентации координатной системы позволяет достаточно точно определить оптимальную ориентацию, но все же основным критерием остается анализ результатов решения обратной задачи.

ГЛАВА 5

СВЯЗЬ ПОЛУЧЕННЫХ РАСПРЕДЕЛЕНИЙ СКОРОСТИ С ОСОБЕННОСТЕЙ

СТРОЕНИЯ РАЙОНА, ШШЧЕННЫМИ ПО ДРУГИМ ГЕОФИЗИЧЕСКИМ

ДАННЫМ

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

§ 5.1. Сравнение распределения групповых скоростей поверхностных волн с особенностями строения земной коры и верхней мантии

Распределение групповых скоростей поверхностных волн с периодами 10-30 сек отражает горизонтальные неоднородности в земной коре и частично верхней мантии. Волны с такими периодами формируются полностью в верхних 40-50 км Земли. Поэтому полученные распределения (гл.З, рис.9-18) можно сравнивать с результатами гравиметрических и магнитных измерений и с данными ГСЗ.

Для периодов 10-15 км на распределение групповых скоростей поверхностных волн будет очень сильно влиять осадочный слой, его толщина и скорости в нем. На рис.47 представлена глубина

0 МЮНХЕН

Рис. 47 Схема глубины залегания консолидированного фундамента 1493

-лцгзалегания консолидированного фундамента в Юго-Восточной Европе [49]. Данные о толщине осадочного слоя в Восточном Средиземноморье, Греции и Эгейском море приведены в $ 1.4. Для Малой Азии отсутствуют данные о мощности осадочного слоя, поскольку объем сейсмических исследований в этом районе очень небольшой.

Зона наиболее выраженного минимума групповой скорости для обоих типов поверхностных волн, расположена над Черным морем. Там по результатам сейсмических исследований мощность осадочного слоя, составленного рыхлыми низкоскоростными породами достигает 14 км. Кроме этого, наличие жидкого слоя влияет на дисперсию волн Релея, что приводит к большей выраженности минимума в распределении этого типа волн. Второй минимум, приуроченный к Лредкарпатскому прогибу вызван большой толщиной осадочного слоя, которая по данным [49] достигает до 12 км Перипьенинского линеа-мента (Западные Карпаты) и увеличивается на восток. Область относительного понижения групповых скоростей в Эгейском море связана с небольшим по мощности слоем осадочных пород (до 3-4 км). Она лучше выделяется в случае волн Релея, что опять же вызвано наличием жидкого слоя. Максимумы в распределении групповых скоростей, приуроченные к Родопскому массиву и центральной части Балканского полуострова характеризуются либо полным отсутствием, либо наличием слаборазвитого осадочного слоя (до 2 км) [25,49]. Уменьшение групповых скоростей в северном и северо-восточном направлениях совпадает с увеличением толщины осадочных слоев к Карпатам, которая в Фокшанской депрессии достигает 18 км. Как было указано в § 3.4, максимум в распределении волн Лява смещен на северо-восток по отношению к максимум для волн Релея. Это вызвано большей глубинной проникания волн Лява. В Родопском массиве к поверхности выходит байкальский фундамент, сложенный из гранитных пород. На Мизийской платформе, консолидированная кора расположена на глубине 5-7 км, и образована более высокоскоростными породами (§ 1.4). Волны Лява с периодом 15 сек захватывают" и эти слои, что определяет перемещение максимума в распределении скоростей этого типа волн.

Как видно из рис.47 в области Внутренних Карпат слабо развит осадочный чехол. В некоторых отдельных зонах он достигает до 5 км толщины. В остальной же части к поверхности выходит кристаллический фундамент. Такое строение этой области предполагает существование максимума групповой скорости для обоих типов поверхностных волн. В полученных же распределениях скоростей такого не наблюдается. Эта область является переходом от высокоскоростной зоны в центре Балканского полуострова к низкоскоростной вдоль Предкарпатского прогиба. Простирание изолиний также не оконтуривает Прёдкарпатский прогиб в области Восточных Карпат. Такой факт может быть вызван как недостаточной информативностью исходного набора данных, так и особенностями строения.

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

Особенности строения, которые могут привести к наблюдаемому плавному изменению скоростей, - это существование высокоскоростных осадочных пород вдоль Восточных Карпат, а также более низкие скорости распространения волн в кристаллическом фундаменте внутренней зоны, Необходимо учитывать, что в пределах складчатой зоны Восточных Карпат по данным ГСЗ мезозойско-палеозой-ский осадочный комплекс характеризуется \пА- 4.5-5.5 км/сек. Скорости продольных волн в гранитном слое достигают 6.0-6.2км/с, таким образом, перепад между скоростными характеристиками осадочных слоев и кристаллического фундамента является не слишком резким. К сожалению, сейсмические исследования во внутренней зоне Карпат и Восточных Карпат в Румынии нам не известны кроме области Фокшанской депрессии, так что полученные распределения, вероятно, отражают особенности скоростного строения этих областей.

В § 3.3 было отмечено повышение групповых скоростей волн Релея на северо-западном крае области, охватывающее западную часть Паннонского бассейна. При сравнении с рис.4?, видно, что это отвечает различной толщине осадочных комплексов в восточной и западной части Венгрии, которые достигают соответственно 9км и 3-4 км.

В распределении скоростей волн Лява на этом периоде выделился небольшой минимум западнее Белграда. По [8] имеются исследования, которые свидетельствуют о низких скоростях на границе Мохоровичича в юго-восточном крае Паннонского бассейна. Учитывая небольшую мощность земной коры в этом районе, можно предположить, что волны Лява отражают именно эти особенности в строении.

Слабое понижение групповой скорости волн Лява над Динарида-ми и Адриатическим морем вызвано наличием мощных осадочных слоев (§ 1.4).

Второй максимум в распределении групповой скорости, приуроченный к Малой Азии по своим значениям подобен Центрально- Бал-кановому. Поскольку, сейсмических данных о строении этого района нет, а геологические исследования показали, что к поверхности в некоторых областях выходит кристаллический фундамент, ' логично предположить аналогичные скоростные характеристики массивов Анатолии и Таврид с Рило-Родопским массивом, т.е.

6.0-6.1 ад/сек. Сглаживание наблюдаемое на периодах 20, 25,30 сек. которое выражено для волн Лява, отражает особенности в строении района, а также свойства поверхностных волн. Волны с такими периодами охватывают всю земную кору и верхи мантии. Поэтому распределения, полученные на периодах 20-30 сек, следует сравнить с картой глубины залегания поверхности Мохоровичича [4Й на рис.48. Данные о Греции и Эгейском море, приведены в §1.4. На периоде 20 сек распределение для волн Релея сохраняет свой вид, В то же время, как было показано в § 3.3 и § 3.4, в области Черного моря происходит быстрое увеличение групповой скорости. Это вызвано небольшой толщиной земной коры и высокими скоростями в консолидированном фундаменте, который состоит только из базальтового слоя (§ 1.4). На периодах 25,30 сек была отмечена тенденция к изменению минимума на максимум в распределении, что определяется влиянием высоких скоростей в верхней мантии. Эта зако-~ номерность сильнее выражена в распределении волн Лява.

Продолжает выделяться область минимума в Предкарпатской прогибе. Как видно, эта зона характеризуется очень большой толщиной земной коры (до 50-60 км), что предопределяет существование минимума в распределении скоростей на всех периодах для обоих типов волн. Надо отметить, что с увеличением периода для волн Лява эта зона выявляется лучше, и она уже приблизительно оконту-ривает простирание Карпат (см.рис.17 гл.З). На Т = 30 сек в распределении волн Лява пониженные скорости охватывают всю область Карпат, что может быть вызвано неоднородностями в строении верхней мантии.

Область Пелагонид и Динарид, которая была выделена как зона понижения скоростей волн Лява на Т=15 сек и Релея на Т= 20, 25 сек, характеризуется большой толщиной земной коры - до 45 км.

Рис. 48 Схема глубины залегания поверхности Мохоровичича С493

Для периода Т-25 сек в распределении волн Релея (рис.12, гл.З) был выделен минимум в южной части полуострова Пеллопенес. Относительно пониженные значения в Эгейском море наблюдались на всех периодах, но при решении обратной задачи для выборки данных по групповым скоростям волн Релея на периоде 25 сек оконтурилась зона, выделенная и по данным ГСЗ. В § 1.4 было указано, что в южной части Эгейского моря граница Мохоровичича расположена на глубине - 22-25 км с Ч, = 7.7-7.9 км/сек. Это является причиной уменьшения значений групповых скоростей на этом периоде. В этом распределении вновь выявлена зона повышенных скоростей в западной части Паннонского бассейна, на этот раз ее можно связать с особенностями строения консолидированного фундамента, который по данным магнитометрии и бурения имеет разный возраст и сложен из различных пород. Отмеченное расширение зоны максимума в центральной части Балканского полуострова для волн Релея вызвано теми же причинами, которые привели к смещению максимума для волн Лява на периодах 10,15 сек. Для волн Лява на Т = 25, 30 сек сильнее выражена тенденция к сглаживанию скоростей. Это явление следует объяснить не большей гомогенностью района в глубине, а тем, что в областях, где наблюдаются мощные низкоскоростные осадочные толщи, земная кора имеет небольшую мощность и верхняя мантия характеризуется нормальными и даже повышенными скоростями. Наиболее высокими значениями групповых скоростей выделяется Малая Азия, что объясняется повышенными скоростями в земной коре и верхней мантии в этой зоне.

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

I) Сглаженные групповые скорости для периодов 10,15 сек

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

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

§ 5.2. Сопоставление распределений сглаженных поправок Р-волн с особенностями геофизических полей (район Юго-Восточной Европы, Малой Азии, Восточного Средиземноморья и Закавказья)

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

Как уже отмечалось (см.рис.3, гл.1) здесь выделены крупные положительные изостатические аномалии в Паннонском бассейне, Эгейском море и Закавказье. Они совпадают с определенными при решении обратной задачи зонами отрицательных поправок к разрезу Джеффриса. Положительные региональные изостатические аномалии требуют существования избытка масс на глубинах 200 км. В тоже время на этих глубинах у нас была определена тенденция к расползанию зоны пониженных скоростей в районе Юго-Восточной Европы и Малой Азии.

Центральная часть Балканского полуострова, определяемая положительными, но менее интенсивными региональными изостатически-ми аномалиями отвечает слегка заниженным скоростям на глубине 200-250 км. В районе Малой Азии по данным [94] наблюдаются положительные изостатические аномалии, но, к сожалению, нам неизвестно, влияние каких масс принималось при их определении, и соответственно, выражают ли они региональные отклонения от изоста

-4 5"0зии. Области Северо-Анатолийского разлома и срединных массивов Анатолии в распределении осредненных поправок также отрицательными значениями. Таким образом, можно прийти к выводу, что в рамках исследуемого района областям с отрицательными значениями осредненных поправок к скорости Р-волн соответствуют положительные региональные изостатические аномалии.

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

В тепловом поле (рис.4, гл.1) области Паннонского бассейна, Эгейского моря, западной Малой Азии, полосы северо-восточного -юго-западного простирания в Малой Азии и Закавказье характеризуются повышенными значениями теплового потока. Как уже было указано выше и в § 4.3-4.7 эти зоны характеризуются пониженными значениями скорости Р-волн. Восточное Средиземноморье, Черное море, Пелагониды, Динариды, в которых тепловой поток ниже среднего, являются зонами повышенных скоростей. Однако, отмеченная корреляция распределений поправок к скорости Р-волн с региональными изостатическими аномалиями и тепловым потоком содержит в себе некоторое противоречие. Избыток масс требует увеличения плотности вещества на этих глубинах. В тоже время, при прямой зависимости между скоростями Р-волн и плотностью ^ (закон Берча) \/р = а + , увеличение ^ должно приводить к увеличению скорости. Это несогласие было отмечено в работе СВОИ (см. § 1.4).

Его возможно объяснить следующим образом. В работе Берча Г9] приведены результаты лабораторных экспериментов по исследованию изменений скоростей продольных и поперечных волн при больших давлениях и температурах, аналогичных существующим в верхней мантии. С целью объяснить существование низкоскоростного слоя в верхней мантии было предложено связывать его с изменениями химического состава вещества. При предположении об оливиновом состоянии мантии образование зоны пониженных скоростей было объяснено увеличением содержание фаялита Fa . Этот минерал имеет большую плотность, чем форстерит Fo , скорости упругих волн в нем меньше. Исходный состав вещества мантии приблизительно равен FogoFa10 . При увеличении содержания фаялита скорости Р- и

S -волн падают. Например, при температурах около 1600°, что соответствует глубине 200 км, при содержании фаялита 10% - скорость Р-волн 8.0 км/сек, при 30% - 7.5 км/сек, при 50% -6.99 км/сек. Эти эксперименты позволили объяснить существование зоны пониженных скоростей наличием обогащенного железом слоя Сдо 23% фаялита). В соответствии с этой интерпретацией, зона пониженных скоростей характеризуется более высокой плотностью по сравнению с вышележащими и нижележащими слоями. Вероятная плото ность этого слоя около 3,48 г/см . Это объяснение дает совпадение зон пониженных скоростей для Р- и S -волн по глубине.

При рассмотрении распределений сглаженных поправок к скорости Р-волн необходимо учитывать, что они получены на основе осреднения, так что локальные значения отклонений от стандартного разреза могут быть значительно больше. Поэтому вероятно, что на глубинах 200-250 км в Юго-Восточной Европе расположен слой пониженных скоростей. Тогда положительные региональные изостатические аномалии возможно объяснить увеличением содержания железа на этой глубине. Необходимо учитывать также, что высокие значения теплового потока предполагают неглубокое расположение изотерм 1200°-1600°, вызванное поднятием разогретого мантийного вещества. В § Í.2 было указано, что для объяснения наблюдаемых аномалий Буге в Паннонском бассейне предлагалось увеличение плотности в верхней мантии. Этот факт ставился под сомнение в связи с установлением волновода на глубине 51-97 км. Предположение об изменении химического состава позволяет снять это противоречие. Следовательно, можно заключить, что зоны пониженных скоростей, выделенных в мантии в областях Ланнонского бассейна, Эгейского моря и Закавказья не связаны с разуплотнением вещества, как принято считать в настоящий момент, а наоборот, более вероятно объяснение, что они связаны с поднятием более плотного, богатого железом вещества к поверхности. Это позволяет связать установленные изостатические, тепловые и скоростные аномалии в исследуемом районе.

Зоны повышенных скоростей связаны с низкими значениями теплового потока. В Восточном Средиземноморье кроме этого наблюдаются и интенсивные отрицательные региональные изостатические аномалии. Одно из возможных объяснений существования таких аномалий - это внедрение в литосферу малоплотных пород и выступ к земной поверхности волновода [5]. Очевидно, эта гипотеза не находит подтверждения ни в особенностях распределения теплового потока, ни в сглаженных поправках к скорости Р-волн. Более вероятен другой механизм образования таких аномалий, предложенный в той же работе [5]. Нарушение изостатического равновесия (отрицательные региональные аномалии) связаны с оттоком глубинного вещества. На больших глубинах появление отрицательных аномальных масс связано с погружением границы раздела между астеносферой и нижележащими слоями мантии. Полученные распределения сглаженных поправок Р-волн подтверждают эту гипотезу. На глубинах до 300 км в области Восточного Средиземноморья и Малой Азии не был выделен слой, в котором бы наблюдались отрицательные поправки к скоростному разрезу Джеффриса. Наблюдаемые величины теплового потока

§ 1.3), показывают также, что в этих областях высокотемпературные слои находятся на больших глубинах. Тектонический механизм, приводящий к наблюдаемым отрицательным изостатическим аномалиям, связан с погружением литосферы и сопутствующим этому оттоком глубинных масс. Предположение о химическом изменении вещества в зоне пониженных скоростей объясняет и существование отрицательных региональных изостатических аномалий в Восточном Средиземноморье. Положительные поправки к скорости Р-волн указывают, что зона пониженных скоростей находится, по крайней мере, глубже 300 км. Поскольку вещество в этой зоне имеет большую плотность, чем вышележащие слои, то его отсутствие приведет к созданию отрицательных региональных аномалий. Это объясняет и низкие значения теплового потока. Если допустить существование астеносферы (предполагаем, что она совпадает с зоной пониженных скоростей), то толщина литосферы в районе Восточного Средиземноморья и Малой Азии больше 300 км, что соответствует толщине литосферы в пределах консолидированных платформ и щитов. Весьма вероятно, что такой слой не существует, или по крайней мере не проявляется в скоростном распределении Р-волн.

Изостатические региональные аномалии в Черном море, Пелаго-нидах и Динаридах близки к нулю, сглаженные поправки до глубины 200 км положительны, на 250-300 км они около 0 или отрицательны. Тепловой поток в этих районах ниже среднего (60 мВт/м^). Поэтому можно считать, что эти области изостатически скомпенсированы, скоростной разрез близок к разрезу Джеффриса. Пониженные значения теплового потока могут быть вызваны как глубинными аномалиями (погружением астеносферы на большие глубины), так и строением земной коры (низкой теплогенерацией пород, большой толщиной осадочного комплекса).

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

Все существующие тектонические построения к настоящему времени основаны на геофизических данных о строении земной коры и верхов мантии. ЧВ качестве критерия выделения отдельных плит используется распределение сейсмичности и механизмов очагов землетрясений. В тоже время, различные литосферные блоки должны отличаться и по своим физическим параметрам. Поэтому можно предположить, что выделенные литосферные блоки найдут отражение и в распределении поправок к скоростям Р-волн. На рис.49 представлено деление исследуемого района на основе анализа механизмов землетрясений, проведенное Маккензи Во]. При сравнении с полученными нами распределениями для района Юго-Восточной Европы, Малой Азии, Восточного Средиземноморья и Закавказья видно, что (рис.44-46, гл.4) только часть границ и блоков выражена в верхней мантии изменениями скоростных параметров среды. Границы между блоками будут определяться большими изменениями в скорости, следовательно, на глубинах, соответ

Йдс. 49 Схема границ плит и их движений Е ЗОН . Стрелками указаны движения плит по отношению к Евразии, а - границы, вдоль которых происходит расширение; б -трансформные разломы; в - границы, вдоль которых происходит сжатие. I - Евразиатская плита; 2 - Африканская плита; 3 - Иранская плита; 4 - Южно-Каспийская плита; 5 - Турецкая плита; б - Эгейская плита; 7 - Черноморская плита; 8 - Арабская плита

-Л5 6ствующих верхней мантии, выражена, например, граница между Черноморской и Турецкой плитой - Северо-Анатолийский разлом. Выделенная зона пониженных скоростей в Эгейском море и западной Малой Азии совпадает с границей Эгейской и Турецкой плит.Эллинский желоб - граница между Африканской и Эгейском плитой выражен в скоростном распределении, но только в своей западной части. Восточная часть этого желоба - от острова Родос к Кипру не нашла отражения в глубинном распределении, так как Малая Азия и Восточное Средиземноморье характеризуются приблизительно одинаковыми скоростями Р-волн. Субдукционные зоны в Закавказье выражены как области пониженных скоростей Р-волн. Черноморская плита характеризуется как приблизительно однородный скоростной блок в глубине, но ее западная граница не нашла отражения в распределении скоростей.

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

I) В пределах исследуемого района выявлена корреляция между ' наблюдаемыми региональными изостатическими аномалиями, тепловым потоком и распределением сглаженных поправок к скорости Р-волн. Отрицательные поправки к скорости наблюдаются в зонах интенсивных положительных изостатических аномалий и повышенного теплового потока. Это объясняется поднятием астеносферы к поверхности, которая характеризуется более высокими плотностями, чем вышележащие слои. Для реализации такого процесса необходимо наличие восходящих движений в мантии.

2) Отрицательные региональные изостатические аномалии связаны с областями повышенных значений скорости в Восточном Средиземноморье и пониженными значениями теплового потока. Это объясняется погружением зоны пониженных скоростей в этой районе, или ее полным отсутствием.

3) Области изостатической компенсации связаны с районами7 где наблюдаются положительные и близкие к нулю поправки к скоростному разрезу Джеффриса. Там, где наблюдаются слабые положительные изостатические аномалии, на глубине 200-250 км определены отрицательные поправки к скорости, и вероятно расположена зона пониженных скоростей.

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

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

6) В Восточных Карпатах и во Внешней и во Внутренней зоне граница Мохоровичича расположена на большой глубине.

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

ЗАКЛЮЧЕНИЕ

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

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

3. В результате решения обратной двумерной кинематической задачи для поверхностных волн по данным о средних групповых скоростях на трассах, пересекающих исследуемую область, получены распределения сглаженных значений скоростей волн Релея и Лява в интервале периодов 10-30 сек. Они характеризуются следующими особенностями: а) На периодах 10-20 сек максимальные скорости волн обоих типов приурочены к центральной части Балканского полуострова и

Малой Азии, Минимальными скоростями характеризуются области Черного моря, Предкарпатского прогиба, Эгейского моря. б) С увеличением периода горизонтальные вариации скорости сглаживаются, наиболее выраженными структурами остаются Предкарпатский прогиб (минимум) и Малая Азия (максимум).

-иов) Распределения скоростей волн Лява на периодах 10-20 сек хорошо коррелируют с распределениями скоростей волн Релея на периодах 20-30 сек, что является следствием большей глубины проникания волн Лява.

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

5. Статистический анализ невязок времен пробега Р-волн в районе Юго-Восточной Европы, Малой Азии и Восточного Средиземноморья позволил выявить значительные аномалии, которые находятся в пределах (-5,-бф (+Зу+4 с).

6. На основе полученных данных о невязках времен пробега была решена обратная трехмерная кинематическая задача и получены сглаженные значения поправок к скоростному разрезу Джеффри-са в исследуемом районе в интервале глубин 50-300 км. Распределения поправок характеризуются следующими особенностями: а) Областями пониженных значений скорости по отношению к стандартному разрезу Джеффриса являются Паннонский бассейн, западная часть Малой Азии, Закавказье. Поправки достигают максимальных по абсолютной величине значений на глубине 200 км и рав ны в этих областях соответственно -0,13, -0,14, -0,25 км/с. б) Областями повышенных значений скорости Р-волн являются Восточное Средиземноморье, южная часть Малой Азии (Тавриды), Черное море и Карпаты, Динариды и Пелагониды. Максимальные значения поправок приурочены к глубине 100 км и соответственно равны 0,25, 0,18 и 0,18 км/с. в) Области пониженных значений скорости соединены между собой зоной отрицательных или близких к нулю поправок. Они расположены под центральной частью Балканского полуострова, над Северо-Анатолийским разломом и срединными массивами центральной части Малой Азии. г) На глубине 200-250 км в Юго-Восточной Европе происходит значительное расширение областей, в которых отмечаются отрицательные поправки к скорости, что свидетельствует о наличии на этих глубинах зоны пониженных скоростей.

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

-1GZ

1. Алексеев A.C. и др. Численный метод определения структуры верхней мантии Земли. - В кн.: Математические проблемы геофизики, вып.2, Новосибирск, ВЦ СО АН СССР, 1969, с.143-165.

2. Аки К. »Ричарде П. Количественная сейсмология.М.,Мир,1983.

3. Аниконов Ю.Е.»Пивоварова Н.Б.,Славина Л.Б. Трехмерное поле скоростей фокальной зоны Камчатки. В кн.:Математические проблемы геофизики. Вып.5,ч.1,Новосибирск, ВЦ СО АН СССР, 1974, с.92-117.

4. Алашнин А.М. О распределении гравитационного поля на доставляющие при изостатическом равновесии масс. Изв.АН СССР, Физика Земли, № II, 1984, с.3-10.

5. Артемьев М.Е. Изостазия территории СССР. М., Наука,1975.215с.

6. Артемьев М.Е.и др. Плотностные неоднородности под системой островных дуг. Изв.АН СССР, шизика Земли, № I, 1982, с.3-25.

7. Бармин М.П. и др. Исследование горизонтальных неоднороднос-тей в литосфере при помощи поверхностных волн. Тезисы докладов Евр.Сейсм.Комм.XIX Ген.Ас.Москва, 1984, с.108.

8. Белявский H.A. Строение земной коры континентов по геолого-геофизическим данным. М., Недра, 198I. 432с.

9. Берч Ф. Плотность и состав верхней мантий (первое приближение при предположении оливинового шэстава). В кн.: Земная кора и верхняя мантия. М., Мир, 1972, с.27-56.

10. Бончев Е. Основни геодинамични проблеми на нашето пространство. В кн.:Геодинамика на Балканите. София, Техника, 1980, с.83-94.

11. Брюн Д.Н. Поверхностные волны и строение коры. В кн-¡Земная кора и верхняя мантия. М., Мир,1972, с.186-198.

12. Буачидзе Г.И. Распределение температуры в земной коре в районе Западного Кавказа и Черного моря. В кн.: Тепловое поле Европы. М., Мир, 1982, с.302-311.

13. Бугаевский Г.Н. Сейсмологические исследования неоднороднос-тей мантии Земли. Киев: Наукова думка, 1978, 184с.

14. Велинов Т.,Бояджнева К.,Петков И. Геетермические данные по Болгарии. В кн.: Тепловое поле Европы. М., Мир, 1982.

15. Величу С.»Деметреску К. Тепловой поток в Румынии и его связь с геологическими и геофизическими характеристиками. В кн.: Тепловое поле Европы. М., Мир, 1982, с.276-285.

16. Велчев Цв.,Дачев Хр.,Петков Ив. Морфологические структуры в Болгарии и их отражение в геофизических полях. д-ео*. Acad. Sol. Иино^. ' V. 15, 1971, р.357-364.

17. Вуллард Дж.П. .Региональные аномалии силы тяжести. В кн.: Земная кора и верхняя мантия. М., Мир, 1972, с.264-290.

18. Гобаренко В.С. Решение обратной трехмерной кинематической задачи сейсмологии методом Бэйкуса-Гильберта (на примере региона Центральной Азии). Канд.дисс.Киев, 1983, 167с.

19. Гобаренко В.С.,Яновская Т.Б. Исследование горизонтальных не-однородностей строения верхней мантии в Саяно-Алтайской зоне. Изв.АН СССР, Физика Земли, № 4, 1983, с.21-35.

20. Гончаров В.П0'Непрочное Н).П . Комплексное исследование Черноморской впадины. М., Наука, 1976,

21. Гравиразведка. Справочник геофизика. М.,Недра, 1981, 397с.

22. Грушинский Н.П. Теория фигуры Земли. М.,Наука,1976.512с.

23. Дачев Хр. Геодинамични проблеми на Балканския полуостров в светлината на регионалните геофизични данни. В кн.¡Геодинамика на Балканите. София, Техника, 1980, с.9-25.

24. Дачев Хр.,Петков Ив. Получение сейсмических волн глубинного происхождения на территории НРБ и опыт их интерпретации. В кн.: Геофизические исследования Юго-Восточной Европы, 1967.

25. Добрев Т.Б.,Щукин Ю.К. Геофизические поля и сейсмичность восточной части Карпато-Балканского региона. М.,Наука,1974, 170с.

26. Драгашевич Т. Результаты комплексного исследования земной коры в Югославии. В кн.: Строение земной коры и верхней мантии по данным сейсмических исследований. Киев, Науковадумка, 1977, с.185-194.

27. Драгашевич Т. Некоторые результаты глубинного сейсмическогозондирования в Динаридах. В кн.: Геофизические исследования строения земной коры Юго-Восточной Европы. М., Наука, 1967, с.56-63.

28. Дрейф континентов. М., Мир, 1966, 232с.

29. Жарков В.Н.»Трубицын В.П. Физика планетарных недр. М., Наука, 1980. 448с.

30. Земная кора и история развития Средиземного моря. М., Наука, 1982. 205с.

31. Кунин И.Я.,Шейх-Заде Э.Р. Латеральные неоднородности поверхности мантии в Западной Европе. Изв.АН СССР, Физика Земли, № 6, 1984 , с.21-36.

32. Лаврентьев М.М.»Романов В.Г. 0 трех линеаризованных обратных задачах для гиперболических уравнений. Докл.АН СССР,1966, 171, № 6, с.1279-1281.

33. Ландер А.В. Аномальные явления в поверхностных волнах насеверо-востоке Евразии и их связь с районом Мамского рифта. В кн.: Математическое моделирование и интерпретация геофизических данных. М., Наука, 1984. (Выч.сейсмология.Вып.16).

34. Левшин А.Л. Поверхностные и каналовые волны. М., Наука, 1973,

35. Ле Пишон К. »Франшто Ж.,Бонин Ж. Тектоника плит. М.» Мир, 1977, 288с.

36. Магницкий В.А. Внутреннее строение и физика Земли. М., Недра, 1965. 380с.

37. Морган П. Тепловой поток Кипра и замечания о геотермическом режиме Восточного Средиземноморья. В кн.: Тепловое поле Европы. М., Мир, 1982.

38. Москаленко В.Н. Строение земной коры Средиземного моря по сейсмическим данным. М., Наука, 1981. 112с.

39. Пошгай К. Сейсмические отражающие * границы и распределение скоростей в земной коре и мантии. В кн.: Строение земной коры и верхней мантии по данным сейсмических исследований. Киев, Наукова думка, 1977, с.228-230.

40. Романов М.Е, Численные методы решения обратной кинематической задачи для горизонтально-неоднородных сред. В кн.: Численные методы в сейсмических исследованиях. Новосибирск, Наука, 1983, с.70-106.

41. Реннер С. Штегена А. Ыссле^обешис.пол* е> Паинонскок бассейне. Ас£а . Ас.а<?( . £с<' и^к^г. , V. Ю, 43£>С } р.2:50-2&/{.

42. Рыжикова Сн.,Петков Ив. Строение земной коры в районе Черноморского бассейна но дисперсии групповых скоростей волн Лява и Релея. Киев, Наукова думка, Геофизический сборник, 1977, вып.80, с.24-32.-ли

43. Саваренский Е.Ф. Руководство по обработке поверхностных волн.

44. Саваренский Е.Ф.,Кирнос Д.П. Элементы сейсмологии и сейсмометрии. М., Гос.изд.те^Нг-теор.лит. ,1955, 544с.

45. Садиков Ф.,Косарев Г. Строение Земли по дисперсии поверхностных сейсмических волн. Ташкент, Фан, 1981, 102с.

46. Страхов В.И.,Романгок Т.В. Восстановление плотности земной коры и верхней мантии по данным ГСЗ и гравиметрии. Изв.АН СССР, Физика Земли, № 6, 1984, с.44-64.

47. Структура земной коры Центральной и Восточной Европы по данным геофизических исследований. Киев, Наукова думка, 1980, 208с.

48. Тихонов А.Н.,Арсенин В.Я. Методы решения некорректных задач. М., Наука, 1979. 285с.

49. Тезддан А.К. Геотермические исследования, их современное состояние и возможности картирования темлового потока в Турции. В кн.: Тепловое Поле Европы, М., Мир, 1982, с.300-322.

50. Фитикас М.Д.,Колиос Н.П. Предварительная карта теплового потока в Греции. В кн.: Тепловое поле Европы. М., Мир, 1982, с.200-210.

51. Хорват Ф.,Бодри Л.,0ттлик П. Геометрия Венгрии и тектонофизи-ка "красного пятна" Паннонской впадины. В кн.: Тепловое поле Европы. М., Мир, 1982, с.176-190.54. 1^бои Т. Гравитационное поле Земли. М., Мир, 1982, 288с.

52. Чермак В. Карта теплового потока Европы. В кн.: Тепловое поле Европы. М., Мир, 1982, с.11-54.

53. Яновская Т.Б. Метод решения обратной кинематической задачи сейсмики для горизонтально-неоднородных сред. В кн.: Методы и алгоритмы интерпретации сейсмологических данных. Выч.

54. Сейсмол. вып.13, M., Наука, 1980, с.96-101.

55. Яновская Т.Б. Распределение групповых скоростей поверхностных волн в Северной Атлантике. Изв.АН СССР, Физика Земли,2, 1982, с.3-И.

56. Яновская Т.Б.»Дорохова Л.Н. Обратные задачи геофизики. Ленинград, изд.ЛГУ, 1983.

57. Яновская Т,Б. Методы решения обратной кинематической задачи сейсмологии для двумерно и трехмерно-неоднородных сред.

58. В кн.: Численные методы в сейсмических исследованиях, Новосибирск, Наука, 1983, с.79-90.

59. Яновская Т.Б.»Николова С.Б. Распределение групповых скоростей поверхностных волн Релея и Лява в Юго-Восточной Европе и Малой Азии. Болгарский геофизический журнал, № 4, София, 1984, с.83-94.

60. Sçu stAA. Soc. Л ил i ; v/. 5" 5 } Zf 4 969 ; p. 4 407- — 4 4 . ÉbacJ¿^c, б-. ÍXuM (ritijLrJt F, Uni queues lu tu (.»л S Гоиaj- LU^ccKtaf^. gtoss cla.be,, PIaCI, Trc^s,. •

61. Soc . Lonoiok . Set. fi-V/, У U4SJ 15?0/ p.

62. Uai Me.ch'ire.'zïau&asL, см^ск Ae-^-c^. . Ik.:

63. Jnije-i. ил^ . о и flu. S-friue*. Ц IbioïAfi of- JW. ^fug . , Spl,'ir} 4S7-6. , p. 72>4- * Hi. H. Piuv,cft; X.} MoitilC c. Qe-opLyitcAJl e-Kptotej-ioi* of HeAi'fjLï'Tc^&u,

64. Sесч. ßo^. (^eo^îs. Тес*. A^., V. 4371 f

65. А- АЛ ^fotVc. a.* a. jpa&oyeo^pbScplLWxi'ie. A^pTW. Ik Jut. s«^ uf, oj- fiyjL Sb .of -KvíL M.Qcl , ¿cAibb.Spttir,7Ü, Hg^vcl^t F. S^be^e-ua L . Tkc. Pawkûuicuv ^sin. ; G:

66. UítckU Vu-Ь-г. Хч Гч^. S^ n-p. o^ -Kvc St4?. W»4s+. of Ujl NW- éft&iuA , 49Í-6 ( p . Ъ2>Ъ-1Ч0ко\)£алл<А 3 Ht^s^.^/. E.S. Tk^e. -cOw^^okci MoXodX t-^. Cuucuuc^ I иллас^е. -fí^ uppe/t

67. SOU.4-A, СХЗ^СЛхЬ Хк XiAolenh' j-ГC-Cuf-Oov, oj1

68. Sou^.ols> ELo^-t-k^U-cxAi-e о^ UndtJ^Yto^Kci En^ibv. fíus*^uíux-Kv Eai^t, Eut-oj^«. ТелЛя-в^и.} v. i3g р. Ib^-lSl.

69. Iut^uctt^ho/ seisi^^ícoi A^lUJrC^ . Ме^^л^.

70. Ы-ъиьиъ of -Кл Ай^Ч S^A.f^M* Ttc^ieX ТСgJL^clt+cdlQ íUiol toav/i^ . J к- A-ÊJ»^u^Out-S, , . W^.'oh CoivHv, . ^IM Î<U.

71. Mus ^ g-Ao4-406. SA. 3ôlv'nso^ L.E., frCfciJ- F. Т«луел-&х'оч IU^IÍAÜÍ^LM

72. Se^SrU^c. za^ , . J^I ieut^. PUuj s . , V/. /2-^

73. Ha £ftVb-U M^c/ ; f-^^AUot^ . EÄVH .ft P^.

74. Le-H-. ; v.24 ( /Í9?>, HaMJbCç, 3. Ctu^W S^-vc¿>j- Кл. S^A а^гЦ4Э?Г; v.

75. Mat¿.t.\'S t3 . TU- UtMJs»6=A Sbu^cÁ^U- co^cá ирр,е.г b^a^bÔL of t^ Ac^tAUv -сл-^cou f DBS . T^ctouö^lv.} v. 4 97-8 ^.¿£5-2.8^. Reilos» 3. Iuv^f^a^'ot, o f . Л377-.-м-о

76. SS. Moil£ едЧ ^ . Sovwe. coubCcle-tJktcoub Ok -К*. ^гоДiWwk,

77. Situais U (ri^ce. Ttt+b^k, «w,gg Me-Uü-Utl'e. D. Слл, frtcÁjSL -ЬгЫъи^С. ít^ttíl«. ÛOIa tV«*«-«* ЬЛu^fttíoM . lu'. I«vb. S^M^.OU Ha. M.&UuS.1. Split, ,p. 485-46*.

78. M^tit n>. AcKUJí. of МхД.'^^иаачreo^U . 3 . Ro^. Aib. Soc. , V. So, Л9?2.,/*. 4d3->!8S\tvrtj„Je.íCAb404 6 . TU-Wo^t., V. ; /(9 p. 332-344 . 62. и^ъч A.E.M , H.W. , S-t«-K*-l F. 6-. (£cl) TU. осдли

79. OMOI uuQ^i'uS . TkjL. Ea^+tA^i el .'"fcV? if P^SS,1. KW-Vox^, A A. , 4

80. Ow^i Й. A- Gtc.vVfcty Аи.ои^а£Гел Ил. Au^Ukll«-.

81. Uí-aiA. M.V. . ave/ ef- Тч.ч1оuj. ; i/. /< g > j p5ч. ОьгЛ^Г F. A ou Ka jxjssí s«^ Гсиичсх af

82. Cßu^i'c. ctc-f \ W ^ ^ tUxL f of tAjL EaS-fe/t^

83. UíA4 37-3 , v/. AS , p. 22?-23"?-.

84. OS". РйкЪч G-.P. TLg. рои/¿-"i <»j- S^v^-i'c. Ьч-г^й-се Viíu/«vui'tU^ |э.есЛг 4-c ОъчЛ4- сид.с( и^р^Л. l^a^f-Ct bt"¿M.<^i~u.ieJl1.ûd^S» . IK'. "TV "bo^ti'oM oj- fluu. f^og-ÔL-U,

85. Ua . С ¿ci, ß . Ca<U«» »• P^SS, tUeMj-Yü\.\c

86. P*J,a"bcuLúS. B.C. ; Ро^Ьк M-.'ucU&S M. ûf+Ч

87. Лоа. Po-^ о . 64 u<; fcJ? S, f^xvcjb-t tí of- Ь*

88. Uï-j-^c* iL . PkaS«. vc/oc^-W/Lî. ç^cL + l^ueJL-j-C'^S-S.

89. S.OC. А-^л , v. ¿T6 ( р. 2 Ъ-liZ.

90. ЛОА, Eejör U-o>« f fc.li P. Л-И w. 6. Р. ^-гл.^»' Ou^ci

91. Rjlc. Ht^üf. .ъя-1 p. ¿-<S1оЧ. V. j So^füT. WefcjW (¡V04 o-f<j/ti)i4.|a4о!Г. P • , G-i-^atöVc; E. TU. c.tt-tS'K^ sf-^ofvt^of

92. Kl WeA+^U Бий^с,"ч-vq íco-í сАа^гг.оил.с< ^«-oát*. ^ /S V/. ( р. 3 7-, AOÊ. Wooot tí. ^ С . ö-tcwi'f^j Сs-lj-ts-tcjí s {"гял.сйл.ч^. »"к «-сч-ьдтг^ц lhax,

93. Oj.'oi*. SoC. A*i . j V/. Я ,c, u*.ilc. imjla.G-JLOJA, feiLS. ; \/.¿4{los. YavooVbi^0^^. ть . "H^ ръо^-^-иof" c.cx/^mc&xj^ ia^ltcJlQ^ ^jlcLso, .ел^Ць. ß. Лй"^. Soa.j Лд&Ч, v. J-*} f.ZIZ -иУЧ-т