Kuyuy · el método

Cómo calculamos cada cifra

Kuyuy reporta lo que el IGP y el USGS ya registraron, con contexto estadístico. Esta página explica de dónde sale cada dato, qué es una medición oficial y qué es un estimado nuestro, con las fórmulas y sus referencias. La regla de la casa: cada cifra enlaza a una fuente que la respalda, y lo estimado se dice estimado.

El catálogo

Combinamos dos fuentes oficiales: el catálogo del IGP (CENSIS), que registra sismos de magnitud ~3 o mayor desde 2012, y el del USGS, con eventos de magnitud 4.5 o mayor desde 1901. Los eventos duplicados entre fuentes se reconcilian (90 segundos y 100 km de tolerancia) y cada ficha muestra la magnitud que reportó cada fuente cuando difieren. El catálogo se actualiza cada 2 minutos.

El catálogo es heterogéneo y lo decimos siempre: antes de 2012 solo hay eventos M≥4.5 (USGS); desde 2012 entra el IGP con M≥~3. Por eso toda comparación entre épocas (histogramas por década, tasas históricas) usa solo M≥4.5, el único umbral común. El salto de densidad que se ve en los gráficos alrededor de 2012 es del registro, no de la sismicidad.

Lo que es oficial y lo que es estimado

Oficial (medido o publicado por una institución): magnitud, profundidad, epicentro, intensidad reportada (escala del IGP), mapas oficiales del evento y de aceleraciones (PGA) del IGP, estaciones de la Red Acelerométrica Nacional, zonificación sísmica de suelos (IGP), bordes de placas (modelo PB2002 de Bird 2003), fallas activas (GEM Global Active Faults) y geometría de la placa subducida (Slab2 del USGS).

Estimado (calculado por Kuyuy con fórmulas publicadas): todo lo de la tarjeta “La física del evento”, el radio potencial de percepción, los segundos de aviso que haría posibles el SASPe, las probabilidades del modo Modelo y los percentiles de zona. Cada uno se explica abajo.

La física del evento (ficha)

Son estimados de orden de magnitud derivados solo de la magnitud, con relaciones empíricas estándar. Cuando la magnitud del catálogo no es Mw (momento), la tratamos como Mw: es un proxy razonable en estos rangos y lo declaramos en la propia tarjeta.

• Energía radiada: log₁₀E = 1.5M + 4.8 (E en julios), Gutenberg-Richter 1956. La equivalencia en TNT usa 1 t = 4.184×10⁹ J. Cada punto de magnitud son ~32 veces más energía.
• Momento sísmico: M₀ = 10^(1.5M + 9.1) N·m (Hanks y Kanamori 1979).
• Tamaño de la ruptura: para eventos de interfase de subducción usamos Strasser et al. 2010; para el resto, Wells y Coppersmith 1994. La incertidumbre típica es de un factor ~2, y el círculo del mapa es ilustrativo: la orientación real de la ruptura no se conoce sin más datos.
• Deslizamiento promedio: D = M₀/(μ·A) (Aki 1966), con rigidez μ de 3×10¹⁰ Pa (cortical) o 4×10¹⁰ Pa (interfase).
• Duraciones: la de ruptura escala con M₀^(1/3) (el doble de la half-duration del catálogo GCMT); la de sacudida fuerte usa Dobry et al. 1978, definida en roca: en suelos blandos dura más, por eso enlazamos la capa de suelos del IGP.
• ¿Se suele sentir?: frecuencia histórica de reportes “sentido” del IGP (2012-hoy) por banda de magnitud y profundidad. Mide reportes, no física pura: también refleja dónde vive la gente.

Percentiles y análisis de zona

El análisis de zona cuenta lo que hay en el catálogo dentro del área elegida: totales, máximos, distribución por década y profundidad, y tasas (media histórica M4.5+ vs los últimos 90 días). El “top X%” de una ficha es el ranking del evento entre todos los registrados en 100 km a la redonda. Nada de esto pronostica: describe el pasado registrado, con la advertencia de heterogeneidad del catálogo.

El modo Modelo (Gutenberg-Richter + Poisson)

La ley de Gutenberg-Richter describe cuántos sismos chicos hay por cada grande (en el Perú, b ≈ 0.84 sobre M≥4.5 desde 1950). Con la tasa local de cada celda de ~111 km y un proceso de Poisson calculamos la probabilidad de al menos un sismo de cierta magnitud en una ventana de años. Es estadística de recurrencia, no predicción: el modelo no tiene memoria, no sabe de fechas y una probabilidad no “crece” porque haya pasado tiempo sin sismos.

Réplicas y secuencias

La sección “La secuencia” agrupa el catálogo con las ventanas espacio-temporales clásicas de Gardner y Knopoff (1974): tras un evento de magnitud M, los sismos menores dentro de cierta distancia y cierto plazo se cuentan como sus réplicas. Es una agrupación estadística del pasado, no un pronóstico. Como referencia empírica citamos la ley de Båth: la réplica mayor de una secuencia suele rondar 1.2 grados menos que el principal. Las secuencias anteriores a 2012 se ven incompletas porque el catálogo solo registraba M≥4.5.

El SASPe y los segundos de aviso

La etiqueta “califica para alerta SASPe” aplica la regla pública del sistema (M≥6.0 con epicentro marino). Los segundos de aviso son un estimado referencial (distancia entre ~6 km/s menos la detección), no un dato oficial del evento, y que un evento califique no significa que las sirenas sonaran: la activación real es selectiva y usa la estimación en tiempo real del IGP.

Los productos del USGS (PAGER, ShakeMap, DYFI, mecanismo focal)

Para eventos que el USGS también registra, la ficha muestra sus productos automáticos: el semáforo PAGER (estimación de impacto y población expuesta), los contornos de intensidad ShakeMap (modelo de sacudida, no medición en cada punto), los reportes ciudadanos DYFI (“Did You Feel It?”) y el mecanismo focal (la “beachball”, que resume cómo se movió la falla). Son estimaciones automáticas del USGS, no cifras oficiales peruanas, y cada módulo enlaza a la página oficial del evento.

Lo que no hacemos

No predecimos sismos. Nadie puede decir la fecha, el lugar o la magnitud de un próximo sismo, y desconfía de quien lo afirme. Lo único que la sismología moderna pronostica con honestidad son probabilidades de réplicas tras un evento grande, siempre como rangos con ventana y vencimiento; cuando Kuyuy incorpore ese módulo, se comunicará exactamente así. Mientras tanto, lo que sí está en tus manos: mochila de emergencia, plan familiar y zonas seguras (INDECI).

Referencias

• Gutenberg, B. y Richter, C.F. (1956). Earthquake magnitude, intensity, energy, and acceleration.
• Hanks, T.C. y Kanamori, H. (1979). A moment magnitude scale. JGR 84.
• Wells, D.L. y Coppersmith, K.J. (1994). New empirical relationships… BSSA 84(4).
• Strasser, F.O., Arango, M.C. y Bommer, J.J. (2010). Scaling of the source dimensions of interface and intraslab subduction-zone earthquakes. SRL 81(6).
• Aki, K. (1966). Generation and propagation of G waves from the Niigata earthquake.
• Dobry, R., Idriss, I.M. y Ng, E. (1978). Duration characteristics of horizontal components of strong-motion earthquake records. BSSA 68(5).
• Bird, P. (2003). An updated digital model of plate boundaries. G³ 4(3).
• Hayes, G.P. et al. (2018). Slab2, a comprehensive subduction zone geometry model. Science 362.
• IGP · SASPe · GEM Global Active Faults

¿Encontraste un error o una cifra sin fuente? Es un reporte para nosotros: escríbenos.