Вестник КРАУНЦ. Физ.-мат. науки. 2020. Т. 33. № 4. C. 157-165. ISSN 2079-6641
ФИЗИКА
УДК 551.510.535
Научная статья
Исследование влияния сейсмической активности на коэффициент полупрозрачности спорадического слоя Es над
^мчаткой
Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, с. Паратунка, ул. Мирная, 7 E-mail: vbogd@ikir.ru, pavlov@ikir.ru
Мелкомасштабная турбулизация спорадического слоя Es (десятки метров) оценивается коэффициентом полупрозрачности AfbEs, а крупномасштабная турбулизация (сотни метров) выражается в Es-рассеянии. Перед крупными землетрясениями с магнитудами M > 5.0 может наблюдаться формирование аномальных спорадических слоев в области E: увеличение электронной концентрации и коэффициента полупрозрачности на интервале от 14 до 3 суток до момента землетрясения и последующее резкое уменьшение значений этих параметров за 1 сутки до землетрясения. В данной работе проведён анализ изменений коэффициента полупрозрачности спорадического слоя Es, предшествующих наступлению землетрясений с магнитудой M > 5.0, произошедших в Камчатском регионе за период 2015-2019 гг. Произведена оценка прогностической эффективности параметра AfbEs по методикам А. А. Гусева и Г. М. Молчана.
Ключевые слова: ионосфера, литосфера, вертикальное радиозондирование, землетрясение, предвестник землетрясений.
Для цитирования. Богданов В. В., Павлов А. В. Исследование влияния сейсмической активности на коэффициент полупрозрачности спорадического слоя Es над Камчаткой // Вестник КРАУНЦ. Физ.-мат. науки. 2020. Т. 33. № 4. C. 157-165. DOI: 10.26117/20796641-2020-33-4-157-165
Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.Org/licenses/by/4.0/deed.ru)
Богданов В. В., Павлов А. В, 2020
Введение
Область Е ионосферы расположена на высотах от 90 до 160 км. На ионограммах в области E могут наблюдаться два слоя: регулярный, создающийся солнечным излучением и нерегулярный спорадический Es-слой. Спорадические слои Es - это
Финансирование. Работа выполнялась в рамках государственного задания по теме «Динамика физических процессов в активных зонах ближнего космоса и геосфер» (Рег.№ АААА-А17-117080110043-4).
В. В. Богданов, А. В. Павлов
DOI: 10.26117/2079-6641-2020-33-4-157-165
Поступила в редакцию: 21.10.2020
В окончательном варианте: 23.11.2020
облака повышенной ионизации, расположеные преимущественно на высотах: 96, 103-107, 111-112, 117-118 и 128 км [1], в которых преобладают ионы металлов метеорного происхождения (Ре+, Ыа+ и др.) [2]. Толщина спорадических
слоев может изменяться от нескольких сотен метров до нескольких километров, а их горизонтальные размеры могут составлять несколько сотен километров.
Механизм формирования стандартных спорадических слоев связан с явлением ветрового сдвига [3], а процесс их расплывания происходит за счет турбулентной и амбиполярной диффузии [4]. Кроме стандартных спорадических слоев существуют также аномальные спорадические слои, образование которых вызвано модификацией атмосферного электрического поля. Аномальные спорадические слои возникают когда электрическое поле на поверхности Земли направлено сверху вниз. Образование аномальных спорадических слоёв, в отличие от стандартных, происходит на высоте 120 км [5].
Исследования сейсмоионосферных связей показывают, что акустические и электромагнитные возмущения, возникающие в процессе подготовки землетрясений, могут быть существенной причиной турбулизации ионосферы в диапазоне масштабов от десятков метров до десятков километров и характерных времен существования порядка нескольких десятков минут. Одним из ионосферных параметров, характеризующих нерегулярные турбулентные процессы, является коэффициент полупрозрачности спорадического слоя. Аномальное изменение значения коэффициент полупрозрачности спорадического слоя Еэ в период подготовки сильных землетрясений наблюдались в различных сейсмоактивных регионах Земли [6, 7].
В данной работе будет произведена оценка прогностической эффективности коэффициента полупрозрачности спорадического слоя Еэ для землетрясений с магнитудой М > 5.0, произошедших в Камчатском регионе за период 01.08.201531.12.2019 гг.
Методика исследования
В данной работе использованы данные радиофизических наблюдений, полученные с помощью автоматической ионосферной станции (АИС) вертикального радиозондирования, расположенной в с. Паратунка (ф = 52.97° с.ш., X = 158.24° в.д.). Наблюдения проводились один раз в 15 минут в импульсном режиме на частотах от 1 до 15 МГц. Для анализа использован каталог землетрясений Камчатского филиала Федерального исследовательского центра "Единая геофизическая служба РАН". Ионосферный слой Еэ характеризуется двумя частотами зондирующих радиоволн - частотой экранирования fbEs, которая соответствует максимальной плотности ионизации в спорадическом слое (fbEs ~ V^тах), и предельной частотой отражения foEs, превышающей fbEs. Частотные параметры спорадического слоя Еэ информируют о структуре и содержании концентрации электронов в Еэ. Коэффициент полупрозрачности AfbEs характеризует рассеяние радиоволн на мелкомасштабных (десятки метров) неоднородностях Еэ и определяется по формуле [8]:
А^^ = - fbEs)/fbEs (1)
где foEs, fbEs - соответственно предельная частота и частота экранирования спорадического слоя в (МГц). Увеличение значений коэффициента полупрозрачности
AfbEs свидетельствует об образовании неоднородностей спорадических слоев. Увеличение степени неоднородности свидетельствует о возникновение турбулентности в спорадическом слое Еб. На рис. 1 представлены ионограммы, содержащие спорадический слой Еб при значениях коэффициента полупрозрачности А^Ез = 0 (рис. 1а) и А^Ез = 5.54 (рис. 1б).
5 МГц
а) б)
Рис. 1. Ионограммы с отмеченным спорадическим слоем Еб при различных значений коэффициента полупрозрачности AfbEs: а) А^Ез = 0; б) А^Ез = 5.54.
На рис. 2 представлены часовые значения коэффициента полупрозрачности AfbEs за период 01.08.2015-31.12.2019 гг.
Рис. 2. Временной ряд значения коэффициента полупрозрачности AfbEs за период 01.08.2015-31.12.2019 гг.
Для выделения аномальных возмущений в ионосфере в период подготовки землетрясений расчет коэффициента полупрозрачности Еб-слоя AfbEs выполнялся для значений foEs > fbEs по формуле (1), на следующем этапе полученные значения коэффициента полупрозрачности осреднялись за сутки местного времени по формуле
1 к
А^Езср = - £ (2)
к1=1
где к - количество значений коэффициента полупрозрачности Еб-слоя за сутки. На рис. 3 представлены значения коэффициента полупрозрачности AfbEs за период 01.08.2015-31.12.2019 гг., осредненные за сутки местного времени.
В качестве примера, на рис. 4 представлены значения А^Езср за период 05.11.2015-21.11.2015 гг., предшествовавшие землетрясениям с М > 5.0. Как видно из рисунка, наступлению землетрясений с магнитудами М > 5.0 предшествовало увеличение значений коэффициента полупрозрачности А^Езср.
Рис. 3. Временной ряд средних за сутки значения коэффициента полупрозрачности AfbEs за период 01.08.2015-31.12.2019 гг.
Рис. 4. Значения коэффициента полупрозрачности АfbEsCp за период 05.11.2015— 21.11.2015 гг. перед наступлением сейсмических событий: 13.11.2015 г. (M = 5.4, h = 51 км), 16.11.2015 г. (M = 5.1, h = 50 км) и 19.11.2015 г. (M = 5.5, h = 51 км).
Для оценки прогностической эффективности параметра АfbEs^ использовались следующие характеристики: надёжность R и достоверность V [9], а также эффективности Jq и Jm, вычисленные соответственно по методикам А.А. Гусева [10] и Г.М. Молчана [11].
Надежность предвестника R определяется как отношение количества землетрясений N+, для которых был выделен предвестник, к числу всех землетрясений N:
R = N+ (3)
N V '
Достоверность предвестника V определяется как отношение числа предвестниковых аномалий п(Ае) к общему числу выделенных аномалий n(A):
П(Ае )
V = (4)
n(A)
Эффективность прогностического признака на основе методики А.А. Гусева вычисляется по формуле [10]:
JQ = ^ (5)
q N/Т
где Т - общее время мониторинга сейсмической обстановки, N+ - количество землетрясений, соответствующих успешному прогнозу за время Т, N - общее количество произошедших землетрясений, Ттр - общее время тревоги (суммарная длительность всех промежутков времени, в которых действовал прогноз по
оцениваемому методу в течение общего времени мониторинга). Статистическая значимость а определяется вероятностью получения значений эффективности /о в отсутствие связи "землетрясение-предвестник" исходя из гипергеометрического распределения
Согласно характеристике информативности прогностических признаков [12], при /о < 1.2 прогностический признак считается неинформативным и его использование не приведёт к заметному улучшению прогноза, при 1.2 < /о < 1.5 прогностический признак может оказаться эффективным в сочетании с другими признаками, а при /о > 1.5 прогностический признак считается информативным и может быть использован в алгоритмах прогноза землетрясений.
Эффективность прогностического признака на основе методики Г. М. Молчана вычисляется по формуле [11]:
/м = 1 - V - т (6)
где V = 1 — - доля пропусков цели, т = - относительный объём тревог.
Для случайного прогноза /м = 0, а для идеального (без пропуска цели и с нулевым временем тревоги) - /м = 1. На диаграмме Молчана диагональ т + V = 1, соединяющая точки (0; 1) и (1;0) соответствует "случайному" прогнозу. Для этой диагонали можно построить доверительный интервал с заданной значимостью а. Если экспериментальные точки (т, V) лежат под нижней границей доверительного интервала с заданной значимостью а, то это интерпретируется как высокая степень надежности выявленной связи рассматриваемого предвестника с прогнозируемыми землетрясениями.
Оценка прогностической эффективности рассматриваемого параметра А^Е,уср производилась для четырех периодов ожидания Гож равных 3, 5, 7 и 14 суток. В качестве прогнозируемых землетрясений рассматривались сейсмические события с магнитудой М > 5.0 и глубиной гипоцентра к < 100 км, произошедших на эпицентральных расстояниях г < 300 км от пункта ионосферных наблюдений в с. Паратунка. Период ожидания длительностью Гож объявлялся при условии превышения значений параметра AfbEsCр значения М(А^Е,ср) + о(А^Е,ср) = 0.6 , где М(А^Е,ср) и о(А^Е,ср) соответственно математическое ожидание и среднеквадратическое отклонение параметра А^Е,ср.
В таблице и на рис. 5 представлены результаты оценки прогностической эффективности коэффициента полупрозрачности слоя Еэ.
При увеличении значений периода ожидания землетрясения Гож возрастают значения параметров надёжности Я и достоверности V прогностического признака, но снижаются значения эффективности /о и /м. Для периода ожидания Тож = 3 суток эффективность прогноза /о = 1.78, т.е прогностический признак можно считать информативным. Точка прогноза (т, V) для Тож = 3 суток лежит ниже 99% доверительного интервала (рис. 5а), что можно интерпретировать как достаточно высокую степень надёжности выявленной связи рассматриваемого предвестника с землетрясениями данного диапазона магнитуд. Для периодов ожидания Тож равных 5, 7 и 14 суток эффективность прогноза /о снижается соответственно с 1.47 до 1.25. При таких значениях /о прогностический признак считается не очень полезным, но для увеличения эффективности прогноза землетрясений может быть использован совместно с другими предвестниками. Точка прогноза (т, V) при Тож = 5 суток лежит на границе 99% доверительного интервала (рис. 5б), а для Тож = 7 суток и Тож = 14 суток - на границы 95% доверительного интервала (рис. 5в,г).
в)
г)
Рис. 5. Диаграммы ошибок для коэффициента полупрозрачности А^Езср для периодов ожидания: а) Тож = 3; б) Тож = 5; в) Тож = 7 суток; г) Тож = 14 суток.
При значениях длительности периода ожидания Тож равных 5, 7 и 14 суток связь рассматриваемого ионосферных предвестника с землетрясениями с магнитудами M > 5.0 не возможно однозначно интерпретировать как надёжную.
Таблица
Прогностическая эффективность коэффициента полупрозрачности А/ЬЕяср
для землетрясений с M > 5.0.
Тож (сутки) 3 5 5 14
N+ 20 24 26 36
N 59 59 59 59
n(Äe ) + 19 22 24 29
n(A) 115 115 115 115
Ттр (сутки) 307 449 554 787
Т (сутки) 1614 1614 1614 1614
V 0.17 0.19 0.21 0.25
R 0.34 0.41 0.44 0.61
Jg 1.78 1.46 1.28 1.25
а 0.01 0.04 0.12 0.06
т 0.19 0.28 0.34 0.49
V 0.66 0.59 0.56 0.39
Jm 0.15 0.13 0.1 0.12
Заключение
Проведён анализ изменений коэффициента полупрозрачности спорадического слоя Еэ, предшествующих наступлению землетрясений с магнитудой М > 5.0, произошедших в Камчатском регионе за период 01.08.2015-31.12.2019 гг. Произведенная оценка прогностической эффективности усреднённого за сутки местного времени параметра А^Е, для периодов ожидания землетрясений различной длительности показала, что наибольшая эффективность /о = 1.78 и /м = 0.15 достигается при Тож = 3 суток. При этом появление аномальных значений прогностического параметра предшествовало 34% землетрясений рассматриваемого диапазона магнитуд (надежность предвестника Я = 0.34). При значении периода ожидания Тож = 3 суток прогностический признак А^Е,ср можно считать информативным. Увеличение длительности периода ожидания землетрясения Тож приводит к снижению эффективности прогноза по данному прогностическому признаку.
Конкурирующие интересы. Авторы заявляют, что конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответсвенность. Все авторы участвовали в написании статьи и полностью несут ответственность за предоставление окончательной версии статьи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Список литература/References
[1] Брюнелли Б.Е., Намгаладзе А.А., Физика ионосферы, Наука, М., 1988, 528 с. [Bryunelli B.E., Namgaladze A.A., Fizika ionosfery, Nauka, M., 1988 (in Russian), 528 pp.]
[2] Liperovsky V.A. et al., "Modification of sporadic E-Layers caused by seismic activity", Surveys in geophysics, 21 (2000), 449-486
[3] Гершман Б.Н., Игнатьев Ю.А., Каменецкая Г.Х., Механизмы образования ионосферного спорадического слоя Es на различных широтах, Наука, М., 1976, 108 с. [Gershman B.N., Ignat'ev Yu.A., Kamenetskaya G.Kh., Mekhanizmy obrazovaniya ionos-fernogo sporadicheskogo sloya Es na razlichnykh shirotakh, Nauka, M., 1976 (in Russian), 108 pp.]
[4] Gurevich A.V., Borisov N.D., Zypin K.P., "Ionospheric turbulence induced in the lower part of the E region by the turbulence of the neutral atmosphere", JGR Space Physics, 102 (1997), 379-388
[5] Пулинец С.А. и др., "Атмосферное электрическое поле как источник изменчивости ионосферы", Успехи физических наук, 168:5 (1998), 582-589. [Pulinets S.A. et al., "At-mosfernoe elektricheskoe pole kak istochnik izmenchivosti ionosfery", Uspekhi fizicheskikh nauk, 168:5 (1998), 582-589 (in Russian)].
[6] Liperovskaya E.V. et al., "Variability of sporadic E-layer semi transparency (foEs - fbEs) with magnitude and distance from earthquake epicenters to vertical sounding stations", Natural Hazards and Earth System Sciences, 3 (2003), 279-284.
[7] Liperovsky V.A. et al., "On spread - Es effects in the ionosphere before earthquakes", Natural Hazards and Earth System Sciences, 5 (2005), 59-62.
[8] Липеровский В.А. и др., "Физические модели связей в системе литосфера-атмосфера-ионосфера перед землетрясениями", Геомагнетизм и аэрономия, 48:6 (2008), 831843. [Liperovskiy V.A. et al., "Fizicheskie modeli svyazey v sisteme litosfera-atmosfera-ionosfera pered zemletryaseniyami", Geomagnetizm i aeronomiya, 48:6 (2008), 831-843 (in Russian)].
[9] Салтыков В.А., "О возможности использования приливной модуляции сейсмических шумов в целях прогноза землетрясений", Физика Земли, 2017, №2, 84-96.
[Saltykov V.A., "O vozmozhnosti ispol'zovaniya prilivnoy modulyatsii seysmicheskikh shu-mov v tselyakh prognoza zemletryaseniy", Fizika Zemli, 2017, №2, 84-96 (in Russian)].
[10] Гусев А.А., "Прогноз землетрясений по статистике сейсмичности", Сейсмичность и сейсмический прогноз, свойства верхней мантии и их связь с вулканизмом Камчатки, Наука, Новосибирск, 1974, 109-119. [Gusev A.A., "Prognoz zemletryaseniy po statistike seysmichnosti", Seysmichnost' i seysmicheskie prognoz, svoystva verkhney mantii i ikh svyaz's vulkanizmom Kamchatki, Nauka, Novosibirsk, 1974, 109-119 (in Russian)].
[11] Molchan G.M., "Strategies in strong earthquake prediction", Physics of the Earth and Planetary Interiors, 61 (1990), 84-98.
[12] Завьялов А.Д., Среднесрочный прогноз землетрясений: основы, методика, реализация, Наука, М., 2006, 254 с. [Zav'yalov A.D., Srednesrochnyy prognoz zemletryaseniy: osnovy, metodika, realizatsiya, Nauka, M., 2006 (in Russian), 254 pp.]
Vestnik KRAUNC. Fiz.-Mat. Nauki. 2020. vol. 33. no. 4. P. 157-165. TSSN 2079-6641
PHYSICS
MSC 86A10, 86A15 Research Article
Investigation of the influence of seismic activity on the semi transparency coefficient of the sporadic Es layer over Kamchatka
V. V. Bogdanov, A. V. Pavlov
Institute of Cosmophysical Researches and Radio Wave Propagation Far-Eastern Branch Russian Academy of Sciences, 684034, Kamchatskiy Kray, Paratunka, Mirnaya st., 7 E-mail: vbogd@ikir.ru, pavlov@ikir.ru
Small-scale turbulization of the sporadic Es layer (tens of meters) is estimated by the semi transparency coefficient AfbEs, and large-scale turbulization (hundreds of meters) is expressed in Es scattering. It is known that before large earthquakes with magnitudes M > 5.0, the formation of abnormal sporadic layers in the E region can be observed: an increase in the electron concentration and semi transparency coefficient at the interval from 14 to 3 days before the earthquake and a subsequent sharp decrease in the values of these parameters 1 day before the earthquake. This paper analyzes changes in the semi transparency coefficient of the sporadic Es layer preceding the onset of earthquakes with a magnitude M > 5.0 that occurred in the Kamchatka region during the period 2015-2019. The prognostic efficiency of the parameter AfbEs was evaluated using the methods of A. A. Gusev and G. M. Molchan.
Key words: ionosphere, lithosphere, vertical radio sounding, earthquake, earthquake precursor.
DOT: 10.26117/2079-6641-2020-33-4-157-165
Original article submitted: 21.10.2020 Revision submitted: 23.11.2020
For citation. Bogdanov V. V., Pavlov A. V. Investigation of the influence of seismic activity on the semi transparency coefficient of the sporadic Es layer over Kamchatka. Vestnik KRAUNC. Fiz.-mat. nauki. 2020, 33: 4,157-165. DOT: 10.26117/2079-6641-2020-33-4-157-165
Competing interests. The authors declare that there are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. All authors contributed to this article. Authors are solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by all authors.
The content is published under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/deed.ru)
© Bogdanov V. V., Pavlov A.V., 2020
Funding. The work was carried out within the framework of the state assignment on the topic "Dynamics of physical processes in active zones of near space and geospheres" (Reg. No. AAAA-A17-117080110043-4).