SismoLabMonitoreo + estimación probabilística

Conectando…
Catálogo vacío
Sismicidad global en tiempo real

Panorama

Registro continuo alimentado por el USGS, el EMSC-CSEM y, dentro de Colombia, el SGC. Cada marca del tambor es un sismo real localizado en las últimas 72 h: su posición es la hora exacta en que ocurrió, su color y su altura crecen con la magnitud, de M2,5 a M9,5. La forma del trazo es un símbolo, no la onda que registró un sismómetro.

Tambor sísmico · últimas 72 h —
M < 3
3 – 5
5 – 6
6 – 7
M ≥ 7

Últimos eventos

Orden cronológico inverso. Toca un evento para centrarlo en el mapa.

Centro de alertas

Clasifica cada evento nuevo del feed en vivo por magnitud. El sonido y las notificaciones del sistema se activan sólo cuando tú los pides — ningún navegador los permite sin tu permiso explícito.
Estas alertas avisan de sismos ya ocurridos y ya localizados, con la latencia de las agencias (segundos a minutos). No son alerta temprana: eso exige acelerómetros propios y una red dedicada. Para emergencias en Colombia, la fuente oficial es el SGC y la UNGRD.

Fuentes conectadas

Estado de cada proveedor. Ninguna requiere clave de API.

Contexto ambiental del epicentro

Condiciones atmosféricas modeladas por Open-Meteo en el último evento seleccionado, más la fase lunar calculada para ese instante.

Selecciona un evento del listado.

Lectura correcta: la meteorología y la fase lunar no predicen sismos tectónicos. Este panel es únicamente contexto ambiental y nunca modifica una probabilidad sísmica.
Lugar + fecha → sismo

Buscador de eventos

Escribe un lugar y una fecha en lenguaje natural. El buscador geolocaliza el sitio, consulta el catálogo del USGS alrededor de esa fecha y te devuelve los sismos ordenados por magnitud, con la distancia real a tu punto.

Formatos válidos: «10 agosto 2026», «10/08/2026», «2026-08-10», «agosto 2026», «2026». También acepta un identificador USGS, por ejemplo us6000tjl2.

Resultados

Sin búsqueda todavía.

Secuencia temporal

Magnitud contra tiempo. Deja ver de un vistazo si hubo premonitores y cómo decayeron las réplicas.

Capas del mapa

Ficha del evento

Selecciona un resultado para ver su ficha completa.
Segundos que salvan vidas

Alerta temprana y sismógrafo

Qué sistemas de alerta existen de verdad, cómo activarlos en su teléfono, y un sismógrafo real que usa el acelerómetro del dispositivo para que vea con sus propios ojos cómo funciona la detección.

SismoLab no es un sistema de alerta temprana y no puede serlo. Una página web no puede avisar antes de un sismo: los datos del USGS y del EMSC se publican minutos después, y una pestaña del navegador no corre en segundo plano. La alerta temprana real necesita sensores densos cerca de la falla, procesamiento en menos de un segundo y entrega a nivel de sistema operativo. Lo que sí puede hacer esta página es explicarle cómo activar los sistemas que sí existen y mostrarle cómo detectan.

Alertas de sismos de Android — actívelas hoy

Google convirtió los teléfonos Android en una red de detección: cada equipo usa su acelerómetro como mini-sismómetro. Cuando varios detectan la onda P a la vez, el servidor confirma el sismo y envía la alerta antes de que llegue la onda S, que es la destructiva. En Colombia funciona.
  1. Abra Ajustes y use el buscador interno: escriba «sismos» o «earthquake».
  2. Si no aparece, pruebe Seguridad y emergencia → Alertas de sismos.
  3. Ruta alternativa en algunos fabricantes: Ubicación → Servicios de ubicación → Alertas de sismos.
  4. Active el interruptor y deje la ubicación encendida: sin ella el sistema no puede avisarle.
  5. Revise que el ahorro de batería no esté restringiendo los Servicios de Google Play.
Llegan dos tipos de aviso: uno informativo para sacudida leve, y uno de acción que toma toda la pantalla con sonido fuerte y dice «Agáchate y cúbrete». Sólo se envían para magnitud 4,5 o mayor.

iPhone en Colombia — la situación real

Conviene decirlo sin rodeos para que nadie se confíe.

Las alertas sísmicas nativas de Apple no están habilitadas para Colombia: sólo operan en unos pocos países. Si busca «Alertas sísmicas» en los ajustes de su iPhone y no aparece, no es un fallo del teléfono.

  • Sí conviene activar Ajustes → Notificaciones → Alertas gubernamentales y dejar encendidas las alertas de emergencia: suenan aunque el teléfono esté en silencio.
  • Para avisos sísmicos hace falta una aplicación de terceros de la App Store.
  • Ninguna de esas aplicaciones iguala la cobertura de la red de Android, que se apoya en millones de sensores repartidos por el territorio.
En una casa con teléfonos mixtos, basta con que un Android tenga las alertas activas para que todos oigan el aviso.

Sismógrafo del dispositivo

Este panel lee el acelerómetro real de su teléfono o portátil, igual que hace Android para detectar sismos. Aplica el disparador STA/LTA, el mismo algoritmo que usan las redes sismológicas de verdad. Golpee la mesa suavemente y véalo funcionar.
sensor inactivo
Aceleración del suelo · últimos 20 s —

Detecciones

El STA/LTA compara la energía de una ventana corta (0,5 s) contra una larga (10 s). Si la razón supera el umbral, hay algo que no es ruido de fondo. Así decide un sismómetro —y su teléfono— que empezó a temblar.
Límite honesto: el acelerómetro de un celular tiene un piso de ruido de milésimas de g. Detecta que usted se está moviendo o que la mesa vibra, pero no percibe sismos lejanos ni pequeños. Google lo resuelve cruzando millones de teléfonos a la vez; uno solo no basta. Este panel es demostrativo y educativo: no lo use como alarma.

Sismograma real

Esto ya no es una demostración: descarga la señal miniSEED de una estación sismológica real y la decodifica aquí mismo (Steim-1 y Steim-2 implementados en JavaScript). Verás las ondas del sismo que elijas en la pestaña Buscar.
Selecciona un sismo en la pestaña «Buscar» y vuelve aquí. Si no, se usa el mayor evento del feed en vivo.
Los datos vienen del servicio FDSN de GEOFON (GFZ Potsdam). La cobertura sobre Colombia es escasa —la red nacional del SGC no publica formas de onda abiertas—, así que verás estaciones regionales a cientos o miles de kilómetros. Pedagógicamente eso es una ventaja: a esa distancia las ondas P y S llegan bien separadas y se ven a simple vista.

Instalar SismoLab en el dispositivo

Instalada queda como una aplicación más: icono propio, pantalla completa sin barra del navegador y arranque inmediato. Además guarda una copia para abrirla sin conexión.

Android y escritorio

  1. Toque Instalar arriba, o el aviso que muestra Chrome.
  2. Si no aparece: menú ⋮ → Instalar aplicación (o Añadir a pantalla de inicio).
  3. Requiere HTTPS. Abierta como archivo local no se puede instalar.

iPhone y iPad

  1. Abra la página en Safari (no en Chrome).
  2. Toque Compartir .
  3. Elija Añadir a pantalla de inicio y confirme.

Antes, durante y después

Antes

  • Anclar estantes, televisores y objetos pesados a la pared.
  • Acordar un punto de encuentro familiar y un contacto fuera de la ciudad.
  • Mochila con agua, linterna, radio, botiquín, copia de documentos y algo de efectivo.
  • Saber dónde y cómo cerrar el gas y la electricidad.

Durante

  • Agacharse, cubrirse y agarrarse. No correr hacia las escaleras.
  • Alejarse de ventanas, vidrios y objetos que puedan caer.
  • Si está afuera, ir a un espacio abierto lejos de fachadas y cables.
  • Si va conduciendo, detenerse lejos de puentes y postes.

Después

  • Póngase zapatos: los vidrios rotos causan la mayoría de heridas evitables.
  • Espere réplicas. Son normales y van disminuyendo con los días.
  • No sature las líneas: mensajes de texto en vez de llamadas.
  • Revise fugas de gas antes de encender cualquier llama o interruptor.
Las indicaciones oficiales en Colombia las emiten el Servicio Geológico Colombiano y la UNGRD. Esta aplicación es material educativo: en una emergencia, siga siempre el canal oficial.
Cualquier institución · cualquier municipio · los 33 departamentos

Modo institución

Busca tu institución entre las 22 530 registradas en los datos abiertos del Ministerio de Educación —los 33 departamentos, de La Guajira al Amazonas—, o cualquiera de los 1 122 municipios del país. La aplicación arma un plan de emergencia propio para esa sede: puntos de encuentro, zonas de riesgo, exposición sísmica local y una lista de verificación que funciona sin conexión. Una vez fijada, el catálogo y el pronóstico arrancan allí.

1. Identifica tu institución

Escribe el nombre del colegio, la universidad o el municipio. Los datos vienen del directorio oficial de establecimientos educativos y de la DIVIPOLA del DANE.
Busca por nombre. Si tu sede no aparece, usa el municipio y ajusta el punto en el mapa.

2. Marca los puntos del plan

Elige un tipo y haz clic en el mapa. Todo queda guardado en este equipo: nada viaja a ningún servidor.

3. Exposición sísmica de tu sede

Calculada con el catálogo real alrededor del punto marcado. No es la microzonificación oficial, pero da el orden de magnitud correcto.

4. Lista de verificación post-sismo

Marca lo que se va cumpliendo. El estado se guarda en el equipo y sirve sin conexión.

5. Plan de emergencia imprimible

Una hoja para pegar en cada salón: puntos de encuentro, responsables, teléfonos y qué hacer antes, durante y después.
Este plan es una ayuda organizativa, no un documento oficial. Debe ser revisado y aprobado por directivos y por el comité de gestión del riesgo escolar, y no reemplaza las instrucciones del SGC, la UNGRD, bomberos o defensa civil. En instituciones con menores, no registres ubicaciones personales de estudiantes: marca lugares, no personas.
Inferencia sobre el catálogo cargado

Estadística sísmica

Magnitud de completitud, ley de Gutenberg–Richter por máxima verosimilitud, decaimiento de Omori–Utsu, desagregación de réplicas y liberación de energía. Todo se recalcula sobre los eventos filtrados.

Filtros de análisis

Recorta la muestra antes de estimar. Los cambios se aplican al pulsar «Recalcular».
La comparabilidad de las magnitudes se comprobará al cargar el catálogo.

Ley de Gutenberg–Richter

Distribución acumulada e incremental con el ajuste de máxima verosimilitud (Aki, 1965).

Determinación de Mc

Bondad de ajuste de Wiemer & Wyss (2000): se busca el menor Mc cuyo residual sintético–observado explique ≥ 90 % de los datos.

Evolución temporal de b

Ventana móvil de N eventos con banda de incertidumbre de Shi & Bolt (1982). Caídas sostenidas de b se asocian a aumento del esfuerzo diferencial.

Índice energético y deformación de Benioff

Equivalencia empírica $\log_{10}E = 1{,}5M + 4{,}8$ y su raíz acumulada. Si se mezclan escalas, se muestra sólo como exploración.

Ley de Omori–Utsu

Ajuste por máxima verosimilitud de $n(t)=K/(t+c)^p$ sobre la secuencia posterior al mayor sismo del catálogo.

Tiempos entre eventos

Si el proceso fuese de Poisson puro, los intervalos serían exponenciales con CV = 1. Un CV > 1 delata agrupamiento (réplicas y enjambres).

Distribución de profundidades

La zona de Wadati–Benioff aparece como cola profunda cuando hay subducción.

Cómo se calculó — paso a paso

Trazabilidad completa del cómputo actual, con la convención de colores del aula.
Probabilidades, no predicciones

Pronóstico probabilístico

Modelos didácticos auditables: Poisson estacionario, decaimiento tipo Omori, recurrencia de Paso Browniano y curva de amenaza tipo Cornell. No constituyen un sistema ETAS operativo ni una predicción exacta.

Límite científico, explícito: nadie puede predecir el lugar, la fecha y la magnitud de un sismo futuro. Lo que sí es reproducible es estimar tasas y probabilidades de excedencia en una ventana de tiempo. Todo lo que sigue son tasas condicionadas al catálogo, con sus barras de error.

Región a pronosticar

Escribe cualquier lugar del planeta. La aplicación lo ubica, descarga su catálogo real, estima Mc, b y a, y recalcula todo el pronóstico sobre esos datos.
Sin región analizada. También puedes descargar un catálogo a mano en la pestaña Catálogo.

Parámetros

Se heredan del ajuste de la pestaña Estadística.
7.0
30
$a$ y $b$ no se editan aquí a propósito: salen del ajuste por máxima verosimilitud del catálogo real. Para cambiarlos, cambia el catálogo o los filtros en Estadística — así ningún número queda desconectado de sus datos.

Probabilidad de excedencia vs. horizonte

$P(M\ge m,\;T)=1-e^{-\lambda(m)T}$ con la banda de $\pm 1\sigma_b$ propagada desde el ajuste.

Periodos de retorno

$T_r(m)=1/\lambda(m)$ derivado de la recta de Gutenberg–Richter ajustada.
El «periodo de retorno» es un promedio a la larga, no un turno. Si dice 100 años y ayer hubo uno, mañana puede haber otro: el terreno no lleva calendario.

Secuencia real, día a día

Réplicas acumuladas del sismo mayor reciente, con el catálogo del SGC. Arrastra para mover, rueda para acercar, y arrastra sobre una regla para estirar sólo ese eje.

Réplicas: exploración didáctica tipo Omori

Tasa esperada de réplicas tras el evento principal, con la ley de productividad $K\propto 10^{\alpha(M_m-M_c)}$.
Después de un sismo fuerte, las réplicas son normales y esperables. Son muchas los primeros días y van bajando con el tiempo. Que sigan sintiéndose no significa que venga algo peor.

Recurrencia BPT · sólo para una falla calibrada

Modelo de Paso Browniano (Matthews, Ellsworth & Reasenberg, 2002) para fallas con comportamiento cuasi-periódico: la probabilidad crece con el tiempo transcurrido desde el último evento.
Algunas fallas se comportan como un resorte que se carga poco a poco: mientras más tiempo pasa sin romperse, más probable es que rompa. No todas son así, por eso este modelo se usa sólo en fallas bien estudiadas.

Curva de amenaza (PSHA didáctico)

Integral de Cornell (1968) sobre magnitud y distancia. El modelo de atenuación es editable — se muestra su forma completa para que nada quede oculto.
Lo que le importa a una casa no es la magnitud, es cuánto se sacude el suelo debajo. Esta curva estima con qué frecuencia se supera cada nivel de sacudida. Es el número que usan los ingenieros para diseñar — aquí en versión didáctica.
$\log_{10} \mathrm{PGA}[g] = c_1 + c_2 M - \log_{10} r - c_4 r,\quad r=\sqrt{d^2+h^2}$ con $c_1=-1{,}02$, $c_2=0{,}249$, $c_4=0{,}00255$, $h=7{,}3$ km (forma tipo Joyner–Boore, 1981; roca–suelo genérico, válida aproximadamente para $5{,}0\le M\le 7{,}7$). Verifica los coeficientes contra la fuente original y sustitúyelos por una GMPE regional antes de cualquier uso que no sea docente.

Comparador de regiones

Cada región analizada se puede fijar aquí para contrastarla con las demás bajo exactamente el mismo método. Es la forma más limpia de mostrar por qué el riesgo no es igual en todas partes.
RegiónNAñosMcb ± σ aλ(M)TrP(M, T)
Analiza una región y pulsa «Añadir».

Validación retrospectiva

La pregunta que ningún pronóstico debería esquivar: ¿habría acertado en el pasado? Se ajusta el modelo con los datos anteriores a un año de corte y se comprueba contra lo que realmente ocurrió después. El modelo nunca ve los datos con los que se le evalúa.
Carga un catálogo y elige el año que separa entrenamiento de prueba.
ObjetivoPrevisto GRPersistenciaObservadoError GRIG
Sin ejecutar.

Sin ejecutar.

Fuentes oficiales para Colombia

Lo que sigue son enlaces, no datos incrustados. Abajo explico por qué.
probando SGC · catálogo nacional
Red Sismológica Nacional, vía ArcGIS REST. Se detecta automáticamente. sismo.sgc.gov.co ↗
enlace Amenaza sísmica · SGC
Estudio nacional de amenaza y mapas de zonificación. amenazasismica.sgc.gov.co ↗
enlace UNGRD
Gestión del riesgo de desastres, alertas y reportes de situación. portal.gestiondelriesgo.gov.co ↗
integrada GEOFON · GFZ Potsdam
Catálogo global alemán. Seleccionable como proveedor en la pestaña Catálogo.
integrada INGV
Instituto italiano. Densidad muy alta en el Mediterráneo.
Corrección: una versión anterior de esta aplicación afirmaba que el SGC no publicaba servicios consultables. Era una generalización incorrecta: sólo se había probado la ruta FDSN. El SGC sí expone servicios ArcGIS REST en srvags.sgc.gov.co con el catálogo sísmico nacional, el mapa de amenaza y las estaciones GNSS de GEORED. Esta versión los integra como proveedor con detección automática: si el servicio responde, se usa; si no, se avisa y se sigue con las fuentes internacionales. La UNGRD sí sigue sin catálogo consultable, por eso permanece como enlace.
Catálogos sintéticos

Simulación de Montecarlo

Se generan miles de catálogos artificiales coherentes con los parámetros estimados ($\lambda$, $b$, $M_c$) y se mide la dispersión real de los resultados. Es la forma honesta de poner barras de error a un pronóstico.

Configuración

Las magnitudes se muestrean por inversión de la exponencial truncada $M = M_c - \frac{1}{\beta}\ln\!\left[1-u\left(1-e^{-\beta(M_{max}-M_c)}\right)\right]$, con $\beta = b\ln 10$. El conteo anual se muestrea de una Poisson de media $\lambda$.

Distribución de la magnitud máxima simulada

Histograma de $M_{max}$ en el horizonte, con los percentiles 5, 50 y 95.

Convergencia del estimador

Probabilidad estimada frente al número de simulaciones, con el intervalo de Wald al 95 %. Sirve para decidir cuántas corridas son suficientes.

Resultado e interpretación

Ejecuta la simulación.

Cartografía real y geometría de la sismicidad

Mapas avanzados

Teselas satelitales y de relieve reales, límites de placas de Bird (2003), densidad por kernel, mapa espacial del valor b y corte vertical de la zona de Wadati–Benioff.

M < 3
3 – 5
5 – 6
6 – 7
M ≥ 7

Capas

Modo corte: haz clic en dos puntos del mapa para definir el perfil. El panel inferior mostrará los hipocentros proyectados sobre esa línea.

Corte vertical (Wadati–Benioff)

Distancia a lo largo del perfil vs. profundidad. El buzamiento de la nube revela la geometría de la losa subducida.

Hipocentros en 3D

Longitud, latitud y profundidad. Arrastra para rotar; la rueda hace zoom.
ClimaLab · tiempo, clima y amenazas hidrometeorológicas

Clima

Las mismas herramientas de profenestor.com/clima, aquí dentro y centradas en el municipio que elegiste. Datos de Open-Meteo, IDEAM, GloFAS, CAMS y NOAA.

Abrir aparte

Cargando ClimaLab…

Pyodide · CPython compilado a WebAssembly

Laboratorio Python

NumPy, SciPy, pandas y Matplotlib corriendo dentro del navegador, con el catálogo cargado expuesto como df y como arreglos. Nada sale del equipo.

Editor

Variables disponibles: df (DataFrame), mag, dep, lat, lon, t (días desde el primer evento), más mc_app y b_app con el ajuste vigente de la pestaña Estadística. Usa show(fig) para graficar.

Salida

Entorno no iniciado. La primera carga descarga ~25 MB.
> esperando…
Método, ecuaciones y límites

La ciencia detrás

Cada número de esta aplicación sale de una ecuación publicada. Aquí está cada una, con su origen y su rango de validez.

Lo primero, porque importa: la predicción determinista de sismos —decir «el día D, en el lugar L, de magnitud M»— no existe. No hay ningún precursor validado que funcione de forma reproducible. Lo que la sismología sí hace, y muy bien, es pronóstico probabilístico: cuantificar tasas de ocurrencia y probabilidades de excedencia. Confundir ambas cosas frente a un grupo de estudiantes es el error más caro que se puede cometer en este tema.

Ley de Gutenberg–Richter

Distribución de tamaños. Es una ley de potencias: los sismos no tienen tamaño «típico».

$$\log_{10} N(\ge M) = a - bM$$

  • $N(\ge M)$: número de eventos con magnitud igual o mayor que $M$ en el periodo del catálogo.
  • $a$: productividad sísmica (depende del volumen y del tiempo observado).
  • $b$: pendiente. Cerca de $1{,}0$ en corteza continental. Un $b$ bajo indica predominio de eventos grandes, asociado a mayor esfuerzo diferencial.

La estimación por mínimos cuadrados sobre la curva acumulada está sesgada. Se usa el estimador de máxima verosimilitud de Aki (1965):

$$\hat b = \frac{\log_{10} e}{\bar M - \left(M_c - \tfrac{\Delta M}{2}\right)}$$

con la incertidumbre de Shi & Bolt (1982):

$$\sigma_b = 2{,}30\,\hat b^{\,2}\sqrt{\frac{\sum_{i}(M_i-\bar M)^2}{n(n-1)}}$$

Magnitud de completitud $M_c$

El umbral por encima del cual la red registra todos los sismos. Ignorarlo arruina cualquier estimación de $b$.

Máxima curvatura: $M_c$ es la magnitud del bin más poblado del histograma no acumulado. Rápido, pero subestima.

Bondad de ajuste (Wiemer & Wyss, 2000): para cada candidato $M_c$ se sintetiza la distribución con los $a$, $b$ ajustados y se mide el residual

$$R(M_c) = 100\left(1 - \frac{\sum_i \left|N_i^{obs}-N_i^{sin}\right|}{\sum_i N_i^{obs}}\right)$$

Se adopta el menor $M_c$ con $R\ge 90\,\%$. Si ningún candidato lo alcanza, se relaja a $95\,\%$ del máximo observado y se reporta.

Ley de Omori–Utsu

El decaimiento de las réplicas es el fenómeno más regular de toda la sismología.

$$n(t) = \frac{K}{(t+c)^p}$$

Los parámetros se estiman maximizando la log-verosimilitud de Ogata (1983) para un proceso de Poisson no homogéneo observado en $[t_s, t_e]$:

$$\ln L = N\ln K - p\sum_{i=1}^{N}\ln(t_i+c) - K\,A(p,c)$$

donde $A = \dfrac{(t_e+c)^{1-p}-(t_s+c)^{1-p}}{1-p}$ para $p\neq1$, y $A=\ln\frac{t_e+c}{t_s+c}$ si $p=1$. Típicamente $p\approx1{,}1$ y $c$ es del orden de horas.

Combinada con la productividad $K\propto10^{\alpha(M_m-M_c)}$ y con Gutenberg–Richter para los tamaños, se obtiene el modelo ETAS (Ogata, 1988), que es hoy el estándar operacional para pronóstico de réplicas.

Energía y momento sísmico

$$\log_{10}E\,[\mathrm{J}] = 1{,}5M + 4{,}8 \qquad \log_{10}M_0\,[\mathrm{N\,m}] = 1{,}5M_w + 9{,}1$$

Un punto de magnitud multiplica la energía por $10^{1,5}\approx31{,}6$. Dos puntos, por $\approx1000$.

La deformación de Benioff acumulada, $\sum\sqrt{E_i}$, se usa para visualizar la liberación progresiva de deformación elástica: un tramo recto sugiere régimen estacionario.

Que la energía sea proporcional a $10^{1,5M}$ es la razón física de que un solo sismo M8 libere más energía que todos los M5 de una década en la misma región.

Desagregación de réplicas

Los modelos de Poisson exigen independencia. Las réplicas la violan, así que se retiran.

Ventanas espacio-temporales de Gardner & Knopoff (1974): para un evento de magnitud $M$,

$$L = 10^{\,0{,}1238M+0{,}983}\ \mathrm{km}$$

$$T = \begin{cases}10^{\,0{,}032M+2{,}7389} & M\ge6{,}5\\[2pt] 10^{\,0{,}5409M-0{,}547} & M<6{,}5\end{cases}\ \mathrm{días}$$

Todo evento menor dentro de esa ventana se marca como dependiente. El catálogo restante —de fondo— es el que alimenta el modelo de Poisson.

Probabilidad de excedencia

Con tasa anual $\lambda(m) = 10^{\,a-bm}/\Delta t_{años}$, un proceso de Poisson da

$$P(N\ge1\ \text{en } T) = 1-e^{-\lambda(m)\,T}$$

y el periodo de retorno es $T_r = 1/\lambda$. Cuidado con la lectura: un periodo de retorno de 100 años no significa «uno cada 100 años exactos»; significa que la probabilidad anual es del 1 %, y que en 100 años la probabilidad de al menos uno es $1-e^{-1}=63\,\%$.

La incertidumbre se propaga desde $\sigma_b$ recalculando $\lambda$ con $b\pm\sigma_b$: por eso las curvas llevan banda.

Paso Browniano (BPT)

Para fallas con ciclo sísmico cuasi-periódico, el modelo de Poisson es demasiado pesimista al principio del ciclo y demasiado optimista al final.

$$f(t)=\sqrt{\frac{\mu}{2\pi\alpha^2t^3}}\;\exp\!\left[-\frac{(t-\mu)^2}{2\mu\alpha^2t}\right]$$

La probabilidad condicional en los próximos $\Delta T$ años, habiendo transcurrido $t_e$ desde el último evento, es

$$P = \frac{\int_{t_e}^{t_e+\Delta T} f(t)\,dt}{\int_{t_e}^{\infty} f(t)\,dt}$$

Matthews, Ellsworth & Reasenberg (2002). $\alpha$ es la aperiodicidad; $\alpha\to0$ es un reloj perfecto, $\alpha\approx0{,}5$ es lo observado en fallas reales.

Amenaza sísmica (Cornell, 1968)

La tasa anual con que un parámetro del movimiento fuerte $y$ (aquí PGA) excede un nivel $y^*$ en un sitio es

$$\lambda(y>y^*)=\sum_k \nu_k \iint P\!\left(Y>y^*\mid m,r\right) f_M(m)\,f_R(r)\,dm\,dr$$

La aplicación discretiza el catálogo como fuente puntual distribuida: cada celda aporta su tasa propia, $f_M$ es la GR truncada y la atenuación es la mostrada en el panel de pronóstico.

Un PSHA real usa modelos de falla con geometría, tasas de deslizamiento geodésicas, varias GMPE ponderadas en árbol lógico y efectos de sitio. Este módulo es una maqueta didáctica del razonamiento, no un insumo de diseño estructural. Para eso está la NSR-10 y el estudio de amenaza del Servicio Geológico Colombiano.

Referencias

Aki, K. (1965). Maximum likelihood estimate of b in the formula log N = a − bM. Bull. Earthq. Res. Inst., 43, 237–239.

Bird, P. (2003). An updated digital model of plate boundaries. Geochem. Geophys. Geosyst., 4(3).

Cornell, C. A. (1968). Engineering seismic risk analysis. BSSA, 58(5), 1583–1606.

Gardner, J. K. & Knopoff, L. (1974). Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian? BSSA, 64(5).

Gutenberg, B. & Richter, C. F. (1944). Frequency of earthquakes in California. BSSA, 34(4), 185–188.

Joyner, W. B. & Boore, D. M. (1981). Peak horizontal acceleration and velocity from strong-motion records. BSSA, 71(6).

Kanamori, H. (1977). The energy release in great earthquakes. JGR, 82(20), 2981–2987.

Matthews, M. V., Ellsworth, W. L. & Reasenberg, P. A. (2002). A Brownian model for recurrent earthquakes. BSSA, 92(6).

Ogata, Y. (1988). Statistical models for earthquake occurrences and residual analysis for point processes. JASA, 83(401).

Shi, Y. & Bolt, B. A. (1982). The standard error of the magnitude–frequency b value. BSSA, 72(5).

Utsu, T., Ogata, Y. & Matsu'ura, R. (1995). The centenary of the Omori formula for a decay law of aftershock activity. J. Phys. Earth, 43.

Wiemer, S. & Wyss, M. (2000). Minimum magnitude of completeness in earthquake catalogs. BSSA, 90(4).

Datos: USGS Earthquake Hazards Program — servicio FDSNWS earthquake.usgs.gov · EMSC-CSEM seismicportal.eu · límites de placas PB2002 (Bird, 2003) · meteorología Open-Meteo · teselas Esri World Imagery / Shaded Relief / Ocean y OpenStreetMap. Todos de acceso abierto.
Informe de análisis
Antes de usarlo en clase

Qué es y qué no es SismoLab

No predice sismos. Ningún método validado lo hace. SismoLab estima tasas y probabilidades a partir de catálogos reales, con las ecuaciones a la vista y las barras de error puestas.

Sí hace

No hace

Cuando el catálogo cambia, todos los números cambian. Esa dependencia es el mensaje pedagógico central: en sismología estadística el resultado nunca es un número, es un número con su intervalo y su modelo declarado.