Archivi Blog

Sequenza sismica in Italia centrale: rapporto di sintesi sul terremoto del 30 ottobre

E’ stato pubblicato sul sito web dell’Istituto Nazionale di Geofisica e Vulcanologia (INGV), all’indirizzo www.ingv.it, il Rapporto di aggiornamento sul terremoto, di magnitudo 6.5, che ha colpito l’Italia centrale il 30 ottobre scorso (in italiano e in inglese). Nelle 49 pagine del rapporto sono descritti sia gli studi in corso sia i risultati preliminari basati sui dati dell’evento sismico in questione e sui rapporti tra questo e i precedenti terremoti del 24 agosto e del 26 ottobre.

map

Mappa della sismicità della regione dal 1985 ad oggi: in blu gli epicentri dal 1985 al luglio 2016, in giallo gli epicentri dal 24 agosto al 25 ottobre 2016, in rosso quelli dal 26 ottobre al 3 novembre 2016. Le stelle sono gli eventi di magnitudo uguale o maggiore di 5.0.

Nel rapporto vengono descritte: le analisi dei dati sismologici, con mappe e sezioni verticali attraverso la zona epicentrale; i modelli di faglia basati sui dati sismometrici e accelerometrici, sui dati geodetici (GPS e da SAR – interferometria radar da satellite) con le prime indicazioni della distribuzione del movimento di dettaglio delle varie faglie; l’impatto del terremoto sul territorio, “visto” dai dati sismici e stimato in base alle Shake maps  e alle analisi sul terreno; la fagliazione, osservata in superficie in tutta l’area interessata dai terremoti più forti dal 24 agosto al 30 ottobre; e, infine, una stima delle probabilità di accadimento delle future repliche (aftershocks).

gps-30ottobre

Mappa degli spostamenti co-sismici GPS – orizzontali (frecce rosse) e verticali (frecce blu) – ottenuti dalla combinazione di tre soluzioni geodetiche indipendenti per l’evento del 30 Ottobre 2016. I quadrati bianchi mostrano la posizione delle stazioni GPS permanenti e quelli color magenta delle stazioni RING-INGV (doi:10.13127/RING). I quadrati arancioni mostrano le stazioni GPS permanenti gestite da DPC e ISPRA. I quadrati verdi e blu mostrano le stazioni della rete CaGeoNet e della rete IGM, rispettivamente, ri-occupati dopo il 24 Agosto.

Dall’analisi di tutti i dati analizzati finora stanno emergendo risultati interessanti sul sistema di faglie che attraversa la regione e che si è attivato in questa sequenza sismica che, lo ricordiamo, è tuttora in corso. Sono infatti ancora diverse centinaia le repliche che vengono localizzate ogni giorno dalla Rete Sismica Nazionale dell’INGV.

L’analisi dei dati geologici, di quelli geodetici e sismologici sono coerenti nell’individuazione delle faglie e della loro cinematica, sebbene siano stati osservati alcuni elementi che sono in corso di approfondimenti e che sono il sintomo della complessità del sistema.

E’ comunque ormai chiaro che le faglie responsabili dei terremoti sono quelle note in letteratura come faglia del Monte Vettore-Bove, faglie normali (ossia a carattere estensionale) orientate in senso NNW-SSE e immergenti verso ovest. Le faglie si sono attivate per l’intero spessore dello strato sismogenetico, da circa 10 km di profondità fino alla superficie, producendo rotture ben visibili in affioramento, con “rigetti” (ossia scalini) che raggiungono in qualche punto i due metri.

faglia-1

Rottura cosismica del 30 Ottobre 2016: dislocazione lungo il piano di faglia in roccia.

Queste rotture osservate in superficie rappresentano la prosecuzione del movimento profondo sul piano di faglia, che si è originato intorno agli 8-10 km (l’ipocentro) e si è propagata lateralmente e verso l’alto. Sia i dati accelerometrici che quelli geodetici sono coerenti nell’individuazione delle zone di massimo spostamento della faglia del 30 ottobre, che si colloca tra le precedenti rotture del 24 agosto a sud e del 26 ottobre a nord, interessando principalmente il settore centrale del sistema di faglie e la sua parte più superficiale, dove vengono individuati spostamenti superiori ai 2.5 metri sul piano di faglia. È infatti evidente che lo spostamento lungo i diversi segmenti di faglia attivi in questi mesi non è avvenuto in maniera omogenea, ma ha avuto forti eterogeneità: spostamenti da pochi decimetri a 1-2 metri sul piano della stessa faglia. Questo potrebbe spiegare la complessità della sequenza, con l’attivazione successiva di segmenti di faglia di grandezza diversa e con spostamenti dei due lembi della faglia anche molto diversi. Sono visibili inoltre altre faglie “minori” che si sono mosse durante la sequenza.

sar-model3faglie

Modello preliminare della sorgente sismica del terremoto del 30 ottobre ottenuta da dati InSAR.

È in corso l’elaborazione di modelli più raffinati per identificare i dettagli di questa geometria e cinematica, confrontando e analizzando congiuntamente tutti i dati disponibili.

L’analisi dei dati accelerometrici del terremoto del 30 ottobre ad Amatrice, dove era stata installata una rete sismica temporanea per effettuare indagini propedeutiche alle attività di microzonazione sismica, ha evidenziato delle forti variazioni delle accelerazioni del suolo a distanze molto brevi, con amplificazioni fino a un fattore 5 rispetto a siti su roccia, dovute principalmente alla struttura geologica superficiale.

amatrice-emersito

Componente Z delle accelerazioni registrate nel paese di Amatrice – San Cipriano durante il terremoto Mw 6.5 del 30 ottobre 2016; per ogni stazione sono anche indicati i valori di PGA verticale.

Le analisi proseguono per seguire attentamente l’andamento delle repliche (il cui numero ha ormai superato quota 26000), per una mappatura di dettaglio degli effetti di superficie, per realizzare dei modelli di faglie che riescano a tener conto di tutti gli elementi osservati sul terreno e dal satellite.

Il report completo è disponibile al seguente link “RAPPORTO DI SINTESI SUL TERREMOTO IN CENTRO ITALIA MW 6.5 DEL 30 OTTOBRE 2016” a cura del Gruppo di Lavoro INGV sul Terremoto in centro Italia (2016).

Citare come: Gruppo di Lavoro INGV sul terremoto in centro Italia (2016). Rapporto di sintesi sul Terremoto in centro Italia Mw 6.5 del 30 ottobre 2016, doi: 10.5281/zenodo.166019


Licenza

Licenza Creative Commons
Quest’opera è distribuita con Licenza Creative Commons Attribuzione – Non opere derivate 4.0 Internazionale.

La sequenza sismica in Italia centrale a un mese dal suo inizio: un aggiornamento sugli studi in corso

Un mese dopo l’inizio della sequenza dell’Italia centrale, attivata dal terremoto di magnitudo 6.0 avvenuto nella notte tra il 23 e il 24 agosto, continuano gli studi dei sismologi e dei geologi per delineare sempre meglio il quadro di quanto è accaduto e fare degli scenari sulla possibile evoluzione.

rieti_23set_ore_12

La sequenza in Italia centrale dal 24 agosto al 23 settembre ore 12.

Anzitutto va detto che la sequenza è ancora in pieno svolgimento, pur con un numero minore di repliche (aftershocks) rispetto alle prime due settimane. Al momento la Rete Sismica Nazionale dell’INGV ha localizzato complessivamente circa 11500 repliche, in un’area che si estende per circa 40 chilometri in direzione NNO-SSE lungo la catena appenninica: 200 terremoti di magnitudo compresa tra 3 e 4, 14 quelli di magnitudo compresa tra 4 e 5 e uno di magnitudo maggiore di 5, oltre naturalmente a quello principale di magnitudo 6 del 24 agosto.

grafico_1mese_data

Andamento temporale della sequenza sismica dal 24 agosto al 23 settembre. In blu i primi dieci giorni, in giallo i dieci giorni successivi, in arancio gli ultimi dieci giorni.

In questo articolo forniamo un aggiornamento sintetico su quanto si è capito finora su questo evento sismico, segnalando che gli studi sono ancora in pieno svolgimento, con i nuovi dati registrati dalle reti di monitoraggio e con quelli che continuano a essere raccolti dalle squadre sul terreno, con le immagini dal satellite e le modellazioni al computer e altro ancora. Un rapporto più esaustivo con maggiori dettagli è disponibile QUI.

Tutti i dati acquisiti dalle reti sismiche e accelerometriche sono stati messi a disposizione della comunità scientifica per studi e ricerche. I dati registrati dalle stazioni della Rete Sismica Nazionale (RSN) e dalle stazioni temporanee installate dopo il terremoto sono disponibili nell’archivio EIDA (European Integrated Data Archive). Tutte le registrazioni accelerometriche della RSN e delle stazioni di emergenza per tutti i terremoti da magnitudo 3 in su avvenuti durante la recente sequenza (oltre 21.000 forme d’onda dal 24 agosto ad oggi) sono pubblicate a scaricabili dal portale web ISMD, ovvero dalla banca dati accelerometrica real time dell’INGV. Al pari, una selezione di forme d’onde processate manualmente, circa 200, è inoltre disponibile su ESM (Engineering Strong Motion Database). Per quanto riguarda le stazioni della Rete Accelerometrica Nazionale (RAN) del DPC, i dati originali sono disponibili sul sito web della RAN.

L’evento principale di ML 6.0 del 24 agosto è stato causato dallo scorrimento di una faglia distensiva (o faglia normale) orientata in direzione NNO-SSE (azimut ~156°) inclinata verso SO con una pendenza di circa 50°. La lunghezza della faglia che si è attivata con questo terremoto è di 20-25 km. A partire dal punto di nucleazione del terremoto, localizzato a circa 8 km di profondità in prossimità di Accumoli, la rottura della faglia è stata bilaterale (propagandosi sia verso NO sia verso SE, ossia verso Norcia e verso Amatrice). Sia i dati accelerometrici che quelli SAR indicano che la faglia non si è dislocata in maniera omogenea lungo la sua estensione ma è caratterizzata da due zone principali di concentrazione dello spostamento sul piano di faglia stesso, con valori massimi di circa un metro.  La durata della rottura, ricavata dai dati accelerometrici, è stata di circa 6 secondi.

Il risultato della modellazione mostra due aree di rottura principali sul piano di faglia uno sopra l’ipocentro ed uno a NW

Il piano di faglia del terremoto del 24 agosto in una vista frontale. Il risultato della modellazione mostra due aree di rottura principali sul piano di faglia, uno sopra l’ipocentro ed uno verso nordovest. La stella rossa indica la posizione del punto di origine della rottura sulla faglia (l’ipocentro), le frecce lo spostamento sulla faglia.

Le aree della faglia caratterizzate dai maggiori valori di slip (movimento lungo il piano di faglia) sono quelle meno popolate dagli aftershocks della sequenza sismica (in particolare questo è vero per l’area con una maggiore concentrazione di slip nel parte più meridionale, vedi figura sotto, pannello a sinistra). Si nota anche come gli aftershocks di magnitudo più elevata (ML>3.0) tendano a concentrarsi lungo i bordi delle aree di maggior slip (figura sotto, pannello a destra).

. Confronto tra la distribuzione degli aftershock e la distribuzione di slip sul piano di faglia principale. Il modello di sorgente riportato in figura è quello derivante dall’inversione dei dati accelerometrici (doi 105281/zenodo.61460) le localizzazioni degli aftershcoks sono mostrate a sinistra tutte, a destra solo quelle con magnitudo superiore a 3.

Confronto tra la distribuzione degli aftershocks e la distribuzione di slip sul piano di faglia principale (proiezione sulla superficie). Il modello di sorgente riportato in figura è quello derivante dall’inversione dei dati accelerometrici (doi 105281/zenodo.61460). Nel pannello di sinistra sono riportate le localizzazioni di tutti gli aftershocks – indipendentemente dalla magnitudo – mentre nel pannello di destra sono graficate solo quelle con magnitudo superiore a 3.

I dati SAR e GPS hanno evidenziato le deformazioni superficiali dell’area dovute al terremoto del 24 agosto, come già descritto in altri articoli del blog. I dati di altri satelliti analizzati in questi giorni hanno permesso di modellare le deformazioni visibili dal SAR con una faglia a slip eterogeneo o con due faglie adiacenti con geometria leggermente differente e analoga distribuzione di slip.

Rappresentazione 3D del modello di slip a due faglie. La linea arancione rappresenta l'intersezione della faglia Nord con la superficie; la linea viola è l'intersezione della faglia Sud. Il simbolo rosso è il mainshock, i simboli verdi sono gli altri eventi maggiori fino al 29/8.

Rappresentazione 3D del modello di slip a due faglie ricavato dai dati SAR. La linea arancione rappresenta l’intersezione della faglia nord con la superficie; la linea viola è l’intersezione della faglia sud. Il simbolo rosso è l’ipocentro dell’evento principale, i simboli verdi sono gli altri eventi maggiori fino al 29 agosto.

Le deformazioni del terreno visibili dai dati SAR analizzati finora sembrano suggerire che il movimento sulla faglia non sia arrivato a interessare direttamente la superficie. Gli spostamenti del terreno osservati dal satellite sembrano piuttosto legati a due fenomeni deformativi a una scala diversa: uno “profondo” legato al movimento sulla faglia a profondità tra i 10 e i 5 km (le aree colorate in blu e celeste nella figura sopra), e un altro, più circoscritto, che produce delle frange caratteristiche negli interferogrammi e che sarebbe legato a deformazioni più superficiali, forse attribuibili a fenomeni gravitativi. Sono in corso verifiche e confronti tra questi dati e quelli di terreno (v. sotto) per discriminare la natura delle rotture evidenziate sul terreno, se queste siano cioè diretta espressione in superficie della faglia responsabile del terremoto del 24 agosto, o, alternativamente, siano degli elementi deformativi secondari.

Come noto, molte “rotture superficiali” del terreno sono state rilevate e mappate lungo la faglia del monte Vettore: si tratta di una zona di fratturazione superficiale continua che si estende per circa 5.2 km lungo il fianco sud-occidentale del monte. Le squadre dei rilevatori stanno tuttora investigando l’area epicentrale per un’estensione totale di circa 40 km tra Castelluccio di Norcia, a Nord, e la Località Ortolano posta a sud del lago artificiale di Campotosto. Al momento, sono state catalogate informazioni geologiche su oltre 3000 punti di osservazione. In generale, nell’intera area investigata sono state segnalate numerose fratture lungo i versanti montuosi e i campi coltivati (questi elementi deformativi risultano spesso ben visibili sul manto stradale), insieme a frane, scoscendimenti e crolli di massi, di piccole-medie dimensioni. Le caratteristiche geometriche degli elementi deformativi osservati sono descritte sinteticamente nella figura sotto.

diagramma a rosa delle direzioni delle rotture cosismiche sui depositi non consolidati (in rosso) e lungo il piano di faglia in roccia del Monte Vettore e del Monte Vettoretto

Diagramma a rosa delle direzioni delle rotture cosismiche sui depositi non consolidati (in rosso) e lungo il piano di faglia in roccia del Monte Vettore e del Monte Vettoretto.

 

Rilievi macrosismici

A seguito del terremoto del 24 agosto, seguito poco dopo dalla forte replica di magnitudo Ml 5.4 (Mw 5.3), le squadre del DPC, CNR‐IGAG e INGV sono partite per effettuare il rilievo macrosismico dell’area epicentrale. Le diverse squadre hanno intrapreso il rilievo nell’area di Norcia, Arquata e Antrodoco, frazione per frazione, talvolta anche nelle stesse località ed in tempi diversi, al fine di confrontare e tarare le osservazioni macrosismiche. Le squadre si sono quindi incontrate più volte sul campo per scambiarsi dati, impressioni e pareri. In remoto hanno operato altri colleghi che raccoglievano notizie e segnalazioni, comunicandole alle squadre.

È opportuno ricordare che il rilievo macrosismico in scala MCS in fase di emergenza è di carattere speditivo, e ha l’obiettivo precipuo di definire al meglio la distribuzione del danneggiamento medio‐grave nell’immediato (2 giorni) e ‐ più in generale ‐ di fornire i limiti dell’area di danno in un tempo ragionevolmente breve (15 giorni) e con un criterio omogeneo e di immediata applicazione, quale quello fornito dalla scala MCS. In particolare, il rilievo viene eseguito utilizzando le percentuali di danneggiamento della scala MCS (1930) quantificate da Molin (2009) in progressione crescente per i gradi di intensità ≥ V MCS e secondo i cinque livelli di danno previsti dalla scala originaria.

D’altro canto, l’applicazione della scala macrosismica europea EMS (Grünthal, 1998), presentando difficoltà collegate alla suddivisione degli edifici di un centro abitato secondo le classi di vulnerabilità proposte dalla scala stessa, è inapplicabile in fase di emergenza. Anche in occasione del terremoto dell’Aquila nel 2009, non è risultato possibile assegnare con sicurezza e con la rapidità generalmente richiesta dai rilievi macrosismici in emergenza la classe di vulnerabilità agli edifici, se non a seguito di sopralluoghi prolungati nei mesi successivi (p.e., Molin et al., 2010).

Dal presente rapporto sono al momento escluse le informazioni di effetti di risentimento che non siano stati direttamente verificati dalle squadre. Alle ore 20 del 15 Settembre sono state rilevate direttamente un totale di 283 località distribuite in 76 comuni. Nella figura sotto sono mostrate le intensità attribuite.

Figura 7. Distribuzione delle intensità rilevate per località in termini di scala macrosismica MCS (cerchi bianchi proporzionali alla intensità di sito). Isosisme dal VI al IX MCS a tratteggio nero. Stelle rosse, eventi con Ml>4 (INGV). Le campiture colorate indicano approssimativamente la distribuzione areale dell'intensità. Sullo sfondo DTM a 20 metri. Le linee blu sono le faglie del Vettore a nord e della Laga a sud. In rosso la rottura superficiale del 2016

Distribuzione delle intensità rilevate per località in termini di scala macrosismica MCS (cerchi bianchi proporzionali alla intensità di sito). Isosisme dal VI al IX MCS a tratteggio nero. Stelle rosse, eventi con Ml>4 (INGV). Le campiture colorate indicano approssimativamente la distribuzione areale dell’intensità. Sullo sfondo DTM a 20 metri. Le linee blu sono le faglie del Vettore a nord e della Laga a sud. In rosso la rottura superficiale del 2016

Queste qui riportate sono solo alcune delle informazioni che si trovano nel “SECONDO RAPPORTO DI SINTESI SUL TERREMOTO DI AMATRICE ML 6.0 DEL 24 AGOSTO 2016 (ITALIA CENTRALE)” a cura del Gruppo di Lavoro INGV sul terremoto di Amatrice (2016).

Citare come: Gruppo di Lavoro INGV sul terremoto di Amatrice (2016). Secondo rapporto di sintesi sul Terremoto di Amatrice Ml 6.0 del 24 Agosto 2016 (Italia Centrale), doi: 10.5281/zenodo.154400


Licenza

Licenza Creative Commons
Quest’opera è distribuita con Licenza Creative Commons Attribuzione – Non opere derivate 4.0 Internazionale.

Terremoto in Italia centrale: Modellazione della sorgente sismica e trasferimento di stress sulle faglie limitrofe

Secondo articolo del Gruppo di Lavoro INGV/IREA-CNR, (Istituto Nazionale di Geofisica e Vulcanologia (INGV) e Consiglio Nazionale delle Ricerche – Istituto per il Rilevamento Elettromagnetico dell’Ambiente, CNR-IREA di Napoli) sull’analisi di dati satellitari volta allo studio delle sorgenti sismiche.


Partendo dalle mappe di spostamento della superficie terrestre è possibile effettuare una modellazione della sorgente sismica del terremoto.

Ma cosa vuol dire “modellazione di una sorgente sismica”? La modellazione consiste nel riprodurre il campo di spostamento misurato dal satellite con un campo di spostamento teorico, facendo variare, di volta in volta, i parametri che caratterizzano la sorgente sismica (faglia sismogenetica), sino ad ottenere una sorgente sismica teorica il cui campo di spostamento in superficie risulta essere molto simile a quello osservato dal satellite. Tale approccio alla modellazione è chiamato “inversione del dato”, proprio perché partendo dal dato osservato si cerca di identificare e caratterizzare la sorgente della deformazione.

L’obiettivo è quello di individuare, tra le faglie note o ancora non note, quella che probabilmente ha generato il terremoto e capire la distribuzione del movimento (slip) lungo il piano della faglia stessa. Naturalmente, individuata la faglia che ha generato il terremoto, geologi e sismologi confrontano i dati di tale modellazione con i dati ottenuti dai rilievi di campagna e le faglie già note in letteratura al fine di caratterizzare i movimenti attuali e possibilmente anche la storia sismica della faglia stessa. Oltre alla geometria del piano di faglia è importante simulare (e quindi ipotizzare il più realisticamente possibile) il movimento della faglia stessa andando a determinare quali porzioni di faglia si sono attivate durante il terremoto e di quanto si sono mosse.

La geologia ci insegna che le faglie non sono dei piani regolari che tagliano la crosta terrestre, ma piuttosto delle superfici articolate che mostrano anche grandi variazione della loro orientazione nello spazio e nel loro angolo di inclinazione rispetto alla verticale. Tuttavia, per poter risolvere in tempi brevi (quelle di un’emergenza sismica) un problema numerico, ovvero quello di riprodurre, con simulazioni al computer, il campo di spostamento prodotto dal movimento di una faglia, c’è bisogno introdurre delle semplificazioni. In particolare, nei modelli presentati in seguito, le faglie vengono semplificate come dei piani e la crosta terrestre come un mezzo perfettamente elastico, tralasciando complicazioni geometriche e comportamenti anelastici della crosta.

Risultati della modellazione della sorgente sismica:

I dati utilizzati nel procedimento di inversione sono: un totale di circa 20 mila valori di spostamento del suolo in linea di vista ricavati dallo srotolamento di 5 interferogrammi (due ALOS, due Sentinel e uno COSMO-SkyMed) e 107 valori di spostamento provenienti da stazioni CGPS (Gruppo di Lavoro INGV-CNT Centro Analisi Dati GPS, 2016).

Il processo di inversione segue due step principali. Il primo è quello di invertire i dati SAR e GPS per ottenere i parametri geometrici della faglia stessa (lunghezza, larghezza, profondità, angolo di inclinazione e orientazione nello spazio), ipotizzando uno scorrimento costante e uniforme su tutto il piano di faglia. Il secondo step consiste nel simulare la distribuzione del movimento (slip) nelle singole porzioni di faglia, tenendo fissi i parametri geometrici ottenuti con il primo step.

Ipotizzando l’attivazione di un singolo piano di faglia, con il primo e il secondo step è stato ottenuto il modello di faglia mostrato in Figura 1.

Figura 1: Distribuzione del movimento (slip) lungo il piano della faglia per il modello a faglia singola. Sono mostrati anche i due eventi maggiori della sequenza al 25 agosto (pentagoni rossi) e le repliche rilocalizzate fino al 29 agosto (puntini neri). Localizzazione delle repliche fornita da INGV AMA_LOC Working Group – DOI: 10.5281/zenodo.61404

Per simulare le possibili variazioni nello spazio (inclinazione e orientazione rispetto al nord) del piano di faglia o la non completa continuità di questo in tutta la sua estensione, è stata effettuata una modellazione con due piani di faglia contigui, caratterizzati da geometrie leggermente differenti l’uno dall’altro. Il risultato è mostrato in Figura 2.

Figura 2: Distribuzione del movimento (slip) lungo il piano della faglia per il modello a due faglie. Sono mostrati i due eventi maggiori della sequenza al 25 agosto (pentagoni rossi) e le repliche rilocalizzate fino al 29 agosto (puntini neri). Localizzazione delle repliche fornita da INGV AMA_LOC Working Group – DOI: 10.5281/zenodo.61404

Per avere una visione più chiara del piano di faglia ottenuto dall’inversione dei dati di deformazione, nelle figure 3 e 4, viene mostrata una rappresentazione 3D dei due modelli a singola e doppia faglia.

Figura 8: Rappresentazione 3D del modello di slip a una faglia. La linea gialla rappresenta l'intersezione del piano del modello con la superficie. Le linee rosse sono tracce di faglie dalla bibliografia geologica. Il simbolo rosso è l’evento principale, i simboli verdi sono gli altri eventi maggiori fino al 29/8.

Figura 3: Rappresentazione 3D del modello di slip a una faglia. La linea gialla rappresenta l’intersezione del piano del modello con la superficie. Le linee rosse sono tracce di faglie dalla bibliografia geologica. Il simbolo rosso è l’evento principale, i simboli verdi sono gli altri eventi maggiori fino al 29 agosto. Ogni quadratino colorato ha il lato di 1 km.

 

Figura 9: Rappresentazione 3D del modello di slip a due faglie. La linea arancione rappresenta l'intersezione della faglia Nord con la superficie; la linea viola è l'intersezione della faglia Sud. Il simbolo rosso è l’evento principale, i simboli verdi sono gli altri eventi maggiori fino al 29/8.

Figura 4: Rappresentazione 3D del modello di slip a due faglie. La linea arancione rappresenta l’intersezione della faglia nord con la superficie; la linea viola è l’intersezione della faglia sud. Il simbolo rosso è l’evento principale, i simboli verdi sono gli altri eventi maggiori fino al 29 agosto. Ogni quadratino colorato ha il lato di 1 km.

I parametri dei due modelli sono riportati nella tabella sottostante:

Modello Lunghezza Larghezza Profondità del tetto della faglia Strike Dip Rake Slip max
Singola faglia ~21 km ~9 km ~1500 m 164° 46° -73° 120 cm
Doppia faglia -Nord ~8 km ~8 km ~3000 m 175° 39° -65° 140 cm
Doppia faglia -Sud ~12 km ~5 km ~2500 m 165° 51° -70° 130 cm

Come spiegato precedentemente, l’utilizzo di un modello con una o due faglie serve a riprodurre (in maniera semplificata) eventuali variazioni nello spazio dei parametri geometrici della faglia e non necessariamente sta a significare l’attivazione di due piani di faglia distinti. I due modelli sono sostanzialmente equivalenti, se confrontati con il dato di spostamento misurato in superficie.

E’ possibile tentare di attribuire la rottura a piani di faglia visibili in superficie, prolungando i piani dei modelli fino ad intersecare la topografia. La traccia superficiale del modello a singola faglia corre parallela e quasi coincidente (entro ±800 m) con la traccia del sistema di faglie Gorzano-Laga-Vettore (Figura 3). Nel modello a due faglie si nota un andamento identico per la faglia sud, mentre la traccia della faglia nord ruota verso NE e dovrebbe emergere circa 3 km ad est del Monte Vettore.

Entrambi i modelli sono compatibili con gli andamenti delle repliche, nei limiti delle incertezze, con un leggero miglior accordo nella parte nord del modello a due faglie rispetto a quello con faglia singola.

La tabella sopra riporta la dimensione della rottura ricavabile dai valori di slip per i due modelli. In entrambi i modelli, andando verso sud dall’ipocentro, la rottura si ferma poco prima di Amatrice, mentre verso nord i valori di slip del modello a doppia faglia vanno a zero circa a metà della Piana di Castelluccio (o poco più a nord per la faglia singola). Le due maggiori concentrazioni di slip sono localizzate a profondità intorno a 6 km nella parte nord e 4 km a sud. A seconda dei modelli, i valori di slip sono circa uguali a zero a profondità minori di 1.5-3 km (a nord) o 700 m (a sud). Nelle parti più superficiali delle faglie vengono modellati in maniera discontinua valori di slip intorno a pochi cm.

Calcolo del trasferimento di stress alla faglie limitrofe

Partendo dal movimento calcolato sui piani di faglia attivatisi durante il terremoto del 24 agosto 2016, è possibile calcolare la variazione di stress statico, noto come Coulomb Failure Function (CFF) sulle faglie note ubicate nell’intorno di quella del 24 agosto.

E’ importante sottolineare che il calcolo delle CFF può dare informazioni utili e interessanti per comprendere quali faglie, nell’intorno dell’ipocentro, possono aver subito delle variazioni di stress e quindi capire meglio l’evoluzione della sequenza. È evidente però che tale calcolo non ha nulla a che vedere con la previsione dei terremoti o con l’individuazione di aree più pericolose rispetto ad altre. Infatti il calcolo delle CFF viene effettuato sulla base di numerose assunzioni circa le caratteristiche reologiche della crosta, i meccanismi di fagliazione e le geometrie dei piani stessi. Assunzioni che possono essere anche molto diverse dalla situazione reale che non possiamo conoscere, ma solo immaginare. Inoltre, lo stesso risultato finale del calcolo delle CFF fornisce la sola variazione dei valori di stress; al contrario non è possibile conoscere il valore di partenza dello stress cui il singolo piano di faglia è sottoposto. Risulta perciò evidente che tale calcolo non può fornirci una previsione sull’evoluzione della sequenza sismica, ma il suo valore è puramente scientifico, come ulteriore elemento utile agli scienziati per studiare i terremoti.

Figura 10: Distribuzione delle differenze della Coulomb Failure Function (CFF), causate dalla dislocazione del modello a due faglie sulle strutture circostanti. Piani di faglia delle strutture vicine forniti dal Gruppo Emergeo.

Figura 5: Distribuzione delle differenze della Coulomb Failure Function (CFF), causate dalla dislocazione del modello a due faglie sulle strutture circostanti. Piani di faglia delle strutture vicine forniti dal Gruppo Emergeo.

Nella simulazione è stato utilizzato il modello a due faglie. La Figura 5 mostra il risultato di tale calcolo sui piani delle sei faglie più vicine (EMERGEO Working Group 2016; Terremoto di Amatrice del 24 agosto 2016: Effetti Cosismici, doi: 10.5281/zenodo.61566). Valori positivi considerevoli di variazione di CFF, fino a 0.5-0.6 MPa si riscontrano al limite NW della faglia del Monte Gorzano, interessata dopo l’evento principale da una intensa attività di bassa magnitudo su una struttura esterna ai modelli presentati. A nord i valori sulla faglia del Vettore-Bove sembrano molto più bassi, tuttavia si evidenzia come sia possibile che il piano utilizzato per il calcolo non sia del tutto rappresentativo, in quanto non è in continuità con la faglia del Monte Vettore. Se il piano proseguisse a SE si avrebbero concentrazioni di CFF probabilmente analoghe a quelle calcolate per la faglia del Gorzano.

a cura del Gruppo di Lavoro INGV/IREA-CNRGiuseppe Pezzo, Christian Bignami, Cristiano Tolomei, Simone Atzori, Andrea Antonioli, Salvatore Stramondo, Stefano Salvi (INGV-CNT); Manuela Bonano, Raffaele Castaldo, Francesco Casu, Claudio De Luca, Vincenzo De Novellis, Riccardo Lanari, Mariarosaria Manzo, Michele Manunta, Antonio Pepe, Susi Pepe, Pietro Tizzani, Ivana Zinno (IREA-CNR)


Articolo estratto da: Gruppo di lavoro IREA-CNR & INGV, 2016, Sequenza sismica di Amatrice: aggiornamento delle analisi interferometriche satellitari e modelli di sorgente, DOI: 10.5281/zenodo.61682


Licenza

Licenza Creative Commons
Quest’opera è distribuita con Licenza Creative Commons Attribuzione – Non opere derivate 4.0 Internazionale.

Terremoto in Italia centrale: Misure dello spostamento del suolo tramite interferometria satellitare

Sin dalle primissime ore dopo il terremoto, il Gruppo di Lavoro INGV/IREA-CNR, costituito dall’Istituto Nazionale di Geofisica e Vulcanologia (INGV) e dal Consiglio Nazionale delle Ricerche (Istituto per il Rilevamento Elettromagnetico dell’Ambiente, CNR-IREA di Napoli) come centri di competenza del Dipartimento della Protezione Civile nei settori della sismologia e dell’elaborazione dei dati radar satellitari,  si è attivato per un’analisi di dati satellitari volta alla misura dei movimenti permanenti del suolo prodotti dal terremoto ed allo studio delle sorgenti sismiche.

Verranno quindi pubblicati due articoli che descrivono i risultati ottenuti dal Gruppo.


L’utilizzo di immagini satellitari radar ad apertura sintetica (Syntetic Aperture Radar, SAR) permette di misurare la deformazione permanente del suolo causata dal terremoto del 24 agosto 2016.

La tecnica utilizzata per effettuare tali tipi di misure è quella nota interferometria differenziale DInSAR, e già utilizzata per effettuare le misure di spostamenti in molti altri terremoti, come quello dell’Emilia del 2012, della Lunigiana del 2013 e dell’Aquila del 2009.

La tecnica

Utilizzando una coppia di immagini radar acquisite dal satellite dalla medesima posizione, prima e dopo un evento sismico è possibile produrre un interferogramma (Figura 1a), cioè una mappa di spostamento espressa in termini di differenze di fase, tra le due immagini radar pre e post evento sismico, del segnale elettromagnetico emesso dal satellite, “riflesso” dalla superficie terrestre e catturato dal sensore (rappresentato con diversi cicli di colore). Ogni ciclo, o frangia, rappresenta uno spostamento della superficie terrestre lungo la linea di vista del satellite di mezza lunghezza d’onda del segnale emesso.

Figura 1: esempio di interferogramma (A) ottenuto dall'elaborazione di due immagini SAR del satellite Sentinel-1. Nel pannello B si riporta la corrispondente mappa delle spostamento del suolo dovuto al terremoto del 24 agosto 2016. Le immagini utilizzate sono relative alle date 21 e 27 agosto 2016.

Figura 1: esempio di interferogramma (A) ottenuto dall’elaborazione di due immagini SAR del satellite Sentinel-1. Nel pannello B si riporta la corrispondente mappa delle spostamento del suolo dovuto al terremoto del 24 agosto 2016. Le immagini utilizzate sono relative alle date 21 e 27 agosto 2016.

Rappresentando ogni frangia un incremento noto di spostamento rispetto la frangia precedente è possibile misurare lo spostamento totale dell’area partendo dalle zone in cui non sono state misurate deformazioni (zone dove non compaiono frange intefirferometriche) sino ad arrivare alla zona più deformata (parte centrale dell’ellissi disegnata dalle frange stesse). Questa operazione è detta “srotolamento” delle frange interferometriche e permette di ottenere una mappa dello spostamento del suolo avvenuto durante il sisma (Figura 1b).

I dati disponibili

Durante l’emergenza sismica iniziata con il terremoto del 24 agosto, l’analisi interferometrica satellitare ha beneficiato di un gran numero d’immagini SAR acquisite, a cavallo dell’evento sismico, dal satellite giapponese ALOS-2 (Japan Aerospace Esploration Agency, JAXA), operante in banda L (lunghezza d’onda 23.6 cm), dai sensori in banda C (lunghezza d’onda 5.6 cm) della costellazione Sentinel-1 del Programma Europeo Copernicus (Agenzia Spaziale Eurpoea, ESA) e da quelli in banda X (lunghezza d’onda 3.1 cm) della costellazione italiana COSMO-SkyMed, sviluppata dall’Agenzia Spaziale Italiana (ASI) in cooperazione con il Ministero della Difesa.

I risultati

In Figura 2 sono riportate alcune delle mappe di deformazione co-sismica in linea di vista generate a partire dalle coppie disponibili sfruttando la tecnica dell’Interferometria SAR Differenziale (DInSAR).

Le mappe di spostamento sono espresse in cm e calcolate lungo la linea di vista del satellite (LoS), inclinata di circa 35° rispetto alla verticale.

Figura 2: Mappe di spostamento del suolo co-sismico (nelle linee di vista del radar) ottenute a partire dai dati radar ALOS 2 e Sentinel-1 acquisiti a cavallo dell’evento sismico del 24 agosto da orbite sia ascendenti sia discendenti. Le mappe di spostamento sono espresse in cm e calcolate lungo la linea di vista del satellite (LoS), inclinata di circa 35° rispetto alla verticale.

Figura 2: Mappe di spostamento del suolo co-sismico (nelle linee di vista del radar) ottenute a partire dai dati radar ALOS 2 e Sentinel-1 acquisiti a cavallo dell’evento sismico del 24 agosto da orbite sia ascendenti sia discendenti. Le mappe di spostamento sono espresse in cm e calcolate lungo la linea di vista del satellite (LoS), inclinata di circa 35° rispetto alla verticale.

Relativamente ai dati Sentinel-1, grazie alla disponibilità dell’ESA, si è avuta la possibilità di ottenere i dati del satellite Sentinel-1B nonostante lo stesso sia ancora in fase di commissioning e quindi non pienamente operativo. Ciò ha permesso di sfruttare in maniera congiunta i dati acquisiti dai due satelliti e di generare coppie interferometriche “cross-sensor” (cioè con un’immagine acquisita da Sentinel-1A e una acquisita da Sentinel-1B) con un intervallo temporale fra le acquisizioni pre- e post- evento di soli 6 giorni (rispetto ai 12 possibili con l’utilizzo di un solo sensore), riducendo i tempi di risposta.

Inoltre, grazie all’utilizzo congiunto delle immagini SAR acquisite lungo orbite ascendenti (da sud a nord) e discendenti (da nord a sud), è stato possibile ricavare le componenti verticale ed est-ovest degli spostamenti del suolo (Figura 3).

Figura 3: Mappe delle componenti verticale (sinistra) ed est-ovest (destra) dello spostamento del suolo, ottenute sfruttando congiuntamente i passaggi ascendenti (sud-nord) e discendenti (nord-sud) dei dati radar ALOS 2 e Sentinel-1 acquisiti a cavallo dell’evento sismico. Nella mappa relativa alla componente est-ovest (destra) sono indicate due aree in cui sono stati rilevati degli spostamenti localizzati e di cui si fornisce maggiore dettaglio in Figura 4.

Figura 3: Mappe delle componenti verticale (sinistra) ed est-ovest (destra) dello spostamento del suolo, ottenute sfruttando congiuntamente i passaggi ascendenti (sud-nord) e discendenti (nord-sud) dei dati radar ALOS 2 e Sentinel-1 acquisiti a cavallo dell’evento sismico. Nella mappa relativa alla componente est-ovest (destra) sono indicate due aree in cui sono stati rilevati degli spostamenti localizzati e di cui si fornisce maggiore dettaglio in Figura 4.

I risultati ottenuti mostrano un abbassamento del suolo dalla caratteristica forma a “cucchiaio” che si estende per circa 20 km in direzione NNW. Il valore massimo dello spostamento del suolo verticale è di circa 20 cm in corrispondenza dell’area di Accumoli. Si nota inoltre che la zona deformata visibile sulla componente est-ovest interessa un’area più estesa rispetto a quella verticale (circa 20 x 25 km2) ed è caratterizzata dalla presenza di quattro aree alternate di spostamento, con valori massimi di deformazione di circa 16 cm verso ovest.

Figura 4: Dettaglio della mappa della componente est-ovest (Figura 2) per due aree interessate da spostamenti localizzati. A sinistra si evidenzia uno spostamento lungo est localizzato su un versante ad ovest dell’abitato di Castelluccio. A destra, invece è indicato il movimento verso ovest del versante ovest del Monte Vettore (si veda anche l’analisi dei dati del satellite COSMO-SkyMed in Figura 5).

Figura 4: Dettaglio della mappa della componente est-ovest (Figura 2) per due aree interessate da spostamenti localizzati. A sinistra si evidenzia uno spostamento lungo est localizzato su un versante ad ovest dell’abitato di Castelluccio. A destra, invece è indicato il movimento verso ovest del versante ovest del Monte Vettore (si veda anche l’analisi dei dati del satellite COSMO-SkyMed in Figura 5).

Dall’analisi delle mappe di deformazione lungo la linea di vista del satellite e delle componenti verticale ed est-ovest dello spostamento, è possibile ricavare informazioni più di dettaglio su alcuni effetti localizzati, probabilmente legati a fenomeni di instabilità di versante (Figura 4). È infatti possibile che lo scuotimento del suolo dovuto al terremoto abbia innescato il movimento di alcune porzioni di versante interessate da instabilità. A seconda dell’entità del movimento e dell’estensione dell’area interessata, il movimento può essere osservato e misurato attraverso l’interferometria.

Figura 5: Mappa della deformazione co-sismica (in alto) ottenuta a partire dai dati radar satellitari acquisiti dalla costellazione COSMO-SkyMed, il 20 agosto 2016 (immagine pre-evento) e il 28 agosto 2016 (immagine post-evento). La freccia indica una deformazione localizzata in corrispondenza del fianco del Monte Vettore, probabilmente legata a un fenomeno d’instabilità di versante. In basso è riportata la mappa di deformazione ottenuta da dati CSK acquisiti precedentemente al sisma (coppia interferometrica 03072016CSK-20082016CSK), dove non si notano segnali deformativi significativi.

Figura 5: Mappa della deformazione co-sismica (in alto) ottenuta a partire dai dati radar satellitari acquisiti dalla costellazione COSMO-SkyMed, il 20 agosto 2016 (immagine pre-evento) e il 28 agosto 2016 (immagine post-evento). La freccia indica una deformazione localizzata in corrispondenza del fianco del Monte Vettore, probabilmente legata a un fenomeno d’instabilità di versante. In basso è riportata la mappa di deformazione ottenuta da dati CSK acquisiti precedentemente al sisma (coppia interferometrica 03072016CSK-20082016CSK), dove non si notano segnali deformativi significativi.

Tali fenomeni sono, inoltre, visibili anche nella mappa di deformazione co-sismica relativa ai dati SAR COSMO-SkyMed (CSK) che, grazie alle elevate risoluzioni spaziali del sistema, permette di stimare gli spostamenti del suolo con grande dettaglio spaziale. A tal riguardo, nell’immagine in alto della Figura 5 è riportata la mappa di deformazione co-sismica CSK in LOS, relativa alla zona che si estende tra Tufo e Pescara del Tronto fino all’area di Castelluccio. Si noti come la migliore risoluzione spaziale consenta di analizzare nel dettaglio anche lo spostamento del suolo di un’area relativamente piccola che si estende per circa 800 m x 600 m sul fianco del Monte Vettore, verosimilmente legato ad un fenomeno di instabilità di versante. Quest’ultimo è probabilmente legato a frane sismo-indotte o a faglie riattivate, visto che l’analisi dei dati CSK acquisiti nel periodo precedente all’evento sismico (immagine in basso della Figura 5) non mostra segnali deformativi significativi.

a cura del Gruppo di Lavoro INGV/IREA-CNRGiuseppe Pezzo, Christian Bignami, Cristiano Tolomei, Simone Atzori, Andrea Antonioli, Salvatore Stramondo, Stefano Salvi (INGV-CNT); Manuela Bonano, Raffaele Castaldo, Francesco Casu, Claudio De Luca, Vincenzo De Novellis, Riccardo Lanari, Mariarosaria Manzo, Michele Manunta, Antonio Pepe, Susi Pepe, Pietro Tizzani, Ivana Zinno (IREA-CNR)


Articolo estratto da: Gruppo di lavoro IREA-CNR & INGV, 2016, Sequenza sismica di Amatrice: aggiornamento delle analisi interferometriche satellitari e modelli di sorgente, DOI: 10.5281/zenodo.61682


Licenza

Licenza Creative Commons
Quest’opera è distribuita con Licenza Creative Commons Attribuzione – Non opere derivate 4.0 Internazionale.

Terremoto in Italia centrale: Effetti in superficie rilevati dal gruppo EMERGEO-INGV

Il terremoto è l’effetto dell’improvviso scorrimento relativo delle masse rocciose che compongono la crosta terrestre lungo zone di rottura dette faglie. I grandi terremoti, con magnitudo generalmente maggiore di 6, producono effetti non solo sul costruito, ma anche deformazioni permanenti sulla superficie terrestre in un’area ampia anche alcune centinaia di km2. Tali deformazioni sono la diretta conseguenza dello spostamento relativo delle porzioni di crosta terrestre lungo il piano di faglia (effetti diretti o primari), ma anche dovuti al passaggio dell’onda sismica (effetti indiretti o secondari).

Schema semplificato del movimento sul piano di faglia in occasione di un terremoto con meccanismo estensionale (in risposta ad una estensione della crosta – vedi frecce arancioni). La deformazione che ha luogo sul piano di faglia alle profondità ipocentrali si propaga fino alla superficie e forma una scarpata di faglia che separa, nel caso di una faglia normale, un settore sollevato [+] da uno ribassato [-].

Schema semplificato del movimento sul piano di faglia in occasione di un terremoto con meccanismo estensionale (in risposta ad una estensione della crosta rappresentata dalle frecce arancioni). La deformazione che ha luogo sul piano di faglia alle profondità ipocentrali si propaga fino alla superficie e forma una scarpata di faglia che separa, nel caso di una faglia normale, un settore sollevato da uno ribassato.

Gli effetti diretti sono quelli che hanno luogo in corrispondenza dell’intersezione del piano di faglia, che si è mosso, con la superficie terrestre. Qui infatti, lo spostamento avvenuto in profondità si evidenzia con la formazione di fratture allineate con il piano di faglia e che spesso formano un vero e proprio scalino (scarpata di faglia). Questi elementi rispecchiano il movimento sulla faglia in profondità sia come geometria, entità e direzione dello spostamento, e essi generalmente sono associati a faglie geologiche segnate nelle carte tematiche dedicate come faglie attive.

Altri effetti diretti del terremoto sono il ribassamento/sollevamento di porzioni della superficie terrestre (frecce gialle nella figura sopra) che possono essere osservate con metodi di studio sia terrestri che satellitari (geodesia e telerilevamento). Quando queste deformazioni sono di entità ingente possono produrre effetto superficiali visibili anche a occhio nudo come ad esempio fenomeni di impaludamento, inondazione, o emersione di porzioni importanti del territorio.

Tra gli effetti indiretti, quindi non necessariamente legati al movimento sul piano di faglia ma solo come risposta allo scuotimento prodotto dal terremoto, la fratturazione senza organizzazione geometrica regolare, la liquefazione e le frane sono i più comuni.

La liquefazione ha luogo in presenza di depositi sciolti, a prevalenza sabbiosa e saturi in acqua, quindi la osserveremo preferibilmente in piane alluvionali e costiere. Il passaggio dell’onda sismica crea delle sovrappressioni nei livelli saturi che porta alla loro liquefazione (perdita di coesione del terreno fino ad assumere un comportamento tipico dei fluidi) e alla conseguente migrazione del materiale liquefatto verso la superficie. Le frane vengono innescate dalle accelerazioni del terreno su versanti acclivi e generalmente instabili anche in condizioni normali.

Ad eccezione del fenomeno della liquefazione di cui non sono ancora pervenute segnalazioni, tutti gli altri effetti sembrano essere stati prodotti nell’area epicentrale del terremoto di Amatrice del 24 agosto 2016.

Foto 1– Frattura superficiale con evidenza di spostamento verticale osservata lungo una faglia attiva.

Frattura superficiale con evidenza di spostamento verticale osservata lungo una faglia attiva.

A partire dal 24 Agosto, le squadre di geologi di EMERGEO, uno dei gruppi operativi di emergenza sismica dell’INGV, si sono attivate per effettuare le prime indagini sugli effetti dei terremoti sul territorio, concentrandosi sia sugli effetti cosismici primari (ovvero direttamente legati alla rottura del terremoto come fagliazione e fratturazione superficiale) che su quelli secondari (ovvero riconducibili allo scuotimento provocato dalle onde sismiche come frane, scoscendimenti, crolli di massi, liquefazioni, etc).

Foto 2 – Frana e crollo di massi sulla strada asfaltata.

Frana e crollo di massi sulla strada asfaltata.

Le squadre dei rilevatori stanno investigando l’area epicentrale per una lunghezza totale di circa 40 km tra Castelluccio di Norcia, a Nord, e la Località Ortolano posta a sud del Lago artificiale di Campotosto. Al momento, sono state catalogate informazioni geologiche su circa 2400 punti di osservazione. In generale, nell’intera area investigata sono state segnalate numerose fratture lungo i versanti montuosi e i campi coltivati (questi elementi deformativi risultano spesso ben visibili sul manto stradale), insieme a frane, scoscendimenti e crolli di massi, di piccole-medie dimensioni.

Maggiori informazioni si trovano nel Primo rapporto preliminare del Gruppo operativo EMERGEO sugli effetti cosismici, disponibile in italiano e in inglese.

E’ possibile contribuire alla attività del Gruppo EMERGEO inviando segnalazioni di effetti geologici legati al terremoto tramite il questionario  oppure pubblicando un Tweet con l’hashtag #emergeoINGV indicando data e località precisa (oppure le coordinate geografiche, latitudine e longitudine).

a cura del Gruppo operativo EMERGEO (2016).


Licenza

Licenza Creative Commons
Quest’opera è distribuita con Licenza Creative Commons Attribuzione – Non opere derivate 4.0 Internazionale.

%d blogger cliccano Mi Piace per questo: