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

Resumo Técnico: Monitoramento da Fase de Topologia da Preparação da Erupção do Kilauea

Definição do Problema
A previsão de erupções vulcânicas de alta precisão continua sendo um desafio significativo porque a preparação para a erupção raramente é sinalizada por uma única variável escalar (por exemplo, um pico na contagem de terremotos). Em vez disso, envolve uma reorganização distribuída da sismicidade, deformação, transporte de gases e conectividade magma-fratura. O monitoramento tradicional muitas vezes depende de canais isolados, o que pode ser ambíguo; por exemplo, enxames de terremotos nem sempre levam a erupções, e o aumento na contagem de eventos pode refletir o ajuste de tensão em vez de atividade iminente. O artigo postula que a preparação da erupção pode ser melhor compreendida como uma transição no "estado-topológico" do campo sísmico — mudando de uma atividade de fundo difusa para um estado conectado, migratório e anisotrópico.

Metodologia
O estudo realiza uma auditoria reprodutível, baseada apenas em sismicidade, do sistema vulcânico Kilauea utilizando dados públicos do catálogo de terremotos FDSN da USGS.

• Escopo de Dados: A análise cobre um raio de 80 km ao redor do Kilauea (19,4069°N, 155,2834°W) para terremotos com magnitude $M \ge 1,0$.
• Janelas Temporais: Os dados foram agregados para janelas de mesmo calendário (1º de janeiro a 5 de maio) de 2014 a 2019. As janelas de 2014–2017 serviram como a distribuição de referência de treinamento, enquanto a janela de 2018 foi o caso de teste alvo para a erupção da zona de rift inferior do leste (LERZ) (início: 3 de maio de 2018).
• Vetor de Estado-Topológico: Uma janela móvel de sete dias foi convertida em um vetor de sete características:
  1. Contagem de eventos ($N$).
  2. Proxy de energia de magnitude ($E$).
  3. Fração da maior componente conectada ($C_{max}$).
  4. Densidade de arestas do grafo ($\rho_G$).
  5. Entropia espacial normalizada ($S_{sp}$).
  6. Anisotropia da nuvem de eventos ($A$).
  7. Migração do centroide ($\Delta r_c$).
  • Construção do Grafo: Os eventos foram conectados se separados por $\le 8$ km em espaço projetado e $\le 2$ dias em tempo.
• Construção do Indicador: Um indicador composto de Fase-Topológica (TPS) foi calculado como a média dos componentes padronizados. A entropia espacial foi negada para refletir a transição de uma organização difusa para uma concentrada. Os limiares foram definidos no percentil 99 da distribuição de treinamento de 2014–2017.

Principais Resultados
A auditoria, utilizando parâmetros travados treinados em 2014–2017, gerou os seguintes resultados em relação à erupção da LERZ em 2018:

• Alerta Antecipado: O indicador composto de fase-topológica cruzou pela primeira vez o limiar do percentil 99 do treinamento em 12 de abril de 2018, proporcionando um tempo de antecedência de 21 dias antes do início da erupção em 3 de maio.
• Comparação com Métricas Escalares:
  • Contagem de Eventos: Cruzou seu percentil 99 apenas em 1 de maio de 2018 (2 dias antes do início).
  • Proxy de Energia: Não cruzou seu limiar antes do início.
• Análise de Componentes: O alerta antecipado foi impulsionado por uma reorganização estrutural coordenada. A densidade do grafo ultrapassou o limiar em 12 de abril, enquanto a fração da maior componente conectada e a entropia espacial ultrapassaram seus respectivos limiares em 17 de abril. A migração mostrou um cruzamento precoce em março, mas não excedeu o máximo do treinamento até 1 de maio. A anisotropia ultrapassou precocemente em janeiro, mas foi considerada insuficiente como um preditor isolado.
• Valores de Pico: O indicador de topologia atingiu um valor de pico de 4,0411 em 1 de maio de 2018, excedendo tanto o percentil 99 do treinamento (2,0462) quanto o máximo do treinamento (2,8166).

Principais Contribuições

1. Detecção de Estado-Topológico: O estudo demonstra que a preparação de uma erupção em um sistema basáltico bem monitorado pode ser detectada como uma transição topológica no campo sísmico, distinta de simples aumentos na contagem de eventos ou energia.
2. Auditoria Reprodutível: Os autores fornecem um fluxo de trabalho totalmente reprodutível, incluindo scripts específicos e arquivos de dados, permitindo a verificação independente da auditoria de "parâmetros travados".
3. Melhoria no Tempo de Antecedência: Os resultados sugerem que a combinação de métricas de conectividade, densidade e organização espacial pode estender significativamente o tempo de antecedência para a detecção de estados de preparação de erupção em comparação com as linhas de base escalares tradicionais.

Significância e Alegações
O artigo faz uma alegação modesta e específica: ele não afirma que a física do estado-topológico pode prever deterministicamente a data exata, a localização da ventilação ou a magnitude de todas as erupções vulcânicas.

• Previsão Probabilística: A principal contribuição é a validação de uma detecção de múltiplos canais probabilística de um "estado de preparação de erupção". O método identifica quando um sistema entrou em um estado de topologia incomumente organizado em relação às bases históricas.
• Alvo Operacional: O alvo de curto prazo apropriado não é a previsão determinística, mas a melhoria do tempo de antecedência e a redução de falsos alarmes escalares através do monitoramento do estado-topológico.
• Limitações e Próximos Passos: Os autores reconhecem que a auditoria atual é baseada apenas em sismicidade e depende de um único caso forte (Kilauea 2018). Eles afirmam que, para o método avançar para uma alegação operacional de alta precisão, requer-se:
  • Validação em intervalos de agitação não eruptivos independentes para estimar as taxas de falsos alarmes.
  • Integração de dados multicanais (deformação, gases, térmicos, inclinação).
  • Comparação contra modelos padrão (ex: b-value, ETAS).
  • Bloqueio estrito de parâmetros antes de testar em novos eventos.

Em conclusão, o artigo apoia a hipótese de que a preparação da erupção de 2018 do Kilauea se manifestou como uma reorganização topológica detectável da sismicidade, oferecendo um caminho potencial para a detecção de estado probabilística mais precoce.