Teoria di base e applicazioni della tomografia sismica
La tomografia sismica è un metodo geofisico utilizzato per "fotografare" l'interno della Terra tramite onde sismiche. Analogamente alla TAC in medicina, che crea un'immagine tridimensionale del corpo a partire dai raggi X, la tomografia sismica modella le strutture del sottosuolo utilizzando i dati relativi al tempo di propagazione, all'ampiezza o alla forma d'onda registrati da una rete di sismometri. Questo metodo è fondamentale perché l'interno della Terra non può essere osservato direttamente, mentre dinamiche come la tettonica a placche, l'attività vulcanica e i terremoti sono fortemente influenzate dalle variazioni delle proprietà fisiche delle rocce in profondità.
1. Concetto base delle onde sismiche
Le onde sismiche sono onde elastiche che si propagano attraverso la roccia. Le onde sismiche si dividono generalmente in:
1. Onde del corpo
– Onde P (primarie/di compressione): si propagano più velocemente, possono attraversare solidi e liquidi e sono sensibili alle variazioni di velocità di compressione e densità.
– Onde S (secondarie/di taglio): più lente delle onde P, non possono propagarsi nei fluidi, sensibili alla rigidità (modulo di taglio) della roccia.
2. Onde superficiali
Ad esempio, le onde di Rayleigh e di Love, che sono generalmente dominanti nelle registrazioni di terremoti distanti, presentano dispersione e sono molto informative per la struttura della litosfera fino all'astenosfera superficiale.
La velocità delle onde sismiche dipende da parametri elastici (modulo di elasticità) e densità. Variazioni di temperatura, composizione minerale, pressione, porosità, fratture e presenza di fluidi o fusi parziali influenzano la velocità. Questa è la base fisica per la mappatura delle variazioni di velocità, che fornisce indizi sulle condizioni geologiche del sottosuolo.
2. Principi di tomografia: problema in avanti e problema all'indietro
La tomografia sismica si basa su due concetti di calcolo:
a) Problema in avanti (problema in avanti)
Dato un modello della Terra (ad esempio, una distribuzione della velocità delle onde), possiamo calcolare una previsione dei dati sismici: il tempo di percorrenza di un'onda da una sorgente (terremoto o sorgente artificiale) a una stazione. In un semplice approccio basato sulla teoria dei raggi, si assume che il percorso dell'onda sia un raggio che segue il principio di Fermat: l'onda sceglie il percorso con il tempo di percorrenza minimo.
Matematicamente il tempo di percorrenza \(T\) può essere scritto come:
\[
T = \int_{\text{raggio}} \frac{ds}{v(\mathbf{x})}
\]
dove \(v(\mathbf{x})\) è la velocità dell'onda nella posizione \(\mathbf{x}\), e \(ds\) è l'elemento di lunghezza del percorso.
b) Problema inverso
Il contrario è vero: disponiamo di dati osservazionali (tempi di percorrenza, residui dei tempi di percorrenza, dispersioni delle onde superficiali o forme d'onda) e vogliamo stimare il modello di velocità che meglio si adatta ai dati. Il problema inverso è generalmente non univoco e mal posto: più modelli possono adattarsi ugualmente bene e i dati contengono incertezza/rumore. Pertanto, è necessaria una regolarizzazione, come ad esempio la levigatura, lo smorzamento o i vincoli geologici.
In pratica, l'inversione viene spesso eseguita in modo linearizzato: il modello iniziale viene leggermente perturbato e quindi la differenza di tempo di percorrenza viene calcolata come:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
dove \(\delta s = \delta(1/v)\) è la perturbazione di lentezza. Queste equazioni vengono organizzate in un grande sistema lineare \( \mathbf{d} = \mathbf{Gm} \), quindi risolto con il metodo dei minimi quadrati smorzati o sue varianti.
3. Tipi di tomografia sismica
1) Tomografia del tempo di percorrenza
Il metodo più comune utilizza gli arrivi delle fasi P e S di numerosi terremoti. I dati sono residui dei tempi di propagazione rispetto a un modello di riferimento (ad esempio, IASP91 o ak135 su scala globale). Adatto per la modellazione della velocità 3D nella crosta e nel mantello superiore, soprattutto se la rete di stazioni è densa.
2) Tomografia a onde superficiali
Sfruttando la dispersione (variazione della velocità di fase/gruppo in funzione del periodo) delle onde di Rayleigh/Love, questa tomografia, sensibile alle strutture superficiali e intermedie, è molto utile per mappare lo spessore della litosfera, le zone a bassa velocità (LVZ) e le variazioni di temperatura.
3) Tomografia telesismica (tomografia telesismica)
Si utilizzano terremoti distanti (telesismici) le cui onde si propagano attraverso il mantello e vengono poi registrate da una rete locale. Il vantaggio: hanno sorgenti multiple e provengono da diverse direzioni, contribuendo così a "illuminare" il volume sotterraneo dell'area di studio, ad esempio al di sotto di vulcani o zone di subduzione.
4) Tomografia della forma d'onda / inversione completa della forma d'onda (FWI)
L'utilizzo dell'intera forma d'onda, e non solo del tempo di arrivo, consente teoricamente di ottenere un'elevata risoluzione perché sfrutta appieno le informazioni di ampiezza e fase, ma richiede notevoli risorse computazionali e un buon modello iniziale per evitare trappole dovute a minimi locali.
4. Fasi generali di elaborazione e inversione
1. Acquisizione dei dati
Raccolta di registrazioni provenienti da sismometri permanenti o temporanei. Su scala esplorativa, la sorgente può essere un'esplosione o un fenomeno vibrosismico; su scala regionale-globale, la sorgente primaria è un terremoto.
2. Identificazione e selezione della fase
Determinazione del tempo di arrivo delle onde P, S o di superficie. La qualità della selezione determina la qualità del modello.
3. Correzione e modellazione iniziali
Correzione temporale (deriva dell'orologio), correzione dell'elevazione e selezione del modello di riferimento. Il modello iniziale può essere monodimensionale o tridimensionale semplice.
4. Tracciamento dei raggi / simulazione delle onde
Calcola il percorso dei raggi o delle forme d'onda sintetiche per costruire la matrice di sensibilità.
5. Inversione e regolarizzazione
Risolvere il sistema per ottenere un modello di velocità. La regolarizzazione viene scelta in base all'obiettivo: evidenziare anomalie marcate o smussare le tendenze.
6. Valutazione della risoluzione e test di affidabilità
Ad esempio, con il test a scacchiera, il test a picco o l'analisi della funzione di covarianza/dispersione del punto per vedere quali parti sono effettivamente risolte dai dati.
5. Interpretazione del modello di velocità
Il risultato della tomografia è solitamente una mappa delle anomalie di velocità relativa: le zone ad alta velocità sono spesso interpretate come rocce più fredde, più dense o più rigide (ad esempio, placche in subduzione). Le zone a bassa velocità possono indicare temperature elevate, rocce alterate, zone di frattura riempite di fluidi o fusione parziale, spesso associate a vulcani attivi o astenosfera calda.
Tuttavia, l'interpretazione deve essere cauta perché la velocità è influenzata da molti fattori. Idealmente, la tomografia sismica viene combinata con altri dati come quelli gravimetrici, magnetotellurici (MT), geodetici, geologici di superficie e petrologici.
6. Applicazioni della tomografia sismica
a) Studio delle zone di subduzione e della dinamica delle placche
La tomografia globale e regionale permette di mappare la placca del mantello in subduzione, inclusi la sua geometria, la profondità e la segmentazione. Queste informazioni sono fondamentali per comprendere le sorgenti dei grandi terremoti, i meccanismi di accoppiamento delle placche e l'evoluzione tettonica di una regione.
b) Sistemi vulcanici e mitigazione dei disastri
Nelle aree vulcaniche, la tomografia può identificare zone a bassa velocità associate a camere magmatiche, vie di risalita dei fluidi o rocce alterate. Con un monitoraggio regolare, le variazioni di velocità (tomografia time-lapse) possono fornire indicatori del movimento di fluidi/magmi rilevanti per l'allerta precoce delle eruzioni.
c) Esplorazione di energia e risorse
A livello della crosta terrestre superficiale, la tomografia viene utilizzata nell'esplorazione di idrocarburi, geotermia e miniere. Le variazioni di velocità aiutano a mappare la litologia, le strutture di faglia, le zone di frattura e i giacimenti. In ambito geotermico, l'integrazione della tomografia con la tomografia magnetotellurica (MT) è spesso efficace: la tomografia fornisce informazioni sull'elasticità, mentre la MT fornisce informazioni sulla conduttività dei fluidi.
d) Caratterizzazione delle faglie attive e dei rischi sismici
La tomografia locale può rivelare zone di debolezza, zone danneggiate ed eterogeneità intorno alle faglie. Questo aiuta a comprendere la segmentazione delle faglie, il potenziale blocco e le variazioni di velocità che influenzano l'amplificazione delle onde (effetti di sito).
e) Struttura regionale della crosta e della litosfera
Utilizzando la tomografia a onde superficiali, i ricercatori possono mappare lo spessore della litosfera, il confine di Moho e le anomalie astenosferiche. I risultati costituiscono la base per i modelli geodinamici, tra cui la formazione dei bacini, l'orogenesi e l'evoluzione continentale.
7. Limitazioni e sfide
La tomografia sismica dipende fortemente dalla distribuzione delle sorgenti e delle stazioni. Le aree con pochi terremoti o stazioni sparse avranno una scarsa illuminazione, con conseguente bassa risoluzione. Inoltre, le ipotesi della teoria dei raggi possono essere meno accurate in presenza di forte eterogeneità o per determinate frequenze. Anche il rumore, gli errori di rilevamento e l'incertezza nella localizzazione dell'ipocentro possono influenzare i risultati dell'inversione. Pertanto, la valutazione della risoluzione e l'integrazione multimetodo sono cruciali per evitare interpretazioni errate.
Chiusura
La tomografia sismica è uno strumento fondamentale per studiare la struttura interna della Terra, dalla scala locale a quella globale. Il suo fondamento teorico si basa sulla relazione tra le proprietà elastiche delle rocce e la propagazione delle onde, nonché sulla risoluzione di problemi inversi che richiedono regolarizzazione e test di risoluzione. Le sue applicazioni sono molteplici: dalla mappatura delle placche in subduzione e dei sistemi magmatici vulcanici all'esplorazione energetica e alla valutazione del rischio sismico. Grazie ai progressi nelle reti sismiche, nell'informatica e nei metodi di inversione (inclusa l'inversione completa della forma d'onda), la tomografia sismica è sempre più in grado di produrre immagini più nitide e informative, sia per scopi scientifici che per la mitigazione dei disastri.
Se lo desideri, posso adattare questo articolo in una versione più accademica (con citazioni e bibliografia) oppure concentrarmi su una delle applicazioni (ad esempio, la tomografia per i vulcani o per l'energia geotermica).
