Seismicity of Morocco and simulation of 08.09.2023 earthuake consequnces with “Extremum” system application

Full Text

ВВЕДЕНИЕ

Альпийско-Гималайский сейсмоактивный пояс (АГСП) (Средиземноморско-Трансазиатский по [5]) чрезвычайно сейсмически активен в 2023 г. Начало было положено разрушительными Турецкими землетрясениями 06.02.2023 г. с Мw = 7.8 и 7.5 на Восточно-Анатолийском разломе, где несколько веков не было столько мощного проявления сейсмичности с катастрофическими последствиями [3, 17]. В противоположных от Турецких очагов направлениях в АГСП произошли землетрясения с разрушительными эффектами: через 215 сут – на западе, в северо-западной части Африки 8 сентября 2023 г. с М = 6.8 (NEIC) в Марокко, через 246 сут – на востоке, в Афганистане, в западной части страны серия землетрясений близкой магнитуды 7 октября 2023 г. (два толчка с разницей в 30 мин. по времени возникновения), 11 и 15 октября 2023 г. с М = 6.3–6.4 и глубинами в верхней части земной коры (оперативные данные ЕГС ФИЦ РАН) (рис. 1).

 

Рис. 1. Карта сейсмичности в Альпийско-Гималайском поясе, по данным ССД ФИЦ ЕГС РАН (выборка сделана от эпицентра Турецкого землетрясения с радиусом 5500 км).

 

Можно высказать предположение о триггерном влиянии Турецких землетрясений 6 февраля на сейсмичность значительной части АГСП. Так называемым “деформационным” волнам, в том числе возникшим после сильного землетрясения, как одному из “триггерных” механизмов подвижности разломно-блоковой структуры континентальной литосферы и генерации сейсмической активности [5] уделяется большое внимание в работах плеяды советских и российских ученых. Довольно полное описание этих работ приводится в монографии С.И. Шермана, посвященной тектонофизической концепции сейсмического процесса и прогнозу землетрясений.

Предметом настоящей статьи является исследование сильнейшего Марокканского землетрясения 08.09.2023 г. с М = 6.8 [13] с целью возможной калибровки ГИС “Экстремум” для этого района.

Авторы этой работы задались целью создания “однородного по магнитуде каталога”. Предполагается провести исследование коэффициентов макросейсмического уравнения для этого района. Для этого создана выборка по различным источникам [7, 10, 12] для Марокко и сопредельных территорий. Анализ результатов моделирования последствий и сравнение с наблюденными данными по интенсивности позволяет сделать вывод о хорошей сходимости между наблюденными и расчетными интенсивностями I_i при использовании уравнений для зоны Верхнего и Среднего Атласа [10], для Алжира и Атласских гор [7] и для территории Турции [12]. Сравнение полученных оценок ущерба для населенных пунктов Adassil (Адассиль) и Targua (Тарга), расположенных в сильно пострадавшей провинции Chichaoua (Чичауа), с космическими снимками показывает удовлетворительную сходимость в оценках расчетных и наблюденных средних степеней повреждений застройки населенных пунктов.

СЕЙСМИЧНОСТЬ АФРИКАНСКОГО КОНТИНЕНТА

Долгое время огромный Африканский континент считался асейсмичным, вплоть до начала ХХ в. были известны лишь незначительные проявления сейсмичности в районе гор Атласа, Абиссинского нагорья, в районе о. Мадагаскар и локально в некоторых других территориях на юге и востоке континента. Известно, что исследованию сейсмичности этого континента уделяли внимание, кроме зарубежных, и советские ученые. Среди них особое место занимает работа “О сейсмичности Африки” выдающего ученого геолога-тектониста-сейсмолога Г.П. Горшкова [1, 14, 15]. Как пишет сам автор, эта работа была выполнена по инициативе и заданию Департамента естественных наук ЮНЕСКО на кафедре динамической геологии МГУ им. М.В. Ломоносова в содружестве с Институтом физики Земли АН СССР и Советом по сейсмологии АН СССР. Г.П. Горшков посчитал своим долгом широко обсудить созданный труд с советскими и зарубежными коллегами, среди которых были, например, Ч. Рихтер и Б. Гутенберг, которые в 1949 и 1954 г. в монографии о сейсмичности Земли [16] опубликовали каталог землетрясений мира, в том числе и Африки [1]. Данные этого каталога послужили Г.П. Горшкову основой для написания монографии “О сейсмичности Африки”. Были составлены списки землетрясений различных сейсмотектонических зон Африки, построен ряд карт, характеризующих распределение сейсмичности и сейсмотектонические зоны (рис. 2 и 3).

 

Рис. 2. Копия рисунка “Тектоническая карта Африки” из [1] с обозначением эпицентра землетрясения 08.09.2023 г.

1 – докембрийский фундамент Африканской платформы, выходящий на поверхность; 2 – области неглубокого погружения докембрийского фундамента; 3 – области глубокого погружения докембрийского фундамента; 4 – внутриплатформенные каледонские складки; 5 – герцинские складки; 6 – капские складки; 7 – предгорный прогиб Карру; 8 – альпийские складки; 9 – Месопотамско-Сирийский предгорный прогиб; 10 – основные разломы.

 

Рис. 3. Сейсмотектонические зоны Африки (частичная копия рисунка карты “Эпицентров землетрясений Африки” из [1].

1 – Восточно-африканские грабены; 2 – границы сейсмотектонических зон I–V: I – область альпийской складчатости (система горных сооружений Атласа), II – область разломов “Великого Рифта Африки”, III – область капских складок, IV – область неглубокого залегания докембрийского кристаллического фундамента, V – то же, что и IV, но по восточному побережью Африки.

 

Установлено, что большую часть территории Африканского континента занимает Африканская платформа. Проявление сейсмичности связано в основном с областью альпийской складчатости (система горных сооружений Атласа) и областью разломов “Великого Рифта Африки” [1, 14, 15].

В более современных исследованиях, например в [8, 11], также подтверждается, что землетрясения в Марокко не превышали М = 6.0 и в основном происходили на Средиземноморском побережье и в районе гор Атлас (рис. 4). Известно разрушительное землетрясение 29.02.1960 г. в Агадире (координаты 30.45° с.ш., 9.62° з.д., Md = 5.9), эпицентр которого был удален примерно на 200 км к юго-западу от эпицентра 08.09.2023 г. По [11] “ужасное землетрясение опустошило Агадир и его окрестности, ущерб катастрофичен: более 12 000 человек погибли, тысячи получили ранения, было разрушено более 75% зданий. В населенных пунктах Ячеч, Касба, Адуар, Фаунти Талборжт более 90% зданий были разрушены или повреждены. Несмотря на небольшое значение магнитуды (Md = 5.9), масштаб негативных последствий можно объяснить низким качеством строительства и небольшой глубиной гипоцентра (3 км), а также близостью к городу”.

 

Рис. 4. Сейсмичность Марокко и его окружения (1901–2010 гг.) [11].

 

В [18] опубликован каталог основных марокканских землетрясений с 1045 по 2005 г. Авторы этой работы задались целью создания каталога однородных по магнитуде землетрясений для дальнейшего его использования в исследовании сейсмической опасности. В представленном каталоге лишь одно 9-балльное землетрясение 1731 г. имело Mw = 6.4, его эпицентр находился примерно в 200 км к юго-западу от эпицентра землетрясения 08.09.2023 на побережье Атлантического океана (видимо в близкой зоне к эпицентру землетрясения 29.02.1960 г.). Эпицентры всех других сильных землетрясений Марокко были приурочены к рифтовым структурам в Средиземноморье.

Из вышеизложенного следует, что землетрясение 08.09.2023 г. было самым сильным в горах Атласа за всю историю, и, безусловно, может быть предметом для нашего исследования в целях калибровки системы “Экстремум” для районов северо-западной части Африки.

МОДЕЛИРОВАНИЕ ПОСЛЕДСТВИЙ ЗЕМЛЕТРЯСЕНИЯ 08.09.2023 г. С М = 6.8 В МАРОККО

Расчеты моделирования последствий землетрясения проводятся с помощью системы “Экстремум” в первые минуты после определения параметров землетрясения для принятия решения о проведении спасательных и других экстренных операций. Точность оперативного прогноза последствий существенно повышается, если для данной сейсмоопасной территории заблаговременно определены рейтинг сейсмологических служб и параметры зависимостей, описывающих особенности затухания сейсмической интенсивности I_i, и региональные особенности поведения элементов риска.

В качестве региональных зависимостей затухания сейсмической интенсивности в данной статье были использованы результаты исследований, полученные в работах [7, 10] для Марокко и сопредельных территорий. Региональные оценки затухания интенсивности для зоны Атласских гор и Алжира [7] оказались близкими к зависимостям N.N. Ambraseys [6] для Балканского региона. Этот факт объясняется принадлежностью рассматриваемых территорий к одной и той же тектонической зоне − зоне альпийской складчатости [1, 14], и позволяет тестировать применимость аналогичных зависимостей для моделирования последствий землетрясения 08.09.2023 г. в Марокко.

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

• для зоны Верхнего и Среднего Атласа (HM) [10], h_ср = 9.7 км:

$I\text{\hspace{0.33em}=\hspace{0.33em}1.72}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+1.\text{51,}$ (1)

$I\text{\hspace{0.33em}=\hspace{0.33em}1.4}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+2.5;$ (2)

• для зоны Западного Высокого Атласа (HO) [10], h_ср = 1.4 км:

$I\text{\hspace{0.33em}=\hspace{0.33em}1.72}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-1.\text{02,}$ (3)

$I\text{\hspace{0.33em}=\hspace{0.33em}1.4}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+0.62;$ (4)

• для зоны Верхнего и Среднего Атласа (HM) [10], h_usgs = 26 км:

$I\text{\hspace{0.33em}=\hspace{0.33em}1.72}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+2.\text{81,}$ (5)

$I\text{\hspace{0.33em}=\hspace{0.33em}1.4}M-1.\text{3ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0013}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+3.5;$ (6)

• для Алжира [7]:

h_ср = 6.82 км,

$I\text{\hspace{0.33em}=\hspace{0.33em}1.43}M-2.\text{28ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0004}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+6.29\pm 0.\text{24}P\text{,}$ (7)

h_usgs = 26 км,

$I\text{\hspace{0.33em}=\hspace{0.33em}1.43}M-2.\text{28ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{0004}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+6.29\pm 0.\text{24}P;$ (7а)

• для Атласских гор [7]:

h_ср = 4.82 км,

$I\text{\hspace{0.33em}=\hspace{0.33em}1.48}M-2.\text{05ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{00074}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+5.\text{16}\pm 0.\text{35}P\text{,}$ (8)

h_usgs = 26 км,

$I\text{\hspace{0.33em}=\hspace{0.33em}1.48}M-2.\text{05ln}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)-0.\text{00074}\left(\sqrt{{\Delta }^2\text{+}{h}^2}\right)+5.\text{16}\pm 0.\text{35}P;$ (8а)

• для юго-восточной Турции [12]:

$I \text{\hspace{0.33em}=} 1.\text{82}M-1.\text{38ln}\left(\Delta \text{\hspace{0.33em}+\hspace{0.33em}}7\right)-0.\text{99;}$ (9)

• для территории Турции [12]:

$I \text{=} 2.\text{12}M-1.\text{38ln}\left(\Delta \text{\hspace{0.33em}+\hspace{0.33em}}7\right)-2.\text{72,}$ (10)

где Δ − эпицентральное расстояние, км; h − глубина очага, км; М − магнитуда землетрясения (инструментально определенное значение Ms или пересчетное из других оценок).

Моделирование возможного распределения интенсивностей для события 08.09.2023 г. было выполнено с использованием уравнений (1) – (10) для коэффициента сжатия эллипса k = 1.5 с учетом ориентации эллипса макросейсмического поля под углом 255° в соответствии с решением механизма очага.

Учитывая небольшой разброс в определении координат эпицентра сейсмологическими службами в оперативном режиме (минимальное расстояние 1.1 км, среднее – 14.1 км) (рис. 5) и то, что плотность населения провинции Марракеш-Сафи по данным http://www.statoids.com/uma.html на 2014 г. составляла 117.6 чел/км², что в 2.6 раза выше, чем в среднем по стране, расчеты проводились для параметров USGS φ = 31.064° с.ш.; λ = 8.391° з.д. Использовалась глубина очага h_usgs = 26 км по данным USGS и средние глубины для зоны h_ср в соответствии с рекомендациями авторов зависимостей [7, 10].

 

Рис. 5. Координаты эпицентра землетрясения 08.09.2023 г. в Марокко по данным разных сейсмологических служб в оперативном режиме.

 

Уравнения (1)–(6) получены преобразованием уравнения Шпонхоера-Ковеслигети [19]:

$I_o-I=3log\left(\frac{\sqrt{R^2+h^2}}{h}\right)+1.3a\left(\sqrt{R^2+h^2}-h\right)$ (11)

к виду

$I=aM-bln\left(\sqrt{R^2+h^2}\right)+c\left(\sqrt{R^2+h^2}\right)+d,$ (12)

где Δ − эпицентральное расстояние, км; h − глубина очага, км; М − магнитуда землетрясения, I₀ – интенсивность в эпицентре, балл; a, b, c, d – коэффициенты.

Для пересчета магнитуд в (1), (3), (5) использовалась зависимость M(I₀) вида:

$M=0.58I_0+\mathrm{0.85,}$ (13)

в (2), (4), (6):

$M=0.716I_0+0.517\text{\hspace{0.33em}log}\left(h\right)-0.207$, (14)

где h − глубина очага, км; М − магнитуда землетрясения, I₀ – интенсивность в эпицентре, балл.

Для выбора уравнения затуханий интенсивности, которые в случае повторения события могут быть использованы в оперативной оценке потерь, расчетные значения интенсивности I_i сравнивались с наблюденными по данным GDACS [https://www.gdacs.org/].

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

На рис. 6 приведены результаты сравнения расчетных и наблюденных значений интенсивностей сотрясений I (точки и аппроксимирующая их кривая) при использовании уравнений затухания (1)–(6) [10].

 

Рис. 6. Сравнение наблюденной и расчетных интенсивностей для уравнений (1)–(6) [10].

 

Анализ рис. 6 показывает систематическое завышение расчетной интенсивности I при использовании уравнений (1), (5) и (6). При использовании уравнения (5) разница составляет несколько баллов. Уравнения (3) и (4) дают систематическое занижение расчетной интенсивности I. На расстояниях до 20–30 км эта разница незначительна, порядка 0.2–0.3 балла, с увеличением эпицентрального расстояния разница возрастает и на расстояниях 80–120 км достигает 1.5–2 балла. Наилучшая сходимость получена для уравнения (2) для зоны Верхнего и Среднего Атласа. На расстояниях до 50 км разница составляет 0.2 балла, на расстояниях свыше 160 км она увеличивается на 0.4 балла.

На рис. 7 приведены результаты сравнения расчетных и наблюденных значений интенсивностей сотрясений I (точки и аппроксимирующая их кривая) при использовании уравнений затухания (7, 7а, 8, 8а), а также (9) и (10).

 

Рис. 7. Сравнение наблюденной и расчетных интенсивностей для уравнений (7), (7а), (8), (8а) по [7]; (9), (10) по [12].

 

Анализ рис. 7 показывает наилучшую сходимость расчетных и наблюденных значений интенсивностей I для уравнения (10) для территории Турции]. Уравнения (7), (7а) для Алжира и (8), (8а) для Атласских гор дают завышенные оценки расчетной интенсивности I_i от 0.3 до 1 балла на близких расстояниях, на расстояниях более 120 км расчетные значения I_i приближаются к наблюденным. Но уравнение (7а) дает почти совпадающие оценки расчетных и наблюденных I_i на эпицентральных расстояниях ∆ > 70 км.

В табл. 1 приведены значения отклонений наблюденных и расчетных интенсивностей ∆I при использовании разных уравнений затухания.

 

Таблица 1. Отклонения наблюденных и расчетных интенсивностей ∆I для разных баллов при использовании уравнений (1)–(10)

I, балл	Уравнения затухания интенсивности
	(1)	(2)	(3)	(4)	(5)	(6)	(7)	(7a)	(8)	(8a)	(9)	(10)
5	2.0	0.8	–0.9	–1.6	3.2	1.8	0.3	0.3	0.4	0.4	0.0	0.2
6	1.4	0.3	–1.2	–2.0	2.7	1.2	0.1	–0.1	0.3	0.2	–0.6	–0.3
7	1.2	0.1	–1.5	–1.9	2.3	0.9	0.3	0.1	0.4	0.1	–0.9	–0.7
8	1.2	0.0	–1.5	–1.7	1.8	0.3	1.0	0.0	1.0	–0.2	–0.8	–0.7
Iср	1.6	0.5	–1.1	–1.7	2.8	1.4	0.3	0.1	0.4	0.3	–0.4	–0.1

 

Анализ рис. 6, 7 и табл. 1 позволяет сделать вывод о хорошей сходимости между наблюденными и расчетными интенсивностями I_i при использовании уравнений (2), (7а), (8а) и (10).

На рис. 8 приведены результаты моделирования последствий землетрясения 08.09.2023 г. с М = 6.8 в Марокко с помощью системы “Экстремум” для уравнений затухания интенсивности (2) для зоны Верхнего и Среднего Атласа и (8а) для Атласских гор. Степени повреждений зданий приведены по шкале MMSK-86 [4], где d = 1 – легкие повреждения, d = 2 – умеренные повреждения, d = 3 – тяжелые повреждения, d = 4 – частичные разрушения, d = 5 – обвалы. Для расчета степеней повреждения d использовалась матрица повреждений по [9].

 

Рис. 8. Моделирование последствий землетрясения 08.09.2023 г. с М = 6.8 в Марокко с помощью системы “Экстремум” для уравнения: а − (8a), б − (2). Средние степени повреждений: d0 – нет повреждений, d1 – легкие, d2 – умеренные, d3 – тяжелые, d4 – частичные разрушения.

 

При использовании уравнений (7а) для Алжира и (10) для территории Турции были получены близкие результаты в определении средних степеней повреждений d_ср зданий в населенных пунктах при использовании матрицы повреждений [9].

Для изучения влияния матрицы повреждений на оценки потерь были выполнены расчеты последствий для уравнения (7а) с использованием матрицы [9] и матрицы по шкале EMS-92.

Анализ рис. 9 показывает, что оценки вероятности отдельных расчетных степеней повреждения зданий при использовании разных матриц повреждений имеют значительный разброс. Для населенного пункта Эминдуни оценки вероятности степени d = 1 для разных матриц отличаются на 27%. Для населенных пунктов Марракеш и Эминдуни оценки вероятности возможных степеней d = 2 не превышает 14%. А для Адассиль разница в оценках вероятности степени повреждения d = 3 при использовании разных матриц составила 60%.

 

Рис. 9. Диаграммы средних степеней повреждений для отдельных населенных пунктов для (7а) для Алжира [7]: а − матрица повреждений по [9]; б − матрицы повреждений по шкале EMS-92.

 

В табл. 2 приведены значения вероятностей разных степеней повреждений для отдельных населенных пунктов при использовании уравнения (7а) и двух матриц повреждений.

 

Таблица 2. Вероятности степеней повреждений для отдельных населенных пунктов при использовании уравнения (7а) для Алжира

Населенный пункт	Вероятность степеней повреждений при использовании матриц по [٩] (в числителе) и по шкале EMS 92 (в знаменателе)
	d = 1	d = 2	d = 3	d = 4	d = 5	dср
Adassil (Адассиль)	0.05/0.12	0.20/0.13	0.35/0.14	0.30/0.45	0.10/0.11	3.20/3.15
Targua (Тарга)	0.26/0.14	0.36/0.14	0.22/0.45	0.09/0.11	0.01/0.0	2.05/2.21
Imindounite (Эминдуни)	0.38/0.28	0.26/0.30	0.11/0.14	0.03/0.02	0.0/0.0	1.37/1.39
Marrakech (Марракеш)	0.38/0.51	0.16/0.14	0.05/0.0	0.0/0.0	0.0/0.0	0.87/0.80

 

Разница в оценке средних степеней повреждений d_ср варьирует от 1.4% для Эминдуни до 8% для Марракеша (см. табл. 2). Т.е. использование разных матриц повреждения не оказывает сильного влияния на оценки средних степеней повреждения d_ср для населенных пунктов в целом с учетом застройки разного типа.

Сравнение полученных оценок ущерба для Адассиль и Тарга, расположенных в сильно пострадавшей провинции Chichaoua (Чичауа), с космическими снимками показывает удовлетворительную сходимость в оценках расчетных и наблюденных потерь. На рис. 10 приведены космические снимки UNOSAT из предварительного отчета EQ20230909MAR от 10.09.2023 г., на которых наблюдается соответствие расчетных и реальных оценок степеней повреждения d = 4 и d = 5 в этих населенных пунктах.

 

Рис. 10. Снимки населенных пунктов: а – Адассиль, б – Тарга.[UNOSAT. Prlimnary_Assessment_Report_EQ20230909MAR_Adassil_10Sep2023-1].

 

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

В процессе калибровки модели макросейсмического поля системы “Экстремум” для территории Марокко выполнен анализ применимости уравнений затухания сейсмической интенсивности, полученных исследователями в разные годы для зоны Верхнего и Среднего Атласа [10], Алжира и Атласских гор [7] и территории Турции [12].

Для оценки возможных последствий землетрясения 08.09.2023 г. применялись методы имитационного компьютерного моделирования с использованием системы “Экстремум”, разработанной с участием авторов данной статьи. Результаты моделирования последствий землетрясения показали хорошую сходимость наблюденных и расчетных значений интенсивности I_i при использовании уравнений для зоны Верхнего и Среднего Атласа, Алжира и Tурции.

При использовании уравнения для зоны Верхнего и Среднего Атласа [10] на расстояниях до 50 км разница ∆I составляет 0.2 балла, на расстояниях более 160 км − 0.4 балла. Наилучшая сходимость наблюденных и расчетных значений интенсивностей I_i получена для уравнения для территории Турции [12], на всех эпицентральных расстояниях ∆I не превышает 0.1 балла. Уравнения для Алжира и Атласских гор [7] дают завышенные оценки расчетной интенсивности I_i от 0.3 до 1 балла на близких расстояниях, на расстояниях более 120 км расчетные значения I_i приближаются к наблюденным. На эпицентральных расстояниях ∆ > 70 км уравнение для Алжира [7] дает почти совпадающие оценки расчетных I_i и наблюденных значений интенсивности.

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

В настоящем исследовании расчетные значения интенсивности I_i сравнивались с наблюденными значениями интенсивности для населенных пунктов в соответствии с информацией, размещенной на сайте GDACS [gdacs.org], которые, судя по расхождениям с информацией, размещенной на сайте DYFI Геологической службы США, также является неокончательной.

Информация на сайте DYFI дает завышенные оценки наблюденной интенсивности по сравнению данными GDACS. Разница достигает 0.7–0.8 балла (рис. 11). На сайте DYFI приведены данные об ощущениях для населенных пунктов на расстояниях 50 км и более от эпицентра.

 

Рис. 11. Сравнение наблюденных значений интенсивности в населенных пунктах по данным GDACS и DYFI USGS.

 

Несмотря на то, что полученные оценки расчетных интенсивностей носят предварительный характер использование уравнения затухания сейсмической интенсивности для Алжира [7] позволяет получить расчетные средние степени повреждения d_ср застройки для наиболее сильно пострадавших населенных пунктов близкие к наблюденным на 10 января 2023 г. [UNOSAT_Preliminary_Assessment_Report_EQ20230909MAR_Adassil_10Sep2023-1].

ВЫВОДЫ

Процедура калибровки региональных моделей затухания сейсмической интенсивности предполагает выделение границ зон с квазистабильными параметрами поля. В дальнейшем для уточнения калибруемых параметров поля и определения границ зон для их применения необходимо проведение работ по сбору информации о макросейсмических проявлениях прошлых событий в рассматриваемом регионе.

Для успешной калибровки моделей макросейсмического поля системы “Экстремум” и других глобальных систем, используемых для оперативной оценки последствий землетрясений, необходимо создание и своевременное пополнение базы знаний о макросейсмических и инженерно-сейсмологических последствиях землетрясений. В этой связи представляется важным международное сотрудничество в области обработки “больших данных” (Big Data) и создание распределенной базы знаний о физических и социально- экономических последствиях прошлых сильных землетрясений в рамках проектов ООН/КОДАТА с указанием границ калибровочных зон.

Авторы выражают благодарность своим коллегам за постоянную поддержку, обсуждение методических вопросов по калибровке моделей для оценки последствий землетрясений в оперативном режиме. Особая благодарность сотрудникам ООО ЦИЭКС за плодотворное сотрудничество и вклад в развитие среды системы “Экстремум”

Статья подготовлена в рамках выполнения государственного задания ИГЭ РАН по теме НИР № 122022400105-9 “Прогноз, моделирование и мониторинг эндогенных и экзогенных геологических процессов для снижения уровня их негативных последствий” и государственного задания ФИЦ ЕГС РАН № 075-01271-23.

About the authors

N. I. Frolova

Sergeev Institute of Environmental Geoscience, Russian Academy of Sciences

Author for correspondence.
Email: frolovanina7@gmail.com
Russian Federation, Ulanskii per., 13, bld. 2, Moscow, 101000

I. P. Gabsatarova

Geophysical Survey, Russian Academy of Sciences

Email: ira@gsras.ru
Russian Federation, Lenina av., 189, Obninsk, Kaluga region, 249035

S. P. Suchshev

Moscow State Technical University named after Bauman

Email: Sersan150@mail.ru
Russian Federation, 2-ya Baumanskaya str., 5, Moscow, 105005

N. S. Malaeva

Sergeev Institute of Environmental Geoscience, Russian Academy of Sciences

Email: frolovanina7@gmail.com
Russian Federation, Ulanskii per., 13, bld. 2, Moscow, 101000

References

Supplementary files