В последние годы активно развиваются системы предупреждения о цунами (СПЦ), основанные на использовании методов Глобальных навигационных спутниковых систем (ГНСС). В отличие от СПЦ, использующих сейсмические методы, системы нового типа нацелены на определение параметров очага землетрясения и, далее, оценку начальных движений морского дна и характеристик цунами. Подобный принцип используется в системах REGARD (Япония) и G-larmS и G-FAST (США) [1–3]. Схема работы этих систем схожа:
1) детектирование косейсмических смещений, порождённых землетрясением, ГНСС-методами (либо их комбинированием, например, с обработкой сейсмологических данных);
2) определение параметров землетрясения путём решения обратной задачи (инверсия косейсмических смещений);
3) моделирование смещений морского дна с использованием полученных параметров землетрясения;
4) расчёт характеристик цунами, заплесков волн на побережье и определение зоны затопления.
От системы требуется не только точность оценки цунамиопасности, но и оперативность её получения. Одним из факторов, влияющих на точность оценки параметров землетрясения и смещений морского дна, является используемая системой модель внутреннего строения и физических свойств Земли. Более подробные модели, участвующие в инверсии и моделировании смещений морского дна, теоретически должны способствовать более точной оценке параметров очага и цунами. Однако использование сложных моделей неизбежно приведёт к дополнительным временным затратам при расчётах и негативно отразится на быстродействии системы. Следовательно, необходимо соблюсти разумное соотношение скорости обработки данных и точности получаемых результатов. Даже в наиболее совершенной на данный момент СПЦ REGARD [1] модель Земли представлена однородным полупространством (инверсия происходит по методу Мацууры и Хасегавы [4]), т.е. не учитывает ни слоистости земной коры, ни сферичности её поверхности.
В работе [5] уже рассматривалось влияние слоистости и сферичности на величины горизонтальных косейсмических смещений, смоделированных с использованием параметров Суматра-Андаманского землетрясения Mw 9,2. В отличие от [5], в настоящей работе исследуется необходимость учёта слоистости и сферичности Земли в работе СПЦ, поскольку одним из ключевых аспектов её работы является оперативное моделирование косейсмических движений морского дна. Зависимость быстродействия гипотетической системы от задействованных моделей Земли рассматриваться не будет, поскольку здесь большую роль играют вычислительные возможности центра обработки системы.
Расчёт косейсмических смещений
Наиболее цунамиопасными регионами России являются Камчатка и Курильские острова, и наличие здесь СПЦ, основанной на ГНСС-методах, крайне желательно [6]. Для простоты рассмотрены только гипотетические землетрясения: три очага, с моментными магнитудами Mw 7,5; 8,0 и 8,5. Очаги локализованы в Камчатской зоне субдукции у берегов п-ова Камчатка, приблизительно в районе генерации землетрясения 1952 г. Mw 9,0. Данные события потенциально способны вызвать волны цунами различной интенсивности. Очаги представлены одной дислокационной плоскостью с равномерным распределением смещения. Длина L и ширина W плоскости дислокаций и величина смещения в очаге D были рассчитаны по формулам (1)–(4) [7].
, (1)
(2)
, (3)
. (4)
Полученные с использованием μ = 3,5×1010 Н/м параметры представлены в табл. 1.
Таблица 1
Параметры очагов гипотетических землетрясений
|
магнитуда |
L (км) |
W (км) |
D (м) |
|
Mw 7,5 |
70,00 |
37,58 |
2,16 |
|
Mw 8,0 |
131,82 |
63,20 |
3,85 |
|
Mw 8,5 |
248,31 |
105,92 |
6,85 |
Для всех плоскостей глубина залегания верхней грани равна 5 км, угол падения 15 ° и направление простирания сейсморазрыва 215 °, что в целом соответствует параметрам землетрясений Камчатской зоны субдукции.
Для расчёта косейсмических смещений использованы три модели Земли: I) однородное полупространство с μ = 3,5•1010 (программный пакет Coulomb3.0), II) сферическая однородная и III) сферическая слоистая модели (программный пакет Static1D). В качестве сферической слоистой была использована модель Preliminary Reference Earth Model (PREM), в то время как в сферической однородной земная кора представлена слоем толщиной 61 км с μ = 3,5•1010 Н/м, коэффициентом Пуассона ν = 0,25 и плотностью 3300 кг/м3.
Для всех очагов и моделей среды были рассчитаны соответствующие им вертикальные косейсмические смещения в узлах сетки с шагом 0,2 ° (рис. 1) и максимальные вертикальные сдвиги морского дна (табл. 2). Разность косейсмических смещений дна моря, рассчитанных с использованием сферической однородной и модели однородного полупространства, показывает влияние сферичности, а разность смещений, полученных с использованием слоистой сферической и однородной сферической моделей, отображает влияние слоистости земной коры.
Таблица 2
Максимальные вертикальные косейсмические смещения, соответствующие рассмотренным гипотетическим землетрясениям и моделям Земли
|
Магнитуда |
Максимальные вертикальные движения морского дна (см) |
СОМЗ ОПМЗ |
ССМЗ СОМЗ |
ССМЗ ОПМЗ |
||
|
Однородное полупространство (ОПМЗ) |
Сферическая однородная модель Земли (СОМЗ) |
Слоистая сферическая модель Земли (ССМЗ) |
||||
|
Mw 7,5 |
87,5 |
85,9 |
77,5 |
0,98 |
0,90 |
0,88 |
|
Mw 8,0 |
171,4 |
162,4 |
148,1 |
0,95 |
0,91 |
0,86 |
|
Mw 8,5 |
307,4 |
297,7 |
276,7 |
0,97 |
0,93 |
0,90 |
Рис. 1. Вертикальные косейсмические смещения, рассчитанные с использованием различных моделей Земли и магнитуд гипотетических землетрясений: а) Mw 7,5, б) Mw 8,0, в) Mw 8,5. Для Mw 7,5 изолинии проведены через 20 см, для Mw 8,0 и Mw 8,5 – через 50 см. Красной линией обозначены смещения, рассчитанные с помощью однородного полупространства, синей линией – по однородной сферической модели, зелёной – по слоистой сферической модели. Проекции плоскостей очагов отображены чёрными прямоугольниками, верхняя грань отображена прерывистой линией
В соответствии с табл. 2 и рис. 1 учёт слоистости и сферичности Земли выражается в получении меньших величин вертикальных смещений. Полученная разница достигает 30 см только при самом сильном рассматриваемом землетрясении, при котором дислокации морского дна составляют порядка 3 м.
Далее, были рассчитаны горизонтальные косейсмические смещения в местах расположения станций постоянных ГНСС-наблюдений на Дальнем Востоке России (рис. 2). Как и в случае с вертикальными, наименьшие смещения соответствуют слоистой сферической модели, несколько большей величины – однородной сферической, самые большие – однородному полупространству. Влияние слоистости и сферичности значительнее на горизонтальные смещения, чем на вертикальные. Также визуально виден различный вклад сферичности и слоистости в уменьшение величины смещений в зависимости от магнитуды события. Так, для очага Mw 7,5 эффект слоистости земной коры ведёт к более значимому преуменьшению смещений, чем сферичность поверхности. Для Mw 8,0, заметно примерно равное влияние факторов, а при Mw 8,5 слоистость практически не влияет на смещения, и к уменьшению их величин, согласно расчётам, ведёт сферичность поверхности. При этом следует заметить, что по отношению к очагам землетрясений рассматриваемые ГНСС-пункты принадлежат к разным зонам – ближней и дальней, в зависимости от магнитуды события. Для события Mw 7,5 все ГНСС-станции будут принадлежать дальней зоне, для Mw 8,0 в неё попадают только ближайшие станции, расположенные на побережье Камчатского п-ова, а для землетрясения Mw 8,5 ближняя зона охватывает уже все станции южной части Камчатки.
Рис. 2. Горизонтальные косейсмические смещения, рассчитанные на пунктах региональных ГНСС-сетей с использованием моделей – сферической слоистой (зелёные векторы), однородной сферической (синие векторы) и однородного полупространства (красные векторы). Использовались гипотетические очаги: а) Mw 7,5, б) Mw 8,0, в) Mw 8,5
Таким образом, можно предположить, что чем сильнее событие и ближе наблюдаемые пункты находятся к очагу, тем больший вклад в уменьшение смещений оказывает сферичность поверхности по сравнению со слоистостью. И наоборот, чем слабее очаг и чем дальше от него расположены пункты наблюдения, тем заметнее влияние слоистости земной коры. В любом случае, оба рассматриваемых фактора приводят к уменьшению косейсмических смещений, причём в большей степени этот эффект сказывается на горизонтальной компоненте сдвигов.
Инверсия
Слоистая сферическая модель намного точнее описывает действительное строение Земли, нежели однородное полупространство. Поэтому для дальнейших вычислений мы приняли косейсмические смещения, рассчитанные с помощью модели PREM в качестве реальных подвижек. Как было показано выше, использование однородного полупространства в качестве модели Земли при моделировании косейсмических смещений приводит к переоценке их значений. И наоборот, инвертирование реальных смещений без учёта слоистости и сферичности должно привести к недооценке магнитуды землетрясения и, соответственно, цунамиопасности.
Для наглядной демонстрации этого эффекта магнитуды гипотетических землетрясений были переопределены путём инверсии по методу [4] рассчитанных с использованием сферической слоистой модели косейсмических смещений. При инверсии пространственные характеристики очагов были зафиксированы, и определялась только величина смещения в сейсмофокальной плоскости (Dinverted). Косейсмическим смещениям присвоены ошибки величиной 1 мм, по аналогии с точностями реальных измерений. После определения Dinverted по формулам (1)–(4) были определены магнитуды инвертированных событий Mwinverted (табл. 3).
Таблица 3
Опорные и инвертированные значения магнитуд землетрясений и смещения в очаге
|
магнитуда |
Dreference (см) |
Dinverted(см) |
Dinverted / Dreference |
Mwinverted |
|
Mw 7,5 |
216,0 |
105,5 |
0,49 |
7,29 |
|
Mw 8,0 |
385,0 |
231,4 |
0,60 |
7,85 |
|
Mw 8,5 |
685,4 |
585,4 |
0,85 |
8,49 |
Инвертированные значения смещения в очаге Dinverted меньше опорных (исходных) Dreference. В абсолютных значениях разница варьируется от 1 м при Mw 8,5 до 1,53 м при Mw 8,0. Однако при относительном сравнении Dreference и Dinverted становится очевидным, что расхождение уменьшается с ростом магнитуды. При Mw 7,5 недооценка составляет около 50 %, при Mw 8,5–15 % от опорных значений. Для большей наглядности были смоделированы поля вертикальных косейсмических смещений (инвертированные поля вертикальных смещений) с использованием инвертированных параметров землетрясений и модели Земли PREM (рис. 3). Максимальные косейсмические смещения согласно опорной dispreference и инвертированной dispinverted моделям землетрясений представлены в табл. 4. Относительное соотношение dispmax inverted и dispmaxreference идентично отношениям Dreference и Dinverted. При этом абсолютное максимальное расхождение в 60 см наблюдается при Mw 8.0.
Таблица 4
Опорные и инвертированные максимальные вертикальные смещения морского дна
|
Магнитуда |
Максимальные вертикальные движения морского дна (см) |
dispinverted (однородное полупространство)/dispreference |
||
|
dispreference |
dispinverted (PREM) |
dispinverted (однородное полупространство) |
||
|
Mw 7,5 |
77,5 |
37,9 |
42,7 |
0,55 |
|
Mw 8,0 |
148,1 |
89,0 |
103,0 |
0,70 |
|
Mw 8,5 |
276,7 |
236,4 |
262,5 |
0,95 |
С другой стороны, если в СПЦ на всех этапах обработки адаптировано однородное полупространство в качестве модели Земли, то недооценка дислокации в очаге может быть частично компенсирована переоценкой смещений морского дна в ходе последующего моделирования. В связи с этим предположением также были рассчитаны поля смещений морского дна с использованием инвертированных параметров землетрясений и однородного полупространства как модели Земли (табл. 4, рис. 3). В данном случае, действительно, инвертированные поля имеют меньшее расхождение с опорными. Для Mw 8,0 и Mw 8,5 улучшение составило 10 % от опорных значений по сравнению с применением слоистой сферической модели. Максимальная недооценка наблюдается при Mw 8,0 и равна 45 см. Для события Mw 8,5 недооценка минимальна как в абсолютном (14 см), так и в относительном (5 %) отношениях.
Рис. 3. Вертикальные поля смещений, соответствующие опорным (зелёные изолинии) и инвертированным очагам (красные изолинии – смещения рассчитаны в однородном полупространстве, синие изолинии – с использованием слоистой сферической модели Земли). Магнитуды гипотетических землетрясений: а) Mw 7,5, б) Mw 8,0, в) Mw 8,5
Для проверки применённой методологии, косейсмические смещения, рассчитанные в однородном полупространстве на ГНСС-станциях, были инвертированы, и полученные магнитуды достигли значений исходных, что говорит о корректности применяемого подхода.
Заключение
На примере трёх гипотетических землетрясений различной мощности было рассмотрено влияние используемой модели Земли на величины моделируемых косейсмических смещений, а также на результаты инверсии. Показано, что модель Земли существенным образом влияет на величины смещений, причём данное влияние более значимо по горизонтальной компоненте. Кроме того, учёт сферичности и слоистости влияет на смоделированные величины смещений неравномерно и варьируется в зависимости от силы землетрясения и дальности расположения наблюдаемых пунктов от очага. Тем не менее некоторые закономерности очевидны: влияние слоистости коры усиливается при удалении от очага, а сферичности поверхности – при увеличении магнитуды землетрясения и приближении к нему. Игнорирование обоих факторов ведёт к переоценке величин смещений в ходе моделирования.
Для целей СПЦ важно хорошее соответствие реальных смещений морского дна (в нашем случае – опорных смещений) с инвертированными. При использовании однородного полупространства, несмотря на некоторую недооценку величины дислокации в очаге и, соответственно, магнитуды землетрясения, недооценка инвертированных смещений в худшем случае не превышает 50 см. Наилучшее соответствие получено при наиболее сильном событии Mw 8,5.
Если учесть, что одноплоскостная модель землетрясения не способна точным образом описать реально возникшие смещения в очаге землетрясения, то полученные расхождения можно считать приемлемыми, по крайней мере для наиболее сильных землетрясений. В случае событий Mw ≤ 8,0 недооценка максимальных смещений морского дна достигает половины от опорных значений, что может отразиться на выявлении самого факта цунами при его моделировании.
Исходя из вышесказанного, можно заключить, что в случае сильных землетрясений Mw ≤ 8,5, происходящих в зоне субдукции на небольшом удалении от побережья и ГНСС-станций, пренебрежение слоистостью и сферичностью не является критическим. С другой стороны, для событий Mw ≤ 8,0 всё же желательно использование более точной модели Земли.
Таким образом, при возникновении сильных землетрясений задействование наиболее простой модели Земли (однородное полупространство) в СПЦ способствует довольно точному и наиболее быстрому на данный момент определению начальных параметров волны цунами.
Работа выполнена при поддержке Японско-Российского центра молодёжных обменов, а также Института прикладной математики Дальневосточного отделения Российской академии наук.