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

В.И. Горельчик, А.В. Сторчеус

В нижних горизонтах земной коры и переходном от коры к мантии слое (интервал глубин 20-40 км) под Ключевским вулканом выявлена зона с аномальными физическими свойствами, имеющая эффективную вязкость на несколько порядков меньше, чем в астеносфере. Многолетними (1977-1996 гг.) сейсмологическими наблюдениями установлено, что в этой зоне существуют длительно действующие источники генерации глубоких длиннопериодных (ГДП) землетрясений, которые имеют иную генетическую природу, чем вулкано-тектонические и тектонические землетрясения в исследуемом районе. В качестве гипотезы, наиболее удовлетворительно объясняющей наблюдаемые кинематические и динамические характеристики ГДП событий под Ключевским вулканом, предлагается рассматривать возникновение импульсов давления в магме. Причиной появления таких импульсов могут быть быстропротекающие фазовые переходы в ограниченных объемах метастабильной магмы. Таким процессом, по нашему мнению, может быть отделение газа в следствие спонтанной (взрывной) полимеризации низкомерных форм силикатов в расплаве.

Введение

В настоящее время все разнообразие сейсмических сигналов, регистрируемых на вулканах, ряд исследователей [3,31,32 и др.] предлагают разделить на две основные группы в зависимости от физической природы генерирующих источников. 

Первая группа включает сейсмические события, возникающие в твердой среде в результате хрупкого разрушения пород под действием сдвиговых и растягивающих напряжений, создаваемых активными магматическими процессами, и распределенных в достаточно большом объеме, окружающем магматические каналы и очаги. Эти землетрясения называются вулкано-тектоническими (ВТ), поскольку по механизму очага [7,21] и виду записи они практически неотличимы от чисто тектонических землетрясений. ВТ землетрясения характеризуются четкими вступлениями волн P и S и быстро спадающей кодой, преобладающие частоты лежат в широком диапазоне от 1 до 15 Гц [40]. Термин "вулкано-тектонические землетрясения" впервые был предложен Г.С.Горшковым при изучении сейсмичности Северной группы вулканов на Камчатке [4]. ВТ землетрясения возникают в районах активных и потухших вулканов как единичные события, так и роями длительностью от часов и дней до недель или месяцев. Часто эти землетрясения являются первыми признаками возобновившейся вулканической активности (долгосрочные предвестники) [22,29,41,45 и др.]. Известно множество примеров значительного усиления ВТ сейсмичности непосредственно перед началом извержений (краткосрочные предвестники) [2,8,22 и др.]. Вторая группа сейсмических сигналов состоит из источников, в которых активную роль в генерации сейсмических волн играют флюиды. К этой группе землетрясений относятся так называемые длиннопериодные (ДП) события: неглубокие ДП землетрясения, взрывные землетрясения, вулканическое дрожание (в.д.), сейсмические сигналы, сопровождающие дегазационные процессы в магме, и другие эруптивные явления [31]. По сравнению с ВТ, ДП землетрясения характеризуются более низкочастотной квазимонохроматической кодой, четко выраженной пиковой формой спектра, диапазон преобладающих частот лежит в интервале от 0,5 до 5 Гц [45]. Как показало теоретическое моделирование [4,32-36], характеристики природных сигналов, зарегистрированных от ДП событий на вулканах, подобны записям синтетических сигналов, генерируемых трещинами или каналами, заполненными жидкостями и флюидами и возбужденными меняющимся давлением. В этом смысле ДП землетрясения могут быть индикаторами пульсаций или перестройки давления в динамически активных магматических системах и, следовательно, содержать ценную прогностическую информацию как о начале, так и о развитии извержений, в частности, что особенно важно, о приближении катастрофической стадии извержения, как это было в 1991 г. на вулкане Пинатубо [37,45]. Так же как ВТ землетрясения, ДП землетрясения происходят в виде единичных событий, роев и продолжительных серий накладывающихся друг на друга сигналов, которые в ряде случаев удается связать с движением магмы [2,38,39]. 

Важной особенностью ДП сейсмичности является подобие формы записи индивидуальных сигналов, составляющих серии и рои, что свидетельствует о неоднократном возбуждении стационарного источника при недеструктивном процессе [34]. Наблюдается сходство пространственных, временных и спектральных характеристик вулканического дрожания и ДП сигналов, подтверждающее общность физических процессов в источниках этих событий [32,39,43]. ДП землетрясения и ДП в.д. были зарегистрированы на многих вулканах (см., например, обзор литературы в [43]). Большинство этих событий возникают на глубинах менее 3-5 км. Это поверхностные длиннопериодные (ПДП) землетрясения, непосредственно связанные с вулканической активностью. ДП землетрясения на глубинах более 10 км относятся к глубоким длиннопериодным (ГДП) землетрясениям, они наблюдались не так часто [45], и их связь с извержениями менее ясна [30].

Настоящая статья посвящена исследованию глубоких длиннопериодных (ГДП) землетрясений под базальтовым Ключевским вулканом - одним из наиболее активных вулканов Северной группы на Камчатке. В статье приводятся результаты изучения сейсмичности слоя, в котором происходят ГДП землетрясения, частично публикуемые в настоящем сборнике в работе [6]. На основании особенностей характеристик ГДП землетрясений делается попытка обсудить возможные причины их возникновения. 

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

Рис. 1

 Использованные данные и методы исследования. В работе использованы каталоги землетрясений района Северной группы вулканов за 1974-1996 гг., составленные по данным сети стационарных сейсмических станций (с/ст) с однотипной амплитудно-частотной характеристикой сейсмических каналов [23] и сети радиотелеметрических станций (РТС). Подробнее о развитии сейсмических наблюдений в этом районе можно прочитать в работе [8]. Наиболее однородный каталог получен авторами [6] для периода 1977-1996 гг. Содержащиеся в нем землетрясения обработаны по единой программе локализации гипоцентров с учетом местного годографа сейсмических волн для земной коры [1,21]. Количественные оценки параметров сейсмического режима Ключевского вулкана проводились по выборке данных, заключенных в прямоугольнике, показанном на рис.1 в статье [6] настоящего сборника. Этот прямоугольник схематично представляет Центральную сейсмоактивную зону Ключевского вулкана, в которой в последние десятилетия была сосредоточена его основная магматическая деятельность. Для выявления пространственно-временных особенностей сейсмичности анализировались как полный каталог, включающий все события, координаты которых было возможно определить, так и выборки землетрясений, локализованных с заданной точностью. Кроме данных, полученных при обработке аналоговых записей смещения грунта (1977-сентябрь 1996 гг.), при исследовании кинематических и динамических характеристик ВТ и ГДП землетрясений были использованы (начиная с октября 1996 г.) цифровые записи скорости смещений, любезно предоставленные лабораторией телеметрических станций КОМСП Геофизической Службы РАН (заведующий лабораторией В.В.Ящук). На рис.1,А представлен район Центральной сейсмоактивной зоны, в котором проводились исследования землетрясений, и его вертикальный разрез (рис.1 ,Б) с гипоцентрами землетрясений, зарегистрированных в 1977-1996 гг. не менее чем четырьмя сейсмическими станциями с погрешностью не более 3 км по глубине и эпицентру. Темными точками на рис.1 ,А обозначены эпицентры исследованных нами ГДП землетрясений, зарегистрированных сетью РТС в 1996 г.

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

Энергетический класс землетрясений Ks = log E (Дж) оценивался по S-волнам по номограмме С.А.Федотова [24]. Ks = 4,6 + 1,5 M, где M - магнитуда, определенная по поверхностным волнам.

Для исследования спектральных характеристик и распределения землетрясений по энергетическим классам использовался пакет статистического программного обеспечения STATISTICA , разработанный фирмой StatSoft Inc., USA . 

Краткая характеристика сейсмичности Ключевского вулкана, связанной с его современной магматической деятельностью. Как было установлено в результате многолетних сейсмологических наблюдений [5,6,9,26], сейсмичность Ключевского вулкана, связанная с его магматической деятельностью в течение последних десятилетий, проявлялась, главным образом, в Центральной сейсмоактивной зоне, охватывающей область вокруг вершинного кратера и круто падающей вниз почти до 40 км. На основании распределения числа и уровня сейсмической энергии землетрясений по глубине в земной коре под вулканом Ключевской были выделены четыре сейсмоактивных горизонта: -4-5 км; 5-2 км; 12-20км и 20-35 км. Подробные характеристики слоев даны в работе [6] настоящего сборника. 

В слоях 1-3 происходят преимущественно ВТ землетрясения, возникающие в упруго-деформированной среде вокруг магматических каналов и магматических очагов под действием меняющегося поля напряжений. Самый высокий уровень сейсмичности, отмечается на глубинах менее 5 км (слой 1), где происходят многочисленные рои землетрясений, связанные с внедрением даек и силлов в вулканическую постройку. Землетрясения в интервале глубин 5-20 км (слои 2 и 3) являются индикаторами подвижек, возникающих на путях миграции магмы в средних горизонтах земной коры, через которые проходит магмоподводящий канал или система каналов. 

В нижних горизонтах земной коры и переходном от коры к мантии слое 4 (интервал глубин 20-40 км) под Ключевским вулканом выявлена зона, в которой генерируется множество ГДП землетрясений [6]. Слой 4 характеризуется максимальным количеством происходящих в нем землетрясений и относительно низкой суммарной сейсмической энергией (более чем на порядок ниже, чем в слое 1). Максимальный энергетический класс землетрясений, происходящих в этом слое, не превышает 6,5-6,7. По приближенным оценкам, величина коэффициента в распределении числа землетрясений по энергии аномально высока: ~ 2,6 [6]. Существенные отличия землетрясений слоя 4 от ВТ землетрясений, происходящих в трех вышележащих горизонтах: низкий максимальный энергетический класс, однотипная форма записи сигналов, более длинные преобладающие периоды сейсмических колебаний [6], позволяют отнести землетрясения, происходящие в низах коры и переходном от коры к мантии слое под Ключевским вулканом, к разряду глубоких длиннопериодных (ГДП) событий.

В работах [6,10] отмечалось, что увеличение числа ГДП землетрясений, как правило, происходит до или сразу после значительных вулканических событий (начало вершинных или побочных извержений, усиление эффузивной и/или эксплозивной деятельности в ходе извержений), опережая на срок от нескольких дней или недель до нескольких месяцев возникновение неглубоких ВТ и ПДП землетрясений. Во многих случаях наблюдалось два типа характерных временных последовательностей сейсмических и вулканических событий:

 группы и рои длиннопериодных землетрясений на глубинах 20-35 (40) км единичные землетрясения в интервале 5-20 км рои ВТ землетрясений в вулканогенно-осадочном слое и вулканической постройке (-4-5 км) появление и/или усиление в.д. и других ПДП сигналов начало или усиление вулканической активности;

 окончание извержения ослабление и/или исчезновение в.д. уменьшение числа поверхностных ВТ и ДП землетрясений появление единичных землетрясений на глубинах 5-20 км возобновление сейсмической активности в интервале глубин 20-40 км. 

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

Рис. 2

 Пространственное распределение землетрясений в средних и нижних горизонтах земной коры на разных стадиях активности Ключевского вулкана. Нами были исследованы особенности пространственного распределения очагов землетрясений на глубинах 5-20 и 20-40 км в Центральной сейсмоактивной зоне относительно кратера Ключевского вулкана и центров побочных извержений в 1974-1990 гг. в зависимости от усиления или ослабления вулканической деятельности. 

Было обнаружено, что во время подготовки или усиления извержений очаги землетрясений в средних горизонтах земной коры Центральной сейсмоактивной зоны под Ключевским вулканом располагаются в непосредственной близости от вершинного кратера, охватывая цепочку "свежих" побочных эруптивных центров или хотя бы часть из них (рис.2 , А).

Рис. 3

 В периоды образования провалов в вершинном кратере Ключевского вулкана и ослабления или прекращения извержений, одновременно с перемещением сейсмичности в низы коры на глубинах 5-20 км отмечалось удаление очагов землетрясений от вершинного кратера и зоны побочных конусов в периферические участки Центральной зоны, на расстояния до 10-15 км и далее (рис3 ,А). Особенно ярко этот эффект проявился после пароксизмального вершинного извержения 7.09.-2.10.1994 г., когда в кратере вулкана в результате значительного опустошения магматического канала начал формироваться огромный провал, дважды зафиксированный аэрофотосъемкой (в августе 1995 г. и в сентябре 2000 г.) [12]. В течение 1995-1999 гг. происходило смещение сейсмической активности на глубинах 5-20 км к северу, западу и северо-западу от кратера Ключевского вулкана к вулканам Ушковский и Крестовский. Наиболее сильные землетрясения были зарегистрированы в районе РТС PDK и под вулканом Крестовский. 

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

Из выше сказанного следует, что:

 ВТ землетрясения на глубинах 5-20 км возникают, по-видимому, в основном при перераспределении напряжений в среде, вмещающей питающую магматическую систему Ключевского вулкана, давление в которой меняется в зависимости от условий накопления или поступления магмы.

 Зона генерации ГДП землетрясений более устойчива в пространстве и не зависит от перераспределения давления в магматической системе.

Рис. 4

 Исследование характеристик ГДП землетрясений. По данным аналоговых записей был проведен сравнительный анализ следующих характеристик ГДП землетрясений и землетрясений, происходящих на глубинах -4-20 км: формы записи, распределения числа землетрясений N по энергетическим классам Ks, периода Ts в максимальной амплитуде колебаний поперечных волн, коэффициента корреляции зависимости Ts = Ts(Ks).

На рис.4 ,А приведены образцы аналоговых записей землетрясений по смещению с энергетическим классом Ks ~ 6, происходящих в выделенных слоях Центральной сейсмоактивной зоны под Ключевским вулканом и вне этой зоны. ГДП землетрясения (рис.4 ,Б) преимущественно однотипны по форме записи, которая в значительной степени отличается как от записей землетрясений в вышележащих слоях (рис.4 ,А), так и от записей тектонических землетрясений (рис.4 ,В) на глубине более 20 км, в которых в ближней зоне регистрации амплитуда коды S-волн спадает по экспоненциальному закону. Средний видимый период колебаний Ts (в S- волнах максимальной амплитуды) для землетрясений слоя 4 - наибольший по сравнению с аналогичными периодами землетрясений, которые происходят в других горизонтах под Ключевским вулканом и в других сейсмоактивных зонах, существующих в районе вулканов Ключевской группы [9]. Для землетрясений слоя 4 с энергетическим классом Кs ~ 6 периоды Ts составляют на эпицентральных расстояниях от 20 км до 45 км ~ 0,6-0,9 с. Для землетрясений слоев 2 и 3 под Ключевским вулканом Тs = 0,3-0,5 с, а в слое 1~ 0,5-0,6 с.

Рис. 5

 Нами были построены гистограммы распределения числа землетрясений по энергетическим классам для 1-3 слоев (рис.5 ,А) и в слое 4 (рис.5 ,Б). При построении использовались данные за 1978-1995гг по 1351 землетрясению в интервале глубин -4-20 км и данные по 2015 землетрясениям с глубиной 20-40 км. Очевидно, что распределение числа землетрясений по Кs для глубин 20-40 км близко к нормальному закону и отличается от распределения для верхних горизонтов глубин, сходного с распределением для тектонических землетрясений (см., например, [28]). 

В последние годы (1997-2000 гг.) в районе Ключевского вулкана увеличилось число сейсмостанций сети РТС и повысилась точность уверенной регистрации слабых землетрясений. Но вид гистограмм распределения землетрясений на глубинах -4-20 км и 20-40 км, происходивших в эти годы, не отличается от приведенных на рис.5 ,А,Б.

Не было выявлено статистической зависимости энергетического класса Кs ГДП землетрясений от периода колебаний в поперечных волнах Ts, что не характерно для тектонических [16] и вулкано-тектонических землетрясений, связанных с образованием трещин при движении магмы. Коэффициент корреляции, рассчитанный для ГДП землетрясений, равен -0,57.

Рис. 6

 Кроме аналоговых записей нами были использованы цифровые данные сети РТС в 1996 г. Для 25-ти ГДП событий, обработанных по цифровым записям скоростей смещения, были построены энергетические спектры 3-х компонент: С-Ю, В-З и Z волн P и S для сравнения со спектрами Р и S волн тектонических землетрясений такого же энергетического класса (Ks ~ 6) и происходящих на таких же глубине и эпицентральном расстоянии от пункта регистрации, что и ГДП землетрясения, но за пределами Центральной сейсмоактивной зоны на удалении ~35 км от кратера вулкана Ключевской. Характерный вид графиков нормированной спектральной плотности вертикальной компоненты Z для ГДП и тектонических землетрясений, записи которых приведены на рис.4 ,Г,Д, представлен на рис.6 . Из рис.6 следует, что спектры Р и S волн ГДП землетрясений подобны и имеют один резко выраженный пик на более низких частотах, чем у тектонического землетрясения, удаленного от вулкана и имеющего, по-видимому, обычный для тектонических землетрясений сдвиговый механизм. На рис.6 видно также, что для ГДП землетрясений частота колебаний в S-волнах приблизительно равна частоте колебаний P-волн. Аналогичная картина наблюдается и для спектров горизонтальных каналов. 

Таким образом, по итогам работы [6] и в результате исследований ГДП землетрясений, выполненных в настоящей статье, было выявлено, что:

Эпицентры ГДП землетрясений расположены на глубине 20 - 35 км и сконцентрированы в районе центрального кратера, независимо от проявлений вулканической активности.

ГДП землетрясения чаще происходят роями продолжительностью до 9 месяцев, а иногда и дольше.

Для ГДП землетрясений характерна однотипная форма записи.

Энергетические спектры P и S волн ГДП землетрясений имеют один резко выраженный пик на более низких частотах, чем для P и S волн тектонических землетрясений такого же энергетического класса.

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

Наиболее характерные особенности динамических и кинематических параметров ГДП землетрясений: а) равенство периодов колебаний в P и S волнах; б) на Z-канале амплитуды колебаний в P и S волнах примерно равны; в) как правило, в пункт регистрации первой вступает волна сжатия.

Обсуждение результатов 

Из анализа сейсмичности слоя 4 совместно с проявлениями вулканической активности можно предположить непосредственную связь ГДП землетрясений с магматической деятельностью. Об этом свидельствует появление роев ГДП землетрясений перед активизацией вершинного кратера или латеральными извержениями [6], а также особенности пространственного распределения очагов землетрясений на глубинах 5-20 и 20-40 км на разных стадиях состояния вулкана. 

Пространственное распределение гипоцентров землетрясений в нижних горизонтах земной коры под Ключевским вулканом обрисовывает систему трещин, по которым магма может продвигаться к поверхности. На аналоговых и цифровых записях ГДП землетрясений четко прописываются S-волны, из чего следует, что на глубине 20-40 км располагается не сплошная магматическая колонна, а среда с системой трещин, в которой могут образовываться и через которую могут проходить поперечные сейсмические волны. По-видимому, на глубине 20-40 км мы имеем дело с системой заполненных магмой мелких трещин, являющихся верхней частью диапира, расположенного глубже 40 км в асейсмичной области между земной корой и верхней мантией, где под Ключевской группой вулканов ранее были обнаружены аномалии затухания короткопериодных сейсмических волн, связываемые с существованием магматических очагов [27]. Относительно небольшие энергетические параметры ГДП землетрясений в сравнении с более мелкофокусными ВТ землетрясениями свидетельствуют о том, что в нижних слоях земной коры за счет свойств среды не происходит интенсивного накопления напряжений. Следовательно, эта среда более пластична (по-видимому, за счет теплового потока от магматического очага).

О механизме генерации ГДП землетрясений. На основании полученных данных нами были рассмотрены две возможные гипотезы о механизме генерации ГДП землетрясений под Ключевским вулканом:

Возникновение мелких трещин в упруго-пластичной среде при движении магмы.

Импульсы давления в магме, заполняющей уже существующие трещины. 

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

Исходя из гипотезы 1, средний размер источника r ГДП землетрясений можно оценить по формуле Брюна [30]:

r = 0,35 VP / f0 , где

VP - cкорость продольных волн, VP 7,5 км/с [1] и

f0 - угловая частота, f0 = 1,5-2,0 Гц. 

Таким образом, мы получим, что r = 1,3 ....1,8 км.

Учитывая размеры источников, рассчитанных по формуле Брюна, справедливой для землетрясений, имеющих сдвиговый механизм, трудно предположить, что в небольшом объеме среды, ограниченном гипоцентрами землетрясений, способно появиться такое количество трещин (сотни) в течение короткого интервала времени (продолжительность одного роя) длиной 1-2 километра. Определенные трудности возникают при интерпретации на основе гипотезы 1 появления роев ГДП землетрясений после снижения внешних проявлений вулканической активности. Кроме того, на основе гипотезы 1 достаточно трудно интерпретировать пункты 3-6 из выше перечисленных характерных особенностей ГДП землетрясений, т.к. параметры землетрясений, связанных с движением магмы в интервале глубин -4-20 км под вулканами, также как и тектонических землетрясений, отличаются от параметров ГДП землетрясений. 

Обсудим подробнее эти различия.

Одна из основных особенностей ГДП землетрясений, выявленных при обработке аналоговых записей - это различные законы распределения ГДП землетрясений и землетрясений на глубинах до 20 км. Протекающие в природе случайные процессы, зависящие от многих факторов, описываются суперпозицией простых законов распределения. Так, для тектонических землетрясений закон повторяемости Гуттенберга-Рихтера, хорошо аппроксимируемый правой ветвью распределения Пуассона, определяется такими основными факторами, как линейный размер блоков и характерное время процесса (время между землетрясениями) [19]. Землетрясения, происходящие на глубинах до 20 км в исследуемом нами районе, также достаточно хорошо укладываются в распределение Пуассона. Коэффициент , рассчитанный для правой ветви распределения ВТ землетрясений по энергетическим классам, примерно равен для тектонических землетрясений Камчатки [25]. Это свидетельствует о том, что механизм очагов подавляющего большинства ВТ землетрясений не отличается от механизма очагов тектонических землетрясений.

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

Важной особенностью ГДП землетрясений является равенство периодов колебаний в P и S волнах, что отличает их от тектонических землетрясений. Это присуще и длиннопериодным землетрясениям, зарегистрированным во время извержения вулкана Пинатубо [45]. Заметим, что для тектонических землетрясений частота колебаний в S-волнах, как правило, ниже, чем в P-волнах (см., например, [16]). Следует также отметить, что у всех исследованных нами ГДП землетрясений с четкими вступлениями P-волн первой проходит волна сжатия, что также не характерно для тектонических и ВТ землетрясений. 

Все это позволяет предположить иную генетическую природу ГДП землетрясений, отличную от природы обычных ВТ землетрясений. Об этом косвенно свидетельствует идентичность формы записей ГДП землетрясений в отличие от землетрясений других слоев, что было отмечено и для длиннопериодных землетрясений при извержении вулкана Пинатубо [41].

Обсудим далее гипотезу 2, впервые предложенную в работе [34], в рамках которой неплохо объясняются перечисленные выше особенности ГДП землетрясений (пункты 1-6) в слое 4 в Центральной зоне под Ключевским вулканом. В этой модели ДП землетрясения генерируются при колебаниях стенок заполненных магмой трещин в результате импульсов давления в магме; волны давления в магме, как предполагают авторы [31,32,34], возникают при движении магмы со сверхзвуковой скоростью. Последнее предположение, по нашему мнению, является самым уязвимым звеном модели, так как это означает движение потока магмы со скоростью более 1 км/с, что представляется нам весьма проблематичным.

С нашей точки зрения наиболее характерные особенности (пункты 3-6) ГДП землетрясений лучше объяснить неизменными в течение длительного времени физико-химическими свойствами магмы, заполняющей трещины, и характерными особенностями процессов, протекающих в магматическом расплаве и являющихся причиной возникновения в нем импульсов давления. Форма записи ГДП землетрясений также свидетельствует об импульсной природе их источника [13]. 

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

2 [ SiO3 (OH)] 3- распл. = [ SiO2 O7] 6-распл. + H2O пар.

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

Из экспериментальных данных [43] известно, что в базальтовых расплавах при внешнем давлении P = 1 атм (105 Па) газовый пузырек может появиться и вырасти до микронных размеров за время, измеряемое миллисекундами. Уровень метастабильности при этом может составлять несколько атмосфер [43]. Но процесс газоотделения на глубине 20-40 км (внешнее давление P = (8-16).108 Па) и при P = 1 атм (105 Па) должен иметь некоторые различия. Во-первых, процесс газоотделения на большой глубине из-за большого гидростатического давления можно рассматривать как происходящий в замкнутом объеме. Во-вторых, отличаются физико-химические свойства магматического расплава, определяющие скорость роста пузырьков при различном давлении, - газосодержаниие, вязкость, коэффициент диффузии. Эти факторы влияют на скорость появления и роста пузырьков. 

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

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

Исходя из предлагаемой гипотезы, попытаемся связать период сейсмических колебаний ГДП землетрясений с физическими характеристиками магматического расплава и сравнить с экспериментальными данными. Для этого определим основные факторы, влияющие на рост пузырьков. 

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

Решение задачи о динамике газовых пузырьков основывается на уравнении Рэлея-Тейлора c учетом вязкости [42]:

(1)

где R - радиус пузырьков, m - плотность магмы, - кинематическая вязкость магмы, P - давление в магме, P1 - давление в газовом пузырьке.

Для того чтобы оценить время релаксации избыточного давления, возникающего при расширении пузырьков, пренебрегаем первыми двумя членами в уравнении (1) вследствие их малости:

(2)

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

(3)

где - коэффициент поверхностного натяжения.

Подставляя выражение (3) в (2), получим уравнение:

,

из которого интегрированием по времени можно определить длительность импульса:

=. (4) 

Для водонасыщенных базальтовых магм на глубине 20 - 40 км по данным [18] вязкость = m ~ 101 -102 Пуаз (1-10 Па/с), а 102 дин/см (10-1 Па/c) [17]. Как следует из формулы (4), для ~ 1с радиус R возникших в магме пузырьков должен составлять 10-1-10-2см, что на несколько порядков превосходит расчетные [20] и экспериментальные [43] данные. Следовательно, время в нашей модели должно определяться другими параметрами.

Рассмотрим теперь процесс роста пузырьков газа исходя из уравнения диффузии:

(5)

где D - коэффициент диффузии, (для базальтовой магмы Ключевского вулкана при P = 5000 атм ( 5.108 Па ) D ~ 10-9 м2/c [16]); dM - масса газа, диффундирующего в пузырек за время dt через поверхность пузырька ds; d g/dx - градиент концентрации свободного газа в магме, равный:

где W(P) - весовая концентрация растворенного в магме газа, dW(P)/dP