Quando iniziò lo sciame sismico all’Aquila?

La sismicità dopo il 6 aprile 2009 in un profilo verticale attraverso la faglia del 6 aprile 2009 (da Valoroso et al., 2013)

Tra pochi giorni verrà emessa la sentenza del processo di appello. Dall’esame dei documenti del primo grado di giudizio, in particolare delle motivazioni della sentenza, rese pubbliche a gennaio 2013 e dai commenti che sono seguiti, è emersa una certa confusione su elementi scientifici di base (per un’analisi dettagliata vedi qui). Ne affrontiamo qui due, legati tra loro: a) quando è iniziata la sequenza sismica sfociata poi nel terremoto del 6 aprile 2009? b) è lecito dal punto di vista scientifico definire quella sequenza uno “sciame” sismico?

L’inizio della sequenza

Per stabilire il momento d’inizio di una sequenza sismica esistono degli algoritmi semplici e accettati dalla comunità scientifica che permettono di distinguere una sequenza (o uno sciame) dalla cosiddetta sismicità di fondo. Quest’ultima è la sismicità che avviene sempre nelle zone sismiche, con tasso di attività variabile da zona a zona ma sostanzialmente costante al suo interno. All’interno della sismicità di fondo possono avvenire saltuariamente terremoti più forti ma, se non cambia il tasso di rilascio sismico, non si può dire che ci sia una sequenza sismica in atto.

Secondo le motivazioni della sentenza la sequenza sarebbe iniziata a giugno 2008. Ad esempio, già a pag. 22, si legge: La scossa in questione (main shock) si verificava nell’ambito di uno sciame sismico che durava già da diversi mesi (dal giugno 2008) e che, prima del 6.4.09, aveva registrato come evento maggiore la scossa di magnitudo 4.1 delle ore 15.38 del 30.3.09.

 La stessa data di giugno 2008 come inizio dello sciame viene citata molte altre volte (alle pagg. 148, 203, 261, 266, 268, 332, 343, 356, 357, 377, 391, 408, 417, 422, 430, 437, 449, 465, 466. 474, 475, 482, 487, 495, 502, 503, 511, 519, 529, 534, 634, 661, 664, 670). In una sola occasione viene fatta iniziare a dicembre 2008 (pag. 647), mentre altre volte viene detto genericamente “iniziato nel 2008”.

Vediamo come stanno le cose nella realtà, applicando l’algoritmo che si usa generalmente per questo tipo di analisi, quello di Paul Reasenberg (1985), sismologo del Geological Survey USA (USGS). L’algoritmo di Reasenberg  analizza le localizzazioni degli eventi sismici avvenuti e calcola il tempo che separa un evento da un altro e la distanza tra di essi: se questi avvengono entro un prefissato intervallo di tempo (per esempio 5 giorni) e in un’area definita (per esempio entro un raggio di 10 km), allora si può ipotizzare che essi appartengano a una sequenza.  Abbiamo chiesto a Franco Mele, sismologo del nostro Istituto esperto in questo tipo di analisi, di verificare se fosse possibile stabilire in maniera chiara il momento d’inizio della sequenza aquilana. Il risultato è che la sequenza è iniziata a gennaio 2009. In particolare, dice l’algortimo di Reasenberg, intorno al 7 gennaio 2009.

 

Distribuzione temporale dei terremoti nella zona aquilana (raggio 30 km intorno alla città) dall’inizio del 2008 al 29 marzo 2009. Ciascun simbolo rosso indica un terremoto la cui magnitudo è leggibile nell’asse verticale a destra. La linea rappresenta la cumulata del numero di eventi sismici (N sull’asse di sinistra è il numero totale): se la pendenza è costante significa un numero costante di terremoti, corrispondente alla sismicità di fondo. Dove si osserva un aumento della pendenza (gennaio 2009) si ha una sequenza o sciame.

La figura sopra mostra chiaramente questo risultato: il tasso di terremoti del 2008 e dei primi giorni del 2009 è stazionario, come mostrato dalla pendenza praticamente costante della retta che indica il numero (cumulato) di eventi che avvengono nel periodo esaminato (trattandosi di un grafico cumulato, una pendenza uniforme vuol dire rilascio costante; in caso di nessun terremoto, la retta sarebbe orizzontale). Il cambio di pendenza (ossia la variazione nel numero medio di terremoti che avvengono nell’area) è abbastanza brusco e si verifica intorno alla metà di gennaio del 2009. Mele ha anche verificato che al variare dei parametri di input (ad esempio della finestra temporale entro la quale considerare due terremoti associati tra loro, in genere 5 giorni) il risultato non cambia. Se si considera un’area più ridotta, di 10 km intorno all’Aquila, la cosa appare ancora più evidente come vediamo nei grafici sotto.

Andamento nel tempo della sismicità dal 1/1/2008 al 30 aprile 2009 in un’area di 30 km intorno a L’Aquila. Ogni punto rappresenta un terremoto di magnitudo come nella scala a sinistra. La linea rossa indica l’inizio della sequenza. Si nota bene che l’andamento prima della linea rossa è variabile ma senza particolari addensamenti (ossia sequenze): è la sismicità di fondo.

 

Andamento nel tempo della sismicità dal 1/1/2008 al 30 aprile 2009 in un’area di 10 km intorno a L’Aquila. Ogni punto rappresenta un terremoto di magnitudo come nella scala a sinistra. La linea rossa indica l’inizio della sequenza. Si vede ancora meglio rispetto alla figura precedente l’inizio della sequenza a gennaio 2009

 

Sciame o sequenza?

Il termine sciame sismico proviene dall’inglese seismic swarm e fu definito probabilmente per la prima volta dai sismologi giapponesi negli anni ’60 del ‘900 (Mogi, 1963; Utsu, 2002). Come per l’inglese swarm, anche in italiano il termine sciame viene usato generalmente per indicare un folto gruppo di insetti o più generalmente una moltitudine di individui o cose in movimento.

Premettiamo che sciame (sismico) e sequenza (sismica) non sono due termini alternativi: il primo è un tipo particolare della seconda, categoria più ampia usata per indicare un addensamento spazio-temporale di terremoti. Risulta sempre difficile stabilire se una sequenza può essere definita uno sciame perché, come scrive Utsu (2002) “non esiste una definizione esatta universalmente accettata di aftershocks, foreshocks, and seismic swarms”. Lo stesso Utsu fornisce la seguente definizione: “uno sciame sismico è una concentrazione (cluster) di terremoti in cui non c’è un singolo terremoto di magnitudo predominante (predominantly large)”. Uno sciame si distingue da una classica sequenza main shock – aftershocks per avere molti terremoti di diverse magnitudo distribuiti irregolarmente nel tempo.

Se guardiamo ai grafici sopra, notiamo come le magnitudo del periodo gennaio – marzo 2009 seguano un andamento di questo tipo, con valori compresi tra 1 e 3 praticamente in tutto il periodo. Quindi, secondo la classificazione di Utsu (2002), nel caso della sismicità aquilana prima del 6 aprile si tratta proprio di uno sciame. Anche i terremoti del periodo tra 30 marzo e 5 aprile rientrano nella tipologia di uno sciame, durante il quale la magnitudo dei terremoti può anche aumentare. Diverso il caso della sequenza dal 6 aprile in poi, in cui a un evento maggiore (ML 5.9, Mw 6.3) sono seguite migliaia di repliche tutte più piccole, che sono diminuite nel tempo in numero e magnitudo seguendo un andamento tipico definito come legge di Omori (1894): una classica sequenza main shock – aftershocks.

Sfatiamo quindi questa credenza, emersa anch’essa durante il processo di primo grado: il termine sciame sismico connota una tipologia di sequenza sismica ben nota da molti decenni e non contiene in sé nessuna valenza rassicurante.

Riferimenti bibliografici

Mogi K. (1963). Some Discussions on Aftershocks, Foreshocks and Earthquake Swarms – the Fracture of a Semi-infinite Body caused by an Inner Stress Origin and its Relation to the Earthquake Phenomena (third paper). Bull. Earthquake Res. Institute, 41, pp. 615-658.

Omori, F. (1894) On the aftershocks of earthquakes. Journal of the College of Science, Imperial University of Tokyo 7: 111–200

Reasenberg, P. (1985) Second-Order Moment of Central California Seismicity, 1969-1982. Journal of Geophysical Research, 90, 5479-5495.

Utsu T. (2002) Statistical features of seismicity. In: International Handbook of Earthquake and Engineering Seismology. Editors: W.H.K. Lee, H. Kanamori, P.C. Jennings, C. Kisslinger. International Geophysics, Volume 81, Part A. Elsevier, ISBN: 978-0-12-440652-0, p. 719-732.

Valoroso, L., L. Chiaraluce, D. Piccinini, R. Di Stefano, D. Schaff, and F. Waldhauser (2013), Radiography of a normal fault system by 64,000 high-precision earthquake locations: The 2009 L’Aquila (central Italy) case study, J. Geophys. Res., 118, doi:10.1029/2012JB009927.