Полное описание геофизических каталогов, пространственной дискретизации Uber H3, физико-математических формул напряжения земной коры и архитектуры калиброванного ансамбля машинного обучения EarthquakeLive.org.
Иерархическая гексагональная сетка для геостатистического анализа без сферических искажений
Традиционные прямоугольные сетки (широта/долгота) страдают от критических искажений площади при приближении к полюсам и имеют два типа соседства (по ребру и по диагонали). В EarthquakeLive.org вся поверхность планеты дискретизируется с помощью дискретной глобальной сетки Uber H3.
Диаметр ячейки: ~110 км
Площадь: ~11 000 км²
Региональный сейсмический пояс
Диаметр ячейки: ~40 км
Площадь: ~2 500 км²
Локальные разломные сегменты
Каждая ячейка аккумулирует локальные исторические толчки, плотность разломов, расчетное b-value и приливные параметры.
Многомерные геофизические индикаторы напряжений и кинематики земной коры
В норме в спокойной земной коре b ≈ 1.0 (на одно событие M=5 приходится 10 событий M=4 и 100 событий M=3). Когда тектоническое сдвиговое напряжение в разломной зоне критически возрастает, происходит эффект «запирания» шероховатостей (asperities). При этом количество мелких микротолчков падает, и значение b опускается ниже 0.8.
Сейсмическая энергия растет как E ~ 10^(1.5 · M). Упругая деформация пропорциональна корню из энергии: √E ~ 10^(0.75 · M). График накопления деформации позволяет выявлять сейсмическое затишье (quiescence) и последующее ускорение высвобождения момента (Accelerating Moment Release — AMR).
Отношение краткосрочной недельной активности к 30-дневному фону. Значения > 2.0 сигнализируют об активизации разлома, форшоковой последовательности или сейсмическом рое.
Отклонения наклона частотно-магнитудного распределения, накопление упругой деформации, ускорение форшоковых роёв и длительное затишье в запертых сегментах разлома.
Геодезическая близость к структурам глобальной базы активных разломов GEM, исторические скорости тектонического смещения и категориальная кинематика разломов (субдукция, сдвиг, надвиг).
Астрономическое эфемеридное моделирование лунного и солнечного приливного потенциала, циклы сизигийного выравнивания (сизигийные и квадратурные приливы) и фаза лунного перигея.
Профилирование глубины гипоцентра (мелкофокусные vs глубокие зоны Беньоффа-Вадати) в сочетании с планетарными геомагнитными индексами возмущения NOAA SWPC (Kp).
Землетрясения в зонах субдукции с малой глубиной очага (h < 50 км) и высокой магнитудой (M ≥ 7.0) вызывают резкое вертикальное косейсмическое смещение морского дна, передавая импульс толще воды. В глубоком открытом океане (средняя глубина d ≈ 4000 м) цунами распространяются по уравнению мелкой воды с фазовой скоростью c = √(g·d) ≈ 200 м/с (~720 км/ч). Прием телеметрии в реальном времени от центров NOAA PTWC / NTWC через протокол CAP XML v1.2 формирует концентрические изохроны добегания волн на 15 мин, 30 мин и 1 час с учетом прибрежного усиления по закону Грина (H₂/H₁ = (d₁/d₂)^(1/4)).
Подъем магмы через хрупкую кору вызывает микротрещиноватость и рои мелкофокусных землетрясений перед извержением. EarthquakeLive сопоставляет каждое событие с каталогом 1 214 голоценовых вулканов Смитсоновского института (GVP) в радиусе 25 км от кратера и на глубине до 20 км, сравнивая плотность событий с 30-дневным базовым уровнем для надежного отличия магматической активности от обычных тектонических афтершоков.
Шлейфы вулканического пепла представляют смертельную угрозу для турбин реактивной авиации из-за плавления силикатного стекла на лопатках компрессора. EarthquakeLive непрерывно опрашивает 9 международных консультативных центров ICAO VAAC, извлекая пространственные полигоны по высотным эшелонам вплоть до FL600+ (60 000 футов / 18 км) с прогнозом дрейфа на 0ч, +6ч, +12ч и +18ч.
Калиброванный ансамбль моделей машинного обучения на прогнозных горизонтах 12ч, 24ч, 72ч и 7 дней
Окно непосредственной угрозы, форшоковых всплесков и экстренного реагирования
Оперативная оценка непосредственной угрозы и афтершоков
Основной оперативный горизонт для систем оповещения
Среднесрочный тренд накопления деформации
| Компонент | Алгоритм | Целевая переменная | Оптимизация |
|---|---|---|---|
| Вероятностный классификатор | Градиентный бустинг решающих деревьев | Бинарное событие: M ≥ 4.5 | Взвешенные классы, logloss + Brier Loss |
| Регрессор магнитуды | Робастная регрессия магнитуды | Ожидаемая M_max в ячейке H3 | Huber Loss (устойчивость к выбросам) |
| Пост-калибратор вероятностей | Монотонный эмпирический калибратор | P_calibrated ∈ [0, 1] | Минимизация Brier Score, калибровка Reliability |
Сырые предсказания градиентного бустинга часто переоценивают или недооценивают экстремальные вероятности. Изотоническая калибровка гарантирует: если модель сообщает вероятность 40%, то ровно в 40 из 100 таких прогнозных окон землетрясение действительно происходит на практике.
Теоретико-игровое разложение прогноза каждой ячейки на физические факторы
Для каждой гексагональной ячейки H3 на карте EarthquakeLive.org рассчитывает вектор Шепли с помощью алгоритма Lundberg & Lee TreeSHAP. Вы можете нажать на любую ячейку карты и увидеть точный процентный вклад каждого фактора в повышение или снижение риска.
Сейсмические каталоги обновляются непрерывно. Модели автоматически дообучаются по скользящему временному окну с верификацией по метрикам ROC-AUC и Brier Skill Score (BSS).
0.864
0.078
Непрерывная автоматическая оценка относительно каталогов землетрясений по протоколу walk-forward.
Все прогнозы, опубликованные на EarthquakeLive, навсегда фиксируются по времени и проверяются по авторитетным каталогам (USGS ComCat и EMSC) после истечения каждого горизонта прогноза. Изотоническая калибровка вероятностей гарантирует, что если региону присвоена вероятность 60%, событие происходит в 60 из 100 таких прогнозных окон.
Инстансинг геометрии Deck.gl v9, нативный векторный фолбэк MapLibre, ограничение видеопамяти DPR 1.0 и сферические траектории камеры с безопасным пересечением антимеридиана
Для мощных браузеров на ПК и смартфонах EarthquakeLive использует Deck.gl v9 с 3D-инстансингом гексагональных призм (PolygonLayer) и пульсирующими эпицентральными кольцами (ScatterplotLayer) при 60 FPS. При работе на Smart TV с жесткими ограничениями памяти (Samsung Tizen WebAPI, LG webOS, Android TV WebView) или слабом оборудовании система динамически переключается на нативный векторный рендеринг MapLibre через проверку shouldUseNativeMapLibre.
Расход видеопамяти строго контролируется ограничением devicePixelRatio до 1.0 на экранах 4K/8K, что снижает объём буфера кадра GPU до 75% без заметной потери четкости при просмотре с расстояния 3 метров.
10-футовый интерфейс полностью управляется физическим пультом ДУ телевизора (ВВЕРХ, ВНИЗ, ВЛЕВО, ВПРАВО, OK, НАЗАД, ВЫХОД). На Samsung Smart TV модуль tvTizenBridge.ts регистрирует коды клавиш в Samsung Tizen WebAPI. На пультах LG webOS и Android TV используются нативные обработчики клавиатурных событий.
Автономный тур автопилота циклически переключается между главными мировыми сейсмокластерами с помощью модуля CameraFlightController. Для предотвращения резкого разворота камеры на 360° при пересечении антимеридиана 180° (например, между Японией и Аляской) траектория рассчитывается через нормализованную угловую дельту:
Траектории выполняются с плавным кубическим сглаживанием, динамическим наклоном высоты и автоматическими паузами, адаптированными для общественных экранов.
При поступлении данных о сильном землетрясении (M ≥ 5.5) от любого из 18 каталогов тур автопилота и экранная заставка мгновенно прерываются. Камера плавно переводится на новый эпицентр, выводится критическая карточка события и включается звуковая сирена.
Строгие теоретические основы, опирающиеся на мировую сейсмологическую литературу
"Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a - bM and its confidence limits. Bulletin of the Earthquake Research Institute, Tokyo University, 43, 237–239."
Роль в системе: Оценивает наклон распределения повторяемости магнитуд на ячейках Uber H3 Res 4; выявляет аномалии концентрации сдвигового напряжения (b < 0.8).
"Gutenberg, B., & Richter, C. F. (1944). Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4), 185–188."
Роль в системе: Фундаментальное физическое соотношение частоты и магнитуды землетрясений в региональных каталогах.
"Benioff, H. (1951). Earthquakes and rock creep (Part I: Creep characteristics of rocks and the origin of aftershocks). Bulletin of the Seismological Society of America, 41(1), 31–62."
Роль в системе: Рассчитывает суммарную акустическую энергию и деформационную отдачу S = ∑10^(0.75 M_i) для фиксации сейсмического затишья или роевой активизации.
"Wald, D. J., Jaiswal, K. S., Allen, T. I., & Hearne, M. G. (2008). USGS Prompt Assessment of Global Earthquakes for Response (PAGER). U.S. Geological Survey Fact Sheet 2008-3093."
Роль в системе: Лежит в основе калиброванных логнормальных функций оценки разрушений зданий и вероятных потерь.
"Esteva, L. (1974). Geology and predicting extreme ground motion for earthquake-resistant design. In Proceedings of the 2nd International Conference on Seismic Risk."
Роль в системе: Расчет радиусов изосейсм для пространственного картирования интенсивности сотрясений (MMI V–VIII).
"Lundberg, S. M., & Lee, S.-I. (2017). A Unified Approach to Interpreting Model Predictions. Advances in Neural Information Processing Systems (NeurIPS 30), 4765–4774."
Роль в системе: Раскладывает прогнозные оценки ансамбля на точные значения Шепли, обеспечивая 100% прозрачность моделей.
"Bird, P. (2003). An updated digital model of plate boundaries. Geochemistry, Geophysics, Geosystems, 4(3), 1027."
Роль в системе: Расчет расстояния до границ тектонических плит, скоростей их относительного движения и кинематики конвергенции.
"Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell, Richmond, VA."
Роль в системе: Вычисление лунно-солнечного гравитационного приливного потенциала и циклов сизигийного выравнивания.
"Titov, V. V., & Synolakis, C. E. (1998). Numerical modeling of tidal wave runup. Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4), 157–171."
Роль в системе: Underpins NOAA MOST tsunami simulation and deep-ocean DART buoy wave inversion models.
"Kajiura, K. (1970). Tsunami source, energy and the ratio of wave energy to seismic energy. Bulletin of the Earthquake Research Institute, Tokyo University, 48, 835–869."
Роль в системе: Formulates the ratio of seismic radiated energy to tsunami potential energy in shallow subduction megathrust ruptures.
"OASIS Standard. (2010). Common Alerting Protocol Version 1.2. OASIS Open Standard."
Роль в системе: International interoperable XML standard for real-time ingestion of NOAA PTWC/NTWC multi-agency hazard bulletins.
Строгая научная воспроизводимость в соответствии с международными стандартами сейсмологии
100% первичных сейсмологических данных и параметров событий поступают в реальном времени из открытых государственных институтов: USGS, EMSC FDSNWS, ISC, BMKG, GeoNet NZ, INGV и GA. Исходные данные никогда не модифицируются и не фабрикуются.
Все модели машинного обучения оцениваются на замороженных тестовых фолдах без заглядывания в будущее (стандарт CSEP). Полный реестр верификации с Brier Skill Score (+0.354) и ROC-AUC (0.814) общедоступен через открытый API метрик.
Мы отвергаем антинаучные заявления о точном детерминированном предсказании. Зарождение очага подчиняется нелинейной механике сплошных сред. EarthquakeLive рассчитывает строго откалиброванные вероятности P(M ≥ 4.5) для поддержки готовности к ЧС.
Верифицированные научные, операционные и экстренные протоколы
Перейдите на интерактивную карту с наложением тектонических разломов, очагов землетрясений и сетки риска в реальном времени.