Исследование сейсмичности в районе Нанкайского глубоководного желоба методом LURR с целью ретроспективной оценки прогноза землетрясений

А.С. Закупин Д.В. Костылев Н.В. Костылева

Журнал: Вестник геонаук @vestnik-geo

Рубрика: Научные статьи

Статья в выпуске: 9 (381), 2026 года.

Бесплатный доступ

Нанкайский глубоководный желоб представляет собой геологический объект с высокой сейсмической опасностью, особенно в южной части, в районе японского острова Кюсю, где в 2016 и 2024 годах. произошли землетрясения магнитудой более 7.0. Целью исследования являлся ретроспективный прогноз сильных землетрясений в южном сегменте Нанкайского глубоководного желоба за период 2010–2025 годы. с использованием метода среднесрочного прогноза Load/Unload Response Ratio (LURR). Расчёты параметра LURR выполнены в девяти круговых областях, равномерно покрывающих очаговые зоны землетрясений магнитудой 7.3 в 2016 году и 7.1 в 2024 году. Выявленные аномалии LURR сопоставлялись во временной и пространственной областях с целевыми землетрясениями магнитудой более 6.0. Установлено, что близко расположенные сейсмические очаги характеризуются высокой чувствительностью к землетрясениям M > 6.0, вызывающим резкие изменения параметра LURR. Такие наведённые аномалии составили до половины всех ложных признаков. После исключения этих событий из статистического анализа результаты продемонстрировали приемлемую эффективность метода для оценки сейсмической опасности исследуемого района: все восемь землетрясений M > 6.0 имели прогностические признаки, составившие 72 % от общего числа выявленных аномалий. Средняя продолжительность периода тревоги составила около шести месяцев. Полученные результаты подтверждают перспективность применения метода LURR для среднесрочного прогноза сильных землетрясений в зонах субдукции.

сейсмичность \ землетрясение \ LURR \ аномалия \ прогноз

Короткий адрес: https://sciup.org/149152353

IDS: 149152353   |   УДК: 550.343.6   |   DOI: 10.19110/geov.2026.9.4

Study of seismicity in the Nankai trough area using the LURR method for retrospective evaluation of earthquake prediction

The Nankai Trough is a geological feature characterized by a high level of seismic hazard, particularly in its southern part near the Japanese island of Kyushu, where earthquakes with magnitudes greater than 7.0 occurred in 2016 and 2024. The aim of this study is to conduct a retrospective forecast of strong earthquakes in the southern segment of the Nankai Trough during the period 2010–2025 using the Load/Unload Response Ratio (LURR) method for medium-term earthquake prediction. LURR calculations were performed in nine circular areas uniformly covering the source zones of earthquakes with magnitudes of 7.3 in 2016 and 7.1 in 2024. The detected LURR anomalies were compared with target earthquakes with magnitudes greater than 6.0 in both the temporal and spatial domains. The results showed that closely located seismic sources exhibit high sensitivity to earthquakes with M > 6.0, causing sharp changes in the LURR parameter. These induced anomalies accounted for up to half of all false precursors. After excluding these events from the statistical analysis, the results demonstrated an acceptable efficiency of the method for assessing seismic hazard in the study area: all eight earthquakes with M > 6.0 were preceded by precursory signals, which accounted for 72 % of the total number of detected anomalies. The average duration of the alarm period was approximately six months. The results confirm the potential of the LURR method for medium-term forecasting of strong earthquakes in subduction zones.

Текст научной статьи Исследование сейсмичности в районе Нанкайского глубоководного желоба методом LURR с целью ретроспективной оценки прогноза землетрясений

Метод среднесрочного прогноза землетрясений LURR известен с 1991 года, когда появились первые публикации китайских ученых (Yin X., Yin C., 1991).

Последние работы в данном направлении были сфокусированы на разработке метода ретроспективных прогнозов землетрясений (Yin X.C. et al., 2006), оползней (Zhang W.J. et al., 2006), а также экспериментах на образцах горных пород для изучения их физико-

механических свойств (Gong F.-Q. et al., 2006; Zhu F.C. et al., 2016; Zhao F., He M.C., 2017).

Одно из новейших исследований посвящено ком-плексированию отклика LURR в деформационном и ге-оэлектрическом полях (Yu C. et al., 2022). Результаты многочисленных работ рассматривались в целом позитивно с точки зрения прогноза, и это, несомненно, вызвало интерес и за пределами Китая. Так, с 2016 года метод LURR активно применяется в России, в Институте морской геологии и геофизики Дальневосточного отделения РАН (ИМГиГ ДВО РАН). Специально разработанное программное обеспечение (Zakupin, 2016) позволило проводить многочисленные расчёты для тестирования методики на различных сейсмических каталогах и оптимизировать настройки для получения лучших результатов.

В 2016 и 2017 годах удалось применить полученные наработки в реальном прогнозе. Прогнозы рассматривались на заседаниях Сахалинского отделения Российского экспертного совета по чрезвычайным ситуациям (СФ РЭС) (протокол № 3 от 11.05.2016 г., протокол № 2 от 16.03.2017 г.). Оба прогноза были признаны полностью реализованными после Онорского (Мw = 5.8, 14.08.2016 г.) и Крильонского (Мw = 5.0, 23.03.2017 г.) землетрясений (Zakupin et al., 2018; Закупин, Семенова, 2018). При этом в активе сахалинских исследователей оказались многочисленные работы по ретроспективному прогнозу землетрясений в различных сейсмоактивных регионах мира (Закупин, 2026).

Продолжением исследований стал запуск эксперимента по мониторингу территории о. Сахалин в режиме реального времени для своевременного выявления аномалий и подготовки прогнозов. Эксперимент был запущен в 2022 году ИМГиГ ДВО РАН и Сахалинским филиалом Федерального исследовательского центра «Единая геофизическая служба РАН» (СФ ФИЦ ЕГС РАН). Подобным в практике российских исследований является продолжительный эксперимент в реальном времени с методами M8 и MSc (Kossobokov, 2011). Как утверждает автор, данные мировой сети обрабатывались каждые 6 месяцев, а результаты в режиме закрытого доступа публиковались на внутреннем сайте Института теории прогноза землетрясений и математической геофизики РАН. При этом даже при относительно скромной статистике целевых землетрясений (17 землетрясений с магнитудами >8.0) с 1992 по 2010 год в работе отсутствуют данные по всем землетрясениям, а приведены лишь некоторые примеры. Судя по итоговым результатам данного моделирования, представленных позже (Kossobokov, 2013), последними в цепочке успешных прогнозов стали землетрясения в Чили (2010 г.) и в Японии (2011 г.).

Известно, что данное исследование продолжается и в настоящее время. Эффективность эксперимента на о. Сахалин к 2026 году выражается в шести оформленных прогнозах, четырех реализованных, одном неуспешном и одном действующем. Результаты проведенной работы опубликованы в ряде статей (Закупин и др., 2023; Закупин, 2025). При этом авторы отмечают, что основной проблемой остаётся процент ложных тревог, который неизбежно будет сохраняться ввиду низкого уровня фоновой сейсмической активности в рабочем диапазоне магнитуд (от 3.3 до 5). В регионах с высокой сейсмической активностью, таких как

Япония или Турция, результаты ретроспективных исследований оказались лучше, поэтому возникает желание продолжать апробацию метода на таких каталогах, пусть и в ретроспективе. При выборе очередной цели для исследования принималось во внимание не только качество данных и уровень активности в регионе, но и наличие фактора ожидания сильных землетрясений по данным долгосрочных оценок. Таковой стала зона на стыке Филлипинской и Евразийской плит в районе о. Кюсю и о. Хонсю (Япония), где одним из интереснейших объектов является Нанкайский глубоководный желоб.

Нанкайский желоб — это подводная впадина, расположенная к югу от региона Нанкайдо на о. Хонсю примерно на 900 км от берега. Подстилающий её разлом Нанкай является источником разрушительных землетрясений. Вдоль Нанкайской впадины на юго-западе Японии Филиппинская морская плита субду-цирует под Евразийскую плиту, и здесь неоднократно происходили крупные межплитные землетрясения. Последним крупным надвиговым землетрясением в центральной части Нанкайской впадины у полуострова Кии было землетрясение Тонанкай 1944 года (Mw = 8.2) (Tanioka et al., 2004). Сейсмичность в зоне субдукции, вероятно, связана с балансом напряжений во входящей и субдуцирующей плитах, а также со взаимодействием между субдуцирующей и перекрывающей плитами (Obana et al., 2005). Исторические данные о крупных межплитных землетрясениях, зарегистрированных с VII века, указывают на интервал повторяемости в 100–200 лет. Эти крупные землетрясения представляли собой межплитные надвиговые землетрясения, вызванные субдукцией Филиппинской морской плиты под Юго-Западную Японию.

В августе 2024 года, после землетрясения с магнитудой Mw = 7.1, произошедшего у восточного побережья острова Кюсю, Японское метеорологическое агентство (Japan Meteorological Agency — JMA) заявило, что вероятность крупномасштабного землетрясения в районе разлома Нанкай выше обычного (URL: https://www. , дата обращения: 27.04.2026). Агентство отметило, что нельзя определить конкретно период, когда может случиться землетрясение, тем не менее предупреждение о возможном катаклизме было разослано в 707 муниципалитетов в 29 префектурах, которые могут пострадать в случае сильных подземных толчков. В этом же месяце Япония отменила состояние повышенной готовности к сильному землетрясению в районе Нанкай-ского жёлоба, но сохранила прогноз вероятности сильнейшего землетрясения в районе разлома Нанкай на ближайшие 30 лет.

Целью работы является изучение и оценка сейсмической опасности в районе Нанкайского глубоководного желоба методом LURR.

Данные и методы

Метод LURR (Yin X., Yin C., 1991; Yin X.C. et al., 1995; Yin X.C. et al., 2006) хорошо описан в монографии (Закупин, 2026), а потому сфокусируемся на нашей практике его использования в методологическом плане. Проводимые расчёты методом LURR позволяют делать среднесрочные оценки сейсмической опас- 39

ности. По ранее полученным данным для о. Сахалин и ряда других регионов, более 85 % случаев успешных ретроспективных прогнозов реализовывались в период до 2 лет с момента появления аномалий и определения зоны прогноза (Закупин, 2026). В этой связи традиционно используем данный период для нашего исследования.

В методике расчетов параметра LURR важную роль играет выбор подхода для определения расчетных областей, а также дискриминация аномалий и формирование зон прогноза. И в данном случае все эти принципы полностью соответствуют тем, что определены на каталоге достаточной глубины (за период 1997– 2019 гг.) сейсмического режима о. Сахалин (Закупин, 2026). Отметим, что в качестве базового параметра для расчета LURR используется корень из сейсмической энергии (деформация Беньоффа), как и в большинстве расчётов, проводимых китайскими коллегами (Yin X.C. et al., 1995).

При вычислении параметра используется покрытие территории расчетными областями (окружности радиусом в один градус с шагом через полградуса). Параметры для математической обработки — скользящее окно в 360 суток и шаг в 5 дней. Итоговая рабочая выборка землетрясений в круговой области расчета дополнительно ограничена в диапазоне магнитуд от 3.3 до 5.0 («магнитудный cutoff», согласно Yin X.C. et al., 1995). Все перечисленные параметры обработки ранее прошли апробацию в многочисленных экспериментальных исследованиях и доказали свою состоя- тельность (Закупин, 2026) не только для сахалинских каталогов. В качестве инструмента для выделения аномалий применяется порог дискриминации, который равен 3 и при превышении которого параметр находится в области аномальных значений. Если в период тревоги в зоне прогноза происходит землетрясение с заявленной магнитудой (M ≥ 5.0), то прогноз считается реализовавшимся и наблюдение завершается.

Исходным материалом для проведенных исследований послужили каталоги JMA с 2010 по 2025 год включительно (URL: php?LANG=en, дата обращения: 15.05.2026). Метод LURR в первоисточнике (Yin X., Yin C., 1991) связывает появление аномалии с подготовкой в расчетной зоне землетрясения с магнитудой выше верхней границы расчетной выборки. При этом нижняя граница расчетной выборки устанавливается произвольно для отсечения слабой (фоновой) сейсмичности. Верхняя же граница показывает минимальную магнитуду ожидаемого землетрясения (главного толчка). Оба этих порога подбираются с учетом характера региональной сейсмичности. Нами экспериментально установлено, что приемлемым вариантом при расчетах нижней границы может быть выбрана М = 3.5. Выборка землетрясений расположена в координатах 30–34° северной широты и 130–134° восточной долготы. В данную область попадают 9 расчетных областей, показанных на рис. 1 в виде пронумерованных кругов. Всего количество землетрясений, попавших в выборку, составило 1056, из них с магнитудой М > 6.0 — 10.

Рис. 1. Карта района исследования: красные круги — расчетные области, желтыми маркерами показаны землетрясения с М ≥ 3.5, синими маркерами — с М > 6.0 за период 2010–2025 гг. (после декластеризации)

Fig. 1. Map of the study area: red circles represent the calculated areas, yellow markers show earthquakes with M ≥ 3.5, and blue markers show earthquakes with M > 6.0 for the period 2010–2025 (after declustering)

Для определения целей прогноза необходимо было провести декластеризацию. Декластеризация выполнена с использованием кластерного метода (Reasen-berg, 1985), позволяющего корректировать различные параметры идентификации афтершоков путем связывания землетрясений с кластерами в соответствии с пространственными и временными зонами взаимодействия. Использовался программный комплекс ZMAP (Wiemer, 2001), в котором данный алгоритм реализован. После декластеризации исходного каталога в нём осталось 873 события, из них с магнитудой М > 6.0 — 8. Эти землетрясения также показаны на карте (рис. 1) в виде синих кружков.

Результаты и обсуждение

Перед тем как приступить к работе, в произвольных ячейках необходимо выполнить расчёты в круговых областях (радиус, как и во всех наших работах, равен одному градусу) с центрами в эпицентрах наиболее сильных землетрясений. Таких у нас два, в 2016 и в 2024 годах (рис. 2, a, b).

Как можно видеть на графиках, в эпицентре первого землетрясения (15.04.2016 г., М = 7.3) параметр LURR 4 раз выходит в область аномальных значений (рис. 2, a), причём аномалия в 2015 году появилась за 10 месяцев до землетрясения и явилась его предвестником. Две аномалии параметра в эпицентре второго землетрясения (08.08.2024 г., М = 7.1) в июле 2016 года и декабре 2021 года (рис. 2, b) не являются предвестниками данного землетрясения. В итоге предвестник зафиксирован только для землетрясения 2016 года.

Результаты расчётов в девяти произвольных областях, равномерно покрывающих исследуемую область, должны показать прежде всего минимальные отклонения от эпицентральных расчётов, то есть выделенные на рис. 2, a, b аномалии должны повторяться хотя бы в некоторых из произвольных расчётных областей. Результаты расчётов представлены в таблице 2, а графики — на рис. 2, c–i. В двух из девяти расчётных областей количество данных оказалось недостаточным для проведения расчётов (области 3 и 6). В последнем столбце таблицы 2 (информация по аномалиям внутри отдельных расчётных областей) зеленым цветом отмечены аномалии, которые подходят в качестве предвестников (произошли в период до двух лет перед землетрясением) для тех из восьми землетрясений из таблицы 1 (13.03.2014 г., 15.04.2016 г., 09.05.2019 г., 21.01.2022 г., 17.04.2024 г., 08.08.2024 г., 13.01.2025 г., 02.04.2025 г.), которые входят в расчётную область и указаны в столбце 3 таблицы 2.

Фиолетовым цветом отмечены аномалии, произошедшие перед землетрясениями, которые не попали в область расчёта, но находятся в зоне нашего исследования (из списка в таблице 1) или имеют значительное влияние на региональную сейсмичность (как, например, землетрясение Тохоку в марте 2011 года). Красным цветом отмечены ложные аномалии, не име-

Рис. 2. Графики LURR-параметра для эпицентральных областей землетрясений 15.04.2016 г. (a) и 08.08.2024 г. (b), а также семи произвольных областей в районе исследования (c–i)

Fig. 2. LURR parameter graphs for the epicentral regions of the earthquakes of 15.04.2016 (a) and 08.08.2024 (b), as well as 7 arbitrary regions in the study area (c–i)

ющие «целевых» землетрясений. В столбце 3 таблицы 2 указаны даты землетрясений, которые попали в расчётную область, а также те, которые вне её пределов на расстоянии не более 10 % от линейного размера геометрической фигуры (круга), то есть 0.2°. Зелёным цветом отмечены те землетрясения, которые имели свои аномалии в данной расчётной области, а красным — которые не имели, это так называемые «пропущенные цели». Относительно красных записей в обоих столбцах есть интересные нюансы, которые обсудим позже, а пока, следуя чёткому формализму, будем их считать ложными аномалиями (фиолетовые также попадают в список ложных) и «пропущенными целями».

На основании данных таблицы 2 можно предложить самую простую статистическую оценку. Из 27 аномалий 11 (41 %) попадают в периоды тревоги (два года) для 8 землетрясений. Соответственно, 59 % (16 ано- малий) — ложные. Такая статистика помещает результат работы на уровень среднестатистический и ничем особо не выделяющийся. При этом из 16 ложных аномалий немалое количество (половина) по времени подходит к землетрясениям из списка (фиолетовый цвет), но в пространстве расчётные области, в которых они фиксируются, не состыковываются со своими возможными целями.

Надо отметить, что из восьми таких аномалий в пяти случаях эпицентры выпадают из зоны расчёта не более чем на 0.3–0.8° (за исключением области 1 и 2, где значительно отклоняются три фиолетовые аномалии) и при прогнозе по данным аномалиям в реальных условиях могли бы быть признаны предвестниками. При этом возможная цель для двух из трёх фиолетовых аномалий в первой и второй расчётных областях находится вне области нашего исследования, но выглядит очень интересно. Аномалии в феврале 2011 года

Таблица 1. Независимые землетрясения в районе исследования с М > 6.0

Table 1. Independent earthquakes in the study area with M > 6.0

Год

Year

Месяц |Month

Число Day

Час

Hour

Минута Minute

Широта, °с. ш.

Latitude, °N

Долгота, °в. д. Longitude, °E

Глубина, км Depth, km

Магнитуда, M Magnitude, M

2014

3

13

17

6

33.69

131.89

78

6.2

2016

4

15

16

25

32.75

130.76

12

7.3

2019

5

9

23

48

31.80

131.97

25

6.3

2022

1

21

16

8

32.72

132.07

45

6.6

2024

4

17

14

14

33.20

132.41

39

6.6

2024

8

8

7

42

31.74

131.72

31

7.1

2025

1

13

12

19

31.83

131.57

36

6.6

2025

4

2

14

3

31.05

131.52

36

6.1

Таблица 2. Аномалии параметра LURR в девяти расчётных областях

Table 2. LURR parameter anomalies in nine calculation areas

№ Центр расчётной области Center of the calculation area Землетрясения из табл. 1, попадающие в область прогноза (внутри круга ±10 % от диаметра) Earthquakes from Table 1 that fall within the forecast area (inside the circle ±10 % of the diameter) Даты выхода параметра в область аномальных значений (аномалии) Dates of the parameter entering the area of anomalous values (anomalies) 1 31° с. ш. _131° в. д. 31° N _131° E 02.04.25 (6.1), 13.1.25 (6.6), 9.5.19 (6.3), 8.8.24 (7.1) 20.02.2011, 21.03.2013, 31.02.2017, 10.11.2018, 30.06.2024 2 31° с. ш._132° в. д. 31° N_132° E 02.04.25 (6.1), 13.1.25 (6.6), 9.5.19 (6.3), 8.8.24 (7.1) 11.02.2011, 20.02.2017, 11.05.2021, 16.06.2024 3 31° с. ш._133° в. д. 31° N_133° E n/a n/a 4 32° с. ш._131° в. д. 32° N_131° E 15.4.16 (7.3), 02.04.25 (6.1), 13.1.25 (6.6), 9.5.19 (6.3), 8.8.24 (7.1) 25.10.2013, 07.03.2016, 25.06.2019 5 32° с. ш._132° в. д. 32° N_132° E 9.5.19 (6.3), 8.8.24 (7.1) 13.1.25 (6.6), 02.04.25 (6.1), 21.1.22 (6.6), 17.4.24 (6.6) 27.10.2013, 17.07.2016, 1.07.2019, 19.12.2021, 10.11.2023 6 32° с. ш._133° в. д. 32° N_133° E n/a n/a 7 33° с. ш._131° в. д. 33° N_131° E 15.4.16 (7.3), 13.3.14 (6.2), 21.1.22 (6.6) 30.02.2014, 11.01.2016, 03.03.2018, 22.05.2019, 01.03.2025 8 33° с. ш._132° в. д. 33° N_132° E 13.3.14 (6.2), 21.1.22 (6.6), 17.4.24 (6.6) 18.08.2023 9 33° с. ш._133° в. д. 33° N_133° E 17.4.24 (6.6), 21.1.22 (6.6) 19.08.2013, 19.08.2016, 05.04.2021, 22.08.2023 появились менее чем за месяц до мегаземлетрясения Тохоку 11 марта 2011 года (М = 9). И в этом случае очень короткий временной фактор и масштаб события перекрыли бы любые расстояния от расчётной области до эпицентра (1345 км). Впрочем, пока для чистоты результата по-прежнему считаем все фиолетовые аномалии ложными.

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

В случае землетрясения 13 марта 2014 года единственная аномалия появилась 30 февраля, за две недели до землетрясения (область 7). Первая аномалия до самого сильного землетрясения 15 апреля 2016 года (М = 7.3) была зафиксирована 11 января 2016 года (область 7). Здесь период тревоги составил 3 месяца. В случае землетрясения 9 мая 2019 года единственная аномалия (область 1) соответствует периоду тревоги 6 месяцев. Землетрясение 21 января 2022 года имеет два предвестника (области 5 и 9), первый из которых появился 5 апреля 2021 года, то есть за 9 месяцев. Для землетрясения 17 апреля 2024 года самая первая аномалия появилась 18 августа 2023 года и период тревоги в данном случае составляет 8 месяцев. Другое землетрясение в 2024 году произошло 8 августа, но находилось от апрельского на значительном удалении, а из двух аномалий в областях расчёта 5 и 8 преимуществом пользуется первая, которая находится в одной области с эпицентром (без допустимых отклонений). В этом случае после появления аномалии 10 ноября 2023 года до землетрясения прошло 9 месяцев. Перед землетрясением 13 января 2025 года появляется аномалия 16 июня 2024 года, которая была выделена в области 2 и произошла за 5 месяцев до события. И последнее землетрясение, 2 апреля 2025 года, имеет в областях 1 и 2 схожие аномалии в июне 2024 года, которые хорошо территориально подходят к данному событию. Здесь период тревоги составляет 10 месяцев. То есть для всех землетрясений период тревоги составил менее одного года, что вполовину меньше нами установленного предела. Средний же период по 8 землетрясениям и вовсе составляет около 6.3 месяца, (максимальный — 10, а минимальный — 0.5).

На рис. 3 показаны карты и представлены землетрясения, для которых в конкретной расчётной области выделены аномалии, включая те, на которые указали аномалии, отмеченные фиолетовым цветом. На представленных картах для каждой расчётной области показаны землетрясения (красные маркеры), на которые в расчетных областях фиксировались аномалии (предвестники), и преимущественно эти землетрясения в зоне или поблизости ретроспективно имели прогноз. Наиболее проблемная в данной связи область 4, где целая группа землетрясений не имела предвестников (желтые маркеры на картах). Однако, если посмотреть на карте сейсмичности района (рис. 1), в 4-ю расчётную область оказались вовлечены практически полностью две совершенно разные очаговые области (с эпицентрами сильнейших землетрясений 2016 и 2024 гг.). Результат такого расчёта вряд ли может быть очевидным, и то, что в этой области всё-таки выявлена аномалия на землетрясение 2016 года, на- верное, лучшее, что можно было ожидать. Также интересно выглядит землетрясение 2019 года, на которое в ряде областей, которые включали в себя его эпицентр, аномалий обнаружено не было.

А теперь обсудим отложенный вопрос по ложным тревогам и «пропущенным целям» (красные записи в таблице 2). Из 8 ложных аномалий 6 произошли практически сразу же после землетрясений 2016 года (через 3 месяца), 2019 года (без паузы), 2024 года (через два месяца) и 2025 года (меньше 2 месяцев), что фактически исключает их из разряда ложных и переводит в разряд наведенных. Такие периоды после сильных землетрясений в ближней зоне вполне приемлемы для изменения напряженно-деформированного состояния на соседних разломах, что, конечно же, повлияет на LURR, который напрямую от этого зависит. То есть в реальной ситуации аномалии, которая появилась бы через 1–3 месяца после землетрясения с М > 6.0, внимания уделено не было бы.

По «пропущенным целям» опять же формальное выполнение требований для оценки прогноза многое оставляет за кадром. Например, в области 8 пропущены землетрясения внутри неё в 2014 и 2022 годах. Однако в соседних областях 7 и 9, для которых их эпицентры пограничные и входят в оценку, предвестники для этих землетрясений есть (в 7 — для 2014 года, в 9 — для 2022 года). Получается, что эти землетрясения в нашем исследовании «пропущенными» назвать довольно проблематично с точки зрения здравого смысла. Такая же история с землетрясением 2019 и других годов, где отрабатывали разные (но не все, куда можно отнести эпицентры) области расчёта.

И если с «пропущенными целями» всё относительно радужно, то что по итогу с ложными аномалиями? Если признать сам факт наведенности и наличия триггерного эффекта в природе как такового, то 6 аномалий надо исключить из рассмотрения. Таким образом, останется всего две ложные аномалии (без учёта фиолетовых) из 27, что фактически приближает этот параметр к ничтожному. При этом, согласно картам на рис. 3, из 7 «фиолетовых» аномалий — прогнозов на землетрясения вне области расчёта — только 3 землетрясения, включая Тохоку 2011 года, находятся на большом расстоянии от своих областей расчёта. Остальные 4 вполне можно считать также успешными прогнозами.

Если сравнить результаты расчётов с аномалиями в эпицентральных зонах, то в целом они имеют соответствие, а аномалия в июне 2015 года (рис. 2, a) произошла раньше, чем единственная аномалия в январе 2016 года, что для эпицентральной области вполне приемлемо. Эпицентры землетрясений 2022 и 2024 годов находятся недалеко друг от друга, поэтому аномалию на рис. 2, б в декабре 2021 года можно отнести к землетрясению в 2022 году.

Таким образом, если принять соответствующие допущения, то есть признать реализованными 4 аномалии перед несильно отдаленными землетрясениями (из числа фиолетовых), а также исключить 6 аномалий из числа ложных, определив их как наведенные, то можно констатировать, что статистика сильно меняется. Из 21 аномалии 15 (72 %) попадают в периоды тревоги для 8 землетрясений, а ложными являются всего 6 (28 %). Пропущенных событий нет вовсе, так 43

Рис. 3. Карта распределения ретроспективных прогнозов землетрясений в девяти расчётных областях: красные маркеры — эпицентры землетрясений, для которых в расчётной области есть аномалия (успешный прогноз), желтые маркеры — эпицентры землетрясений в области расчёта без предвестников («пропущенные цели»), синие маркеры — оставшиеся (вне прогнозных категорий), по данным таблицы 2

Fig. 3. The map of the distribution of retrospective earthquake forecasts for nine calculation areas: red markers — earthquake epicenters for which there is an anomaly in the calculation area (successful forecast), yellow markers — earthquake epicenters in the calculation area without precursors («missed targets»), blue markers — remaining according to the data in Table 2 (outside the forecast categories)

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

Заключение

В работе представлены результаты анализа сейсмичности в районе о. Кюсю (Япония, стык Филиппинской и Евразийской плит) методом LURR c 2010 по 2025 год. с целью ретроспективной оценки прогноза землетрясений с М > 6.0. С учётом предложенных допущений (обоснованных на качественном уровне) результат сле- 44

дующий: 72 % успешных прогнозов, 28 % ложных тревог и отсутствие «пропущенных целей». Отмечено, что результат достигается в том числе достаточным количеством расчётных областей. Важным итогом также стало полученное значение для среднего периода тревоги (около полугода), которое в четыре раза меньше, чем традиционно применяемое авторами для оценок данным методом.

Работа выполнена при поддержке Минобрнауки России в рамках государственного задания № 075-00609-26 Института морской геологии и геофизики ДВО РАН и с использованием данных, полученных на уникальной научной установке «Сейсмоинфразвуковой комплекс мониторинга арктической криолитозоны и комплекс непрерывного сейсмического мониторинга Российской Федерации, сопредельных территорий и мира» , .