Seismic Compression-Quiescence Anomalies as Earthquake Precursors: Prospective Prediction for the Kashmir Himalayan Fault Belt and Validation Across Four Tectonic Regimes

1. Planteamiento del problema

El arco de Cachemira-Himalaya representa una de las regiones sísmicamente más peligrosas del mundo, impulsada por la convergencia hacia el norte de las placas India y Euroasiática. Esta convergencia almacena una energía de deformación elástica significativa a lo largo de sistemas de fallas paralelos (Falla Frontal Principal, Falla Límite Principal, Falla de Panjal y Falla Central Principal). A pesar de catástrofes históricas como el terremoto de Cachemira de M 7.6 en 2005, los métodos actuales para la predicción de terremotos en esta región enfrentan tres limitaciones específicas:

• Falta de un umbral cuantitativo único que sea válido en diferentes entornos tectónicos.
• Ausencia de registros de predicción prospectivamente verificables específicamente para Cachemira.
• Incapacidad para atribuir las señales de carga a filamentos de falla específicos en lugar de regiones geográficas amplias.

Los métodos existentes, como la quiescencia sísmica (estadístico Z), el mapeo del valor-b y la liberación acelerada de momento, no han logrado alcanzar estos objetivos simultáneamente.

2. Metodología

El estudio introduce un marco basado en datos centrado en el Índice de Compresión de Magnitud ($C$) y la señal de caída de C (C-drop), validado a través de cuatro regímenes tectónicos distintos (Cachemira, Japón, Turquía, Italia).

Métricas principales:

• Índice de Compresión de Magnitud ($C$): Una medida libre de parámetros definida como $C(k) = \max(0, 1 - \sigma[W] / \sigma_0)$. Cuantifica el estrechamiento de las distribuciones de magnitud de los terremotos a medida que una falla acumula tensión. A medida que la carga de tensión homogeneiza la sismicidad local, la desviación estándar ($\sigma$) de las magnitudes disminuye, causando que $C$ aumente hacia 1.
• La caída de C (C-drop): La señal precursora primaria. Se define como un descenso sistemático y sostenido post-pico en $C$ (al menos 0.050 desde un pico $\ge$ 0.25). Esta caída ocurre cuando la falla comienza a fallar, permitiendo que terremotos de mayor magnitud entren en la ventana de observación y aumentando la varianza de la magnitud.
• Índice de Quiescencia Sísmica ($Q$): Una medida secundaria de reducción de tasa ($Q = \max(0, 1 - \lambda[W] / \lambda_0$). Una $Q$ elevada combinada con una $C$ elevada indica una falla críticamente bloqueada.

Procesamiento de datos y identificación de zonas:

• Ventana rodante: Los cálculos utilizan una ventana rodante de los $W=20$ terremotos más recientes en una zona, suavizada por un promedio móvil de 12 eventos.
• Línea de base: Los valores de referencia ($\sigma_0, \lambda_0$) se derivan del periodo 2000–2009.
• Zonificación basada en datos: A diferencia de estudios previos que utilizan coordenadas arbitrarias, los centros de las zonas se ubican en picos naturales de densidad sísmica identificados mediante suavizado gaussiano del catálogo ISC (2000–2026). Los radios de las zonas capturan el 90% de los cúmulos de sismicidad local.
• Sistema de puntuación de seis criterios: Para filtrar falsos positivos (por ejemplo, secuencias de réplicas), una zona se clasifica como CRÍTICA solo si cumple al menos 3 de 6 criterios simultáneamente:
  1. Compresión: $C \ge 0.25$.
  2. Quiescencia: $Q \ge 0.20$.
  3. Caída de C: Descenso confirmado desde el pico.
  4. Valor-b: $b < 0.75$ (indicando concentración de tensión).
  5. Baja actividad: $<15$ eventos en los últimos 90 días.
  6. Frescura: El pico ocurrió dentro de los últimos 180 días.

3. Contribuciones clave

• La señal de caída de C (C-drop): Establece el descenso post-pico en la compresión de magnitud como un precursor primario, físicamente fundamentado, para terremotos significativos inminentes, distinto de la simple quiescencia basada en la tasa.
• Descomposición de zonas basada en datos: Demuestra que las zonas sísmicas compuestas pueden descomponerse para identificar filamentos de falla activos específicos. Por ejemplo, el estudio aisla la señal de carga activa en la región de Anantnag a la Falla de Panjal (Ramban–Banihal) en lugar del núcleo descargado de la Falla Central Principal.
• Verificación prospectiva: El artículo presenta una predicción formal archivada en Zenodo (18 de junio de 2026) antes de la ocurrencia de eventos activadores, proporcionando un registro verificable de la aplicación del método.
• Validación trans-tectónica: El método se prueba retrospectivamente en cuatro tipos de límites de placa (colisión, subducción, desgarre y arco de compresión), demostrando la universalidad del mecanismo de la caída de C.

4. Resultados

Validación retrospectiva:
A través de 53 ciclos confirmados de caída de C en cuatro regímenes tectónicos, un terremoto de magnitud objetivo ocurrió dentro de los 90 días en el 83.0% de los casos (44/53).

• Japón (Subducción): Tasa de detección del 88.9%.
• Turquía (Transformante): Tasa de detección del 91.7%.
• Italia (Compresión): Tasa de detección del 75.0%.
• India/Cachemira (Colisión): Tasa de detección del 73.3%.
• Subconjunto de alta confianza: Para ciclos con $C_{pico} \ge 0.30$ y caída $\ge 0.055$, la tasa de detección alcanzó el 100% en todas las regiones.

Estado de la zona de Cachemira (al 20 de junio de 2026):

• Zona Z21 (Anantnag): Puntuación 5/6. $C_{drop} = 0.129$. Dos eventos de M 4.8 ocurrieron el 11 de abril y el 5 de junio de 2026, durante la fase activa de la caída de C. El ciclo permanece abierto con $C=0.260$ (67% de la carga máxima retenida).
• Zona Z15b (Norte de Cachemira): Puntuación 5/6. $C_{drop} = 0.092$. Mayor carga en 26 años ($C_{pico}=0.585$). La tasa de detección histórica para M $\ge$ 4.5 es del 96.3%. Predicción de M 5.0–6.0 para julio–septiembre de 2026.
• Zona Z15 (Muzaffarabad): Puntuación 4/6. $C=0.342$ con una quiescencia excepcionalmente profunda ($Q=0.795$). Se espera la confirmación de la caída de C (umbral 0.050) dentro de 2–4 semanas.
• Zona FZ-NPB (Gilgit): Puntuación 4/6. Valor de $C$ más alto del estudio (0.530). La caída de C acaba de comenzar; se ha marcado una predicción emergente para agosto–octubre de 2026.

Eventos de validación:
Tres eventos validaron la ventana de predicción prospectiva:

1. M 3.3 (19 de junio de 2026): Ocurrió 18 horas después del archivo de Zenodo, rompiendo un silencio de 13 días en la Zona Z21.
2. M 4.8 (11 de abril de 2026) y M 4.8 (5 de junio de 2026): Ambos ocurrieron durante la fase activa de la caída de C de la Zona Z21, representando liberaciones parciales de tensión.

5. Significado y Reivindicaciones

El artículo afirma proporcionar una predicción de terremotos prospectiva y prospectivamente verificable para el cinturón de fallas de Cachemira-Himalaya, abordando una brecha crítica en la práctica sismológica actual.

• Mecanismo físico: Los autores argumentan que la caída de C refleja un proceso universal de mecánica de fallas donde la varianza de la magnitud disminuye bajo carga y aumenta tras la falla, independientemente del entorno tectónico.
• Precisión: Al utilizar la descomposición de zonas basada en datos, el estudio refina la atribución espacial de la tensión, identificando la Falla de Panjal como la fuente activa para la señal de Anantnag en lugar de la anteriormente asumida Falla Central Principal.
• Evaluación de riesgos: El estudio destaca que aproximadamente 10 millones de personas viven en las zonas afectadas a través de Jammu y Cachemira, Cachemira administrada por Pakistán y Gilgit-Baltistán.
• Llamado a la acción: Los autores invitan a la comunidad sismológica internacional (incluyendo NCS y CSEP) a monitorear y evaluar independientemente la predicción, enfatizando que el método es autoconsistente para distinguir entre estados de falla cargados (predictibles) y descargados (no predictibles).

El artículo concluye que el método de la caída de C, combinado con la zonificación basada en datos y un filtro de seis criterios, ofrece un marco robusto para identificar estados críticos de carga de fallas y emitir pronósticos prospectivos con alta confianza estadística.