Topology-Phase Monitoring of Kilauea Eruption Preparation: A Reproducible Seismic Audit and Forecasting Potential Assessment

Resumen Técnico: Monitoreo de la Fase de Topología de la Preparación de la Erupción del Kilauea

Planteamiento del Problema
El pronóstico de erupciones volcánicas de alta precisión sigue siendo un desafío significativo porque la preparación de una erupción rara vez es señalada por una única variable escalar (por ejemplo, un pico en el conteo de terremotos). En su lugar, implica una reorganización distribuida de la sismicidad, la deformación, el transporte de gases y la conectividad magma-fractura. El monitoreo tradicional suele depender de canales aislados, lo cual puede ser ambiguo; por ejemplo, los enjambres de terremotos no siempre conducen a erupciones, y un aumento en el número de eventos puede reflejar un ajuste de estrés en lugar de actividad inminente. El artículo postula que la preparación de una erupción puede entenderse mejor como una transición en el "estado-topológico" del campo sísmico: pasando de una actividad de fondo difusa a un estado conectado, migratorio y anisotrópico.

Metodología
El estudio realiza una auditoría reproducible, basada únicamente en sismicidad, del sistema volcánico Kilauea utilizando datos públicos del catálogo de terremotos FDSN de la USGS.

• Alcance de los Datos: El análisis cubre un radio de 80 km alrededor del Kilauea (19.4069°N, 155.2834°W) para terremotos con magnitud $M \ge 1.0$.
• Ventanas Temporales: Los datos se agregaron para ventanas de mismo calendario (del 1 de enero al 5 de mayo) de 2014 a 2019. Las ventanas de 2014–2017 sirvieron como la distribución de referencia de entrenamiento, mientras que la ventana de 2018 fue el caso de prueba objetivo para la erupción de la zona de rift este inferior (LERZ) (inicio: 3 de mayo de 2018).
• Vector de Estado-Topológico: Una ventana rodante de siete días se convirtió en un vector de siete características:
  1. Conteo de eventos ($N$).
  2. Proxy de energía de magnitud ($E$).
  3. Fracción de la componente conectada más grande ($C_{max}$).
  4. Densidad de aristas del grafo ($\rho_G$).
  5. Entropía espacial normalizada ($S_{sp}$).
  6. Anisotropía de la nube de eventos ($A$).
  7. Migración del centroide ($\Delta r_c$).
  • Construcción del Grafo: Los eventos se conectaron si estaban separados por $\le 8$ km en el espacio proyectado y $\le 2$ días en el tiempo.
• Construcción del Indicador: Se calculó un Indicador de Fase de Topología (TPS) compuesto como la media de los componentes estandarizados. La entropía espacial se negó para reflejar la transición de una organización difusa a una concentrada. Los umbrales se fijaron en el percentil 99 de la distribución de entrenamiento de 2014–2017.

Resultados Clave
La auditoría, utilizando parámetros bloqueados entrenados con los datos de 2014–2017, arrojó los siguientes hallazgos respecto a la erupción de la LERZ en 2018:

• Alerta Temprana: El indicador compuesto de fase de topología cruzó por primera vez el umbral del percentil 99 de entrenamiento el 12 de abril de 2018, proporcionando un tiempo de adelanto de 21 días antes del inicio de la erupción del 3 de mayo.
• Comparación con Métricas Escalares:
  • Conteo de Eventos: Cruzó su percentil 99 solo el 1 de mayo de 2018 (2 días antes del inicio).
  • Proxy de Energía: No cruzó su umbral antes del inicio.
• Análisis de Componentes: La alerta temprana fue impulsada por una reorganización estructural coordinada. La densidad del grafo superó el umbral el 12 de abril, mientras que la fracción de la componente conectada más grande y la entropía espacial superaron sus respectivos umbrales el 17 de abril. La migración mostró un cruce temprano en marzo, pero no excedió el máximo de entrenamiento hasta el 1 de mayo. La anisotropía superó el umbral temprano en enero, pero se consideró insuficiente como predictor independiente.
• Valores Pico: El indicador de topología alcanzó un valor máximo de 4.0411 el 1 de mayo de 2018, excediendo tanto el percentil 99 de entrenamiento (2.0462) como el máximo de entrenamiento (2.8166).

Contribuciones Clave

1. Detección del Estado-Topológico: El estudio demuestra que la preparación de una erupción en un sistema basáltico bien monitoreado puede detectarse como una transición topológica en el campo sísmico, distinta de un simple aumento en el conteo de eventos o energía.
2. Auditoría Reproducible: Los autores proporcionan un flujo de trabajo totalmente reproducible, incluyendo scripts específicos y archivos de datos, lo que permite la verificación independiente de la auditoría de "parámetros bloqueados".
3. Mejora del Tiempo de Adelanto: Los resultados sugieren que la combinación de métricas de conectividad, densidad y organización espacial puede extender significativamente el tiempo de adelanto para la detección de estados de preparación de erupciones en comparación con las líneas base escalares tradicionales.

Significado y Reivindicaciones
El artículo realiza una reivindicación modesta y específica: no afirma que la física del estado-topológico pueda predecir de manera determinista la fecha exacta, la ubicación de la chimenea o la magnitud de todas las erupciones volcánicas.

• Preaviso Probabilístico: La principal contribución es la validación de una detección de múltiples canales y probabilística de un "estado de preparación de erupción". El método identifica cuándo un sistema ha entrado en un estado de topología inusualmente organizado en relación con las líneas de base históricas.
• Objetivo Operativo: El objetivo a corto plazo apropiado no es la predicción determinística, sino la mejora del tiempo de adelanto y la reducción de las falsas alarmas escalares mediante el monitoreo del estado-topológico.
• Limitaciones y Pasos Futuros: Los autores reconocen que la auditoría actual se basa únicamente en sismicidad y depende de un único caso fuerte (Kilauea 2018). Establecen que, para que el método avance hacia una reivindicación operativa de alta precisión, requiere:
  • Validación en intervalos de inquietud no eruptivos independientes para estimar las tasas de falsas alarmas.
  • Integración de datos multicanal (deformación, gases, térmicos, inclinación).
  • Comparación contra modelos estándar (por ejemplo, valor-b, ETAS).
  • Bloqueo estricto de parámetros antes de probar en nuevos eventos.

En conclusión, el artículo respalda la hipótesis de que la preparación de la erupción del Kilauea en 2018 se manifestó como una reorganización topológica detectable de la sismicidad, ofreciendo una vía potencial para la detección de estados de forma probabilística y temprana.