In questo secondo contributo alla conoscenza del terremoto del 2009 vediamo come sia stato possibile ricostruire il processo di rottura associato alla scossa principale del 6 aprile alle ore 3:32, attraverso l’utilizzo di dati geodetici e sismologici.
Deformazioni e modelli di faglia geodetici
La deformazione della crosta terrestre causata dalla sequenza sismica dell’Aquila del 2009 – spiegano Daniele Cheloni (INGV-ONT) e Antonella Cirella (INGV-Rm1) – è stata misurata sia da stazioni GPS collocate a terra in un’ampia regione dell’Italia centrale (Anzidei et al., 2009; Cheloni et al., 2010; D’Agostino et al., 2012; Gualandi et al., 2014), sia dai satelliti con le tecniche radar (Atzori et al., 2009; Walters et al., 2009), e sia da tecniche di geodesia terrestre (Cheloni et al., 2014).
Spostamenti cosismici orizzontali misurati (frecce nere) e simulati (frecce rosse) ai caposaldi GPS presenti nella zona epicentrale. Il rettangolo blu rappresenta il modello di faglia (proiezione sulla superficie) ricavato da tali spostamenti. Il riquadro in basso a destra mostra invece gli spostamenti verticali osservati e modellati (Anzidei et al., 2009.)
Tali tecniche hanno permesso di evidenziare gli spostamenti della superficie terrestre e realizzare così un modello di faglia (posizione, estensione, spostamento dei due lembi della faglia) per l’evento principale del 6 aprile del 2009 (Anzidei et al., 2009; Atzori et al., 2009, Cheloni et al., 2010; Walters et al 2009).
In particolare i massimi spostamenti orizzontali e verticali osservati furono dell’ordine dei 10 e 15 cm alla stazione GPS denominata CADO.
Gli interferogrammi (ottenuti dai sensori radar satellitari in passaggi prima e dopo il terremoto) mostrano un campo di deformazione caratterizzato da frange concentriche (Figura 2) indicanti che il terreno si era allontanato dal satellite (lungo la linea di vista – Light Of Sight – del satellite) su un’area di circa 480 km2 estesa verso sud-est a partire dall’epicentro del terremoto del 6 aprile. I massimi spostamenti osservati lungo la linea di vista dei satelliti furono dell’ordine dei 20-28 cm, misurati tra la città dell’Aquila e Fossa (AQ).
Interferogrammi calcolati con i satelliti (a) COSMO-Sky-Med e Envisat (b and c) per la scossa del 6 aprile. La stella rossa indica l’epicentro del terremoto del 6 aprile (Atzori et al., 2009).
Per la prima volta in Italia, sono stati inoltre osservati dei movimenti lenti della superficie terrestre nell’area circostante l’epicentro, dovuti al movimento post-terremoto che avviene sul piano di faglia successivamente al terremoto (Cheloni et al., 2010; D’Agostino et al., 2012; Cheloni et al., 2014; Gualandi et al., 2014). Definiamo “lento” questo spostamento perché avviene in un arco temporale di molte settimane o alcuni mesi, mentre durante il terremoto lo spostamento di tutta la faglia avviene in pochi secondi, come vedremo più avanti. Questo movimento lento viene definito afterslip (scivolamento post-sismico). Le registrazioni giornaliere GPS disponibili mostrano chiaramente il lento movimento avvenuto nei giorni successivi alla scossa principale (Figura 3).
Spostamento misurato alle stazioni GPS dell’Aquila (AQUI) e di Paganica (PAGA) durante (freccia rossa) e nelle settimane successive (freccia blu) al terremoto del 6 aprile 2009. Sono mostrate le tre componenti del movimento (North, East, Up) (Cheloni et al., 2010).
In particolare, i dati GPS misurati unitamente agli interferogrammi mostrati sopra, misurati prima e dopo il terremoto principale del 6 aprile, hanno permesso di calcolare lo spostamento del terreno e ricavare quindi un modello di faglia per la sequenza del 2009. I principali modelli di faglia sono stati proposti da Anzidei et al. (2009), Atzori et al. (2009), Walters et al. (2009), Cheloni et al. (2010) e Gualandi et al. (2014). Altri modelli vennero proposti da Balestra et al. (2015) e Castaldo et al. (2018).
Gli spostamenti cosismici medi sull’intero piano di faglia ottenuti dall’inversione dei dati GPS (Anzidei et al., 2009; Cheloni et al., 2010; Gualandi et al., 2014), sono di circa 50-60 cm, in accordo con l’inversione di dati interferometrici (Atzori et al., 2009, Walters et al., 2009), con massimi movimenti di circa 1 metro (Figura 4).
Modello di faglia e distribuzione di movimento (slip) sul piano di faglia stimato da misure di spostamento GPS. Le frecce nere rappresentano gli spostamenti osservati, mentre quelle bianche gli spostamenti previsti dal modello. La scala di colori rappresenta l’entità di movimento stimata sul piano di faglia (Gualandi et al., 2014).
In generale,- spiegano Daniele Cheloni (INGV-ONT) e Antonella Cirella (INGV-Rm1) – i vari modelli di faglia proposti per la scossa del 6 aprile 2009, concordano nel definire come sorgente sismogenetica della sequenza dell’Aquila del 2009 la faglia di Paganica. Tutti i modelli la caratterizzano come una faglia con geometria planare con un angolo di immersione (dip) verso SW di circa 50° fino ad una profondità di circa 10 km, attivatasi per una lunghezza di circa 16 km. La geometria della faglia identificata dai dati sopra descritti è in accordo con quanto mostrato dalla distribuzione in mappa e in profondità delle repliche (o aftershocks).
Ulteriori dettagli sul processo di fagliazione, in particolare quelli legati all’evoluzione temporale della rottura durante il terremoto del 6 aprile, possono essere ricavati dai dati accelerometrici registrati dalle stazioni poste in area epicentrale, come illustrato nel seguito.
Modellazione congiunta sismologica / geodetica
Il terremoto del 6 aprile 2009 e i principali eventi della sequenza sismica ad esso associati sono stati registrati da diverse stazioni digitali appartenenti alla “Rete Accelerometrica Nazionale” (RAN) gestita dal Dipartimento della Protezione Civile, da diverse stazioni accelerometriche a larga banda della Rete MedNet e dalle stazioni sismiche permanenti digitali della Rete Sismica Nazionale Italiana dell’INGV . Subito dopo l’evento principale, l’INGV in collaborazione con il Laboratoire de Géophysique Interne et Tectonophysique (LGIT) di Grenoble ha provveduto all’installazione di una fitta rete temporanea composta di ulteriori 40 stazioni sismiche digitali (Chiaraluce et al., 2011), che ha permesso la registrazione dell’intera sequenza. I dati accelerometrici registrati dalle stazioni dislocate nella regione epicentrale (Figura 5) durante il terremoto dell’Aquila del 2009, costituiscono per la comunità scientifica un set di osservazioni senza precedenti per un evento con meccanismo di faglia normale.
Mappa del terremoto dell’Aquila: – il rettangolo in rosso rappresenta la proiezione in superficie del piano di faglia; – i triangoli bianchi indicano le stazioni accelerometriche strong-motion e i punti in viola i siti GPS selezionati nello studio di Cirella et al., (2012) per la modellazione della sorgente sismica. L’epicentro del terremoto è identificato dalla stella rossa mentre in giallo vengono riportate le posizioni della città de L’Aquila e di Paganica.
L’analisi e la modellazione della radiazione sismica associata al terremoto dell’Aquila, eseguite attraverso l’applicazione di metodologie che tengono conto degli effetti dovuti alla vicinanza della sorgente sismogenetica, ha permesso di ricostruire l’evoluzione spazio-temporale della rottura co-sismica avvenuta sulla faglia di Paganica, responsabile del terremoto. In particolare, l’inversione congiunta di dati geodetici (GPS e DInSAR) e dati sismologici (Cirella et al., 2009; Yano et al., 2009; Cirella et al., 2012; Gallovič et al., 2015; Del Gaudio et al., 2015) ha consentito di ottenere una descrizione dettagliata del processo di sorgente sismica, in termini di distribuzioni dei parametri cinematici (picco della velocità di dislocazione, velocità del fronte di rottura, durata e direzione della dislocazione) sul piano di faglia.
La Figura 6a mostra il modello di rottura ottenuto per il terremoto dell’Aquila del 6 aprile 2009, da Cirella et al. (2012). I pannelli in alto, al centro e in basso mostrano, rispettivamente, le distribuzioni di dislocazione, la durata ed il picco della velocità di dislocazione sul piano di faglia. Le isolinee in bianco rappresentano i tempi di rottura e i vettori in nero corrispondono alla direzione di dislocazione. La stella rossa identifica la posizione dell’ipocentro. In Figura 6b si ha il confronto tra i sismogrammi osservati (in blu) e modellati (in rosso) alle stazioni riportate in Figura 5.
a) Modello della sorgente sismica responsabile del terremoto dell’Aquila, descritto in termini di distribuzioni della dislocazione sul piano di faglia (in alto), durata (centro) e picco (in basso) della velocità di dislocazione sul piano di faglia. b) Confronto tra le forme d’onda osservate (blu) e modellate (rosso). I numeri indicano i valori di picco, in cm/s, osservati su ciascuna stazione e per ogni componente del moto.
Il video mostra l’evoluzione temporale della velocità di dislocazione (in m/s) sul piano di faglia proiettata sulla superficie terrestre. I punti in rosso identificano i siti della città dell’Aquila e dei villaggi di Paganica ed Onna. Si nota come l’intero processo di rottura della faglia duri meno di 10 secondi. Questa durata non va confusa con quello dello scuotimento, che è molto maggiore per il propagarsi delle onde sismiche nella crosta, con riflessioni e rifrazioni multiple, come si vede nel video della propagazione delle onde in Italia centrale.
Questo tipo di indagini fornisce uno strumento essenziale per ottenere una descrizione della sorgente sismica che sia il più possibile rappresentativa dei reali processi sismogenetici, nell’ottica di migliorare la conoscenza dei meccanismi che sono alla base della generazione di un terremoto. Conoscenza indispensabile per la prevenzione e la mitigazione del rischio sismico.
