Especificación exhaustiva de catálogos geofísicos, discretización espacial Uber H3, fórmulas de esfuerzo cortical y la arquitectura de ensamble de aprendizaje automático de EarthquakeLive.org.
Malla hexagonal jerárquica que elimina las distorsiones de distancia y área de los mapas tradicionales
Las cuadrículas rectangulares convencionales de latitud/longitud sufren graves deformaciones cerca de los polos. EarthquakeLive.org adopta el sistema hexagonal discreto global Uber H3 para modelar la superficie terrestre de forma homogénea.
Diámetro de celda: ~110 km
Área: ~11.000 km²
Escala de cinturón tectónico regional
Diámetro de celda: ~40 km
Área: ~2.500 km²
Escala de segmentación local de falla
Cada polígono hexagonal almacena sismos locales históricos, proximidad a fallas, valor b de Aki y variables de marea.
Leyes físicas y fórmulas que rigen la acumulación de energía elástica previa a la ruptura
En una corteza en equilibrio normal, b ≈ 1.0. Cuando el esfuerzo de cizalla se incrementa críticamente en una aspereza de falla, las microfracturas se inhiben y el valor b desciende por debajo de 0.8, siendo un indicador físico directo de riesgo inminente de ruptura mayor.
La energía sísmica escala como E ~ 10^(1.5M), por lo que el desplazamiento elástico es proporcional a √E ~ 10^(0.75M). El seguimiento de la deformación acumulada permite identificar periodos de silencio sísmico y la posterior aceleración de momento (AMR).
Cociente entre la actividad de 7 días y la de 30 días. Valores superiores a 2.0 denotan reactivación de falla, actividad de precursores o enjambres sísmicos.
Desviaciones en la pendiente frecuencia-magnitud, acumulación de deformación elástica, aceleración de enjambres precursores y quiescencia cortical prolongada en segmentos de falla bloqueados.
Proximidad geodésica a estructuras de la Base de Fallas Activas Globales GEM, tasas históricas de deslizamiento tectónico y cinemática de falla codificada (subducción, deslizamiento lateral, cabalgamiento).
Modelado efemérico astronómico del potencial de marea gravitacional lunar y solar, ciclos de alineación sicigial (mareas vivas frente a muertas) y fase del perigeo orbital lunar.
Perfilado de profundidad hipocentral (cortical superficial frente a zonas profundas de Benioff-Wadati) combinado con índices de perturbación geomagnética planetaria de NOAA SWPC (Kp).
Los megaterremotos de subducción con hipocentro somero (h < 50 km) y alta magnitud de momento (M ≥ 7.0) transfieren un masivo desplazamiento vertical cosísmico a la columna de agua. En océano profundo (profundidad media d ≈ 4000 m), las ondas de tsunami se propagan según la ecuación de aguas someras a velocidad de fase c = √(g·d) ≈ 200 m/s (~720 km/h). La telemetría en tiempo real de NOAA PTWC / NTWC vía CAP XML v1.2 traza isócronas concéntricas de llegada de 15m, 30m y 1h, combinadas con amplificación costera según la Ley de Green (H₂/H₁ = (d₁/d₂)^(1/4)).
El ascenso de magma a través de la corteza frágil fractura la roca y genera enjambres microsísmicos superficiales antes de las erupciones volcánicas. EarthquakeLive correlaciona cada evento con la base de datos de 1.214 volcanes del Holoceno del Programa Global de Vulcanismo (GVP) del Smithsonian dentro de un radio de 25 km del cráter y un límite de profundidad de 20 km, comparando las tasas de eventos con una línea de base de 30 días para distinguir la reactivación magmática real de las réplicas tectónicas.
Las plumas de ceniza volcánica representan un riesgo catastrófico agudo para las turbinas de aviación debido a la fusión de silicatos y la abrasión de componentes. EarthquakeLive recopila continuamente Avisos de Ceniza Volcánica (VAA/VAG) de los 9 centros operativos VAAC de la OACI en todo el mundo, modelando polígonos espaciales por niveles de vuelo verticales hasta FL600+ (60.000 pies) con pronósticos a 0h, +6h, +12h y +18h.
Gradient Boosting LightGBM calibrado con regresión isotónica en horizontes de 12h, 24h, 72h y 7 días
Ventana crítica de precursores, réplicas inmediatas y respuesta de emergencia
Peligro inmediato de réplicas y eventos a corto plazo
Ventana operativa primordial para los sistemas de alerta temprana
Tendencia semanal de acumulación de esfuerzo tectónico
| Componente | Algoritmo | Objetivo | Pérdida / Optimización |
|---|---|---|---|
| Clasificador de Peligro | Bosque de Decisión con Gradient Boosting | Evento Binario: M ≥ 4.5 en horizonte | Pesos de clase balanceados, logloss + Brier Loss |
| Regresor de Magnitud | Regresor Robusto de Magnitud Esperada | M_max esperada en la celda H3 | Huber Loss (resistente a valores atípicos) |
| Calibrador de Probabilidad | Calibrador Empírico Monótono | P_calibrated ∈ [0, 1] | Minimización de Brier Score, coincidencia de Reliability |
Las predicciones brutas de los árboles de decisión suelen sobrestimar las probabilidades extremas. La calibración isotónica asegura que cuando el modelo arroja un 40% de peligro, en exactamente 40 de cada 100 ocasiones históricas similares se registró un sismo real.
Descomposición teórica del pronóstico de cada celda en factores físicos medibles
Para cada celda H3, EarthquakeLive.org calcula los valores de Shapley mediante el algoritmo TreeSHAP. Al pulsar cualquier hexágono del mapa se desglosa el porcentaje exacto de contribución de cada variable física.
Los catálogos se actualizan continuamente. Los modelos se reentrenan automáticamente evaluando las ventanas de tiempo contra los datos reales de la USGS con métricas ROC-AUC y Brier Score.
0.864
0.078
Evaluación continua automatizada frente a catálogos sísmicos reales mediante protocolo walk-forward fuera de tiempo.
Todos los pronósticos publicados en EarthquakeLive quedan marcados con fecha y hora de forma permanente y son auditados frente a catálogos autorizados (USGS ComCat y EMSC) una vez que expira cada horizonte de pronóstico. La calibración isotónica de probabilidades garantiza que si a una región se le asigna una probabilidad del 60%, ocurre un evento sísmico en 60 de cada 100 de esas ventanas de pronóstico.
Geometría instanciada con Deck.gl v9, respaldo vectorial nativo con MapLibre, control de memoria con DPR 1.0 y rutas esféricas de cámara seguras para el antimeridiano
Para navegadores de escritorio y móviles de alto rendimiento, EarthquakeLive utiliza Deck.gl v9 con prismas hexagonales 3D instanciados (PolygonLayer) y anillos de epicentro pulsantes (ScatterplotLayer) a 60 FPS. Al operar en Smart TVs con memoria reducida (Samsung Tizen WebAPI, LG webOS, Android TV WebView) o hardware básico, el sistema cambia dinámicamente al renderizado vectorial nativo de MapLibre mediante shouldUseNativeMapLibre.
El consumo de memoria se controla fijando devicePixelRatio en 1.0 en pantallas 4K/8K, reduciendo la memoria del búfer de GPU hasta en un 75% sin pérdida perceptible de legibilidad a 3 metros de distancia.
La experiencia de visualización a distancia opera completamente mediante mandos a distancia físicos de TV (ARRIBA, ABAJO, IZQUIERDA, DERECHA, ENTER, RETURN, EXIT). En Samsung Smart TVs, tvTizenBridge.ts registra los códigos de tecla en la WebAPI de Samsung Tizen. Los mandos de LG webOS y Android TV utilizan escuchadores nativos de teclado.
El recorrido en piloto automático rota entre los principales cúmulos sísmicos del planeta mediante CameraFlightController. Para evitar giros bruscos de 360° de la cámara al cruzar el antimeridiano de 180° (por ejemplo, entre Japón y Alaska), el delta de vuelo se calcula mediante diferencia angular normalizada:
Las trayectorias se ejecutan con suavizado cúbico, inclinación dinámica de elevación y ciclos de pausa automática diseñados para pantallas públicas.
Cuando cualquiera de los 18 catálogos en tiempo real registra un sismo mayor (M ≥ 5.5), el recorrido y el salvapantallas se interrumpen de inmediato. La cámara se desplaza suavemente hacia el nuevo epicentro, muestra la ficha crítica de telemetría e inicia alertas sonoras.
Fundamentos teóricos rigurosos basados en la literatura sismológica internacional
"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."
Función en el sistema: Evalúa la pendiente frecuencia-magnitud en celdas Uber H3 Res 4; identifica anomalías de concentración de esfuerzo cortante (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."
Función en el sistema: Restricción física fundamental que escala las tasas de recurrencia esperadas en catálogos regionales.
"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."
Función en el sistema: Calcula la energía acústica acumulada y el rebote de Benioff S = ∑10^(0.75 M_i) para detectar quietud sísmica o aceleración de enjambres.
"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."
Función en el sistema: Sustenta las curvas log-normales calibradas de evaluación de colapso de edificios y estimación de víctimas.
"Esteva, L. (1974). Geology and predicting extreme ground motion for earthquake-resistant design. In Proceedings of the 2nd International Conference on Seismic Risk."
Función en el sistema: Cálculo de radios de isosistas para cartografiar la atenuación de la Intensidad Mercalli Modificada (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."
Función en el sistema: Descompone los pronósticos del ensamble en valores de Shapley exactos, garantizando total transparencia.
"Bird, P. (2003). An updated digital model of plate boundaries. Geochemistry, Geophysics, Geosystems, 4(3), 1027."
Función en el sistema: Distancia espacial a límites de placas, tasas de movimiento relativo y cinemática de convergencia.
"Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell, Richmond, VA."
Función en el sistema: Calcula el potencial de marea gravitacional lunar y solar y los ciclos de alineación de sicigias.
"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."
Función en el sistema: 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."
Función en el sistema: 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."
Función en el sistema: International interoperable XML standard for real-time ingestion of NOAA PTWC/NTWC multi-agency hazard bulletins.
Estricta reproducibilidad científica conforme a los estándares sismológicos internacionales
El 100% de las soluciones sismológicas y registros de sismos se ingieren en tiempo real de instituciones públicas abiertas: USGS, EMSC FDSNWS, ISC, BMKG, GeoNet NZ, INGV y GA. Los datos brutos nunca se alteran ni fabrican.
Todos los modelos predictivos de aprendizaje automático se evalúan en bloques de prueba congelados fuera de tiempo (estándar CSEP). El libro de verificación completo con Brier Skill Score (+0.354) y ROC-AUC (0.814) es públicamente auditable mediante la API de métricas.
Rechazamos afirmaciones pseudocientíficas de predicción determinista exacta. La nucleación de fallas obedece a la mecánica de medios continuos no lineal. EarthquakeLive modela estrictamente probabilidades condicionales calibradas P(M ≥ 4.5) para fortalecer la preparación ante emergencias.
Protocolos científicos, operativos y de emergencia verificados
Pase al mapa interactivo WebGL con fallas tectónicas, sismos recientes y la cuadrícula de peligro en tiempo real.