Historia real
Se usa el catálogo público de USGS (ComCat), con 74.384 sismos en y cerca de Chile desde 1964. Incluye el de Maule de 2010.
Un proyecto de investigación · Chile
Este mapa compara unas zonas con otras usando el historial real de terremotos. Es investigación, no un aviso de emergencia.
En cuatro pasos
Sin magia y sin inteligencia artificial. Solo historia, un umbral y un modelo clásico.
Se usa el catálogo público de USGS (ComCat), con 74.384 sismos en y cerca de Chile desde 1964. Incluye el de Maule de 2010.
El modelo solo mira sismos de magnitud 5 o más (tipo mwc). Debajo de ese umbral el catálogo global se ve incompleto.
Cada zona tiene un ritmo de fondo. Después de un sismo hay más réplicas cerca, y luego se calma. Eso se junta en un modelo ETAS.
Chile se parte en celdas. En cada una se estima la chance de al menos un sismo, por rango de magnitud, en los próximos 7 días.
Los próximos 7 días
El color más cálido marca más actividad esperada, comparada con el resto de Chile. No es un porcentaje de peligro.
Lo que ya ocurrió
Los puntos del gráfico son observaciones de USGS, no del modelo.
Cargando sismos recientes de USGS…
Cada cálculo queda guardado. El de arriba es el que se usa. El de magnitud mb se conserva como archivo y no se borra.
Método
CHILE-OEF es un experimento de pronóstico operacional de sismos (Operational Earthquake Forecasting, Jordan 2011). No busca adivinar un terremoto. Busca estimar, con un modelo publicado, dónde hay más actividad esperada en una ventana corta.
El catálogo vivo es USGS ComCat (FDSN), no la red del Centro Sismológico
Nacional. Por eso los sismos chicos se ven peor que en un catálogo local.
La fila científica con magnitud mb se conservó: ese tipo se satura
cerca de 6,2 y da un valor b artificialmente alto (~2,13). El mapa público
usa magnitud momento mwc, con b ≈ 1,12.
Primero se estima la magnitud de completez (Mc) con el método EMR de Ogata y Katsura. Luego la pendiente Gutenberg–Richter con el estimador de Aki y la corrección de binning de Utsu. El motor espacial es un ETAS espacio-temporal (Ogata 1988/1998) con fondo homogéneo y γ = 0. Cada celda de 0,1° y cada rango de magnitud reciben un recuento esperado; la cifra del mapa es la probabilidad de Poisson de al menos un evento, P(N ≥ 1) = 1 − exp(−λ).
Las ciudades suman las λ de las celdas a 40 km y aplican la misma fórmula. Eso no es amenaza sísmica, ni Intensidad de Mercalli, ni un aviso para tu casa.
No hay puntaje CSEP publicado. El catálogo se cargó de una vez, así que
available_at no reproduce la latencia histórica. Una prueba hacia
el pasado vería un catálogo vacío o “hacer trampa”. La prueba justa
empieza desde ahora hacia adelante. Tampoco hay intervalos de
incertidumbre de los parámetros, ni aprendizaje automático, ni datos
del CSN en el modelo.
El USGS publica outlooks de réplicas (OAF). Japón tiene outlooks de réplicas de la JMA. Italia (INGV) opera OEF semanal para autoridades, no como mapa público de este tipo. CHILE-OEF no es un equivalente operacional de ninguno de ellos.
Con honestidad
Aún no comparamos este mapa con sismos que ocurrieron después de publicarlo. No vamos a inventar un puntaje.
El modelo no cambia el ritmo de fondo según la fosa, la placa o la ciudad. γ = 0: la extensión espacial no crece con la magnitud.
USGS ve peor los sismos chicos. Por eso el mapa parte en magnitud 5. Un catálogo chileno local sería otra cosa.
El porcentaje urbano suma celdas a 40 km. No es riesgo de tu edificio, ni de tsunami, ni de que se caiga algo.
No hay intervalo de los parámetros. El estado del pronóstico es uncalibrated_point_forecast: un punto, sin calibrar.
No hay semáforo, no hay “salgan de la casa”, no hay fecha. Si hay emergencia: CSN y SENAPRED.
Para saber si el mapa sirve hay que dejarlo funcionar y, después, comparar con lo que realmente ocurrió. No sirve “probarlo” con terremotos que ya conocemos, porque eso hace ver cualquier modelo mejor de lo que es.
Recién cargamos toda la historia. Así que la prueba justa empieza desde ahora hacia adelante. Cuando haya resultados reales, los vamos a mostrar aquí.
El motor de evaluación existe (pruebas CSEP N/M/S/L en caminata sintética). Cada sismo del backfill quedó disponible al momento de la carga, no en la fecha histórica del evento. Un cálculo hacia el pasado vería un catálogo vacío o usaría información que entonces no existía. Por eso no hay un backtest publicado ni un puntaje CSEP en esta página.
Preguntas cortas
No. Nadie puede dar fecha, hora, lugar y tamaño con certeza. Esto solo compara actividad esperada.
Nada distinto de lo de siempre: kit, punto de encuentro, y seguir al CSN y a SENAPRED si hay un sismo fuerte.
Porque un sismo de magnitud 5 a 6, en 7 días y a 40 km, es raro incluso en Chile. Un 0,14% es “1 en 700”, no “va a temblar”.
No. El mapa público es ETAS + Gutenberg–Richter. Si algún día se prueba otra cosa, se va a decir con claridad, y solo si gana en una prueba hacia adelante.