Previsione dei Terremoti con ETAS.inlabru

Informazioni Fondamentali

• ID Articolo: 2510.13930
• Titolo: Earthquake Forecasting with ETAS.inlabru
• Autore: Ziwen Zhong
• Classificazione: stat.AP (Statistica - Applicazioni)
• Data di Pubblicazione: 15 ottobre 2025 (preprint arXiv)
• Link Articolo: https://arxiv.org/abs/2510.13930v1

Riassunto

Questo articolo affronta la previsione dei terremoti utilizzando il modello ETAS, proponendo il metodo "inlabru" basato su INLA per l'inferenza bayesiana al fine di risolvere i problemi di elevato costo computazionale e forte correlazione parametrica dei metodi MCMC tradizionali. L'articolo convalida l'efficacia del metodo attraverso tre esperimenti: (1) esplorazione dell'influenza dei parametri fissi sulle distribuzioni a posteriori di altri parametri, rivelando una chiara relazione di mutua influenza tra i parametri α e K; (2) confronto della distribuzione dei tempi normalizzati tra eventi da dati reali e cataloghi sintetici; (3) previsione degli eventi sismici entro dieci settimane dopo la scossa principale utilizzando eventi precedenti, valutando l'accuratezza predittiva mediante test di conteggio e punteggio di probabilità ordinale continua (CRPS).

Contesto di Ricerca e Motivazione

Definizione del Problema

I terremoti, come uno dei disastri naturali che minacciano la vita e i beni umani, presentano caratteristiche di improvvisità e difficoltà di evacuazione tempestiva. Sebbene la previsione della scossa principale rimanga l'ideale massimo della ricerca, sono stati compiuti progressi significativi nella previsione delle scosse di assestamento. La previsione affidabile delle scosse di assestamento aiuta i governi e i soccorritori a pianificare meglio le operazioni di soccorso post-sisma.

Importanza

L'accuratezza della previsione delle scosse di assestamento è direttamente correlata a:

1. Garantire maggiori opportunità di sopravvivenza alle persone colpite
2. Proteggere meglio la vita dei soccorritori
3. Ottimizzare l'allocazione delle risorse di soccorso post-sisma

Limitazioni dei Metodi Esistenti

1. Stima di Massima Verosimiglianza (MLE): il metodo è diretto e semplice, ma non può ottenere l'incertezza della stima parametrica
2. Catena di Markov Monte Carlo (MCMC): come metodo bayesiano può ottenere la distribuzione a posteriori dei parametri, ma presenta i seguenti problemi:
   • Elevato costo computazionale
   • Alta correlazione tra i parametri del modello
   • Difficoltà nel garantire una buona distribuzione a posteriori dei parametri

Motivazione della Ricerca

Il metodo 'inlabru' basato su Approssimazione Laplaciana Integrata Annidata (INLA) è in grado di:

• Ridurre significativamente il tempo di adattamento del modello mediante approssimazione della distribuzione a posteriori
• Prestazioni eccellenti quando la relazione di dipendenza tra parametri di covarianza e variabili latenti è approssimativamente lineare
• Stimare internamente l'approssimazione gaussiana, migliorando l'effetto dell'inferenza bayesiana

Contributi Principali

1. Contributo Metodologico: applicazione del metodo inlabru basato su INLA al modello ETAS per la previsione dei terremoti, risolvendo il problema dell'efficienza computazionale dei metodi MCMC tradizionali
2. Scoperta delle Relazioni Parametriche: rivelazione sistematica attraverso esperimenti delle relazioni di mutua influenza tra i parametri del modello ETAS, in particolare la forte correlazione tra i parametri α e K
3. Miglioramento dell'Efficienza Computazionale: dimostrazione che il fissaggio dei parametri α o K ai valori reali riduce il tempo di adattamento del modello di oltre il 50%
4. Verifica della Capacità Predittiva: creazione di un framework di valutazione predittiva basato su dati sismici reali, determinazione dei requisiti minimi di dati per previsioni affidabili
5. Intuizioni Teoriche: chiarimento delle condizioni necessarie per ottenere buone distribuzioni a posteriori dei parametri, fornendo orientamenti per l'applicazione pratica

Spiegazione Dettagliata del Metodo

Definizione del Compito

Input: sequenza di eventi sismici storici, contenente tempo di occorrenza t, posizione s e magnitudo m Output: distribuzione di probabilità predittiva degli eventi sismici nella finestra temporale futura Vincoli: framework del modello ETAS basato su processi puntuali auto-eccitati

Architettura del Modello

Fondamenti del Modello ETAS

Il modello ETAS è basato sul processo di Hawkes, la cui funzione di intensità è:

λ_ETAS(t,m|H_t) = (μ + Σ K·e^(α(m_h-M_0))·(t-t_h/c + 1)^(-p))·π(m)

Dove:

• μ: tasso di occorrenza di fondo
• K: tasso di produzione, controlla l'intensità della funzione di innesco
• α: scala di magnitudo, controlla l'influenza della magnitudo degli eventi storici
• c: offset temporale
• p: parametro di decadimento delle scosse di assestamento
• π(m): distribuzione di magnitudo che segue la legge di Gutenberg-Richter

Approssimazione della Verosimiglianza Logaritmica

Utilizzo dell'espansione di Taylor e strategia di binning per approssimare la verosimiglianza logaritmica:

L(θ|H) = -Λ_0(T_1,T_2) - Σ Σ Λ_i(t_j^(bi), t_{j+1}^(bi)) + Σ log λ(t_i|H_{t_i})

Procedura di Inferenza Bayesiana

1. Utilizzo dei parametri iniziali θ_0 per ottenere l'approssimazione della serie di Taylor della funzione di verosimiglianza logaritmica
2. Aggiunta al logaritmo della priorità, calcolo della moda della distribuzione a posteriori θ_1*
3. Aggiornamento dei parametri mediante ricerca lineare: θ* = αθ_0 + (1-α)θ_1*
4. Iterazione fino alla convergenza (differenza inferiore all'1%)

Punti di Innovazione Tecnica

1. Ottimizzazione della Strategia di Binning: utilizzo di binning adattivo, bin stretti sul lato sinistro per aumentare la precisione, bin larghi sul lato destro per ridurre il carico computazionale
2. Trasformazione Integrale di Probabilità Inversa: conversione di campioni da distribuzione gaussiana a campioni da distribuzione prioritaria target
3. Approssimazione Lineare: realizzazione dell'approssimazione efficiente della funzione di verosimiglianza mediante espansione di Taylor

Impostazione Sperimentale

Dataset

1. Terremoto dell'Aquila: utilizzato per generare cataloghi sintetici e analisi parametrica
2. Terremoto di Amatrice (2016-2017, Italia): sequenza sismica reale contenente tre scosse principali
3. Catalogo Sintetico: generato sulla base dei parametri del terremoto dell'Aquila, utilizzato per il controllo sperimentale

Metriche di Valutazione

1. Test di Conteggio (N-test): valutazione dell'accuratezza del numero di eventi predetti

   δ_2 = |{N_j|N_j ≤ N_obs, j=1,...,m}|/m
   

2. Punteggio di Probabilità Ordinale Continua (CRPS): valutazione simultanea di accuratezza e precisione

   S(F,N_true) = Σ(F(k) - I(N_true ≤ k))²
   

Metodi di Confronto

• Modello ETAS completo vs modelli con singolo parametro fisso
• Confronto di modelli con diverso numero di scosse principali storiche
• Confronto della distribuzione tra dati reali e dati sintetici

Dettagli di Implementazione

Distribuzioni Prioritarie:

• μ ~ Gamma(0.3, 0.6)
• K, α, c ~ U(0, 10)
• p ~ U(1, 10)

Valori Iniziali: μ₀=0.5, K₀=0.1, α₀=1, c₀=0.1, p₀=1.1

Risultati Sperimentali

Risultati Principali

Esperimento 1: Effetto del Fissaggio dei Parametri

Scoperte Chiave:

• Chiara relazione di mutua influenza tra α e K: il fissaggio di uno al valore reale migliora la distribuzione a posteriori dell'altro
• Relazione di mutua influenza anche tra c e p
• μ è relativamente indipendente, il fissaggio di altri parametri non influenza la distribuzione a posteriori di μ

Miglioramento dell'Efficienza Computazionale:

Esperimento 2: Distribuzione dei Tempi tra Eventi

• La distribuzione dei tempi normalizzati tra eventi da dati reali e cataloghi sintetici è sostanzialmente coerente
• Rispetto alla distribuzione esponenziale Exp(1), esiste un numero maggiore di eventi a intervalli brevi e lunghi, indicando l'esistenza di clustering
• L'aumento di μ rende la distribuzione più vicina a Exp(1) (processo di Poisson uniforme)
• L'aumento di p aumenta significativamente il numero di eventi a intervalli brevi

Esperimento 3: Valutazione delle Prestazioni Predittive

Risultati della Previsione del Terremoto di Amatrice:

• Utilizzo della seconda scossa principale: sovrapprevisione in tutte le dieci settimane, ma il valore reale rientra nell'intervallo di confidenza del 95%
• Utilizzo della terza scossa principale: valore δ₂ vicino a 0.5, previsione più accurata
• CRPS mostra che la previsione della terza scossa principale è significativamente migliore di quella della seconda

Esperimenti di Ablazione

Influenza del Numero di Scosse Principali Storiche:

1. Nessuna scossa principale: distribuzioni a posteriori estremamente scadenti per tutti i parametri eccetto μ
2. Una scossa principale: previsione di base utilizzabile, ma precisione limitata
3. Due scosse principali: miglioramento significativo delle prestazioni predittive, in particolare della distribuzione a posteriori del parametro α

Analisi di Casi

Ordine di Difficoltà nella Stima dei Parametri:

1. μ (Più Facile): stimabile con pochi eventi
2. c, p (Medio): richiedono più eventi sismici, stimabili in cataloghi tranquilli
3. α, K (Più Difficile): richiedono almeno una sequenza completa scossa principale-assestamento; più scosse principali migliorano significativamente la distribuzione a posteriori di α

Scoperte Sperimentali

1. Correlazione Parametrica: forte correlazione tra α-K e c-p, richiedono considerazione congiunta
2. Requisiti di Dati: previsione affidabile richiede almeno una sequenza completa scossa principale-assestamento
3. Miglioramento Predittivo: l'aumento del numero di scosse principali storiche migliora significativamente le prestazioni predittive
4. Ottimizzazione Computazionale: l'utilizzo ragionevole di informazioni prioritarie può migliorare significativamente l'efficienza computazionale

Lavori Correlati

Principali Direzioni di Ricerca

1. Sviluppo del Modello ETAS: dall'originale modello di Ogata (1988) a varie versioni estese
2. Metodi di Stima Parametrica: evoluzione da MLE a MCMC a metodi basati su INLA
3. Valutazione della Previsione Sismica: sviluppo di metriche di valutazione come test di conteggio e CRPS

Relazione di Questo Articolo con Lavori Correlati

• Costruito sulla base della teoria dei processi di Hawkes
• Eredita i concetti fondamentali del modello ETAS
• Approfondisce la ricerca sulla base dell'applicazione di inlabru di Serafini et al. (2023) e Naylor et al. (2023)

Vantaggi Relativi

1. Efficienza computazionale significativamente superiore ai metodi MCMC tradizionali
2. Analisi sistematica delle relazioni tra parametri
3. Fornisce un framework pratico di valutazione predittiva

Conclusioni e Discussione

Conclusioni Principali

1. Validità del Metodo: il metodo inlabru mostra prestazioni eccellenti nel modello ETAS, migliorando significativamente l'efficienza computazionale
2. Relazione Parametrica: la relazione di mutua dipendenza tra α e K è il fattore chiave che influenza le prestazioni del modello
3. Capacità Predittiva: il modello possiede affidabile capacità di previsione delle scosse di assestamento, con qualità predittiva che migliora con l'arricchimento dei dati storici
4. Requisiti di Dati: è necessaria almeno una sequenza scossa principale-assestamento per previsioni affidabili

Limitazioni

1. Limitazioni dei Dati: i cataloghi sintetici contengono solo una singola scossa principale, limitando la generalità delle conclusioni
2. Definizione dell'Intervallo tra Eventi: considera solo gli intervalli tra eventi adiacenti, non la struttura interna dei cluster di scosse principali
3. Sfide Computazionali: il tempo di adattamento del modello aumenta significativamente con l'aumento del numero di scosse principali, con difficoltà di convergenza
4. Analisi Insufficiente delle Caratteristiche: l'influenza di altre caratteristiche sismiche sulla distribuzione a posteriori, oltre al numero di scosse principali, non è stata sufficientemente esplorata

Direzioni Future

1. Strategie di Aggiornamento Parametrico: sviluppo di metodi di apprendimento incrementale, utilizzando la distribuzione a posteriori della scossa principale precedente come priorità per la successiva
2. Esperimenti Multi-Scossa Principale: ripetizione degli esperimenti su cataloghi sintetici contenenti più scosse principali
3. Analisi del Clustering degli Eventi: ricerca della distribuzione degli intervalli tra eventi all'interno e tra i cluster di scosse principali
4. Ricerca sull'Influenza delle Caratteristiche: analisi sistematica dell'influenza di varie caratteristiche sismiche sulle prestazioni del modello

Valutazione Approfondita

Punti di Forza

1. Innovazione Metodologica: applicazione riuscita del metodo INLA al modello ETAS, risolvendo problemi computazionali pratici
2. Progettazione Sperimentale Completa: tre esperimenti valutano completamente le prestazioni del metodo da diverse prospettive
3. Analisi Parametrica Approfondita: rivelazione sistematica delle relazioni di mutua influenza tra i parametri del modello ETAS
4. Alto Valore Pratico: fornisce orientamenti chiari sui requisiti di dati per l'applicazione pratica della previsione sismica
5. Alta Credibilità dei Risultati: verifica con dati sismici reali, conclusioni di significato pratico

Insufficienze

1. Analisi Teorica Insufficiente: mancanza di spiegazione teorica della correlazione parametrica
2. Ambito Sperimentale Limitato: principalmente basato su dati sismici italiani, rappresentatività geografica insufficiente
3. Confronto Metodologico Incompleto: mancanza di confronto dettagliato con altri metodi bayesiani moderni
4. Quantificazione dell'Incertezza: analisi dell'incertezza predittiva non sufficientemente approfondita

Impatto

1. Contributo Accademico: fornisce nuova metodologia al campo interdisciplinare della sismologia e della statistica
2. Valore Pratico: applicabile direttamente alla previsione sismica pratica e alla valutazione del rischio
3. Riproducibilità: basato sul pacchetto R open-source inlabru, buona riproducibilità
4. Potenziale di Diffusione: il metodo è estendibile ad altre applicazioni di processi puntuali auto-eccitati

Scenari Applicabili

1. Previsione delle Scosse di Assestamento: previsione dell'attività di assestamento a breve termine dopo il terremoto
2. Valutazione del Rischio Sismico: valutazione del rischio sismico basata su dati storici
3. Pianificazione della Risposta di Emergenza: supporto alle decisioni scientifiche per le operazioni di soccorso post-sisma
4. Ricerca Sismologica: comprensione della distribuzione spazio-temporale dell'attività sismica

Bibliografia

La bibliografia principale include:

• Ogata (1988): proposta originale del modello ETAS
• Hawkes (1971): fondamenti teorici del processo di Hawkes
• Rue et al. (2017): framework teorico del metodo INLA
• Serafini et al. (2023): applicazione di inlabru nei processi di Hawkes
• Naylor et al. (2023): modellazione bayesiana del modello ETAS

Valutazione Complessiva: questo è un articolo di importante valore pratico nel campo della previsione sismica, che risolve con successo il problema dell'efficienza computazionale del modello ETAS e fornisce attraverso esperimenti sistematici analisi parametrica e valutazione predittiva di valore. Sebbene vi sia spazio per miglioramenti nella profondità teorica e nell'ampiezza sperimentale, la praticità del metodo e l'operabilità delle conclusioni lo rendono un contributo importante in questo campo.