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

УДК 550.4

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

А.В. Викулин, Д.Р. Акманова, Н.А. Осипова

С использованием большого объема данных выявлены периодичности планетарного сейсмического процесса. В результате спектрального Фурье-анализа получены периоды T0 ~ 116 ± 1, 2T0 ~ 195 ± 6, 4T0 ~ 388 ± 4, 8T0 ~ 786 ± 9. «Четность» этого ряда, по-видимому, может быть связана с симметрией процесса - «замкнутостью» сейсмических поясов. Анализ пространственно-временных закономерностей сейсмического процесса окраины Тихого океана показал, что землетрясения с M > 7 имеют тенденцию мигрировать по часовой стрелке в направлении от Новой Зеландии до Южной Америки со скоростью около 300 км/год.

With the help of great scope of data periodicity ofplanetary seismic process has been revealed.

As a result of spectral Fourier’s analysis the periods equal T0 ~ 116 ± 1, 2T0 ~ 195 ± 6, 4T0 ~ 388 ± 4 and 8T0 ~ 786 ± 9 were obtained. «Parity» of these periods can be connected with symmetry of process - «isolation» of seismic belts. The analysis of existential laws of seismic process of suburb of Pacific Ocean has shown, that earthquakes with M > 7 tend to migrate clockwise in a direction from New Zealand to the South America with a speed about 300 km/year.

Первые результаты исследования сейсмичности позволили выявить и продемонстрировать на достаточно большом материале такие свойства сейсмичности, как повторяемость землетрясений и их миграция - закономерное перемещение во времени и пространстве вдоль всего сейсмического пояса некоторой его части или отдельно взятых очагов сильнейших землетрясений (Тамразян, 1962; Рихтер, 1963; Duda, 1963; Mogi, 1968, 1969). Обзор таких данных на начало XXI в. содержится в работах (Викулин, 2001, 2003).

На основании анализа мировых баз данных (Соловьев, Го, 1974, 1975; Новый каталог ..., 1975; Викулин, Ким, 1983; Catalogue ..., 1958, 1966, 1968; Duda, 1965, 1992; Centennial Earthquake Catalog; ISC; JMA; NEIC) был составлен электронный каталог, включающий все известные (датированные) землетрясения (2150 г. до н. э. - август 2008 г.) «до начала инструментального периода» (до XX в.) и с M > 6,0, произошедшие в течение «инструментального периода» (с XX в. до настоящего времени). Сейсмические события в базе характеризуются следующими параметрами: время в очаге (год, месяц, день, час, минута, секунда), координаты очага (долгота и широта в градусах), глубина очага (в км) и магнитуда. Всего в базе приведены сведения о 12 145 землетрясениях (Акманова, Осипова, 2007; Викулин, Водинчар, Мелекесцев и др., 2007).

Повторяемость землетрясений

Известно много примеров, свидетельствующих о повторяемости сильных землетрясений в одном месте через определенные интервалы времени. Для сильнейших (М > 8) тихоокеанских землетрясений выделяется период их повторения в одном месте: Т = 100 ± 50 лет.

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

Таблица 1

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

Фурье-анализ

Мировой, Тихоокеанский и Альпийский каталоги Мировой каталог Тихий океан Альпийский разлом

1 2 3 4

116±1 110-130 100-120 110-130

195±6 Ш-200 Ш-220 Ш-210

250-270

1 2 3 4

388±4 300-400 380-410 350-430

430-470 510-530

786±9 780-800 760-790 770-790

Из данных, представленных в табл. 1, следует, что выявленные периоды расположены в большом диапазоне значений: от 100 до 800 лет. Приведенные цифры подтверждают результаты, полученные ранее в работах (Кириллова, 1957; Федотов, 1965; Филипас, 1965; Shimazaki, Nakata, 1980; Jacob, 1984; Викулин, Викулина, 1989). Как видим, выявленные периоды для выборки 1 (табл. 1) попадают в пределы трех других выборок (2-4).

Между выявленными периодами (в годах) существует очевидная взаимосвязь:

T й 116 ± 1;

2T0 й 195 ± 6;

4T0 й 388 ± 4;

8T0 й 786 ± 9.

Большая часть из выявленных периодов может быть интерпретирована «гармоническим» рядом Т0, 2Т0, 4Т0, 8Т0. «Четность» этого ряда, по-видимому, может быть связана с симметрией сейсмического процесса - «замкнутостью» поясов.

Полученные данные можно рассматривать как подтверждение выявленной структуры периодичности сейсмического процесса.

Миграция очагов землетрясений

Изучение миграционных процессов позволяет выявлять районы и интервал времени (в годах) возникновения потенциальных очагов землетрясений соответствующих магнитуд.

С целью изучения свойств миграции сейсмического процесса и проверки рассчитанных ранее значений скоростей миграции очагов землетрясений исследовались закономерности распределения в пространстве и времени очагов сильнейших (М > 7) землетрясений окраины Тихого океана с использованием разработанного авторами математического алгоритма (Mogi, 1968, 1974; Вилькович, Шнирман, 1982; Кузнецов, Кейлис-Борок, 1997; Акманова, Осипова, 2007; Викулин, Водинчар, Мелекесцев и др., 2007; Акманова, Викулин, Осипова, 2008). На основе данных, представленных в базе, формировались исходные выборки: землетрясения с эпицентрами вдоль окраины Тихого океана с глубинами очага от 0 до 100 км и с магнитудамиM > 7; 7,5; 8; 8,5.

Исследование распределения тихоокеанских землетрясений в пространстве и времени проводилось на плоскости с осями «расстояние по дуге 0 < L < 45 400 км - время извержений t». В качестве расстояния вдоль дуги была выбрана линия, примерно совпадающая с осями глубоководных желобов и границами тектонических плит. Общая протяженность окраины Тихого океана от острова Окленд в Новой Зеландии (Lmin = 0 км) до оконечности Южной Америки составляет Lmax = 45 400 км. Математическое определение этой линии приведено в работах (Викулин, Водинчар, 2005; Викулин, Водинчар, 2006).

Алгоритм построения миграционной цепочки на плоскости L - t в каждой из четырех сформированных совокупностей землетрясений сводился к следующему. Для каждого i-го землетрясения с координатой Li и временем в очаге ti в каталоге искали такое землетрясение i + 1, в очаге которого Li+i > Li, ti+1 > ti. Землетрясение, попавшее хоть один раз в миграционную цепочку, исключалось из дальнейшего анализа. Примеры и параметры таких миграционных цепочек приведены на рис. 1-4 и в табл. 2-3.

— - линия, вдоль которой определялись пространственные координаты эпицентров землетрясений

О, А , П - очаги землетрясений (М > 7; 7,5; 8,5) в миграционной цепочке

------- линия, указывающая направление миграции очагов сильнейших землетрясений

Ь, км 40 000

30 000

20 000

10 000 0

Рис. 1. Положение эпицентров тихоокеанских землетрясений

Ь, км 25 000

20 000

15 000

10 000

5000

Т, год

|*_ 1750 1800 1850 1900 1950 2000

0

1940

Т, год

1960____1980

б

2000 2020

Ь, км

Рис. 2. Пример миграционных цепочек очагов сильнейших землетрясений с М > 7,0

Таблица 2

Параметры выявленных миграционных цепочек

М > 7,0 М > 7,5 М > 8,5

р М Т, год Ь, км М Т, год Ь, км М Т, год Ь, км

1 2 3 4 5 6 7 8 9 10

1 7,4 1763 21 677 7,9 1943 1225 9,0 1687 40 312

2 7,0 1812 32 289 7,5 1943 7846 8,5 1730 42 896

3 7,7 1837 34 400 7,7 1946 10 624 8,7 1730 42 911

4 7,5 1875 34 527 7,5 1957 12 501 8,5 1751 43 333

5 8,0 1882 37 698 7,9 1969 13 621

6 7,6 1906 37 933 7,5 1972 17 732

а

1 2 3 4 5 6 7 8 9 10

7 7,3 1928 39 511 7,5 1977 18 053

8 7,6 1933 41 390 8.0 1978 18 446

9 7,3 1940 41 670 7.7 1999 18 477

10 7,0 1956 42 089

11 7,0 1965 42 100

а, км/год 90 а, км/год 265 а, км/год 50

АТ, год 202 АТ, год 56 АТ, год 64

АЬ, км 20 423 АЬ, км 17 252 АЬ, км 3021

Примечание. Символ р - порядковый номер землетрясения в цепочке; М - магнитуда; Т - год землетрясения; Ь - расстояние вдоль дуги от Ьт1п; а - скорость миграции очагов; АТ - продолжительность цепочки во времени; АЬ -протяженность цепочки в пространстве.

— - линия, вдоль которой определялись пространственные координаты эпицентров землетрясений

[О", И, очаги землетрясений (М > 7; 7,5; 8,5) в миграционной цепочке

-----линия, указывающая направление миграции очагов сильнейших землетрясений

Ь, км 15 000

10 000 500

Рис. 3. Положение эпицентров тихоокеанских землетрясений

40 000 30 000'

20 0001

10 000

б

Ь, км

Т, год

1949 1959 1969 1979 1989 1999 2009

а

Ь, км

Рис. 4. Пример миграционных цепочек очагов сильнейших землетрясений с М > 7,0

Параметры выявленных миграционных цепочек

М > 7,0 М > 7,5 М > 8,0 М > 8,5

п М Т, год Ь, км М Т, год Ь, км М Т, год Ь, км М Т, год Ь, км

1 7,0 1949 4655 8,1 1897 13 853 8,3 1986 4667 9,0 1687 40312

2 7,7 1959 4608 7,9 1897 13 853 8,0 1995 6134 8,5 1730 42896

3 7,5 1961 4679 7,8 1899 19 902 8,2 2000 10783 8,7 1730 42911

4 7,4 1963 5003 7,8 1900 20 261 8,5 1751 43333

5 7,5 1970 8741 7,6 1904 21 982

6 7,8 1980 8929 7,8 1905 24 433

7 7,5 1982 8941 8,0 1906 25 435

8 7,1 1989 9405 7,8 1917 26 823

9 7,5 1996 10 340 7,6 1922 31 342

10 7,3 2000 10 742 8,1 1932 34 491

11 8,2 2000 10 783 7,9 1932 34 552

12 7,0 2001 10 932 7,7 1937 35 193

13 7,2 2004 12 646 7,5 1950 35 219

а, км/год 142 а, км/год 408 а, км/год 404 а, км/год 50

АТ, год 55 АТ, год 53 АТ, год 14 АТ, год 64

АЬ, км 7992 АЬ, км 21 366 АЬ, км 6116 АЬ, км 3021

Примечание. Символ р - порядковый номер землетрясения в цепочке; М - магнитуда; Т - год землетрясения; Ь - расстояние вдоль дуги от Ьтіп; а - скорость миграции очагов; АТ - продолжительность цепочки во времени; АЬ - протяженность цепочки в пространстве.

В результате применения приведенного алгоритма в каталоге тихоокеанских землетрясений (684-2007 гг.) было выявлено достаточно большое количество миграционных цепочек очагов землетрясений.

Приведем количество миграционных цепочек землетрясений (^ разных магнитуд (М), выявленных в каталоге тихоокеанских землетрясений (684-2007 гг.):

М N

М > 7 86

М > 7,5 56

М > 8 25

М > 8,5 9

Расчет средней миграционной цепочки

Средняя миграционная цепочка (СМЦ) - это цепочка со средними параметрами: скоростью

(V, км/год), количеством событий в цепочке (п), продолжительностью цепочки (АТ , год) и ее

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

Таблица 4

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

Параметры СМЦ Магнитудные диапазоны Среднее значение

М > 7 М > 7,5 М > 8 М > 8,5

Скорость миграции V , км/год 203 298 238 546 320 ± 112

Продолжительность миграции АТ , год 162 168 239 256 206 ± 41

Протяженность АЬ, км 18 214 17 716 18 850 13219 17 000 ± 1890

Число событий в цепочке п 19 12 8 4 12 ± 5

Как видим, параметры (V, АТ , АЬ , п) СМЦ не слишком зависят от магнитуд мигрирующих землетрясений. При этом среднюю скорость миграции очагов землетрясений можно принять за V~ 320 км/год.

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

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

появления одной миграционной цепочки с параметрами V, АТ , АЬ , п необходимо рассчитывать путем определения значений соответствующих условных вероятностей. Однако представляется, что наименьшее значение вероятности существования одной СМЦ будет соответствовать случаю, когда все параметры цепочки рассматриваются как независимые.

Расчет вероятностей (табл. 5) осуществлялся с использованием следующих параметров: количество событий (И), протяженность вдоль окраины Тихого океана (Ь, км), продолжительность временного интервала между последним и первым событием в выборке (Т, лет), количество выявленных цепочек (Ис). Тогда вероятность появления одной СМЦ определяется по формуле

Р = Р1 Р2 Рз,

где Р1 - вероятность появления цепочки с п событиями в цепочке из общего числа событий

N каталога:

Р = N

Р2 - вероятность появления цепочки длительностью АТ во всем временном интервале Т каталога:

Р2 =АТТ;

Р3 - вероятность появления цепочки с пространственной протяженностью АЬ в пределах

окраины Тихого океана протяженностью Ь:

Рз =Т.

Вероятность появления Ис таких цепочек рассчитывается как

РМс = РИс .

Таблица 5

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

Параметры Магнитудный диапазон

M > 7 M > 7,5 M > 8 M > 8,5

N 1600 652 187 37

AL 45 300 45 300 43 697 43 697

AT 1359 1359 1359 1133

Nc 86 56 25 9

Примечание. Символ N - количество событий; AL - протяженность цепочки вдоль окраины Тихого океана; AT - продолжительность временного интервала выборки; Nc - количество выявленных цепочек.

Используя данные табл. 5, для вероятностей, рассчитанных с помощью соотношений (1)-(5), получили значения, приведенные в табл. 6. Как следует из табл. 5, значения вероятностей P, Pn, PдТ и PAL различаются в зависимости от магнитуды землетрясения не более чем на порядок. Однако для всех магнитудных диапазонов значение PN остается неизменным: PN ~ 0,1.

Значения, определяющие вероятности существования одной средней миграционной цепочки в случае,

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

- Магнитудный диапазон

М > 7 М > 7,5 М > 8 М > 8,5 Все магнитуды

Рі 0,01 0,02 0,04 0,11 0,01-0,1

Р2 0,12 0,12 0,18 0,23 0,1-0,2

Рз 0,40 0,39 0,43 0,30 0,3-0,4

Р 0,0006 0,0008 0,003 0,007 0,001-0,01

РИ 0,05 0,05 0,08 0,06 = 0,1

Примечание. Символ Р - вероятность появления одной цепочки; р - вероятность появления цепочки с п событиями в цепочке из общего числа событий N каталога; Р2 - вероятность появления цепочки длительностью АТ во всем временном интервале каталога Т; Р3 - вероятность появления цепочки с пространственной протяженностью АЬ в пределах окраины Тихого океана протяженностью Ь; Р^ - вероятность появления N0 таких цепочек.

Как показывает анализ, в случае снятия ограничения (одно значение может только один раз попасть в миграционную цепочку) число выделяемых цепочек увеличивается на порядок (табл. 7). При этом параметры цепочек не слишком отличаются от параметров СМЦ, представленных в табл. 4. Это означает, что вероятность существования совокупности цепочек увеличивается в пределе PN ^ 1. Другими словами, объединение в цепочки очагов землетрясений с магнитудами М > М0 в широком диапазоне магнитуд М0 > 7 является не случайным, а закономерным явлением сейсмического процесса в пределах Тихого океана. При этом средняя скорость миграции очагов землетрясений с М > 7,0 должна составлять около 300 км/год, или 1 см/с, что подтверждает вывод, сделанный в работе (Викулин, 2003).

Таблица 7

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

Магнитудный диапазон N1 N2

М > 7 1546 86

М > 7,5 631 56

М > 8 181 25

М > 8,5 34 9

Всего 2392 176

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

Выводы

1. В результате проведенного Фурье-анализа данных мирового, Тихоокеанского и Альпийского каталогов землетрясений с использованием разработанного программного комплекса было показано существование периодичностей сейсмического процесса (в годах): Т0 «116 ± 1 лет, 2Т0 «195 ± 6 лет, 4Т0 « 388 ± 4 лет, 8Т0 « 786 ± 9 лет.

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

3. Анализ пространственно-временных закономерностей сейсмического процесса окраины Тихого океана показал, что землетрясения с М > 7 имеют тенденцию мигрировать по часовой стрелке в направлении от Новой Зеландии до Южной Америки со скоростью около 300 км/год.

Расчеты сделаны студентами КамчатГТУ во время написания ими курсовых работ и прохождения преддипломной практики в лаборатории геодинамики переходных зон ИВиС ДВО РАН.

Литература

1. Акманова Д.Р. Особенности вулканической активности окраины Тихого океана за последние 12 тыс. лет // Пробл. комплексного геофизического мониторинга Дальнего Востока России: Тез. докл. региональной науч.-техн. конф., Петропавловск-Камчатский, 11-17 ноября 2007 г. -Петропавловск-Камчатский: КФ ГС РАН, 2007. - С. 98.

2. Акманова Д.Р., Осипова Н.А. О взаимосвязи сейсмического и вулканического процессов на примере окраин Тихого океана // Вестник КРАУНЦ. Науки о Земле. - 2007. - № 10. - С. 144-155.

3. Акманова Д.Р., Викулин А.В., Осипова Н.А. Миграция сейсмической и вулканической активности как тектонофизический процесс // Тектонофизика и актуальные вопросы наук о Земле: Тез. докл. Всерос. конф., Москва, 13-17 октября 2008 г. - М.: ИФЗ РАН, 2008. - С. 205-207.

4. Викулин А.В. Миграции и осцилляции сейсмической активности и волновые движения земной коры // Пробл. геодинамики и прогноза землятресений. - Хабаровск: ИТиГ ДВО РАН, 2001. -С. 205-224.

5. Викулин А.В. Физика волнового сейсмического процесса. - Петропавловск-Камчатский: КГПУ, 2003. - 151 с.

6. Викулин А.В., Викулина С.А. Закономерности размещения очаговых областей сильнейших землетрясений в районе желоба Нанкай. - Петропавловск-Камчатский: КГС ИФЗ АН СССР, 1989. - Препринт № 5.

7. Викулин А.В., Водинчар Г.М. Спектр энергии потока сейсмичности // Материалы ежегод. конф., посвященной Дню вулканолога. - Петропавловск-Камчаткий: ИВиС ДВО РАН, 2005. -С. 167-174.

8. Викулин А.В., Водинчар Г.М. Волны миграции сейсмической энергии // Тектоника, глубинное строение и минерагения востока Азии: Тез. докл. V Косыгинских чтений, Хабаровск, 24-27 января 2006. - Хабаровск, 2006. - С. 206-209.

9. Моделирование геодинамических процессов окраины Тихого океана / А.В. Викулин, Г.М. Водинчар, И.В. Мелекесцев и др. // Солнечно-земные связи и предвестники землетрясений: Сб. докл. IV междунар. конф., Петропавловск-Камчатский, 14-17 августа 2007 г. - Петропав-ловск-Камчатский: ИКИР ДВО РАН, 2007. - С. 275-280.

10. Викулин А.В., Ким Ч.У. Курило-Камчатские землетрясения: данные наблюдений. 1911-1952 гг. - М.: МЦДБ, 1983. - 82 с.

11. Викулин А.В., Ким Ч.У. Курило-Камчатские землетрясения: данные наблюдений. 1953-1961 гг. - М.: МЦДБ, 1983. - 84 с.

12. Вилькович Е.В., Шнирман М.Г. Математические модели строения Земли и прогноз землетрясений / Вычислительная сейсмология: Сб. науч. ст. - М.: Наука, 1982. - С. 27-37.

13. Кириллова И.В. О периодичности разрушительных землетрясений Кавказа и Турции // Докл. АН СССР. - 1957. - Т. 115. - № 4. - С. 771-773.

14. Кузнецов И.В., Кейлис-Борок В.И. Взаимосвязь землетрясений Тихоокеанского сейсмического пояса // Докл. РАН. - 1997. - Т. 355. - № 3. - С. 389-393.

15. Новый каталог сильных землетрясений на территории СССР с древнейших времен до 1975 г. - М.: Наука, 1977. - 536 с.

16. Рихтер Ч. Элементарная сейсмология. - М.: Изд-во иностр. лит., 1963. - 670 с.

17. Соловьёв С.Л., Го Ч.Н. Каталог цунами на западном побережье Тихого океана. - М.: Наука, 1974. - 309 с.

18. Соловьёв С.Л., Го Ч.Н. Каталог цунами на восточном побережье Тихого океана. - М.: Наука, 1975. - 203 с.

19. Тамразян Г.П. О периодичности сейсмической активности в течение последних полутора двух тысяч лет (на примере Армении) // Изв. АН СССР. Сер. «Геофизика». - 1962. - № 1. -С. 76-85.

20. Федотов С.А. О закономерностях распределения сильных землетрясений Камчатки, Курильских островов и Северо-Восточной Японии // Тр. ИФЗ АН СССР. - 1965. - № 203/36. -С.66-93.

21. Филипас С. Ф. О некоторых закономерностях сейсмической активности северо-западной части Тихоокеанского пояса // Современные движения земной коры. - 1965. - № 2. - С. 366-375.

22. Catalogue of Major Earthquakes which occurred in and near Japan (1926-1956) // Seismol. Bull. Japan Meteorolog. Agency. Suppl. 1. - Tokyo, 1958. - 91 p.

23. Catalogue of Major Earthquakes which occurred in and near Japan (1957-1962) // Seismol. Bull. Japan Meteorolog. Agency. Suppl. 2. - Tokyo, 1966. - 47 p.

24. Catalogue of Major Earthquakes which occurred in and near Japan (1963-1967) // Seismol. Bull. Japan Meteorolog. Agency. Suppl. 3. - Tokyo, 1968. - 61 p.

25. Centennial Earthquake Catalog: (http://earthquake.usgs.gov/research/data/centennial.php).

26. Duda S.J. Strain release in the Circum-Pacific belt, Chile 1960 // J. Geophys. Res. - 1963. -№ 68. - P. 5531-5544.

27. Duda S.J. Secular seismic energy release circum-pacific belt // Tectonophysics. - 1965. -№ 2(5). - P. 409-452.

28. Duda S.J. Global earthquakes 1903-1985. - Hamburg F.R.Germany: NEIC, 1992. - 183 p.

29. ISC: (http://www.isc.ac.uk/Bulletin/rectang.htm).

30. Jacob K.H. Estimates of long-term probabilities for future great earthquakes in the Aleutians // Geophys. Res. Lett. - 1984. - № 4. - P. 295-298.

31. JMA: (http://www.jma.go.jp/jma/index.html).

32.Mogi K. Migration of Seismic Activity // Bulletin of the Earthquake Research Institute. - 1968. - Vol. 46. - P. 53-74.

33.Mogi K. Some features of recent seismic activity in and near Japan. Activity before and after great earthquakes // Bull. of the Eathquake Res. Inst. - 1969. - Vol. 47. - P. 395-417.

34.Mogi K. Active periods in the world’s shief belts // Tectonophysics. - 1974. - Vol. 22. -№ 3-4.- P. 265-282.

35. NEIC: (http://neic.usgs.gov/neis/epic/epicglobal.html).

36. Shimazaki K., Nakata T. Time-predictable recurrence model for large earthquakes // Geophys. Res. Lett. - 1980. - № 4. - P. 279-282.