Научный журнал
Успехи современного естествознания
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,823

О ПРИМЕНЕНИИ МЕТОДОВ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СПУТНИКОВЫХ СИСТЕМ ДЛЯ ЦЕЛЕЙ РАННЕГО ПРЕДУПРЕЖДЕНИЯ О ЦУНАМИ В КУРИЛО-КАМЧАТСКОМ РЕГИОНЕ

Нечаев Г.В. 1, 2 Шестаков Н.В. 1, 2 Такахаши Х. 3 Герасименко М.Д. 2 Сысоев Д.В. 1, 2
1 ФГАОУ ВО Дальневосточный федеральный университет
2 ФГБУН Институт прикладной математики Дальневосточного отделения Российской академии наук
3 Институт сейсмологии и вулканологии Хоккайдского университета
Действующая система предупреждения о цунами (СПЦ) в Курило-Камчатском регионе основана на применении сейсмических методов и не позволяет получить в режиме реального времени точные параметры очага землетрясения, необходимые для надежного прогнозирования цунами. Кроме того, сейсмическим методам свойственна недооценка магнитуды события при возникновении сильных землетрясений (Mw > 8,0), нередко происходящих в Курило-Камчатской зоне субдукции. В последние годы для решения вышеописанных проблем рядом зарубежных СПЦ широко применяются методы Глобальных навигационных спутниковых систем (ГНСС-методы). Суть данной методологии заключается в инвертировании детектированных ГНСС-станциями косейсмических смещений в параметры очага землетрясения и последующем расчете начальных смещений морского дна и параметров цунами. В данной работе на примере Камчатского землетрясения 1952 г. Mw8,7–9,0 рассматривается возможность применения ГНСС-методов в региональной СПЦ. Косейсмические смещения на местах расположения региональных ГНСС-станций, смоделированные программой Coulomb3.0 по опорной модели, были инвертированы с использованием различной точности их определения. По инвертированным плоскостям землетрясений далее были рассчитаны косейсмические смещения дна и сравнены со смещениями, соответствующими опорной модели. В результате определено, что инвертирование данных даже простой одноплоскостной моделью очага способно удовлетворительно описать реальный очаг, а соответствующие инвертированным очагам начальные движения морского дна довольно точно повторяют реальные смещения. Точности определения смещений в 3 см по горизонтали и 5 см вертикали достаточно для адекватной оценки параметров землетрясения. Таким образом, подтверждено, что применение ГНСС-методов и использование существующей региональной ГНСС-сети способны существенно улучшить надежность СПЦ при нынешних возможностях ГНСС-технологий.
ГНСС-методы
предупреждение о цунами
моделирование землетрясений
инверсия
косейсмические смещения
1. Пинегина Т.К., Базанова Л.И. Новые данные о параметрах исторических цунами на побережье Авачинского залива (Камчатка) // Вестник КРАУНЦ. Науки о Земле. 2016. № 3 (31). С. 5–17.
2. Kawamoto S., Ohta Y., Hiyama Y., Todoriki M., Nishimura T., Furuya T., Sato Y., Yahagi T., Miyagawa K. REGARD: A new GNSS-based real-time finite fault modeling system for GEONET. Journal of Geophysical Research: Solid Earth. 2017. vol. 122. P. 1–26. DOI:10.1002/2016JB013485.
3. MacInnes B., Weiss R., Bourgeois J., and Pinegina T.K. Slip distribution of the 1952 Kamchatka Great earthquake based on near-field tsunami deposits and historical records. Bulletin of the Seismological Society of America. 2010. vol. 100. no. 4. P. 1695–1709. DOI: 10.1785/0120090376.
4. Matsuura M., Hasegawa Y. A maximum likelihood approach to nonlinear inversion under constraints. Physics of the Earth and planetary interiors. 1987. vol. 47. P. 179–187.
5. Fratarcangeli F., Savastano G., Chiara D’Achille M., Mazzoni A., Crespi M., Riguzzi F., Devoti R., Pietrantonio G. VADASE reliability and accuracy of real-time displacement estimation: application to the central Italy 2016 earthquakes. Remote sensing. 2018. 10(8). 1201. P. 1–20. DOI: 10.3390/rs10081201.

Полуостров Камчатка и Курильские острова расположены в районе высокой сейсмической активности. Регулярные землетрясения происходят здесь в результате сложного взаимодействия Тихоокеанской и Северо-Американской плит и Берингийской и Охотской микроплит. Наиболее сильные землетрясения с магнитудой до Mw9,0 генерируются в Курило-Камчатской зоне субдукции, в зоне контакта Тихоокеанской и Охотской плит и являются наиболее частым источником цунами. Также регион подвержен воздействию телецунами, возникающих в результате катастрофических землетрясений, происходящих у берегов Северной и Южной Америки. Так, в 1960 г. высота заплесков Чилийского цунами на побережье Камчатки достигала 6 м [1]. Однако наиболее сильные и разрушительные для Курило-Камчатского побережья цунами генерируются именно в близлежащей одноименной зоне субдукции. Наиболее разрушительными в XX в. стали Камчатское землетрясение Mw8,7–9,0 и цунами 1952 г., когда г. Северо-Курильск был разрушен и погибли более 2000 человек. В связи с этим случаем к 1958 г. была разработана региональная система предупреждения о цунами (СПЦ-1958). После трагических событий 2004 г., связанных с Суматра-Андаманским землетрясением и цунами, система была усовершенствована. Однако и сейчас она основана на использовании сейсмических методов, и ее главной задачей является определение магнитуды и местоположения сейсмического события. Подобный подход имеет значительные недостатки, главным из которых является недооценка магнитуды сразу после сильного землетрясения (Mw > 8,0) в связи с насыщением сигнала. Другим несовершенством таких систем является невозможность определения геометрических параметров сейсморазрыва, а также величины смещения в очаге.

В последние годы развивается принципиально новая методика определения параметров очага землетрясения в реальном времени, основанная на применении технологий Глобальных навигационных спутниковых систем (ГНСС). Суть метода заключается в использовании вызванных землетрясением косейсмических смещений, зафиксированных ГНСС-станциями, для определения магнитуды события и параметров очага, которые далее применяются для моделирования смещений морского дна и расчета характеристик цунами. Данный способ свободен от проблемы насыщения сигнала и потому хорошо работает в случае сильнейших землетрясений, которые происходят в зонах субдукции. Особенно активно системы, использующие ГНСС-методы, начали разрабатываться после 2004 и 2011 гг., когда действовавшие на тот момент СПЦ существенно недооценили магнитуду землетрясений, что привело к огромным человеческим жертвам. На сегодняшний день существует несколько развивающихся и работающих систем раннего предупреждения о землетрясениях и цунами, основанных на ГНСС-методах: GITEWS (Индонезия), G-larmS и G-FAST (США), REGARD (Япония), из которых последняя является наиболее развитой и находится в опытной эксплуатации с 2017 г. [2]. Ретроспективная оценка магнитуды землетрясения Тохоку 2011 г. системой REGARD составила Mw8,7, что намного ближе к реальной магнитуде Mw9,0, чем оценка, полученная Японским Метеорологическим Агентством (Mw7,9) сразу после события. Подобный опыт позволяет предположить перспективность использования ГНСС-методов для предупреждения о цунами в Курило-Камчатском регионе, что даст возможность сократить число ложных тревог и повысить надежность и скорость оценки опасности цунами.

В данной работе на примере Камчатского землетрясения 1952 г. Mw8,7–9,0 рассмотрены возможность применения ГНСС-методов для целей предупреждения о цунами, а также влияние точности определения косейсмических смещений на оценку параметров землетрясения и движений морского дна.

Методология

Искомые параметры очага землетрясения получают посредством инверсии наблюденных косейсмических смещений. Таким образом, смещения на станциях наблюдений являются наиболее важной информацией для рассматриваемого типа систем. К сожалению, для события 1952 г. подобных данных не существует. Поэтому косейсмические смещения были смоделированы на местах положения действующих (либо прекративших работу по различным причинам) станций постоянных ГНСС-наблюдений региональных сетей (ДВО РАН, KamNET). Расчеты производились с помощью программного пакета Coulomb3.0, в котором модель Земли представлена однородным полупространством с жесткостью μ = 3,5×1010 [3]. В качестве опорной (эталонной) была использована многоплоскостная модель очага, определенная MacInnes [3] (рис. 1).

nechaev1.tif

Рис. 1. Региональные ГНСС-станции (красно-белые квадраты) и соответствующие им косейсмические смещения (синие стрелки) рассчитанные по модели очага землетрясения MacInnes [3] (голубые квадраты). Цифрами обозначены смещения в каждой подплоскости (в метрах)

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

M0 = μLWD, (1)

nechaev01.wmf (2)

Для инверсии рассчитанных смещений применялся метод, предложенный Matsuura & Hasegawa [4], в котором так же, как и в Coulomb3.0, модель Земли представлена однородным полупространством. Данный способ является одним из наиболее распространенных; кроме того, он позволяет довольно быстро определить параметры землетрясения, что важно для целей раннего предупреждения о цунами. В качестве начальных параметров для инверсии была задана сейсмофокальная плоскость, соответствующая событию Mw7,5 (длина L = 70 км, ширина W = 37,6 км, смещение D = 2,17 м, формулы 1, 2), центр которой соответствует положению эпицентра землетрясения, представленному в каталоге USGS (рис. 3). Остальные параметры заданы в соответствии с типичными значениями для зоны субдукции: направление простирания плоскости str = 180 °, угол падения dip = 15 ° и направление дислокации в очаге rake = 90 °.

Поскольку указанные смещения в реальности отягощены ошибками измерений, необходимо присвоить смоделированным сдвигам точность их определения. Величины ошибок зависят от разных факторов, наиболее значимым из них является способ обработки ГНСС-данных. К примеру, используемый системой REGARD программный пакет RTK-lib позволяет определить горизонтальные смещения с точностью 5 см, вертикальные – 8 см. Однако существуют методы, например VADASE [5], способные достичь точности 5 мм в горизонтальной плоскости и 1,5 см по вертикали. Принимая во внимание существующие возможности ГНСС-технологий, в работе мы использовали три варианта точности определения косейсмических смещений: 1) 2 см и 4 см (условие I), 2) 3 см и 5 см (условие II), 3) 5 см и 8 см (условие III) по горизонтали и вертикали соответственно. Применяя разные величины ошибок измерений при инверсии одного и того же набора косейсмических смещений, мы выяснили, каким образом точность определения смещений влияет на результат инверсии.

Результаты инверсии

В результате инверсии были получены три сейсмофокальные плоскости (рис. 2), соответствующие каждому условию точности ГНСС-наблюдений. Их параметры (h – глубина верхнего края, L – длина, W – ширина, str – направление простирания, dip – угол падения, D – смещение в очаге) и моментные магнитуды, вычисленные по формулам 1, 2, приведены в табл. 1.

Таблица 1

Параметры инвертированных плоскостей землетрясений

Условия

h, км

L, км

W, км

str, °

dip, °

D, м

Mw

I

17,76663

362,5566

151,0067

204,2511

16,21113

8,65518

8,78

II

18,547

343,8774

144,2041

201,7034

15,80829

8,99657

8,76

III

18,79407

188,2681

110,3264

185,9245

15,52989

42,08663

8,95

Согласно данным табл. 1 и рис. 2, а параметры инвертированных плоскостей I и II хорошо совпадают, их положение согласуется с положением подплоскостей референсной модели с максимальными смещениями. Однако положение плоскости III (рис. 2, б) не коррелирует с положением опорных плоскостей и плоскостей I и II. Таким образом, магнитуды всех инвертированных событий превосходят опорное значение.

nechaev2a.tif

а)

nechaev2a.tif

б)

Рис. 2. Смоделированные косейсмические смещения и инвертированные по ним плоскости очагов землетрясений. а) эталонная (красный контур) и инвертированные по I и II условиям (зеленый и синий прямоугольники соответственно) плоскости землетрясений. Эталонные (красные стрелки) и инвертированные по I и II условиям (зеленые и синие стрелки соответственно) косейсмические смещения; б) эталонная (красный контур) и инвертированная по III условию (синий прямоугольник) плоскости землетрясений. Эталонные (красные стрелки) и инвертированные по III условию (синие стрелки) косейсмические смещения

Очевидно, что третье условие точности привело к значительной переоценке как магнитуды, так и величины смещения в очаге. Тем не менее при всех условиях опорное значение магнитуды было достигнуто путем инверсии косейсмических значений региональных ГНСС-станций, что свидетельствует об уверенном преодолении проблемы недооценки магнитуды, свойственной действующей СПЦ, с помощью ГНСС-методов. Стоит отметить, что моментная магнитуда, соответствующая опорной модели (опорная магнитуда), согласно формулам 1, 2 равна Mwref 8,72. Магнитуда события, инвертированного по второму условию MwII 8,76, является наиболее приближенной к Mwref, несмотря на бо́льшие ошибки определений смещений, чем по условию I (3 и 2 см по горизонтали соответственно).

Расчет мгновенных смещений морского дна

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

Для этих целей в ближней зоне землетрясения были рассчитаны поля косейсмических смещений, соответствующие опорной и инвертированным плоскостям сейсморазрывов (рис. 2). Согласно рис. 2, а условия I и II показывают примерно одинаковое распределение и величины горизонтальных смещений, поэтому далее из этих двух вариантов рассматриваться будет только условие II как менее требовательное к возможностям ГНСС-методов. Инвертированное по второму условию и опорное поля смещений частично соответствует. Наилучшее согласие горизонтальных смещений наблюдается в центральной части эталонной модели в местах наибольших смещений в очаге, а также в южной части п-ова Камчатка и северных Курильских островов. Значительные расхождения приходятся на места, где опорная и инвертированная модели не перекрывают друг друга: в северной части поля и к юго-востоку от опорной модели.

Данное разногласие легко объясняется невозможностью точного повторения сложной референсной модели одной плоскостью. Вертикальные смещения в целом неплохо соответствуют (рис. 3, а). Максимальные смещения практически идентичны (табл. 2). Инвертированное поле не описывает локальный максимум северной части опорного поля и расходится по долготе с южным максимумом за счет разности направлений простирания эталонной и инвертированных плоскостей.

nechaev3a.tif nechaev3b.tif

а) б)

Рис. 3. Вертикальные косейсмические смещения: а) эталонные (красный контур) и инвертированные по II условию (синий контур), горизонтали проведены через 1 м; б) эталонные (красный контур, горизонтали проведены через 1 м) и инвертированные по III условию (синий контур, горизонтали проведены через 2 м)

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

Таблица 2

Максимальные горизонтальные и вертикальные косейсмические смещения, соответствующие опорной и инвертированным моделям очага

Модель землетрясения

Горизонтальные смещения, см

Вертикальные смещения, см

условие I

376,5

381,9

условие II

378,6

391,6

условие III

1483,6

1738,4

опорная модель

360,3

391,1

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

Заключение

Рассчитанные косейсмические смещения значительно превосходят по точности смещения, определенные ГНСС-методами, что делает их легко детектируемыми в режиме реального времени. Однако ошибки определения смещений значительным образом влияют на результаты инверсии. Расчеты показали, что точности определения в реальном времени косейсмических смещений в 3 см и 5 см в горизонтальной и вертикальной плоскостях достаточно для удовлетворительной оценки параметров землетрясения, а главное, величины смещения в очаге.

Следует отметить, что данная работа опирается на реально существующие и существовавшие ГНСС-станции, которые могут быть использованы в перспективе для модернизации существующей СПЦ. Несмотря на простоту модели очага и методологии, использованной для инверсии, косейсмические смещения морского дна, посчитанные по инвертированной (условие I и II) и референсной моделям, очень близки, что позволяет получить близкие к реальным параметры цунами в источнике.

Таким образом, можно уверенно заключить, что уже сейчас ГНСС-методы могут быть использованы для целей раннего предупреждения о цунами в Курило-Камчатском регионе и решения ныне существующей проблемы недооценки магнитуды при сильных землетрясениях. Мы подтвердили, что существующее состояние ГНСС-технологий и ГНСС-сетей в регионе способно значительно улучшить надежность работы нынешней СПЦ.

Данная работа выполнена при поддержке Японско-Российского центра молодежных обменов, а также Института прикладной математики Дальневосточного отделения Российской академии наук.


Библиографическая ссылка

Нечаев Г.В., Шестаков Н.В., Такахаши Х., Герасименко М.Д., Сысоев Д.В. О ПРИМЕНЕНИИ МЕТОДОВ ГЛОБАЛЬНЫХ НАВИГАЦИОННЫХ СПУТНИКОВЫХ СИСТЕМ ДЛЯ ЦЕЛЕЙ РАННЕГО ПРЕДУПРЕЖДЕНИЯ О ЦУНАМИ В КУРИЛО-КАМЧАТСКОМ РЕГИОНЕ // Успехи современного естествознания. – 2020. – № 3. – С. 92-98;
URL: http://www.natural-sciences.ru/ru/article/view?id=37351 (дата обращения: 20.10.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074