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

Sintesi Tecnica: Monitoraggio della Fase Topologica della Preparazione dell'Eruzione del Kilauea

Definizione del Problema
La previsione ad alta precisione delle eruzioni vulcaniche rimane una sfida significativa perché la preparazione eruttiva raramente è segnalata da una singola variabile scalare (ad esempio, un picco nel conteggio dei terremoti). Al contrario, essa comporta una riorganizzazione distribuita della sismicità, della deformazione, del trasporto di gas e della connettività tra magma e fratture. Il monitoraggio tradizionale spesso si basa su canali isolati, che possono risultare ambigui; ad esempio, gli sciami sismici non portano sempre a eruzioni, e l'aumento del numero di eventi può riflettere un assestamento dello stress piuttosto che un'attività imminente. Il documento postula che la preparazione dell'eruzione possa essere meglio compresa come una transizione nello "stato-topologico" del campo sismico — passando da un'attività di fondo diffusa a uno stato connesso, migratorio e anisotropo.

Metodologia
Lo studio esegue un audit riproducibile, basato esclusivamente sulla sismicità, del sistema vulcanico del Kilauea utilizzando i dati pubblici del catalogo sismico USGS FDSN.

• Ambito dei Dati: L'analisi copre un raggio di 80 km attorno al Kilauea (19.4069°N, 155.2834°O) per terremoti con magnitudo $M \ge 1.0$.
• Finestre Temporali: I dati sono stati aggregati per finestre dello stesso calendario (dal 1° gennaio al 5 maggio) dal 2014 al 2019. Le finestre 2014–2017 hanno servito come distribuzione di riferimento per l'addestramento, mentre la finestra 2018 è stata il caso di test target per l'eruzione della zona di rift orientale inferiore (LERZ) (inizio: 3 maggio 2018).
• Vettore dello Stato-Topologico: Una finestra mobile di sette giorni è stata convertita in un vettore di sette caratteristiche:
  1. Conteggio degli eventi ($N$).
  2. Prossimità dell'energia della magnitudo ($E$).
  3. Frazione della componente connessa più grande ($C_{max}$).
  4. Densità dei bordi del grafo ($\rho_G$).
  5. Entropia spaziale normalizzata ($S_{sp}$).
  6. Anisotropia della nuvola di eventi ($A$).
  7. Migrazione del baricentro ($\Delta r_c$).
  • Costruzione del Grafo: Gli eventi sono stati collegati se separati da $\le 8$ km nello spazio proiettato e $\le 2$ giorni nel tempo.
• Costruzione dell'Indicatore: È stato calcolato un Indicatore di Fase-Topologica (TPS) composito come la media delle componenti standardizzate. L'entropia spaziale è stata negata per riflettere la transizione da un'organizzazione diffusa a una concentrata. Le soglie sono state impostate al 99° percentile della distribuzione di addestramento 2014–2017.

Risultati Chiave
L'audit, che utilizza parametri bloccati addestrati sui dati 2014–2017, ha prodotto i seguenti risultati riguardanti l'eruzione LERZ del 2018:

• Preavviso Precoce: L'indicatore composito della fase topologica ha superato per la prima volta la soglia del 99° percentile di addestramento il 12 aprile 2018, fornendo un tempo di anticipo di 21 giorni prima dell'inizio dell'eruzione del 3 maggio.
• Confronto con Metriche Scalari:
  • Conteggio degli Eventi: Ha superato il suo 99° percentile solo il 1° maggio 2018 (2 giorni prima dell'inizio).
  • Prossimità dell'Energia: Non ha superato la sua soglia prima dell'inizio.
• Analisi delle Componenti: L'avviso precoce è stato guidato da una riorganizzazione strutturale coordinata. La densità del grafo ha violato la soglia il 12 aprile, mentre la frazione della componente connessa più grande e l'entropia spaziale hanno violato le rispettive soglie il 17 aprile. La migrazione ha mostrato una violazione precoce a marzo ma non ha superato il massimo di addestramento fino al 1° maggio. L'anisotropia ha violato la soglia precocemente a gennaio ma è stata ritenuta insufficiente come predittore autonomo.
• Valori di Picco: L'indicatore topologico ha raggiunto un valore di picco di 4.0411 il 1° maggio 2018, superando sia il 99° percentile di addestramento (2.0462) che il massimo di addestramento (2.8166).

Contributi Chiave

1. Rilevamento dello Stato-Topologico: Lo studio dimostra che la preparazione dell'eruzione in un sistema basaltico ben monitorato può essere rilevata come una transizione topologica nel campo sismico, distinta da un semplice aumento del conteggio degli eventi o dell'energia.
2. Audit Riproducibile: Gli autori forniscono un flusso di lavoro completamente riproducibile, inclusi script e file di dati specifici, che consentono la verifica indipendente dell'audit a "parametri bloccati".
3. Miglioramento del Tempo di Anticipo: I risultati suggeriscono che la combinazione di metriche di connettività, densità e organizzazione spaziale può estendere significativamente il tempo di anticipo per il rilevamento degli stati di preparazione dell'eruzione rispetto alle tradizionali basi scalari.

Significato e Rivendicazioni
Il documento esprime una rivendicazione modesta e specifica: non afferma che la fisica dello stato-topologico possa prevedere deterministicamente la data esatta, la posizione del condotto o la magnitudo di tutte le eruzioni vulcaniche.

• Preavviso Probabilistico: Il contributo primario è la validazione di un rilevamento probabilistico e multicanale di uno "stato di preparazione all'eruzione". Il metodo identifica quando un sistema è entrato in uno stato topologico insolitamente organizzato rispetto alle linee di base storiche.
• Obiettivo Operativo: L'obiettivo a breve termine non è la previsione deterministica, ma il miglioramento del tempo di anticipo e la riduzione dei falsi allarmi scalari attraverso il monitoraggio dello stato-topologico.
• Limitazioni e Passi Futuri: Gli autori riconoscono che l'attuale audit è basato solo sulla sismicità e si basa su un singolo caso forte (Kilauea 2018). Affermano che, affinché il metodo progredisca verso una rivendicazione operativa ad alta precisione, è necessario:
  • Validare su intervalli di inquietudine non eruttivi indipendenti per stimare i tassi di falsi allarmi.
  • Integrare dati multicanale (deformazione, gas, termico, inclinazione).
  • Confrontarsi con modelli standard (es. b-value, ETAS).
  • Applicare un rigoroso blocco dei parametri prima di testarli su nuovi eventi.

In conclusione, il documento supporta l'ipotesi che la preparazione dell'eruzione del Kilauea del 2018 si sia manifestata come una riorganizzazione topologica rilevabile della sismicità, offrendo una potenziale via per il rilevamento precoce e probabilistico dello stato.