CC BY
35
УДК 550.34
ПЕРИОДИЧНОСТЬ КАТАСТРОФИЧЕСКИХ ИЗВЕРЖЕНИЙ И ИХ МИГРАЦИЯ ВДОЛЬ ОКРАИНЫ ТИХОГО ОКЕАНА
А.В. Викулин, Д.Р. Акманова, Н.А. Осипова (Институт вулканологии и сейсмологии ДВО РАН), С.В. Чебанюк, Т.В. Залипина, Ю.О. Фирсова, И.А. Ященко (КамчатГТУ)
С помощью созданной базы данных извержений вулканов мира выявлены периодичности вулканического процесса как в мировом масштабе, так и в локальном: для отдельных районов окраины Тихого океана и отдельно взятого вулкана. Оказалось, что выявленные периоды близки периодам сейсмического процесса (в годах): T1 ~ 198 ± 17, T2 ~ 2T1 ~ 376 ± 12 и Т3 ~ 4Т1 ~ 762 ± 17. « Четность » этих периодов мо^жет свидетельствовать о «замкнутости» окраины Тихого океана. Кроме того, изучалось пространственно-временное распределение вулканической активности в пределах окраины Тихого океана. Показано, что вулканические извержения с W > 4 (объем изверженного материала составляет 0,01 км3 и более) за последние 12 тыс. лет имеют тенденцию мигрировать по часовой стрелке в направлении от Новой Зеландии до Южной Америки со скоростью около 10 км/год, и это подтверждает ранее полученные результаты.
By means of world’s volcanoes eruptions database periodicity of volcanic process in separate areas of suburb of Pacific Ocean and for separately taken volcano is revealed. It has appeared that the revealed periods are close to the periods of seismic process: T1 ~ 198 ± 17, T2 ~ 2T1 ~ 376 ± 12 and T3 ~ 4T1 ~ 762 ± 17 years. ”Parity” of these periods can testify to ”isolation” of suburb of Pacific ocean. Existential distribution of volcanic activity within suburb of Pacific Ocean was studied.
It has been shown, that volcanic eruptions with W > 4 (the volume of erupted is material 0.01 км3 and more) for last 12 thousand years tend to migrate clockwise in the direction from New Zealand to the South America with a speed about 10 km/year, confirming earlier received results.
Сейсмичность и вулканизм как процессы планетарного масштаба взаимосвязаны. Большая часть всех землетрясений, в том числе и самые сильные, располагаются в пределах тектонически активных поясов. Вблизи таких поясов расположены и самые активные вулканы мира. Согласно обзору данных, приведенных в работах (Акманова, 2007; Викулин, 2003; Викулин, Во-динчар, Мелекесцев и др., 2007), отсутствуют данные, однозначно указывающие на существование взаимосвязи этих процессов. Тем не менее ряд исследователей склоняется к мысли о существовании такой взаимосвязи.
Ранее авторами (Викулин, 2003; Акманова, 2007; Викулин, Водинчар, Мелекесцев и др., 2007) была сформулирована гипотеза о том, что взаимосвязь имеет не локальный, а региональный характер, т. е. взаимосвязаны не конкретные землетрясения и извержения вулканов, а их совокупности.
Обзор опубликованных данных по миграции вулканической активности показывает, что миграция была выявлена как в пределах отдельно взятых вулканических центров (Kenneth et al, 1986; Clague, Dalrymple, 1987; Lonsdale, 1988; Леонов, 1991), так и в пределах вулканических дуг (Berg, Sutton, 1974; Sauers, 1986; Викулин, 2003). Скорости миграции вулканической активности в пределах отдельно взятых вулканических центров расположены в пределах V = 3 • 10-5 2 • 10-3 км/год; в пре-
делах вулканических дуг V = 100 и 900 км/год. В работе (Викулин, 2003) установлено на региональном уровне существование корреляции между извержениями вулканов и сильнейшими землетрясениями. На этом основании был сделан вывод о миграции вулканической активности.
Данные, полученные ранее авторами, позволили сделать предположение о существовании миграции вулканической активности в пределах окраины Тихого океана, являющимся одним из самых активных регионов планеты (Акманова, 2007; Акманова, Осипова, 2007; Викулин, Во-динчар, Мелекесцев и др., 2007).
С целью изучения свойств миграции вулканического процесса исследовались закономерности распределения в пространстве и времени катастрофических извержений (Ж > 6) окраины Тихого океана, а также велась проверка рассчитанных ранее (Акманова, 2007; Викулин, Водин-чар, Мелекесцев и др., 2007; Акманова, Викулин, Осипова, 2008) значений скоростей миграции. С этой целью был создан электронный каталог, включающий все известные (датированные) данные о вулканических извержениях планеты: временной период Т = 9850 г. до н. э. - 2007 гг., число извергавшихся вулканов N = 630, количество их извержений п = 6415 (Акманова, 2007; Акманова, Осипова, 2007; Викулин, Водинчар, Мелекесцев и др., 2007). В качестве количественной характеристики извержений была использована шкала 1 < Ж < 7, отражающая в логарифмическом масштабе объем извергнутого материала. При этом значениям Ж соответствуют объемы материала, равные 105, 106, ..., 1011 м3.
Периодичность вулканического процесса
Развитие любого глобального процесса Земли происходит неравномерно, и в проявлении его максимумов выделяется периодичность. И.И. Гущенко отмечал наличие определенной цикличности в извержениях магматических вулканов (И.И. Гущенко, 1985). Он выделяет циклы со следующими периодами: 1 год, 5-6 лет, 23 года, 60-90 лет, 180 лет. В работе (Белов, 1986) убедительно доказано, что в истории Земли в течение последних 600 млн лет в глобальном проявлении вулканизма имели место периодические пульсационные колебания вулканической активности, происходившие на фоне общего увеличения масштабов вулканизма. По данным количественных расчетов, в глобальном проявлении древнего вулканизма Земли также устанавливается четкая периодичность, при которой максимумы проявления вулканизма отмечаются через каждые 190-200 млн лет.
Анализ данных, характеризующих частоту вулканических извержений за последние 400 лет, показал экспоненциальное нарастание числа событий от года к году. Однако это связано не с природным нарастанием активности процесса, а с возрастанием числа сообщений о них по мере расширения границ географических знаний и количества наблюдателей в средневековье. На фоне этого роста отличаются локальные временные всплески в 1630-1700, 1730-1810, 1840-1910, 1920-1930 и 1950-1960 гг., которые возможно, но не очевидно отражают активизацию процесса как такового. Наиболее вероятные гармоники, характеризующие частоту вулканических извержений, имеют продолжительность 20-30, 50, 80, 110-130, 160 и около 200 лет. Наиболее характерные гармоники вариаций энергии вулканических извержений имеют продолжительность 30, 70 и 100 лет (Леви и др., 2003). В работе (Хаин, Халилов, 2008) приводятся данные о наличии 1-2-, 11-, 22-, 50-, 60-, 80-летних циклов.
Методом Фурье-анализа выявлялись периоды вулканического процесса не только в мировом, но и в локальном масштабе, т. е. периоды для вулканов, расположенных в пределах окраины Тихого океана и его отдельно взятых регионов, в том числе восточной и западной окраин. Анализ показал наличие большого количества выявленных периодов, заключенных в широком диапазоне - от 170 до 2700 лет. Значения выявленных периодов оказались различными для каждой выборки. Все выборки содержат периоды, продолжительности которых заключены в следующих диапазонах:
Т к 198 ± 17 лет;
Т2 к 2Т ^ 376 ± 12 лет;
Т3 к 4Т к 762 ± 17 лет.
Кроме этого, методом Фурье-анализа выявлялись периоды для отдельно взятых вулканов мира N = 9), число извержений которых составляет более 70. Анализ показал наличие только четных периодов: 2Т0 к 210 лет, 4Т0 к 390 лет, 6Т0 к 610 лет, 8Т0 к 830 лет, 10Т0 к 1010 лет,
12Т0 к 1 230 лет, где Т0 к 110 ± 10 лет, что может быть связано с «замкнутостью» окраины
Тихого океана.
Полученные периоды наблюдаются в извержениях как одного вулкана, так и нескольких, расположенных в пределах одного региона. Так, извержения вулкана Кракатау (Индонезия) в 416 и 1883 гг. произошли с интервалом в 1467 года, извержения вулкана Бона (Аляска) в 110 и 740 гг. -с периодом в 630 лет. Вулкан Чурчил (Аляска) извергался в 60 и 740 гг. с периодом 640 лет, вулканы Амбрим (Новые Гибриды) в 50 г. и Кувае (Новые Гибриды) в 1452 г. - с периодом в 1402 года, вулканы Кагуяк (Аляска) в 415 г. и Новарупта (Аляска) в 1912 г. - с периодом в 1497 лет. Интервал между извержениями вулканов Ксудач и Опала (Камчатка) составляет 370 лет.
Миграция вулканической активности
Исследование пространственно-временного распределения извержений тихоокеанских вулканов будет проводиться на плоскости с осями «расстояние по дуге 0 < Ь < 45 400 км - время извержений ^». В качестве расстояния вдоль дуги была выбрана линия, примерно совпадающая с осями глубоководных желобов и границами тектонических плит. Определенная таким образом общая протяженность окраины Тихого океана от вулкана Бакл Айленд (Антарктика) до Десепсьон (Южные Шетландские острова) составила 45 400 км (Викулин, Водинчар, 2005; Викулин, Водинчар, 2006).
Из данных об извержениях вулканов окраины Тихого океана, приведенных в каталоге, были сформированы три выборки: извержения с объемом выброшенного материала Ж > 4, 5, 6. Алгоритм выявления и построения миграционной цепочки в каждой из трех сформированных совокупностей извержений вулканов сводился к следующему. Для каждого /-го извержения вулкана с координатой Ь/ и временем ^ в каталоге искалось такое извержение вулкана I + 1, координата и время которого Ь/ + 1 > Ь/, ^ + 1 > и. Извержение вулкана, попавшее хотя бы один раз в миграционную цепочку, исключалось из дальнейшего анализа. В результате применения описанного алгоритма в каталоге тихоокеанских извержений вулканов в период 9850 г. до н. э. - 2007 гг. было выявлено достаточно большое количество миграционных цепочек (N) вулканических извержений с различным объемом выброшенного материала (Ж):
Ж N
Ж > 4 36
Ж > 5 20
Ж > 6 11
Всего 67
Вследствие этого были выявлены все возможные миграционные цепочки вулканических извержений и определены их параметры: количество событий в каждой цепочке, их продолжительность и протяженность вдоль окраины Тихого океана, скорость миграции в цепочке. Примеры выявленных цепочек и их параметры представлены на рис. 1-4 и в табл. 1, 2.
— - линия, вдоль которой определялись пространственные координаты эпицентров землетрясений
О > А , - извержения с Ж > 4,
5, 6 в миграционной цепочке
-------линия, указывающая направление миграции вулканических извержений
Рис. 1. Расположение вулканов окраины Тихого океана (цифры обозначают последовательность вулканических извержений в разных цепочках событий)
Ь, км
Рис. 2. Пример миграционных цепочек очагов сильнейших землетрясений с Ж > 6 (а), Ж > 5 (б), Ж > 4 (
— - линия, вдоль которой определялись пространственные координаты вулканических извержений
о, А . - извержения с Ж> 4,
5, 6 в миграционной цепочке
-------- линия, указывающая направление миграции вулканических извержений
Рис. 3. Расположение вулканов окраины Тихого океана
Параметры выявленных миграционных цепочек
Таблица 1
Ж > 6 Ж > 5 Ж > 4
р Ж Т, год Ь, км Ж Т, год Ь, км Ж Т, год Ь, км
1 2 3 4 5 6 7 8 9 10
1 6 -250 4490 5 1310 3386 6 -7460 17 821
2 6 50 8467 6 1452 8391 4 -7450 22 019
3 6 540 10 769 6 1580 10 362 5 -7050 23 386
4 6 710 10 906 5 1641 13 919 6 -6600 23 720
5 6 800 10 967 7 1815 14 451 5 -6400 23 848
6 6 1660 11 311 5 1963 14 548 5 -5400 23 883
а
1 2 3 4 5 6 7 8 9 10
7 6 1883 15 514 6 1991 17 821 5 -2550 36 580
8 4 -2250 36 846
9 4 -1050 37 861
10 4 -850 37 870
11 5 -690 38 320
12 4 -270 39 111
13 4 1993 41 787
а, км/год 4 а, км/год 18 а, км/год 3
АТ, год 2133 АТ, год 681 АТ, год 9453
АЬ, км 11 024 АЬ, км 14 435 АЬ, км 23 966
Примечание. Символ р - порядковый номер землетрясения в цепочке; Ж - объем изверженного материала; Ь - расстояние вдоль дуги от Ьш;п; а - скорость миграции; АТ - продолжительность цепочки во времени; АЬ - протяженность цепочки в пространстве.
Ь, км
Рис. 4. Пример миграционных цепочек очагов сильнейших извержений с Ж > 6 (а), Ж > 5 (б), Ж > 4 (в)
Параметры выявленных миграционных цепочек
Таблица 2
Ж > 6 Ж > 5 Ж > 4
р Ж Т, год Ь, км Ж Т, год Ь, км Ж Т, год Ь, км
1 2 3 4 5 6 7 8 9 10
1 6 -7460 17 821 5 1310 3386 5 -6150 3386
2 6 -7350 20 142 6 1452 8391 5 -5150 10 776
3 6 -6600 22 019 6 1580 10 362 7 -4350 19 750
4 6 -6600 23 720 5 1641 13 919 4 -4350 19 889
5 6 -5700 23 786 7 1815 14 451 4 -3790 23 620
6 7 -4900 31 149 5 1963 14 548 6 -1645 27 062
7 6 -1890 44 364 6 1991 17 821 5 -800 27 119
8 4 -570 31 727
9 4 -170 36 613
10 4 170 36 846
11 4 1330 37 861
12 4 1350 37 870
13 4 1570 38 274
14 4 1660 38 348
а
1 2 3 4 5 6 7 8 9 10
15 4 1744 38 376
16 4 1916 39 111
а, км/год 5 а, км/год 18 а, км/год 4
АТ, год 5570 АТ, год 681 АТ, год 8066
АЬ, км 26543 АЬ, км 14435 АЬ, км 35 725
Примечание. Символ р - порядковый номер землетрясения в цепочке; Ж - объем изверженного материала; Ь - расстояние вдоль дуги от Ьш;п; а - скорость миграции; АТ - продолжительность цепочки во времени; АЬ - протяженность цепочки в пространстве.
Для каждой из выявленных 67 цепочек определялось значение скорости миграции методом наименьших квадратов. Диапазон скоростей «глобальной» миграции извержений вулканов с различными значениями Ж (Ж > 4, Ж > 5 и Ж > 6) оказался расположенным в пределах 2 < V < 96 км/год (рис. 3-5). Как видим, полученные в работе (Акманова, Викулин, Осипова, 2008) данные о скоростях миграции извержений вулканов с Ж > 4 не противоречат ранее полученным результатам.
Гистограммы, представленные на рис. 5, указывают на наличие некоторых закономерностей в распределении скоростей миграции катастрофических вулканических извержений (Ж > 4, 5, 6).
V, км/год
а
п
V, км/год
б
п
V, км/год
в
Рис. 5. Графики зависимости миграционных цепочек (п) от скорости (V) для извержений вулканов с Ж > 6 (а), Ж > 5 (б), Ж > 4 (в)
Общей закономерностью для всех значений Ж в первом (грубом) приближении можно считать наличие «широкого высокоамплитудного» максимума в диапазоне 2 < V < 24 (км/год) и нескольких «узких малоамплитудных» максимумов, соответствующих более высоким значениям скоростей 31 < V < 96 км/год. Такое распределение скоростей миграции позволяет предположить существование общих закономерностей для всего класса катастрофических извержений в диапазоне Ж > 4, 5, 6. Все «локальные» максимумы скоростей миграции, выявленные на гистограммах (см. рис. 5), приведены в табл. 3.
Таблица 3
«Локальные» максимумы скоростей миграции вулканических извержений
Ж > 6 Ж > 5 Ж > 4
2 2
3
4 4 4
6
7 7
9 9
12
13
15
16
20
22
24
31
41
48
54
69
86
96
Согласно табл. 3 в диапазоне значений 2 < V < 22 (км/год) скорости для всех анализируемых диапазонов Ж имеют тенденцию группироваться в пределах достаточно «узких» диапазонов (ограниченных жирными линиями): 3 ± 1, 8 ± 1, 19 ± 3 км/год. Именно в эти узкие диапазоны, отстоящие друг от друга на AV = 16 км/год, попадают 53 (79%) значения скоростей из 67.
Расчет средней миграционной цепочки
Средняя миграционная цепочка (СМЦ) - это цепочка со средними параметрами: скоростью (V , км/год), количеством событий в цепочке (п), продолжительностью цепочки (АТ, год), протяженностью цепочки вдоль дуги (АЬ, км). Приведенные выше данные позволили рассчитать параметры СМЦ (табл. 4) для различных значений Ж (объем выброшенного материала).
Таблица 4
Параметры средних миграционных цепочек для различных объемов выброшенного материала
Параметры СМЦ Объем выброшенного материала
Ж > 4 Ж > 5 Ж > 6
V 17 11 3
АТ 4449 4808 4684
А1 19 784 19 731 13 496
п 13 8 5
Как следует из табл. 4, параметры (V, АТ , АЬ , п) средней цепочки не слишком зависят от объема выброшенного материала мигрирующих извержений. При этом среднюю скорость миграции можно принять за V ~ 10 км/год.
Вероятность появления средней миграционной цепочки
Известно, что вулканические извержения имеют тенденцию перемещаться в пространстве и времени. Поэтому определение вероятности появления одной миграционной цепочки с параметрами V, АТ , АЬ , п (см. табл. 4) необходимо рассчитывать путем определения значений соответствующих условных вероятностей. Однако представляется, что наименьшее значение вероятности существования одной СМЦ будет соответствовать случаю, когда все (V, АТ , АЬ , п) параметры цепочки рассматриваются как независимые.
Если каждая выборка (Ж > 4, 5, 6) характеризуется своим набором параметров (табл. 5): количеством событий (п), протяженностью вдоль окраины Тихого океана (АЬ, км), продолжительностью временного интервала между последним и первым событием в выборке (АТ, лет), количеством выявленных цепочек (Ас), - тогда вероятность появления одной СМЦ определяется по формуле
Р = Р1 Р2 Рз, (1)
где Р1 - вероятность появления цепочки с п событиями в цепочке из общего числа событий N каталога:
Р=N; (2)
Р2 - вероятность появления цепочки длительностью АТ во всем временном интервале каталога Т:
Р2 =АТ; (3)
Р3 - вероятность появления цепочки с пространственной протяженностью АЬ в пределах окраины Тихого океана протяженностью Ь:
Рз =Т. (4)
Вероятность появления Ас таких цепочек рассчитывается как
Рас = РАС. (5)
Таблица 5
Параметры Объем выброшенного материала
Ж > 4 Ж > 5 Ж > 6
п 460 162 58
АЬ 41 913 41 913 41 061
АТ 11 857 11 841 9471
Ас 36 20 11
Примечание. Символ п - количество событий в цепочке; АЬ - протяженность цепочки вдоль окраины Тихого океана; АТ - продолжительность временного интервала выборки; Ас - количество выявленных цепочек.
С помощью данных табл. 5 были получены значения для вероятностей, рассчитанных с применением соотношений (2)-(6) (табл. 6).
Как можно видеть из данных, приведенных в табл. 5, значения вероятностей Р, Р1, Р2, Р3 различаются в зависимости от объема выброшенного материала не более чем на порядок. Однако для всех Ж значение РАс остается равным: 0,2 < РАс < 0,66. Как показывает анализ, в случае снятия ограничения (одно значение может попасть в миграционную цепочку только один раз) число выделяемых цепочек увеличивается на порядок (см. табл. 6). Это означает, что РАс ^1.
Другими словами, объединение вулканических извержений с объемом выброшенного материала Ж > 4 является не случайным, а закономерным явлением вулканического процесса в пределах Тихого океана. При этом средняя скорость миграции извержений вулканов с Ж > 4 должна составлять около 10 км/год, или к 0,3 см/с, что подтверждает вывод, сделанный в работе (Викулин, 2003).
Таблица 6
Значения, определяющие вероятности существования одной средней миграционной цепочки в случае, когда ее параметры считаются независимыми
Значения вероятностей Объем выброшенного материала
Ж > 4 Ж > 5 Ж > 6 Диапазон Ж
Р1 0,03 0,05 0,4 0,03-0,4
Р2 0,4 0,4 0,5 0,4-0,5
Р3 0,5 0,5 0,3 0,3-0,5
Р 0,006 0,01 0,06 0,006-0,01
РУс 0,2 0,2 0,66 0,2-0,66
Примечание. Символ Р - вероятность появления одной цепочки; Р1 - вероятность появления цепочки с п событиями в цепочке из общего числа событий N каталога; Р2 - вероятность появления цепочки длительностью АТ во всем временном интервале каталога Т; Р3 - вероятность появления цепочки с пространственной протяженностью АЬ в пределах окраины Тихого океана протяженностью Ь; РУс - вероятность появления N0 таких цепочек.
В случае, когда одно извержение вулкана может быть неоднократно составной частью миграционных цепочек, число выявленных цепочек миграции извержений вулканов (/У) с разным объемом выброшенного материала (Ж) представляется в виде следующего соотношения:
Ж N
Ж > 4 454
Ж>5 160
Ж > 6 55
Всего 669
Выводы
1. Анализ пространственно-временных закономерностей вулканического процесса окраины Тихого океана показал, что вулканические извержения с Ж > 4 за последние 12 тыс. лет имеют тенденцию мигрировать по часовой стрелке в направлении от Новой Зеландии до Южной Америки со скоростью около 10 км/год, что подтверждает полученные ранее результаты (Акманова, Викулин, Осипова, 2008).
2. В результате проведенного анализа данных о вулканических извержениях мирового каталога, окраины Тихого океана и его регионов с использованием разработанного программного комплекса было доказано существование периодичностей вулканического процесса (в годах): 2Т0 к 210, 4Т0 к 390, 6Т0 к 610, 8Т0 к 830, 10Т я 1010, 12Т я 1 23, где Т я 110 ± 10 лет.
Расчеты сделаны студентами КамчатГТУ во время написания ими курсовых работ и прохождения преддипломной практики в лаборатории геодинамики переходных зон ИВиС ДВО РАН.
Литература
1. Акманова Д.Р. Особенности вулканической активности окраины Тихого океана за последние 12 тыс. лет // Пробл. комплексного геофизического мониторинга Дальнего Востока России: Тез. докл. I региональной науч.-техн. конф., Петропавловск-Камчатский, 11-17 ноября 2007 г. - Петропавловск-Камчатский: КФ ГС РАН, 2007. - С. 98.
2. Акманова Д.Р., Осипова Н.А. О взаимосвязи сейсмического и вулканического процессов на примере окраин Тихого океана // Вестник КРАУНЦ. Науки о Земле. - 2007. - № 10. - С. 144-155.
3. Акманова Д.Р., Викулин А.В., Осипова Н.А. Миграция сейсмической и вулканической активности как тектонофизический процесс // Тектонофизика и актуальные вопр. наук о Земле: Тез. докл. Всерос. конф., Москва, 13-17 октября 2008 г. - М.: ИФЗ, 2008. - С. 205-207.
4. Белов С.В. О периодичности современного и древнего вулканизма Земли // ДАН СССР. -1986. - Т. 291. - № 2. - С. 421-425.
5. Викулин А.В. Физика волнового сейсмического процесса. - Петропавловск-Камчатский: КГПУ, 2003. - 151 с.
6. Викулин А.В., Водинчар Г.М. Спектр энергии потока сейсмичности // Материалы ежегод. конф., посвященной Дню вулканолога. - Петропавловск-Камчатский: ИВиС ДВО РАН, 2005. -С.167-174.
7. Викулин А.В., Водинчар Г.М. Волны миграции сейсмической энергии // Тектоника, глубинное строение и минерагения востока Азии: Материалы V Косыгинских чтений, Хабаровск, 24-27 января 2006 г. - Хабаровск, 2006. - С. 206-209.
8. Моделирование геодинамических процессов окраины Тихого океана / А.В. Викулин, Г.М. Водинчар, И.В. Мелекесцев и др. // Солнечно-земные связи и предвестники землетрясений: Сб. докл. IV междунар. конф., Петропавловск-Камчатский, 14-17 августа 2007 г. - Петропав-ловск-Камчатский: ИКИР ДВО РАН. - 2007. - С. 275-280.
9. Голицын Г.С. Объяснение зависимости частота - объем извержений вулканов // ДАН. -2003. - Т. 390. - № 3. - С. 394-396.
10. Гущенко И.И. Цикличность извержений вулканов мира // Вулканология и сейсмология. -1985.- № 2. - С. 31-43.
11. Современная геодинамика и гелиогеодинамика: 500-летняя хронология аномальных явлений в природе и социуме Сибири и Монголии: Учеб. пособие для вузов. Кн. 2 / К.Г. Леви, Н.В. Задонина, Н.Е. Бердникова и др. - Иркутск: ИрГТУ, 2003. - 383 с.
12. Леонов В.Л. О некоторых закономерностях развития гидротермальной и вулканической деятельности на Камчатке // Вулканология и сейсмология. - 1991. - № 2. - С. 28-40.
13. Хаин В.Е., Халилов Э.Н. Пространственно-временные закономерности сейсмической и вулканической активности. - Бургас: SWB, 2008. - 304 с.
14. Berg E., Sutton G.H. Dynamic interaction of seismic and volcanic activity of the Nazca plate edges // Phys. of the Earth and Plan. Inter. - 1974. - № 9. - P. 45-68.
15. Clague D.A., Dalrymple G.B. The Hawaiian - Emperor volcanic chain. P. I. Geologic Evolution // Volcanism in Hawaii. Chap. I. U.S. Geol. Survey Prof. Paper 1350. - Washington, 1987. - P. 5-54.
16. Gresta S., Marzocchi W., Mulargia F. Is there a correlation between larger local earthquakes and the end eruptions at mount Etna volcano, Sicili // Geophys. J. Int. - 1994. - № 1. - P. 230-232.
17. Migration of volcanism in the San Francisco volcanic field, Arizona / L. Kenneth, E.M. Tanaka, G. Shoemaker et all // GSA Bulletin. February. - 1986. - Vol. 97. - № 2. - P. 129-141.
18. Lonsdale P. Geography and history of the Louisville hot spot chain in the southwest Pacific // J. Geophys. Res. - 1988. - Vol. 93. - № 34. - P. 3078-3104.
19. Sauers J. The westward migration of geophysical events in the Aleutians, Springs, 1986 // Cycles. - 1986. - № 9. - P. 203-204.
