Сейсмологи научились точечно вычислять места сильнейших землетрясений

30 июля 2025 года у восточного побережья Камчатки произошло одно из мощнейших землетрясений в истории инструментальных наблюдений. Его моментная магнитуда составила 8,84. В этом регионе Тихоокеанская тектоническая плита непрерывно движется на северо-запад и уходит в мантию под Охотоморскую плиту со скоростью около 80 миллиметров в год. Этот процесс называется субдукцией.

Катастрофа 2025 года произошла практически на том же участке границы плит, где в ноябре 1952 года случилось другое гигантское землетрясение с магнитудой 9,0. Однако детальный анализ сейсмических волн и смещений земной поверхности сразу выявил два принципиальных противоречия.

Во-первых, не сходился баланс механических смещений. За 73 года, прошедших между событиями 1952 и 2025 годов, плиты при скорости 80 миллиметров в год должны были накопить взаимный сдвиг величиной около 5,8 метра. Это смещение соответствует магнитуде не более 8,7. Однако расчеты показали: во время главного толчка 2025 года горные породы внутри очага сместились относительно друг друга на 20, а местами на 30 метров.

Во-вторых, землетрясение почти не проявило себя на побережьях океана. В 1952 году аналогичный по силе подземный толчок породил волны цунами высотой от 15 до 18 метров, которые полностью разрушили город Северо-Курильск. В 2025 году глубоководные океанические датчики давления зафиксировали волны высотой всего 0,85 метра, а на побережье Гавайских островов высота подъема воды не превысила двух метров.

Геофизики Аксель Периолла и Гарет Фаннинг из Калифорнийского университета в Риверсайде провели исследование, опубликованное в научном журнале Geophysical Research Letters. Они применили метод численного моделирования, чтобы восстановить реальную геометрию контакта плит и объяснить физические причины обоих феноменов.

Физика контакта плит: скольжение и сцепление

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

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

В других зонах трение между плитами слишком велико. Шероховатости, выступы и повышенное давление заставляют два массива горных пород намертво сцепляться. Смещение плит в этом месте временно прекращается. Возникает заблокированная зона (в сейсмологии ее называют «асперити»).

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

Дополнительную сложность процессу придает эффект тени напряжений. Когда прочный участок заблокирован, его механическое сопротивление принудительно тормозит движение соседних, незаблокированных зон. В результате даже участки со слабым трением замедляют свой ход, если находятся рядом с крупным заблокированным массивом.

Математический поиск без субъективных допущений

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

Периолла и Фаннинг исключили влияние человеческого фактора. Они объединили расчеты методом граничных элементов с алгоритмом статистической выборки Метрополиса — Гастингса.

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

Затем авторы использовали архив измерений камчатской сети высокоточных спутниковых станций геодезического мониторинга (GNSS) за межсейсмический период с 1952 по 2025 год. Эти приборы непрерывно фиксировали, с какой скоростью и в каком направлении перемещаются разные участки суши под действием давления океанической плиты.

Алгоритм работал следующим образом:

1. Создавался случайный вариант распределения зон сцепления по всей площади разлома.
2. Программа рассчитывала, какие смещения на поверхности суши должна вызвать такая конфигурация с учетом упругости горных пород и эффекта тени напряжений.
3. Расчетные смещения сравнивались с реальными данными наземных приборов GNSS.
4. Если случайное изменение приближало расчет к реальности, оно сохранялось. Если расхождение увеличивалось, вариант отбрасывался.

Ученые запустили один миллион независимых расчетных цепочек, каждая из которых содержала по 10 000 последовательных шагов. Усреднив результаты миллионов финальных конфигураций, исследователи получили объективную карту вероятности сцепления для каждого квадратного километра зоны субдукции.

Архитектура очага и физический смысл форшоков

Результаты расчетов показали, что вдоль восточного побережья Камчатки существует единая устойчивая область сцепления площадью около 100 000 квадратных километров. Вероятность блокировки контакта здесь устойчиво превышает 45%, а статистическая неопределенность близка к нулю.

Верхний край этой области залегает на глубине около 20 километров, а нижний опускается до 68 километров. Ниже отметки в 68 километров вероятность сцепления падает практически до нуля. Это согласуется с физикой процессов внутри Земли: на таких глубинах температура превышает предел хрупкости силикатных пород, и контакт плит переходит в состояние постоянного пластического течения без накопления упругих напряжений.

Когда исследователи разделили эту область на отдельные сегменты, они обнаружили прямую связь между зонами наибольшего трения и точками зарождения землетрясений. Внутри структуры выделились компактные ядра, где вероятность сцепления достигала 88-98%.

Оказалось, что гипоцентр (точка начала разрыва) землетрясения 1952 года на глубине 21 километр точно совпал с северным ядром сцепления. Гипоцентр землетрясения 2025 года, расположенный на глубине 35 километров, совпал со смежным центральным ядром. Разрушение коры в обоих случаях начиналось именно там, где плиты удерживали друг друга сильнее всего.

Кроме того, модель объяснила механику предварительных толчков — форшоков. За одиннадцать месяцев до катастрофы, в августе 2024 года, в регионе произошло землетрясение магнитудой 7,1. А всего за несколько дней до основного удара в июле 2025 года последовал второй форшок магнитудой 7,4.

Оба этих толчка произошли точно по внешним границам заблокированного кластера. Согласно закону трения горных пород Дитериха, максимальная скорость концентрации сдвиговых механических напряжений возникает не в центре сцепленного участка, а на его геометрических кромках — там, где неподвижная часть контакта граничит с непрерывно ползущей породой.

Перед главным ударом 2025 года на границах зоны сцепления происходили локальные срывы: разрушались краевые удерживающие участки. Когда эти опоры разрушились, механическое напряжение перераспределилось на центральное ядро зоны контакта, что привело к мгновенному распространению магистрального разрыва.

Гравитационные границы: почему разрыв остановился

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

Для этого авторы сопоставили геометрию разрыва с данными спутниковой гравиметрии — аномалиями силы тяжести в редукции Буге, ориентированными параллельно глубоководному желобу (TPGA). Изменения гравитационного поля отражают плотность пород в земной коре: участки с высокой плотностью создают положительные гравитационные аномалии, а менее плотные и трещиноватые породы — отрицательные.

Сопоставление показало: северная граница разрыва как 1952, так и 2025 года зафиксирована ровно в том месте, где начинается зона сильной положительной гравитационной аномалии со значениями более 50 миллигал.

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

Откуда взялись тридцать метров смещения

Определив границы очага, авторы вернулись к первому парадоксу: как плиты смогли сместиться на 20-30 метров, если за 73 года между землетрясениями накопилось всего 5,8 метра дефицита?

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

1. Неполная разгрузка в предыдущих циклах. Землетрясение 1952 года сняло лишь часть упругой энергии, накопленной за предыдущие столетия. Перед событием 1952 года сильнейший толчок на Камчатке происходил в 1737 году. За те 215 лет дефицит взаимного перемещения плит достиг величины, эквивалентной магнитуде 9,1-9,3. Значительная часть этой упругой энергии осталась в породах после 1952 года и суммировалась со смещением, накопившимся к 2025 году.
2. Постсейсмическое ускорение. Первые годы и десятилетия после сильного землетрясения глубокие слои коры и верхняя мантия деформируются с повышенной скоростью за счет релаксации напряжений. Этот процесс заставляет заблокированный участок коры накапливать деформации значительно быстрее средних 80 миллиметров в год.

Расчеты показали, что слабые и умеренные землетрясения, происходившие в регионе между 1952 и 2025 годами, высвободили незначительную долю общей энергии. Почти весь объем деформаций оставался законсервированным в породах и высвободился концентрированно во время события 2025 года.

Почему не возникло разрушительное цунами

Второй парадокс — умеренная высота волн цунами при магнитуде 8,84 — объясняется тем, на какой именно глубине происходило основное смещение пород.

Волна цунами образуется не от колебаний грунта, а от вертикального вытеснения объема воды. Чтобы в океане возникла опасная волна, дно должно резко приподняться или опуститься на большой площади. Эффективность этого процесса напрямую зависит от глубины залегания разрыва:

• Сценарий 1952 года (приповерхностный сдвиг): разрыв начался на глубине 21 километр и распространился вверх по плоскости разлома, выйдя непосредственно на прижелобный мелководный участок контакта плит. Дно океана на краю континентального склона резко сместилось по вертикали на многие метры прямо под толщей воды, вытолкнув огромную массу жидкости к побережью.
• Сценарий 2025 года (глубокий сдвиг): гипоцентр находился глубже — на отметке 35 километров. Главное смещение пород сосредоточилось на глубинах от 30 до 50 километров и развивалось преимущественно вниз и вбок по простиранию плиты.

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

Спутник океанографической альтиметрии SWOT (Surface Water and Ocean Topography) зафиксировал форму водной поверхности над очагом землетрясения 2025 года из космоса. Радары подтвердили, что вблизи глубоководного желоба возникли короткие дисперсионные волны небольшой амплитуды. Приповерхностная часть разлома испытала лишь остаточное смещение: основная кинетическая энергия землетрясения осталась запертой на глубине.

Слепая зона современной геофизики

Исследование Периоллы и Фаннинга доказало: мегаземлетрясения подчиняются закономерностям, которые можно рассчитывать математически. Участки мощного трения плит сохраняют свое положение на протяжении сотен лет и определяют координаты зарождения будущих катастроф, а плотностные барьеры коры задают предельные границы разрушений.

В то же время работа выявила фундаментальную инструментальную проблему. Все современные расчеты распределения напряжений на Камчатке строятся на данных станций GNSS, установленных на суше. Однако сухопутная часть полуострова удалена от глубоководного желоба на расстояние от 100 до 200 километров.

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

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

Для решения этой задачи требуется переход к донной геодезии. Установка автономных акустических станций и датчиков гидростатического давления непосредственно на дне Тихого океана — единственный способ напрямую измерить степень сцепления плит у кромки желоба и точно прогнозировать последствия неизбежных тектонических срывов.

Источник: Geophysical Research Letters