Retrieving the Seismic Moment Release Rate of Small Seismic Events Using the Adjoint Method - Theory and Numerical Validation

Quando la Terra trema, anche un minuscolo sussulto porta con sé una storia nascosta su come il terreno si sia fratturato. I sismologi sono in grado da tempo di individuare dove è iniziato un terremoto e quanta energia ha rilasciato, ma comprendere esattamente come tale energia sia stata rilasciata nel tempo è rimasto un enigma difficile. Questa storia è raccontata da quella che viene chiamata funzione del tempo sorgente, un registro del battito cardiaco della rottura, dal suo primo scricchiolio alla sua interruzione finale. Conoscere questo ritmo è fondamentale perché rivela la velocità della rottura, la direzione per cui si è propagata e i cambiamenti di stress avvenuti nelle profondità del sottosuolo. Tuttavia, estrarre questo ritmo dalle onde disordinate registrate in superficie è notoriamente difficile. I metodi tradizionali si affidano o alla ricerca di un terremoto minore e simile da utilizzare come riferimento — un processo che può essere soggettivo e impreciso — o tentano di eseguire simulazioni massicce e computazionalmente pesanti che sono troppo lente per un uso di routine.

In uno studio recente, i ricercatori Wojciech Dębski e Kamil Waśkiewicz dell'Istituto di Geofisica dell'Accademia Polacca delle Scienze hanno proposto un modo diverso per ascoltare questi ritmi nascosti. Hanno sviluppato una tecnica basata sull'idea fisica semplice che il suono e le onde sismiche possano viaggiare all'indietro nel tempo proprio come viaggiano in avanti. Prendendo le registrazioni dalla superficie, invertendole nel tempo e inviandole nuovamente nel terreno attraverso un modello computazionale, le onde tornano naturalmente a convergere nel punto in cui l' terremoto è iniziato. Mentre queste onde invertite convergono, ricostruiscono il modello originale di rilascio dell'energia. Il team ha testato questa idea utilizzando simulazioni al computer di piccoli terremoti causati dall'uomo, come quelli innescati dall'attività mineraria, utilizzando dati sintetici privi di rumore. I loro risultati dimostrano che questo metodo può recuperare con successo la tempistica dettagliata della rozione senza necessitare di complesse inversioni matematiche o della ricerca di un terremoto di riferimento. Offre un modo più veloce e robusto per comprendere la meccanica dei piccoli eventi sismici, trasformando un difficile problema matematico in un processo fisico diretto.

Il nucleo di questo nuovo approccio si basa su una proprietà fondamentale della fisica delle onde: la reversibilità temporale. In un mezzo perfettamente elastico, dove l'energia non viene persa in calore o attrito, le equazioni che descrivono il movimento delle onde funzionano allo stesso modo sia che il tempo proceda in avanti che all'indietro. Immaginate un video di una pietra che cade in uno stagno: se riavvolgete il video, le increspature viaggiano verso l'interno e la pietra balza fuori dall'acqua. Sebbene non sia possibile invertire realmente il tempo nel mondo reale, la matematica permette di simulare questo processo. I ricercatori hanno utilizzato questo concetto per creare un algoritmo che chiamano TRSTF. Il processo inizia prendendo le onde sismiche registrate in varie stazioni in superficie. Queste onde vengono poi invertite nel tempo e reintrodotte in un modello computazionale della crosta terrestre. Mentre le onde invertite viaggiano attraverso il modello, rimbalzano sui confini e interagiscono con la complessa geologia, per poi convergere infine nel luogo di origine della sorgente.

Quando queste onde invertite si incontrano nella sorgente, interferiscono tra loro. Nell'esatta posizione del terremoto, le onde si allineano perfettamente, creando un segnale forte e focalizzato. In qualsiasi altra posizione, le onde arrivano fuori fase e si annullano a vicenda, lasciando solo un debole rumore. Questo effetto di focalizzazione è la chiave del metodo. Il segnale che emerge nel punto focale è una ricostruzione diretta della funzione del tempo sorgente, ovvero il ritmo del rilascio di energia del terremoto. I ricercatori hanno scoperto che questa tecnica funziona molto bene nei loro esperimenti sintetici privi di rumore, anche quando il modello dell'interno della Terra non è perfettamente noto. Nei loro esperimenti al computer, hanno simulato tre diversi scenari per mimare le configurazioni di monitoraggio del mondo reale: uno in cui i sensori erano posizionati vicino al livello dell'attività mineraria, uno in cui i sensori erano sulla superficie sopra l'evento, e un terzo con una diversa disposizione superficiale. In ogni caso, il metodo ha recuperato con successo la forma del modello di rottura originale.

Lo studio ha anche esplorato come gli errori nel modello della velocità del suono della Terra influenzino i risultati. Nel mondo reale, i geologi non conoscono l'esatta velocità delle onde sismiche in ogni punto del sottosuolo; essi si affidano a stime. I ricercatori hanno testato il loro metodo introducendo errori casuali e variazioni a blocchi nei loro modelli computazionali, simulando una situazione in cui il modello di velocità era accurato solo per circa il 90 percento. I risultati hanno indicato che l'uso di modelli di velocità inaccurati introduce errori sistematici nella funzione del tempo sorgente ricostruita. Il metodo si è dimostrato più sensibile agli errori organizzati e a blocchi che a quelli casuali e sparsi. Tuttavia, nonostante questi errori, il segnale recuperato rimaneva riconoscibile e forniva informazioni utili, dimostrando che il metodo mantiene la sua utilità anche con modelli imperfetti, sebbene l'accuratezza non sia immune dalle imprecisioni del modello.

Una delle scoperte più sorprendenti è stata come la disposizione dei sensori abbia influenzato la qualità del risultato. I ricercatori hanno scoperto che avere sensori più vicini ai confini dell'area oggetto di studio favoriva la ricostruzione. Questo perché le riflessioni da questi confini agiscono come ulteriori sensori virtuali, fornendo maggiori informazioni al sistema. In uno dei loro casi di test, in cui i sensori erano posizionati vicino al bordo superiore del modello, il segnale ricostruito era più nitido e accurato rispetto alle altre configurazioni, nonostante i dati grezzi apparissero più distorti a causa dell'interferenza. Questo risultato controintuitivo suggerisce che il metodo trae beneficio dalla complessità del campo d'onda, utilizzando le riflessioni a proprio vantaggio invece di trattarle come rumore.

I ricercatori hanno anche misurato caratteristiche specifiche dei segnali ricostruiti, come la rapidità con cui la rottura è iniziata e quanto è durata. Hanno scoperto che, sebbene i segnali ricostruiti fossero leggermente più ampi degli originali, probabilmente a causa del numero limitato di sensori utilizzati nella simulazione, la forma generale e la tempistica venivano preservate. La capacità di stimare la durata della rottura è vitale perché aiuta gli scienziati a calcolare lo stress drop dinamico, una misura di quanto stress sia stato rilasciato durante l'evento. Questa informazione è essenziale per comprendere la meccanica della sismicità indotta, come i terremoti causati dall'attività mineraria o dall'iniezione di fluidi. Lo studio ha dimostrato che il nuovo metodo può stimare questi parametri con un'accuratezza comparabile, o migliore, rispetto alle tecniche esistenti, ma con un impegno computazionale molto minore.

A differenza della inversione dell'intero campo d'onda (full-waveform inversion), che richiede l'esecuzione di migliaia di simulazioni e complessi cicli di ottimizzazione per trovare il miglior adattamento, questo nuovo metodo richiede un unico passaggio delle onde invertite attraverso il modello. Ciò lo rende computazionalmente efficiente, capace di girare su normali computer desktop per modelli piccoli o su grandi supercomputer per scale regionali. Evita inoltre la necessità di selezionare un terremoto di riferimento, un passaggio che può introdurre pregiudizi umani e incertezza nei metodi tradizionali. I ricercatori sottolineano che il loro approccio è particolarmente adatto per i piccoli eventi sismici, dove la sorgente è abbastanza piccola da poter essere trattata come un singolo punto. Essi riconoscono che il metodo assume che il terreno si comporti come un materiale elastico perfetto, il che significa che non assorbe energia, e che la sorgente è piccola rispetto alla lunghezza d'onda delle onde. Queste sono ipotesi ragionevoli per molti eventi legati all'attività mineraria, ma i ricercatori notano che sono necessari ulteriori studi per comprendere come il metodo si comporti in ambienti più complessi e reali con una significativa perdita di energia.

L'articolo conclude evidenziando il potenziale di questa tecnica nel diventare uno strumento standard nella sismologia. Fornendo un modo veloce, affidabile e fisicamente intuitivo per recuperare la funzione del tempo sorgente, essa apre la strada all'analisi di routine dei piccoli eventi sismici che prima erano troppo difficili da caratterizzare in dettaglio. I ricercatori suggeriscono che questo metodo potrebbe essere integrato nei sistemi di monitoraggio esistenti per fornire intuizioni immediate sulla natura della sismicità indotta. Sebbene l'attuale studio sia stato limitato alle simulazioni al computer, la base teorica è solida e i risultati sono promettenti. Il passo successivo, secondo gli autori, sarà testare il metodo su dati sismici reali e affinare la comprensione dei suoi limiti e delle sue incertezze. Per ora, lo studio rappresenta una prova di concetto che il tempo reversibile può essere una potente lente per guardare nel passato della Terra, rivelando la storia nascosta di come il terreno si rompe.