Тонкая структура косейсмического электромагнитного отклика по данным геомагнитных и сейсмологических наблюдений

Полный текст

ВВЕДЕНИЕ

Основная проблема в изучении геомагнитного отклика от состоявшихся землетрясений связана со сложностью получения высококачественных данных из пунктов наблюдений, обеспечивающих мониторинг одновременно геомагнитной и сейсмической активностей в достаточной близи от очагов сильных землетрясений. Для возможности детального исследования такого отклика должны быть выполнены три условия: экстремально большая магнитуда землетрясения (Mw ≥ 7), наличие качественных 1-секундных наблюдений геомагнитного поля на относительно небольших (не более 25° или 3000 км) эпицентральных расстояниях, наличие сейсмической станции рядом с геомагнитной обсерваторией [Соловьев, 2023]. Регистрация секундных наблюдений геомагнитного поля приемлемого качества началась примерно с 2020 г. На сегодня ее обеспечивают несколько десятков обсерваторий, расположенных, в основном, в европейском секторе [Thomson, Flower, 2021]. В итоге лишь единичные сейсмические события удовлетворяют всем перечисленным условиям. К ним, в частности, относятся несколько землетрясений, произошедших в 2020–2023 гг.

30.10.2020 г. в 10 км к северу от греческого острова Самос в Эгейском море произошло землетрясение с магнитудой Mw = 7.0 на глубине 21 км (по данным Геологической службы США (USGS), https://earthquake.usgs.gov/earthquakes/). Данная территория на севере граничит с зоной северо-анатолийского разлома, на юге – с Эллинским желобом [Sakkas, 2021]. 06.02.2023 г. в провинции Кахраманмараш, расположенной на юге Турции, произошли разрушительные землетрясения, наиболее сильные из которых характеризовались магнитудами Mw = 7.5 и 7.8. Согласно USGS глубины гипоцентров составили 7.4 и 10 км соответственно. Сейсмичность данного региона, представляющего собой сочленение крупных Африканской, Анатолийской и Аравийской тектонических плит, изучена весьма подробно [Kadirioğlu, Kartal, 2016; Ahmad et al., 2017; Akkar et al., 2018; Alver et al., 2019; Güvercin et al., 2022]. Журнал “Физика Земли” № 6, 2023 г. целиком посвящен всестороннему анализу серии турецких землетрясений 2023 г. В основном он касается изучения процессов в очаге, тектонических аспектов, реконструкции сопутствующих геодинамических процессов и др. (см.введение к выпуску [Смирнов,Шебалин, 2023 ] )

Вместе с тем, по-прежнему остаются недостаточно исследованными механизмы воздействия землетрясений и связанных с ними процессов в литосфере на геомагнитное поле. Их изучению посвящена отдельная дисциплина – сейсмоэлектромагнетизм [Molchanov, Hayakawa, 2008]. На сегодня преобладают два описания возможного механизма генерации геомагнитного отклика на сейсмическое событие. Первое связано с так называемым “косейсмическим сигналом”, предполагающим генерацию электромагнитного сигнала сейсмической волной. Генерация таких сигналов обусловлена электрокинетическим или геомагнитно-индуктивным эффектом. Электрокинетический эффект возникает при разделении зарядов во флюиде, содержащемся в земной коре. Деформации коры, вызванные сейсмическими волнами, приводят в движение ионизированный флюид, что создает электрокинетический ток. Геомагнитно-индуктивные возмущения при распространении сейсмических волн могут возникать в результате эффекта индукции за счет колебаний проводящей земной коры, находящейся в постоянном геомагнитном поле. На основе численного моделирования косейсмического эффекта было показано, что длительность сигналов должна увеличиваться с расстоянием из-за дисперсии поверхностных сейсмических волн [Surkov et al., 2018; Yamazaki, 2024]. Второй механизм предполагает распространение акустико-гравитационных волн от очага с последующим возбуждением ионосферы и генерацией геомагнитных пульсаций. Очевидно, в этом случае речь идет о более длиннопериодных сейсмогенных вариациях геомагнитного поля. Так, например, в работе [Рябова, Шалимов, 2022] геомагнитные вариации с периодами больше 13 мин интерпретируются как результат распространения медленных магнитогидродинамических волн, возбужденных в ионосфере акустическим импульсом после землетрясений, или как результат прохождения перемещающихся ионосферных возмущений.

Явление возбуждения электромагнитного поля при распространении механических колебаний в пористых влагонасыщенных горных породах (также называемое сейсмоэлектрическим эффектом второго рода) экспериментально было открыто в 1939 г. [Иванов, 1939]. Интерес был обусловлен тем, что возбуждаемые поля несут информацию о глубине залегания влагонасыщенной породы, ее флюидосодержании, пористости, проницаемости и других петрофизических параметрах. Данный механизм расценивался как перспективный при решении разнообразных прикладных задач инженерной и нефтегазовой геологии, геофизики и геодинамики. Позже, в 1960-х годах, совершались попытки использовать сейсмоэлектрический эффект для детектирования подземных ядерных испытаний. Тем самым, регистрация сейсмогенных электромагнитных возмущений выполнялась на нерегулярной основе в интересах сугубо прикладных задач.

В работе [Соловьев, 2023] сообщалось о детектировании геомагнитного отклика от природных землетрясений в высокоточных данных регистрации магнитного поля с 1-секундной частотой дискретизации. Было показано, что при магнитуде Mw ≥ 7.0 сигнал надежно детектируется в радиусе 3000 км от эпицентра. При этом, геомагнитный эффект наилучшим образом выделяется в скорости изменения поля, достигая аномальных амплитуд >10 нТл/с. Для каждой из рассмотренных геомагнитных обсерваторий была определена задержка прихода сигнала от землетрясений. Естественным продолжением данного исследования является более точная оценка скоростей распространения электромагнитного сигнала для установления приоритетного механизма его передачи. Часть обсерваторий, данные которых рассматривались в работе [Соловьев, 2023], обеспечивают комплексный набор геофизических измерений, включающих не только геомагнитные, но и сейсмологические наблюдения. Таким образом, появляется возможность комплексирования геомагнитных и сейсмических данных, синхронно регистрируемых на едином пункте наблюдений. Именно этому исследованию и посвящена настоящая работа.

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

ИСХОДНЫЕ ДАННЫЕ

Основные затруднения в изучении обсуждаемых геомагнитных эффектов связаны с недостатком 1-секундных наблюдений геомагнитного поля и их неполным географическим охватом. Как было показано в работе [Соловьев, 2023], искомый эффект в исходных данных представляет собой квазипериодический сигнал с характерным периодом 5–20 с, небольшой амплитудой порядка нескольких нТл и длительностью от 5 до 15 мин. Таким образом, в 1-минутных наблюдениях, регистрация которых была начата существенно раньше и которых накоплено значительно больше по всему земному шару, исследуемые эффекты пропадают за счет сглаживания данных при усреднении. Кроме того, согласно работе [Соловьев, 2023] на расстоянии более 3000 км от эпицентра сильного землетрясения сейсмогенный геомагнитный эффект рассеивается и перестает детектироваться в скорости изменения магнитного поля, поэтому выбор данных ограничивается обсерваториями, расположенными в указанном радиусе от очагов изучаемых землетрясений.

Учитывая эти обстоятельства, были отобраны обсерватории, обеспечивающие мониторинг геомагнитного поля с 1-секундным разрешением. Из них мы оставили те, у которых на территории или в непосредственной близости имелись постоянно действующие широкополосные сейсмические станции. Отобранные пункты наблюдений представлены на карте на рис. 1. Из них комплексный мониторинг геомагнитной и сейсмической активности обеспечивают обсерватории “Михнево” (Московская обл., IAGA-код магнитной обсерватории MHVg, код сейсмической станции MHVs) [Адушкин и др., 2016; Soloviev et al., 2023], “Борок” (Ярославская обл., IAGA-код магнитной обсерватории BOX, код сейсмической станции BROK) [Гоев и др., 2024] и “Климовская” (Архангельская обл., IAGA-код магнитной обсерватории KLI, код сейсмической станции KLM) [Соловьев и др., 2016; 2022; Soloviev et al., 2015; Антоновская и др., 2022]. Расстояния между парными пунктами наблюдений составляют от 30 до 350 км.

В результате, за период 01.01.2020–15.02.2024 гг. нашлось три землетрясения с магнитудой Mw ≥ 7 в 3000-километровой окрестности рассматриваемых обсерваторий (см. рис. 1):

1. 30.10.2020 11:51 UT Mw = 7.0 H = 21 км, вблизи о. Самос (Греция);
2. 06.02.2023 01:17 UT Mw = 7.8 H = 10 км, провинция Караманмараш (Турция);
3. 06.02.2023 10:24 UT Mw = 7.5 H = 7.4 км, там же.

Обозначим их условно ЗТ-1, ЗТ-2 и ЗТ-3. Их эпицентры отмечены на карте серыми кружками (см. рис. 1).

 

Рис. 1. Карта эпицентров землетрясений с магнитудой Mw ≥ 7 за период 01.01.2020–15.02.2024 гг. по данным USGS (серые круги), отобранные магнитные обсерватории, предоставляющие 1-секундные данные (черные звезды – обсерватории сети ИНТЕРМАГНЕТ, полые звезды – не входящие в ИНТЕРМАГНЕТ обсерватории), и сейсмические станции (обозначены треугольниками).

 

В табл. 1 перечислены пары рассматриваемых пунктов геомагнитных и сейсмических наблюдений, указаны расстояния от эпицентров до пунктов сейсмических наблюдений и время прихода на них P-волн. Для тех сейсмических станций, которые предоставляют данные в Международный сейсмологический центр ISC, информация о времени первого вступления была загружена с сайта ISC [International…, 2024a; 2024b]. Для остальных станций соответствующие временные отметки определялись вручную.

 

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

Сейсм. станция	Геомагнитная обсерватория	Расстояние от сейсм. станции до эпицентра (D, км)			Задержка прихода P-волны (dP, с)
		ЗТ-1	ЗТ-2	ЗТ-3	ЗТ-1	ЗТ-2	ЗТ-3
TASB	GLK	1530	700	650	258	94	89
KHC	BDV	1630	2300	2250	206	282	277
GKP	HLP	1850	2340	2280	235	284	275
MHVs	MHVg	2060	1970	1880	253	252	235
BROK	BOX	2390	2320	2230		284	272
PUL	SPG	2440	2550	2460	293	304	296
VJF	SPG	2510	2670	2590	300	317	306
KLM	KLI	2700	2630	2540	318
ADZR	KLI	2990	2970	2880	340	343	335
SHIP	KLI	2950	2780	2690	339	327	317

 

Все российские геомагнитные обсерватории (MHVg, BOX, SPG, KLI) и армянская обсерватория GLK входят в межрегиональный сегмент геомагнитных наблюдений, ядром которого является ЦКП “Аналитический центр геомагнитных данных” на базе Геофизического центра РАН [Гвишиани и др., 2018]. Данные этих обсерваторий по мере поступления в центр непрерывно обрабатываются средствами аппаратно-программного комплекса (АПК) МАГНУС на предмет распознавания техногенных и природных аномалий магнитного поля [Кудин и др., 2021; Kudin et al., 2023] и публикуются онлайн с присуждением DOI [Soloviev et al., 2015; 2020; 2021].

ОЦЕНКА СТАБИЛЬНОСТИ ПОСТАМЕНТОВ, ГЕОМАГНИТНОЙ ОБСТАНОВКИ И ШУМОВЫХ ХАРАКТЕРИСТИК

В обсерваторской практике калибровочная кривая dF служит важным индикатором корректной и согласованной работы векторного и скалярного магнитометров, которые непрерывно измеряют вариации трех компонент и модуль вектора магнитной индукции соответственно. Согласно стандартам ИНТЕРМАГНЕТ [St-Louis, 2020] dF рассчитывается по формуле:

dF(i) = F_s(i) – sqrt(X(i)² + Y(i)² + Z(i)²),

где F_s – измеренный модуль вектора магнитной индукции; X, Y, Z – северная, восточная и вертикальная компоненты вектора магнитной индукции; i – текущий отсчет по времени.

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

 

Рис. 2. Графики Fs (сверху) и dF (снизу) за 30.10.2020 г., построенные по секундным данным обсерваторий BDV (а) и KLI (б), и за 5–6 февраля 2023 г. по секундным данным обсерваторий GLK (в) и KLI (г). Моменты первого вступления обозначены вертикальными пунктирными линиями.

 

В анализе дополнительно использовались данные индексов геомагнитной активности, которые служат наглядным индикатором состояния космической погоды. Данные доступны на сайте Международной службы геомагнитных индексов ISGI (https://isgi.unistra.fr/). В течение 29–30 октября 2020 г. и 5–6 февраля 2023 г. наблюдался пониженный уровень геомагнитной активности: геомагнитные K-индексы на рассматриваемых обсерваториях и индекс Kp, отражающий общепланетарную обстановку, не превышали 4 баллов, а индекс Dst, характеризующий интенсивность магнитосферного кольцевого тока, не опускался ниже –25 нТл. Графики индексов приведены на рис. 3. Это обстоятельство исключает ложную трактовку выделяемого сигнала как эффекта космической погоды вместо эффекта сейсмического события, а также позволяет более отчетливо выделить геомагнитный отклик на общем фоне.

 

Рис. 3. Графики геомагнитных индексов Kp, ap (сверху) и Dst (снизу) за 29–30 октября 2020 г. (а) и 5–6 февраля 2023 г. (б).

 

На рис. 4 приведены характерные шумовые характеристики 1-секундных векторных измерений магнитного поля, которые используются в настоящей работе. На рисунке приведены результаты, полученные по вертикальной Z-компоненте, наиболее чувствительной к внешним электромагнитным шумам. Из графиков видно, что характерный нижний порог амплитуды шумов составляет 0.1–0.2 нТл. К сожалению, далеко не все обсерватории, на которых осуществляется регистрация 1-секундных наблюдений, обеспечивают такое качество измерений, что не позволяет их использовать в изучении геомагнитного отклика от землетрясений. Например, если геомагнитная обсерватория находится в зоне относительно сильных электромагнитных шумов, то это вынуждает срезать частоты интересующего диапазона (>0.05 Гц) в 1-секундных геомагнитных измерениях на уровне АЦП.

 

Рис. 4. Шумовые характеристики данных вариометров обсерваторий KLI (а) и MHVg (б) по 1-секундным наблюдениям вертикальной Z-компоненты.

 

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Для корректного качественного сопоставления 1-секундных данных геомагнитных обсерваторий с более высокочастотными сейсмическими записями, исходная частота дискретизации которых составляет от 20 до 100 Гц, последние предварительно были усреднены до 1 Гц. Для каждой пары наблюдений, приведенных в табл. 1, были построены суточные графики по каждой из трех компонент смещений и dB/dt, синхронизированные по времени, в исходных единицах и в нормированном виде за вычетом среднего. Примеры по отдельным парам наблюдений исследуемых землетрясений ЗТ-1, ЗТ-2 и ЗТ-3 приведены на рис. 5 (на верхнем графике рис. 5а сейсмическая запись отсутствует ввиду того, что в это время соответствующий канал на станции не работал). На графиках сейсмических записей, зарегистрированных на нескольких станциях, отчетливо видны вступления основных сейсмических фаз, как объемных – первичной (P) и вторичной (S), так и поверхностных – Рэлея (RL) и Лява (LV). На сейсмических и магнитных станциях ведется запись вертикальной (Z) и двух горизонтальных (N, E) компонент, ориентированных в плоскостях Север–Юг и Восток–Запад соответственно. Для удобства представления и дальнейшего анализа следует перейти к системе координат, в которой одна из горизонтальных осей (радиальная R) лежит в плоскости падения, а другая (тангенциальная T) – перпендикулярна этой плоскости. Технически, это соответствует повороту вокруг вертикальной оси на обратный азимут. В этом случае, вступление P-волны и SV-волны проявляется на радиальной компоненте, а поперечная волна с поляризацией SH – на тангенциальной. То же относится к поверхностным волнам: запись волны Рэлея содержится на компонентах Z и R, волны Лява – на T-компоненте. Следует отметить, что для данной конфигурации взаимного расположения эпицентров землетрясений и станций в почти половине случаев обратный азимут близок к 180°. Очевидно, что в этом случае описанное выше преобразование не требуется. Далее будет рассмотрен электромагнитный отклик среды объемных и рэлеевских волн.

 

Рис. 5. Примеры совмещенных графиков северной (верхний), восточной (средний) и вертикальной (нижний) компонент смещений (синий) и dB/dt (оранжевый) в нормированном виде за вычетом среднего по данным пары наблюдений KLM-KLI для ЗТ-1 (а); VJF-SPG для ЗТ-2 (б); KHC-BDV для ЗТ-2 (в); MHVs-MHVg для ЗТ-3 (г); TASB-GLK для ЗТ-3 (д). Вертикальной пунктирной линией обозначен момент прихода P-волны на соответствующую сейсмическую станцию.

 

Наилучшим образом геомагнитный сигнал от землетрясений выделяется на записях скорости изменения геомагнитного поля dB/dt в виде сильных изолированных всплесков с последующим затуханием. При этом, чем ближе геомагнитная обсерватория к эпицентру, тем сигнал более локализован во времени (рис. 5д) и его амплитуда выше. Напомним, что при горизонтальной однородности геоэлектрических свойств подстилающей среды ориентация вектора dB/dt соответствует направлению возбуждаемого теллурического поля Е и токов в поверхностных слоях Земли. Наиболее отчетливо всплески dB/dt детектируются на горизонтальных компонентах магнитного поля (рис. 5). Пиковые амплитуды dB/dt, вызванные тремя землетрясениями, приведены в табл. 2. Так, на ближайшей к эпицентрам ЗТ-2 и ЗТ-3 обсерватории GLK (Армения) амплитуда возмущения |dX/dt|, вызванного землетрясением ЗТ-2, достигла почти 17 нТл/с. Такое значение скорости изменения поля является экстремальным даже при очень сильных геомагнитных бурях на высоких широтах, где вариабельность поля максимальна. Сейсмогенный геомагнитный эффект в данных обсерватории BOX оказался слишком слабый из-за вынужденного среза частот > 0.05 Гц.

 

Таблица 2. Пиковые значения амплитуд dB/dt (нТл/с), вызванные землетрясениями (пустые ячейки свидетельствуют об отсутствии исходных данных либо о наличии в них сильных шумов)

Геомаг. обсерв.	ЗТ-1			ЗТ-2			ЗТ-3
	dX/dt	dY/dt	dZ/dt	dX/dt	dY/dt	dZ/dt	dX/dt	dY/dt	dZ/dt
GLK				–16.9	16.615	9.26	8.855	10.125	–5
BDV	–0.82	0.69	0.395	–0.85	–0.77	0.455	–0.93	–0.835	0.43
HLP	–1.185	1.66
MHVg	–0.36	–0.41		–0.865	1.195	0.33	–0.895	–1.295	–0.32
SPG				0.765	1.145		–0.85	0.99
KLI	–0.425	0.25	0.11	0.91	–1.375	–0.235	0.745	1.4	0.245

 

Для всех событий по всем парам геомагнитных и сейсмических наблюдений начало геомагнитного возмущения совпадает с приходом P-волны (см. табл. 1 и рис. 5). С приходом S-волны, чья амплитуда в несколько раз больше, геомагнитный сигнал усиливается. С приходом рэлеевской волны, имеющей максимальную амплитуду, геомагнитный отклик также становится максимальным. Такая динамика возмущений вполне очевидна, если рассматривать косейсмический эффект как линейный электромагнитный отклик на механическое воздействие сейсмических волн. Это полностью соответствует теоретической модели косейсмического эффекта, описанной в статье [Surkov et al., 2018], что говорит в пользу “косейсмического” механизма возбуждения регистрируемого электромагнитного сигнала.

Для сопоставления частотного состава сейсмических и геомагнитных записей были построены периодограммы для каждого события на временном интервале существования геомагнитных возмущений. Использовались нормированные данные для сопоставимости спектральной плотности мощности разнородных наблюдений. При этом рассматривался 15-минутный отрезок с начала вступления P-волны, являющийся характерной продолжительностью всплеска dB/dt. Геомагнитный сигнал преимущественно занимает диапазон частот от 0.05 до 0.2 Гц; сейсмический сигнал тяготеет в область более низких частот (меньше 0.1 Гц). Примеры периодограмм, построенных по трем компонентам сейсмических колебаний и двум горизонтальным компонентам dB/dt, приведены на рис. 6.

 

Рис. 6. Примеры периодограмм по нормированным данным для трех компонент смещений (синий) и двух горизонтальных компонент скорости изменения магнитного поля (оранжевый), зарегистрированных за 15-минутный интервал с момента первого вступления: пары пунктов наблюдений GKP-HLP и MHVs-MHVg для события ЗТ-1 (а); KHC-BDV и ADZR-KLI для события ЗТ-2 (б); TASB-GLK и PUL-SPG для события ЗТ-3 (в).

 

Из рис. 6 видно, что спектральные составы сейсмического и электромагнитного сигналов в целом согласуются друг с другом, но имеют существенное различие в деталях. Поэтому нами были рассчитаны спектры отдельно для объемных P- и S-волн и поверхностных волн Рэлея. Предварительно была выполнена фильтрация сейсмических и электромагнитных наблюдений в диапазоне частот от 20 до 125 мГц фильтром Баттерворта 6 порядка. При выборе границ частотного интервала мы исходили из свойств сейсмического сигнала. При распространении объемных волн происходит рассеяние высокочастотной составляющей. На телесейсмических расстояниях продольные волны на частотах выше 1 Гц не превышают шума, а спектр поперечных волн от удаленных событий ограничен частотой порядка 0.2 Гц. Характерные частоты поверхностных волн еще ниже и составляют 0.025–0.05 Гц. Если рассматривать косейсмический сигнал как линейный электромагнитный отклик среды на сейсмическое возмущение, то подавленные фильтром возмущения следует рассматривать как помеху.

На рис. 7 приведены периодограммы трех пар сейсмических и геомагнитных наблюдений, построенных на отрезке прохождения поверхностных волн. На верхних панелях приведена амплитуда спектра Z-компоненты сейсмического сигнала, содержащего ряд спектральных линий рэлеевских мод. Аналогичную структуру имеет спектр R-компоненты (не приведен на рисунках). Согласно работе [Surkov et al., 2018], косейсмический эффект имеет электромагнитную природу и проявляется на X- и Z-компонентах. Наличие многочисленных обертонов говорит о достаточно высокой добротности эффективного резонатора под обсерваторией.

 

Рис. 7. Периодограммы Z-компоненты сейсмического сигнала при прохождении волны Релея (вверху) и трех компонент (X, Y, Z) электромагнитного отклика dB/dt на трех парах пунктов наблюдений: VJF-SPG (ЗТ-1) (а); MHVs-MHVg (ЗТ-3) (б); PUL-SPG (ЗТ-3) (в).

 

На рис. 8 приведен спектр модуля сейсмического сигнала объемных волн и три компоненты электромагнитного отклика. Высокочастотная область спектра соответствует продольной волне. Спектр на низких частотах определяется поперечной волной. Здесь основной косейсмический эффект является электрокинетическим и проявляется на Y-компоненте.

 

Рис. 8. Периодограммы модуля сейсмического сигнала |(N, Z)| при прохождении объемных P- и S-волн (вверху) и трех компонент (X, Y, Z) электромагнитного отклика dB/dt на паре пунктов наблюдений MHVs-MHVg для события ЗТ-2.

 

ЗАКЛЮЧЕНИЕ

Обнаружение квазипериодических сигналов с частотой 0.05–0.2 Гц на геомагнитных записях подтверждает возможность генерации низкочастотных электромагнитных полей под воздействием на среду под обсерваторией сейсмических волн. Для сильных землетрясений с магнитудой Mw ≥ 7 эффект проявляется на значительных эпицентральных расстояниях вплоть до 27° (около 3000 км по поверхности). Наилучшим образом сигнал детектируется в записях скорости изменения магнитного поля dB/dt в виде изолированных всплесков с последующим затуханием, что свидетельствует об индуцировании геоэлектрического поля E ~ dB/dt в поверхностных слоях Земли. Чем ближе геомагнитная обсерватория к эпицентру, тем всплеск имеет бо́льшую амплитуду и тем более он локализован во времени. На бо́льших расстояниях УНЧ волны затухают при распространении в литосфере [Molchanov, Hayakawa, 2008; Surkov et al., 2018], и сигнал в скорости изменения магнитного поля становится соизмерим с его шумовой составляющей. Наблюдаемая пиковая амплитуда |dX/dt|~ ~17 нТл/с (обсерватория GLK) является экстремальной для любых широт и едва достижима даже во время очень сильных магнитных бурь. Таким образом, из всех природных процессов, возбуждающих быстрые изменения геомагнитного поля, одними из самых мощных являются сейсмические события на небольшом расстоянии от пункта наблюдений. Очевидно, величина магнитного эффекта зависит от параметров среды. Скорость распространения электромагнитного сигнала определяется скоростью распространения сейсмических волн. На близкорасположенных друг к другу пунктах высокоточного геомагнитного и сейсмического мониторинга начало возмущения скорости изменения магнитного поля совпадает с приходом P-волны. Это говорит в пользу механизма генерации электромагнитного сигнала по мере распространения сейсмической волны (т.н. “косейсмический электромагнитный эффект”) [Surkov et al., 2018; Yamazaki, 2024] против механизма его передачи посредством акустико-гравитационных волн, воздействующих на ионосферу [Рябова, Шалимов, 2022; Hayakawa et al., 2010; Канониди, 2014]. Длительность геомагнитного сигнала увеличивается с расстоянием, что также согласуется с результатами численного моделирования косейсмического эффекта [Molchanov et al., 2001].

Таким образом, нами впервые представлены экспериментальные наблюдения, которые отражают данный эффект. Кроме фактической регистрации эффекта, была выявлена его тонкая структура – вклады, соответствующие разным типам сейсмических волн. Показано, что объемные волны возбуждают, главным образом, геомагнитные возмущения, возникающие в результате градиента механического давления. В статье [Surkov et al., 2018] такой эффект назван электрокинетическим, в отличие от электромагнитных возмущений, возбуждаемых поверхностной волной. Последние нами также были идентифицированы.

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

БЛАГОДАРНОСТИ

Авторы выражают благодарность ФИЦКИА УрО РАН за доступ к данным УНУ “Архангельская сейсмическая сеть” и Международной федерации цифровых сейсмических сетей (International Federation of Digital Seismograph Networks, https://www.fdsn.org/) за предоставление доступа к данным зарубежных сейсмических станций.

В исследовании использовались данные, регистрируемые геомагнитными обсерваториями. Мы благодарим национальные институты, которые поддерживают обсерватории, сеть ИНТЕРМАГНЕТ за продвижение высоких стандартов данных (http://www.intermagnet.org) и Межрегиональный центр геомагнитных данных (http://geomag.gcras.ru) за свободное распространение данных в режиме онлайн. В работе также использовались данные и сервисы ЦКП “Аналитический центр геомагнитных данных” Геофизического центра РАН (https://ckp.gcras.ru/) и УНУ “Среднеширотный комплекс геофизических наблюдений “Михнево” Института динамики геосфер им. академика М.А. Садовского РАН (https://idg.ras.ru/tsentry/unu-skgn-mikhnevo/). Авторы признательны двум рецензентам за полезные замечания.

ФИНАНСИРОВАНИЕ РАБОТЫ

Исследование выполнено в рамках государственных заданий Геофизического центра РАН, Института физики Земли им. О.Ю. Шмидта РАН и Института динамики геосфер им. академика М.А. Садовского РАН (№ 122040400015-5), утвержденных Министерством науки и высшего образования Российской Федерации.

Об авторах

А. А. Соловьев

Геофизический центр РАН; Институт физики Земли им. О.Ю. Шмидта РАН

Автор, ответственный за переписку.
Email: a.soloviev@gcras.ru
Россия, г. Москва; г. Москва

И. М. Алешин

Геофизический центр РАН; Институт физики Земли им. О.Ю. Шмидта РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва; г. Москва

С. В. Анисимов

Институт физики Земли им. О.Ю. Шмидта РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва

А. Г. Гоев

Институт физики Земли им. О.Ю. Шмидта РАН; Институт динамики геосфер им. академика М.А. Садовского РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва; г. Москва

А. Н. Морозов

Геофизический центр РАН; Институт физики Земли им. О.Ю. Шмидта РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва; г. Москва

Д. С. Сапронов

Геофизический центр РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва

Е. Н. Соловьева

Геофизический центр РАН; Институт физики Земли им. О.Ю. Шмидта РАН

Email: a.soloviev@gcras.ru
Россия, г. Москва; г. Москва

Список литературы

Дополнительные файлы