Le glaciazioni tra fine Devoniano e inizio del Carbonifero che precedono la grande fase glaciale del Paleozoic superiore e il ruolo svolto in proposito dal tenore di CO2 atmosferico

Questo post è in qualche modo il seguito del post sulle estinzioni di massa del Devoniano superiore, perché parla di una fase cruciale nella storia del nostro pianeta, il passaggio avvenuto tra Devoniano e Carbonifero fra l’ambiente caldo del Paleozoico inferiore e quello decisamente più freddo del Paleozoico inferiore, caratterizzato dalle grandi glaciazioni del Gondwana del Permo – Carbonifero. Queste glaciazioni hanno una differenza fondamentale con quelle attuali: se dalla fine del Terziario inferiore ad oggi la presenza di ghiacci è in buona parte attribuibile ad un importante cambio nella paleogeografia e cioè all’approssimarsi verso le alte latitudini di entrambi gli emisferi di importanti masse continentali (Antartide, America Settentrionale ed Eurasia), nel Paleozoico sostanzialmente il Gondwana risiedeva “da un bel pò” nei dintorni del Polo Sud e verso le alte latitudini settentrionali c’erano ben poche terre emerse, per cui l’avvento delle glaciazioni non è  legato all’arrivo di estese masse continentali nelle aree polari. Questa circostanza fa esplorare altre soluzioni per spiegare l’enorme differenza climatica fra Paleozoico inferiore e superiore, in particolare il possibile ruolo giocato dal tenore atmosferico di CO2 in tutti questi avvenimenti. Con gli ovvi riflessi per la situazione attuale.

Il Tiktaalik rosae, un classico antenato dei tetrapodi,
vissuto nel Devoniano superiore della Groenlandia,
all'epoca in posizione equatoriale nel nord di Euramerica

IL DEVONIANO SUPERIORE: UNO DEI MOMENTI CRUCIALI DELLA STORIA DELLA TERRA. Nel Devoniano accadono degli avvenimenti particolarmente importanti. Dal punto di vista biologico, oltre alla definitiva affermazione delle piante terrestri, sulle terre emerse arrivano anche i vertebrati, con i primi anfibi; questo passaggio è sempre più conosciuto e le ultime ricerche hanno evidenziato una diffusione degli antenati dei tetrapodi da Euroamerica al Gondwana prima del previsto (Gess e Ahlberg, 2018), dimostrando che tutto sommato i due continenti fossero già piuttosto vicini in quest'epoca. Dal punto di vista geologico registriamo l’accelerazione degli eventi che porteranno alla formazione della Pangea (la quale, checché se ne dica, non comprendeva tutte le terre emerse, ne parlerò, spero presto, in un prossimo post) e alle glaciazioni del Paleozoico superiore. 

Tutto questo non avviene in modo tranquillo, ma anche attraverso due estinzioni di massa piuttosto importanti di cui ho parlato qui: la prima, al passaggio Frasniano – Famenniano è una delle “big five” (le maggiori estinzioni di massa) riconosciute già da Sepkoski (1996), la seconda alla fine del Famenniano e quindi in corrispondenza del passaggio Devoniano  - Carbonifero. Come al solito c’è una sincronia fra estinzioni di massa e messa in posto di Large Igneous Provinces: nel primo caso si tratta della LIP della Yacuzia, nel secondo, come ho già scritto, ci sono ancora dubbi sulla collocazione dell’evento (ma non sulla sua esistenza!): per alcuni autori si potrebbe addirittura trattare di una seconda fase della stessa LIP asiatica (Ricci et al 2013). 

Paleogeografia tra Devoniano superiore e carbonifero inferiore
Le stelle indicano i depositi glaciali della fine del Famenniano (Modificata da Larkin et al 2016)
N.Am = Nord america, S.Am = Sud America. Af = Africa 

LA PALEOGEOGRAFIA TRA DEVONIANO SUPERIORE E CARBONIFERO INFERIORE. La caratteristica dominante della paleogeografia del Paleozoico (e quindi anche di questi momenti) è la distribuzione sbilanciata delle masse continentali verso l’emisfero meridionale, come lo è quella attuale verso l’emisfero settentrionale: posizionato alle alte latitudini meridionali, il grande continente di Gondwana comprendeva oltre alle parti che si sarebbero frammentate nel Mesozoico formando Africa, Sud America, Australia, India, Antartide, Arabia, qualche altro frammento che si stava staccando in quel momento; anche la maggior parte di Euramerica (o Laurussia) si trovava nell’emisfero meridionale: solo le sue propaggini più settentrionali arrivavano fino all’equatore e un po' oltre come dimostrano i fossili delle Svalbard. Nell’emisfero settentrionale invece troviamo tra latitudini medio – alte e l’equatore solo alcune masse minori isolate fra loro, quelle che di lì a poco formeranno l’ossatura primaria dell’Asia con la fascia orogenica dell’Asia Centrale (ne ho parlato qui): Siberia, Amuria, l’arco dei Verhojansk, il Kazakhstan e le due Cine.

DA UNA TERRA CALDA A UNA TERRA FREDDA. Un aspetto ancora poco conosciuto di questi momenti è come sia avvenuto il passaggio fra la Terra calda del Paleozoico inferiore e le glaciazioni del Permocarbonifero (in inglese si dice da una Terra a “greenhouse” ad una Terra a “icehouse”) e non si sa neanche di preciso quando poter collocare questo limite. La presenza di depositi glaciali diventa continua dalla fine del Mississipiano, intorno a 320 milioni di anni fa, ma prima dell’avvento delle calotte più massicce e permanenti del Carbonifero superiore, si registrano almeno tre fasi glaciali minori, appunto fra Devoniano superiore e Carbonifero inferiore da 355 milioni di anni fa. A questo si aggiunge un po' di confusione nella nomenclatura della cronologia del Carbonifero in cui fra Europa e America settentrionale sono stati usati termini diversi, che necessitano di esser quindi confrontati con la cronologia attualmente vigente, che si vede qui accanto (International Chronostratigraphic chart, versione 2017).  

LE TRE FASI GLACIALI ISOLATE DI QUEI TEMPI. Non ci sono ancora molti dati su questi eventi (e non so quanti altri se ne potranno trovare); una ottima rassegna in materia è contenuta in Larkin et al (2016), testo a cui faccio riferimento dove non specificato diversamente quando elenco le zone interessate. Ad oggi sono stati riconosciuti tre eventi glaciali: uno alla fine del Devoniano, che fa parte della catena di avvenimenti dell’evento Hangenberg, e due nel Carbonifero inferiore, nel Tournasiano e uno nel Viseano. 

Le tracce fondamentali che permettono di riconoscere un evento glaciale sono diverse: la presenza di diamictiti (morene glaciali fossili) è sicuramente l’aspetto più diagnostico, al pari di una diminuzione del livello marino a scala mondiale che provoca l'interruzione della sedimentazione un pò dappertutto nelle zone precdentemete ricoperte da mari di bassa profondità; inoltre si riscontra una migrazione verso latitudini più basse di forme di vita di clima più fresco.

Il ciclo glaciale della fine del Devoniano (Famenniano finale) è particolarmente intrigante perché, appunto, si colloca vicino all’evento di Hangenberg in senso stretto. 

Variazioni del livello marino e stratigrafia
in Germania tra fine Devoniano e inizio Carbonifero
da Bless et al (1992)

Anzi, lo precede: il Famenniano terminale è stato interessato da una glaciazione le cui diamictiti si trovano in Africa e nelle Americhe, e la loro presenza nella fascia appalachiana dell’America Settentrionale significa che siamo davanti ad un intervallo glaciale che è arrivato fino a basse latitudini; inoltre è attestato a livello mondiale un generale ed improvviso abbassamento del livello marino di un centinaio di metri, evidentissimo in Germania nel Massiccio Renano (Kaiser et al 2015), ma anche tra Francia e Belgio (Kumpan et al, 2014), America settentrionale e in tutta l’area dell’Orogene dell’Asia Centrale. In particolare nel massiccio renano il canale di Seiler è una incisione profonda 100 metri scavata dall'erosione nella piattaforma continentale di Euramerica durante la fase a mare basso del Famenniano superiore e nella quale si sono sedimentati, quando il livello marino è risalito durante la fase di deglaciazione, il conglomerasto di Seiler e gli scisti di Hangenberg che registrano il successivo evento omonimo (Bless et al, 1992). Il ciclo di ritiro del mare e avanzamento dei ghiacci seguito dal ritorno alle condizioni precedenti sembra essere durato non più di 100.000 anni (Streel et al, 2013). La deposizione di scisti scuri perché ad alto contenuto di materia organica (che quindi si sono deposti in una fase in cui di ossigeno nelle acque ce n’era ben poco) è un evento a scala mondiale che si associa ad una perturbazione nel ciclo del Carbonio. 

La catena di eventi dell’Hangenberg è quella classica che si verifica in corrispondenza della messa in posto di una Large Igneous Province: un raffreddamento associato ad una regressione marina, successivamente seguito da una trasgressione con deposizione di sedimenti ricchi in materia organica, estinzione di massa e perturbazioni geochimiche, fra le quali la più nota è quella del rapporto fra gli isotopi del Carbonio, ma ce ne sono diverse altre.

Una valle incisa da una regressione marina,
che viene successivamente riempita da sedimenti
durante la successiva trasgressone (da Moxness et al, 2018)

Gli altri due episodi glaciali sono decisamente di minore portata e i depositi glaciali si ritrovano solo nelle alte latitudini meridionali dell'epoca: quello del medio Tournaisiano  (circa 350 MA) è attestato in Brasile, Falkland e Sudafrica, quello del Viseano in Brasile, Argentina, Falkland e Malesia. Le modalità  dei due eventi sono molto diverse fra loro; il Medio Tournasiano ricorda il Famenniano: un solo ciclo di regressione marina con  formazione di ghiacciai seguita da un riscaldamento contrassegnato da una deposizione di sedimenti anossici in Germania e altrove, accompagnato da anomalie geochimiche; Insomma, si tratta di due eventi a scala mondiale ma temporanei e monociclici, decisamente diversi da quelli classici successivi del Permo-Carbonifero e del Terziario, caratterizzati da continue alternanze di trasgressioni e regrssioni marine. Quindi sono portato a pensare che anche l’episodio del Tournasiano sia legato alla messa in posto di una Large Igneous Province. Sulle motivazioni del raffreddamento che precede la messa in posto di una Large Igneous Province ci sono ancora poche ipotesi. Ricordo che negli anni ‘80, quando ancora i flood basalts continentali e i plateau oceanici non erano ancora stati uniti nelle Large Igneous Provinces e le prove di collegamento fra queste grandi eruzioni e le estinzioni non erano così stringenti come oggi, era stata avanzata una correlaizone fra le estinzioni di massa e le variaizoni del livello marino (Hallam, 1989)

Nel Viseano, poco più di 350 milioni di anni fa, le cose appaiono un pò diverse: a differenza dei precedenti, questo intervallo sembra essere policiclico e non presenta una trasgressione contrassegnata da black shales e anomalie geochimiche e non c’è una particolare accelerazione del tasso di estinzione (Stanley, 2016). Insomma, a differenza degli episodi precedenti, quello del viseano appare più come una fase glaciale normale, un annuncio di quello che sarebbe successo da lì a poco dopo.

C’è solo una cosa che mi stupisce: non capisco come mai due eventi singoli come Hangenberg e medio Tournasiano abbiano lasciato così tante tracce.

GLACIAZIONE Sì E GLACIAZIONE NO A PARITÀ DI PALEOGEOGRAFIA. C’è poi da capire perché le glaciazioni si verificano solo nel Paleozoico superiore quando già in quello inferiore le condizioni dal punto di vista paleogeografico per averle, e cioè la presenza di una ampia area continentale nella zona polare meridionale, c'erano di già. Per alcuni Autori uno dei fattori che lo avrebbero determinato è stato la riorganizzaione della circolazione oceanica a causa dello scontro fra Euramerica e Gondwana e la chiusura del passaggio fra Panthalassa e Tetide, ma per altri Autori si tratta di un evento avvenuto più tardi (per una rassegna in merito: Davydov e Kozar, 2017). Noto invece che per la glaciazione attuale uno dei maggiori eventi che l'hanno innescata potrebbe essere stato, al contrario, la nascita della corrente circumpolare quando si è aperta la soglia di Ross fra Argentina e Penisola Antartica (Katz et al, 2011).

Il tenore di CO2 atmosferico dal Devoniano ad oggi (Franks et al, 2014)
Si vede come le fasi glaciali iniziano quando questo tenore si abbassa

Qui entra in gioco, però, un'altra importante forzante: la figura qui accanto, rielaborata da Franks et al (2014) descrive il contenuto medio di CO2 atmosferico dal Devoniano a noi, calcolato analizzando gli stomi delle foglie fossili, la cui densità è correlabile appunto al tenore atmosferico di CO2. Si vede che dopo un massimo nel Devoniano inferiore, il tenore atmosferico di CO2 inizia a diminuire. Questa diminuzione significa che nel sistema - Terra il CO2 prodotto non era sufficiente a compensare la sua richiesta da parte dei fenomeni che lo consumano (essenzialmente fotosintesi, alterazione delle rocce silicatiche e formazione di rocce carbonatiche). È molto probabile che la colonizzazione dei continenti da parte delle piante sia stata la nuova variabile che ha provocato una nuova e massiccia richiesta di CO2; si deve notare inoltre come il calo si interrompe bruscamente (e momentaneamente) nel Famenniano, l‘ultimo piano del periodo, per proseguire poi per tutto il Carbonifero in quanto una grande quantità di gas verrà stoccato nelle calotte glaciali. Si rialzerà solo alla fine del Permiano. Ho già ipotizzato che questo momentaneo stop del calo sia dovuto alle emissioni di CO2 delle due importanti Large Igneous Provinces della fine del Devoniano (Piombino, 2017). Questo calo ha diminuito l’effetto – serra, provocando appunto l’inizio delle glaciazioni, mentre il tenore elevato di CO2 dell’inizio del Devoniano potrebbe essere la chiave del perché con una situazione paleogeografica sostanzialmente analoga non ci sinao tracce di evidenti di glaciazioni.

Sul raffreddamento pre-LIP e pre estinzioni di massa ci sono ancora poche ipotesi. La mia idea è che la fase principale dell’attività magmaticasia preceduta da una drammatica emissione di tufi, cioè che avvenga all’ennesima potenza quello che succede dopo una massicica eruzione esplosiva di un vulcano: le particelle coprono la Terra e bloccano la radiazione solare e se in una Terra “greenhouse” il calore latente è ancora piuttosto elevato e quindi non troviamo glaciazioni importanti, quando il tenore atmosferico di CO2 diminuisce il blocco della radiazione solare provoca un raffreddamento. Questo raffreddamento però si conclude molto presto perché, appena inizia l’attività effusiva, gli immensi quantitativi di CO2 da essa emessi provvedono a ristabilire l’effetto serra.

Insomma, una ennesima dimostrazione nel passato di quanto conti per le temperature globali il tenore atmosferico dei gas-serra.


Bless et al (1992) Eustatic cycles around the Devonian-Carboniferous boundary and the sedimentary and fossil record in Sauerland (Federal Republic of Germany). Annales de la Societe Geologique de Belgique, 115, 689-702

Davydov e Kozar (2017) The formation of the Alleghenian Isthmus triggered the Bashkirian glaciation: Constraints from warm-water benthic foraminifera Palaeogeography, Palaeoclimatology, Palaeoecology in press

Gess e Ahlberg (2018) A tetrapod fauna from within the Devonian Antarctic CircleScience 360, 1120–1124  

Hallam (1989) The case for sea-level change as a dominant causal factor in mass extinction of marine invertebrates. Philos. Trans. R. Soc. London, Ser. B 325, 437–455.

Katz et al (2011). Impact of Antarctic Circumpolar Current Development on Late Paleogene Ocean Structure, Science 332,1076–1079

Kumpan et al (2014) Sea-level and environmental changes around the Devonian–Carboniferous boundary in the Namur–Dinant Basin (S Belgium, NE France): A multi-proxy stratigraphic analysis of carbonate ramp archives and its use in regional and interregional correlations Sedimentary Geology 311,43-59

Lakin et al (2016) Greenhouse to icehouse: a biostratigraphic review of latest Devonian–Mississippian glaciations and their global effects Geological Society, London, Special Publications, 423, 439-464, 15 April 2016

Moxness et al (2018) Sedimentology of the mid-Carboniferous fill of the Olta paleovalley, eastern Paganzo Basin, Argentina: Implications for glaciation and controls on diachronous deglaciation in western Gondwana during the late Paleozoic Ice Age journal Of South American Earth Sciences 84, 127-148

Piombino (2016) The Heavy Links between Geological Events and Vascular Plants Evolution: A Brief Outline International Journal of Evolutionary Biology volume 2016, Article ID 9264357 http://dx.doi.org/10.1155/2016/9264357

Ricci et al (2013) New 40Ar/39Ar and K–Ar ages of the Viluy traps (Eastern Siberia): Further evidence for a relationship with the Frasnian–Famennian mass extinction Palaeogeography, Palaeoclimatology, Palaeoecology 386 (2013) 531–540 

Sepkoski (1996) Patterns of Phanerozoic extinction: A perspective from global databases. Global Events and Event Stratigraphy in the Phanerozoic, ed Walliser OH (Springer, Berlin), pp 35–51

Stanley (2016) Estimates of the magnitudes of major marine mass extinctions in earth history PNAS 113, E6325-E6334

Streel et al (2013) What do latest Famennian and Mississippian miospores from South American diamictites tell us? Palaeobiodiversity and Palaeoenvironments, 93, 299-316

I mammiferi terrestri siciliani, la loro storia e un istmo nello stretto di Messina

Prima della fine del Miocene i territori oggi facenti parte della Sicilia appartenevano a due piattaforme continentali diverse in cui le isole erano di ridotte dimensioni e sporadiche nello spazio e nel tempo. Da quel momento i primi lembi dell’isola sono emersi stabilmente e, complici sia gli eventi tettonici che, nel quaternario, le oscillazioni del livello marino, il suo territorio ha variato spesso configurazione, da quella a più isole (tre grandi, la catena settentrionale, gli Iblei e Malta) a una terra molto vasta che arrivava oltre Malta. Per la sua posizione al centro del Mediterraneo, la Sicilia può essere vista come una via principale per le migrazioni tra le coste meridionali e settentrionali del bacino, ma questo è ciò che vediamo ora: in un passato anche molto recente le comunicazioni fra l'isola (o, meglio, le isole che componevano l’arcipelago siculo – maltese) e le aree circostanti sono state sporadiche. La colonizzazione della Sicilia da parte dei vertebrati terrestri è stata quindi molto complessa, come anche la colonizzazione umana che sembra essere avvenuta molto in ritardo rispetto alle altre sponde del Mediterraneo. Ci sono diverse questioni in campo che sono ancora aperte e irrisolte o forniscono risultati contraddittori sulle provenienze faunistiche, che potevano sfruttare il cosiddetto istmo siculo-tunisino o la penisola Italiana (ma possono essere vere entrambe...). Su queste istanze prettamente zoologiche si innesta quella antropologica: quando è avvenuta la colonizzazione umana e da dove provenivano questi uomini? 

L'Italia nel Pliocene: tante isole separate.

I RARI FOSSILI DEL MIOCENE E L’ASSENZA DI TESTIMONIANZE NEL PLIOCENE. La Sicilia è un'area in cui due placche si sono scontrate tra l’inizio del Miocene e l’inizio del Pliocene (grossolanamente tra 20 e 5 milioni di anni fa) ed è composta da tre domini principali: la catena settentrionale, l’avampaese Ibleo e, in mezzo, l’avanfossa di Gela e Caltanissetta. La zona degli Iblei, il settore sud-orientale dell'isola, prosegue nel canale di Sicilia fino all'arcipelago di Malta, e rappresenta il margine settentrionale della placca africana.

Prima del Miocene, le aree ora comprese nella Sicilia appartenevano a piattaforme continentali in cui le isole erano saltuarie nel tempo e nello spazio, ma dall’inizio di questo periodo, circa 20 milioni di anni fa, le cose iniziano a cambiare perchè le due piattaforme si scontrano e nella zona degli attuali Iblei troviamo i primi segni della futura e definitiva emersione: una serie di isole, nelle quali si trovano denti di una forma nana di Mastodonte di chiara affinità africana (provenienza logica, dato che era una continuazione della piattaforma di quel conitnente) e che rappresentano le più antiche tracce di mammiferi terrestri in Sicilia. Più tardi, verso la fine del Miocene (ma non in quello terminale) è stata ritrovata una fauna che include taxa euroasiatici e africani, senza particolari endemismi, che però è stata distrutta durante il terremoto del 1908. La duplice provenienza potrebbe essere spiegata da un mix di forme di provenienza iblea e da nord (Rook et al, 2006) 

Se i fossili del Miocene sono frammentari, quelli del periodo successivo, il Pliocene, sono praticamente assenti: è stato un periodo di alto mare in cui l'Italia e il Nordafrica erano come le attuali Filippine: una serie di isole (molto piccole…) separate da bracci di mare più o meno ampi. Non certo un posto ideale per i mammiferi terrestri.

I COMPLESSI FAUNISTICI FOSSILI DEL PLEISTOCENE IN SICILIA. Nel Pleistocene la situazione cambia notevolmente: specialmente durante le fasi a basso livello marino dei massimi glaciali, al posto dell'arcipelago Siculo - Maltese si forma una vasta area emersa che comrende la Sicilia, Malta e vaste aree del canale di Sicilia. Da questa situauzone hanno tratto giovamento i mammiferi terrestri, di cui troviamo numerose testimonianze, specialmente nei tempi più recenti. Sono stati individuati diversi complessi faunistici che cambiano nel tempo (Bonfiglio et al, 2002), ben descritti anche in Sineo et al (2015)

1. Complesso faunistico di Monte Pellegrino. È il più antico: si tratta di una fauna scarsamente diversificata, composta principalmente da piccoli mammiferi probabilmente discendenti della popolazione della fine del Miocene che hanno eroicamente resistito alle difficoltà del Pliocene. Sono forme estremamente particolari, con un endemismo pronunciato che riflette il lungo isolamento. Da questo momento l'attuale area della Sicilia emerge con una certa continuità sia nel tempo che nello spazio e cominciano ad arrivare cose diverse.

La regione italiana nei momenti di basso livello marino
durante i massimi glaciali: Sicilia e Malta unite

2. Complesso faunistico di Elephas (Palaeoloxodon) falconeri. È dominato, come attesta il nome da questo elefante nano. Si tratta di un mix fra elementi di origini africane e elementi locali antichi come quelli del complesso faunistico precedente. È comunque una fauna povera, che riflette una situazione di scarsa continuità tra le varie isole che componevano l’arcipelago siculo – maltese e dove tranne il proboscidato tutti gli altri esponenti erano di dimensioni estremamente ridotte. Le poche datazioni parlano di circa 400.000 anni fa

3. Complesso faunistico di Palaeoloxodon mnaidriensis. Anche questo complesso prende il nome da un proboscidato, il quale è sempre una forma abbastanza piccola, ma adesso compaiono altri mammiferi di dimensioni medie come bisonti, uri, ippopotami, cervi rossi, daini, cinghiali, iene maculate, lupi e volpi. Con la rilevante eccezione di un ippopotamo (nano) si tratta di forme chiaramente europee. Il cambiamento non ha riguardato i mammiferi di piccole dimensioni, che continuano più o meno ad essere come quelli di prima. Le datazioni sui reperti forniscono delle età comprese fra 200 e 145 mila anni fa.

Entrambi questi complessi sono noti pure a Malta (Hunt & Schembri, 1999)

4. Complesso faunistico San Teodoro Cave–Pianetti. Abbastanza simile al precedente per i macromammiferi, a parte la perdita dell’ippopotamo e la comparsa di un cavallo, si discosta fortemente dal precedente per i piccoli mammiferi, che vengono sostituiti completamente. La comparsa del cavallo dimostra la scomparsa dei boschi in favore delle steppe (ma anche qualcos’altro, ne leggerete subito dopo) ed in effetti siamo esattamente intorno all’ultimo massimo glaciale di 20.000 anni fa

5. Complesso faunistico di Castello. Si colloca dopo il massimo glaciale ed è caratterizzato dalla fine degli endemismi residui: la fauna è praticamente la stessa dell’Italia continentale.

Quindi nel Pleistocene superiore, da 300.000 anni fa in poi, si assiste al ribaltamento delle provenienze faunistiche,  dall’Africa all’Europa. Immagino che la cosa rifletta importanti cambiamenti paleogeografici: il ritiro dei mari interni della penisola italiana e l’approfondirsi nel Canale di Sicilia della soglia fra Mediterraneo Occidentale e Orientale, con la definitiva perdita della continuità della piattaforma continentale tra Italia e Africa.

Però il massiccio arrivi di piccoli mammiferi carnivori, dei cavalli e dei leoni dimostrerebbe che anche lo stretto di Messina non esisteva più durante l’ultimo massimo glaciale (LGM)

Andiamo con ordine.

Confronto fra uomo, elefanti nani della Sicilia e Elephas antiquuus
(tipica forma europea quaternaria) da Benton (2010) 

TEORIA DEL POPOLAMENTO DELLE ISOLE. È indubbio che i mammiferi terrestri in qualche caso abbiano colonizzato delle isole a relativa distanza dalla terraferma. Il caso classico è quello del Madagascar, dove l’arrivo dei mammiferi placentati è posteriore alla separazione dell’isola dall’Africa. Ma ci sono altri casi, come per esempio i Primati e i roditori che nel Paleocene sono passati dall’Europa all’Africa e nell’Eocene dall’Africa in America meridionale e, sempre nel continente latino – americano, l’invasione degli ungulati del nordamerica fino a formare i gruppi tipici ora estinti come i Notoungulati e i Litopterni nel Paleocene. 

Nelle isole spesso le dinamiche sono diverse da quelle della terraferma e difatti un metodo valido per riconoscere una fauna insulare è il nanismo delle forme in teoria più grandi, dovuto alla mancanza di predatori e ad una minore abbondanza di risorse, mentre quelle più piccole tendono ad aumentare di dimensione.

Per conquistare un’isola ci vogliono diverse qualità. La prima è intuitiva: la capacità, più o meno forzata, di passare un braccio di mare; la seconda è la capacità di adattamento, che in generale significa “riesco a mangiare tutto quello che trovo”. 

Una parte delle colonizzazioni sono avvenute in modo talmente improbabile  da essere chiamato la "vittoria nella lotteria", che può essere “a nuoto” o “a cavallo di una zattera”: in concomitanza di forti eventi atmosferici la corrente ha trascinato in mare degli individui che nuotando o attaccati a delle zattere vegetali sono riusciti a sbarcare in una nuova isola anziché annegare o morire di sete. La cosa è stata dimostrata recentemente nei Caraibi con delle iguane che dopo un uragano sono arrivate in un’isola dove non c’erano loro conspecifici (Dawkins, 2004).

Altre situazioni sono la dispersione attraverso corridoi moderatamente ampi formati da aree a basso fondo marino o sistemi di lagune parzialmente emerse, con il caso estremo della formazione temporanea di vere connessioni completamente terrestri.

Dal punto di vista della "lotteria", mammiferi di piccole dimensioni hanno la possibilità di resistere abbastanza a lungo su delle zattere a patto che le foglie forniscano acqua dolce, mentre quelli di grandi dimensioni hanno a disposizione solo la capacità di nuotare (e infatti a parte, forse, gli antenati degli ungulati americani e gli elefanti, si tratta sempre di forme leggere e di dimensioni ridotte come piccoli Primati, Roditori o piccoli insettivori come i Tenrec).

Elefanti e cervidi sono buoni nuotatori (gli elefanti anche in mare, anche se suona strano ...) e difatti i loro fossili sono diffusi in molti siti che erano isole ad una certa distanza dalla terraferma non solo nel Mediterraneo, ma anche in Cina, Indonesia e California (Johnson, 1980). Viceversa gli ippopotami non nuotano, ma camminano sul fondo del corso d’acqua.

Ma oltre alle qualità del nuoto occorrono anche le capacità di opportunismo alimentare: i cervidi richiedono un foraggio altamente nutriente, tuttavia possono passare facilmente da un tipo di cibo all'altro. I bovini invece evolvono specializzazioni sofisticate per un certo tipo di cibo. Pertanto, se i cervidi non possono competere con i bovidi negli ecosistemi maturi, la loro strategia opportunistica li rende buoni colonizzatori che possono sopravvivere in nuove aree, dove i bovidi specializzati non riescono a nutrirsi. Pertanto, l'opportunismo ecologico e l'assenza di specializzazioni alimentari danno ai cervidi delle buone possibilità di sopravvivenza. Anche ippopotami e proboscidati non sono molto schizzinosi in fatto di cibo.

Insomma, mentre gli elefanti e i cervidi sono buoni nuotatori e si adattano a tutto (e quindi sono dei perfetti saltatori), gli ippopotami nuotano molto male ma si adattano a tutto e i bovidi magari nuotano ma sono molto limitati nel tipo di cibo. Invece i cavalli sono la negazione totale di questa specialità: non nuotano e hanno bisogno di cibo specifico e abbondante (difficile in un'isola secialmente la seconda..). Quindi trovare un cavallo in Sicilia significa che il passaggio via terra c’era davvero..

lo stretto di Messina con la soglia tra Villa San Giovanni e Ganzirri, da Antonioli et al (2014)

IL PASSAGGIO SULLO STRETTO DI MESSINA. Abbiamo visto che il complesso faunistico di San Teodoro – Cave Pianetti differisce da quello precedente (quello del Palaeoloxodon mnaidriensis) per l’arrivo di leoni, cavalli e diversi piccoli mammiferi. Un piccolo mammifero è un conto, può essere un caso di vincita alla lotteria, ma quando sono diverse specie la cosa “puzza” un pochino. Se poi ci mettiamo anche leoni e cavalli, che a fare salti in mezzo al mare non ci pensano assolutamente, e che dove c’è cavallo c’è steppa (e quindi un clima decisamente più freddo) si può ipotizzare che il tutto sia avvenuto in coincidenza di un momento di freddo. L'ipotesi è giusta: questo complesso faunistico si colloca esattamente intorno all’ultimo massimo glaciale di 20.000 anni fa. Ma oltre al freddo c’è anche l’abbassamento del livello marino. Abbassamento che, come dimostrano Antonioli et al (2014) ha portato all’emersione di una soglia posta più o meno tra Ganzirri e Villa San Giovanni nel periodo dell’LGM, tra 22 e 19 mila anni fa,  quando solo 100.000 anni prima il braccio di mare era largo circa il doppio di adesso.

Nell’immagine si vede come già con un livello marino di 108 metri più basso di quello attuale, un istmo collegherebbe attualmente le due sponde. Il lavoro prende in esame diversi parametri, a partire dall’erosione delle correnti, per modellizzare l’altezza di questa soglia. 

La presenza di questo ponte di terra è quindi provata sia dal lato geologico che da quello biogeografico. È inoltre possibile (e le datazioni sono coerenti con questo) che anche il passaggio fra i due complessi faunistici ad elefanti sia avvenuto durante il periodo glaciale precedente nella stessa zona, vista la presenza di ippopotami, bovini, iene, lupi e volpi, animali non particolarmente noti per avventure acquatiche, e che a parte l’ippopotamo si tratta di apporti tipicamente europei.

LA COLONIZZAZIONE UMANA DELLA SICILIA. Veniamo adesso all’ultima questione aperta, la colonizzazione umana. Ci sono possibili evidenze di presenza umana, forse sporadica, tra 32.000 e 27.000 anni fa; in particolare qualche Autore riporta dei reperti risalenti all'Aurignaziano (Chilardi et al, 1996). Tuttavia, questi resti, estremamente frammentari, hanno uno scarso valore diagnostico perché non è neanche sicura la loro collocazione stratigrafica. Reperti sicuri esistono solo negli ultimi 15.000 anni,  . In questo caso la soglia nello stretto di Messina potrebbe aver giocato un ruolo fondamentale per questa prima migrazione nell’isola. Morfologia e genetica dimostrano che questi colonizzatori hanno una chiara affinità europea, escludendo così un ampio contributo dal lato opposto del Mediterraneo (Sineo et al 2015). Resta la curiosità di uno stanziamento nell’isola estremamente recente, addirittura più recente che nella più lontana dal continente Sardegna.

Antonioli et al (2014) Timing of the emergence of the Europe–Sicily bridge (40–17 cal ka BP) and its implications for the spread of modern Humans 
Benton et al (2010) Dinosaurs and the island rule: Dinosaurs from Hateg Island, Palaeogeography, Palaeoclimatology, Palaeoecology (2010), doi:10.1016/j.palaeo.2010.01.026.
Bonfiglio et al (2002) Pleistocene calabrian and sicilian bioprovinces Geobios, mem.spec. 24

Chilardi et al (1996). Fontana Nuova di Ragusa (Sicily, Italy): southernmost Aurignacian site in Europe. Antiquity, 70, 553–563.

Dawkins (2004). The Ancestor Tale Houghton Mifflin Harcourt publishing, Boston 
Hunt & Schembri (1999) Quaternary environments and biogeography of the Maltese Islands.
In: Mifsud, A. & Savona Ventura, C. [eds] Facets of Maltese prehistory. pp. 41-75; Malta: The Prehistoric Society of Malta; vii + 243pp
Johnson 1980. Problems in the land vertebrate zoo- geography of certain islands and the swimming powers of elephants. Journal of Biogeography, 7, 383–398.
Rook et al (2006) Lands and endemic mammals in the Late Miocene of Italy: Constrains for paleogeographic outlines of Tyrrhenian area. Palaeogeography, Palaeoclimatology, Palaeoecology, 238, pp. 263-269
Sineo et al (2015) Human Peopling of Sicily During Quaternary  - The Quaternary Period AcademyPublish.org 

La deformazione dell'area del monte Mantap in Corea del Nord a seguito dell'esperimento termonucleare del 3 settembre

Le applicazioni del radar interferometrico satellitare (InSAR) sono numerose. In questa occasione i dati del satellite tedesco TerraSAR-X hanno determinato le deformazioni seguite all'ultimo esperimento nucleare in Corea del Nord, dimostrando come questa esplosione abbia pesantemente modificato la topografia del sito, rendendolo inutiilizzabile per attività future e costringendo il regime nordcoreano a cambiare la sua strategia. È stato inoltre possibile determinare con ls deformazioned ella superficie topografica un valore di Mw che oltre all'energia inziale comprende anche le deformazion post-evento e che sostanzialmente è in accordo con quello determinato esclusivamente con metodi geofisici.

Il 3 Settembre 2017 i sismografi di tutto il mondo hanno registrato un evento sismico localizzato nel poligono nucleare di Punggye-ri nella Corea del Nord, indicando che nel sito era stato appena eseguito un altro esperimento nucleare, e, soprattutto, che questo ultimo esperimento è stato ben più potente dei precedenti. Poco dopo l’agenzia ufficiale dello Stato ha confermato quello che molti esperti temevano: il Paese asiatico aveva fatto esplodere una bomba termonucleare a due stadi (una normale bomba a fissione ha innescato una bomba a fusione, la classica trafila di una bomba H). Già a novembre in questo post ho descritto la geologia del sito nucleare e fatto notare che a seguito dell’esperimento di settembre qualcosa non aveva funzionato a dovere: il primo allarme è stato costituito da una replica M 4.1 registrata appena 8 minuti e mezzo dopo l’esperimento. A questa prima replica ne sono seguite altre due, il 23 settembre e il 12 ottobre; secondo fonti solitamente ben informate in questa ultima occasione numerosi addetti ai lavori sarebbero morti per il crollo di un tunnel che stavano scavando. Insomma, già ai primi di novembre iniziavano i sospetti che la nuova strategia del regime coreano sulle armi nucleari fosse dovuta aulla impossibilità di proseguire l’attività a  Punggye-ri, ipotesi che è diventata certezza qualche tempo dopo che ho scritto quel post. Oggi i radar montati su satellite consentono grazie alla tecnica InSAR di vedere chiaramente le componenti verticale ed E-W della deformazione associata all’evento e dare un’altra valutazione indipendente da quella geofisica sull’energia rilasciata nell’esperimento.

Le deformazioni dopo gli esperimenti
del 1992 da Vincent et al (2003)

LE DEFORMAZIONI RILEVATE DA InSAR A SEGUITO DEGLI ULTIMI ESPERIMENTI NUCLEARI IN USA NEGLI ANNI '90. Vincent et al (2003) dimostrarono 15 anni fa che con la tecnica InSAR un radar satellitare riesce a determinare le deformazioni superficiali causate da test nucleari sotterranei anche nei casi in cui l’esplosione non abbia provocato effetti visibili, per esempio la formazione di un cratere o delle fratturazioni. È evidente invece che una leggera subsidenza del terreno non sia invece facilmente rilevabile (o, quantomeno, quantificabile), con mezzi tradizionali.

La ricerca è stata effettuata nei mesi che seguirono gli ultimi esperimenti nucleari nel poligono del Nevada, alla fine del 1992. Le serie di dati InSAR hanno evidenziato la subsidenza dell’area sovrastante i siti degli esperimenti e che dopo una forte componente immediata (subsidenza cosismica), il movimento proseguiva in un’area di circa 1 km di raggio per giorni, se non per mesi o addirittura anni, dopo l’evento. La subsidenza avviene sia nel caso di formazione di un cratere sia nel caso in cui questo cratere non si formava. Il valore finale dell’abbassamento della superficie topografica va da uno a qualche centimetro.

In questa immagine, tratta da quel lavoro, si osservano i segnali della deformazione cosismica catturati in 14 mesi di osservazione di tre test sotterranei del 1992. I punti rossi rappresentano crateri di test precedenti.

La deformazione della superficie topografica
secondo Wang et al (2018)

LE DEFORMAZIONI A PUNGGYE-RI DOPO IL 3 SETTEMBRE 2017. Su Science è uscito in questi giorni (Wang et al, 2018) un articolo di un gruppo di ricerca internazionale che applica la stessa tecnica al poligono di Punggye-ri: usando i dati provenienti dal satellite tedesco TerraSAR-X, è stata misurata la deformazione causata dall’esperimento del 3 settembre, visibile nella figura qui accanto: il movimento è stato scomposto nella componente orizzontale e in quella verticale. La componente orizzontale è rappresentata dalle frecce, che mostrano una evidente dislocazione orizzontale che diverge a raggiera dal punto dell’esplosione, il cui valore arriva addirittura a 3.5 m. 

La componente verticale invece è rappresentata dai colori sulla carta.

La deformazione indica che la sorgente dell’esplosione è localizzata dalle coordinate 129.078°E e 41.300°N±50 m, ed è avvenuta a 1750 ± 100 m di quota sopra il livello medio del mare, cioè 450 ± 100 m sotto la cima del monte Mantap.

Si osserva molto bene che un buona parte della superficie del monte Mantap si è abbassata in un’area di qualche centinaio di metri dal punto dell’esplosione, mentre un movimento orizzontale ha innalzato quelle intorno; inoltre lo spostamento orizzontale è maggiore nei fianchi meridionale e occidentale della montagna, il che evidenzia uno stretto controllo della sua entità da parte della topografia: più elevata è la pendenza, maggiore è lo spostamento centrifugo. Anche l’abbassamento è legato alla pendenza, in particolare l’area a massima subsidenza si prolunga a distanza maggiore dal sito dell’esplosione nei versanti più scoscesi a W e a S. É invece curioso vedere che i valori di innalzamento siano molto superiori nella parte nord che in quella a sud.

Un aspetto interessante è che se la topografia esercita un controllo sulla entità del movimento, non lo esercita invece sulla sua direzione, per cui non si tratta di movimenti franosi, che, ovviamente, si sarebbero dovuti muovere nella direzione del pendio. Sono comunque distinguibili anche i segnali di alcune frane innescate dallo scuotimento, che sono completamente diversi e localizzati: dei debris flows localizzati in avvallamenti preestenti non in grado si produrre i moviklenti orizzontali a larga scala evidenziati dal satellite. 

Quindi la maggior parte della deformazione della superficie topografica non è il risultato delle frane ma è proprio dettato dall’esplosione e tale quadro deformativo è connesso all’espansione di una cavità dovuta all’esplosione e al suo successivo collasso: l’espansione ha provocato il rapido allontanamento dal centro delle rocce intorno all’ordigno, le quali hanno spinto le zone adiacenti che hanno reagito sollevandosi. Successivamente il collasso di quanto stava sopra alla cavità ha provocato l’abbassamento. Il raggio della cavità è calcolato sui 300 metri. Nella figura qui sotto la serie degli eventi come descritta da Wang et al (2018)

I movimenti secondo Wang et al 2018

Il cerchio bianco a sud del punto dell'esplosione indica la posizione della replica M. 4.1 avvenuta 8 minuti e mezzo dopo, localizzata circa 700 m verso sud; non a caso questo secondo evento si colloca precisamente sotto una delle zone a massima subsidenza e di movimento centrifugo nel più scosceso fianco meridionale del monte Mantap, tra il sito dell’esplosione e il portale di accesso al poligono.

In questo lavoro non viene preso in considerazione l’evento del 23 ottobre, che poi è quello che ha confermato i dubbi sulle condizioni del poligono dopo il test, che la replica M 4.1 aveva innescato. Penso che la sua posizione possa coincidere con quell’area di leggera subsidenza a NE della carta.

LA DETERMINAZIONE DELLA MAGNITUDO DELL'EVENTO CON METODO GEODETICO. Oltre alla precisa collocazione dell’evento, l’interferometria radar è capace di definire anche una Magnitudo Momento (non quella locale) dell’esperimento, combinando la profondità ricavata dalla geodesia e l’energia ricavata dai sismogrammi: il risultato è Mw = 5.5, un po' più alto di quello ricavato dai soli dati sismici (Mw = 5.24), perché include la deformazione lenta successiva all’evento che non ha prodotto onde sismiche; il cambiamento di volume è stato stimato in 0.01 km3. 

DEFORMAZIONI ED ABBANDONO DEL POLIGONO. Questo lavoro è interessante perché dimostra come ci sia un altro metodo – satellitare – che si affianca a quello geofisico per determinare sito e potenza di un esperimento nucleare. 

Le deformazioni del sito di Punggye-ri rilevate con l’interferometria radar da un lato confermano l’energia e la posizione dell’esperimento del 3 settembre ricavate dai dati geofisici, ma vanno ben oltre, perché ci danno un quadro completo della deformazione che ne è seguita, dimostrando che la chiusura del poligono non è stata volontaria ma, come già si paventava, è una conseguenza di una esplosione troppo potente per quanto potevano resistere le rocce al suo intorno. Per fortuna, trattandosi di una esplosione termonucleare a due stadi, l’ordigno a fissione capace di produrre alte ricadute radioattive era molto piccolo; mi chiedo comunque come sia possibile che in tutto questo macello non ci sia stata una fuga di radioattività.

Vincent et al (2003) New signatures of underground nuclear tests revealed by satellite radar interferometry Geophysical Research Letters, 30 / 22, 2141 
Wang et al. (2018), Science 10.1126/science.aar7230 (2018)
   

La recrudescenza dell'attività sismica nella parte settentrionale della zona interessata dalla sequenza iniziata nel 2016

Con questo post vorrei fare un pò il punto della situazione a Muccia e Pievo Torina, dove in questi ultimi mesi c'è stata una recrudescenza dell'attività sismica: si ricordano in particolare gli eventi dall'inizio di aprile, ma era da gennaio che era stato notato un aumento della sismicità nella zona a nord del bacino di Caselluccio, quella più settentrionale dell'areale interessato dalla sequenza iniziata il 24 agosto 2016, in particolare quelli con M>3. Ci si domanda poi se questa sequenza sia diversa da quelle che hanno sconvolto l'Appennino centrale nei secoli passati e, soprattutto, quanto ancora durerà e se in un futuro prossimo aree limitrofe saranno colpite. Purtroppo nessuno ha la sfera di cristallo per rispondere a queste domande e continuo a sostenere che l'unica soluzione per essere sicuri è quella di vivere, lavorare, studiare, passare il tempo libero in edifici dalle caratteristiche compatibili con il rischio sismico specifico di ogni area e che la ricostruzione delle zone terremotate deve essere preceduta da una zonazione sismica capace di individuare dove si possa ricostruire.

In questo momento la situazione nella zona di Muccia e Pieve Torina è più delicata dal punto di vista umano che da quello geofisico e sismologico: il problema fondamentale oggi è sicuramente lo stress a cui la popolazione è sottoposta da quasi due anni e ne ho avute dimostrazioni pratiche anche grazie a delle domande che mi sono state fatte personalmente. Per questo ho già perso la pazienza diverse volte contro i ciarlatani che creano allarmismi e con quelli delle postvisioni del tipo “avevo l’apparecchiatura (che peraltro non legge nulla, ndr) in manutenzione”. Il drammatico è quando, come adesso, anche esponenti del mondo scientifico si mettono a far polemica: per dirla con Niels Bohr, “è difficile fare delle previsioni, specialmente per il futuro”;  Bohr non parlava dei terremoti, ma è un fatto che proprio i terremoti siano ancora, allo stato attuale, estremamente imprevedibili.  In questo grafico vediamo in rosso il numero degli eventi con M>2 e in blu quello degli eventi con M>3. I numeri sull’asse x sono i mesi dal 24 agosto. È chiaro che i mesi non sono tutti lunghi uguale, ma insomma, si vede come nel mese successivo al 24 agosto e in quello dopo il 30 ottobre ci sono stati moltissimi eventi (e giustamente, il secondo terremoto, essendo ben più forte, ha provocato più scosse del primo); poi la sismicità è scesa a valori molto più bassi.

Il meccanismo degli eventi dell'Appennino Centrale:
una estensione 

ALCUNI CONCETTI FONDAMENTALI DA TENERE A MENTE. voglio qui riassumere alcuni aspetti che non tutti hanno ben chiari:

1. ad eccezione degli sciami sismici, in cui ci sono diverse scosse ravvicinate nel tempo con una Magnitudo simile, in una sequenza sismica normale un singolo evento principale è seguito da una serie di repliche più deboli rispetto al primo, secondo la legge di Omori per la quale in media la frequenza delle repliche diminuisce iperbolicamente con il tempo dopo un forte terremoto. Naturalmente ogni terremoto fa storia a se e quindi i coefficienti della legge di Omori variano da sequenza a sequenza e quindi, e inotre durante l’attenuazione della sequenza non si possono escludere colpi di coda un po' sopra le righe, una situazione che sta diventando molto pesante a Muccia

2. c’è poi una seconda legge, quella di Gutenberg-Richter che integra la legge di Omori: spannometricamente si può dire che il numero di repliche di un evento principale aumenta di 10 volte scendendo di una unità nel valore della Magnitudo: ad esempio un evento M5 sarà seguito da 10 eventi di M 4, 100 di M3, 1000 di M 2 e via discorrendo, fermo restando che sotto un certo limite che mettiamo convenzionalmente a M = 1 gli eventi sismici diventano indistinguibili dal rumore di fondo

Il meccanismo compressivo tipico
dei terremoti emiliani del 2012

3. tra la sismicità dovuta a un quadro tettonico compressivo e la sismicità dovuta a un quadro tettonico distensivo  c’è un abisso nella modalità dei meccanismi focali e in quello che succede dopo: possiamo dire che in un evento compressivo i due lati del piano di faglia si “uniscono” sempre di più; in un evento distensivo invece si crea un vuoto

4. è per questo che la sequenza che segue un evento distensivo dura molto di più rispetto ad una crisi dovuta ad un evento compressivo. Lo possiamo facilmente notare confrontando la durata di pochi mesi delle repliche dei terremoti emiliani del 2012 con quella dell’Appennino centrale attuale; ma anche con quella, distensiva come quella attuale, del 1980 in Irpinia, che ancora nel 1982 ha visto alcuni eventi con M > 4 nella zona epicentrale 

5. è indubbio che la legge di Omori sia rispettata: nell’area interessata dagli eventi principali del 2016 le scosse stanno diminuendo in intensità e frequenza praticamente dappertutto, anche se nella zona di Muccia e di Pieve Torina c’è, come dire, una recrudescenza dei fenomeni dall’inizio del 2018 e in questi giorni è in corso la sequenza di Omori al seguito dell’evento del 10 aprile (che per la cronaca ha un meccansismo focale e una direzione del piano di faglia perfettamente coerente con le scosse principali dell’ottobre 2016).

MAGNITUDO, DISTRIBUZIONE E NUMERO DEGLI EVENTI. Osservando esclusivamente gli eventi con M uguale o superiore a 3, la sequenza iniziata il 24 agosto 2016 è caratterizzata da 3 fasi, come si osserva nelle tre carte qui sotto.

Qui sopra, nella prima fase, dal 24 agosto al 25 ottobre, la sismicità è essenzialmente annidata fra Preci a nord e il lago di Campotosto a sud, ma con gli eventi più meridionali “staccati” dagli altri. Questo denota che il settore del Monte Vettore ancora non ha iniziato a muoversi.

Qui vediamo invece come il 25 ottobre 2016 si attiva anche la parte settentrionale tra Preci e Camerino, e fino al 15 gennaio 2017 le replice coprono tutta l'area tra Muccia a nord e Pizzoli a sud, a parte una leggera discontinuità che separa gli eventi della zona di Amatrice da quelli di Campotosto. 

Dopo il 15 gennaio 2017 non esiste più il gap tra Amatrice e Campotosto e la distribuzione delle repliche è praticamente continua.

Quello che si può notare è che applicando la legge di Gutemberg – Richter i conti ancora non tornano: ci dovrebbero essere stati, dopo il 30 ottobre, 10 eventi con M compresa tra 5 e 6, 100 con M compresa fra 4 e 5 e 1000 con M compresa fra 3 e 4. Invece siamo a 5 (e tutti fuori zona, il 30 ottobre è stata interessata l’area più a nord, e queste repliche intense sono annidati a Campotosto, esattamente dalla parte più a sud…), 65 e poco più di 1000 (il valore di 1132 per questioni grafiche comprende anche gli eventi dello spoletino e qualcun altro fuori dalla zona “calda”). Insomma, ce ne sono stati meno nelle classi di Magnitudo superiore (la metà di quelli aspettati) e un numero abbastanza compatibile in quella di M tra 3 e 3.9.

Dal punto di vista temporale in tutta l'area interessata gli eventi complessivi con M3+ sono 260 entro il 23 ottobre 2016, e 957 entro il 14 gennaio 2017; consoderando solo la zona del M 6.5 del 30 ottobre 2016 i 3 eventi più forti che hanno seguito la scossa principale sono:

M 4.8 1 novembre 2016

M 4.7 3 novembre 2016

M 4.6 10 aprile 2018

Naturalmente questo non vuole dire che ci saranno per forza altri eventi forti per colmare il gap nelle magnitudo più alte secondo la legge di Gutemberg -  Richter, ma il rischio esiste.

UN CONFRONTO FRA OGGI E IL PASSATO. La domanda fondamentale è se precedenti terremoti di simile entità nell’Appennino Centrale abbiano avuto la stessa evoluzione di quella attuale. In parte, consultando e cronache storiche, la risposta può essere considerata “affermativa”, perché numerose fonti attestano la presenza di repliche prolungate nel tempo; dobbiamo inoltre notare come molti terremoti siano incerti non solo nella magnitudo, ma persino nell'essere avvenuti o meno; questo vale sia per gli eventi sismici principali, ma soprattutto per le repliche. Spesso si legge che le scosse sono durate mesi (o anni), ma non è dato sapere quanti, anche perché in genere gli eventi che restano sono quelli distruttivi e l’evento principale spesso distruggeva praticamente tutto quello che c’era da distruggere: in genere le repliche sono riportate dalle cronache solo se particolarmente avvertibili e/o causa di frane o ulteriori rovine negli edifici. E sicuramente una loro buona parte non è stata citata. Ne segue che il catalogo è sicuramente lacunoso.

Il post 2016 si differenzia dalle precedenti sequenze per alcune condizioni particolari dal punto di vista umano:

• le cronache del passato non possono citare le scosse strumentali, e, come ho detto, sono abbastanza lacunose anche su eventi sismici “molto sopra le righe” (quando, appunto, non hanno provocato effetti pratici).
• fino al 1915 non ci sono registrazioni sismometriche e, come fa notare Alessandro Amato, mai nella storia italiana la rete sismografica è stata particolarmente sofisticata nel numero e nella qualità dei sensori (una condizione ben più dettagliata non solo rispetto al 1915, ma anche rispetto al 1997 e al 2009)
• la Magnitudo decisamente più alta rispetto agli eventi principali recenti (1997 e 2009)

Quindi non è possibile dire se questa sequenza si differenzia da altre di Magnitudo simili documentate dalle fonti storiche. Tantomeno è possibile prevedere quando il tutto si esaurirà.

MUCCIA, PIEVE TORINA E DINTORNI. Veniamo ora ai possibili scenari futuri, argomento del contendere di questi giorni con la parola d’ordine “cosa succede a Muccia e a Pieve Torina?”. La persistenza dell’attività nella zona settentrionale è probabilmente dovuta alla maggiore intensità del terremoto del 30 ottobre nei confronti di quello del 24 agosto, ma è innegabile che dopo una fase di stasi nella seconda metà del 2017, la “coda” si sia nuovamente intensificata. Non chiedetemi però il perché….

Venendo agli utlimi 2 mesi, dal 4 aprile ci sono stati fino alle 10.00 del 23 maggio 2018, nell’area molto localizzata evidenziata nella carta, 172 eventi con M > 2.1, dei quali 19 con M 3+ con le punte di M 4.0 il 4 aprile e 4.6 il 10 aprile. Non solo: in quest’area dall'inizio del 2018 si registra un’attività superiore a quella dei mesi precedenti, come si vede dal grafico qui accanto: da gennaio 2018 ben 36 dei 38 eventi con M >3 di tutta l'area tra Muccia e Pozzoli sono stati registrati qui.

L’evento del 4, ma soprattuto quello del 10, hanno innescato a loro volta una sequenza di Omori e si vede che in questo caso la legge di Gutemberg-Richter sottostima un po' il numero delle repliche.

Dobbiamo poi registrare una leggera ripresa negli ultimi giorni, dopo un paio di settimane fondamentalmente calme; in particolare un evento del 21 maggio è stato il più forte dopo il 10 aprile.  Per fare un confronto con le sequenze storiche ci si dovrebbe domandare quanti di questi eventi non solo sarebbero stati avvertiti dalla popolazione, ma soprattutto di quanti resterebbe nelle cornache il ricordo. Probabilmente ben pochi… senza la rete di INGV la maggior parte sarebbe passata inosservata e, se fosse successo anche solo 200 anni fa, le cronache storiche probabilmente non avrebbero registrato il ricordo dell’evento del 10 aprile 2018.

I POSSIBILI SCENARI FUTURI. Purtroppo, come ho detto, non si possono fare previsioni, ma è possibile delineare degli scenari futuri, a scala locale (la zona interessata dal 24 agosto 2016 in poi) e a scala regionale (l’Appennino centrale).

A scala locale lo scenario ottimale dal punto di vista umano sarebbe “finisce tutto subito”. Purtroppo non è realistico, nel senso che l’attività della sequenza sismica innescatasi nel 2016 sta ancora continuando. In questo quadro di generale diminuzione della frequenza e dell’intensità delle repliche sono purtroppo possibili  ancora “colpi di coda” come questo di Muccia, con il dubbio  sulla scarsa quantità delle repliche di M più alta e del perché della ripresa dopo una stasi. Insomma.. ancora “non ci siamo”.

La distribuzione irregolare nel tempo degli
eventi sismici maggiori dell'appennino centrale
 è chiaramente evidente in questo grafico
tratto da Tondi e Cello (2003)

A scala dell’Appennino centrale, torniamo al gennaio 2017 e al comunicato della Commissione Grandi Rischi, che ho commentato qui: gli eventi sismici nell’Appennino centrale non sono distribuiti casualmente nel tempo, ma tendono a raggrupparsi in cluster temporali ben definiti.

Esaminando il “catalogo parametrico dei terremoti italiani” dell’INGV, si notano dei momenti in cui il settore a cavallo fra Lazio, Umbria e Marche è stato colpito da una serie di eventi con M superiore a 5.5 ravvicinati nel tempo: limitandosi al periodo tra il XIII e il XVIII secolo, è successo per esempio tra il 1269 e il 1279 e negli anni 1348 – 49, e tra il 1458 e il 1466, mentre in tutto il 1500 si segnala una attività scarsissima e nel 1600 si contano “appena” 3 eventi maggiori, tutti in un lasso temporale ristretto; nella prima metà del XVIII secolo, invece, si sono verificate scosse intense circa ogni 10 anni (1703, 1719, 1730, 1741, 1747 e 1751), poi dopo qualche decennio di calma abbiamo terremoti nel 1781, 1785, 1791 e 1799. Vediamo la distribuzione nel tempo dei terremoti più importanti in questa figura tratta da Tondi e Cello (2003).

Comunque, ci sono anche degli eventi “isolati”, come nel 1298, 1328, 1599, per cui l'asserzione "se abbiamo un evento forte, allora ne verranno altri" è stata per fortuna varie volte contraddetta.  Adesso abbiamo avuto in 20 anni il 1997, il 2009 e il 2016 (più qualche alto evento un po' fortino). Quindi è indubbio che siamo di fronte ad una crisi sismica importante. Crisi che potrebbe essere conclusa con solo queste 3 sequenze principali, o forse 4 se gli eventi del 24 agosto e del 30 ottobre possano essere considerati distinti oppure appartenenti alla stessa sequenza. Comunque il problema è che attualmente la pericolosità sismica nelle aree dell’Appennino Centrale adiacenti a quelle interessate dalle sequenze dgli ulimi 20 anni è teoricamente più elevato di quello che c’era 50 anni fa.

Le strutture sismogenetiche principali sono suddivise in due sistemi grossomodo paralleli, uno orientale che va dal Vettore alla Maiella e uno occidentale che va da Colfiorito all’alta valle del Sangro (Galadini e Galli, 2000). C’è poi un terzo allineamento, ancora più occidentale, non compreso nella carta  e obliquo rispetto ai due precedenti, che dalla Valtiberina e dalla valle Umbra arriva a Leonessa e alla valle del Salto unendosi a quello occidentale nell’area del Fucino (Boncio et al, 2004). Ne ho parlato in dettaglio qui.

Queste faglie accomodano la deformazione provocata dalla estensione (di oltre 1,5 mm /anno) che sta subendo la crosta, perché il settore adriatico si muove verso Est più di quanto non lo faccia il settore tirrenico (Farolfi e Delventisette 2016, (ne ho parlato qui). Dei due sistemi, negli ultimi secoli quello occidentale è stato più attivo di quello orientale e gli eventi del 1979 (Valnerina), 1997 e 2009 appartengono a questo, come quasi tutti gli altri forti terrmeoti storici, mentre il sistema orientale ha originato i terremoti del 2016, e quelli del Fabrianese e quelli di Amatrice e Laga del XVII secolo.

Secondo la Commissione Grandi Rischi i settori più a rischio attualmente sono quello di Norcia, che fa parte del sistema occidentale, perché si trova in mezzo fra le zone interessate nel 1997 e nel 2009 ed è di fronte al settore attualmente interessato dalla sequenza, quello della Laga, perche è quello immediatamente a sud del settore del Vettore (ma forse gli eventi di Campotosto sono già una risposta) e l’Alta Valle d’Esino, fino a Fabriano, dove il terremoto Mw 6.2 del 24 aprile 1741, che seguì e precedette il gli eventi del 1747 e 1751 della zona di Gualdo Tadino.

Ripeto nuovamente, comunque, il concetto: NON non si tratta di una previsione ma si tratta di una valutazione scientifica di possibili scenari. Le previsioni, checchè ne pensino i soliti personaggi, NON sono attualmente possibili. Punto e basta.

Dal punto di vista pratico, la cosa principale che emerge dallo scenario peggiore è che per la ricostruzione si dovrà tenero conto del rischio sismico elevato e che sono numeriosi gli edifici nelle zone vicine che vanno urgentemente adeguati.

Boncio et al 2004. Defining a model of 3D seismogenic sources for seismic hazard assessment applications: the case of central Apennines (Italy). J. Seismol. 8, 407–425

Galadini e Galli (2000) Active tectonics in the central Apennines (Italy) — input data for seismic hazard assessment. Natural Hazards 22: 225–270

Tondi e Cello (2003) Spatiotemporal evolution of the Central Apennines fault system (Italy) Journal of Geodynamics 36, 113–128

Estinzione del Devoniano superiore e Large Igneous Province della Yacuzia

L’evento Kellwasser è una delle “Big 5” di Sepkosky, cioè uno dei principali eventi di estinzione di massa. È avvenuto nel Devoniano superiore, ma non alla fine del periodo, che peraltro corrisponde ad un’altra fase di estinzione, nota come Evento Hangenberg. Sulle sue cause si è discusso parecchio, anche per la presenza in quel tempo di un impatto meteoritico; oggi appare abbastanza ovvio che si tratti delle conseguenze della messa in posto della Large Igneous Province dei basalti della Yacuzia e anche le ultime ricerche puntano il dito contro il vulcanismo: durante il Kellwasser notiamo una forte anomalia positiva di mercurio nei sedimenti, che si accompagna a quella dell’iridio, già nota da tempo. Il problema è che le datazioni dei basalti della Yacuzia tornano poco a causa del pessimo stato di conservazione di quella serie vulcanica,É piuttosto interessante notare una circostanza fondamentale nella storia della vita e cioè la comparsa dei primi vertebrati terrestri proprio nel tempo che intercorre fra il Kellwasser e l’ Hangenberg.

Un tipico pesce senza mascelle del Devoniano

Nel Devoniano superiore ci sono stati due diversi eventi di estinzione di massa  a poca distanza l’uno dall’altro. Il primo al limite Frasniano – Famenniano (FFB, conosciuto come evento Kellwasser) e che corrisponde ad uno dei “big five”, i 5 maggiori eventi di estinzione di massa della storia (Sepkosky, 1996); il secondo è di poco successivo e corrisponde alla fine del Famenniano (evento Hangenberg) e siccome il Famenniano è l’ultimo stadio del Devoniano, il secondo evento corrisponde anche al limite Devoniano – Carbonifero. Quindi come succederà in seguito nel Permiano superiore, abbiamo due estinzioni di massa separate da poco più di 10 milioni di anni; la differenza fra le due coppie è che nel Devoniano il primo è più importante del secondo, mentre nel Permiano, l’evento più importante è il secondo.  

Venendo specificamente al Kellwasser, si tratta di una estinzione di massa che ha coinvolto, come sempre succede, una vasta gamma di forme viventi appartenenti a gruppi molto diversi: coralli, stromatoporoidi (animali all’epoca molto diffusi, ormai estinti e ancora di incerta classificazione), trilobiti, cefalopodi, brachiopodi, conodonti e pesci (in particolare quelli senza mascelle). L’estinzione è stata veloce, ma graduale (Becker, 1993); fondamentalmente si è trattato di una sostituzione progressiva della fauna precedente con organismi sempre più adatti a vivere in condizioni di scarsa ossigenazione delle acque (guarda caso come al K/T…). 

Le cause dell’evento Kellwasser sono state variamente addebitate: cambiamenti nella circolazione oceanica, provocati alternativamente da riscaldamento o raffreddamento climatico globali (o altro non specificato), cambiamenti del livello marino, impatto di un meteorite (immancabile...).

Il cratere di Siljan, in Svezia.
 Probabilmente precede di poco l'evento Kellwasser

Al di là delle caratteristiche biotiche, ci sono molte circostanze geologiche interessanti, che poi sono comuni anche agli altri eventi di estinzione di massa (in particolare al K/T):

• la più evidente è la deposizione dei classici sedimenti finissimi e scuri perché ricchi di materia organica, sintomo di condizioni di mancanza (o quantomeno di forte scarsità) di ossigeno. Le serie meglio studiate sono quelle dello stato di New York, dove troviamo più livelli del genere, che precedono e seguono quello principale posto al limite fra Frasniano e Famenniano. Questi livelli assomigliano al Fiskeler, il famoso sedimento che troviamo in Danimarca depositato esattamente al limite K/T (Christensen et al 1973). Come in esso troviamo cristalli di pirite e laminazioni sottilissime che possono restare evidenti solo in caso di assenza di vita nel fondo marino, perchè altrimenti questa finissima stratificazione verrebbe completamente distrutta dagli animali che vivono sul fondo e, soprattutto, da quelli che lo scavano
• un forte disturbo, a livello globale, del rapporto isotopico del Carbonio, il δ13C: fra gli eventi di estinzione globale alcuni presentano alternativamente oscillazioni positive o negative di questo rapporto; in questo caso la situazione è un pò complessa del solito, in quanto al Kellwasser troviamo sia oscillazioni positive che negative a seconda della regione
• forti oscillazioni del livello marino; in particolare una fase a basso livello coincide esattamente con la parte terminale del Frasniano, documentata da interruzioni della sedimentazione in molte aree (Nordamerica, Europa, Iran etc etc)
• una anomalia positiva nella concentrazione di Iridio nelle argilliti scure, con livelli molto elevati anche se non tanto quanto al K/T (Over et al 1997) e, è arrivata la notizia in questi giorni, anche di Mercurio (Racki et al, 2018)
• in Belgio troviamo anche delle microtectiti, tipiche di un impatto meteoritico (Claeys and Casier, 1994) e non è una coincidenza casuale in quanto nei dintorni del Kellwasser si colloca un grande impatto meteoritico, quello di Siljian, in Svezia. Si tratta di un cratere piuttosto grande (il diametro  varia a seconda degli Autori fra 65 e 90 km)
• sempre in Belgio sono stati anche trovati alti valori di molecole organiche, che indicano diffusi incendi di foreste ed anche evidenze di una accelerazione dell’erosione del suolo (Kaiho et al, 2013)
• e, soprattutto, la messa in posto di una Grande Provincia Magmatica (LIP, Large Igneous Province), i basalti della Yacuzia 

La stratigrafia dei sedimenti del Frasniano nello Stato di New York, che evidenzia
i numerosi episodi di sedimentazione di black shales con anomalia dell'Iridio (Over et al, 1997)

L’IRIDIO: voglio soffermarmi in particolare sulla questione dell’Iridio. La concentrazione di elementi del gruppo del Platino (PGE) nei sedimenti dello stato di New York è stata oggetto di un dettagliato studio (Over et al 1997), dove vari orizzonti di scisti neri, sia del Fammenniano superiore che quello corrispondente propriamente al Kellwasser, contengono concentrazioni anomale di Iridio; i rapporti quantitativi fra gli elementi PGE non sempre sono simili a quelli delle condriti. I campioni presi nei sedimenti “normali” invece non presentano arricchimento in iridio, platino e rutenio. Inoltre, come succede al K/T per le tectiti dello Yucatan, i sedimenti dell’impatto meteoritico di Siljan trovati in Belgio non presentano anomalia dell’iridio. Insomma, sia i rapporti quantitativi fra gli elementi del gruppo del platino, sia la presenza di diversi orizzonti arricchiti separati nel tempo non indicano una singola fonte extraterrestre, come quella rappresentata da un impatto. Anche nell’evento della fine del Devoniano troviamo una anomalia dell’Iridio in Cina (Wang et al, 1993) e in Nordafrica (Kaiser et al, 2011)

DATAZIONE DELL’IMPATTO DI SILJAN. Sull’onda dell’ipotesi del meteorite killer dei dinosauri (connessone che ormai è stata dimostrata falsa, ma che continua ad essere molto popolare) sono state ricercate connessioni fra l’impatto di Siljan e il Kellwasser. Il problema è ancora aperto in quanto la forbice cronologica degli eventi è piuttosto grande e, a complicare le cose, ci si mette lo stesso limite Frasniano – Famenniano, che negli ultimi 40 anni è oscillato tra 370 e 376 Ma; nell’ultima versione della carta cronostratigrafica ufficiale (2017) è posto a 372,2 ± 1.6 milioni di anni, quindi fra 370,6 e 373,8 Ma. Quanto al cratere, le età ricavate dalla bibliografia sono diverse, ma in sostanza tendono ad essere un pò più vecchie dell’intervallo in cui si colloca il Kellwasser; quelle ricavate più di recente che ho trovato parlano di 380.9 ± 4.6 Ma (Jourdan et al., 2012). Precedentemente due campioni avevano fornito un’età rispettivamente di 377.2 ± 2.5 Ma e 376.1 ± 2.8 (Reimold et al. 2005) per cui una connessione diretta fra i due eventi è piuttosto difficile.

Carta della giounzione tripla della Yacuzia
da Ricci et al (2013)

VULCANISMO E EVENTO KELLWASSER. I fenomeni registrati alla fine del Frasniano sono sintomi della presenza di un vulcanismo importante. Ma, come ho accennato, proprio in questi giorni ne è venuto fuori uno ancora più stringente: nei sedimenti deposti al FFB c’è una vistosa anomalia di mercurio (Racki et al 2018), anomalia che si riscontra anche in corrispondenza di altri eventi di estinzione di massa, come alla fine di Permiano, Triassico e Cretaceo (Bergquist, 2017). 

Il vulcanismo nel Frasniano è stato molto intenso, per esempio negli Stati Uniti ed in Europa, dove si stava perfezionandosi la collisione fra Euromerica e Gondwana, in Sudan e in Cina. Ci sono inoltre delle grandi eruzioni basaltiche nel cratone est europeo:  i basalti di Kola e quelli del rift di Pripyat–Dniepr – Donets tra Ucraina e Bielorussia (Kravchinsky, 2012).

Ma, soprattutto, in quei tempi avveniva la messa in posto di una delle più massicce coperture basaltiche degli ultimi 500 milioni di anni, la Large Igneous Province della Yacuzia, nota anche come Trappi di Viluy. 

I TRAPPI DELLA YAKUZIA. Come dice il nome della regione, siamo nella parte NE della Siberia, in una regione venuta alla ribalta con il Risiko, dove nel Devoniano superiore si è formata una giunzione tripla con la formazione di 3 rift (Ernst and Buchan, 1997), in una situazione che ricorda molto quella attuale dell’Afar. Dei tre rami del rift, due formano gli attuali margini della piattaforma siberiana ma sono abortiti. Il terzo invece ha avuto una storia diversa: dopo una fase di divergenza che ha formato un bacino oceanico il regime è tornato presto compressivo e il conseguente scontro fra la Siberia e il microcontinente di Kolyma -  Omlon ha formato nel mesozoico l’orogene dei monti Verkhoyansk (Oxman, 2003).

Nel rift di Viluy, il ramo abortito di SW, che attualmente è lungo 800 km e largo fino a 450 e in altre aree intorno alla giunzione tripla troviamo la Large Igneous Province della Yacuzia, di cui sono visibili alcuni resti lungo i lati dei rift, perchè la maggior parte è stata coperta dai sedimenti che si sono accumulati nei rift stessi; le stime parlano di quasi 2 milioni di km cubi fra lave e il loro sistema di alimentazione (Kiselev et al., 2006). 

La domanda a questo punto è se proprio la LIP di Yacuzia sia il killer del Kellwasser. Diciamo che la cosa avrebbe senso, dato che anche gli altri eventi di estinzione di massa, con solo probabilmente una eccezione, il passaggio Eocene – Oligocene, sono associati a Large Igneous Provinces, e agli altri fenomeni descritti all’inizio del post, come si vede dall’immagine. 

DATAZIONE DEI TRAPPI DELLA YACUZIA. Le indicazioni stratigrafiche indicavano il Devoniano superiore come età di queste rocce ignee, ma le datazioni assolute erano poche, e poco affidabili, anche perchè  la maggior parte di queste rocce sono molto alterate. Nel 2010 con il metodo 40Ar/39Ar furono trovate due età di 360.3 ± 0.9 and 370.0 ± 0.7 Ma, il che ha suggerito che due impulsi diversi di questi basalti abbiano provocato sia la grande estinzione al liimte Frasniano – Famenniano che quella successiva alla fine del Devoniano (l’ Hangenberg) (Courtillot et al, 2010). 

Più recentemente, anche Ricci et al (2015) esaminando diversi campioni hanno notato che con il metodo 40Ar/39Ar questi si raggruppano in due gruppi diversi con età medie ponderate di 364,4 ± 1,7 e 376,7 ± 1,7 Ma. 

Il quadro ricavato supporta quindi in entrambi gli studi la presenza di due impulsi principali (come succede spesso nelle LIP) e tutto sommato l’intervallo che si ricava è lungo più o meno come il Famenniano. Per cui è possibile (ed intrigante) che sia il Kellwasse che l’Hangenberg siano dovuti a due diverse fasi parossistiche dei Trappi di Viluy.  

Il problema è che soprattutto a causa della alterazione (che ha probabilmente interessato il sistema dell’argon) questi risultati mi paiono un pò aleatori e poco in accordo con la cronologia attualmente in vigore: insomma, oggi le stime per il Kellwasser vanno da 373,8 a 370,6 MA, quelle di Ricci e soci per la prima fase dei trappi di Viluy vanno da 374,0 a 378,4 MA. Quindi le due età si toccano quasi. Per quanto riguarda invece l’Hangenberg, è un limite abbastanza ben definito dal punto di vista cronologico, in quanto ha una incertezza nella sua datazione assoluta di poche centinaia di migliaia di anni: al massimo può essere avvenuto a 359.3, ben dopo il limite più recente della forbice cronologica della Yacuzia, che è post a 362,7. 

Insomma, il fatto che in entrambi i casi le vulcaniti risultino più vecchie degli eventi biotici può essere dovuta alla pessima qualità dei campioni dal punto di vista della loro alterazione, ma se la correlazione fra eruzioni della Yacuzia e Kellwasser è “estremamente probabile” dal punto di vista della cronologia relativa, le datazioni assolute sono ancora poco chiare. Invece sull’Hangenberg i dubbi sono ancora estremamente alti. La tabella qui accanto mostra le attuali incertezze nelle datazioni assolute.

In conclusione si può dire che la geochimica delle serie sedimentarie evidenzia sicuramente al limite Frasniano – Famenniano degli importanti fenomeni vulcanici, e i vari episodi conosciuti di quell’epoca non paiono dal punto di vista quantitativo essere minimamente paragonabili a quelli della Siberia orientale. Per cui è piuttosto verosimile che il killer della fine del Frasniano sia ancora una volta una LIP, quella della Yacuzia, anche se i meccanismi che portano una LIP ad essere uno spietato killer non sono ancora chiarissimi e la datazione assoluta delle eruzioni è ancora incerta. 

Becker (1993) Anoxia, eustatic changes, and Upper Devonian to lowermost Carboniferous global ammonoid diversity. In: M.R. House (Editor), The Ammonoidea. Environment, Ecology, and Evolutionary Change. Syst. Assoc. Spec. 47, 115–163
Bergquist (2017) Mercury, volcanism, and mass extinctions: PNAS 114, 8675–8677 
Courtillot et al (2010) Preliminary dating of the Viluy traps (Eastern Siberia): eruption at the time of Late Devonian extinction events? Earth and Planetary Science Letters 300, 239–245. 
Christensen et al (1973) Sedimentology and depositional environment of Lower Danian fish clay from Stevns Klint, Denmark. Bull. Geol. Soc. Denmark, 22, 193–212
Ernst e Buchan (1997) Giant radiating dyke swarms: their use in identifying pre-Mesozoic large igneous provinces and mantle plumes. In: Mahoney, J., Coffin, M. (Eds.), Large Igneous Provinces: Continental, Oceanic, and Planetary Volcanism: Am. Geophys. U. Geophysical Monograph Series, 100, 297–333. 
Golonka and Gawęda (2012) Plate Tectonic Evolution of the Southern Margin of Laurussia in the Paleozoic, Tectonics Evgenii Sharkov, IntechOpen
Kaiho et al (2013) A forest fire and soil erosion event during the Late Devonian mass extinction Palaeogeography, Palaeoclimatology, Palaeoecology 392, 272–280
Kaiser et al (2011) Climate-controlled mass extinctions, facies, and sea-level changes around the Devonian–Carboniferous boundary in the eastern Anti-Atlas (SE Morocco). Palaeogeography, Palaeoclimatology, Palaeoecology 310, 340–364 
Kiselev et al (2006) Middle Paleozoic basic magmatism of the northwestern Vilyui Rift: composition, sources, and geodynamics. Petrology 14, 588–608
Kravchinsky (2012) Paleozoic large igneous provinces of Northern Eurasia: Correlation with mass extinction events Global and Planetary Change 86-87, 31–36
Over et al (1997) Platinum group element enrichments and possible chondritic Ru:Ir across the Frasnian–Famennian boundary, western New York State Palaeogeography, Palaeoclimatology, Palaeoecology 132,399–410 
Over (2002) The Frasnian/Famennian boundary in central and eastern United States Paleogeography, Palaeoclimatology, Palaeoecology 181,153-169 
Oxman (2003) Tectonic evolution of the Mesozoic Verkhoyansk–Kolyma belt (NE Asia) Tectonophysics 365,45–76
Racki et al (2018) Mercury enrichments and the Frasnian-Famennian biotic crisis: A volcanic trigger proved? Geology  DOI: https://doi.org/10.1130/G40233.1  
Ricci et al 2013 New 40Ar/39Ar and K–Ar ages of the Viluy traps (Eastern Siberia): Further evidence for a relationship with the Frasnian–Famennian mass extinction Palaeogeography, Palaeoclimatology, Palaeoecology 386 (2013) 531–540 
Sepkoski JJ (1996) Patterns of Phanerozoic extinction: A perspective from global databases. Global Events and Event Stratigraphy in the Phanerozoic, ed Walliser OH (Springer, Berlin), pp 35–51
Wang et al (1993) Global Iridium anomaly, mass extinctions and redox change at the Devonian-Carboniferous boundary. Geology 21, 1071 – 1074