CC BY
16УДК 550.348.435
СПЕКТРАЛЬНЫЕ и ОЧАГОВЫЕ ПАРАМЕТРЫ ЗЕМЛЕТРЯСЕНИЙ КРЫМСКО- ЧЕРНОМОРСКОГО РЕГИОНА Б.Г. Пустовитенко, Э.Э. Эреджепов
Институт сейсмологии и геодинамики ФГАОУВО «Крымский федеральный университет имени В.И. Вернадского», г. Симферополь, [email protected], [email protected]
Аннотация. Представлены спектральные и динамические параметры очагов (Мо, го, До, 8, по, О г , и, Eu и Ым>) тринадцати землетрясений Крыма с Кп=7.2-11.2, Ым>рег=2.6-4.1, восстановленные по 123 амплитудным спектрам продольных и поперечных сейсмических волн, зарегистрированных цифровыми региональными сейсмическими станциями. Аппроксимация спектров и расчет по их характеристикам динамических параметров очагов землетрясений выполнены в рамках дислокационной модели Бруна. В осреднении очаговых параметров участвовало от четырех до семи станционных определений, что обеспечило малую величину стандартного отклонения - показателя степени рассеяния индивидуальных оценок. Наилучшая сходимость станционных определений отмечена для радиуса круговой дислокации. Наибольшие значения динамических параметров получены для ощутимого землетрясения 2 марта с Кп= 11.2, а наименьшие - для его слабого афтершока 6 апреля с Кп=7.2. Во всем диапазоне энергий величина сброшенных напряжений не превысила Да=106 Па (10 бар), а кажущихся напряжений по<3-105 Па (3 бар). Средние значения сейсмических моментов и радиусов круговой дислокации в пределах погрешностей их определения удовлетворяют средним долговременным зависимостям параметров от энергетического уровня землетрясений, полученным по аналоговым записям. Для наиболее сильных землетрясений 2 марта с Кп=11.2, М^рег=4.1 и 18 октября с Кп=10.1, М^рег=3.8 получены решения механизмов очагов. Землетрясение 2 марта произошло в центральной части региона вблизи Южного берега Крыма под действием горизонтальных сил растяжения, ориентированных близширотно. Тип подвижки в очаге - сбросо-сдвиг, с преобладанием сбросовой компоненты над сдвиговой. Землетрясение 18 октября возникло в центральной части Черноморской впадины в обстановке субширотного сжатия и субмеридионального растяжения. Тип подвижки в очаге - чистый сдвиг. Для главного толчка 2 марта и его шести афтершоков по данным станции «Алушта» приведены энергетические спектры и их основные параметры. Дан анализ полученных результатов.
Ключевые слова: механизм очага, амплитудный спектр, модель Бруна, сейсмический момент, радиус дислокации, сброшенное и кажущееся напряжения, подвижка по разрыву, радиационное трение, энергетический спектр.
DOI: 10.35540/1818-6254.2020.23.25
Для цитирования: Пустовитенко Б.Г., Эреджепов Э.Э. Спектральные и очаговые параметры землетрясений Крымско-Черноморского региона // Землетрясения Северной Евразии. - Вып. 23 (2014 г.). -Обнинск: ФИЦ ЕГС РАН, 2020. - С. 250-262 doi: 10.35540/1818-6254.2020.23.25
Исходные данные. Большинство (около 70 %) из всех зарегистрированных в 2014 г. землетрясений [1] относится к слабым (КП<7) сейсмическим событиям, непредставительным для региона [2]. Для таких землетрясений основные параметры, особенно глубины очагов, определены с меньшей надежностью, только по данным ограниченного числа ближайших к очаговой зоне сейсмических станций. В связи с этим их записи не были привлечены для восстановления очаговых параметров. Для анализа выбраны только все наиболее значимые представительные сейсмические события года с Кп>9.5, а также более слабые толчки из зоны ощутимого землетрясения 2 марта с КП=11.2 и районов недостаточно изученных территорий северо-западной части Крыма и Азовского моря (районы № 6, № 7). Для предварительно отобранных землетрясений был проведен анализ качества сейсмограмм на сейсмических станциях с надежными амплитудно-частотными характеристиками приборов. Выбраны записи с четкими волновыми формами объемных P- и «-волн, не осложненные микросейсмическим шумом, с превышением полезного сигнала над фоном помех в два и более раза. По такому критерию для расчета спектральных и динамических параметров очагов отобрано 13 землетрясений с энергетическими классами Кп=7.2-11.2 (рис. 1, табл. 1).
Для двух наиболее сильных землетрясений: первого - 2 марта, произошедшего в центральной части региона вблизи Южного берега Крыма с Кп=11.2 и интенсивностью до трех баллов в городах Ялта и Алушта [2], и второго - 18 октября с Кп=10.1, возникшего в центральной части
Черноморской впадины, получено решение механизма очага [3]. Для этого использованы записи продольных Р-волн по станциям Крымской, Турецкой и Кавказской сетей, а также данные о знаках первых вступлений Р-волн, приведенные в бюллетенях международных сейсмологических агентств [4-6]. Углы выхода сейсмических волн определены для удаленных станций по годографу «ак135» [7], а для близких - по региональному годографу [8].
Рис. 1. Карта эпицентров землетрясений Крымско-Черноморского региона за 2014 г., для которых
восстановлены очаговые параметры
1 - энергетический класс Кп [9]; 2 - номер района; 3 - сейсмическая станция; 4 - граница района; 5 - номер землетрясения по табл. 1; 6 - глубина очага h, км.
Максимальное количество изученных землетрясений (семь событий) относится к Ялтинской зоне региона (район № 2), два - к Азово-Кубанской (район № 7) и по одному -к Алуштинской, Керченско-Анапской зоне, Степному Крыму и Черноморской впадине
(районы № 3, 5, 6 и 9 соответственно).
Таблица 1. Основные параметры землетрясений Крымско-Черноморского региона за 2014 г., для которых определены очаговые параметры
№ Дата, к, Эпицентр К, Магнитуда Кп Район
д м ч мин с ф°, N Е км ть Мм'рег/Л (табл. 4)
1 2 3 4 5 6 7 8 9 10 11
1 04 01 07 48 00.3 46.14 36.57 6 3.6 3.52/7 3.3 9.4 Азово-Кубанский (№ 7)
2 02 03 03 34 28 3 44.30 34.33 28 3.4 4.12/7 3.9 11.2 Ялтинский (№ 2)
3 02 03 18 10 25.9 44.26 34.34 26 - 2.88/4 - 8.5 Ялтинский (№ 2)
4 17 03 23 13 36.0 44.43 34.34 19 - 3.31/7 3.0 9.8 Ялтинский (№ 2)
5 25 03 17 04 40.5 44.69 36.53 13 3.3 3.32/7 3.2 9.6 Керченско-Анапский (№ 5)
6 29 03 05 41 21.1 44.44 34.33 19 - 3.05/7 - 8.6 Ялтинский (№ 2)
7 30 03 07 55 59.1 45.51 32.98 21 - 3.13/4 - 8.2 Степной Крым (№ 6)
8 06 04 20 47 13.2 44.30 34.28 23 - 2.63/5 - 7.2 Ялтинский (№ 2)
9 27 04 02 51 33.6 44.23 34.33 22 - 2.9/6 - 8.5 Ялтинский (№ 2)
10 28 04 06 30 35.8 44.62 34.47 22 - 3.44/5 2.9 9.3 Алуштинский (№ 3)
11 20 08 17 05 27.2 44.49 34.25 15 - 2.83/7 - 7.9 Ялтинский (№ 2)
12 18 10 15 34 54.0 42.74 35.11 5 3.7 3.83/7 3.7 10.1 Черноморская впадина (№ 9)
13 09 11 04 08 13.7 45.51 37.07 25 - 3.05/4 - 8.5 Азово-Кубанский (№ 7)
Примечание. Параметры землетрясений в графах 2-7, 9-11 соответствуют таковым в [1]; значения Мм>/п - из табл. 4.
Наибольшее количество спектров получено по записям станции «Судак». Процент участия этой станции в общей оценке динамических параметров отдельных землетрясений максимальный и составляет 100 %, для остальных также высокий - 85 %.
Пока не использованы для анализа записи сейсмических станций «Феодосия» и «Керчь», которые расположены, соответственно, на мергелистых глинах и мшанковых известняках с пониженными прочностными свойствами и повышенным затуханием высоких частот. Сейсмические колебания от местных землетрясений на этих станциях не только осложнены микросейсмическим шумом, но и сложным образом искажены амплитудно-частотной характеристикой (АЧХ) среды, которую в настоящее время без специальных исследований корректно учесть не представляется возможным. Записи сейсмической станции «Тарханкут», открытой в 2012 г. [10], также не привлечены для анализа из-за отсутствия надежных АЧХ сейсмографов и АЧХ локальных инженерно-геологических условий места расположения станции. На момент расчета спектров и оценки динамических параметров не были уточнены основные параметры сейсмографов станции «Ялта» за 2014 г., в связи с этим сейсмограммы этой станции также не использованы в данной работе.
Механизмы очагов землетрясений [3] получены с использованием теории дислокации по стандартной методике в рамках двойного диполя [11].
В соответствии с полученными решениями (табл. 2, рис. 2) землетрясение 2 марта (№ 2 по табл. 1) произошло под действием горизонтальных сил растяжения (PZT=8°), ориентированных близширотно (AZMt=1 17°). Тип подвижки в очаге - сбросо-сдвиг с преобладанием сбросовой компоненты над сдвиговой (табл. 2, рис. 2). Первая плоскость разрыва NP1 имела диагональную ориентацию (<S7Xnpi=44°) с наклонным падением (DPnpi=57°) на юго-восток, а другая NP2 - близмеридиональную (<STKnp2=184°), с пологим падением (DPnp2=41°) на запад.
Таблица 2. Параметры механизма очагов землетрясений 2 марта 2014 г. в 03ь34т с Кп=11.2 и 18 октября в 15ь34т с Кп=10. 1
Дата, д м to, ч мин с h, км Магнитуды Оси главных напряжений Нодальные плоскости Ис-
Mw ть Кп T N P NP1 NP2 точ-
(MLwsg) PL AZM PL AZM PL AZM STK DP SLIP STK DP SLIP ник
02.03 03 34 28.3 28 4.1 3.4(3.9) 11.2 8 117 21 210 68 5 44 57 -65 184 41 -122 [121
18.10 15 34 54.0 16* 3.8 3.7(3.7) 10.1 28 167 60 326 9 72 302 77 27 206 64 166 [12]
Примечание. При расчетах углов выхода использована глубина очага по данным ISC [5].
Землетрясение 18 октября (№ 12 по табл. 1) возникло в обстановке восток-северо-восточного сжатия (А2Мр=72°, PLp=90) и субмеридионального растяжения (42Мр=167°, PLт=280). Тип подвижки в очаге - практически чистый сдвиг. Обе нодальные плоскости близвертикальны (ОРт\ =77° и £РОТ2=64°) и простираются диагонально (БТКтх =302° и «ТК^ =206°).
0 0
Рис. 2. Диаграммы механизма очагов землетрясений в проекции нижней
полусферы: а - 2 марта в 03ь34т с Кп=11.2; б - 18 октября в 15ь34т с Кп=10.1
1 - нодальные линии; 2, 3 - оси главных напряжений растяжения и сжатия соответственно; зачернена область волн сжатия.
По параметрам механизма очагов определены поправки за направленность излучения Reф из очага на станцию регистрации (табл. 3), которые далее использованы при оценке скалярного сейсмического момента М0 для этих землетрясений по формуле (1).
270
90 270
180 а
180 б
Таблица 3. Значения станционных поправок для землетрясений 2 марта и 18 октября 2014 г.
Дата, Станция $ЕУ SUDU SUDU ALU ALU
д м Тип волны 5 5 Р 5 Р 5 5
02.03 Rеф 0.54 0.32 0.21 0.39 0.56 0.27 0.54
18.10 Rеф 0.19 0.18 0.73 0.28 0.65 0.23 0.19
Методика обработки и интерпретации амплитудных спектров. Амплитудные спектры поперечных волн (5) во всех случаях рассчитывались по двум горизонтальным составляющим (Ы^) и (Е^) и вычислялся полный вектор колебаний, а продольных (Р) - по вертикальной составляющей (2). Относительная длительность т исследуемого участка записи принята равной интервалу времени от момента вступления 5- и Р-волн до времени спада их максимальных амплитуд Атах на уровень 1/з Атах [12].
Станционные амплитудные спектры рассчитаны стандартным методом быстрого преобразования Фурье. Для интерпретации выбрана теоретическая дислокационная модель Бруна (ш-2) [13], согласно которой амплитудный спектр объемной волны характеризуется тремя основными параметрами: спектральной плотностью О о(х, / в длиннопериодной части (при /^0),
угловой частотой /0(© о) и углом наклона у--2. Полученные амплитудные спектры записей
землетрясений Крыма за 2014 г. уверенно аппроксимируются двумя прямыми: в длиннопериодной части спектра - прямой, параллельной оси частот и наклонной, с угловым коэффициентом у—2 в области коротких периодов (рис. 3), что соответствует выбранной теоретической модели.
О %я.п
О лм.п
10
10н
10"7
10"
10"
1—1—........ 8иои, Е |-1—гт
101 10° КУ/, Гц 101 10° КУ/, Гц 101 10° Ю^Гц Ю1 10° Ю1/,Гц
4 января в 07М8т0.38(№ 1) 2 марта в 03ь34т28.38(№ 2) 30 марта в 071г55т59.18(№ 7) 28 апреля в 6ь30т35.88(№ 10)
Рис. 3. Примеры амплитудных спектров объемных сейсмических волн землетрясений Крыма за 2014 г. по записям региональных цифровых сейсмостанций и их аппроксимации в рамках теоретической модели Бруна (номера и даты землетрясений соответствуют таковым в табл. 1)
Расчет динамических параметров очага. Согласно теоретической дислокационной модели Бруна [13], спектральная плотность О 0 пропорциональна скалярному сейсмическому моменту М0, а угловая частота /0 прямо связана с размером дислокации Г0.
Формула для расчета М0 по станционному спектру имеет вид:
М 0=О 0(4пр •V3)/Д9ф • G (А, h) С (©) • 5т (/), (1)
где О0 - максимальное значение спектральной плотности при V - скорость распространения
объемной волны; р - плотность пород в окрестности очага, й0ф - направленность излучения из очага на станцию регистрации; G(А, К) - поправка за геометрическое расхождение; С(©) - частотная характеристика среды под станцией; 5т(/) - поправка за неупругое затухание в мантии.
Методика учета всех этих факторов при переходе от станционного спектра к спектру источника подробно изложена в работах [12-15].
Геометрические размеры очага определялись по угловой частоте /0 амплитудного спектра смещения. Для дислокационной модели Бруна с разрывом в виде круга, радиус дислокации Г0 вычисляется по формуле:
г = 2.34
Ур
/0
По найденным значениям сейсмического момента Мо и размеров разрыва го с использованием формул из работ [12-17] определены другие динамические параметры очагов: сброшенное напряжение Ас, величина деформации сдвига в, кажущееся напряжение , величина радиационного трения с r, сейсмическая энергия Es , средняя подвижка по разрыву u (или величина
дислокации), энергия образования дислокации в очаге Ей и моментная магнитуда Мнуг (табл. 4). Моментная магнитуда рассчитана по индивидуальным оценкам скалярных сейсмических моментов Мо с использованием формул из [18]:
Mw=2/3-lgMo-10.7, (3)
Mw=2/3-(lgM0+7)-10.7, (4)
где Mo имеет размерность [дин^см] в системе СГС и [Нм] - в системе СИ.
Большинство динамических параметров получено по группе независимых станционных определений.
Осреднение станционных динамических параметров (кроме радиационного трения Асг) и вычисление их стандартных отклонений проведены, как и прежде, с учетом логнормального закона распределения величин [12]. Станционные значения Ас г получились знакопеременными, поэтому их среднее значение для данного очага вычислено по среднегеометрическим напряжениям Ас и ^с по формуле Асг=1/2 Ас - ^с из [16]. В этом случае графа 5S в табл. 4 не заполнялась. Среднее значение моментной магнитуды Mw^ определено как среднее арифметическое с соответствующей погрешностью.
Результаты расчета станционных и средних для землетрясения динамических параметров представлены в табл. 4. Для каждой станции указаны эпицентральное расстояние А, км, тип использованной волны (Р, S) и составляющая записи, где (N+E) - полный вектор колебаний по N-S и E-W.
Таблица 4. Спектральные и динамические параметры очагов землетрясений Крыма за 2014 г.
Код ст Тип волны Сост. А, км Ü0-10-6, м-с SÜ0-10-6, м-с А Гц M0-1013, H-м Г0, км Ао-105, Па е-10-6 й-10-2, м г|о-105, Па Аог-105, Па Eu-108, Дж MWper
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Землетрясение (1) 4 января: ¿0=07h48m0.3s; ф=46.14°, Х=36.57°; h=6 км; Кп=9.4
ALU P Z 234 0.041 - 3.2 9.43 0.58 2.1 7.0 0.3 1.59 -0.55 3.3 3.25
ALU S N+E 234 - 0.41 2.1 20.4 0.53 5.9 19.7 0.76 0.74 2.2 20.1 3.48
SEV S N+E 287 - 0.39 2.1 23.7 0.53 6.9 23.0 0.89 0.63 2.8 27.3 3.52
SIM P Z 233 0.1 - 3.25 11.4 0.57 2.66 8.9 0.37 1.31 0.02 5.1 3.31
SIM S N+E 233 - 1.41 2.1 34.9 0.53 10.1 33.8 1.31 0.43 4.63 58.9 3.63
SUDU P Z 186 0.22 - 3.1 40.2 0.60 8.11 27.0 1.18 0.37 3.68 54.2 3.67
SUDU S N+E 186 - 1.44 2.0 56.8 0.56 14.2 47.5 1.93 0.26 6.86 135.0 3.77
s 23.65 0.56 5.95 19.87 0.81 0.63 23.65 23.52 3.52
SS 0.11 0.01 0.11 0.11 0.11 0.11 0.11 0.22 0.14
Землетрясение (2) ^2 марта: fo=03h34m28.3s; ф=44.30°, Х=34.33°; h=28 км;^Кп=11.2
ALU P Z 43 2.0 - 2.7 148.0 0.91 8.6 28.6 1.9 3.21 1.08 212.0 4.05
ALU S N+E 43 - 6 1.7 182 0.76 13.8 46.0 2.8 2.61 4.29 419 4.11
SEV S N+E 59 - 3 1.65 123 0.86 8.5 28.4 1.77 3.87 0.39 174 4.0
SIM P Z 75 0.6 - 2.8 280 0.88 18.1 60.5 3.86 1.7 0.74 848 4.23
SIM S N+E 75 - 12 1.52 189 0.93 10.2 34.1 2.31 2.52 2.6 322.0 4.12
SUDU P Z 84 0.5 - 2.5 220.0 0.98 10.1 33.7 2.41 2.17 2.89 370.0 4.16
SUDU S N+E 84 - 4.95 1.5 223.0 0.94 11.6 38.8 2.66 2.13 3.69 433.0 4.17
s 188.96 0.89 11.18 37.32 2.46 2.52 1.68 352.5 4.12
SS 0.04 0.01 0.04 0.04 0.04 0.05 0.15 0.09 0.04
Землетрясение (3) 2 марта: ¿0=18h10m25.9s; ф=44.26°, X=34.34°; h=26 км; Кп=8.5
ALU P Z 47 0.01 - 5.0 0.94 0.49 0.35 1.16 0.04 2.01 -1.84 0.05 2.59
ALU S N+E 47 - 0.15 2.9 2.69 0.49 1.01 3.38 0.12 0.7 -0.19 0.46 2.89
SIM S N+E 79 - 0.15 2.5 3.13 0.57 0.76 2.52 0.1 0.6 -0.22 0.39 2.93
SUDU S N+E 87 - 0.12 2.8 5.48 0.51 1.86 6.19 0.23 0.35 0.58 1.69 3.1
s 2.57 0.51 0.84 2.8 0.1 0.74 -0.32 0.35 2.88
SS 0.16 0.02 0.15 0.15 0.16 0.16 - 0.32 0.14
Код ст Тип волны Сост. Д, км 00-10-6, м^с И00-10-6, мс /0, Гц Мг1013, Нм Г0, км Да-105, Па е-10-6 й-10-2, м г|а-105, Па Даг-105, Па &-108, Дж MWрег
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Землетрясение (4) 17 марта. ¿0=23ь13т36.08; ф=44.43°, Х=34.34°; й=19, км; Кп=9.8
ALU Р ъ 28 0.15 - 3.6 7.06 0.65 1.14 3.81 0.18 4.25 -3.68 1.35 3.17
ALU 5 28 - 0.7 2.3 6.3 0.58 1.39 4.64 0.2 4.76 -4.06 1.46 3.14
SEV 5 54 - 0.4 2.85 9.05 0.47 3.8 12.7 0.43 3.31 -1.41 5.74 3.24
SIM Р ъ 60 0.14 - 3.9 9.09 0.6 1.87 6.23 0.27 3.3 -2.37 2.83 3.24
SIM 5 60 - 1.64 2.1 20.2 0.64 3.4 11.3 0.53 1.48 -0.21 11.5 3.47
SUDU Р ъ 73 0.1 - 3.8 15.7 0.61 3.0 9.97 0.44 1.91 -0.41 7.84 3.4
SUDU 5 73 - 0.8 2.15 24.0 0.62 4.3 14.5 0.65 1.25 0.92 17.4 3.52
£ 11.58 0.59 2.43 8.1 0.35 2.59 -1.38 4.69 3.31
55" 0.09 0.02 0.09 0.09 0.08 0.09 - 0.16 0.13
Землетрясение (5) 25 марта: /0=17ь4т40.58; ф=44.69°, Х=36.53°; h=13, км; Кп=9.6
ALU Р ъ 169 0.012 - 3.5 4.21 0.67 0.63 2.1 0.1 5.05 -4.73 0.44 3.02
ALU 5 169 - 0.15 2.3 10.0 0.58 2.2 7.4 0.31 2.1 -1.0 3.87 3.27
SEV 5 227 - 0.11 2.6 9.86 0.52 3.15 10.5 0.39 2.15 -0.58 5.15 3.27
SIM Р ъ 193 0.08 - 3.2 15.9 0.73 1.8 6.0 0.32 1.33 -0.43 4.82 3.4
SIM 5 193 - 0.32 2.05 12.2 0.65 1.9 6.3 0.30 1.74 -0.79 3.87 3.33
SUDU Р ъ 123 0.06 - 4.0 15.3 0.58 3.4 11.3 0.48 1.39 0.31 8.67 3.39
SUDU 5 123 - 0.6 2.2 29.3 0.61 5.65 18.8 0.84 0.73 -2.1 27.6 3.58
£ 12.01 0.62 2.24 7.48 0.33 1.77 -0.65 4.53 3.32
55 0.1 0.02 0.11 0.11 0.11 0.1 - 0.2 0.12
Землетрясение (6) 29 марта: ?0=05ь41т21.Р; ф=44.44°, Х=34.33°; h=19, км; Кп=8.6
ALU Р ъ 28 0.06 - 4.0 2.82 0.58 0.63 2.1 0.09 0.95 -0.63 0.29 2.9
ALU 5 N+E 28 - 0.46 2.5 4.14 0.54 1.17 3.9 0.15 0.65 -0.06 0.81 3.0
SEV 5 53 - 0.16 2.3 3.55 0.58 0.78 2.6 0.11 0.75 -0.36 0.47 2.97
SIM Р ъ 59 0.1 - 3.5 6.61 0.67 0.98 3.28 0.16 0.4 -0.09 1.08 3.15
SIM 5 N+E 59 - 0.43 2.15 5.43 0.62 0.98 3.26 0.15 0.49 -0.003 0.89 3.09
SUDU Р ъ 73 0.03 - 4.0 4.65 0.58 1.03 3.44 0.15 0.57 -0.06 0.8 3.05
SUDU 5 N+E 73 - 0.24 2.25 7.1 0.60 1.74 4.9 0.21 0.38 0.36 1.74 3.17
£ 4.68 0.59 1 3.25 0.14 0.57 -0.07 0.76 3.05
55 0.05 0.01 0.05 0.04 0.04 0.05 - 0.09 0.08
Землетрясение (7) 30 марта: ?0=07ь55т59.Р; ф=45.51°, Х=32.98°; h=21, км; Кп=8.2
SEV 5 121 - 0.07 2.9 3.35 0.46 1.48 4.95 0.17 0.4 0.3 0.83 2.95
SIM Р ъ 108 0.06 - 3.0 6.7 0.78 0.63 2.1 0.12 0.2 0.11 0.7 3.15
SIM 5 N+E 108 - 0.4 2.0 8.53 0.67 1.24 4.13 0.2 0.16 0.46 1.76 3.22
SUDU 5 N+E 173 - 0.12 2.25 8.2 0.6 1.7 5.65 0.25 0.16 0.68 2.32 3.21
£ 6.29 0.62 1.18 3.95 0.18 0.21 0.38 1.24 3.13
55 0.09 0.05 0.1 0.1 0.07 0.09 0.17 0.13 0.09
Землетрясение (8) 6 апреля: /0=20ь47т13.28; ф=44.30°, Х=34.28°; h=23, км; Кп=7.2
ALU Р ъ 44 0.01 - 4.2 0.68 0.55 0.18 0.58 0.02 0.25 -0.16 0.02 2.49
ALU 5 44 - 0.075 2.8 0.98 0.48 0.39 1.3 0.04 0.17 0.02 0.06 2.6
SEV 5 N+E 55 - 0.04 2.4 0.94 0.56 0.24 0.79 0.03 0.18 -0.06 0.04 2.6
SUDU Р ъ 87 0.008 - 4.0 1.49 0.58 0.33 1.11 0.05 0.11 0.05 0.08 2.72
SUDU 5 N+E 87 - 0.045 2.3 1.61 0.58 0.35 1.18 0.05 0.1 0.07 0.09 2.74
£ 1.08 0.55 0.29 0.95 0.04 0.15 -0.01 0.05 2.63
55 0.07 0.02 0.06 0.06 0.08 0.07 - 0.12 0.08
Землет рясение (9) 27 апреля: ¿0=02ь51т33.68; ф=44.23°, Х=34.33°; h=22, км; Кп=8.5
ALU Р ъ 51 0.02 - 5.0 1.52 0.47 0.66 2.2 0.07 1.25 -0.9 0.17 2.72
ALU 5 51 - 0.15 2.75 2.18 0.49 0.82 2.7 0.09 0.87 -0.46 0.3 2.83
SIM Р ъ 82 0.02 - 4.0 1.75 0.58 0.39 1.3 0.06 1.1 -0.88 0.11 2.77
SIM 5 N+E 82 - 0.25 2.85 4.19 0.47 1.76 5.9 0.2 0.45 0.43 1.23 3.02
SUDU Р ъ 90 0.018 - 5.0 3.44 0.47 1.49 5.0 0.17 0.55 0.2 0.86 2.96
SUDU 5 N+E 90 - 0.14 2.8 5.12 0.48 2.04 6.8 0.24 0.37 0.65 1.74 3.08
£ 2.75 0.49 1.02 3.4 0.12 0.69 -0.18 0.47 2.9
55 0.09 0.01 0.11 0.11 0.1 0.09 - 0.2 0.12
Землетрясение (10) 28 апреля: ¿0=06ь30т35.88; с р=44.62°, Х=34.47°, h=22, км; Кп=9.3
SEV 5 63 - 0.64 2.15 17.0 0.62 3.07 10.2 0.46 0.63 0.9 8.69 3.42
SIM Р ъ 46 0.39 - 3.8 20.4 0.61 3.83 13.0 0.58 0.53 1.42 13.2 3.48
SIM 5 N+E 46 - 1.56 2.1 16.6 0.64 2.63 8.75 0.41 0.68 0.63 6.84 3.4
Код ст Тип волны Сост. Д, км Пй-Ю"6, м-с ШтЮ-6, мс А Гц М0-1013, Нм Г0, км Да-105, Па е-10-6 й-10-2, м г|а-105, Па ДаЛ05, Па Еи108, Дж MWрег
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
SUDU Р ъ 51 0.1 - 4.2 11.4 0.55 2.93 9.76 0.39 0.93 0.53 5.56 3.31
SUDU 5 ^Е 51 - 1.4 2.2 30.5 0.61 5.89 19.6 0.87 0.35 2.6 29.9 3.59
£ 18.21 0.61 3.51 11.73 0.52 0.59 1.17 10.55 3.44
55 0.07 0.01 0.06 0.06 0.06 0.07 0.13 0.13 0.08
Землетрясение (11) 20 августа: /0=17ь05т27.28; ф=44.49°, Х=34.25°; h=15 км; Кп=7.9
АШ Р ъ 24 0.02 - 5.1 1.18 0.46 0.54 1.86 0.06 0.57 -0.29 0.11 2.65
АШ 5 ^Е 24 - 0.15 3.1 1.7 0.43 0.92 3.06 0.09 0.4 0.06 0.26 2.76
SEV 5 ^Е 46 - 0.1 3.15 1.92 0.43 1.1 3.6 0.11 0.35 0.19 0.35 2.79
SIM Р ъ 52 0.05 - 3.9 2.79 0.6 0.57 1.91 0.08 0.41 0.46 0.27 2.9
SIM 5 ^Е 52 - 0.2 2.5 2.13 0.54 0.6 2.02 0.08 0.31 -0.01 0.21 2.82
SUDU Р ъ 73 0.01 - 5.05 1.55 0.46 0.69 2.31 0.08 0.43 -0.08 0.18 2.73
SUDU 5 ^Е 73 - 0.21 3.0 6.22 0.45 3.05 10.2 0.33 0.11 1.42 3.2 3.13
£ 2.16 0.48 0.88 2.93 0.1 0.34 0.1 0.32 2.83
55 0.09 0.02 0.1 0.1 0.09 0.09 - 0.18 0.11
Землетрясение (12) 18 октября: ^=15ь34т54.08; ф=42.74°, Х=35.11°; h=15 км; Кп=10.1
АШ Р ъ 224 0.11 - 2.95 33.5 0.79 3.0 9.9 0.57 3.57 2.08 16.6 3.62
АШ 5 N+E 224 - 0.48 2.0 62.5 0.63 10.8 35.9 1.65 1.91 3.47 112 3.8
SEV 5 N+E 232 - 0.45 1.82 77.4 0.78 4.57 33.5 1.7 1.54 3.48 130 3.86
SIM Р ъ 255 0.3 - 3.0 48.8 0.62 6.88 15.2 0.86 2.45 -0.16 37.2 3.73
SIM 5 ^Е 255 - 1.2 1.8 107.0 0.70 13.5 45.0 2.3 1.11 5.64 242 3.96
SUDU Р ъ 240 0.26 - 2.95 73.8 0.79 6.57 21.9 1.26 1.62 1.67 80.8 3.85
SUDU 5 N+E 240 - 1.5 1.75 129 0.68 17.9 59.6 2.95 0.93 8.01 384 4.01
£ 69.78 0.71 7.71 27.06 1.42 1.71 2.15 94.54 3.83
55 0.08 0.02 0.1 0.1 0.09 0.08 - 0.18 0.08
Землетрясение (13) 09 ноября: ^=04ь08т13.78; ф=45.51°, Х=37.07°; h=25 км; Кп=8.5
АШ 5 N+E 228 - 0.04 3.3 2.44 0.41 1.59 5.3 0.16 0.78 0.02 0.65 2.86
SEV 5 N+E 287 - 0.04 3.1 4.53 0.43 2.45 8.2 0.26 0.42 0.8 1.85 3.04
SUDU Р ъ 176 0.02 - 4.1 7.29 0.57 1.74 5.8 0.24 0.26 0.61 2.12 3.18
SUDU 5 ^Е 176 - 0.08 3.0 5.57 0.45 2.73 9.1 0.3 0.34 1.03 2.54 3.1
£ 4.6 0.46 2.07 6.92 0.23 0.41 0.63 1.6 3.05
55 0.1 0.03 0.06 0.06 0.06 0.1 0.4 0.13 0.1
Обсуждение результатов. Количественные параметры амплитудного станционного спектра существенно зависят от энергии в очаге, направленности излучения на станцию регистрации, свойств глубинной среды на пути распространения сейсмических волн и под станцией регистрации, а также от ряда факторов, которые трудно учесть. Поэтому для надежных оценок очаговых параметров в настоящей работе использованы не только записи нескольких станций, но и разные типы волн. В осреднении очаговых параметров (5) участвовало от четырех до семи станционных определений, что обеспечило малую величину стандартного отклонения 55 показателя степени рассеяния индивидуальных оценок. Наилучшая сходимость станционных определений получена для радиуса круговой дислокации (5го=0.01-0.05), величина которого (см. формулу 2) в явном виде не зависит от направленности излучения энергии из очага и условий среды под станцией регистрации. Стандартные отклонения 55 по другим параметрам для большинства землетрясений менее <0.1. Уверенно определены средние значения моментных магнитуд Мнуг с погрешностью не выше ±0.14 или округленно ±0.1. В целом, восстановленные значения динамических параметров очагов землетрясений Крымско-Черноморского региона за 2014 г. уверенно относятся к категории надежных.
Наибольшие значения динамических параметров (М0, Г0, Аа, в, ^а , й, Ей иMwрег) получены для очага ощутимого землетрясения 2 марта с Кп=11.2, произошедшего в центральной части региона (район 2) на глубине h=28 км (№ 2), а наименьшие - для его слабого афтершока 6 апреля с Кп=7.2 (№ 8). При этом отметим, что радиационное трение сг для всех исследованных очагов зоны ощутимого землетрясения 2 марта имело отрицательное значение, указывая на сложное
скольжение разрыва в очаге за счет неравномерного распределения прочностных свойств глубинной среды в данной очаговой зоне [16]. Во всем диапазоне энергий величина сброшенных напряжений не превысила Да=1-106 Па (10 бар), а кажущихся напряжений - по<3-105 Па (3 бар).
Проведем сравнение полученных в 2014 г. динамических параметров очагов со средними их долговременными величинами, установленными по аналоговым записям [12, 19] (рис. 4). Сравнение, как и в предыдущие годы [20-22], проведено только для сейсмического момента Мо и радиуса круговой дислокации го, поскольку они являются главными, определяемыми непосредственно по спектрам записей сейсмических колебаний, а остальные параметры вычисляются по их значениям. Для сравнения использованы зависимости М0 (КП) и г0 (КП) из [19].
Рис. 4. Сравнение динамических параметров очагов землетрясений за 2014 г.: сейсмического момента Мо (а) и радиуса круговой дислокации г0 (б) с долговременными зависимостями М0(КП) и г0(Кп) из [19]
1 - экспериментальные значения за 2014 г; 2 - значения для землетрясений с известными механизмами очагов: пунктиром обозначены пределы погрешностей.
Как видно из рис. 4, для большинства землетрясений 2014 г. средние значения М0 и г0 находятся в пределах доверительных областей долговременных зависимостей М0(Кп), Г0(Кп). При этом значения Мо равномерно распределены относительно регрессии М0(Кп), а - в двух случаях имели наибольшие положительные отклонения 5г0 от зависимости г0(КП). Это относится к сравнительно слабым землетрясениям 30 марта в 07ь55т с Кп=8.2, h=21 км (№ 7) из Северо-Западного района (район № 6) и 6 апреля в 20ь47т с Кп=7.2, h=23 км (№ 8) из центральной зоны региона (район № 2). Завышение значений Г0 для сейсмических событий низкого энергетического уровня может быть связано не только с естественными особенностями очагов, но с меньшей точностью определения глубин очагов и, соответственно, выбором скоростной модели среды для этих глубин. Наилучшее соответствие долговременным параметрам М0 (рис. 4) отмечено для землетрясений с известными механизмами очагов, для которых учтена реальная направленность излучения в соответствии с табл. 3.
Энергетические спектры землетрясений. Спектральные свойства очагов землетрясений 2014 г. изучены также по энергетическим спектрам. Переход от спектральной плотности амплитудного спектра к энергетическому спектру осуществлялся по формуле из [23]:
Ч = ~Р— ш3^2(ю), (5)
2л^ е
где р - плотность пород в окрестности очага, с — скорость распространения объемных волн (Р или 5), П(ю) - модуль спектральной плотности из амплитудного спектра. Методика и алгоритм расчета энергетических спектров, их обработка описаны в [12, 23].
На рис. 5 приведены примеры сглаженных энергетических спектров для двух компактных групп землетрясений центральной зоны региона по данным станции «Алушта» с эпицентраль-ными расстояниями Д=(46±3) км и Д=(27±2) км.
Сглаживание проводилось скользящим осреднением в пределах октавной ширины по оси частот. Для всех энергетических спектров (рис. 5) получена характерная резонансная форма с выраженным максимумом плотности энергии Чтах в ограниченной области частот, как и для других землетрясений Крымско-Черноморского региона [12].
N
1А
4х/ \\
гМ1
12 11 /
1 Ш)
-2 -3 -4 -5 -6 -1 -8
-ш у Е
А
1 х2
г* 1 1
П§(/
б
Рис. 5. Сглаженные энергетические спектры двух компактных групп землетрясений по данным станции «Алушта»: а - Д=(46±3) км; б - Д=(27±2) км (номера спектров (1-7) соответствуют таковым в табл. 5)
Рассмотрены основные характеристики энергетических спектров: преобладающие частоты, на которые приходится максимум плотности энергии (^шах) и ширина спектра / на уровне 0.9 от qшax в высокочастотной fq\) и низкочастотной (fq2) части спектра:
/=/4/2. (6)
Результаты расчета fqшax и / для двух групп очагов по записям станции «Алушта»: эпи-центральные расстояния (46±3) км и (27±2) км по продольным волнам (составляющая 2) и поперечным Е) представлены в табл. 5.
Таблица 5. Значения максимума fqшax и ширины Ъfq энергетических спектров землетрясений за 2014 г. по станции «Алушта» для двух групп очагов
№ Дата, Кп Д, км /qшax Г ц / С%-/?2), Гц
дн мес г 2 N Е 2 N Е
Первая группа, Д=(46±3) км
1 02.03.2014 11.2 43 5.69 4.56 2.83 4.35 5.79 4.43
2 02.03.2014 8.5 47 5.69 4.51 2.88 5.55 4.95 4.58
3 06.04.2014 7.2 44 5.71 4.53 2.26 5.42 5.67 4.63
4 27.04.2014 8.5 51 4.51 2.85 2.27 6.62 4.23 4.54
Вторая группа, Д=(27±2) км
5 17.03.2014 9.8 28 5.76 - 2.29 5.58 - 4.39
6 29.03.2014 8.6 28 5.76 - 2.27 5.52 - 3.02
7 20.08.2014 7.9 24 5.69 3.59 2.86 6.74 5.34 5.43
а
Из табл. 5 видно, что, несмотря на вариации индивидуальных энергетических спектров, значения qmax по всем станциям и разным составляющим записи находились в сравнительно ограниченном диапазоне частот:fmax=2.26,...5.76 Гц. При этом для продольных волн (Z) преобладали максимальные значения плотности сейсмической энергии вблизи f=5.7 Гц, а для поперечных - 4.5 Гц (по N) и <3 Гц (по Е). Отмечен ожидаемый сдвиг qmax в сторону более высоких частот для продольных волн (Z-составляющая) по сравнению с поперечными волнами. Прямой зависимости значений fqmax и f от энергетического уровня землетрясения и эпицентрального расстояния из-за ограниченного числа статистических данных не установлено. В то же время ранее по представительной выборке данных (112 значений) по аналоговым записям поперечных S-волн на сейсмической станции «Алушта» для компактной группы с Д=(20±5) км была установлена прямая корреляционная связь между периодом qmax и энергетическим классом [12]:
№тах=(0.023±0.007)К"и - (0.79±0.05), р=0.8.
Расчеты по этой формуле дают для волн S как близкие (в пределах погрешностей), так и несколько отличающиеся значенияfqmax=1/ T qmax, в сравнении с теми, что приведены в табл. 5. Только наращивание статистических данных по энергетическим спектрам за длительный срок наблюдений и расширение динамического диапазона исследуемых землетрясений позволит уточнить ранее полученные зависимости по аналоговым записям [12], исследовать причины индивидуальных отклонений и общие пространственно-временные и станционные свойства для их использования в задачах инженерной сейсмологии.
Литература
1. Козиненко Н.М., Свидлова В.А., Сыкчина З.Н. (отв. сост.). Каталог землетрясений Крымско-Черноморского региона за 2014 г. // Землетрясения Северной Евразии. - Вып. 23 (2014 г.). - Обнинск: фИц ЕГС РАН, 2020. - Приложение на CD-ROM.
2. Пустовитенко Б.Г., Свидлова В.А., Бондарь М.Н. Крымско-Черноморский регион // Землетрясения Северной Евразии. - Вып. 23 (2014 г.). - Обнинск: ФИЦ ЕГС РАН, 2020. - С. 38-48. doi: 10.35540/1818-6254.2020.23.03
3. Пустовитенко Б.Г. (отв. сост.) Каталог механизмов очагов землетрясений Крымско-Черноморского региона за 2014 г. // Землетрясения Северной Евразии. - Вып. 23 (2014 г.). - Обнинск: ФИЦ ЕГС РАН, 2020. - Приложение на CD-ROM.
4. Сейсмологический каталог (сеть телесейсмических станций), 2014 (2020) // ФИЦ ЕГС РАН [сайт]. -URL: ftp://ftp.gsras.ru/pub/Teleseismic_Catalog/2014
5. International Seismological Centre (2020), On-line Bulletin, https://doi.org/10.31905/D808B830
6. DDA http://www.isc.ac.uk/cgi-bin/serve_contribution?repid=600387869
7. Kennet B.L.N. Seismological Tables: ak135 // Research School of earth Sciences Australian national University. Australia, Canberra, ACT0200. - 2005. - P. 80.
8. Кульчицкий В.Е., Сафонова Г.П., Свидлова В.А. Годографы сейсмических волн Крымско-Черноморских землетрясений // Сейсмологический бюллетень Западной территориальной зоны ЕССН СССР (Крым-Карпаты) за 1983 г. - Киев: Наукова думка, 1986. - С. 94-103.
9. Пустовитенко Б.Г., Кульчицкий В.Е. Об энергетической оценке землетрясений Крымско-Черноморского региона // Магнитуда и энергетическая классификация землетрясений. - М.: ИФЗ АН СССР, 1974. - Т. 2. - С. 113-125.
10. Пустовитенко Б.Г, Свидлова В.А Крымско-Черноморский регион // Землетрясения Северной Евразии. - Вып. 21 (2012). - Обнинск: ФИЦ ЕГС РАН, 2018. - С. 49-58.
11. Балакина Л.А., Введенская А.В., Голубева Н.В., Мишарина Л.А., Широкова Е.И. Поле упругих напряжений Земли и механизм очагов землетрясений. - М.: Наука, 1972. - 198 с.
12. Пустовитенко Б.Г., Пантелеева Т.А. Спектральные и очаговые параметры землетрясений Крыма. -Киев: Наукова думка, 1990. - 251 c.
13. Brune J.N. Tectonic stress and the spectra of seismic shear waves from earthquakes // J. Geophys. Res. -1970. - 75. - N 26. - P. 4997-5009.
14. Аптекман Ж.Я., Белавина Ю.Ф., Захарова А.И., Зобин В.М., Коган С.Я., Корчагина О.А., Москвина А.Г., Поликарпова Л.А., Чепкунас Л.С. Спектры Р-волн в задаче определения динамических параметров очагов землетрясений. Переход от станционного спектра к очаговому и расчет динамических параметров очага // Вулканология и сейсмология. - 1989. - № 2. - C. 66-79.
15. Аптекман Ж.Я., Дараган С.К., Долгополов В.В., Захарова А.И., Зобин В.М., Коган С.Я., Корчагина О.А., Москвина А.Г., Поликарпова Л. А., Чепкунас Л.С. Спектры Р-волн в задаче определения динамических параметров очагов землетрясений. Унификация исходных данных и процедуры расчета амплитудных спектров // Вулканология и сейсмология. - 1985. - № 2. - С. 60-70.
16. Костров Б.В. Механика очага тектонического землетрясения. - М.: Наука, 1975. - 179 с.
17. Кейлис-Борок В.И. Исследование источников, приближенно эквивалентных очагам землетрясений // Труды Геофизического Института АН СССР. -1959. - № 9 (136). - C. 20-42.
18. Hanks T.C., Kanamori H. A moment magnitude scale // J. Geophys. Res. - 1979. - 84. - N 35. -P. 2348-2350.
19. Пустовитенко Б.Г., Пустовитенко А.А., Капитанова С.А., Поречнова Е.И. Пространственные особенности очаговых параметров землетрясений Крыма // Сейсмичность Северной Евразии - Обнинск: ГС РАН, 2008. - С. 238-242.
20. Пустовитенко Б.Г., Калинюк И.В., Мержей Е.А. Динамические параметры очагов землетрясений Крымско-Черноморского региона // Землетрясения Северной Евразии, 2010 год. - Обнинск: ГС РАН,
2016. - С. 296-304.
21. Пустовитенко Б.Г., Калинюк И.В., Пустовитенко А.А. Очаговые параметры землетрясений Крымско-Черноморского региона // Землетрясения Северной Евразии в 2011 г.- Обнинск: ФИЦ ЕГС РАН,
2017. - С. 282-293.
22. Пустовитенко Б.Г., Калинюк И.В., Пустовитенко А.А. Очаговые параметры землетрясений Крымско-Черноморского региона // Землетрясения Северной Евразии. - Вып. 21 (2012 г.). - Обнинск: ФИЦ ЕГС РАН, 2018. - С. 286-293.
23. Ризниченко Ю.В., Сейдузова С.С. Спектрально-временная характеристика сейсмической опасности. - М.: Наука, 1984. - 180 с.
SPECTRAL and SOURCE PARAMETERS of CRIMEAN-BLACK SEA EARTHQUAKES
B.G. Pustovitenko, E.E. Eredzhepov
Institute of Seismology and Geodynamics, Crimean Federal University V.I. Vernadsky, Simferopol, [email protected]; [email protected]
Abstract. The spectral and dynamic source parameters (Mo, ro, Act, 8, t|ct, CT r , ü, Eu and Mw) of 13 Crimean earthquakes with Kn = 7.2-11.2 (Mw=2.6-4.1), restored by 123 amplitude spectra of longitudinal and transverse seismic waves recorded by digital regional seismic stations is shown. Approximation of the spectra and the source dynamic parameters calculation based on their parameters is performed in the framework of the Brune dislocation model. Four to seven station definitions participated in the averaging of focal parameters, which ensured a small standard deviation, the scattering degree index of the individual estimates. The best convergence of the station definitions is obtained for the radius of a circular dislocation. The highest values of dynamic parameters have been obtained for the perceptible earthquake on March 2 with Kn=11.2, and the lowest values have been obtained for its weak aftershock on April 6 with Kn = 7.2. Within the total range of energies the value of the stress drop does not exceed Act= 106 Pa (10 bar), and the apparent stress drop does not exceed t|ct <3105 Pa (3 bar). The average values of seismic moments and circular dislocation radius within the errors of their determination match the average long-term correlations of the parameters on the earthquake energy level obtained by the analog recording. For the strongest earthquakes on March 2 with Mwreg=4.1 and October 18 with Mwreg= 3. 8 solutions of the focal mechanisms were obtained. The March 2 earthquake occurred in the central part of the region near the South Coast of Crimea under the action of horizontal tensile forces oriented near the latitude. The type of movement in the focus is a obligue slip with predominance of a normal component over a strike-slip. The earthquake on October 18 occurred in the central part of the Black Sea basin in condition of sublatitudinal compression and submeridional extension. Type of movement in the source is a pure strike-slip. For the main shock on March 2 and its six aftershocks, energy spectra according to the data of Alushta station and their main parameters are given. The analysis of obtained results is given.
Keywords: source mechanism, amplitude spectrum, Brune model, spectral density, the angular frequency, seismic moment, dislocation radius, apparent stress drop, radiation friction, energy spectra.
DOI: 10.35540/1818-6254.2020.23.25
For citation: Pustovitenko, B.G., Eredzhepov E.E. (2020). [Spectral and Source Parameters of Crimean-Black Sea Earthquakes]. Zemletriaseniia SevernoiEvrazii [Earthquakes in Northern Eurasia], 25(2014), 250-262. (In Russ.). doi: 10.35540/1818-6254.2020.23.25
References
1. Kozinenko, N.M.,Svidlova, V.A., & Sykchina, Z.N. (2020). Catalog of earthquakes in the Crimean Black Sea region for 2014. Zemletriaseniia Severnoi Evrazii [Earthquakes in Northern Eurasia], 23(2014), Appendix on CD. (In Russ.).
2. Pustovitenko, B.G., Svidlova, V.A., & Bondar, M.N. (2020). Crimean Black sea region. Zemletriaseniia Severnoi Evrazii [Earthquakes in Northern Eurasia], 23(2014), 38-48. (In Russ.). doi: 10.35540/18186254.2020.23.03
3. Pustovitenko, B.G. (2020). Catalog of earthquake foci mechanisms in the Crimean Black Sea region for 2014. Zemletriaseniia Severnoi Evrazii [Earthquakes in Northern Eurasia], 23(2014), Appendix on CD. (In Russ.).
4. GS RAS, Catalogues of Teleseismic Stations, 2014 (2020). Retrieved from ftp://ftp.gsras.ru/pub/Teleseis-mic_Catalog/2014
5. International Seismological Centre (2020), On-line Bulletin, https://doi.org/10.31905/D808B830
6. DDA. (2020). Retrieved from http://www.isc.ac.uk/cgi-bin/serve_contribution?repid=600387869
7. Kennet, B.L.N. (2005). Seismological Tables: ak135. Research School of Earth Sciences. Australian national University. Australia, Canberra, ACT0200, 80 p.
8. Kul'chickij, V.E., Safonova, G.P., & Svidlova, V.A. (1986). [Hodographs of seismic waves of the Crimean-Black sea earthquakes]. In Seismologicheskii biulleten' Zapadnoi territorial'noi zony ESSN SSSR (Krym-Karpaty) za 1983 g. [Seismological Bulletin of the Western Territorial Zone of the ESSN of the USSR (Crimea-Carpathians) for 1983] (pp. 94-103). Kyiv, Ukraine: Naukova Dumka Publ. (In Russ.).
9. Pustovitenko, B.G., & Kul'chickij, V.E. (1974). [On the energy assessment of the earthquakes of the Crimean Black Sea region] In Magnituda i energeticheskaia klassifikatsiia zemletriasenii, tom 2 [Magnitude and energy classification of earthquakes, Pt. 2] (pp. 113-125). Moskow, Russia: IEP AS USSR Pabl. (In Russ.).
10. Pustovitenko, B.G, & Svidlova, V.A. (2018). [Crimean Black Sea region]. In Zemletriaseniia Severnoi Evrazii v 2012 [Earthquakes in Northern Eurasia, 2012] (pp. 49-58). Obninsk, Russia: GS RAS Publ. (In Russ.).
11. Balakina, L.A., Vvedenskaya, A.V., Golubeva, N.V., Misharina, L.A., & Shirokova, E.I. (1972). Pole upru-gikh napriazhenii Zemli i mekhanizm ochagov zemletriasenii [Field of elastic stresses of the Earth and the mechanism of earthquake foci]. Moscow, Russia: Nauka Publ. 198 p. (In Russ.).
12. Pustovitenko, B.G., & Panteleeva, T.A. (1990). Spektral'nye i ochagovye parametry zemletriasenii Kryma [Spectral and focal parameters of earthquakes of the Crimea]. Kiev, Ukraine: Naukova Dumka Publ., 251 p. (In Russ.).
13. Brune, J.N. (1970). Tectonic stress and the spectra of seismic shear waves from earthquakes. Journal of geophysical research, 75(26), 4997-5009.
14. Aptekman, Zh.Ya., Belavina, Yu.F., Zakharova, A.I., Zobin, V.M., Kogan, S.Ya., Korchagina, O.A., Mos-kvina, A.G., Polikarpova, L.A., & Chepkunas, L. (1989). [Spectra of P-waves in the problem of determining the dynamic parameters of earthquake sources. The transition from station to focal range and the calculation of dynamic source parameters]. Vulkanologiia i seismologiia [Journal of Volcanology and Seismology], 2, 66-79. (In Russ.).
15. Aptekman, Zh.Ya., Daragan, S.K., Dolgopolov, V.V., Zakharova, A.I., Zobin, V.M., Kogan, S.Ya., Korchagina, O.A., Moskvina, A.G., Polikarpova, L.A., & Chepkunas, L.S. (1985). [Spectra of P-waves in the problem of determining the dynamic parameters of earthquake sources. Harmonization of raw data and the procedure of calculating the amplitude spectra]. Vulkanologiia i seismologiia [Journal of Volcanology and Seismology], 2, 60-70. (In Russ.).
16. Kostrov, B.V. (1975). Mekhanika ochaga tektonicheskogo zemletriaseniia [Mechanics of the tectonic earthquake focus]. Moscow, Russia: Nauka Publ., 179 p. (In Russ.).
17. Kejlis-Borok, V. I. (1959). [Study of the sources, approximately equivalent to the foci of earthquakes]. Trudy Geofizicheskogo InstitutaANSSSR. [Proceedings of the Geophysical Institute of the USSR Academy of Sciences], P(136), 20-42 (In Russ.).
18. Hanks, T.C., & Kanamori, H. (1979). A moment magnitude scale. Journal of Geophysical Research, 84, 2348-2350.
19. Pustovitenko, B.G., Pustovitenko, A.A., Kapitanova, S.A., & Porechnova, E.I. (2008). [Spatial features of focal parameters of earthquakes in the Crimea]. In Seismichnost'Severnoi Evrazii [Seismicity of Northern Eurasia] (pp. 238-242). Obninsk, Russia: GS RAS Publ. (In Russ.).
20. Pustovitenko, B.G., Kalinyuk, I.V., & Merzhej, E.A. (2016). [Dynamic parameters of the foci of earthquakes of the Crimea]. In Zemletriaseniia Severnoi Evrazii, 2010 god. [Earthquakes in Northern Eurasia, 2010] Obninsk, Russia: GS RAS Publ., 296-304. (In Russ.).
21. Pustovitenko, B. G.. Kalinyuk, I. V. & Pustovitenko, A.A. (2017). [Focal parameters of earthquakes of the Crimea in 2011]. In Zemletriaseniia Severnoi Evrazii, 2011 god. [Earthquakes in Northern Eurasia, 2011] Obninsk, Russia: GS RAS Publ., 282-293. (In Russ.).
22. Pustovitenko, B.G., Kalinyuk, I.V., & Pustovitenko, A.A. (2018). [Focal parameters of earthquakes of the Crimea in 2012] Zemletriaseniia Severnoi Evrazii [Earthquakes in Northern Eurasia], 21(2012), 286-293. (In Russ.).
23. Riznichenko, Yu.V., & Sejduzova, S.S. (1984) Spektral'no-vremennaia kharakteristika seismicheskoi opas-nosti [Spectral-temporal characteristics of seismic hazard]. Moscow, Russia: Nauka Publ., 180 p. (In Russ.).
