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.

Sequenza sismica in Italia centrale: forte aftershock M4.1 stanotte

Ancora una forte replica del terremoto del 24 agosto in provincia di Rieti, questa notte all’1:34 (le 23:34 UTC di ieri 19/8). L’evento, di magnitudo (Richter) 4.1, è stato localizzato nella zona dell’evento principale del 24 agosto, a 4 km da Accumoli. Nelle ore successive sono stati registrati numerosi aftershocks, soprattutto nella zona a nord, verso Norcia, con un evento di magnitudo 3.4 alle 5:30. Sia l’evento di M4,1 che quelli delle ore successive hanno interessato l’area già attiva dal 24 agosto. Tutti i dati sugli eventi recenti si trovano qui.

 

Epicentro del terremoto di questa notte all'1:34, di magnitudo (Richter) 4.1

Epicentro del terremoto di questa notte all’1:34, di magnitudo (Richter) 4.1. I colori si riferiscono alle ore e ai giorni trascorsi precedenti l’ora di questo terremoto, il cui epicentro è rappresentato dalla stella.

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.

Sequenza sismica in Italia centrale: aggiornamento 16 settembre, ore 11:00

Dall’inizio della sequenza, con il terremoto di magnitudo 6.0 avvenuto alle ore 03:36 italiane del 24 agosto, la Rete Sismica Nazionale dell’INGV ha localizzato complessivamente circa 9500 eventi: 186 i terremoti di magnitudo compresa tra 3.0 e 4.0, 15 quelli localizzati di magnitudo compresa tra 4.0 e 5.0 ed uno di magnitudo maggiore di 5.0 (quello di magnitudo 5.4 (Mw 5.3) avvenuto il 24 agosto alle ore 04:33 italiane nella zona di Norcia (PG). 

La mappa della sequenza sismica in Italia centrale aggiornata alle ore 11:00 del 16 settembre.

La mappa della sequenza sismica in Italia centrale aggiornata alle ore 11:00 del 16 settembre.

Dalle ore 11.00 del 12 settembre (ora dell’ultimo aggiornamento) sono stati registrati 10 terremoti di magnitudo maggiore o uguale di 3.0, come riportato nella tabella qui sotto.

Gli eventi più forti, 2 di magnitudo 3.9 ed uno di magnitudo 3.7,  si sono verificati nella giornata di ieri a distanza di pochi  minuti tra le 16:40 e le 16.44 (ora italiana), localizzati in provincia di Perugia, pochi chilometri a sud-est di Norcia.

Data e Ora (UTC) Magnitudo Provincia/Zona
2016-09-15 17:27:54 3.1 Perugia
2016-09-15 16:51:44 3.0 Perugia
2016-09-15 14:44:24 3.9 Perugia
2016-09-15 14:43:51 3.7 Perugia
2016-09-15 14:40:52 3.9 Perugia
2016-09-15 14:40:47 3.1 Perugia
2016-09-15 13:51:29 3.1 Perugia
2016-09-14 16:06:49 3.0 Ascoli Piceno
2016-09-14 03:07:44 3.1 Rieti
2016-09-13 11:04:20 3.4 Rieti
%d blogger cliccano Mi Piace per questo: