Jedi Simon Linkshaa Solution
 
Jedi Simon Foundation
 Resonant Vibration Department Research - Cymatics Research -

Approccio Sinergico Integrato: Linkshaa Science
per risolvere quello che una sola scienza non puo',
al servizio della vita. Multivariable Earthquake precursors.
Si applichi nel caso di vita o di morte il Diritto condiviso della conoscenza, al servizio dell'uminatià.
Nel nome della conoscenza, della consapevolezza nell'interesse della comunità.

L’Aquila 30 Novembre 2006 Cronistoria dell’attività di Ricerca
della S C S sui “Precursori di eventi sismici”.Anni 2000-2001

Un piccolo gruppo di Fisici e Tecnici che lavorano presso i Laboratori Nazionali del Gran Sasso
(CNR-INFN) decidono, per comune interesse e curiosità scientifica, di iniziare una ricerca basata su
un nuovo metodo di misura delle emissioni di 222Radon.
Questa ricerca e’ interamente privata e viene fatta con strumentazione riadattata e installata in uno
scantinato di una civile abitazione alla periferia dell’Aquila.

Le persone coinvolte sono:
Sig. Gioacchino Giuliani – CTER - IFSI - LNGS
Sig. Roberto Giuliani – Tecnico Informatico INFN – LNGS
Prof. Victor Alekseenko - Ricercatore del Baksan Neutrino Astronomy (Russia)
L’interesse per questa ricerca, era rappresentato dalla possibilità di poter determinare con certezza, se
le variazioni di concentrazione di Radon emesse dalla crosta terrestre, in prossimità di un terremoto,
avvenissero prima dell’evento, durante o dopo.
Va segnalato che ricerche su emissioni di 222Radon mirate alla possibilità di osservare Precursori
sismici, sono state effettuate nel corso degli ultimi 40 anni, da ricercatori e università di tutto il
mondo, ma senza ottenere mai risultati significativi.
A tutt’oggi ricerche di questo tipo sono ancora in corso.
Anni 2001-2002

Il primo periodo di lavoro, è stato speso per realizzare un sismografo a pendolo verticale.
Essenzialmente il sismografo avrebbe fornito informazioni sulla dinamicità nell’area di sensibilità del
rivelatore di Radon, e tracciare una mappa sismogenetica locale, per ottenere un rapporto di
correlazione tra flusso medio di Radon e l’attività dinamica prodotta dagli eventi sismici osservati.
L’importanza di un sismografo per le correlazioni delle variazioni di concentrazione di Radon ed
eventi sismici, risulta essenziale per la determinazione della sismogeneticità della zona e la media
statistica di eventi per anno, per stagione, per mese.
 

L’immagine a sinistra mostra il sismografo a pendolo verticale.
Esso è stato posizionato a circa tre metri di distanza dal prototipo rivelatore
di Radon. Entrambi i sensori, sono tutt’ora localizzati in una stanza interrata
a circa tre metri di profondità, in assenza di ventilazione.

Il secondo periodo di attività è stato dedicato alla realizzazione di diversi tipi di contatori, per
ricercare un metodo ed una tecnica attendibili.
E’ naturale, credo, che dovendo trovare correlazioni tra un elemento ed un fenomeno tra loro
associabili, si puntualizzi lo studio o sull’elemento o sul fenomeno.
Parve quindi naturale, realizzare degli strumenti artigianali, per il rilevamento di particelle alfa
prodotte dal decadimento del Radon attraverso delle semplici camere a ionizzazione o tubi
proporzionali.

Ricordiamo, per i meno informati, alcuni elementi caratteristici del Radon.
Esso è un elemento mobile, radioattivo e chimicamente inerte, questa doppia proprietà, di essere
elemento estremamente mobile e chimicamente non reattivo, gli danno l’appellativo di elemento
nobile.

In natura si trovano tre isotopi del Radon, il 222Rn prodotto dall’ 238U, il 220Rn prodotto dal 232Th ed il
219Rn prodotto dal 235U.
Dalle diverse catene di decadimento radioattivo, 238U, 232Th e 235U, tutti gli isotopi del Radon vengono
prodotti per emissione di particella alfa dall’elemento 226Radio.
Sono state così realizzati due tipi di contatori, uno a camera di ionizzazione ed uno a tubi, così come
mostrano le immagini seguenti.

L’immagine a destra mostra la sezione digitale del sismografo. Essa è costituita dal sistema di taratura, dal datario e dalla
stampante termica per i sismogrammi. E’ tra l’altro prevista la possibilità di inviare il segnale sismico direttamente su un Pc,
che permette di avere un allarme sonoro in coincidenza dell’evento sismico.


L’ immagine a sinistra mostra 2 camere a ionizzazione, costituita da due piastre contrapposte, (Catodo), attraversate nella parte centrale e lungo la
diagonale da un cavo di circa 40 μm (Anodo) sul quale viene applicata una tensione di circa 1700 V.
Il campo elettrico prodotto all’interno della camera, permette di osservare la cascata elettronica prodotta dal passaggio degli ioni di Radon e dei gas presenti
nell’ambiente.

Il secondo periodo di attività è stato dedicato alla realizzazione di diversi tipi di contatori, per
ricercare un metodo ed una tecnica attendibili.
E’ naturale, credo, che dovendo trovare correlazioni tra un elemento ed un fenomeno tra loro
associabili, si puntualizzi lo studio o sull’elemento o sul fenomeno.
Parve quindi naturale, realizzare degli strumenti artigianali, per il rilevamento di particelle alfa
prodotte dal decadimento del Radon attraverso delle semplici camere a ionizzazione o tubi
proporzionali.

Ricordiamo, per i meno informati, alcuni elementi caratteristici del Radon.
Esso è un elemento mobile, radioattivo e chimicamente inerte, questa doppia proprietà, di essere
elemento estremamente mobile e chimicamente non reattivo, gli danno l’appellativo di elemento
nobile.

In natura si trovano tre isotopi del Radon, il 222Rn prodotto dall’ 238U, il 220Rn prodotto dal 232Th ed il
219Rn prodotto dal 235U.
Dalle diverse catene di decadimento radioattivo, 238U, 232Th e 235U, tutti gli isotopi del Radon vengono
prodotti per emissione di particella alfa dall’elemento 226Radio.
Sono state così realizzati due tipi di contatori, uno a camera di ionizzazione ed uno a tubi, così come
mostrano le immagini seguenti.
 

L’ immagine in alto mostra 2 camere a ionizzazione, costituita da due piastre contrapposte, (Catodo), attraversate nella parte centrale e lungo la
diagonale da un cavo di circa 40 μm (Anodo) sul quale viene applicata una tensione di circa 1700 V.
Il campo elettrico prodotto all’interno della camera, permette di osservare la cascata elettronica prodotta
dal passaggio degli ioni di Radon e dei gas presenti nell’ambiente.
 

L’immagine in basso mostra invece un rivelatore costituita da 3 tubi, il cui funzionamento si basa sullo
stesso principio delle camere precedenti. Anche in questo caso l’interno del cilindro è attraversato da un cavo di 40 μm,
sul quale è applicata una tensione continua negativa, pari a circa 1400 – 1600 V.

La figura in alto mostra lo schema elettrico utilizzato sia per le camere a ionizzazione che per i tubi.

Il sistema sopra descritto, dopo un paio di mesi di test, si è rivelato poco pratico ed anche piuttosto
impreciso per le informazioni che volevamo ottenere.
Le camere a ionizzazione non hanno una buona velocità di risposta, quando il flusso del gas risulta
piuttosto elevato, per cui si incontrano difficoltà sostanziali nell’utilizzare questi strumenti come
rivelatori di impulsi singoli.
Nel periodo Gennaio, Febbraio 2001, risultò subito evidente che bisognava cambiare metodo di
ricerca, la strada che si stava perseguendo, avrebbe dato informazioni esclusivamente sulla
concentrazione relativa del Radon nell’ambiente.
La stessa risposta si sarebbe potuta ottenere acquistando un normale Radometro, ma l’incidenza del
costo, per un discreto strumento, e le risposte scientifiche che avremmo ottenuto, sarebbero state le
stesse acquisite fino ad oggi dalla scienza ufficiale; perciò si decise di seguire una strategia di ricerca
diversa da quella tradizionale.

Il grafico in basso, mostra la catena di decadimento dell’ 238Uranio ed è da questa catena che si è

intuito, su quale ramo di decadimento del Radon si dovesse focalizzare l’attenzione per realizzare un
nuovo rivelatore che permettesse misure correlabili con eventi sismici.
Attraverso un’attenta osservazione della catena di decadimento del 222Rn, in modo particolare nella
fase in cui vengono prodotti gli isotopi 214Pb, 214Bi e 214Po, si osserva che questi elementi, a differenza
degli altri, anziché decadere con processo alfa, decadono con processo beta, con emissione di fotoni
gamma.
(L’immagine seguente mostra la sezione della catena, nella fase di decadimento beta)

Disponendo quindi di alcuni Fotomoltiplicatori e di scintillatore plastico, è stato realizzato un primo
prototipo che ci avrebbe permesso di contare flash di particelle gamma prodotte dal decadimento beta
del 214Pb e del 214Bi.

La rivelazione di questi elementi, con questo metodo, ci avrebbe permesso di rilevare impulsi singoli
con buona precisione; inoltre essendo essi diretti isotopi del Radon, il valore della loro variazione di
concentrazione, equivaleva alla variazione di concentrazione del Radon emesso dalla crosta terrestre.
Nasce così il primo prototipo di particelle gamma, costituito da 4 fotomoltiplicatori ed uno
scintillatore plastico.

Esso viene posto in acquisizione in uno scantinato ad una profondità di tre metri dalla superficie
stradale, nella zona di Coppito (AQ), nel mese di Maggio 2002.
 

L’immagine a sinistra mostra il primo rivelatore chiamato PM4, durante la fase di preparazione e prima di entrare in
funzione, presso la stazione di rilevamento di Coppito (AQ)

Il grafico a destra rappresenta il primo counting rate ottenuto dopo i primi periodi di acquisizione del rivelatore
PM4 di Coppito (AQ). Sulle ascisse sono riportate le ore relative ai giorni di presa dati, dal 27/06/02 al 01/07/02
Sulle ordinate è riportato il numero degli eventi ogni 7200 secondi, (Counting rate).

Già dalle prime settimane di presa dati, i risultati apparvero più che soddisfacenti.
Si era centrato l’obbiettivo! Per la prima volta, si osservava un segnale continuo che rappresentava la
variazione di concentrazione di Radon, ottenuta dall’osservazione di suoi diretti figli di decadimento:
gli isotopi di 214Pb e 214Bi.

Dopo alcuni mesi, dedicati alla taratura dell’elettronica ed allo studio del segnale prodotto, vengono
ottenute le prime correlazioni tra variazione di concentrazione di Radon ed eventi sismici.
Nel periodo Luglio-Ottobre 2002, i dati osservati vengono analizzati con il metodo della Varianza.
Ogni settimana, veniva controllata la media di flusso del Radon, cercando anomalie nel counting rate
acquisito; contemporaneamente il sismografo monitorava l’attività sismica relativa a piccoli eventi
registrati in prossimità delle coordinate del rivelatore PM4, (Lat. +42° 21’N, Long. +13° 20’E).
L’immagine a sinistra mostra il primo rivelatore chiamato PM4, durante la fase di preparazione
e prima di entrare in funzione, presso la stazione di rilevamento di Coppito (AQ)
 

Il grafico a destra rappresenta il primo counting rate ottenuto dopo i primi periodi di
acquisizione del rivelatore PM4 di Coppito (AQ).
Sulle ascisse sono riportate le ore relative ai giorni di presa dati, dal 27/06/02 al 01/07/02
Sulle ordinate è riportato il numero degli eventi ogni 7200 secondi, (Counting rate).

Giorno Lat. Long. Area Magnitudo
03/08/2002 42°.6 13°.0 C. Italia M 3.2
06/09/2002 41°.9 12°.5 S. Italia M 4.5
07/09/2002 41°.6 15°.8 S. Italia M 2.7
09/09/2002 42°.4 11°.8 C. Italia M 4.5
24/09/2002 41°.7 13°.3 S. Italia M 2.5
06/10/2002 39°.6 12°.4 Tirreno M 3.2
19/10/2002 40°.7 12°.7 Tirreno M 2.9
23/10/2002 42°.7 17°.3 Adriatico M 4.0

Ii- counting rate di due ore
La tabella di seguito mostra i primi allarmi prodotti dall’analisi della media del flusso di Radon
osservato e correlato con gli eventi sismici, registrati e dal nostro sismografo e dalla Sala sismica
centrale dell’INGV di Roma.(www.ingv.it)
PRIME CORRELAZIONI TRA VARIAZIONI DI
CONCENTRAZIONE DI RADON ED EVENTI SISMICI
DATA LAT. LOMG. ORA Previsto MAG.

Ev.
03/08/2002 42°.6 13°.0 19.21 02/08/2002 M. 2.7
06/09/2002 41°.9 12°.5 13.45.00 05/09/2002 M. 4.5
07/09/2002 41°.6 15°.8 7.38.16 06/09/2002 M. 2.7
09/09/2002 42°.4 11°.8 2.53.46 08/09/2002 M. 4.5
24/09/2002 41°.7 13°.3 15.22.09 23/09/2002 M. 2.5
06/10/2002 39°.6 12°.4 1.22.41 05/10/2002 M. 3.2
19/10/2002 40°.7 12°.7 16.02.51 18/10/2002 M. 2.9
23/10/2002 42°.7 17°.3 11.01.25 22/10/2002 M. 4.0
29/10/2002 serie eventi S. Giuliano 31/10/02
30/10/2002 serie eventi S. Giuliano 01/11/02
43°.6 12°.8 11.32 25/03/2003
26/03/2003 19.01 25/03/2003 M. 2.5
30/03/2003 41°.6 14°.7 16.42 29/03/2003 M. 3.3
11/04/2004 42°.30 13°.40 11.04 10/04/2003 M. 2.4

In questo stesso periodo era stata presentata alla Regione Abruzzo e Protezione Civile Abruzzo, una
richiesta di finanziamento per la nostra ricerca, nella speranza di poter ammortizzare le spese del
nostro progetto, che seppur condotto in maniera artigianale, incideva pesantemente sul bilancio
famigliare.

In cambio del finanziamento promesso, la Protezione Civile dell’Aquila, avrebbe ricevuto dalla sala
sismica della S.C.S., gli allarmi > del 3° Richter, previsti il giorno prima dell’evento.
Per motivi che qui non stiamo a sottolineare, la S.C.S. non ha mai ricevuto nessun genere di aiuto.
La tabella in alto mostra gran parte degli eventi allarmati nel periodo Agosto – Ottobre 2002 alla
Protezione Civile di L’Aquila.

In particolare, le date segnate in rosso nella tabella, mostrano gli allarmi trasmessi, all’allora
Assessore Regionale ai Lavori Pubblici, con delega per la Protezione Civile, Dr. Giorgio De Matteis.
Il 20 Dicembre 2002 viene depositata domanda di brevetto sul sistema di “Precursore di eventi
sismici”

Anno 2003
Sulle ali dell’euforia, il gruppo di ricerca della S.C.S., Roberto Giuliani, Victor Alekseenko e
Giampaolo Giuliani, per i risultati ottenuti sulle correlazioni Radon - terremoti e per l’individuazione
del segnale chiamato precursore sismico, decise di incrementare lo sforzo della ricerca in corso.
Il gruppo era convinto che la scoperta del sistema per rilevare allarmi sismici da 6 a 24 ore di anticipo
sull’evento, sarebbe stata apprezzata sia nell’ambiente scientifico italiano, che negli ambienti della
Protezione Civile nazionale.

Confortati dalla certezza, fu deciso di costruire un secondo rivelatore da posizionare ad una certa
distanza da quello in funzione a Coppito (AQ).
Con l’elettronica e la meccanica ancora disponibile, si dette quindi inizio alla realizzazione di un
secondo detector, con caratteristiche identiche al primo.
Per questo secondo rivelatore furono utilizzati 2 Fotomoltiplicatori anziché 4 come il precedente.
Questa scelta avrebbe permesso di verificare la riproducibilità del sistema, l’efficienza, la sensibilità e
l’eventuale correlazione con siti tra loro distanti e diversi.

Bisognava nello stesso tempo incrementare le analisi sui dati forniti dal PM4.
Il dr. V. Alekseenko propose un algoritmo sulla Varianza delle 24 ore del 222Rn, per mezzo del quale
sarebbero state effettuate correlazioni con le variazioni di concentrazione di Radon ed un maggior
approfondimento sui diversi fenomeni, quale ad esempio le armoniche, osservate nel continuo dei dati
rilevati.

Il Sig, G. Giuliani propose un algoritmo sulla media di flusso del 222Rn in ambiente chiuso e senza
ventilazione, rilevato ogni 2 ore.
Questo algoritmo sarebbe stato applicato al programma automatico di allarmi sismici sulla rete
internet, con dei messaggi via e-mail.

Formula dell’algoritmo applicato:
| CI – CI+n | > 3 < C > n = 1,…..,12
Dove: < C > = å12
1
Ci /12

Il Sig. R. Giuliani avrebbe organizzato invece un software per il controllo dei dati acquisiti, adattabile
ad un diverso numero di sistemi tra loro interallacciati e posti sulla rete internet, per ottenere in tempo
reale preallarmi ed allarmi attraverso le soglie dei precursori sismici connessi in modo nodale.
Esempio di allarme sismico del 29 dicembre 2005, via e-mail e relativo al sisma del 30 dicembre
2005 Alta_Val_Tiberina delle ore 10:37U.T. pari a 2,5M:
pm4
ALLARME
Massimo: 29/12/2005 02:00:00 150117
Minimo: 29/12/2005 16:00:00 146780
Differenza: 3337
Arco Tangente: -1,5686986373335159

La necessità di un secondo rivelatore, posto ad una certa distanza dal primo, avrebbe definitivamente
risposto ai diversi quesiti che assillavano un po’ tutti.
Fino ad allora non si aveva notizia nell’ambiente scientifico, di osservazioni effettuate sull’emissione
di Radon dalla crosta terrestre, controllata on-line, in tempo reale e da diversi siti di osservazione.
Gli interrogativi che chiedevano risposte erano:
- Quale andamento dei flussi sarebbe stato osservato dai rivelatori in siti diversi?
- Quale andamento avrebbero avuto i flussi monitorati, in prossimità di eventi sismici, rilevati
dai diversi detectors?
- Il flusso e l’andamento del Radon misurato dipendeva anche dalla natura del terreno dove
venivano situati i rivelatori?
Naturalmente l’interrogativo più importante era costituito dalla necessità di sapere se due o tre
rivelatori potessero correlarsi tra loro e fornire, con una certa precisione, epicentro, grado sismico ed
ora dell’evento.

Con un solo rivelatore purtroppo, non era possibile dare una stima sull’epicentro dell’eventuale sisma
allarmato, se non nel raggio d’azione del rivelatore; mentre era possibile avere una stima del grado
magnitudo con un errore di circa 0.5°-0.6° di Grado.
Verso Maggio-Giugno 2003, il secondo rivelatore con 2 Fotomoltiplicatori, era pronto per entrare in
funzione.
Per un periodo di qualche mese, sarebbe stato tarato, al fianco del fratello chiamato PM4.
Pur lavorando nel campo della ricerca dei precursori sismici, con molta discrezione ed avendo avuto
relazioni solo con dirigenti della Protezione Civile dell’epoca, fummo contattati nei mesi Giugno -
Luglio 2003 dall’imprenditore Vincenzo Passarelli.
L’imprenditore si disse interessato alla richiesta di brevetto da noi presentata in Italia, sul metodo e
tecnica per la rivelazione dei precursori sismici, garantendo la possibilità di commercializzare
all’estero i nostri rivelatori.
Si giunse ad un accordo definitivo, che prevedeva la realizzazione di due reti di monitoraggio anti
sismiche, una in Algeria ed una in Turchia, finanziate dai due governi.
Vincenzo Passatelli avrebbe finanziato un test ufficiale da effettuarsi a Reggio Calabria, sotto l’egida
di un gruppo di Geofisica dell’Università di quella città.

Purtroppo la trattativa si interruppe proprio il giorno della firma per l’accordo, presso l’ufficio del
legale preposto alla stesura del contratto.
Nel mese di Ottobre 2003, la Caen S.p.a. di Marcello Givoletti ci propose di subentrare al posto
dell’imprenditore Passatelli.
La Caen rappresenta a tutt’oggi la più importante società italiana per la produzione di apparecchiature
elettroniche per la fisica nucleare.
La possibilità di essere sostenuti da una così importante società, ci spinse ad accettare l’offerta del
Presidente Marcello Givoletti.

Data la delicatezza dell’argomento e pur con tutta la “prudenza” necessaria, d’accordo con Marcello
Givoletti (Caen), sentimmo il dovere di parlare della scoperta con qualche Istituzione Scientifica.
Il problema era: “cosa fare dei dati che erano stati rilevato” e soprattutto “cosa sarebbe stato giusto
ancora fare” per tentare di dare, se possibile, una risposta scientifica e rigorosa sull’argomento.
Più in particolare si pensava di :
- Aprire una o più collaborazioni con Istituzioni o Università?
- Collaborare con qualche Istituto per migliorare l’efficienza del rivelatore di particelle finora usato?
- Creare una rete di almeno 10-15 sistemi di misura distanti fra loro circa 50 Km col fine di
interpolare le misure di ogni rivelatore (per ottenere una maggiore accuratezza sulla stima
dell’epicentro) come dimostrano i dati fino ad ora analizzati?
- Studiare un sistema di calibrazione per ogni singolo sistema in funzione della diversa struttura
geologica del terreno di misura.
- Integrare le misure con altre informazioni (campo elettrico e magnetico terrestre ??)
- Reperite risorse umane e finanziarie per portare avanti un simile programma (spese vive valutate
intorno a 500 K€).

Nella seconda metà del 2003 entrammo in contatto con il Sen. Zamberletti, presidente dell’ISPRO
(www.ispro.it) e notoriamente da sempre coinvolto col fenomeno Terremoti.
Il Sen. Zamberletti si dichiarò molto interessato alla scoperta ed organizzò a Roma presso i locali
dell’ISPRO, una riunione allargata a una decina di persone interessate all’argomento fra le quali
citiamo il Prof. Enzo Boschi (Presidente INGV), il Dott.Galanti, Direttore Servizio Sismico
nazionale, il Gen. Mollicone (Aeronautica Militare), ecc….
Gioacchino Giuliani espose nella circostanza quanto era stato fatto fino a quel momento,
sottolineando, che tutte le spese erano state sostenute da un auto finanziamento famigliare.
Espose inoltre, quanto scientificamente fosse stato svolto fino ad allora , senza alcuna pretesa, propose
inoltre se poteva essere presa in considerazione la formazione di un gruppo di lavoro congiunto per
verificare e ampliare il programma di ricerca.

Le reazioni dei singoli partecipanti sono state molto diverse fra loro (in particolare quella del Prof.
Boschi e’ stata a dir poco “scortese”) , ma complessivamente la riunione non ha portato a nessuna
conclusione positiva. Va segnalato che nessuno allora (e fino a tutt’oggi) ha trovato errori o inesattezze sulle misure
effettuate! Al fine di dimostrare ufficialmente che i nostri dati erano attendibili, fu accettata la proposta dell’On.
Zamberletti, di sottoporre il sistema ad un test di funzionamento. Il Sistema avrebbe inviato dal 20 Gennaio 2003
all’8 Gennaio 2004, gli “allarmi e preallarmi segnalati” come sismi.

Gli allarmi sarebbero dovuti pervenire un “notaio” di Viareggio, (Al fine di certificarne la data e l’ora
certa), alla Caen S.p.a. ed alla ISPRO dell’On. Zamberletti.
Gli allarmi e preallarmi spediti, hanno anticipavano gli eventi sismici, con un intervallo di tempo
variabile dalle 6 e le 24ore. L’operazione riuscì perfettamente.

Nell’incontro successivo effettuato alla ISPRO, il Sen. Zamberletti avanzò l’ipotesi che un
programma di ricerca congiunto, sarebbe stato finanziato dalla società che coordina la costruzione del
“Ponte sullo stretto di Messina”
Nel mese di Dicembre 2003, venne sottoscritto in Viareggio, alla presenza di un notaio, l’accordo con
la Caen S.p.a. di Marcello Givoletti, per la costituzione di una società in L’Aquila, la “Caen-
Geo”entro il 31 Dicembre 2004.

La Caen S.p.a. avrebbe sostenuto tutte le spese di industrializzazione per la realizzazione dei rivelatori
PM4 e PM2 ed alla copertura della spesa per l’estensione del brevetto internazionale, a nome del Sig.
Gioacchino Giuliani; in cambio avrebbe ottenuto l’esclusiva per 20 anni, della commercializzazione
del brevetto stesso.

Anno 2004
Nel mese di Gennaio 2004, i dati del secondo rivelatore (PM2) mostrava risultati eccezionali, il suo
counting rate, era perfettamente sovrapponibile al counting rate del PM4.
Questo significava che il primo prototipo realizzato era perfettamente riproducibile, come mostra
l’immagine seguente:

L’istogramma nella figura mostra il counting rate del PM4 e del PM2, ottenuto ponendo i due
rivelatori a 2 metri di distanza tra loro e nello stesso ambiente.
Per il test successivo sarebbe stato indispensabile trovare un luogo dove mettere in acquisizione il
PM2, possibilmente ad una distanza dal PM4, > di 20 Km, < di 80 Km.
Verso l metà del 2004 il Prof. Alessandro Bettini, all’epoca direttore dei Laboratori Gran Sasso, visto
il clamore creato da alcuni giornalisti locali, che a gran voce gridavano che un ricercatore aquilano
aveva scoperto un sistema per prevedere i terremoti, ci propose di far rianalizzare tutti i dati fino
allora prodotti.

Fu accetta la proposta e furono consegnati i dati per la loro rielaborazione.
Dopo tre mesi fu consegnata una relazione al prof. A. Bettini, in cui veniva specificato che la
rielaborazione dei dati, conduceva alle stesse deduzioni del Sig. Giuliani. Quel metodo permetteva di
prevedere eventi sismici con 6 – 24 ore di anticipo.
Da un colloquio con la Senatrice Maria Claudia Ioannucci, riuscimmo ad ottenere dall’Onorevole
Borghini, responsabile del nucleo industriale “SI Sviluppo Italia”, l’opportunità di installare
nell’incubatore di Avezzano (AQ) il PM2 testato.

Nello stesso periodo la Caen mise a disposizione della S.C.S. 2 ricercatori, il Dr Nicola Zaccheo ed il
Dr Claudio Raffo, per rianalizzare tutti i dati fino a quel momento prodotti.
L’acquisizione, i dati, gli algoritmi di calcolo, le correlazioni, mostrano ripetitività ed attendibilità,
suffragate dai test effettuati dal dr N. Zaccheo, attraverso le misure effettuate per ottenere lo spettro di
energia nella regione di misura dei prototipi realizzati.

Vengono così ottenere “correlazioni” attendibili ed abbastanza ripetitive fra le misure effettuate
(preallarmi ed allarmi) e gli eventi sismici effettivamente avvenuti e riscontrati, oltre che dallo staff
della S.C.S., anche dai ricercatori della Caen.

Da marzo 2004 ad ottobre 2004, con la messa in funzione del secondo detector nella città di
Avezzano, in modo nodale ed in rete con quello dell’Aquila, attraverso l’analisi dei dati ottenuti, è
stato possibile non solo, perfezionare tutto il sistema di acquisizione ma anche dare una risposta
definitiva sul comportamento del radon misurato da punti di monitoraggio diversi ed in tempo reale.
Nella prima metà del 2004, è stata realizzata in l’Aquila, l’analisi sui dati acquisiti dal detector
chiamato “PM-4”, dal luglio 2002 al settembre 2003, estrapolando dai 4256 terremoti avvenuti in
Italia nello stesso periodo, i 381 terremoti registrati in un raggio di circa 1° geografico nel raggio
d’azione del “PM-4” applicando ad essi il metodo della Varianza ed ottenendo un’efficienza di
correlazione del 86% su tutti gli eventi correlati da 2.0° M a 5,6° M.

Tra dicembre 2004 e gennaio 2005, è stata portata a termine l’analisi e le correlazioni dei dati
acquisiti dai due detector posti tra L’Aquila ed Avezzano, applicando ad essi l’algoritmo sulla
variazione di flusso del 222Radon a 12 ore ottenendo una efficienza di circa il 92% - 94%, rispetto ad
un solo detector e rispetto alla Varianza applicata sulle 24 ore.
La stessa analisi dei dati acquisiti tra i due detector Pm-4 e PM-2 ha permesso di elaborare
l’algoritmo, che triangolando con più detector tra loro interconnessi, permette di ottenere l’epicentro
dell’evento allertato, all’interno di un raggio di 2000m – 3000m.

Gli istogrammi sopra e quello di sotto, mostrano le prime correlazioni tra le stazioni di monitoraggio
di Coppito (AQ) ed Avezzano (AQ), tra loro distanti circa 40 Km.
Le risposte fino ad allora teorizzate si concretizzano subito, con le analisi dei primi istogrammi, che
mostrano l’andamento correlato del flusso di Radon, i precursori sismici correlati dalle diverse
stazioni, la possibilità di correlare tra loro epicentro e grado sismico degli eventi.
 

Alla luce di questi risultati la Caen ed il Sig. Giuliani, si adoperano allora per trovare altri gruppi di
ricerca interessati a lavorare sull’argomento e vengono firmate due convenzioni, con le seguenti
 

Università :
Università della Calabria (Capo gruppo Prof. Ignazio Guerra)
Università di Bari – Dipart. Di fisica (Capo gruppo Prof Franco Romano)
In ambedue i casi le Università sono prive di fondi da dedicare a questa ricerca.
Venne quindi presentata una proposta operativa alla Società “Stretto di Messina SpA”, con l’analisi
dei costi ripartiti fra i partecipanti al progetto (CAEN e l’Università della Calabria, l’Università della
Basilicata, l’Università di Messina e l’Università di Palermo)
Nessuna notizia dalla Soc. del Ponte sullo Stretto di Messina !!


Dopo mesi di stop al programma e di estenuante attesa, Caen ed il Sig. G. Giuliani discutono della
cosa con alcuni ricercatori USA ricevuti in Viareggio, insieme al Prof. Edward Nixon (Geologo e
fratello dell’ex Presidente USA Richard Nixon).
Il Prof E. Nixon e’ interessato al programma ed incontra gli altri “Partner” delle due Universita’
Italiane e si decide di avviare le procedure per la richiesta di fondi USA al DOE (Dipartimento
dell’Energia) e DOI (Dipartimento degli Interni).
Anno 2005


Nessuna notizia dalla Soc. del Ponte sullo Stretto di Messina !!
Il Prof. E.Nixon e’ in continuo contatto con CAEN, ha gia coinvolto altri ricercatori e Università USA
e hanno gia individuato in California un sito che sembra “ideale” per la ricerca e spera di avere i fondi
necessari all’inizio del 2006.
Nel mese di Aprile la regione Toscana ha manifestato interesse a monitorare alcune aree a rischio
sismico utilizzando il Sistema CAEN.
I contatti sono tuttora in corso.
Nell’attesa che si concretizzassero i contatti con gli enti sopra citati, il Sig. G. Giuliani propose al
Presidente della Caen, M. Givoletti, di collaborare, alla realizzazione di un prototipo sottomarino, per
verificare se gli stessi fenomeni sul flusso di Radon, osservati sulla superficie terrestre, fosse possibile
osservarli anche sul fondo del mare.

La Caen rifiutò la possibilità di collaborare alla ricerca dei precursori sismici sottomarini, adducendo
che la scoperta fin qui realizzata era più che sufficiente e che sarebbe stato meglio spendere tutte le
risorse per cercare di convincere enti di ricerca e governativi ad interessarsi di ciò che fino ad allora
era stato realizzato.
Dal 20 Aprile 2004 i Sig. R. Giuliani e G. Giuliani, ritenendo indispensabile per le ricerche in corso,
dare una risposta alla possibilità di misurare di Radon dal fondo del mare, decidono di realizzare un
prototipo rivelatore di Radon sottomarino.
G. Giuliani acquista una vecchia barca da pesca destinata alla demolizione e da tutti considerata non
più adatta alla navigazione.

Sotto la guida di un maestro d’ascia, dopo 6 mesi la vecchia barca viene trasformata in un natante
laboratorio, per effettuare ricerche sul fondo marino sui medesimi precursori sismici osservati sulla
terra ferma.

Come si presentava la barca prima del recupero, abbandonata in secca da 6 anni, in un piccolo
cantiere di Roseto Degli Abruzzi. 20 Aprile 2005. Stazza 4 Ton. L.f.t. 8m Semi cab.
Motore BMV 3 Cil. 45CV.

Dopo un anno di preparativi e test effettuati al largo della costa Adriatica, il 29 Aprile 2005, sono stati
registrati i primi dati acquisiti a 15 metri di profondità, a circa 5 miglia dalla costa.

Di fianco la barca laboratorio dopo il restauro.
L’immagine mostra la barca prima della partenza, per gli ultimi test effettuati al largo della costa
adriatica, tra Roseto e Giulianova.

Il prototipo rivelatore sottomarino utilizza la stessa tecnica dei rivelatori di terra, l’unica differenza
riguarda l’ingegnerizzazione del detector, che prevede una tenuta stagna, per operare su fondali marini
da 15m a 40m di profondità.
Il grafico sotto mostra le prime misure di Radon effettuate a 18 m di profondità a circa 5 mgl dalla
costa tra Roseto (Te) e Giulianova (Te).

Dopo questi risultati si decise quindi di istallare una stazione di rilevamento nella cittadina di Pineto
(Te) per correlare i dati ottenuti dalle stazioni di Coppito (Aq) e quella dei Laboratori Nazionali Gran
Sasso in Assergi (Aq).
La stazione di Pineto, avrebbe effettuato correlazioni con i dati rilevati dalla barca laboratorio al largo
della costa. Di fianco la barca laboratorio dopo il restauro.
L’immagine mostra la barca prima della partenza, per gli ultimi test effettuati al largo della costa
adriatica, tra Roseto e Giulianova.
 

Questa ricerca sulle correlazioni tra il fondo del mare e la costa è a tutt’oggi in corso, anche se, a
causa di mancanza di fondi, al momento subisce un forte rallentamento.
G. Giuliani.
 

Dal 16 gennaio scorso, chi abita nell'aquilano ha un pò la tremarella:. Questo è il prospetto stralciato da gentili indicazioni dell'Istituto Nazionale
di Geofisica e Vulcanologia circa le scosse registrate nel mese di Febbraio 2009.

Data Ora Prof(Km) Mag Distretto Sismico
2009/02/17 18:13:06 9.6 Ml:2.4 Aquilano
2009/02/17 13:13:03 10.4 Ml:2.3 Aquilano
2009/02/17 06:08:34 9 Ml:2.6 Aquilano
2009/02/16 23:16:39 10.0 Mw:5.3 Ionian Sea
2009/02/16 18:43:08 124.8 Ml:2.7 Isole_Lipari
2009/02/16 12:41:53 10.4 Ml:1.8 Aquilano
2009/02/15 19:16:38 10.2 Ml:2.5 Aquilano
2009/02/15 07:41:40 11.5 Ml:2 Aquilano
2009/02/13 10:29:57 10 Ml:2.5 Aquilano
2009/02/13 06:39:13 8.6 Ml:3 Alta_Val_Tiberina
2009/02/12 12:12:58 68 Ml:2.7 Costa_siciliana_settetrionale
2009/02/11 17:34:53 33 Mw:7.5 Talaud Islands, Indonesia
2009/02/10 07:57:25 11.1 Ml:2.4 Prealpi_lombarde
2009/02/09 14:18:16 8.8 Ml:2.6 France
2009/02/07 08:02:35 6 Ml:3.4 Madonie
2009/02/05 22:29:48 5 Ml:2.4 Piana_del_Volturno
2009/02/05 17:43:55 10 Ml:2.8 Tirreno_centrale
2009/02/05 14:50:13 30.7 Ml:3.1 Mar_Ionio
2009/02/05 07:43:12 10 Ml:1.8 Irpinia
2009/02/05 03:55:22 22.8 Ml:3.2 Mar_Ionio
2009/02/04 16:57:32 10 Ml:2 Valle_del_Topino
2009/02/03 20:17:18 6.9 Ml:1.4 Garfagnana
2009/02/01 14:52:00 5.8 Ml:3.8 Mar_Ligure
2009/02/01 05:52:04 9.1 Ml:2.2 Aquilano
2009/02/01 00:57:09 35 Mb:4.8 Greece
+ altre sette trettecate nel mese di gennaio, degne di nota.

Sul sito www.ilcapoluogo.it, si dice però che l'Aquila dispone, "unica al mondo", di un sofisticato sistema di "previsione dell'attività sismica",
brevettato nel 2003! da un geofisico aquilano.

"Il "PM4", è una scatola pesante mezzo quintale, semplice come tutte le opere di genialità, realizzata e brevettata da Giampaolo Giuliani, un tecnico aquilano in servizio presso i Laboratori Nazionali del Gran Sasso, e prevede i terremoti in un ampio raggio d'azione e con un anticipo compreso tra le 6 e le 24 ore, con una precisione molto alta. I precursori sismici vengono incrociati con altre 4 stazioni così attrezzate che si trovano nella frazione di Coppito, nei laboratori del Gran Sasso, nel vicino Comune di Fagnano e a Pineto (Teramo). La triangolazione dei dati consente di individuare con una discreta precisione epicentro e intensità di un evento sismico. Il sistema è stato collaudato dal Dipartimento di Ingegneria delle strutture, delle acque e del terreno (Disat) dell'Università dell'Aquila. I rivelatori - a differenza dei sismografi, che non consentono di osservare un precursore sismico - percepiscono l'energia di una faglia e si basano sull'analisi del decadimento del radon, gas nobile con vita media 3,8 giorni: i 2 isotopi Piombo 214 e Bismuto 214 emettono fotoni percepiti da uno scintillatore plastico e visualizzabili grazie a fotomoltiplicatori, collocati in un cubo di piombo del peso di 6 quintali, con pareti spesse 7 centimetri. Molte cose sono cambiate nel breve volgere di qualche mese. I risultati raggiunti da Giuliani, dopo anni di studi, test e lunghi esperimenti svolti a L'Aquila, erano regolarmente snobbati dalla scienza ufficiale e dai geofisici. Secondo questi ultimi, prevedere i terremoti era impossibile, ma la caparbia (e la precisione dei dati) del tecnico aquilano sta gradualmente intaccando il muro di gomma eretto dagli scettici. L'ideatore del "PM4" è entusiasta per i risultati registrati: "Abbiamo previsto tutti gli eventi sismici con un'esattezza indiscutibile. La protezione civile riceve regolarmente i nostri preallarmi e allarmi, emessi dai nostri rilevatori. Qui, per adesso, vengono solamente analizzati dal tecnico Roberto Ianni il quale ne controlla l'intensità. La scossa di ieri mattina (alle 7:08, magnitudo 2.6), ad esempio, era stata preannunciata con un allarme della stazione di Coppito alle 18 di ieri. Nell'ultimo mese le scosse avvertite nel territorio aquilano - una o due al giorno - sono state puntualmente previste dalle stazioni gestite da Giuliani il quale tranquillizza: "Il fatto che ci siano piccoli terremoti ogni giorno non consente che si accumuli energia tale da scatenare un evento più forte". Ma lo sciame sismico dura ininterrottamente da un mese, cosa sta accadendo? "Il primo evento, ampiamente previsto dai nostri rilevatori, c'è stato il 16 gennaio" spiega Giuliani."Da allora, escludendo un paio di giorni di pausa sismica, ha scaricato praticamente tutti i giorni. Gli epicentri e ipocentri più importanti sono sempre gli stessi e cioè Roio (località aquilana) e la bassa valle dell'Aterno. Si tratta di uno sciame tellurico del tutto normale - prosegue Giuliani - in quanto eventi simili si sono verificati anche negli anni passati, ma quest'anno alcuni terremoti sono stati percepiti dalla popolazione". Nel corso dei suoi studi sui terremoti, Giuliani e il suo staff hanno monitorato 4.500 scosse sismiche registrate in tutta Italia, dalla Valle d'Aosta a Pantelleria, giungendo a un'altra clamorosa scoperta: "Abbiamo notato che la maggior parte dei terremoti si verificano durante i mesi invernali, o meglio, quando il sistema terra-luna è nel perielio, quindi più vicino al sole. In inverno, quando la terra subisce uno stress gravitazionale maggiore, si registrano più eventi sismici (60-70%) che in estate. La percentuale si mantiene ancora più alta quando c'è la luna nuova. Il magma che scorre sotto la crosta terrestre risente delle attrazioni gravitazionali, come accade per gli oceani". Giuliani ha fiducia sul fatto che la scoperta del secolo "made in L'Aquila" possa dare lustro alla città: "Spero con tutto il cuore che i benefici che si potranno trarre da questo brevetto restino a L'Aquila. Mi basta che il mondo sappia di questa opportunità e che un domani, magari molto prossimo, si possano divulgare le "previsioni sismiche" come accade per le previsioni metereologiche. Tutto questo ci offre il lato per congedarci da Giuliani con la domanda che tutti vorrebbero fargli: la popolazione civile sarà avvisata se il "PM4" dovesse prevedere una scossa di forte intensità? "La gente - risponde - potrebbe essere avvertita con 24 ore di anticipo e ci sarebbe tutto il tempo necessario per premunirsi. Ma per arrivare a tale fine bisogna innanzitutto educare la popolazione civile, e non solo". (di Luca Ceccarelli) "
 

Angela Baglioni L'AQUILA

Il susseguirsi di scosse, anche se di bassa magnitudo, sta ormai diventando una costante per il territorio aquilano. Nella giornata di ieri i sismografi dell'Istituto nazionale di geofisica e vulcanologia di Roma hanno registrato altri due terremoti. Il primo, di magnitudo 2, si è verificato alle 11,46. Il secondo, pari a 2,4 gradi della scala Richter, è stato invece avvertito alle 13,11. Nella città dell'Aquila, ormai, serpeggia una vera psicosi e il ripetersi frequente di movimenti tellurici non aiuta a tranquillizzare gli animi, nonostante le rassicurazioni da parte della comunità scientifica. Tra l'altro, ai terremoti veri si sono purtroppo aggiunte nei giorni scorsi notizie allarmanti diffuse da un camioncino su presunte scosse distruttive, creando notevole allarme tra la popolazione.

La Digos della Questura dell'Aquila, a questo proposito, ha avviato delle indagini per tentare di risalire all'identità di coloro che, subito dopo i terremoti più importanti, si sono "prodigati" nell'invitare la gente ad abbandonare le rispettive abitazioni preannunciando scenari catastrofici. Scenari, è bene ripeterlo, smentiti in maniera categorica dalla protezione civile che a più riprese ha ribadito come alla luce delle attuali conoscenze scientifiche non sia assolutamente possibile prevedere i terremoti. Fino a questo momento sono state ascoltate due persone ma l'attività investigativa prosegue. Eppure, ancora durante la giornata di ieri, circolavano voci sull'allestimento di una tendopoli a Tornimparte, dove far "convergere" la popolazione in caso di sisma. Anche questa notizia, ovviamente, non risponde al vero ma testimonia il clima di quasi panico che ormai si è diffuso. Gli investigatori non escludono il sospetto che qualcuno abbia interesse ad evacuare le case per motivi non certo legati alla pubblica incolumità, ma più prosaicamente per poter approfittare dell'assenza dei proprietari. Le forze dell'ordine, comunque, non tralasciano alcuna pista, neanche quella della goliardata di pessimo gusto. Nei giorni scorsi all'Aquila si è riunita la commissione grandi rischi presieduta da Franco Barberi. Alla riunione, voluta dal capo del dipartimento della protezione civile, Guido Bertolaso, hanno partecipato tutte le componenti istituzionali coinvolte ai vari livelli. Tra gli altri, erano presenti l'assessore regionale Daniela Stati e il sindaco dell'Aquila, Massimo Cialente. Dopo l'incontro è stato ribadito ancora il messaggio che non è possibile prevedere i terremoti.
Verosimilmente, hanno detto, è possibile ipotizzare lo scenario a partire dalle caratteristiche geo-morfologiche del territorio. Di sicuro, non si potrà sapere quando un evento sismico si verificherà. Lo sciame che interessa il territorio aquilano e la conca di Sulmona va avanti ormai da più di tre mesi, con decine di movimenti tellurici registrati dai sismografi. Nella sola giornata di lunedì scorso, quando la sequenza ha raggiunto l'acme con la scossa di magnitudo 4, la popolazione ha avvertito nell'arco di appena nove ore ben dieci terremoti. La protezione civile mantiene alto il livello di attenzione anche se, come ribadito dagli esperti, si tratta di una sequenza tutto sommato normale per un territorio ad alto grado di sismicità, come quello aquilano. Il danno atteso, ha detto il vice capo del dipartimento operativo nazionale di protezione civile, Bernardo De Bernardinis a margine dell'incontro di martedì, non dovrebbe essere molto diverso da quello che è stato determinato finora e che non ha tuttavia interessato strutture portanti di edifici, salvo qualche rara eccezione dovuta più che al terremoto a carenze gravi già note prima dell'avvio dello sciame sismico. Il consigliere regionale del Pdci, Antonio Saia, nel frattempo, ha presentato un'interrogazione per sapere quali provvedimenti intende assumere la Giunta in relazione ai danni causati dagli eventi sismici e quali iniziative saranno intraprese per ultimare la ricostruzione degli edifici danneggiati dal terremoto di 25 anni fa.

A L'Aquila la prima sala di previsione sismica al mondo

E' esecutiva dall'agosto del 2008 e funziona perfettamente la prima stazione di previsione dei terremoti esistente al mondo.

Lo sciame sismico che dal 16 gennaio scorso sta interessando la città dell'Aquila è la circostanza ideale per testarne l'attendibilità. Lo abbiamo fatto andando a trovare chi ha realizzato e brevettato il "PM4", una magica scatola di piombo
che permette di prevedere un evento sismico dalle 6 alle 24 ore prima dell'evento.

La "sala sismica" per la previsione di terremoti è stata attrezzata all'interno della scuola elementare "De Amicis", a San Bernardino,
grazie alla concessione di alcuni locali da parte del Comune dell'Aquila che, assieme al dirigente provinciale Francesco Bonanni e al direttore dei Laboratori
del Gran Sasso d'Italia, Eugenio Coccia, hanno prestato fede e scommesso sull'invenzione di Giuliani.

Il "PM4", è una scatola pesante mezzo quintale, semplice come tutte le opere di genialità, realizzata e brevettata da Giampaolo Giuliani,
un tecnico aquilano in servizio presso i Laboratori Nazionali del Gran Sasso, e prevede i terremoti in un ampio raggio d'azione e con un anticipo
compreso tra le 6 e le 24 ore, con una precisione molto alta.

I precursori sismici vengono incrociati con altre 4 stazioni così attrezzate che si trovano nella frazione di Coppito, nei laboratori del Gran Sasso,
nel vicino Comune di Fagnano e a Pineto (Teramo). La triangolazione dei dati consente di individuare con una discreta precisione epicentro e intensità di un evento sismico.

Il sistema è stato collaudato dal Dipartimento di Ingegneria delle strutture, delle acque e del terreno (Disat) dell'Università dell'Aquila.
I rivelatori - a differenza dei sismografi, che non consentono di osservare un precursore sismico - percepiscono l'energia di una faglia
e si basano sull'analisi del decadimento del radon, gas nobile con vita media 3,8 giorni: i 2 isotopi Piombo 214 e Bismuto 214 emettono
fotoni percepiti da uno scintillatore plastico e visualizzabili grazie a fotomoltiplicatori, collocati in un cubo di piombo del peso di 6 quintali,
con pareti spesse 7 centimetri. Il brevettatore del sistema (brevettato nel 2003) è l'aquilano Gioacchino Giampaolo Giuliani,
tecnico di ricerca dell'Inaf (Istituto Nazionale di Astrofisica), associato Infn (Istituto Nazionale di Fisica Nucleare),
che dal 2002 ha speso tutte le sue energie sulla realizzazione di un sistema di previsione dei terremoti dopo aver riscontrato
che l'andamento del flusso di radon, un gas naturale, è strettamente legato ad un evento sismico.

Molte cose sono cambiate nel breve volgere di qualche mese. I risultati raggiunti da Giuliani, dopo anni di studi, test e lunghi esperimenti svolti a L'Aquila,
erano regolarmente snobbati dalla scienza ufficiale e dai geofisici. Secondo questi ultimi, prevedere i terremoti era impossibile,
ma la caparbia (e la precisione dei dati) del tecnico aquilano sta gradualmente intaccando il muro di gomma eretto dagli scettici.

L'ideatore del "PM4" è entusiasta per i risultati registrati: "Abbiamo previsto tutti gli eventi sismici con un'esattezza indiscutibile.
La protezione civile riceve regolarmente i nostri preallarmi e allarmi, emessi dai nostri rilevatori. Qui, per adesso, vengono solamente
analizzati dal tecnico Roberto Ianni il quale ne controlla l'intensità. La scossa di ieri mattina (alle 7:08, magnitudo 2.6), ad esempio,
era stata preannunciata con un allarme della stazione di Coppito alle 18 di ieri. Nell'ultimo mese le scosse avvertite nel territorio aquilano
- una o due al giorno - sono state puntualmente previste dalle stazioni gestite da Giuliani il quale tranquillizza:
"Il fatto che ci siano piccoli terremoti ogni giorno non consente che si accumuli energia tale da scatenare un evento più forte".

Ma lo sciame sismico dura ininterrottamente da un mese, cosa sta accadendo?
"Il primo evento, ampiamente previsto dai nostri rilevatori, c'è stato il 16 gennaio? spiega Giuliani.
"Da allora, escludendo un paio di giorni di pausa sismica, ha scaricato praticamente tutti i giorni.
Gli epicentri e ipocentri più importanti sono sempre gli stessi e cioè Roio e la bassa valle dell'Aterno.
Si tratta di uno sciame tellurico del tutto normale - prosegue Giuliani - in quanto eventi simili
si sono verificati anche negli anni passati, ma quest'anno alcuni terremoti sono stati percepiti dalla popolazione".

Nel corso dei suoi studi sui terremoti, Giuliani e il suo staff hanno monitorato 4.500 scosse sismiche registrate in tutta Italia,
dalla Valle d'Aosta a Pantelleria, giungendo a un'altra clamorosa scoperta: "Abbiamo notato che la maggior parte dei terremoti
si verificano durante i mesi invernali, o meglio, quando il sistema terra-luna è nel perielio, quindi più vicino al sole.
In inverno, quando la terra subisce uno stress gravitazionale maggiore, si registrano più eventi sismici (60-70%) che in estate.
La percentuale si mantiene ancora più alta quando c'è la luna nuova. Il magma che scorre sotto la crosta terrestre risente
delle attrazioni gravitazionali, come accade per gli oceani".

Giuliani ha fiducia sul fatto che la scoperta del secolo "made in L'Aquila" possa dare lustro alla città:
"Spero con tutto il cuore che i benefici che si potranno trarre da questo brevetto restino a L'Aquila.
Mi basta che il mondo sappia di questa opportunità e che un domani, magari molto prossimo,
si possano divulgare le "previsioni sismiche" come accade per le previsioni metereologi che.

Tutto questo ci offre il lato per congedarci da Giuliani con la domanda che tutti vorrebbero fargli: la popolazione civile sarà avvisata se il "PM4"
dovesse prevedere una scossa di forte intensità? "La gente - risponde - potrebbe essere avvertita con 24 ore di anticipo e ci sarebbe tutto il tempo
necessario per premunirsi. Ma per arrivare a tale fine bisogna innanzitutto educare la popolazione civile, e non solo".

Altro articolo:
Giampaolo Giuliani utilizza il Radon, un gas nobile emesso dalla crosta terrestre. Lo stesso nome e’ stato dato al suo apparecchio che dovrebbe prevedere i sismi.
Ma la protezione civile smentisce che lui ci riesca. Lui e’ convinto di aver creato l’oggetto giusto per prevedere tante catastrofi. Per questo i suoi cubi ora sono sparsi
per l’Abruzzo: nel 2004 reagirono anche con largo anticipo al sisma che causo lo Tsunami in Indonesia.
 
Domenica Giampaolo Giuliani ha previsto una violenta scossa nel pomeriggio, dopo quella della mattina, a Sulmona. Subito il fisico ha chiamato la polizia e il sindaco,
che si trovava a Roma al congresso del Popolo della Libertà. Peccato che il terremoto non si sia mai verificato. E con la sua falsa previsione ha causato un caos allucinante.
Tanto che ora e’ stato denunciato per aver creato falso allarmismo.
 
Il sindaco commenta cosi’: “Non è possibile che si vada in giro a creare allarmismi, domenica ho passato la giornata più brutta della mia vita perché dopo aver parlato con quel signore, che mi aveva annunciato un sisma devastante, molto più forte di quello che c’era stato in mattinata, mi sono trovato in una situazione difficilissima. Perché non sapevo che cosa fare: far scattare il piano di evacuazione o far finta di niente, con il rischio di portarmi sulla coscienza una possibile tragedia?“.
 
Ma Giampaolo Giuliani va avanti con il suo lavoro: lui e’ convinto che i suoi cubi dicano il vero!

Previsione analitica e non pronostico calcistico o profezia come e' stato scritto da alcuni.


 Scosse sismiche Aquila: Gennaio, Febbraio, Marzo.

Viste le numerose richieste degli utenti riportiamo la tabella delle scosse sismiche rilevate a L’Aquila e provincia, negli ultimi 4 mesi.

Sono dati puramente informativi dei quali non possiamo garantire l’assoluta correttezza.

la prima colonna riporta la data e l’ora in tempo locale.
la seconda colonna il magnitudo
la terza colonna la profondità dell’evento in chilometri
la quarta colonna la provincia interessata

26/03/2009 03:10:38 2 11 L Aquila
24/03/2009 09:27:54 1.8 10 L Aquila
22/03/2009 08:02:55 2.3 10 L Aquila
18/03/2009 10:23:09 2.2 9 L Aquila
13/03/2009 16:29:55 2.4 10 L Aquila
13/03/2009 14:49:14 1.8 9 L Aquila
13/03/2009 14:09:52 1.3 5 L Aquila
11/03/2009 22:55:49 2.9 9 L Aquila
10/03/2009 12:25:27 2.2 10 L Aquila
22/02/2009 08:03:40 1.3 12 L Aquila
21/02/2009 18:04:10 2.2 10 L Aquila
17/02/2009 19:13:06 2.4 10 L Aquila
17/02/2009 14:13:03 2.3 10 L Aquila
17/02/2009 07:08:34 2.6 9 L Aquila
16/02/2009 13:41:53 1.8 10 L Aquila
15/02/2009 20:16:38 2.5 10 L Aquila
15/02/2009 08:41:40 2 12 L Aquila
13/02/2009 11:29:57 2.5 10 L Aquila
01/02/2009 06:52:04 2.2 9 L Aquila
31/01/2009 00:00:01 2.3 10 L Aquila
30/01/2009 22:54:19 2.5 10 L Aquila
30/01/2009 01:39:03 2.5 9 L Aquila
24/01/2009 10:06:07 2.4 10 L Aquila
23/01/2009 19:51:56 2 8 L Aquila
18/01/2009 09:13:43 2.1 9 L Aquila
16/01/2009 23:47:23 2.4 11 L Aquila
14/12/2008 09:00:00 1.8 10 L Aquila

fonte dei dati: European-Mediterranean Seismological Centre


PROGETTO DI MASSIMA DI UN SISTEMA DI MONITORAGGIO PER LA PREVISIONE
SISMICA ED ADATTABILE SU QUALSIASI TERRITORIO AD ALTO RISCHIO SISMICO.
Ogniqualvolta si verifica una calamità naturale, ci si ricorda della fragilità del pianeta su cui
viviamo.
Ci si rende conto della drammaticità dell’evento, nel momento in cui si contano le vittime, si
raccolgono i feriti e si calcolano i danni ingenti.
Trascorsa la crisi, ci si dimentica dei disagi delle popolazioni colpite, dei danni ingenti, delle vite
umane perdute, del gran numero di invalidi che restano.
Alcuni dati degli ultimi forti terremoti:
16 Agosto 1999 11:11 terremoto M = 7.8 a IZMUT (Turchia), morti ≈ 10.000, Feriti 45.000.
26 Gennaio 2001 03:16 terremoto M = 7.9 Gujarat (India), morti > di 20.000.
21 Maggio 2003 18:44 terremoto M = 6.7 ALGERIA morti ≈ 2000, feriti 6782.
26 Dicembre 2003 01:56 terremoto M = 6.3 Bam (Iran), morti ≈ 40.000 (Bam distrutta)
24 Dicembre 2004 terremoto M = 9.2 Indonesia e area asiatica morti > di 180.000.
Costo dei danni incalcolabili.
La Sezione Ricerca della S.C.S. (Società di collaborazioni Scientifiche) ha realizzato nel 2002 un
sistema che permette di registrare allarmi sismici con un margine di anticipo sull’evento da 6 ore
a 24 ore. Il sistema è stato brevettato in Italia con n. 0001333804 a nome del Sig. Gioacchino
Giuliani.
L’installazione di un sistema di monitoraggio per la previsione di eventi sismici, fornito dalla
S.C.S., trova applicazione presso qualsiasi territorio terrestre ad alto rischio Sismico.
Il sistema si avvale dell’applicazione di una rete di rivelatori di gas Radon, distribuita sul territorio
o regione prescelta per il monitoraggio.
La rete di rivelatori del tipo PM4 e PM2, forniscono allarmi sismici con un margine di anticipo
sull’evento da 6 a 24 ore prima di un terremoto.
Il metodo utilizzato dai rivelatori del tipo PM4 e PM2, permette di rilevare il flusso medio del gas
Radon ed attraverso un’analisi del continuo dello stesso gas di rilevare quei segnali considerati
precursori sismici, che vengono prodotti da 6 a 24 ore prima di un evento sismico.
L’intera rete di monitoraggio sarà costituita da una serie di punti nodali distribuiti lungo il territorio
destinato al controllo della sismogeneticità della zona prescelta.
Ogni nodo dell’intera rete di monitoraggio, è costituito da almeno tre rivelatori del tipo PM4 e
PM2, posti ad una distanza fra loro pari ad un minimo di 25 Km ad un massimo di 80 Km.
La distanza tra i rilevatori nodali è determinata dalla natura del territorio e dalla sismogeneticità dei
luoghi da sottoporre a monitoraggio.
Ogni rivelatore di Radon, (PM4 o PM2) costituirà quindi una Stazione di rilevamento ed
acquisizione dati.
Le varie stazioni dei diversi punti nodali, saranno tra loro interallacciate e collegate ad una Sala
Sismica Centrale.
I dati forniti dalle singole stazioni e gli allarmi sismici provenienti dai diversi punti nodali,
convergeranno sui monitors della Sala Sismica Centrale.
L’allarme sismico prodotto dai punti nodali, sarà costituito da un messaggio e-mail contenente
l’informazione sulla stazione o stazioni che hanno inviato l’allarme, le coordinate della località
dove è previsto l’evento ed il valore del grado sismico richter previsto.
Un messaggio telefonico può essere inviato nello stesso tempo ad un utente addetto al controllo
degli allarmi.
L’eventuale allarme di grado richter considerato pericoloso per la popolazione della regione
interessata, sarà gestito da personale dell’ente preposto alla diffusione delle informazioni sugli
allarmi sismici.
Naturalmente il sistema di gestione degli allarmi sismici, prevede un periodo di specializzazione da
parte del personale preposto al controllo della sala sismica, al controllo degli allarmi delle singole
stazioni ed all’analisi dei dati provenienti dai vari punti nodali.
Ogni stazione di rilevazione, che dovrà ospitare un rivelatore di tipo PM4 o PM2, necessita di una
piccola stanza, interrata o semi interrata, senza ventilazione o possibilità di scambi di aria con
ambienti diversi da quello destinato all’alloggiamento del rivelatore.
Il locale in questione necessita di un punto presa per la corrente elettrica 220/250 Volts ed un punto
presa di collegamento alla rete internet.
Ogni rivelatore ha un ingombro medio di 70 cm x 70 cm x 70 cm.
Il peso del rivelatore può variare da 400 Kg a 500 Kg.
L’ingombro dell’elettronica di acquisizione ed analisi per una stazione sarà pari a 70 cm x 70 cm x
100 cm h. Peso dell’elettronica da 40 Kg a 60 Kg.
La potenza elettrica impegnata per l’alimentazione di una stazione, è pari a 1000 W.
La S.C.S. provvederà all’istallazione dei rivelatori, al controllo degli stessi ed all’allestimento della
Sala Sismica di controllo.
A titolo di esempio vengono di seguito riportate immagini del rivelatore, dell’elettronica e delle
mappe cartografiche raffiguranti il tipo di rete di monitoraggio.

 

Esempio di ingombro del tipo di elettronica utilizzato per una stazione di monitoraggio, per il
rilevamento della variazione di concentrazione di Radon.

Rivelatore di tipo PM2, utilizzato in una stazione di monitoraggio per la rivelazione del flusso
medio della concentrazione di Radon in ambiente.

Esempio di piccola rete di monitoraggio in Abruzzo. Costituita da n. 3 rivelatori di Tipo PM4 e
PM2. Le Stazioni operative rappresentate sono quelle di L’Aquila, Gran Sasso (Assergi) e
Pineto (TE)

Rivelatore di tipo Pm4 attualmente operativo presso la stazione di rivelamento Coppito (L’Aquila)

Esempio di una grande rete di monitoraggio, per la previsione sismica.
Studio effettuato dalla S.C.S. per un territorio ad altissimo rischio sismico, utilizzando mappe
sismogenetiche, relative ad una area attraversata da una grande faglia.

La relazione descritta è stata curata dai Sig.ri Gioacchino Giuliani e Roberto Giuliani, per conto
Della S.C.S.
Società di Collaborazioni Scientifiche
Via del Corso, 29/31
67110 L’Aquila – Italy
Tel. 00390862362813
Cell. 00393389798972.

Altro articolo

Due tecnici Aquilani ed un amico Russo, Giampaolo Giuliani, Roberto Giuliani
e Victor Alekseenko potrebbero cambiare il mondo e la nostra città
L'Aquila, 9 nov. - "...e ora passiamo alle previsioni sismiche; prevista una
leggera scossa sismica per domani pomeriggio nella zona di Los Angeles." Non
si tratta di una scena da film futuristico hollywoodiano, ma di una realtà
molto vicina, anzi, a L'Aquila già operativa. Si chiama "PM4", è una scatola
pesante mezzo quintale, semplice come tutte le opere di genialità,
realizzata e brevettata da Giampaolo Giuliani, un tecnico aquilano in
servizio presso i Laboratori Nazionali del Gran Sasso, e prevede i terremoti
in un ampio raggio d'azione e con un anticipo compreso tra le 6 e le 24 ore,
con una precisione molto alta. La eclatante notizia risale al 2003 e fu
anticipata dal nostro direttore Gianfranco Colacito sulle pagine
dell'Opinione, poi scese il silenzio... (segue)

"I tempi erano forse prematuri" - ci spiega Giuliani - "da quando Colacito
pubblicò quel bellissimo articolo, ho dovuto subirne tante". Ma andiamo per
ordine. Siamo andati a trovare Giampaolo Giuliani dopo aver appreso che
l'Istituto di Fisica Nucleare del Gran Sasso, pur non prendendo parte
direttamente alla sperimentazione, ha deciso di installare uno dei
macchinari "MP4" per testarne l'efficacia. Entriamo nel laboratorio privato
di Giuliani, Presso la Ditta S.C.S. di Coppito, "CHIOCCIOLANDIA", ed è come
ci aspettavamo che fosse: una specie di officina artigianale con macchinari
scintillanti di lucette e computer, non certo gli avveniristici laboratori
con gli effetti speciali o scansioni laser dell' iride per accedere nelle
segrete sale della scienza ufficiale. Solo lavoro, passione e genialità. La
prima domanda è logica: "Ha scoperto una macchina del genere e il mondo non
sa nulla? Funziona davvero o è una sperimentazione?" - "Altrochè se
funziona, e già da qualche anno! Non è semplice spiegare tutto. Non sono uno
scienziato, sono soltanto un tecnico che ha scoperto "l'uovo di Colombo". La
scienza ufficiale è convinta che non si possano prevedere i terremoti, e
potete immaginare a quanti incontri-scontri con l'elite scientifica ho
dovuto far fronte. Ma i dati e i rilevamenti che ho sono inconfutabili". -
La solita vecchia storia, dunque, quando per la Scienza non è piacevole
comprendere che si è insegnato per anni qualcosa di inesatto come ad esempio
la teoria tolemaica, e dover cominciare a ridisegnare teorie e studi. L'MP4
è stato ideato e costruito da Giuliani con l'aiuto del figlio Roberto, del
fisico russo Victor Alekseenkor e di pochi altri collaboratori ed è
brevettato a livello nazionale ed internazionale. La tecnica utilizzata si
basa sulla rilevazione del radon 222, un gas radioattivo, da cui il
Rivelatore evidenzia alcuni figli diretti della catena di decadimento
radioattivo (così li chiama Giuliani) e da cui un software specifico
seleziona i così detti Precursori Sismici, quando i precursori vengono
rivelati, un particolare algoritmo sempre via software provvede a lanciare
un Preallarme o un Allarme dalle 6 alle 24 ore prima dell'evento.
Con almeno tre rivelatori distanti fra loro un certo numero di Km. Sarà
possibile dare anche la prossimità dell'epicentro ed il grado sismico.
Abbiamo superato con successo tutti i test e in quattro anni abbiamo
previsto praticamente tutti i terremoti rilevati dall'MP4, con una
efficienza fuori dalle attese. La scoperta di Giuliani è consistente e di
portata mondiale, tanto è vero che Eduard Nixon, fratello dell'ex presidente
degli Stati Uniti, attualmente capo della Protezione civile di San
Francisco, ha chiesto finanziamenti al governatore per l'acquisto di 5
detector. Attualmente vi sono in funzionamento 2 rilevatori, uno a L'Aquila
ed uno a Pineto. Il team di Giuliani è stato inoltre il primo al mondo ad
effettuare rilevamenti di radon sul fondo del mare al largo delle coste
abruzzesi ottenendo risultati inattesi. E ora la nota dolente.
L'Aquila rischia di perdere ancora una volta una possibilità enorme
dal momento che nessuno è stato in grado di afferrarne le
enormi possibilità. In un territorio che perde ogni giorno posti di lavoro
si presenta una specie di miracolo che potrebbe portarci all'attenzione
mondiale. Giuliani da aquilano verace ed un po' campanilista vorrebbe che
sia L'Aquila e l'Abruzzo a beneficiare dei frutti che potrebbero scaturire
dalla scoperta, una volta che avrà ottenuto l'ufficialità anche dal mondo
scientifico. Il rischio di vedere trasportata altrove un'azienda che produca
migliaia di questi macchinari è alto.
In Italia sono già pronte molte città che farebbero salti mortali per
produrre il brevetto, oltre che, naturalmente, le grandi potenze americane,
russe e giapponesi. -

ROMA - Giampaolo Giuliani, il tecnico denunciato per procurato allarme dopo aver previsto che un terremoto di grande entità avrebbe potuto colpire la zona di L'Aquila,
era stato intervistato il 24 marzo anche dal blog DonneDemocratiche.it. Nell'intervista, dopo aver raccontato com'è nato il suo apparecchio, aveva corretto la sua prima previsione,
dicendo: "Mi sento di poter tranquillizzare i miei concittadini". In un'altra intervista, rilasciata ad una tv locale, ora visibile su You Tube e sul sito di Repubblica.it, era stato più determinato
e aveva detto che una rete di "precursori" potrebbe garantire una capacità predittiva di grande rilevanza umana e sociale.

Sulle sue dichiarazioni, che in queste ore fanno il giro dei mezzi di informazione, si sta scatenando una polemica, se è vero che l'Istituto di Geofisica e Vulcanologia ribadisce oggi in un comunicato ufficiale: "Si sottolinea la circostanza secondo la quale, allo stato attuale delle conoscenze, non è possibile realizzare una previsione deterministica dei terremoti (previsione della localizzazione, dell'istante e della forza dell'evento). Ciò è vero anche in presenza di fenomeni quali sequenze o sciami sismici che nella maggior parte dei casi si verificano senza portare al verificarsi di un forte evento. Una scossa quale quella che si è manifestata oggi viene normalmente seguita da numerose repliche, alcune delle quali probabilmente assai sensibili."

Nell'intervista a DonneDemocratiche, Giuliani specifica che il suo "precursore sismico" arriva ad anticipare il verificarsi di un terremoto fino a 6-24 ore prima.

La tecnica è basata sullo studio del "comportamento" dell'elemento chimico chiamato Radon. Ecco il racconto di Giuliani: "Nel 2001 stavamo osservando il misuratore di particelle cosmiche presso l'Istituto quando, in corrispondenza del terremoto in Turchia, rilevammo una quantità straordinaria, rispetto al solito, di radon. Così ho impiegato quasi 2 anni per realizzare da solo uno strumento in grado di rilevare il radon, iniziai ad osservarlo ed a studiarlo, e con l'aiuto di un sismografo mi resi conto che la concentrazione di radon aumentava in corrispondenza di un evento sismico. Nel 2002, ad esempio, in corrispondenza del terremoto di S. Giuliano, registrammo valori 100 volte maggiori alla norma, ma disponendo di 1 solo precursore sismico eravamo in grado di emanare un allarme per un evento sismico che distava più di 50 km da L'Aquila, senza poter fornire altre informazioni circa la collocazione o la direzione dell'evento stesso. Oggi con 5 precursori saremmo in grado di essere molto più precisi, triangolando i dati ed i segnali di concentrazione del radon"

Dopo aver fornito elementi utili sulle caratteristiche del Radon, Giuliani dice che il suo apparecchio non si trova nel campo della teoria, me ne esistono già cinque esemplari:" "I 5 Precursori sismici si trovano a Coppito, nel Laboratorio del Gran Sasso (ospite dell'INFN), presso la scuola De Amicis, a Fagnano e a Pineto; sono tutti a più di 3 metri sotto terra e in corrispondenza di un evento sismico rilevano nello stesso momento, lo stesso segnale creando un grafico perfettamente sovrapponibile."

Poi Giuliani entra nel merito degli allarmi che si sono verificati in Abruzzo e della sua previsione: "I dati ottenuti in questi 9 anni di studi, ci hanno consentito di rilevare un rischio sismico maggiore nel periodo invernale che va da novembre ad aprile. Senza voler banalizzare, ma per semplificare i concetti, posso aggiungere anche che l'attività sismica è strettamente correlata alle fasi lunari. In particolare quest'anno, il sistema Terra-Luna, si è venuto a trovare al Perielio (punto più vicino al Sole, in inverno) con la Luna nello stesso periodo alla minima distanza dalla Terra, e con il Pianeta Venere allineato, in fase di Venere piena anch'essa vicina. L'attrazione gravitazionale delle masse sulla Terra hanno intensificato l'effetto marea sul nostro pianeta, rendendo gli eventi sismici più rilevanti, rispetto agli altri sciami, cui siamo stati interessati negli anni precedenti. Mi sento di poter tranquillizzare i miei concittadini, in quanto lo sciame sismico andrà scemando con la fine di marzo."

Ma evidentemente il metodo non è ancora accettato uniformemente. Giuliani descrive così l'atteggiamento della comunità scientifica: "Mi osserva con interesse ... una parte mi da fiducia, come dimostra anche il Direttore dell'INFN che mi ha messo a disposizione un locale per ospitare uno dei 5 Precursori sismici; una parte è un po' più cauta e scettica."

GIAMPAOLO gIULIANI RIBADISCE: LA DENUNCIA: «SI POTEVA CAPIRE» - Oggi, dopo la tragedia, Giuliani parla con amarezza: «C'è il rischio che domani mi mettano in galera - dice - ma confermo: non è vero, è falso, che i terremoti non si possono preveder. Sono 10 anni che noi riusciamo a prevedere eventi di questo tipo in una distanza di 100-150 chilometri da noi. Da tre giorni - continua- vedevamo un forte aumento di radon, al di fuori della soglia di sicurezza. E forti aumenti di radon segnalano forti terremoti. Questa notte il mio sismografo denunciava una forte scossa di terremoto e ce l'avevamo online. Tutti potevano osservarlo e tanti l'hanno osservato. Poteva essere visto ce ci fosse stato qualcuno a lavorare o si fosse preoccupato. Abbiamo vissuto la notte più terribile della nostra vita, sono sfollato anche io... Questi scienziati canonici, loro lo sapevano che i terremoti possono essere previsti».

Altro Articolo:

  Ecco il tracciato del Sisma che ha colpito l'Aquila, di magnitudo 5.8

La Rete Sismica Nazionale dell’INGV ha registrato un terremoto di Magnitudo 5.8 (Magnitudo Richter) (6.2 Mw=magnitudo momento) nella zona dell'Aquilano, il 6 Aprile 2009 alle 3:32 (ora italiana). Le coordinate epicentrali risultano le seguenti: Lat. 42.33N e Long. 13.33 E. La profondità dell'ipocentro è pari a 8.8 km. Il terremoto è caratterizzato da un meccanismo di tipo estensionale, con piani di faglia orientati NW-SE e direzione di estensione NE-SW (anti-appenninica)".


Nell'articolo Portali-addendum, si legge: "Il Gran Sasso, cui accennai in Portali (prima parte), si trova più o meno sulla stessa longitudine (13 gradi e 34 primi) di Lipsia (12 gradi e 23 primi). In laboratori ricavati nel cuore della montagna, non si studierebbero solo neutrini, ma si compirebbero esperimenti segreti. Alcuni scalatori hanno affermato di aver colà notato anomalie".

La longitudine del sisma coincide quasi perfettamente con quella in cui è ubicato il laboratorio sotterraneo. Ora, con i dati a nostra disposizione, non è possibile stabilire se la scossa sia stata naturale o artificiale. Stando ad alcuni ricercatori, un terremoto generato da sistemi d'arma come H.A.A.R.P. è contraddistinto da una frequenza più alta.

Secondo la scienziata statunitense Rosalie Bertell, armi elettromagnetiche “sono in grado di causare terremoti in siti scelti come bersaglio, sprigionando energie equivalenti alle più forti esplosioni nucleari”.

E' possibile che il fenomeno tellurico sia stato prodotto da movimenti delle faglie: infatti "l’attività di questi giorni si colloca tra la terminazione meridionale della faglia che si è attivata nel terremoto del 1703 (Int. MCS del X grado MCS, pari a Magnitudo circa 6.7) ed i limiti settentrionali della faglia associata nei cataloghi al terremoto del 1349 e di quella denominata 'Ovindoli-Piani di Pezza'”.

Comunque sia, è sintomatico che la "scienza" ufficiale si ostini a dichiarare che i sismi non possono essere in nessun modo preveduti. Giampaolo Giuliani ed altri dissentono. (Si ascolti l'intervista che riporta you tube sull'argomento).

E' inammissibile che tristi personaggi come Guido Bertolaso osino denunciare per procurato allarme (sic) chi avverte circa possibili pericoli. Cassandra andrebbe ascoltata, mentre gli incompetenti (ammesso e non concesso che si tratti solo di incompetenza) meritano inesorabile riprovazione.

Ci chiediamo a che cosa servano i tanto decantati "progressi scientifici", se ancora non si riesce a predire con una buona approssimazione eventi tanto disastrosi. Come spesso avviene, quasi tutte le innovazioni riguardano gli apparati militari. Forse evitare che un terremoto provochi danni e vittime non è nell'interesse dei governi: alcuni geologi d'altronde, invece di adoperarsi per monitorare la situazione sismica, invece di compiere studi utili per la collettività, preferiscono indagare le scie di "condensazione".

Secondo Angelo Ciccarella, la città di Lipsia, nell'ex Germania est si può considerare un portale.
Nell'articolo intitolato La faccia occulta del Comunismo, lo studioso riferisce di audaci e pericolosi esperimenti compiuti dalla Stasi,
la famigerata polizia segreta della Repubblica "democratica" tedesca, nel centro militare di Lipsia.

L'autore scrive: "Ad un certo punto degli esperimenti collettivi, qualcosa avvenne, qualcosa di terribile si rivelò.
Un arco di radiazioni elettromagnetiche in espansione saliva attraverso gli strati delle vecchie e mai perdute trasmissioni radio di tutto il mondo.
Mentre l'arco procedeva in avanti nello spazio, viaggiava indietro nel tempo e tutti i sensitivi, i tecnici, gli agenti Stasi presenti nel laboratorio,
ascoltarono voci di persone morte da tempo, fino a tornare verso la prima trasmissione di Marconi e... verso il silenzio.
Poi, di botto, una cascata inarrestabile di fenomeni: onde scalari, anomali flussi energetici, ampio spettro di frequenze E.M. e nota di fondo dominante,
fragorosa. Uno squarcio nel vuoto ed ecco comparire un grande arazzo sopra le teste dei presenti. Vi era raffigurato un animale, una salamandra in mezzo alle fiamme.
Si era creato un punto cieco, un varco dimensionale. Inavvertitamente, pigiando un tasto, oltre lo specchio si era acceso un led. Le facoltà extrasensoriali
di uomini e donne unite ad un glifo psicoacustico ed ad un codice alieno (?) avevano trapassato le Porte del Tempo. La Quarta dimensione non è impenetrabile.
Quando si lacera la sua struttura connettiva, ne sopravviene un universo tangente, instabile, limitato nella trasmissione. Esisterebbero così infiniti universi, nei quali si svilupperebbero, in maniera indipendente, tutte le alternative possibili, a partire da una situazione che ha la possibilità di evolvere in diversi modi... Ciascuna alternativa diventerebbe così realizzata e praticabile in un mondo a sé, senza essere cancellata dalle altre nelle vicende di un unico mondo... Lipsia era un cancello stellare che attendeva da tempo di essere aperto, un dispositivo in grado di viaggiare nel tempo e nelle dimensioni. I massimi esperti di fisica e di elettronica provarono inutilmente ad usare il portale come arma. Il portale da solo non serve a nulla; per funzionare ha bisogno di un suo omologo di uscita nel segmento spazio-temporale di arrivo. Indovinate dove si trova l'altra porta e soprattutto chi la detiene?"


Ciccarella intravede dietro le quinte della storia una cospirazione di proporzioni gigantesche: il nome in codice dell’Entità che agisce nell’ombra è Salamandra.
La salamandra, animale che, nella tradizione, era in grado di vivere nel fuoco, ricorda le figure create dal plasma.

Nel suo testo criptico, non privo di autocompiacimento e di qualche tono da ierofante, Ciccarella ha comunque il merito di rispolverare gli inquietanti scenari di un conflitto segreto, il conflitto delle psicospie. Durante gli anni della Guerra fredda, agenti psichici statunitensi (famoso è lo scienziato ed artista Ingo Swann) e dell'est europeo si confrontarono in esperimenti di visione a distanza ed in altri ambiti relativi all'uso a fini spionistici dei poteri mentali. L'ipotesi dell'autore pare la seguente: forze extrasensoriali furono all'origine di un'apertura nella sfera spazio-temporale. E' probabile che ai risultati concorsero diaboliche tecnologie. La convergenza tra dimensione medianica e tecnologica sembra essere all'origine delle scoperte che resero celebre Marconi; un ascolto di esseri dell'altrove illuminò la mente di Tesla. Gli scientisti reputano tale approccio ridicolo, poiché essi distinguono tra la scienza e la superstizione. Sappiamo, però, quanto le ricerche di eminenti studiosi anche razionalisti siano venate di magia, di occultismo.

La riflessione sui portali, su questi aditi verso regioni invisibili ci conduce al confine della percezione ordinaria, là dove i concetti razionali si appannano per compenetrarsi con idee controintuitive, con i paradossi delle verità eccentriche ed inaccettabili, appena rasentate da intelletti sublimi. Quando ci accostiamo a certe verità mai attinte prime, esse si sbriciolano e si depauperano, a contatto con il linguaggio e la logica.

Dove potrebbe essere l'altra porta verso l'Ade? Forse nelle viscere del Gran Sasso, cui accennai in Portali, visto che si trova più o meno sulla stessa longitudine (13 gradi e 34 primi) di Lipsia (12 gradi e 23 primi)? In laboratori ricavati nel cuore della montagna, non si studierebbero solo neutrini, ma si compirebbero esperimenti segreti. Alcuni scalatori hanno affermato di aver colà notato anomalie.

Nulla è certo; tutte queste informazioni e voci vanno prese con il beneficio del dubbio, ma è evidente che la realtà è molto più enigmatica e paradossale di quanto si possa immaginare.

Cronaca di una tragedia annunciata

News di Attualità | 7 Aprile 2009 | Commenta l'articolo
Lo spaventoso terremoto che ha colpito l’Abruzzo, sventrando le case e mietendo centinaia di vittime, non è arrivato all’improvviso. Già nelle scorse settimane in tutta la regione c’erano state tante piccole scosse sismiche ripetute: delle avvisaglie che qualcosa di molto più devastante era in arrivo. Proprio in virtù di questo - e basandosi su un suo metodo sperimentale - il tecnico di laboratorio dell’Istituto Nazionale di Fisica Nucleare (Infn) Gioacchino Giampaolo Giuliani aveva dato l’allarme preannunciando una catastrofe. Ma non gli avevano creduto, dandogli della Cassandra. Non solo: era stato anche denunciato per “procurato allarme” e il capo della Protezione Civile Guido Bertolaso lo aveva apostrofato come “uno di quegli imbecilli che si divertono a diffondere notizie false”. Notizie rivelatesi, poi, tragicamente vere.

Le precisazioni dell’Inaf. L’Inaf precisa che “Gioacchino Giampaolo Giuliani è un collaboratore tecnico non laureato dell’Istituto di Fisica dello Spazio Interplanetario di Torino, che è una delle venti strutture Inaf” e che l’ente “ha come missione e scopo lo studio dei fenomeni che avvengono nello spazio e nell’universo e non di terremoti o altri fenomeni legati alla geofisica”. L’Istituto precisa inoltre che “il signor Giuliani lavora, come collaboratore tecnico, presso i Laboratori Nazionali dell’Istituto Nazionale di Fisica Nucleare del Gran Sasso per conto di Ifsi-Inaf, nel quadro di una collaborazione multipartner nell’esperimento Lvd (Large Volume Detector) per la rivelazione di neutrini prodotti da collasso gravitazionale stellare”. Pertanto, conclude la nota, “le attività del signor Giuliani rispetto la presunta possibilità di previsione di terremoti non sono una ricerca Inaf, ma vengono svolte dal Giuliani stesso per scopi personali al di fuori dell’orario di servizio per l’Istituto”. Una cosa, però, l’Inaf non la dice: perché non è stato dato alcun credito (né supporto) alle parole di Giuliani, che pure aveva portato dati circostanziati sugli eventi futuri?

La solidarietà su Facebook. Dopo la tragedia del terremoto in Abruzzo e le nuove polemiche per le dichiarazioni di Gioacchino Giampaolo Giuliani, che avrebbe messo a punto un sistema, ancora sperimentale, per prevedere un sisma sul social network Facebook fioccano i gruppi che solidarizzano col tecnico. In molti addirittura propongono che diventi il nuovo capo della Protezione Civile. Sono nati i gruppi su «Giuliani eroe nazionale», «Chiedere scusa a Giuliani», «Solidarietà per Giuliani», «Giuliani aveva ragione sul terremoto». Mentre conta già 2151 iscritti il gruppo «Sosteniamo Giampaolo Giuliani ed il brevetto per la previsione di terremoti».

Ingv: “Previsioni con radon fondate, ma non pronte”. Prevedere i terremoti basandosi sull’analisi del radon è una possibilità allo studio dell’ Istituto Nazionale di Geofisica e Vulcanologia (Ingv), «ma questo strumento di previsione non è ancora pronto». Federica Quattrocchi, responsabile dell’Unità funzionale di geochimica dei fluidi stoccaggio geologico e geotermia dell’Ingv, lavora da quasi vent’anni sui segnali che annunciano i terremoti. «Ho promosso questa ricerca in Italia nel 1990, ma spesso - rileva - il nostro lavoro è stato affossato proprio da persone come Giuliani, che si muovevano ai margini della comunità scientifica e si rivolgevano direttamente ai sindaci lanciando allarmi per far evacuare i Comuni».

Del sistema di previsione proposto da Gioacchino Giampaolo Giuliani, il collaboratore tecnico dell’Istituto Nazionale di Astrofisica che sostiene di avere previsto il terremoto, Quattrocchi dice: «Non contesto la sua ricerca, ma il metodo: si muove a margine della comunità scientifica e sostiene che le previsioni siano deterministiche». La verità, prosegue, “è che questo sistema di previsione non è ancora pronto e attualmente non è in grado di fare previsioni esatte, ossia di stabilire la data e il luogo esatto in cui avverrà un terremoto. Per avere dati affidabili è necessaria una rete nazionale di rivelatori di radon, con una decina di strumenti per regione e dal costo complessivo di circa un milione di euro. In pratica, se la rete venisse finanziata ora, i primi risultati potrebbero arrivare dopo il 2050″.

Di Pietro: “La tragedia poteva essere contenuta”. «Oggi è soltanto il giorno del dolore, della solidarietà e anche della vicinanza umana e politica verso le istituzioni, Governo compreso, che si stanno prodigando per venire incontro alla popolazione duramente colpita dal sisma», afferma Antonio Di Pietro sul suo blog. «In Italia - aggiunge poi il leader di Idv - ‘ambasciator porta pena’, anzi rappresenta egli stesso la pena. Se denunci sei denigrato, sospeso, denunciato a tua volta e, come avviene negli Stati a diritto compromesso, ne paghi le conseguenze. È accaduto a Giampaolo Giuliani, fisico e tecnico presso il laboratorio nazionale del Gran Sasso.

Quella che è accaduta a L’Aquila - sottolinea Di Pietro - è una tragedia che poteva essere contenuta. Giampaolo Giuliani, che aveva preannunciato questo disastro, è stato denunciato per procurato allarme. Sarebbe stato il caso di verificarne accuratamente la segnalazione, attivando le precauzioni necessarie per evitare l’irreparabile o, comunque, contenerne le conseguenze. L’Italia, come tutti sanno, è un Paese ad altissima attività sismica. Tra gli edifici crollati in Abruzzo ci sono anche scuole, diversi palazzi pubblici ed una Casa dello studente. Mi domando - conclude Di Pietro - che senso abbia investire nel progetto del ponte sullo Stretto di Messina, esso stesso in zona ad alto rischio sismico, o discutere sulla costruzione di centrali nucleari sicure, se non si riescono neppure a garantire i requisiti antisismici nelle abitazioni».

DOMANDA: La cosa che non mi torna è questa: se i terremoti non si possono prevedere (cosa che mi pare realistica) come mai i Bertolaso e i vari altri laureati una settimana fa prevedevano che non ci sarebbe stato? Sarà mica che la prevedibilità funziona solo in una direzione??!
 

Poi c'e' un video sul test per la prevenzione sismica fatto con l'universita dell'aquila, dietro alla scuola di san berardino, dove si procede a trivellare un foro nel terreno sitiato su un angolo della piazza, per trenta metri di profondità. per una prova sismica, dove sarà prodotta una onda sismica artificiale, per vedere come il sottosuolo dell'aquilano si comporta.

Se è vero che un sisma è sempre accompagnato da precursori sismici percepibili da strumenti quali il PM4, SE FOSSE PROVATO CHE UNA SCOSSA SISMICA NON LO FOSSE STATA, RISULTEREBBE EVIDENTE LA REALTA' DI UN'ALTRA MOTIVAZIONE  ESOGENA DEL FENOMENO OSSERVATO.  Nell'evenienza di una tale circostanza, sarebbe opportuno indagare e verificare quale sistema di onde potrebbero produrre il sisma anomalo caratterizzato dall'assenza dei suddetti precursori sismici.

Qualora una forma d'onda "X" fosse all'origine di una relativa scossa "Y no-radon", basterebbe indagare su tre diversi campi ondulatori di frequenze per scoprire quale origine possa avere. Analisi di spettro dei campi Elettrico, Magnetico o Gravitazionale terrestre offrirebbero subito una risposta chiara. Come ad ogni spinta a carica risonante segue una oscillazione relativa proporzionale alla somma delle quantità addizionate, cosi la forza liberata di un'onda smorzata di questo tipo, risuonerebbe come un gong, e sebbene sarebbe caratterizzata dall' assenza di precursori radon, sarebbero facilmente riscontrabili tracce della sua costruzione a livello vlf e xlf, nello spettro magnetico e quello gravitazionale.

MISURE DI RADON E SISMICITÀ NELLA ZONA DI CONFINE TRA ITALIA E SLOVENIA

 

GNGTS – Atti del 23° Convegno Nazionale / 05.10
M. Santulin (1), A. Riggio (1), S. Sancin (1), A. Popit (2) e J. Vaupotic (2)
(1) Istituto Nazionale di Oceanografia e Geofisica Sperimentale - OGS, Trieste, Italia
(2) Jozef Stefan Institute, Ljubljana, Slovenia
MISURE DI RADON E SISMICITÀ NELLA ZONA DI CONFINE
TRA ITALIA E SLOVENIA
Misure della concentrazione di gas radon a scopo sismologico sono state
eseguite in Italia (Friuli) e in Slovenia rispettivamente, dalla fine del 1995 e dal 1998.
Le tecniche di acquisizione, come anche il tipo di strumento sono differenti, ma i dati
sembrano comparabili. Una collaborazione scientifica tra l’OGS e l’Istituto Jozef
Stefan di Lubiana, iniziata nel 2000, è stata ufficializzata con un progetto bilaterale
Italia-Slovenia nel 2003. In Slovenia sono in funzione tre stazioni di radon in acqua
termale, ma in questo lavoro sono stati utilizzati solamente i dati registrati dalla
stazione di Bled perché quest’ultima ha avuto una maggiore continuità di
registrazione nel periodo considerato. Nello stesso sito vengono misurate anche la
temperatura e la conducibilità elettrica dell’acqua. Lo strumento utilizzato è una
sonda Barasol (MC-450, Algade, France) con sensibilità di 50 Bq/m3. L’intervallo di
campionamento è di 1 ora.
Il sito italiano è equipaggiato con uno strumento a registrazione continua del
tipo a cella a scintillazione di Lucas, modello Prassi 5S della Silena. La sensibilità è
di 4 Bq/m3 e l’intervallo di campionamento utilizzato è di 3 ore. La misura è eseguita
nell’aria, all’interno di un pozzo dove il livello medio della falda è di circa 17 metri. A
pochi metri di distanza una stazione meteorologica provvede alla registrazione della
pressione atmosferica, della pluviometria, della temperatura dell’aria e dell’acqua e
del livello della falda.
Nel maggio 2004, a circa 30 metri di distanza dal sito italiano di registrazione, è
stata installata un’altra sonda Barasol per l’acquisizione di radon in suolo. I dati
registrati dai due strumenti di tipo diverso sono stati analizzati e confrontati. I primi
risultati sembrano confermare che dati registrati con metodologia e strumentazione
differenti presentano un andamento simile (Fig .1).
La concentrazione di radon è analizzata contemporaneamente ai dati
sismologici registrati dalla Rete Sismometrica del Friuli, gestita dall’OGS, e dalla
Rete Sismometrica Slovena, gestita dall’Agenzia dell’Ambiente della Repubblica
Slovena. Per i terremoti con epicentro al di fuori dell’area in esame, sono utilizzati i
dati riportati in cataloghi internazionali.
La selezione dei singoli terremoti è stata fatta secondo la relazione di Hauksson
e Goddard (1981), che determina la magnitudo minima necessaria perché venga
registrata una variazione nella concentrazione di radon ad una data distanza
dall’epicentro.
Le variazioni della concentrazione di radon a lungo periodo sono state, invece,
confrontate con le variazioni nel tempo del parametro b della legge di Gutenberg e
Richter (1944), calcolato utilizzando tutti i terremoti con epicentro nell’area
considerata e con magnitudo maggiore o uguale a 1,8.
I risultati hanno mostrato il legame tra l’andamento della concentrazione di
radon e lo stato di stress presente nella zona. Tale relazione è evidente, nei dati
relativi al periodo di preparazione dei due terremoti forti avvenuti nel 1998 e nel 2004
a Bovec, in Slovenia.


GNGTS – Atti del 23° Convegno Nazionale / 05.10
Fig. 1 - Confronto fra le concentrazioni di radon registrate in giugno 2004 dalle due diverse
strumentazioni, nel sito di Cazzaso (Friuli). I valori registrati con il Barasol sono stati divisi per 5 per
una migliore rappresentazione grafica.
BIBLIOGRAFIA
Gutenberg B. and Richter C.F.; 1944: Frequency of earthquakes in California. Bull.Seism.Soc.Am. 34,
185-188.
Hauksson E. and Goddard J.G.; 1981: Radon earthquake Precursor studies in Iceland.
J.Geophys.Res. 86, 7037 – 7054.

Geophysical Research Abstracts,
Vol. 10, EGU2008-A-10334, 2008
SRef-ID: 1607-7962/gra/EGU2008-A-10334
EGU General Assembly 2008
© Author(s) 2008
Radon precursors to earthquakes: a search for physical
mechanisms
H. Woith (1), G. Martinelli (2), J. Heinicke (3) & H. Surbeck (4)
(1) GeoForschungsZentrum Potsdam, Germany, (2) Servizio Cartografico e Geologico,
Regione Emilia-Romagna, Bologna, Italy, (3) TU Bergakademie Freiberg, Germany, (4) Swiss
Center of Hydrogeology, University of Neuchatel, Switzerland (radon@gfz-potsdam.de / Fax:
+49 331 2881204)


We investigated a series of publications where radon anomalies were reported to precede
earthquakes with special emphasis to the proposed physical mechanisms. The
aim of our study is not an in extenso listing of all reported radon anomalies but rather
an attempt to elaborate certain classes of radon anomalies and their physical origin. By
far the most popular attempt to relate radon anomalies with earthquakes was – and still
is – a 1:1 approach, i.e. one radon anomaly is related to one earthquake. Whereas the
anomaly definition is more or less straightforward, the selection of the “corresponding”
earthquake is not. In fact, it is the most critical, sometimes arbitrary procedure.
Other types are 1:n and n:n relations describing one radon anomaly in relation to the
seismicity of an area or a time segment and a series of radon transients in relation
to seismicity in time and/or space, respectively. Finally, we discuss 1:0 types of relations,
i.e. a radon anomaly which is not followed by an earthquake. Needless to say,
that the latter type is rarely published. We will present examples that non-seismically
originated radon anomalies can be strikingly similar to the “real” thing and thus the
underlying physical processes might be strikingly similar as well.

E in Thailandia vediamo cosa scrivono al riguardo:
Soil gas radon as an earthquake precursor: some considerations on data improvement

K. Wattananikorna, M. Kanareea and S. Wiboolsakea

a Department of Physics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand

Received 6 June 1998;
accepted 24 July 1998.
Available online 28 January 1999.

Abstract: A data processing technique based on impulse responses from multivariable time series was applied to radon data. It was found that soil gas radon at 50 cm depth was affected more by meteorological parameters than at 100 cm depth. This fact together with the fact that the radon change at the latter depth was larger when an earthquake occurred, radon measurement at 100 cm depth was tested as an earthquake precursor. The result showed correlation between radon anomalies and earthquake events.

Article Outline

1. Introduction
2. Geology of the study area
3. Meteorological and earthquake data collection
4. Radon measurement
5. Theoretical consideration
5.1. Movement of radon in overburden
5.2. Radon data processing
6. Result and discussion
6.1. Improvement of soil gas radon data
6.2. Radon data as an earthquake precursor
7. Conclusion
Interessante!!! Non lo sanno neppure loro. Vediamo se anche altri ignorano la cosa... Ora andiamo in India...
Lascio l'articolo in Inglese perché e' bene conoscere almeno tre lingue oggigiorno e un poco di sforzo no fa mai male.

COMMUNICATIONS
1288 CURRENT SCIENCE, VOL. 86, NO. 9, 10 MAY 2004
A study of fluctuation in radon concentration behaviour as an earthquake precursor
Kaushik Majumdar
Electronics and Communication Sciences Unit, Indian Statistical
Institute, 203, B.T. Road, Kolkata 700 108, India

Before, during and after many large earthquakes significant
change in concentration of radon has been
observed in China, Japan and India. To establish this
behaviour as a potential earthquake precursor, it is
important to understand the phenomenon in terms of
geodynamics. Here I propose a geodynamic model to
establish a correlation between precursory radon concentration
fluctuations and an impending earthquake.
The model has been validated by field data collected
from Japan and China. Further investigation with this
model also helps to find a possible explanation for
abnormal variation in the geomagnetic field intensity
observed before, during and after some earthquakes.
LARGE earthquakes (M 5.5 or more) are the most devastating
of all natural calamities to human life and property.
Therefore, scientific research leading to successful prediction
of large earthquakes, sufficiently in advance, has
great significance to mankind. Unfortunately, to this date,
no reliable method for successful earthquake prediction
with precise timing, location and intensity has been developed.
Nevertheless, based on long-time studies, some
earthquake precursors like change in ion concentration in
water, variation in concentration of He, Ne, Ar, Rn and
N2 in the environment of the affected zone, abnormality
of behaviour in some animals, occurrence of milder foreshocks
before a large earthquake, sudden water-level change
in some wells, ground deformation, stress build-up in the
ground rocks (which may in turn alter electric resistance
of the rocks)1, etc. have been identified. But most of these
precursors behave erratically and therefore have been
poorly understood so far and earthquake prediction has
become a controversial issue2,3.
To understand the exact relationship of a precursor
with an impending large earthquake, it is essential to
know the geodynamical reason behind the occurrence of
that precursor. Erratic change in radon concentration has
been observed in many earthquake-prone zones a few
months before, during and after a large earthquake1,4–6.
Such behaviour has been observed in deep wells where
inducement of radon concentration fluctuations from the
environment can virtually be ruled out. Therefore, it is
tempting to consider a sudden erratic fluctuation in radon
concentration, for days on end, particularly in deep wells
in an earthquake-prone zone, as a potential omen for an
impending large earthquake. Here an analysis of the precursory
radon concentration behaviour during the Izu-
Oshima-kinkai earthquake (M 7.0) during 1978 in Japan
will be undertaken, and based on this a geodynamic model
will be proposed, which has been validated by data collected
during some other earthquakes.
Let us try to find an explanation for the erratic fluctuations
in radon concentration recorded in an artesian well
in Izu peninsula of Japan, before, during and after the
14 January 1978 Izu-Oshima-kinkai earthquake (M 7.0,
Figure 1)4.
Most of the rocks in the earth’s crust, as well as soil
and alluvium derived from it, contains uranium. About
99.3% of this uranium7 is the isotope U-238. Decay of an
U-238 atom marks the beginning of a series of 14 decays
that ends at the isotope Pb-206. Rn-222 (half-life 3.82
days, directly derived from Ra-226) is an intermediate
radioactive, inert, odourless, colourless, gaseous daughter
isotope in the decay7 of U-238 to Pb-206. Because of the
short half-life of Rn, it is generally believed that most of
the naturally occurring Rn in the environment comes
from only within a few metres of depth in the earth’s
crust7. The continental crust has an average thickness of
30 km, reaching its highest value of 70 km in the Tibetan
plateau8, and minimum at the spreading plate boundaries.
A standard argument for the variation in dissolved gas
concentration in water wells before, during and after an
earthquake in that zone, is the variation in release of the
gases entrapped in crustal rocks due to pore collapse and/
or crack build-up caused by variation in stress5. Being an
inert gas, radon concentration cannot be varied due to
chemical reaction with other substances exposed due to
crack build-up in the rocks. Normally, its source on the
earth is U-238 in the near-surface crustal rocks, whose
typical concentration level is 1 to 4 parts per million (ppm)7.
To explain the erratic precursory radon concentration
behaviour4 in Figure 1, I propose the following model.
e-mail: kaushik_r@isical.ac.in

RESEARCH COMMUNICATIONS
Figure 1. Precursory radon concentration fluctuations of 1978 Izu-
Oshima-kinkai earthquake observed in an artesian well, 25 km away
from the epicentre. The heavy dotted line shows the average annual
variation calculated from data for eight years. (After Wakita et al.4).
RESEARCH COMMUNICATIONS
CURRENT SCIENCE, VOL. 86, NO. 9, 10 MAY 2004 1289
The pattern of convection beneath the earth’s surface
plate, on which the plate tectonics and therefore earthquakes
depend, remains one of the main unsolved problems
in geophysics9. Because both seismic earth models
and mineral physics studies suggest that any present-day
chemical density heterogeneity in the mantle cannot exceed
1–2%, the mantle is often assumed to be homogeneous
for the purposes of dynamic modelling9. Clearly, this is
only an approximation. The upper mantle, near the subduction
zones, where subducted lithospheres with strong
chemical heterogeneity on the scale of tens of kilometres8
are slowly being recycled back into the mantle, must bear
the signature of this heterogeneity. In fact, chemical
segregation is favoured over homogenization by kinematic
mixing8.

In laboratory experiments, under much more simplified
conditions compared to the actual mantle, ‘doming’ regimes
have been observed at a layer of viscosity stratification in
a thermochemical flow model9. The Rayleigh numbers of
the experiments are comparable to those of the actual
mantle flow. Hot plumes come out from a small neighbouring
area centering the uppermost point of the dome.
In a recent study, four high-viscosity stratification layers
in the mantle have been identified, the uppermost being
only 220 km deep10. Chemical heterogeneity near a subduction
zone, coupled with this viscosity stratification,
should also give rise to smaller (only of the order of a
few thousand cubic km), shortlived blobs (existing only
for a few tens of years before being mixed-up in the
mantle), unlike the megablobs in existence for hundreds
of million of years, as inferred by Davaille9. Hot, cylindrical,
axisymmetric plumes of diameter only from a few
hundred metres to a few kilometres will reach up to just
below the plates. There is a region of upper mantle, 50 to
150 km deep, with low viscosity10 and low seismic velocity8
(i.e. with low density) over which plates glide, like
rafts floating on sea. Hot-plume injection in this region
near a subduction zone will expedite plate motion, by
reducing viscosity and/or increasing negative buoyancy
(by density reduction), which may lead to the advancement
of an earthquake. In this regard, it is interesting to
note the excellent correlation observed between seismic
activity of an earthquake-prone zone with artificial injection
of water in the rocks1, although the two mechanisms
and effects are completely different.
Studies carried out10, clearly indicate that the flow inside
the mantle is turbulent, in the sense that (i) the
velocity field is random and (ii) vortices occur chaotically.
Turbulent flow and temperature variations in the
upper mantle may occasionally increase the pressure
inside the blob, leading to a sudden increase in plume
discharge from it. This, in turn, will increase the plate
motion sometimes (but not always) leading to an earthquake.
Alternatively, the conduit through which the plume
is upwelling may get choked for sometime, leading to a
pressure build-up inside the blob. When the pressure
level becomes sufficiently high the conduit may burst,
much like a blocked pressure cooker valve, ejecting
greater amount of plumes than usual for sometime, till
the pressure level inside the blob returns to normal. This
process may expedite the plate motion more than it has
initially slowed down depending on the duration of
choking and the magnitude of the upwelled magma after
the blast.

To explain the precursory erratic radon concentration
behaviour of Izu-Oshima-kinkai earthquake (Figure 1),
let us note that the high annual average radon concentration
cannot be due to a source from near the surface of
the earth, but rather from much deeper underground. It is
actually coming out of the Moho through a conduit
across the plate. The Moho in this region got its supply
(steady at least for several years) of U-238 from a blob
just underneath, through upwelling plumes. Just a few
months before the 14 January 1978 event, the magma
upwelling conduit above the blob got choked-up, leading
to an intermittent release of smaller-than-usual amount of
uranium. (The uranium concentration is diluted fast through
mixing in the low viscosity and low density region,
where the velocity of ejection is much higher than the
usual convection speed in that layer.) The average concentration
of U-238 in the interior of the earth8 is 30.8 ×
10–9 kg/kg and that at the uranium ore7 is 10–3 kg/kg, as
characterized by sudden step-like lowering in the rate of
radon emission in Figure1. When the burst occurred, it
cleared the block at the magma upwelling conduit leading
to a greater-than-usual uranium supply at the Moho just
above, which in turn increased the radon emission rate
just before and after the 14 January event, as shown in
Figure 1. The hypothesis is further bolstered by the data
presented in Figure 2 in conjunction with the abovementioned
radon emission rate (see Figure 2 a). In Figure
2 b, water-temperature variation in a 500 m well, at
an epicentral distance of 30 km, has been recorded. In
Figure 2 c and d water-level variation at the well of
(Figure 2 b) and the ground strain variation (at a borehole
strain-meter placed at an epicentral distance of 50 km)
have respectively been recorded. From the data of Figure
2, it is apparent that the water in the well (Figure 2 b) is
directly being heated by hot gas emanating from the
plume under the plate and reaching the well by diffusion
through small cracks across the plate (unfortunately, no
chemical analysis data of the well-water at that time is
available to further confirm or reject this assumption). As
the plate slowly subducts, the ground strain is relaxed
leading to crack enlargement, which decreases the water
level in the well (by 0.3 m) by consuming some of it. But
since the plume has ceased to upwell due to temporary
choking of the plume conduit over the blob in the upper
asthenosphere region8, the water temperature falls instead
of increasing even after enlargement of the cracks. This
is compatible with a similar fall in the radon emission
rate level at the well (Figures 1 and 2 b).

The subsequent
RESEARCH COMMUNICATIONS
1290 CURRENT SCIENCE, VOL. 86, NO. 9, 10 MAY 2004
rise in radon emission rate level (Figure 2 a) and watertemperature
increase (Figure 2 b) are also compatible. It
is worthwhile to note that during observation of radon
concentration in the well of Figure 1 from May 1977 to
August 1986, no significant change was observed due to
other earthquakes of magnitude greater than M 5 within
an epicentral distance of 100 km and greater than M 6.5
within 600 km4.

A simple mathematical model of elliptic spreading of
hot magma (for clue see Figure 1 b of Manga11), in the
Moho, ejected through the plume conduit on top of a blob
is being proposed. According to this model, the twodimensional
magma flow in the Moho is turbulent at any
point of time. According to some opinion, turbulent flows
cannot be modelled by the Navier–Stokes equations12,13.
Here I propose the following set of equations for modelling
the two-dimensional magma spreading through the
Moho. (Equations (1) and (2) represent a general form of
flow in a plane. By choosing the values of ni, Ai, mi and
the form of the functions fi, any flow in a plane can be
modelled by eqs (1) and (2), as shown here in case of an
elliptical flow.)
( , , )
1( 1 2 ) 1( 1 2 ) 1 , , , , e 1 1 1 2 m f x x t u x x t v x x t A - = + , (1)
.
( , , )
2 ( 1 2 ) 2 ( 1 2 ) 2
u x , x , t = v x , x , t + A e-m2 f2 x1 x2 t (2)
ui is the total velocity component along xi, vi is the linear
velocity along xi, Ai emifi(xi,x2,t) is the curvilinear (spiral)
component of velocity along xi, for i = 1, 2. Ai and mi
are constant numbers whose values lie in some closed,
bounded intervals of the real line. (The absolute value of
mi is an increasing function of Rayleigh number and
determines the angular shape of curving, Ai determines
the radius of the curving), fi is a function giving the
angular span of curving. A1 and A2 are of the order of tens
of kilometres. Moreover, if the flow is approximately
elliptical, then we should take
log(sin )
1
1
f1 t m = and
log(cos )
2
1
f2 m t = ,
where t is of the order of years and therefore sint and cost
can be taken as positive in a few months to few years
range. Since here I am concerned with the curvilinear
spreading of flow only, keeping the centre of propagation
(the hole of plume ejection) fixed, v1 = v2 = 0.

Geochemical heterogeneity within the upper mantle is
not limited to the subduction zones alone. Due to existence
of diverse isotopic suits throughout the upper
mantle, the upper mantle has been modelled like a ‘plum
pudding’14. So blobs can also form below, deeper inside
the plate. Ejection of plumes from such blobs may affect
the intraplate stress build-up, expediting an intraplate
earthquake.

In a two-layer mantle convection model, the concentration
of radiogenic heat sources (U-238, U-235, Th-232
and K-40) in the lower mantle is only about 1.5% of the
average heat-source concentration in the entire mantle8.
In a more detailed model, the upper mantle itself is subdivided,
according to viscosity stratification, into two
zones separated by a high-viscosity layer at 220 km
depth10. The density variation is the highest in the whole
mantle above this layer10, possibly forming the continental
and the oceanic crusts at the surface, like accumulation
of cream over boiled milk after cooling. Mantle flow
through the high viscosity barriers (at least four)10 is
Figure 2. Precursory changes of the 1978 Izu-Oshima-kinkai earthquake.
a, Radon concentration change observed at the well of Figure 1.
b, Water-temperature variation observed at a 500 m well (epicentral
distance 30 km). c, Water-level variation at the well of (b). d, Strain
change measured by a borehole strainmeter at the tip of the peninsula
(distance 50 km). After Wakita et al.4.
a
b
c
d
RESEARCH COMMUNICATIONS
CURRENT SCIENCE, VOL. 86, NO. 9, 10 MAY 2004 1291
maintained through upwelling of plumes in a smaller
scale, compared to the ones in Hawaii, which gradually
depletes the lower parts of the mantle of radiogenic heat
producing materials and enriches the upper parts of the
mantle with their deposit. This process becomes the
richest in radiogenic substances through the uppermost
high-viscosity layer at 220 km depth. The spike-like fluctuation
in the precursory radon concentration behaviour
before several earthquakes of intensity M 5.5 or more, as
recorded in Guza station (located at the intersection of
three major fault systems)4 in China (Figure 3), is
directly related to this phenomenon. The smaller plume
conduits, where plume strength varies in time and possibly
deformed by other motions in the flow, often tend to
get choked-up. Due to radiogenic heating inside the blob,
the pressure keeps building up. When the pressure becomes
sufficiently high, the choke in the conduit is blown
up, releasing a sudden spurt of U-238-rich plume for
some-time before being choked-up again. This may affect
the strain build-up in the plate directly above by expediting
the plate motion, sometimes (but not always)
leading to an intraplate earthquake. Radon, produced by
U-238 decay, might be reaching the station through some
cracks across the plate and after the plume conduit chokes
up, the uranium concentration below the plate is diluted
fast through rapid mixing in the low-viscosity and lowdensity
region.
The simulation given in Figure 4 will help analyse the
precursory change in geomagnetic field intensity during
some moderate-to-large earthquakes in China and Mongolia
(Figure 5).
The propagation of hot magma as shown in Figure 4, is
in general turbulent, due to numerous collisions with the
subsolidus rock. But if the radioisotope-rich magma at
some portion gets an average angular velocity, as shown
in Figure 4, then the emitted a, b and g rays due to
radioactive decay ionize some of the neighbouring materials
and release free electrons. The circulation of free
electron-rich magma in a path indicated by arrows in
Figure 4, will produce a magnetic field. According to its
orientation, it may increase or decrease the geomagnetic
field intensity as experienced on the surface of the earth.

Magma upwelling below a plate near the boundary
may help plate subsistence due to negative buoyancy (upwelled
magma is hotter and therefore less dense than the
surrounding). This may expedite an interplate strike-slip
earthquake, like the one that occurred on 12 June 1897 in
Assam, India (M 8.1)15. Radon concentration fluctuation
monitoring has been going on for about two decades in
the Himalayan region also6,16. Here too, the correlation
between erratic fluctuations in radon concentration and
seismic activity is clear6,16. It has been demonstrated that
Figure 3. Spike-like changes in radon concentration observed at
Guza station, Sichuan Province, China before several earthquakes of
magnitude M 5.5 or more, within an epicentral distance of 320 km.
After Wakita et al.4.
Figure 4. Magma propagation below the plate after upwelling
through an axisymmetric plume at the centre (diameter 2 km). In
reality, the upwelled magma will collide with the subsolidus rock in the
Moho and the velocity vector will be random at each point. A more
complicated calculation of eqs (1) and (2) at a finer resolution will give
a more realistic simulation. But even the present simple simulation will
help explain phenomena like precursory change in geomagnetic field
intensity4. If the radioisotope-rich magma in certain regions gets an
average angular velocity as indicated by the arrows, it will affect the
geomagnetic field intensity.

RESEARCH COMMUNICATIONS
1292 CURRENT SCIENCE, VOL. 86, NO. 9, 10 MAY 2004
even in case of small earthquakes (M 4.5 or less), erratic
fluctuations in radon concentration is clearly observable
before, during and after an earthquake16. Radon concentration
fluctuation data collected from Chamba station (in
Himachal Pradesh, India) during March 1995 have been
presented in time series form (see figure 3, Virk et al.6).
An earthquake of M 5.1 occurred in Chamba on 24 March
1995. Here too, several days before the quake, radon
concentration level went down (possibly because of choking
as described above) and then suddenly shot up
abnormally (because of the sudden burst of the choked
conduit) a few days before the earthquake. Right at the
time of the quake, it reached almost the usual level. But
the earlier events had already expedited the plate motion,
leading to the 24 March earthquake, or at least some
advancement of it.
Validating an earthquake precursor through an acceptable
geodynamic modelling and accounting for its occurrence
is a challenging task. Earthquakes are seen as consequences
of plate tectonics. Plate motions are controlled
by the complex dynamics of mantle flows, which depends
on the thermal states and chemical composition of the
mantle. Upwelling of hot plumes and downwelling of
colder, denser materials keep running the cycles of vertical
mantle flow over geologic time. Therefore, it is only natural
to expect that plume upwelling and propagation below
the plates may have significant effects on the seismic
activities. Based on the findings of some recent research9,10,
here a model has been proposed for upwelling and spreading
of radioisotope-rich magma below the plates, which
leaves its signature in the erratic radon concentration
fluctuations, as observed before, during and after some
earthquakes in China, Japan and India. With this model,
the geomagnetic field intensity change as a precursor to
some earthquakes can also be explained. Based on the
validity of the model, radon concentration and geomagnetic
field intensity variation give a direct measure of the
velocity field and chemical composition of magma
upwelling and propagation below the plate, from which
the prospect of an intermediate-term prediction of an
impending large earthquake can sometimes be improved.
This should be viewed as a complementary and not
competitive effort for earthquake prediction in addition to
the existing ones like the measurement of rupture
propagation3,17. Further testing of this model with more
diverse datasets is needed.
1. The Cambridge Encyclopaedia of Earth Sciences, Cambridge
University Press, Cambridge, 1981.
2. Geller, R. J., Shake-up for earthquake prediction. Nature, 1991,
352, 275–276, Nature, 353, 611–612.
3. Vidale, J. E., Do big and little earthquakes start differently? Science,
1996, 271, pp. 953–954.
4. Wakita, H., Nakamura, Y. and Sano, Y., Short-term and intermediateterm
geochemical precursors. In Intermediate-Term Earthquake
Prediction (eds Stuart, W. D. and Aki, K.), Birkhauser, Basel,
1988, pp. 267–278.
5. Donald, T., Geochemical precursors to seismic activity. In Intermediate-
Term Earthquake Prediction (eds Stuart, W. D. and Aki,
K.), Birkhauser, Basel, 1988, pp. 241–266.
6. Virk, H. S., Walia, V. and Sharma, A. K., Radon precursory signals
of Chamba earthquake. Curr. Sci., 1995, 69, pp. 452–454.
7. http://www.azgs.az.gov/RADONDTE.htm, Spencer, J. E., Radon
gas – a geologic hazard in Arizona. In Arizona Geological Survey:
Down-to-Earth Series 2, 1992.
8. Schubert, G., Turcotte, D. L. and Olson, P., Mantle Convection in
the Earth and Planets, Cambridge University Press, Cambridge,
2001.
9. Davaille, A., Simultaneous generation of hotspots and super-wells
by convection in heterogeneous planetary mantle. Nature, 1999,
402, 756–760.
10. Forte, A. M. and Mitrovica, J. X., Deep-mantle high-viscosity
flow and thermochemical structure inferred from seismic and geodynamic
data. Nature, 2001, 410, 1049–1056.
11. Manga, M., Shaken, but not stirred. Nature, 2001, 410, 1041–
1043.
12. Ladyzhenskaya, O. A., The Mathematical Theory of Viscous Incompressible
Flow, Gordon and Breach, New York, 1969.
13. Hardy, J. and Pomeau, Y., Microscopic model for viscous flow in
two dimensions. Phys. Rev. A, 1977, 16, 720–726.
14. Batiza, R., Inverse relationship between Sr isotope diversity and
rate of oceanic volcanism has implications for mantle heterogeneity.
Nature, 1984, 309, 440–441.
15. Bilham, R. and England, P., Plateau ‘popup’ in the great 1897
Assam earthquake. Nature, 2001, 410, 806–809.
16. Virk, H. S., Sharma, A. K. and Walia, V., Correlation of alphalogger
radon data with microseismicity in N-W Himalaya. Curr.
Sci., 1997, 72, 656–663.
17. Kasahara, K., Earthquake Mechanics, Cambridge University Press,
Cambridge, 1981.
18. Hart, S. R., A large-scale isotope anomaly in the southern hemisphere
mantle. Nature, 1984, 309, 753–757.
19. van der Hilst, R. D. and Karason, H., Compositional heterogeneity
in the bottom 1000 kilometers of Earth’s mantle: Toward a hybrid
convection model. Science, 1999, 283, 1885–1888.
ACKNOWLEDGEMENTS. During this research I was supported by
a Senior Research Fellowship of the Indian Statistical Institute and
received some technical help from Mr Subhashis Mitra. Helpful suggestions
made by the anonymous reviewers are acknowledged.
Received 15 July 2003; revised accepted 24 February 2004
Figure 5. Geomagnetic field intensity change observed in Beijing,
during the earthquakes in Holinkoerh, Inner Mongolia (6 April, M 6.3),
Tacheng (22 April, M 4.4) and Tangshan (28 July, M 7.8 and 7.1), both
in China. Taken from Wakita et al.4 after editing out the portions which
are not relevant to this study.

E ora torniamo in Italia.... ma siamo proprio sicuri che non se ne sapeva nulla?

Tracking precursors and degassing by radon monitoring during major eruptions at Stromboli Volcano (Aeolian Islands, Italy)

C. Cigolini DSMP, Università di Torino, Turin, Italy

G. Gervino DFS, Università di Torino, Turin, Italy

R. Bonetti IGGA, Università di Milano, Milan, Italy

F. Conte DSMP, Università di Torino, Turin, Italy

M. Laiolo DSMP, Università di Torino, Turin, Italy

D. Coppola DSMP, Università di Torino, Turin, Italy

A. Manzoni IGGA, Università di Milano, Milan, Italy

Ma sono tutti ricercatori Italiani che scrivono in Inglese!!! Come mai?

We monitored radon anomalies before, during and after the recent paroxysmal eruptions of Stromboli. Major eruptive events are preceded by relative minima in radon emissions at stations located at the base of the cone, whereas three of the summit stations reached peak values (∼20,000 Bq/m3 or higher) 12 to 14 days before major eruptions, and are nearly coeval with earth tides. Radon emissions at/or above such high threshold values are proposed to result from the dynamic response of the fracture network due to degassing processes associated with magma ascent. Relative minima recorded at the base of the cone are related to fractures' self-sealing within the hydrothermal system, eventually coupled with the atmospheric stack-effect. The coexistence of the above anomalies is interpreted as the first precursory signals to major eruptive events.

Received 2 February 2005; accepted 19 May 2005; published 24 June 2005.

Citation: Cigolini, C., G. Gervino, R. Bonetti, F. Conte, M. Laiolo, D. Coppola, and A. Manzoni (2005), Tracking precursors and degassing by radon monitoring during major eruptions at Stromboli Volcano (Aeolian Islands, Italy), Geophys. Res. Lett., 32, L12308, doi:10.1029/2005GL022606.

Questi posso prevedere addirittura le urruzioni vulcaniche con 12 / 14 giorni di anticipo!!! Viva i Profeti del Vulcano !!!

The Smithsonian/NASA Astrophysics Data System

Natural Radioactive Characterization by Gamma Ray Mapping in Relation with Radon-222, Radium-226 Contents in the Atmosphere, Dripping Groundwater and Mineral in a Dead-end Horizontal Tunnel

American Geophysical Union, Fall Meeting 2004, abstract #H21E-1059
The concentration of radon-222 has been monitored since 1995 in the atmosphere of a 2 m transverse dimension, 128 m long, dead-end horizontal tunnel located in the French Alps, at an altitude of 1600 m. The total volume of the tunnel is about 496 m3 with an estimated total wall surface area of 1041 m2. Most of the time, the radon concentration is stable, with an average value ranging from 200 Bq m-3 near the entrance to about 1000 Bq m-3 in the most confined section, and with an equilibrium factor F between radon and its short-lived decay products varying from 0.61 to 0.78. However, radon bursts are repeatedly observed, with amplitudes reaching up to 36 x 103 Bq m-3 and durations varying from one to several weeks, with similar spatial variations along the tunnel as the background concentration. These spatial variations are qualitatively interpreted in terms of natural ventilation. To understand the source term, we realised a comprehensive 3D map of total gamma emission of rock on all surfaces of the tunnel, i.e., floor, wall and roof, representing 700 measurements points with a section of eight points by two-meter steps along the tunnel. Dose rate varies from 100 to 350 nSv h-1 and reveals two major contrasted geological structures, correlated with the radon-222 profile in the air along the tunnel (230 Bq m-3 at the entrance to 770 Bq m-3 at the end, measured during a period of background radon level), and with other profiles measured in previous work. In addition, a correlation is evidenced with the radon-222 contents in dripping groundwater and with the radon precursors - i.e., the daughters of uranium-238 (radium-226, for example) - content in rocks. The radon and radium contents in water and rocks, and the study of flow rate of the dripping groundwater are in direct relation with the mechanism of radon gas transport by water in porous media. By interpreting all available measurements, it is now possible to model the time variation of radon-222 in the tunnel, and to constrain the origin of the radon peaks generated in the gallery. The bursts result from transient radon exhalation upsurges at the tunnel wall, that could be due to meteorological effects, or possibly to combined hydrological and mechanical forcing associated with the water level variations of the nearby Roselend reservoir lake.

Mechanical Forcing suona un po come stress sismico non e' vero? Ma ora andiamo in Antartide....

The peculiarities of cross-correlation between two secondary precursors - radon and magnetic field variations, induced by tectonic activity  ( Ma daiii... )

Abstract: A model of precursor manifestation mechanisms, stimulated by tectonic activity and some peculiarities of observer strategy, whose main task is the effective measurement of precursors in the spatial area of their occurrence on the Earth's daylight, are considered. In particular, the applicability of Dobrovolsky's approximation is analyzed, when an unperturbed medium (characterized by the simple shear state) and the area of tectonic activity (local inhomogeneity caused by the change only of shear modulus) are linearly elastic, and perturbation, in particular, surface displacement is calculated as a difference of the solutions of two independent static problems of the theory of elasticity with the same boundary condition on the surface. Within the framework of this approximation a formula for the spatial distribution (of first component) of magnetic field variations caused by piezomagnetic effect in the case of perturbed regular medium, which is in simple shear state is derived. Cogent arguments in favor of linear dependence between the radon spatial distribution and conditional deformation are obtained. Changes in magnetic field strength and radon concentrations were measured along a tectonomagnetic profile of the total length of 11 km in the surroundings of the "Academician Vernadsky" Station on the Antarctic Peninsula (W 64 16', S 65 15'). Results showed a positive correlation between the annual surface radon concentration and annual changes of magnetic field relative to a base point, and also the good coincidence with theoretical calculation.
Cite as: arXiv:physics/0605244v1 [physics.geo-ph]

 

Allora vediamo se in casa nostra nonlo sapeva proprio nessuno... o e' il solito segreto di pulcinella...

A. Riggio, S. Sancin, G. F. Gentile, P. Zennaro e R. Belletti

Istituto Nazionale di Oceanografia e di Geofisica Sperimentale, Trieste

MISURE DI RADON IN FRIULI E CONFRONTO CON LA SISMICITÀ

 

Riassunto. Scopo del lavoro è quello di contribuire alla conoscenza dello stato di stress presente nelle strutture, verificando eventuali correlazioni tra emissione di radon e attività sismica. A tal fine, nel 1995, è stata installata una stazione per il rilevamento del gas radon in Friuli, in funzione continuativamente dal 1996 a meno di alcune piccole interruzioni dovute a cause tecniche. Nel presente lavoro è considerato il periodo 1996-1997 per uniformità di dati. Contemporaneamente sono stati controllati i parametri meteorologici e l'andamento della falda freatica. È stata riscontrata una scarsa dipendenza della concentrazione di radon dai parametri meteorologici, ma temperatura e radon hanno in comune le variazioni giornaliere e le variazioni stagionali. Tramite la "crosscorrelazione del segnale giorno per giorno", è stato messo in evidenza che le possibili "anomalie" sono costituite da variazioni nell'andamento giornaliero del gas radon. L'accostamento degli eventi sismici, verificatisi durante il periodo preso in esame, mostra i possibili legami tra variazione di radon e terremoti più forti. Le variazioni più significative sono state accompagnate da variazioni del livello di falda non giustificate dalla pluviometria.L'anomalia più evidente, anche nella semplice rappresentazione della concentrazione del gas radon, è quella presente durante il periodo gennaio-marzo 1996. In questo periodo si è verificata la serie sismica di Claut costituita da tre terremoti principali che, con 3,7, 4,0 e 4,3 di magnitudo, sono stati i più forti con epicentro entro il territorio regionale, da gennaio 1996 a dicembre 1997.

 RADON MEASUREMENTS IN FRIULI  AND COMPARISON WITH THE SEISMICITY

 Abstract. The purpose of this work is to contribute to knowledge of the stress status of the structures by the control of the possible correlation between radon emission and the seismic activity.For this end, on 1995, a radon survey station has been set up in Friuli and since 1996 it is recording continuously. In this work has been considered the period 1996-1997 for homogeneity of the radon data. At the same time have been controlled meteorological parameters and the groundwater level. The radon concentration results connected with the meteorological parameters in a quite weak way, but temperature and radon show the same behaviour and the same seasonal variation.Through the "day by day croscorrelation of the signal" has been pointed out that the "anomaly" consists mainly of a daily behaviour variation of the radon concentration. The coupling of the seismic events have showed the correlation between the radon variation and the strongest earthquakes. During more important radon variations, also the groundwater level changes without any rain justification.The clearest anomaly, in the same radon concentration representation too, has been noticed during the period January-March 1996. In this period happened the Claut seismic series constituted by three main earthquakes that with magnitude 3.7, 4.0 and 4.3, have been the strongest verified over whole region since 1996 January to 1997 December.

 INTRODUZIONE

 La Previsione dei terremoti, intendendo con questa, quando, dove e con quale energia si verificherà un evento sismico, come sanno bene tutti gli "addetti ai lavori", rimane ancora un'aspirazione per tutti (Wyss, 1991).

I tentativi portati avanti in passato sia dalle varie scuole (russa, cinese, greca, ecc.) che da dilettanti, hanno trovato spesso l'opposizione della comunità scientifica, perché ritenuti non significativi. Il nocciolo del problema è capire come gli stress che possono precedere un terremoto o provocare lenti movimenti, si propaghino nel sottosuolo e come le strutture reagiscano ad essi. Questo comporta l'analisi delle strutture dal punto di vista dinamico attraverso lo studio delle variazioni nel tempo dei parametri fisici sensibili agli stress (Ren et al., 1984; Kissin et al., 1996; Cosentino et al., 1976; Riggio e Sancin, 1986; Riggio  et al., 1987). Tra questi parametri, la concentrazione di radon ha la caratteristica di poter essere rilevata con un costo relativamente basso e di essere molto sensibile alle variazioni di stress anche profondo e proveniente da strutture più lontane. Il radon si diffonde attraverso le rocce e si disperde nel suolo e nell'acqua in cui ha buona solubilità (Dongarrà e Martinelli, 1995; Martinelli, 1992; Martinelli, 1993).

La misura del radon è utilizzata per scopi tettonici (per individuare faglie sepolte, controllare quelle attive dopo disastrosi eventi sismici, etc.) (Facchini et al., 1995), per individuare dei percorsi idrici e per scopi ambientali (radon indoor, individuazione delle perdite di biogas nelle discariche, inquinamenti da idrocarburi, ecc.).

Le fratture che si formano prima di un terremoto, secondo il modello della dilatanza (Birchard e Libby, 1980) incrementano la diffusione del fluido nei pori, quindi modificano la resistenza e la pressione della roccia cosicché il flusso del radon viene in qualche modo modulato dagli stress. Queste variazioni sono considerate significativi precursori del terremoto (Riggio e Sancin, 1995).

Le prime applicazioni di questo tipo di misure furono fatte in Unione Sovietica, e i primi risultati positivi si ebbero in coincidenza del terremoto del 1966 a Tashkent (Mogi, 1985).  Il graduale incremento nei valori di concentrazione di radon iniziò già sei anni prima del terremoto e si stabilizzò alcuni mesi prima. Le anomalie di radon si presentano in forme differenti: con incrementi o decrementi, con variazioni lente o rapide tipo "spike" (Friedman et al., 1988; Heinicke et al., 1995). Attualmente misure continue di radon vengono eseguite in Cina, Giappone, Stati Uniti e altri paesi. Negli Stati Uniti ebbero occasione di notare anomalie di radon in corrispondenza di quattro sismi con magnitudo tra 4,2 e 4,7.

Anche in Italia dall'inizio degli anni '90 si sono sviluppati gli studi di monitoraggio della geochimica dei fluidi miranti ad ottenere le variazioni dei vari parametri, sia nello spazio che nel tempo (Calcara e Quattrocchi, 1993a, 1993b; Carapezza et al., 1980; Dall'Aglio, 1992).

 IL PROCESSO FISICO

 La misura della concentrazione di radon è eseguita più frequentemente in acqua, (Friedmann e Hernegger, 1978; Heinecke et al., 1995), ma non mancano rilevamenti in suolo (Friedman et al., 1988; Sabbarese et al., 1993). Nonostante il radon abbia un tempo di dimezzamento di soli 3,5 giorni (Montagner, 1994), sfruttando le correnti ascensionali di altri gas riesce per convezione a raggiungere la superficie anche da profondità di vari chilometri. Le variazioni di radon possono, quindi, essere causate da terremoti con epicentro distante anche diversi chilometri ma con energia sufficiente. La formula empirica (Hauksson and Goddard, 1981):

                                                       M = 2,4 log10 (D) - 0,43                                                   (1)

 definisce la magnitudo minima M richiesta per ottenere un'anomalia di radon ad una distanza D espressa in chilometri. Indicativamente un terremoto di magnitudo 4,0 può essere preceduto da un'anomalia di radon fino ad una distanza di 70 km. Questo valore è interessante in quanto dell'ordine delle massime distanze dei confini regionali da Cazzaso, sito di rilevamento del radon.

 LINEAMENTI  SISMOTETTONICI

 Nel 1995 è stata realizzata una stazione per il rilevamento del gas radon in Friuli, in località Cazzaso nel comune di Tolmezzo. Nei primi mesi il suo funzionamento è stato saltuario, e dal 1996 è in funzione quasi continuativamente. 

La zona fa parte dell'area Carnico-Friulana caratterizzata da accavallamenti tagliati da faglie subverticali a direzione N-S e da altre a direzione dinarica, tutte con attività spesso trascorrente. L'insieme costituisce, unitamente all'antistante pianura veneto-friulana, la Catena sudalpina orientale (Bernardis et al., 1996; Slejko et al., 1987).

L'ubicazione della stazione per il rilevamento del gas radon è riportata in Fig. 1 dove è riprodotto anche lo schema tettonico. La zona si trova immediatamente a nord di quella con la più alta attività sismica che comprende anche l'evento del 1976 di magnitudo 6,5 (Riggio, Gentile, 1984). In accordo allo schema riportato, l'abitato di Cazzaso si trova in prossimità della linea Alto Tagliamento e di una faglia subverticale con direzione Dinarica con attività spesso trascorrente.

Il paese è, oltretutto, interessato da un evidente dissesto franoso che interessa la zona delle pendici del monte Diverdalce, la cui nicchia è ben visibile anche da lontano. La frana sembra avere avuto origine nel novembre 1851 a seguito di un periodo d'intense piogge. Nel 1898 il movimento franoso causò la distruzione di alcuni edifici e l'intero paese sarebbe scivolato di circa trenta metri verso valle.

Storicamente in questa zona avvennero terremoti disastrosi come quelli di Tolmezzo del 20 ottobre 1788 (Gentile et al., 1984) e del 27 marzo 1928 dell'VIII-IX grado della scala MCS (Mercalli-Cancani-Sieberg). Poco più ad ovest, il terremoto di Raveo nel 1700 ebbe un'intensità del X grado MCS ma i suoi effetti furono rilevanti in un'area molto limitata. Più recentemente, i terremoti che sono avvenuti il 3 febbraio 1949 a Paularo e il 26 aprile 1959 ad Arta, rispettivamente a circa 17 km e a 7 km di distanza da Cazzaso, hanno avuto un'intensità epicentrale dell'VIII grado MCS.

La sismicità attuale della zona è costituita da terremoti con magnitudo intorno al 3,0 ed epicentri tra Tolmezzo, Cazzaso e Villa Santina.

 DESCRIZIONE DELLA STRUMENTAZIONE

 La stazione di rilevamento del gas radon è costituita da un rilevatore denominato Prassi mod.55 con cella a scintillazione, predisposto sia per la misura continua che campionata. Il rilevatore è costituito da una camera rivestita internamente di solfuro di zinco attivato argento. Il volume di questa cella è di 2 litri ed è accoppiato ad un fotomoltiplicatore a bassa deriva. Lo strumento misura i conteggi in un tempo prefissato in modo ripetitivo. La concentrazione di gas radon in Bq/m3 o pCi/l è calcolata dai conteggi misurati. La larga capacità di memoria permette di memorizzare fino a 6500 misure. Lo strumento è stato installato all'imboccatura di un pozzo profondo 40,50 metri con un livello dell'acqua, dal piano di campagna, di circa 17,5 metri (Fig. 2). La stratigrafia relativa ai primi metri dal piano di campagna è un'alternarsi di limi sabbiosi-argillosi e ghiaie sino alla profondità di circa 17 metri, dove iniziano argille più dure, seguono sabbie, poi siltiti marnose con venature di calcite.

 
 

Fig.  SEQ Fig. \* ARABIC 1 - Ubicazione della stazione di rilevamento del gas radon e schema tettonico del Friuli Occidentale da Bernardis et al., (1996).

 

L'aria viene aspirata, tramite una sonda flessibile, dal pozzo alla profondità di circa 7 metri dal piano di campagna e a circa 6 metri dalla superficie dell'acqua. L'intervallo di campionamento è stato fissato in tre ore e 30 minuti, dall'inizio delle registrazioni sino al 27 maggio 1996, in quattro ore, dal 27 maggio al 17 giugno 1997 e in tre ore, dal 17 giugno in poi. Questi cambiamenti si sono resi necessari per poter correlare la concentrazione del radon con le misure meteorologiche che hanno campionamento orario. Lo strumento è stato protetto, dagli effetti meteorologici, inserendolo in un contenitore in plexiglas con coperchio a chiusura stagna e alloggiando il tutto in un armadietto in alluminio con barre di rinforzo e chiusura con lucchetto. L'alimentazione è realizzata tramite batterie ricaricate con due pannelli solari e regolatore di carica.

Attualmente, circa ogni mese viene fatto il recupero dei dati in loco tramite un PC portatile, ma sarebbe possibile anche trasmetterli direttamente nella sede dell'elaborazione.

La Regione Friuli-Venezia Giulia, negli anni 1995-1997 ha stipulato una convenzione con l'OGS al fine di valutare l'evoluzione del fenomeno franoso. La zona era, quindi, contemporaneamente attrezzata con pozzi dotati di piezometri e inclinometri tarati per i movimenti superficiali delle frane, nonché di termometro per la temperatura dell'acqua e dell'aria, di pluviometro e di barometro.

 

Fig.  SEQ Fig. \* ARABIC 2 - Stazione di rilevamento del gas radon. 

ANALISI DEI DATI

Scopo del lavoro è quello di contribuire alla conoscenza dello stato di stress presente nella struttura. A tal fine si è proceduto a verificare la presenza di una quantità sufficiente di radon, la sensibilità della concentrazione di radon alle variazioni di stress e infine se, in qualche modo, la presenza dello stress viene confermata da terremoti.

Dal momento che durante il 1998 è stata apportata una modifica al sistema di rilevamento, per uniformità di dati è stato considerato il periodo 1996-1997.

La conoscenza dei parametri meteorologici ha permesso una migliore valutazione dei dati radon, e di constatare, che con il sistema d'acquisizione usato, le influenze ambientali sono modeste. Solamente in caso di piogge intense, che possono penetrare all'interno del pozzo, si possono avere variazioni indotte anche notevoli (in due anni si è verificato un solo caso, con 40 mm di pioggia in un giorno).

I dati sismometrici sono quelli registrati dalla Rete Sismometrica del Friuli Venezia-Giulia, gestita dal CRS (Centro Ricerche Sismologiche) dell'OGS di Trieste, dalla stazione Broad Band di Trieste, gestita dall'OGS e dall'Università di Trieste e dal catalogo USGS. I terremoti presi in esame sono tutti quelli con magnitudo maggiore o uguale a 2,0 (valore minimo per poter considerare completo il catalogo), e che soddisfano l'equazione (1). Si tratta in tutto di 95 terremoti, di cui 88 localizzati entro il territorio regionale (Fig. 3).

Fig.  SEQ Fig. \* ARABIC 3 - Mappa degli epicentri dei terremoti locali utilizzati.

 In Fig. 4 è riportato il diagramma relativo al radon, alla temperatura dell'aria, alla pluviometria, al livello della falda nel pozzo e ai terremoti considerati per il 1996 e in Fig. 5 le stesse grandezze relative al 1997.

È evidente la variazione stagionale della concentrazione del radon. In inverno il suo valore si aggira intorno ai 500 Bq/m3, mentre nel periodo estivo il suo valore medio è di circa 5000 Bq/m3 con picchi che raggiungono anche i 24000 Bq/m3. Questo andamento costituisce la componente a bassa frequenza del segnale per noi di poco interesse. Per tale componente a lungo periodo l'analisi spettrale del segnale relativo al radon non è significativa, dal momento che il campione di tempo è troppo piccolo in relazione alla periodicità da evidenziare. Il segnale però è costituito anche dalla componente ad alta frequenza. Per meglio evidenziare quest'ultima, in Fig. 6 è rappresentato l'andamento della concentrazione di radon per un periodo di dieci giorni, a partire dal 17/6/97, con il corrispondente andamento della temperatura. Ambedue i segnali hanno delle periodicità giornaliere. Questo è ben visibile nell'analisi spettrale dello stesso periodo. Sono contemporaneamente presenti, in ambedue i segnali, una forte componente alle 24 ore e una molto debole alle 12 ore, mentre per le altre componenti i due segnali sono indipendenti e senza una relazione causale. Questo sarà confermato dalle successive elaborazioni. 

 

Fig.  SEQ Fig. \* ARABIC 4 - Concentrazione del radon a), temperatura dell'aria b), pluviometria c), quota idrometrica della falda d) e terremoti locali considerati durante il 1996 e).

 

Fig.  SEQ Fig. \* ARABIC 5 - Concentrazione del gas radon a), temperatura dell'aria b), pluviometria c), quota idrometrica della falda d) e terremoti locali considerati durante il 1997 e).

 CONCENTRAZIONE DI RADON E TERREMOTI 

Dall'analisi della concentrazione di radon emerge un evidente e inatteso andamento dei valori, durante il periodo gennaio-marzo 1996. In tale periodo si è verificata la serie sismica di Claut che si trova a circa 35 km dalla stazione di rilevamento del radon. Le tre scosse principali di magnitudo 3,7, 4,0, e 4,3 verificatesi rispettivamente il 27 gennaio, il 27 febbraio e il 13 aprile hanno avuto un meccanismo del tipo faglia inversa come tutti i loro foreshocks, con un'unica eccezione per quello dell'ultimo terremoto più forte (Bernardis et al., 1996). In corrispondenza di quest'ultimo, la concentrazione di radon è ormai sui valori medi, indicando che la fase di massima compressione nelle strutture sottostanti la stazione di rilevamento, si è esaurita. L'anomalia, in bassa frequenza (Fig. 7), della durata di una decina di giorni, dal 20 febbraio al 2 marzo, raggiunge un valore massimo di 7 kBq/m3 superando di ben tre volte il valore medio dei valori nei periodi circostanti. L'andamento anomalo inizia in realtà il 24 gennaio con un rapido incremento dei valori e il loro mantenimento anche durante le ore notturne. Questo fenomeno costituisce la principale caratteristica delle variazioni in alta frequenza che analizzeremo in seguito.  In modo meno evidente la perturbazione continua fino a maggio. Tra il 22 e il 23 febbraio si verifica anche un repentino decremento del livello della falda freatica di circa 30 cm e successivamente, tra il 6 e il 7 marzo, il livello risale fino a raggiungere i valori di partenza, delineando una forma anomala in controfase e in leggero ritardo rispetto a quella del radon. Il tutto senza la presenza di piogge.

Dal momento che le variazioni osservate sembrano consistere in variazioni dell'andamento giornaliero, si è proceduto all'introduzione della correlazione, senza ritardo, dei valori di ogni giorno con quelli del giorno precedente. Nelle Figg. 8 e 9 sono riportate, rispettivamente per gli anni 1996 e 1997, tali correlazioni relative alla concentrazione di radon e della temperatura dell'aria insieme al grafico delle magnitudo dei terremoti presi in considerazione. Risulta evidente che, in coincidenza o prima dei terremoti, ci sono coefficienti di crosscorrelazione con valori tra 0,95 circa e 0,55, ma è difficile trovare una diretta relazione tra il valore basso di crosscorrelazione e il terremoto poiché ci sono molti terremoti di piccola magnitudo ravvicinati nel tempo. Solamente per i valori di coefficienti più bassi è possibile vedere una correlazione con i terremoti con magnitudo maggiore o uguale a 3,2. Da questa analisi si ricava il valore del coefficiente di crosscorrelazione per il quale è possibile parlare di “anomalia”.

Dai dati sperimentali sono state considerate come "anomalie" quelle legate a valori del coefficiente di crosscorrelazione minori o uguali a 0,70. Sulla base di questa definizione, i terremoti legati ai valori anomali della concentrazione di radon durante il 1996 sono: la serie sismica di Claut (Friuli Occidentale), e il terremoto del 15 ottobre con magnitudo 5,8 e distante 250 km dalla stazione radon. Durante il 1997 i terremoti sono: quello del 17 marzo, distante 25 km e con magnitudo 3,2, la serie sismica dell'Umbria, iniziata il 26 settembre con il terremoto di magnitudo 6,4 e, distante 380 km e i terremoti della Slovenia del 13 marzo e del 12 aprile 1998 con magnitudo rispettivamente 5,2 e 6,0 e distanti da Cazzaso 132 km e 25 km. Tutti quanti sono riportati nelle Figg. 8b e 9b.    

Da queste ultime rappresentazioni è evidente che le anomalie di radon non sono causate da corrispondenti variazioni della temperatura. Anche le piogge non sembrano essere una causa diretta, poiché, non sono presenti in tutti i periodi in cui si sono riscontrate le anomalie di radon. Nelle Figg. 10, 11 e 12 sono stati riportati i particolari dei periodi sopra menzionati. Da notare che a meno dell'evento locale (M=3,2), gli eventi citati sono preceduti anche da rapide cadute del livello piezometrico.

 Fig.  SEQ Fig. \* ARABIC 6 - Concentrazione del gas radon a), temperatura dell'aria b), spettro di potenza del gas radon c) e spettro di potenza della temperatura d).


Fig.  SEQ Fig. \* ARABIC 7 - Concentrazione del gas radon a) e quota idrometrica della falda b) durante il periodo 11 gennaio - 20 marzo 1996. Le tacche indicano i tempi dei terremoti principali e ne riportano la magnitudo.

 Fig.  SEQ Fig. \* ARABIC 8 - Coefficienti, senza ritardo, della crosscorrelazione dei valori di un giorno con quelli del giorno precedente relativi alla concentrazione del radon a) e della temperatura dell'ambiente esterno c). Magnitudo dei terremoti presi in considerazione b). Anno 1996.

 

Fig.  SEQ Fig. \* ARABIC 9 - Coefficienti, senza ritardo, della crosscorrelazione dei valori di un giorno con quelli del giorno precedente relativi alla concentrazione del radon
 a) e della temperatura dell'ambiente esterno c). Magnitudo dei terremoti presi in considerazione b). Anno 1997.

 Pur non essendoci una dipendenza lineare, a valori di magnitudo più alta corrispondono i valori di coefficienti di correlazione più bassi, inoltre a parità di magnitudo
          i tempi precursori diminuiscono con l'aumentare della distanza, mentre a parità di distanza aumentano con la magnitudo.

 Nonostante ci sia una sola stazione, si è costruito il diagramma magnitudo-distanza (Fig. 13) che evidenzia la buona sensibilità del radon allo stress.

La linea continua tangente il dominio degli eventi capaci di produrre effetti sulla emissione di radon, sta ad indicare che le magnitudo richieste per eventi vicini sono più alte       rispetto alla relazione di Haksson e Goddard (1981), ma restano ancora accettabili.

L'equazione (1) per i nostri dati locali sperimentali diventa:

 

                                               Mmin   =  2,013 log10 (D) + 0,565                                          (2).

 

Fig.  SEQ Fig. \* ARABIC 10 - Concentrazione del gas radon a), quota idrometrica della falda b). Le tacche indicano i tempi dei terremoti principali e ne riportano le magnitudo. Anno 1996.


 

Fig.  SEQ Fig. \* ARABIC 11 - Concentrazione del gas radon a), quota idrometrica della falda b). Le tacche indicano i tempi dei terremoti principali e ne riportano le magnitudo. Anno 1997.

 

Fig.  SEQ Fig. \* ARABIC 12 - Concentrazione del gas radon a), quota idrometrica della falda b). Le tacche indicano i tempi dei terremoti principali e ne riportano le magnitudo. Anno 1997.

Fig.  SEQ Fig. \* ARABIC 13 - Diagramma magnitudo-distanza degli eventi principali (non dipendenti). * eventi associabili ad anomalie della concentrazione di radon, + eventi non accompagnati da anomalie, x eventi troppo vicini nel tempo ad altri di magnitudo maggiore, relazione empirica locale, relazione empirica Haksson e Goddard (1981).

 CONCLUSIONI

 Pur nei limiti degli eventi sismici verificatisi in regione, l'ubicazione molto favorevole della stazione di rilevamento ha permesso di evidenziare la validità della misura di radon per la segnalazione delle situazioni di stress presismico delle strutture tettoniche. Nell'unica situazione locale di maggior sviluppo di energia sismica, costituita dalla serie di Claut, si è avuta un'indicazione molto chiara di alterazioni dell'emissione di radon.

Ultimamente è stata modificata la tecnica di aspirazione dell'aria dal pozzo per cui la misura della concentrazione di radon sembra essere molto più indipendente dalle condizioni ambientali e climatiche. Con tale metodo di misura dovrebbe essere possibile individuare la relazione tra piccole anomalie e terremoti di bassa magnitudo.

La presenza di più punti di rilevamento permetterebbe di valutare l'influenza della distanza e della direzione della sorgente di stress sull'emissione di radon e correlare i risultati fra loro e con quelli delle altre misure passive per delineare i contorni delle strutture sismogenetiche.

 

 

Ringraziamenti. I dati sismologici, registrati dalla Rete Sismometrica del Fiuli, sono di proprietà del Dipartimento della Protezione Civile della Regione Friuli-Venezia Giulia.
I dati meteorologici e piezometrici sono di proprietà del Dipartimento Ambiente della Regione Friuli-Venezia Giulia. Si ringraziano il Sig. Medeot per lo sviluppo della strumentazione, il Sig. Jungwirth per la preparazione dei dati meteorologici, il Sig. Perini per i lavori in campagna e lo staff del Dipartimento CRS di Udine dell'OGS di Trieste per la preparazione dei dati sismologici.

Quindi la protezione civile della Regione Friuli Venezia Giulia lo sapeva già dal 1995!!

BIBLIOGRAFIA

 Bernardis G., Poli M.E., Renner G., Snidarcig A. e Zanferrari A.; 1996: Le tre sequenze sismiche del 1996 a Claut (Prealpi Carniche). Atti del 15° Convegno del Gruppo Nazionale di Geofisica della Terra Solida. Esagrafica srl, Roma. 343-348.

Birchard G.F. e Libby W.F.; 1980: Soil Radon concentration changes preceding and following four magnitude 4,2-4,7 earthquakes on the San Jacinto fault in Southern California. J.Geophys res. Vol. 85, n. B6, 3100-3106.

Calcara M. e Quattrocchi F.; 1993 a: Una nuova stazione di monitoraggio in continuo del 222Rn nell'aria dei Colli Albani: primi risultati. Atti del 12° Convegno del Gruppo Nazionale di Geofisica della Terra Solida. Esagrafica srl, Roma. 449-462.

Calcara M. e Quattrocchi F.; 1993 b:  Sulla scelta dei siti idonei al monitoraggio geochimico ai fini della sorveglianza sismica della Calabria settentrionale: Valle del Crati, Piana di Sibari. Atti del 12° Convegno del Gruppo Nazionale di Geofisica della Terra Solida. Esagrafica, Roma. 463-479.

Carapezza M., Nuccio P.M. e Valenza M.; 1980: Geochemical Precursors of Earthquakes. High pressure science and technology, Pergamon Press. 90-103.

Cosentino P., Luzio D. e Riggio A.M.; 1976: Geomagnetic surveys for siting a magnetic observatory in western Sicily. Bollettino di Geofisica Teorica e Applicata Vol. XIX n. 71, 143-150.

Dongarrà G. e Martinelli G.; 1995: Migration processes of  Radon towards the earth surface: Implications for the prediction of seismic and volcanic events In: Valerio Spagna & Enrico Schiavon (ed), Proceedings of Scientific Meeting on the seismic protection, Regione del Veneto, Dipartimento per la Geologia e le Attività estrattive. 114-118.

Dall’Aglio M. ; 1992:  Geochemical forerunners of earthquakes and volcanic eruptions. In:. Enzo Boschi and Michele Dragoni (ed), Earthquake Prediction. Il Cigno Galileo Galilei Edizioni di Arte e Scienza, Roma. 537-561.

Facchini U., Magnoni S., Ferlin A., Garavaglia M. e Verdiani G.; 1995: Emission de radon et d'hélium dans une faille dans les Préalpes de Bergamo. Gas geochemistry. Supplement to Volume 16 (1994) of Environmental Geochemistry and Health. Edit by Claude Dubois. Science Reviews. 281-294.

Friedmann H. e Hernegger F.; 1978: A method for continuous measurement of radon in water of spring for earthquake prediction. Geoph. Res. Lett., Vol. 5, n. 7, 565-568.

Friedmann H., Aric K., Gutdeutsch R., King C.-Y., Altay C. e Sav H.; 1988: Radon measurements for earthquake prediction along the North  Anatolian Fault Zone: a progress report.  Tectonophysics, Vol.152, n. 3\4, 209-214. 

Gentile G.F., Renner G., Riggio A.M., Slejko D. e Zacchigna M.; 1984: Problematiche riguardanti lo studio dei terremoti storici: tre esempi significativi. Atti del convegno Finalità ed esperienze della Rete Sismometrica del Friuli-Venezia Giulia. Ufficio Stampa e pubbliche relazioni della Regione Friuli-Venezia Giulia. 211-229.

Hauksson E. e Goddard J.G.; 1981: Radon Earthquake Precursor studies in Iceland  J.Geophys.Res. Vol. 86, 7037-7054.

Heinecke J., Kock U. e Martinelli G.; 1995: CO2 and Radon measurements in the Vogtland area (Germany) - a contribution to earthquake prediction research. Geoph. Res.Lett., Vol. 22, n.7, 771-774.

Kissin I.G., Belikov V.M. e Ishankuliev G.A. ; 1996: Short -term graundwater level variations in a seismic region as an indicator of the geodynamic regime. Tectonophysics, Vol.265, n. 3-4, 313-326.

Martinelli G.; 1992: Geochemistry and fluid dynamics of  subterrestrial precursory gases. In: Earthquake prediction. International School of Solid Earth Geophysics "Ettore Majorana Centre for  Scientific Culture", Erice. 563-592.

Martinelli G.; 1993: Radon geochemistry and geophysics in deep fluids in Italy. In Procedings of the Second Workshop on radon monitoring in radioprotection, enviromental and/or Earth Sciences. International Centre for theoretical Physics, Trieste, Nov, 25-Dec. 6, 1991. World Scientific, Singapore.

Mogi K; 1985:  Earthquake Prediction. Academic Press. 355.

Montagner D.; 1994: Anomalie di concentrazione di Rn 222 in acque sorgive delle dolomiti bellunesi. Tesi di Laurea. Università di Padova. Relatore Prof. C.Manduchi. 148.

Ren X., Qi G. e Zhan Z.; 1984: The change of geomagnetic total field in Beijing area before and after the 1976 Tangshan earthquakes. Acta Seismologica Sinica, Vol.6, n.3, 271-286.

Riggio A.M. e Gentile G.F.; 1984: Sismicità in Friuli dal maggio 1977 al dicembre 1982: distribuzione spazio-temporale e rilascio di energia. Atti del convegno: Finalità ed esperienze della Rete Sismometrica del Friuli-Venezia Giulia. Ufficio stampa e pubbliche relazioni della Regione Autonoma Friuli-Venezia Giulia. 93-157.

Riggio A.M. e Sancin S.; 1986. Variazione nel tempo del parametro b quale precursore. Atti del 5° Convegno Annuale del Gruppo Nazionale di Geofisica della Terra Solida. Esagrafica srl, Roma. 407-419.

Riggio A.M., Sancin S. e Slejko D.; 1987: Variazioni del rapporto delle velocità delle onde P ed S in relazione al verificarsi di eventi sismici. Atti del 6° Convegno Annuale del Gruppo Nazionale di Geofisica della Terra Solida. Esagrafica srl, Roma. 3-16.

Riggio A. e Sancin S.; 1995: Analysis of some physical and seismological parameters as earthquakes precursors. In: Valerio Spagna & Enrico Schiavon (ed), Proceedings of Scientific Meeting on the seismic protection, Regione del Veneto, Dipartimento per la Geologia e le Attività estrattive. 114-118.

Sabbarese C., De Martino S., Del Pezzo E., Maresca R., Cantore L. e De Cillis A.; 1993: Misure di diffusione di Radon nei suoli della valle dell’Irno.- Atti del 12° Convegno del Gruppo del Gruppo Nazionale della Geofisica della Terra Solida , Roma. 437-448.

Slejko D., Carulli G.B., Carraro F., Castaldini D., Cavallin A., Doglioni C., Iliceto V., Nicolich R., Rebez A., Semenza E., Zanferrari A. e Zanolla C.; 1987: Modello Sismotettonico dell’Italia Nord-Orientale. CNR - GNDT  Rendiconto n.1. Ricci, Trieste. 82

Wyss M.; 1991: Evaluation of proposed earthquake precursors. Geophys. Monogr. American Geophys. Union, Washington, M.Wyss editor. 94.

Ed ora possiamo infine associare il fenomeno con un secondo tipo di precursori sismici di ordine Elettro-Magnetico

The peculiarities of cross-correlation between two secondary precursors - radon and magnetic field variations, induced by tectonic activity
Rusov, V. D.
Maksymchuk, V. Yu.
Ilic, R.
Pavlovych, V. M.
Bakhmutov, V. G.
Saranuk, D. N.
Vaschenko, V. M.
Skvarc, J.
Hanzic, L.
Kosenko, S. I.
Location: http://arxiv.org/abs/physics/0605244

 
A model of precursor manifestation mechanisms, stimulated by tectonic activity and some peculiarities of observer strategy, whose main task is the effective measurement of precursors in the spatial area of their occurrence on the Earth's daylight, are considered. In particular, the applicability of Dobrovolsky's approximation is analyzed, when an unperturbed medium (characterized by the simple shear state) and the area of tectonic activity (local inhomogeneity caused by the change only of shear modulus) are linearly elastic, and perturbation, in particular, surface displacement is calculated as a difference of the solutions of two independent static problems of the theory of elasticity with the same boundary condition on the surface. Within the framework of this approximation a formula for the spatial distribution (of first component) of magnetic field variations caused by piezomagnetic effect in the case of perturbed regular medium, which is in simple shear state is derived. Cogent arguments in favor of linear dependence between the radon spatial distribution and conditional deformation are obtained. Changes in magnetic field strength and radon concentrations were measured along a tectonomagnetic profile of the total length of 11 km in the surroundings of the "Academician Vernadsky" Station on the Antarctic Peninsula (W 64 16', S 65 15'). Results showed a positive correlation between the annual surface radon concentration and annual changes of magnetic field relative to a base point, and also the good coincidence with theoretical calculation. 

Vediamo in Taiwan come fanno:

RADON MONITORING FOR EARTHQUAKE PREDICTION IN SOUTH CENTRAL TAIWAN CHYI, L. L.1, CHOU, C. Y.2, YANG, F. T.2, and CHEN, C. H.2, (1) Dept. of Geology, The Univ. of Akron, Akron, OH 44325-4101, lchyi@uakron.edu, (2) Dept. of Geosciences, National Taiwan Univ, 245 Chousan Rd, Taipei, 106, Taiwan

An integrated radon monitoring method is employed in Chuko fault zone in south central Taiwan for earthquake prediction.

Soil gas radon is monitored continuously with a solid-state detector and recorded with a data logger. The detector assembly is housed in a PVC pipe to reduce the influence of environmental factors. The fault zone is known to have deep source gases and sensitive to earthquake activities. It has spring upwelling and gases released were measured to have high 3He/4He ratios. The quantities of spring upwelling and natural gas releases are known to vary with earthquake activities.

By placing the detector assembly housing in a structured soil gas-collecting ditch of much larger size, radon level is raised but fluctuation reduced. Data retrieval from a remote laboratory site at the National Taiwan University, from the end of October 2000 to the end of February 2001, showed that spike-like radon anomaly precursors occurred essentially before every major earthquake with a ML > 4.0. The strong correlation between spike-like anomalies and major earthquakes suggests that these precursors are indeed seismic.

The precursors of the 15 major earthquakes during the recording period with a few service interruptions indicated that there are 0.49 to 7.40 days leading periods before the onset of earthquakes occurred within a 30 by 100-km radii ellipse. The trending of the long axis of the ellipse appears to reflect the direction of the current stress and the old structure trend of Taiwan. The new integrated approach in soil gas radon monitoring reveals that the spike-like anomalies are superior precursors and could be used to delineate time and the place of an earthquake. However, nature of basement geology could be another critical factor in determining these. Additional units of this type are being installed nearby with different basement geology for the testing.

Assai precisi anche nelle profezie.... e specialmente nella possibilita' che avvenga in un tempo ben determinato..... da 10 ore a 7 giorni e mezzo!!!

E questo libro allora!!!... Anomalia associata a micro terremoto...

A clear case of radon anomaly associated with a micro-earthquake event in a Stable Continental Region

Auteur(s) / Author(s)

REDDY D. V. (1) ; SUKHIJA B. S. (1) ; NAGABHUSHANAM P. (1) ; KUMAR Devender (1) ;

Affiliation(s) du ou des auteurs / Author(s) Affiliation(s)

(1) National Geophysical Research Institute, Hyderabad, INDE
Résumé / Abstract

Though many worldwide attempts to establish radon anomalies in soil and groundwater as precursors for major earthquakes are made, yet practically no such work has been reported for micro-earthquakes in a Stable Continental Region (SCR). From the continuous radon measurements for five months at a depth of 1 m in fractured granitic area, we report for the first time that enhanced radon concentration preceded seventeen hours the occurrence of micro earthquake. The peak radon level (twice the normal concentration) was reached thirteen hours before the occurrence of the event and decayed exponentially during the next 5 days to come to normal level. This radon precursor pertains to a dry period, unaffected by meteorological factors, thus exemplifies clearly its association with micro earthquake.

Risposta esatta!!!

Revue / Journal Title Geophysical research letters   ISSN 0094-8276   CODEN GPRLAJ 
Source / Source 2004, vol. 31, no10, pp. L10609.1-L10609.4 (20 ref.)
Langue / Language Anglais
Editeur / Publisher American Geophysical Union, Washington, DC, ETATS-UNIS  (1974) (Revue)
Mots-clés anglais / English Keywords

Dry spell ; concentration ; depth ; Precursor ; precursors ; ground water ; soils ; earthquakes ; anomalies ; radon ;

Mots-clés français / French Keywords

Période sèche ; Concentration ; Profondeur ; Précurseur ; Phénomène précurseur ; Eau souterraine ; Sol ; Séisme ; Anomalie ; Radon ;

Mots-clés espagnols / Spanish Keywords

Periodo seco ; Concentración ; Profundidad ; Precursor ; Fenómeno precursor ; Agua subterránea ; Suelo ; Sismo ; Anomalía ; Radón ;

Ma vediamo cosa e' il RADON...

Radon is an odorless, colorless, radioactive, though chemically unreactive gas. It has an atomic number of eighty-six, which corresponds to the number of protons found in the nucleus of any isotope of radon. There are more than thirty known isotopes of radon, and each one emits some combination

of alpha, beta, and gamma radiation when undergoing radioactive transformation, commonly referred to as "decay."

 

Radon gas is ubiquitous in the natural environment. This is because the precursors to radon, such as the aforementioned radium isotopes, and others such as radium, thorium, and uranium isotopes, are present in some rock formations. Radon is also found in the man-made environment because many of the materials, consumer products, and foodstuffs of everyday life come from the naturally radioactive environment.

Radon is one of the few examples in nature of a gaseous element that results from the decay of a solid element and then decays into another solid element. This increases its potentially harmful effect in humans. For example, radon-222, the most common isotope of radon, is a product of the alpha decay of radium-226 atoms, found in rocks. Radon-222 atoms subsequently produce polonium-218 in a similar alpha-decay process, and it is this solid substance that can lodge in human tissue.

Solid-state radionuclides remain where created by decay processes unless they are redistributed by dissolving in groundwater or by becoming airborne. Given the chemically inert nature of radon, there are no known compounds that include this element. Thus isotopes of radon may diffuse away from their place of origin and usually end up dissolved in ground water or mixed with air above the soil and rocks that bear their solid precursors.

People's exposure to radon primarily occurs when radon seeps out of air spaces above soil or rocks and into surrounding indoor or outdoor air, such as the basements of houses built over radium-bearing rocks. It is not exposure to radon gas that actually may lead to harm, but exposure to the decay products of radon, specifically the ones with short half-lives that emit alpha radiation. Radon-222 offspring, like polonium-218 and polonium-214, become attached to dust particles that may be breathed in by people exposed to the gas and become lodged in the respiratory tract. Decay of the radon progeny while in the lungs is the means by which the radiation dose is delivered to the lungs. This dose, which is the energy of alpha particles absorbed by cells that line the lungs, is what gives rise to the potential for lung cancer associated with exposure to radon.

Radon has been labeled by the Environmental Protection Agency as the second-leading cause of lung cancer in humans (after tobacco smoke), based on mathematical risk estimates derived from many published studies of exposure of subsurface uranium miners to highly elevated levels of the gas, primarily radon-222. Many radiation health scientists have challenged such findings because of the vast difference in exposure levels between homes and buildings on the one hand, and subsurface mines on the other. However, a variety of action levels and exposure limits for radon gas exposure have been recommended or set into law for the protection of the public. The Surgeon General and the EPA recommend that radon levels of four picocuries or more inside homes be reduced. The EPA states that radon levels less than four picocuries still pose a risk, especially for smokers.

Methods to both detect and mitigate indoor radon exposure have been devised as well. Detection and measurement methods usually make use of a device to collect radon gas atoms or the offspring particles. The simplest real-time method would be a "grab sample," in which air is drawn into an evacuated flask that is then taken back to a laboratory for analysis. The most popular short-term measurement device is the activated charcoal canister, a small container of steam-treated charcoal that is opened and left at the sampling location for a prescribed time. Radon is adsorbed by the charcoal, and the decay products of the radon are analyzed after the canister has been resealed and retrieved. The simplest mitigation methods include sealing cracks and penetrations through foundations, as well as diverting the radon away from the slab or out of the ground, with vacuum or ventilation systems.

National Research Council. (1999). Health Effects of Exposure to Radon—BEIR VI. Washington, D.C.: National Academy Press.

National Research Council. (1999). Risk Assessment of Radon in Drinking Water. Washington, D.C.: National Academy Press.

(PRIMAPRESS) ROMA - Giampaolo Giuliani, l'uomo denunciato per procurato allarme perchè aveva predetto un terremoto nell'aquilano di grosse proporzioni sbagliando solo di qualche giorno non è l'unico al mondo a sostenere che attraverso le misurazioni del gas Radon si possa capire se si sta per verificare un evento sismico di grosse proporzioni. Da Taiwan agli Stati Uniti e persino dagli esperti giapponesi tale sistema di predizione dei terremoti da anni è consuetudine. Da numerosi studi compiuti in zone soggette a rischio sismico è emerso che anomalie nella concentrazione di Radon nel sottosuolo o nelle acque provenienti da circolazione profonda (o comunque da sorgenti) sono correlate a variazioni della tensione della crosta terrestre, e quindi ad eventi di natura sismica. Negli ultimi anni la ricerca sulla predizione dei terremoti è basata principalmente sull'osservazione di fenomeni precursori. In ogni caso la correlazione tra terremoti e fenomeni precursori è ancora di difficile comprensione perché le cause e le condizioni che li governano sono complesse. Un fenomeno precursore avviene prima che abbia luogo un terremoto. Negli ultimi venti anni sono stati individuati diversi precursori, tra i quali la variazione anomala di Radon nel suolo o in sorgenti situate presso le faglie attive. Il problema più importante che si presenta nel momento in cui si decide di eseguire una campagna di monitoraggio del Radon è quello di distinguere il segnale anomalo dall'emissione caratteristica del luogo o addirittura dal rumore. A questo punto non si capisce bene se l'allarme, lanciato da Giuliani da diversi giorni, seppur indirizzato alla città di Sulmona, potesse essere preso meno sottogamba. Sta di fatto che oggi si assiste ad un tentativo di minimizzare il precedente e di prevenire le polemiche, dicendo che nessun terremoto si può prevenire. Eppure basta cliccare sulla frase "Radon earthquake precursor" da qualsiasimotore di ricerca per rendersi conto che è una pratica diffusa in tutto il Globo e studiata e pErsino dalle nostre Università. Insomma non basta dire che non si poteva prevedere e punto. Lo dimostrassero Boschi, Barberi e gli altri che Giuliani ha colto nel segno per una coincidenza. Fortuna? Coincidenza? Oppure ben altro. I morti di oggi, se non subito, ma nelle prossime settimane, obbligano alla massima chiarezza su questa vicenda. E di sicuro Bertolaso non si sottrarrà ad una simile e necessaria verifica. (PRIMAPRESS)

E ora in Turchia... anche loro lo sanno...

Geochemical monitoring in the Marmara region (NW Turkey): A search for precursors of seismic activity

Physical and chemical properties of warm and hot spring waters as well as soil radon concentrations were measured continuously during a 3-year period in the Marmara region; following the devastating İzmit earthquake of 17 August 1999 (Mw = 7.4). Promising and encouraging anomalies in ground radon emanation have been recorded and found to be closely related to seismic activity. The temporal and spatial variations in the soil radon data are presented. The earthquakes with magnitude >4 in the region were correlated with positive radon anomalies. Furthermore, during quiescence (absence of seismic activity) the radon data indicate random walk behavior of radon in soil and show Rayleigh-type probability density function (pdf), however, during the earthquake build-up period, the data show deviations from Rayleigh-type pdf. The radon positive anomalies indicate disturbance of the path of gas movement or gas release pattern prior to earthquakes. However, systematic and consistent anomalies in physical and/or chemical properties of the spring waters have not been detected for earthquakes occurring in the observation period (M < 5.3).

Radon variation in Spring water before and after Chamoli earthquake,
Garhwal Himalaya, India
V.M.Choubey1, P.K.Mukherjee1 and R.C.Ramola2
1Wadia Institute of Himalayan Geology,
33-General Mahadeo Singh Road, Dehra Dun 248 001, India
E-mail: vchoubey1@rediffmail.com
2Deptt. of Physics, H.N.B. Garhwal University
Badshahi Thaul Campus, Tehri Garhwal -249 001, India
Abstract: Radon levels measured in the springs of Garhwal Himalaya before Chamoli earthquake (Mb=6.3) of 29th March 1999 have been compared with those measured immediately after the earthquake. Radon levels in springs are found to be affected by the earthquake. In all, 39 springs were surveyed for radon which are divided into three sections across the tectonic structures. The water discharge data were also collected along with the radon measurements. The enhancement/reduction of radon concentration and water discharge rate due to the earthquake were then calculated as Enhancement Factor (EF) by dividing the pre-earthquake data with that of post earthquake measurements.
It was observed that in most of the springs, radon levels were increased after earthquake, particularly in worst effected area around Rudraprayag and Chamoli. Most of the springs that shows very high radon enhancement factor are located in the vicinity of the different thrusts or maximum damaged zones with development of ground cracks. Probably, this may indicate that there is some degree of reactivation of the thrusts and faults present in this area. However, radon concentrations were found to decrease in some springs. A similar variation was also recorded for water discharge from the spring before and after earthquake. The enhancement or depletion of radon and water discharge may be attributed to the combined effect of earthquake and partial mobilization of major structural features and lithological control. The effects of earthquake on the radon levels and water discharge in the springs of Garhwal Himalaya are discussed in details. A significant change in radon concentration and water discharge from springs in relation to earthquake shows that this technique may be exploited as an additional tool in earthquake prediction programme in Himalayan region.
1. Introduction
Himalaya on the whole is experiencing vertical uplift and transverse block movements due to the compression caused by the continued northward movement of Indian plate, even after the collision with Eurasian plate in Paleocene. As a result Himalayan region is being subjected to intense neotectonic movements and seismic activities. During the past decade two major earthquakes struck in the Garhwal region alone, i.e., Uttarkashi earthquake [1] of 20th Oct. 1991 (6.5 Mb) and Chamoli earthquake of 29th March 1999 (Mb 6.3).
Geological and geohydrological signals preceding significant earthquakes have been reported since 1960 and have been successfully used to predict the 1975 Haicheng earthquake in China. Perceptible changes in terms of hydrogeochemical parameters have been reported as a result of noticeable earthquake by several researchers [2,3,4]. During the last few decade, research on seismic precursors have received enormous attention because of its social relevance. Monitoring of radon concentration has been found to have a great potential as a reliable precursor [5]. Several studies suggests that the emanation pattern of radon (222Rn) can be used as a seismic precursor [6,7,8,9,10,11,12,13]. The emanation of radon increases with the increasing stress in the crustal layers, and supposed to be maximum with the release of the stored energy in the form of the earthquakes. Okabe [14] suggested that along the active faults the soil gas is enriched in radon. The concentration of radon in air, groundwater, springs close to the epicentral area increases anomalously just before an earthquake [15]. Deep seated faults and fractures provide passages for migration of the gas and thereby concentration of radon dissolved in water increases. King [7] observed high radon near active faults and temporal variations related to large seismic activity.
The authors have been engaged in radon studies in soil and spring water for last few years in Garhwal
1
region. The data related to radon concentration in spring water were collected from this area in May 1997. This forms the basis for the present study to compare the variation in radon concentration before and immediately after the recent Chamoli earthquake. However, a more meaningful information can be extracted if the radon monitoring along with metrological parameters are made on a continuous basis.
2. Chamoli Earthquake
A moderate earthquake (Mb 6.3) occurred on the early morning hours of March 29, 1999 (0:35:13:59 IST) near the Chamoli town in the higher reaches of Garhwal Himalaya. Apart from the loss of lives and properties, several ground cracks have been observed in the vicinity of the area extending up to 1 km in length with width varying from 10 to 20 cm. Numerous rock falls, landslides and changes in ground water levels have also been observed after the Chamoli earthquake. This earthquake was followed by several aftershocks with magnitude up to 5.3. According to a US Geological Survey report pure thrust faulting has been said to be the cause of the earthquake based on the focal mechanism solution. It is worthwhile to mention that the area is traversed by several thrust/faults. Preliminary reports by the scientists of wadia Institute of Himalayan Geology suggest that the activation of Main Central Thrust (MCT) seems to be the causative factor for this tragic event. It may be noted that the Uttarkashi earthquake of 1991 (Mb 6.5) was also caused by re-activation of MCT.
3. Radon Measurement Technique
750 ml of water sample was taken in radon-tight reagent bottle of 1 litre capacity connected in a close circuit with a ZnS coated detection chamber through a hand operated rubber pump and a glass bulb containing CaCl2 to absorb the moisture. Air was then circulated in close circuit for a period of ten minutes till the radon formes a uniform mixture with the air and the resulting alpha activity was recorded. The details of the measurement technique are discussed elsewhere [16,17]. The calibration factor 1 count/min. = 0.11 Bq/l was used to convert the recorded alpha counts in Bq/l [18].
4. Regional geology and structure
The geological structure and tectonics of the area has been studied and published by several researchers [19,20,21,22,23,24,25], however a brief account of the salient features are presented here.
There are four principal lithotectonic units present in the area affected by the said earthquake. From south to north ther are: a) Deoban-Barinag formation b) Rautagara formation c) Ramgarh formation and d) Munsiari formation (Fig.1). The crystalline rocks of Ramgarh formation forms a large antiformal structure overlying the Lesser Himalayan rocks. Deoban and Berinag formations to the north of Chamoli are highly folded into a antiform referred to as Chamoli Window and inturn is thrusted over by the rocks of the Munsiari formation. The contact between the two is tectonic, called Munsiari Thrust [20] or MCT-1 [26,23].

Figure 1 Geological map of Chamoli area [after 20] showing radon sample locations.
5.

In any comparative earthquake related studies on the variation of radon concentration, it is important to consider the effects due to seasonal and meteorological parameters. In the present case, however, the data collection period is almost the same (i.e. May 1997 and April 1999 ). As such, the effects due to the above factors may be considered as minimum. In order to have a meaningful comparison, the post earthquake data were collected from the same springs, as many of them as possible, at the same spot using the same instrumentation and methodology. In addition, few more springs were also sampled in the vicinity of the maximum damage zone. In all, 39 springs were surveyed which are divided into three sections across the tectonic structures. (Table 1 and Fig.1). The water discharge data were also collected along with the radon measurements (Table 1). The enhancement/reduction of radon concentration and water discharge rate due to the earthquake were then calculated as Enhancement Factor (EF) by dividing the pre-earthquake data with that of post earthquake measurements.
Along Karanprayag–Chamoli-Pakhi section a total of 11 springs were sampled out of which only 7 were repeated in the post earthquake measurements. Pre-earthquake radon concentration varies from 0.9 to 35.9 Bq/l as against 1.2 to 68.4 Bq/l observed in post-earthquake measurements with enhancement factors as high as 26.4 to the lowest of 0.6. It was observed that most of the radon concentration in spring water and its discharge in this section has invariably shown some degree of enhancement (EF>1) or remain nearly same, except for the sample no. eq22. In case of sample no. eq21, however, the discharge rate has drastically decreased where as the radon enhancement factor is highest (26.4 ). It may be noted that these two samples lie on the faults of local extent north of Karanprayag as well as a series of thrusts within a few km radius. The next highest radon value as well as discharge is also located in the same area (sample no. eq23).
The post earthquake radon measurements in 8 springs out of total 12 samples were carried out along Rudraprayag-Tilwara-Okhimath section. The comparison of the data reveal that the springs located in Tilwara (sample nos.eq3 & eq4) and Kund (sample nos.eq8 & eq9) show significant enhancement (EF=3.6 to 26.4) of radon as well as the water discharge (EF=1.2 to 3.0) with the exception of sample no. eq3 where the water discharge has decreased by a factor of 0.25. All these four samples are again related to the same fault extending from north of Karanprayag and a couple of thrusts (Fig.1). The Bhatwari–Garponga area witnessed notable damage and therefore additional measurements were made in the springs of this section (sample nos. eq11 to eq16). Despite of the increase in water discharge from springs (samples eq6, eq7 and eq7a), the radon concentrations were found to decrease. The
reason for the depletion in spite of considerable damage, may be attributed due to escape of radon gas through the ground cracks that has developed extensively in this area during the earthquake. The additional water, transported from non-mineralised area, may also have diluted the radon content in the springs.
In Ghonti-Budhakedar section there are 7 springs that have been sampled in both the sessions. In addition 5 more springs were also studied during the post-earthquake measurements. Strikingly the discharge of the spring have drastically reduced up to 75% (EF= 0.25 to 0.8) after the earthquake. Some of the previously studied springs have even completely dried up. The earlier measurements indicate that radon levels were moderately high, ranging between 4.5 to 159.7 Bq/l. In the later measurements, three of the springs recorded marginal increase in radon contents (EF=1.03 to 3.9), where as rest four springs have recorded the lower radon concentrations (EF= 0.2 to 0.5). The drying up of the springs and reduction in discharge rate suggest that there has been a narrowing down of the fractures and weak planes. This has resulted in the lower radon concentration due to obstruction of the radon pathways in the subsurface strata. Probably the emanation of radon might have increased many fold due to the earthquake in this area, however, it was restricted to a large extent due to the same reason. It is worthwhile to mention here that uranium concentration in the granitic rocks near Budhakedar is high [27] and there are reports of uranium mineralisation also [28]. The anomalous radon concentrations of 622.9 and 106.8 Bq/l, observed in sample nos. eq52 and eq53, respectively, near Budhakedar, are due to the combine effect of the earthquake and high uranium contents of the rocks.
 

Conclusions
The behavior of radon gas built up and its release before and after the earthquake has been well established based on the continuous monitoring. However, in the present case, we have collected data a couple of year prior to the Chamoli earthquake and about a week later. Although we have lost many of the signatures in between, nevertheless there are some valuable information that can be drawn based on the available data:
1. There is an overall enhancement of radon concentration in the post-earthquake measurements of the spring along Karanprayag-Chamoli-Pipalkoti and Rudraprayag-Ukhimath sections. This is also coupled with an increase in water discharge rate of springs.
2. The discharge rate of the springs in Budhakedar-Ghonti section have considerably diminished, even some of the springs studied earlier have dried up. Except for two springs towards north near Munsiari thrust, the radon concentration have decreased after the earthquake.
3. Most of the springs that shows very high radon enhancement factor are located in the vicinity of the different thrusts or maximum damaged zones with development of ground cracks. Probably, this may indicate that there is some degree of reactivation of the thrusts and faults present in this area.
The above mentioned results and discussion reveal that there is a significant relationship between ground water radon concentration and earthquakes. Therefore, continuous measurements of radon both in soil-gas and spring/groundwater at several sites in a well organised grid pattern are necessary to make it as useful precursor in earthquake prediction programme. In addition, the combination of radon measurements with the measurements of CO2, CH4, N2, H2, He, Cl and F along with the temperature and electric conductivity of the spring water may prove a useful hydrogeochemical precursor for any destructive earthquake.
References
1. Thakur, V.C. and Kumar S. 1994. Seismotectonics of the 20 October 1991 Uttarkashi earthquake in Garhwal, Himalaya, north India. Terra Nova 6, 90-94.
2. Roeloffs, E.A., 1988, Hydrologic precursors to earthquakes–A review: Pure and Applied Geophysics, v. 126, p. 177-209.
3. Sharma, S.C. 1993. Isotope and geochemical studies on the thermal discharge in the Alaknanda and Bhagirathi Valleys, Chamoli and Uttarkashi district. Geol. Surv. Indi. Recd., 126, pt. 8, 211-213.
4. Sharma, S.C. 1997. Thermal springs as gas monitoring sites for earthquake prediction. In: Rare gas Geochemistry (Ed. H.S. Virk), GND Univ., Amritsar, pp. 193-199.
5. Braitenberg, C., Garavaglia, M., Magnoni, S. and Zadro, M. 1997. Radon and crustal deformation monitoring in seismic area of Friuli, north-east of Italy. In: Rare gas Geochemistry (Ed. H.S.Virk), GND University, Amritsar, pp. 143-157.
6. Fleischer, R.L. and Mogro-Campero 1979. Radon enhancement in the earth, evidence for the intermittant upflow. Geophys. Res. Lett. 6, 361-364.
7. King, C.Y. 1978. Radon emanation on San Andreas Fault. Nature 271, 516-519.
8. King, Y.C., Koizumi, N. and Kitagawa, Y. 1995. Hydrogeochemical anomalies and the 1995 Kobe earthquake. Science, 269, 33-34.
9. H. Friedmann, F. Hernegger: A Method for Continuous Measurement of Radon in Water of Springs for Earthquake Prediction. Geophys. Research Letters 5, No.7, p.565-568, 1978.
10. Wakita, H., Nakamura, Y., Notsu, K., Noguchi, M. & Asada, T. 1980. Radon anomaly: A possible precursor of the 1978 Izu-Oshima-Kinkai earthquake. Science, 207, 4433: 882-883.
11. Ramola, R.C., Singh, M., Sandhu, A.S., Singh, S. and Virk, H.S. 1990. The use of radon as an earthquake precursor. Nucl. Geophys. 4, 275-287.
12. Virk, H.S. 1993. Radon and earthquake prediction in India : Present Status. Nucl. Tracks Radiat. Meas. 22(1-4), 483-494.
13. Virk, H.S. 1995. Radon monitoring of microseismicity in the Kangra and Chamba Valleys of Himachal Pradesh, India. Nucl. Geophys. 9(2), 141-146.
14. Okabe, S. 1956. Time variation of atmospheric radon content near the ground surface with relation to some geophysical phenomena. Memoir Sci. Coll. Univ. Kyoto Sev. A28, 99-115.
15. Duarah, R., Baruah, S., Yadav, D. and Sitaram, M.V.D. 1997. Anomalies in radon emanation as earthquake precursor phenomena. In: . In: Rare gas Geochemistry (Ed. H.S.Virk), GND University, Amritsar, pp.175-187.
16. Ramola, R.C., Rawat, R.B.S., Kandari, M.S. and Choubey, V.M. 1997. Measurement of Radon in drinking water and indoor air. Radiation Protection Dosimetry, 74(1/2), 103-105.
17. Choubey, V.M. 1998. Radon measurements in soil and water and its relation with geology, Garhwal Himalaya. In Proceeding volume of the 7th Tohwa International Symposium “Radon and thoron in the human environment” Japan, (Eds. A.Katase and M. Shimo) World Scientific Publisher, pp. 193-198.
18. Singh, M., Singh, N.P., Singh, S and Virk, H.S. 1986. Calibration of radon detectors. Nucl. Tracks-Radiat. Meas., 12, 739-742.
19. Gansser, A. 1964. Geology of the Himalaya. Inter Science Publication. John Wiley & Sons, 289p.
20. Valdiya, K.S. 1980.Geology of Kumaun Lesser Himalaya. WIHG Publ. Dehra Dun.291p.
21. Saklani, P.S. and Bahuguna, V.K. 1986. Thrust tectonics of the Garhwal Himalaya. In: Himalayan thrust and associated rocks (Ed. Saklani. P.S.), Today & Tommorow Pub. New Delhi, 1-25.
22. Roy, A.B. and Valdiya, K.S. 1988. Tectonometamorphic evolution of the Great Himalayan Thrusts sheets in Garhwal region, Kumaun Himalaya. Jour. Geol. Soc. India. 32, 106-124.
23. Bist, K.S. and Sinha, A.K. 1990. Imbrication of Main Central Thrust and its significance in the central sector of the Himalaya. In: Crustal evolution and orogeny, Ed. S.P.H.
Sychanthavong., Oxford & IBH Publ. Delhi. 293-302.
24. Thakur, V.C. 1992. Geology of Western Himalaya. Physics and Chemistry of the earth. Pergamon press, 19, (I-VI), 1-363.
25. Metcalfe, R.P. 1993. Pressure, temprature and time constraints on metamorphism across the Main Central Thrust Zone and High Himalayan slab in the Garhwal Himalaya. Himalayan Tectonics. Geol. Soc. Special Publication , 74, 485-509.
26. Fuchs, G. and Sinha, A.K. 1978. The tectonics of the Garhwal Kumaun Lesser Himalaya. Jaharb. Geol. B-A ,121, 219-241.
27. Choubey, V.M., Bist K.S., Saini N.K. and Ramola, R.C. 1999. Relation between soil-gas radon variation and the different lithotectonic units of Garhwal Himalaya, India. Applied Radiation and Isotopes. 51 (5), 587-592.
28. Nashine, S.K., Dhanaraju, R., Bhatnagar, G.S. and Narayan Das, G.R. 1982. Uranium occurrences close to Main Central Thrust around Sileth, Dhargaon, Chamyala Balganga Valley, Tehri Garhwal (U.P.). Him. Geol. 12., 305-316.

Table 1. Insitu measurements of the Radon concentrations and Water discharge data during post (April, 1999) and pre(May, 1997)
Chamoli earthquake of Garhwal Himalaya (March 29, 1999)
Radon
Water Discharge
Post-EQ
Bq/l
Pre-EQ
Bq/l
Pre-EQ
Lt/min.
Post-EQ
Lt/min.
Rn-EF
Dis-EF
1. Karanprayag-Chamoli-Pakhi
eq21 68.4 2.6 6 2.5 26.38 0.42
eq22 1.9 3.2 1 0.25 0.59 0.25
eq23 4.5 0.9 2.5 6.0 5.17 2.40
eq24 8.0
eq25 4.2
eq26 27.8 24.3 12 20.0 1.15 1.67
eq28 59.4 35.9 3 5.0 1.66 1.67
eq29 11.3 4.2 5 6.5 2.71 1.30
eq30 1.2
eq31 2.1
eq32 3.0 2.3 6 8.0 1.28 1.33
2.a) Rudraprayag-Tilwara-Okhimath
eq1 0.6
eq3 15.0 3.5 2 0.5 4.27 0.25
eq4 22.8 0.9 15 40 26.44 2.67
eq5 2.1
eq6 56.9 60.0 3 4 0.95 1.33
eq7 6.1 10.2 1.5 3 0.61 2.00
eq7a 5.0 39.7 6 10 0.13 1.67
eq8 219.2 60.4 3 9 3.63 3.00
6
eq9 28.2 7.8 3 3.5 3.63 1.17
eq10 3.4
eq17 25.5 25.8 16 20 0.99 1.25
eq18 4.5
2.b) Subsection: Bhatwari-Garponga
eq11 19.6
eq12 14.6
eq13 97.8
eq14 3.9
eq15 162.4
eq16 61.6
3. Ghonti-Ghansali-Budhakedar
eq41 2.2 4.5 3 2.5 0.50 0.83
eq42 18.5 18.0 4.5 3 1.03 0.67
eq45 5.4 19.5 4.5 1.5 0.27 0.33
eq47 1.9 6.5 0.5 0.29
eq48 1.3 1
eq49 31.9 87.3 3 0.37
eq50 22.6
eq51 10.3
eq52 622.9 159.7 10 2.5 3.90 0.25
eq53 106.8 66.7 6 2 1.60 0.33
EQ=Earthquake, EF=Enhancement factor, Dis.=Spring discharge, Rn=Radon

Misurazioni sull'acqua. Ottimo.

Earthquake Precursors 129 * E-mail: chaddha@ngri.res.in
Precursors for certain earthquakes indicate that, there is reason enough to believe that precursors do exist and
can be observed for forecasting earthquakes. The concept that catastrophic damage is the culmination of smaller
events provides a physical basis for the presence of precursors. In case of an earthquake rupture, certain
precursory activity can be expected, if the observations are made in the near vicinity of the causative rupture. Even
if they are genuine, the precursors will remain circumstantial, unless linked to earthquake rupture through a sound theory.
And for this to happen, a study of precursors in seismically active regions is needed. Application of chaotic dynamics
approach to certain seismically active regions offers some optimism.


Globally, efforts to predict earthquakes were started about a century ago and peaked during 1970s, when there
was a cautious optimism that reliable earthquake prediction might be possible. This optimism was mainly due to a few
successful attempts, like the one for August 3, 1973, Blue Mountain Lake earthquake of M 2.6 in New York based on
the precursory change in the Vp/Vs ratio. The observations appeared to be consistent with the Dilatancy Hypothesis,
which explains the behaviour of crustal rocks under nearcritical stress levels. This attempt was followed by another
spectacular and accurate prediction of the February 4, 1975 Haicheng earthquake of M 7.3 in China based on several
precursors including immediate changes in foreshock escalation and, clustering and abnormal animal behaviour
in the region. This timely prediction saved about 1,00,000 human lives. But, the occurrence of M 7.6 earthquake in
Tangshan, China, 14 months later, for which no prediction could be made, attested to the complexity of earth processes
prior to the occurrence of large earthquakes. Later, when the predicted earthquake did not occur within the specified
time frame (1985-1993) in the Parkfield prediction experiment in the United States on San Andreas Fault, the collective
optimism of the scientific community turned into collective pessimism. This pessimism was further compounded when
a similar effort in Japan, where a prediction had been made for the Tokai region did not yield positive results. An
earthquake did occur in 2003 in the Tokai region but no definite precursors were reported prior to this event. Since
then, there are some success stories of definite precursory changes prior to several moderate to large earthquakes.
These are region specific, for example, shear-wave splitting and in-situ strain studies in Iceland. Russian seismologists
have been working for several years on statistical and other related approaches, but have not succeeded in short-term
forecasts.

With the increase in the reporting of precursors, the International Association of Seismology and Physics of the
Earth’s Interior (IASPEI) appointed a sub-commission on Earthquake Precursors
R.K. CHADDHA*
National Geophysical Research Institute, Uppal Road, Hyderabad - 500 007
Earthquake Precursors for their validation. The subcommission defined a precursor as a quantitatively
measurable change in an environmental parameter that occurs before main shock and is thought to be linked to
the preparation of the main shock. Enhanced foreshocks, seismicity pattern, hydrological and radon gas have been
considered as few of the notable precursors in the recent times. Subsequently, there has been a debate amongst
seismologists on earthquake prediction. During the 22nd General Assembly of the International Union of Geodesy
and Geophysics in 1999 at Birmingham, UK, the IASPEI sponsored a symposium on “Nature of Seismic Sources and
Prediction of Earthquakes”. The proceedings, giving the status of earthquake prediction related research were
published in 2001 in Tectonophysics. In a joint UNESCOEOLSS publication on “Encyclopedia of Life Support System”
in 2003, the history on earthquake prediction and provided a glimpse on the present status of earthquake prediction
(Gupta, 2003). Earlier, in a Special Issue of Current Science in 1992 on “Seismology in India – An overview” the status
of earthquake prediction in India was discussed in one of the papers. All these attempts point towards a strong desire
of the scientific community to document and understand the physical phenomena of earthquakes.
INDIAN SCENARIO
In India, earthquake precursory studies are being attempted at a few research/academic institutions. DST has
been continuously supporting research projects related to the precursory studies based on Radon and Helium
measurements, ULF/VLF radio emissions; GPS-based deformation studies, Borehole observations such as water
level, temperature and other physical parameters in Koyna region, etc. The notable have been the geomagnetic field
precursory studies in Himalaya and in Koyna by the Indian Institute of Geomagnetism, Mumbai. Geo-electrical studies
for precursors were studied in the Shillong region. ULF/ VLF magnetic field emissions associated with earthquakes
are in progress since February 1998 at the R.B.S. College Agra. The Saha Institute of Physics, Kolkata has been
regularly monitoring the geochemical constituents He Rn, CH4 etc., at Bakreswar spring and observed anomalies
related to the 26th December, 2004 Sumatra earthquake. However, two systematic and notable studies in
earthquake precursors are the Northeast India in Himalaya and Koyna region in Maharashtra. The first study deals with
a medium-term forecast of an M 8 ±½ based on seismicity pattern, which was considered a moderate success. The
IASPEI Sub-commission on earthquake prediction also reviewed this study. In the Koyna region various
geophysical, geochemical and hydrological precursors are  being continuously monitored, as this region is a natural
laboratory with earthquakes continuously occurring for the last 40 years.

130 Glimpses of Geoscience Research in India
NORTH-EAST INDIA
On the basis of a study in 1986, which examined the time series of earthquakes associated with several of the
moderate earthquakes that occurred in the area between 200 N to 320 N latitude and 870 E to 1000 E a medium-term
forecast was proposed for a possible future earthquake of magnitude 8 + ½ earthquake in the NE region (210 N to
25 ½0 N and 930 E to 960 E) before 1990. This forecast was based on the concept that moderate magnitude to great
earthquakes in the north-east India region are found to be preceded, generally by precursory swarms and quiescence
periods. The medium-term forecast came true with a magnitude 7.3 earthquake occurring on August 6, 1988. The
epicentral coordinates for this earthquake were 25.1490 N 95.1270 E and were within the proposed region of forecast.
 

KOYNA REGION EARTHQUAKE PRECURSORY STUDIES
A multi-parametric study has been initiated by National Geophysical Research Institute, Hyderabad to monitor
seismological, geophysical, geochemical and hydrogeological precursors in Koyna region. Seismological
precursors are being monitored by a broadband network.
The geophysical precursory studies like gravity, magnetotellurics, electrical resistivity etc., are also in progress. The
seismological precursory study in Koyna region has led to the proposal of a nucleation method. The proposed method
is based on the premise of detecting and observing a nucleation zone of foreshock cluster, which increases over
a time window of 50 to 100 hours period prior to a M~4 to 5 earthquakes. Studies on the deepening of such a zone
to foresee a hypocentre of a future moderate-size earthquake at the base of the seismogenic layer. Three successful
predictions of moderate-size earthquakes in the Koyna region during 2006-07 have been made by National
Geophysical Research Institute base on the nucleation concept. Precursor studies for monitoring well water levels
in deep bore wells in the Koyna region have resulted in the identification of pre and co-seismic precursors for
earthquakes of magnitudes > 4 occurring in this region.

Radon is being continuously monitored at three sites in the vicinity of the epicentral region in Koyna. The analysis of
data showed a sudden increase in Radon anomaly by a factor of two, about 19 days prior to 14th March 2005
earthquake (M 5.1) in the region. The Ministry of Earth Sciences, Government of India has taken a major initiative
to establish Multi-parametric observatories to monitor precursors in certain seismically areas in the country.
Two such observatories are already in operation at Ghuttu in Himlaya, being operated by the Wadia Institute
of Himalayan Geology, Dehradun and the other at Koyna, operated by the National Geophysical Research Institute,
Hyderabad. Similar, observatories are being planned in Andaman and Nicobar Islands and Northeast India.


REFERENCES
Chadha, R.K., Pandey, Ajeet P. and Kuempel, H.J., 2003, Search
for earthquake precursors in well water levels in alocalized
seismically active area of Reservoir Triggered Earthquakes
in India: Geophys. Res. Lett., v. 30 (7), p. 1416, doi:10.1029/
2002GL016694.
Das, Nitish, K., Bhandari, Rakesh, K., Ghose Debasis, Sen
Prasanta and Sinha Bikash, 2005, Anomalous fluctuation
of radon, gamma dose and helium emanating from a
thermal spring prior to an earthquake: Current Science,
v. 89(8), p. 1399-1404.
Gupta, K.H. and Singh, H.N., 1986, Seismicity of the north-east
India region: Jour. Geol. Soc. India, v. 28, p. 367-406.
Gupta, H.K., 1992, Seismology in India: An overview: Current
Science, v. 62 (1, 2).
Gupta, H.K., Chadha, R.K. and Srinagesh, D., 2001, The Nature
of Sources and Prediction of Earthquakes: Tectonophysics,
v. 338(3,4), p. 357.
Gupta, H.K., 2003, Nature of earthquakes: Publishers Joint
UNESCO-EOLSS, p 58.
Gupta H.K. et al., 2006, An earthquake of M~5 may occur at
Koyna: Curr. Sci., v. 89 (5), p. 747-748.
Gupta, H.K. et al., 2006, Prediction of an M~4 earthquake in
the Koyna region comes true!: Jour. Geol. Soc. India,
v. 68, p. 149-150.
Gupta, H.K. et al., 2007, Earthquake forecast appears feasible
at Koyna India: Curr.Sci., v. 93(6), p. 843-848.
Reddy, D.V., Sukhija, B.S., Nagabhushanam, P. and Devender
Kumar., 2004, A clear case of radon anomaly associated
with a micro-earthquake event in a Stable Continental
Region: Geophy. Res. Letters, v. 31, L10609.
Singh, B., Kushwah, V., Singh, V., Tomar, M. and Hayakawa,
M., 2005, Simultaneous ULF/VLF amplitude anomalies
observed during moderate earthquake in Indian region,
Indian: J. Radio & Space Phys., v. 34, p. 221.
Singh, R.P., Singh, B., Kushwah, V.K. and Chauhan, R.V.S.,
2004a, Attenuation of ULF-VLF seismo- electromagnetic
signals and their propagation to long distances: Indian
J. Radio & Space Phys., v. 33, p. 189.
Singh, B., Kushwah, V., Singh, O.P., Lakshmi, D.R. and Reddy,
B. M., 2004b, Ionospheric perturbations caused by some
major earthquakes in India: Phys. Chem. Earth, v. 29,
p. 537.

 

Jedi Simon

             Ora scusatemi ma c'e un 4,2 in corso e continuo la ricerca appena posso.... vado....

 

Passata per ora...fiuuuu....

Ora associamo all'analisi dinamica del gas, l'analisi del campo elettro-magnetico per consentire un controllo incrociato di dati, comparati
al fine di definire con maggiore certezza e precisione riguardo al luogo e all'ora nella quale avrà luogo l'evento distruttivo. Con grande ringraziamento
a
Dionysios Dimakos (SV1ZJ)

   EARTHQUAKE PRECURSORS
By Dionysios Dimakos (SV1ZJ)

I am 38 years old and I reside in Athens, Greece.  My major is electronic engineering.  Since 1991, I have been involved – as an amateur - in monitoring and recording the magnetic field variations, aiming to trace earthquake precursors.  For this purpose, I have constructed and tested many different circuits and sensors. It is a fact that the major problem during measurements is created by the use itself of electric devices, as well as the power lines.
 


Scope of Measurement

During the past few years, I have focused my study on 0,1 Hz to 40 Hz frequency spectrum, measuring and recording the instant variations of the magnetic field and not its absolute value.
Unfortunately, in the vicinity of big cities, the levels of signals from natural sources are 1.000-30.000 times lower than those of signals transmitted by our technological civilization.  For this reason, the use of filters in the electronic circuit is imperative.
 


The Electronic Circuit

The electronic circuit diagram can be seen on Picture 1 and Picture 2.
It consists of amplifying stages and electronic frequency filters.


 

Its first stage consists of the operational amplifier LT1128. The alternative voltage (in the 50Hz frequency) of the LT1128 Pin 6 should not exceed +-2V (in reference to ground); alternatively, we reduce the value of R1.

The second stage is a filter “Notch” in the frequency of 50Hz and consists of two OP27.

The three variable resistors of the filter are 20-turns “trimmers”, adjusted to the values specified on the diagram, before placed on the circuit. When the circuit is in operation measuring the magnetic field, the “trimmers” are readjusted micro metrically on the original values, in order to achieve the lowest AC voltage in the frequency of 50Hz, at the circuit output.

From “trimmer’ R2, we should have 0 V DC at the circuit output, while from the “trimmer” of the output stage, we adjust the value of the output signal which we forward to the Line-In of our computer sound card.

The capacitors C1, C2, C3, C4 are of “Silver mica” or “foil polystyrene” type in order to maintain stable capacity during temperature variations.

All stable resistances are of “metal film” 1% type.

Having tested many different operational amplifiers combinations, I decided on the ones of the diagram since they offer lower noise levels and better overall performance on the specific circuit.


The Sensor

The sensor consists of a 50cm long ‘PVC’ pipe, with 5cm external diameter. We wrap three layers of 0.04’’ “MuMetal” sheet around the pipe. Over the “MuMetal”, we wrap insulating tape and coil enamel coating wire, 1 mm thick, at pipe's length. We coil ~9 layers of wire (~ 4.500 coils). When the coil is completed, we wrap it up with insulating tape and plate it with aluminum foil or self-adhesive aluminum tape.  The aluminum is connected to the plate of the coaxial cable (RG-213, RG-8…) that we use as signal transfer line.
 

  The end of the external winding of the coil is also connected to the shield of the cable, while the internal end is connected to the internal conductor of the coaxial cable. The coil in inserted into a 60cm long ‘PVC” pipe, with 10cm diameter it is then stabilized and waterproofed.

The coil has pure resistance R = 16 Ohm and self-induction L = 1 H.  The final construction is placed horizontally 10cm above the ground and stabilized.  The ends of the coil are adjusted at 110o and 290o respectively in reference to the magnetic north. . 

 

Just to emphasize that the coil should have low pure resistance (including the coaxial cable), in order to achieve reduced noise in the circuit


Recording Software

For the recording, I use the “spectrogram” with the following settings:


In order to save the images, I use program “2020” with the following settings:


Receiving Samples
 

28_4_2003  Record

Quake 29_4_2003, 5R  230km SW
  8_6_2003   Record

Quake  9_6_2003, 5R  230km NNW

 
6_7_2003   Record

Quake  6_7_2003, 5.4R  320km NE
  13_8_2003  Record

Quake  14_8_2003, 5.9R  275km WNW

 
8_2_2004  Record

Quake  9_2_2004, 4.6R  250km SW
  29_2_2004  Record

Quake  1_3_2004, 5R  175km SW

 
07_3_2004    Record

Quake 11_3_2004, 4.5R  Dodecanese Isl
  14_3_2004  Record

Quake 17_3_2004, 6.1R  South of Crete Isl
 

 


Findings

I have already completed an analysis of the records produced during the period from 01-02-2003 to 01-03-2004.  (During that period, the reception/ recording device has been unavailable for a total of 30 days due to technical reasons.)

According to Athens National Observatory, 12 earthquakes of 4,8R or above have been recorded for the same period. (Aftershocks that follow bigger earthquakes within a few days are excluded; earthquakes that occurred the days from 23-07-2003 till 26-07-2003 are also not recorded due to technical problem.)

Analyzing the remaining 10 earthquakes in relation to the registered precursors, we can draw the following conclusions:

Studying all precursors that have been recorded the past 14 months as well as the earthquakes that followed, we can draw the following conclusions:

1. The time of precursor occurrence is somehow related to the positions of the Sun and the Moon in relation to Earth and the gravitational powers generated. A signal possibly related to a forthcoming earthquake is twice as likely to appear within 72 hours after the actual time of each phase of the Moon (new moon, first quarter, full moon, last quarter), than within the 72 hours before.

2. The average time of occurrence is app. 1 hour before sunset

3. 75% of the earthquakes following the precursors occur within 72 hours after the precursors
 


Earthquake Precursors Pattern

The pattern of earthquake precursors could be described as “pre-earthquake signature” since every reoccurrence maintains common characteristics.  Observing the recorded signal, one realizes that, throughout the spectrum 1-40Hz, there are lines almost parallel to each other at any given time. Any spectral line differs 1-3Hz from the previous one.

The duration of the “pre-earthquake signature” varies from 15 minutes up to 3 hours, while for the specific spectral impressions, the intensity of the magnetic field varies from 2 PT to ~20 PT. It is most likely that the common characteristics of the signal will still appear in frequencies higher than 40 Hz.  Further observation is required in order to verify the following:

What is the correlation of the magnetic field intensity and the duration of the “pre-earthquake signature” with the data of the earthquake to follow?
What is the highest frequency in which the “pre-earthquake signature” is observed, as well as the correlation of the frequency with the earthquake magnitude or the distance from the epicenter to the reception device?


Conclusion

The intensity of the earthquake precursor will be increasing as one approach the epicenter of the imminent earthquake.  Consequently, earthquake prediction through the method described above requires a network of base-stations throughout the country, each station including at least two sensor devices and the respective recording equipment.
I believe that observation of the earthquake “signature” used in parallel with other methods can provide very satisfactory results in earthquake prediction.
Teamwork, selfless interest and publication of the progress and research findings of different scientific groups are the prerequisites for successful earthquake prediction.  Selfish motives can only cause further delay to a possible solution, which means more death and tragedy….
I hope that soon all scientists working in the specific field will cooperate and combine their knowledge in order to achieve a common goal, earthquake prediction, which is essential to humanity.


In the future, “pre-earthquake signatures” will be posted in this site, as soon as possible following occurrence.
We are always at your disposal to provide further clarifications and data to anyone interested in constructing the specific device. More detailed pictures are available at my web site: http://www.ham.gr/ulf/index_en.php  Qui puoi trovare gli spettrogrammi degli ultimi terremoti.
 

EARTHQUAKE  TRIGGERING, EARTHQUAKE  PRECURSOR  SIGNAL  GENERATION,
AND  EARTHQUAKE  SENSITIVITY  THEORY  PICTURES

Latest Update:   May 20, 2004

Earthquake Triggering     Earthquake Precursor Signal Generation
Earthquake Sensitive Humans, Animals, And Plants
Earthquake Triggering Forcess    Reference Information

Variabili da tenere presenti nel calcolo delle previsioni, per l'analisi della somma di tutti i vettori e di tutti gli elementi che possono prendere parte alla costruzione di un'onda sismica. 
 


       Site visitors are invited to make suggestions regarding additions to or corrections to these general theory pictures. Please send the suggestions to: seismic@ix.netcom.com.
       Many of the individual theories outlined here have been proposed by and discussed by researchers around the world over the years. 
       And one goal is to try to locate and include references to as many of those researchers’ publications as possible.


EARTHQUAKE  TRIGGERING

       The following theory pictures propose how earthquakes are being triggered.  Just a few of the most important earthquake triggering forces and phenomena are being discussed.  A number of others are listed in the following report:  Earthquake Triggering Processes.
 

A - 1   The movement of tectonic plates relative to one another and other processes cause strain to build in earthquake fault zones around the world.

A - 2   Various obstacles such as the resistance of rock layers to smoothly slip past one another cause that strain energy to build to dangerously high levels in some fault zones.

A - 3   Enough strain eventually accumulates in the fault zone rock layers to overcome those obstacles. And some of that stored energy is released in the form of a powerful earthquake.

A - 4   Two groups of forces, phenomena, and factors which can affect earthquake triggering are ones related to the gravitational pulls of the sun and the moon on the Earth, and fault zone environmental factors such as it north and south or east and west orientation.

A - 5   External forces or phenomena such as atmospheric storm related wind drag on continents and fault zone rock layer heating associated with solar and geomagnetic storms can temporarily or permanently add or subtract sufficient amounts of strain energy to or from the fault zone to move an earthquake occurrence time forward or backwards hours, days, or weeks in time.

A - 6   Forces or phenomena such as shockwaves from powerful earthquakes occurring around the world or nuclear tests can produce permanent instabilities in a fault zone and move the earthquake occurrence time forward or backwards days, weeks, months, decades, and even centuries in time.

A - 7   The actual earthquake occurrence time window is often finally determined by strain which is temporarily or permanently added to the fault zone by forces or phenomena which are directly and/or indirectly associated with the gravitational pulls of the sun and the moon. Some of those phenomena would be ocean tides and the Solid Earth Tide.

       The following drawing shows how forces related to the gravitational pulls of the sun and the moon on the Earth's crust can trigger an earthquake. A fault zone running east and west (EW) may fracture when the gravitational pulls of the sun and the moon pull it towards the north and the south at the same time. And a fault zone running north and south (NS) may fracture when the gravitational pulls of the sun and the moon pull it towards the east and the west at the same time.

Fault Zone Environment

A - 8   Two specific earthquake triggering force related theories are presently being considered. The first (see D – 1) involves earthquake triggering pressure being directly exerted on the Earth’s crust by the gravitational pulls of the sun and the moon. The second (see D – 2) involves earthquake triggering pressure being exerted on the crust by the Earth’s solid core as it attempts to move away from the sun – Earth – moon barycenter or center of mass which the core, the crust, the Earth as a whole, and the moon rotate around once each month.

A - 9   When forces or phenomena such as atmospheric and earthquake related shockwaves move the earthquake occurrence time forward or backwards in time from when it would have ordinarily occurred then the occurrence time shifts to another distinct time window which is controlled by the gravity related forces. Shift like that do not involve arbitrary amounts of time.

A - 10   Sun and moon gravity related forces and phenomena can reach high and low strength points in a number of different time cycles. Several of the most important cycles are multiples of roughly 6.25 hours (including 24 hours and 50 minutes) and multiples of about 7 days (including 28 days).

A - 11   A number of cycles having the same time length can be active at a given time. For example, one 6 hour cycle may be controlled by the position of the moon in the sky while at the same time another 6 hour cycle is controlled by the location of ocean tide crests and troughs.

A - 12   Which time cycles control occurrence times for earthquakes in a given fault zone system is strongly related to how forces and phenomena related to the cycles cause different rock layer sections in the fault zone system to be bent, stretched, or compressed in relation to one another. Another way of saying that is to state that what might be described as time cycle controlled Pressure Waves within the Earth's crust are actually triggering most earthquakes.

A - 13   It has been proposed by at least one researcher that in most cases earthquake occurrence times are linked with points of high and low strain in the time cycles while a smaller percentage of earthquake occurrence times are linked with the times when time cycle related earthquake triggering forces are experiencing their most rapid change in strength.

A - 14   Earthquake aftershocks are often not random events as far as their occurrence times are concerned. Those times are frequently controlled in part by those sun and moon gravity related time cycles. In fact, the aftershocks are even more likely to be controlled by them than the original earthquake.


EARTHQUAKE PRECURSOR SIGNAL GENERATION

Only a few earthquake precursor signal generation related processes are being discussed here.

B - 1   Phenomena such as the piezoelectric effect cause different types of electromagnetic energy field fluctuations to be generated in fault zones in the days, weeks, months, and years before earthquakes occur.

B - 2   Geomagnetic storm energy (which often has its origins in solar storms) can interact with fault zones and in the process contribute to or amplify certain fault zone activity related electromagnetic energy field fluctuations.

B - 3   Solar storm activity is associated with events taking place in and around the sun.  Interactions between solar storms and the Earth's geomagnetic field and earthquake fault zones are affected by the nature of the solar storm and also the orientation of the sun's magnetic field lines relative to those of the Earth.  About every 7 days the sun magnetic field lines orientation can shift from being in the same direction as those of the Earth to being in the opposite direction.

B - 4   Earthquake fault zone activity related electromagnetic energy fields can be focused in specific directions by a variety of phenomena. As a result before an earthquake a sensor located a certain distance on one side of the fault zone might be detecting strong energy field fluctuations while another one an equal distance away on the fault zone's other side might not detect anything. This is one of the important reasons that some energy field detectors give the appearance of not producing consistent results.

B - 5   The geomagnetic storm - earthquake fault zone interactions take place most often during the days and weeks just before an earthquake occurs when the fault zone’s rock layer's physical, chemical, and electrical properties change as it gets close to fracturing. But they can also take place months and perhaps years before the earthquake occurs.

B - 6   Phenomena such as lighting strikes and I believe geomagnetic storms can cause low frequency radio waves to be generated. I expect that events taking place within fault zones can also cause them to be generated.

       Those radio waves can travel around the world with relatively little loss in signal strength. Air is a reasonably good conductor for them. Water is better. And metal is best.

       As they interact with fault zones and/or when they are generated as the result of fault zone activities, information regarding fault zone events can be superimposed on them and then be carried by them around the world.

B - 7   Some fault zone event related electromagnetic energy field fluctuations range in time from 0.25 seconds to perhaps 30 seconds in duration. Others may persist for hours, days, and even weeks.

B - 8   Because the times when those energy field fluctuations occur can be controlled by forces or phenomena related to the gravitational pulls of the sun and the moon and also forces or phenomena related to geomagnetic storms etc., the fluctuations can occur at high strain level times etc. in the same time cycles which control earthquake occurrence times. Or they can occur at somewhat random looking times which may be days, weeks, or months apart.

B - 9   Other earthquake precursors possibly related to electromagnetic energy field fluctuations have been reported over the years. They would include Earthquake Lights which are a type of temporary, local luminescence, glowing balls of light which may be suspended in the air, and static noise on radio, television, and telephone transmissions.

B - 10   Interactions taking place between earthquake fault zones located even great distances apart can result in signals being generated which are stronger than they might otherwise be. For a hypothetical example, an approaching 8.0 magnitude earthquake in one fault zone might cause a signal associated with a 4.0 magnitude earthquake which is about to occur in another, distant fault zone to be amplified to a point where it appears to be for a 6.0 magnitude earthquake.


EARTHQUAKE SENSITIVE HUMANS, ANIMALS, AND PLANTS

C - 1   Electromagnetic energy field fluctuations linked with fault zone events and in some cases also with low frequency radio waves which are affected by or generated as a result of fault zone events, can produce strong reactions in humans, animals, and perhaps some plants.

C - 2   Various groups of organisms such as dogs and catfish will have different types of and different levels of sensitivity to those energy field fluctuations. And within a single group of organisms such as dogs different members will experience different types of and different levels of sensitivity to those energy field fluctuations.

C - 3   The responses of different organisms to the energy field fluctuations can result from basic energy field - biological system interactions. Those types of interactions would be similar to having water heated by the energy field of a microwave oven.

C - 4   They can result as side effects of energy field - biological system interactions which the evolutionary process has deliberately built into the members of a specific group of organisms for some reason such as to provide them with an ability to detect light or to sense the orientation of the Earth's magnetic field.

C - 5   They can result from energy field - biological system interactions which result from a deliberate effort by the evolutionary process to provide the members of a specific group of organisms with an ability to detect the approach of an earthquake.

C - 6   I have never heard about any research being done on plant sensitivity to those energy field fluctuations and am simply assuming that such effects probably exist.

C - 7   Earthquake sensitive animals may become temporarily agitated and/or disoriented by those energy field fluctuations. Some pet animals may wander away from home and become lost.

C - 8   Animals respond mostly to energy field fluctuations which are linked with events taking place in local earthquake fault zones. However, an occasional animal will be able to sense those fluctuations when they are associated with fault zones which are hundreds of miles away.

C - 9   Some pet animals which might behave as though they are earthquake sensitive are actually responding to fault zone activity related behavioral changes etc. which they are observing in their earthquake sensitive human owners.

C - 10   All humans are probably earthquake sensitive to some degree. Some fraction of the population is sufficiently sensitive to be able to clearly tell when a local earthquake is going to occur. A small fraction are sufficiently sensitive to be able to detect signals associated with earthquakes occurring around the world.

C - 11   Earthquake sensitive humans may experience physical, random or monthly type hormonal, emotional, and brain activity related responses to those energy field fluctuations. Headaches are among the most common responses. Nausea and dizziness can occur. Some people may experience abrupt, strong, and temporary or even somewhat longer lasting mood swings. High or low frequency tones may be heard in one ear or the other or in both ears at the same time. Short duration muscle jolts are common. Physical pain or some other type of sensation in a hand or foot and in extreme cases actual physical injury can occur. Dreams may be influenced by those energy field fluctuations. Especially sensitive people may begin hearing voices or seeing ghostly images. I believe that some people might begin seeing auras of different colors surrounding other people.

C - 12   Some earthquake sensitives can determine from which hand or foot is being affected by those energy field fluctuations etc. roughly where on the planet an approaching earthquake is likely to occur.

C - 13   Because water is a good conductor for the low frequency radio waves, earthquake sensitive humans will sometimes experience stronger reactions to the fault zone activity related energy field fluctuations when they stand near a large body of water such as on the ocean shore. Their reactions may get abruptly stronger when they stand near metal water pipes in a house and even when they begin running tap water into the kitchen sink.

C - 14   People who have water dowsing abilities are probably often detecting those energy field fluctuations as the signals travel through underground sources of water. I understand that dowsers are not usually able to detect the presence of electrically isolated bodies of water such as water present in a plastic swimming pool.

C - 15   Diet modification and dietary supplements can at times be used to control painful headaches or other undesired responses to those fault zone activity related electromagnetic energy field fluctuations. Common pain control medications such as aspirin may help. People who are strongly earthquake sensitive may need to rely on prescription drugs. Certain types of medications and foods can temporarily increase a person's sensitivity to the energy field fluctuations.

       Biofeedback can reportedly help reduce earthquake sensitivity. And I would expect that hypnosis might also work for some people.

       Some earthquake sensitives have reportedly been able to reduce the number and severity of their headaches or mood shifts etc. by moving to some other part of the country where fault zone activity related electromagnetic energy field fluctuations are weaker or at least different. After such a move it might require several weeks for the undesired energy field related responses to be reduced in strength.

C - 16   An inexpensive and effective way to detect approaching earthquakes might involve having a number of people who are earthquake sensitive contact some central facility when they were having earthquake warning experiences such as headaches. If 100 in some area were reporting in and only few in surrounding areas were doing the same then that might be a sign that an earthquake was about to occur near that first area. Additionally, earthquake forecasters and researchers could always ask them for more information when desired. That is obviously not possible with earthquake sensitive animals or conventional radio signal receivers or seismic activity monitors.  See:  A Demonstration Earthquake Prediction Program


EARTHQUAKE TRIGGERING FORCES

D - 1   An External Source Of Earthquake Triggering Pressure On The Earth’s Crust

       As the following drawing shows, the Gravity Point is the location on the surface of the Earth where the combined gravitational pulls of the sun and the moon are strongest.

Sun - Earth - Moon Barycenter

       As the Earth rotates each day on its North and South Pole axis, its crust on the side which is rotating towards the Gravity Point is accelerated by the sun and moon gravity while the side moving away from the Gravity Point is decelerated or slowed. That creates strain within the crust. And earthquakes are triggered at locations where that strain causes the crust to be bent, stretched, or compressed in just the right directions.

D - 2   An Internal Source Of Earthquake Triggering Pressure On The Earth’s Crust

       As that previous drawing shows, once each month the Earth and the moon rotate around a sun – Earth – moon barycenter or center of mass which is a certain distance down beneath the surface of the Earth. The location of that barycenter is constantly moving towards and away from the center of the Earth depending upon the location of the sun relative to the moon. And it is constantly moving around the location of a line drawn between the center of the Earth and the center of the moon.

       The solidified iron believed to be the major component of the Earth’s core is roughly 5 times as heavy or dense as the material in the Earth’s crust. As both the core and the crust rotate around that barycenter once each month the mass of the Earth in general attempts to assume the shape of an egg with its ends pointing away from the barycenter. And the solid core tries to move in a direction almost directly away from both that barycenter and from the location of the moon.

       The Earth’s gravity, rotation, and its somewhat rigid structure prevent it from undergoing rapid, appreciable changes in shape. And the core is kept from moving within the Earth by those factors. That results in the solid core exerting earthquake triggering pressure at different locations or pressure points in the crust from inside the Earth. Where those pressure points are located is determined by factors such as the location of center of mass of the core relative to the sun – Earth – moon rotation barycenter and the rotational orientation of the Earth as it revolves once each day on its North and South Pole axis.


REFERENCE INFORMATION

How Earthquake Sensitives Predict Earthquakes, GeoMonitor: Earthquake Prediction News; April, 1996

Breathing of the seafloor: Tidal correlations of seismicity at Axial volcano, Geology; June 2002; v. 30; no. 6; p. 503-506 (See especially Figure 4A.) (Also, that paper contains a helpful list of references.)

'Restless Earth' May Give Advance Notice Of Large Earthquakes (NASA News)
 

Prediction of Destructive Earthquake by the MDCB-Method, Wang, Wen-xiang and Yang, Wu-yang, Xian Sub-Institute, Coal Science Research Institute of China  MDCB Web Site

Cellphone can warn of earthquake

The 7.4 magnitude devastating earthquake in Pakistan on October 8 has raised several questions and posed new problems for scientists and engineers.

It is for the first time in seismic history that military personnel inside bunkers have been killed.

There have been cases before of soldiers killed in earthquake while overground. But the present case is unique and military engineers will now have to revise their design to make bunkers safe against seismic forces.

But the most frequently asked question is whether it is possible to predict the occurrence of an earthquake.

The answer to this question is, 'yes, to a sufficient degree of accuracy, about a few hours before the occurrence of an earthquake.'

Complete coverage: Tremors across borders

The earthquake precursors discussed here are exclusively for moderate to large earthquakes, meaning earthquakes of magnitude 6.0 or more on the Richter Scale. Earthquakes of lesser magnitude may or may not necessarily exhibit all the precursors.

Earthquake precursors have been divided into three broad categories. a. long term precursors. This means more than five years in advance; b. medium term precursors. This means within about four to six months in advance; and c. short term precursors. This refers to signs within 20 to 30 hours.

Long term seismic precursors: The location and magnitude of the likely epicentre is calculated based primarily on statistical studies, and the prediction is probabilistic. It is reliable when the probability of occurrence in the next five years or less is more than 0.6. However, such predictions are more useful to academic professionals and researchers.

Medium term precursors: Such predictions are able to give the location of the quake to a sufficient degree of accuracy, and the likely time of occurrence and magnitude are somewhat grotesque. The period of six months or so could be used for creating seismic awareness, taking suitable mitigatory measures, imparting seismic training to various forces etc. These are indicative of changes in the geomagnetic field, changes in the gravity field, rise in sub-surface temperature, rise in ground radioactivity etc. All these measurements indicate that something is wrong inside the earth. But it is not possible to arrive at any definite inference. Inputs from such parameters are useful to administrators and academicians.

Short term precursors: These could be seen by the common man and are highly reliable. People in seismically active areas may keep a copy of this article for further use. There are several types of short term precursors which could be seen and observed by different strata of society.

Thermal precursor: A few months before the occurrence of an earthquake the average temperature of the area keeps increasing. Weather report bulletins refer to temperatures above or below average by so many degrees. It is seen that in case the area is heading for an earthquake, the average temperature goes on increasing. On the day of the earthquake it is about 5 to 9 degrees Celsius above the average normal temperature for that day. This could be monitored by the meteorology department and also by thermometers inside homes.

Water precursor: There is a sudden rise or fall in water level in the wells. It could be as high as one metre. Sometimes the well water may turn muddy. At times a fountain appears inside the well. All these changes happen about one or three days before the earthquake. Sometimes a fountain appears in the ground. This normally happens a few hours before the quake.

The Gujarat earthquake of 2001

Another form of water precursor is seen in the flow of riverwater. There is a sudden and rapid increase or decrease of water flow in the river or nallah. This happens about one to two days before the quake.

Seismo-electromagnetic precursor: This is a very reliable precursor. It occurs and is exhibited about 10 to 20 hours before the quake. Before the occurrence of an earthquake the subsurface temperature rises. As a result of this the geomagnetic field is reduced. The reduction in geomagnetic field adversely affects the propagation of electromagnetic waves. This is experienced abundantly on the radio, television and telephone.

If one is tuning the radio at 1000 kHz, then the same station will be received in the potential epicentral area, about 10 to 30 hours before the quake, at higher frequencies. Maybe 1100, 1200, 1300, 1800, 1900, 2000 kHz or so.

Similarly, reception on television is affected by audio, visual and spectral disturbances. The telephone's reception is also adversely affected and there is a continuous disturbance that you can hear.

Apart from this, a recent reliable precursor is the mobile telephone. It has been found that about 100 to 150 minutes before the occurrence of an earthquake, mobile telephones start non-functioning or malfunctioning. The timespan indicated is sufficient to take all necessary mitigatory measures.

A word of caution: all the above precursors are valid only when seen and manifested extensively. Failure of one or two instruments should not be taken as a seismic precursor.

Animal precursor: It is seen that 10 to 20 hours before the occurrence of an earthquake, the entire animal kingdom becomes highly disturbed and restless. They move in a directionless manner and in fear. Birds do not sit on trees but move about at a low height, emitting a shrill noise. Rodents like rats, mongooses etc are in a panic. Domestic animals like cows, dogs, cats etc struggle against being tied up, and even turn on the owner.

Human precursor: Doctors and nurses observe that some sensitive patients in hospitals become highly disturbed. They exhibit a sudden rise in blood pressure, heart trouble, headache, migraine, respiratory disorders etc. Further, these psychosomatic signs are manifested without any provocation. It is also seen that the number of patients in the out-patient department increases by five to seven times, some 10 to 20 hours before the quake.

The best indicator is the number of child deliveries in any hospital. On the penultimate day of the earthquake the number of deliveries goes up about three to five times, while on the day of the earthquake it is as high as seven to eight times the normal.

Arun Bapat is a research seismologist. He was formerly Head, Earthquake Engineering Research Division at the Central Water and Power Research Station, Pune. He can be contacted at: arun_bapat@vsnl.com

Just prior to an earthquake, the birds go to the bottom of the cages and huddle in a corner.

Air-borne ions have reportedly a pronounced biological impact on humans
and animals (Charry, 1984; Krueger & Reed, 1976). Exposure to negative ions
is generally perceived by humans as pleasant, may enhance the cognitive
performance, and has overall beneficial health effects (Baron, 1987). By
contrast, positive ions elicit adverse reactions including headaches, nausea,
depression, and overall ill feelings. Adiabatically compressed, katabatic winds,
blowing downslope, are known in the Alps as Foehn, in the Rocky Mountains as
Chinook, and in Southern California as Santa Ana (Byers, 1974). Such winds are
laden with positive ions and are reputably the source of increased levels of
tension and irritability in the general population.
Combining these known adverse physiological effects of positive ions on
humans with the build-up of positive surface charges at the Earth’s surface, due
to the arrival of p-hole charge clouds, suggests that animals may react to the
positive ions that are most likely emitted from the ground. This would provide
a more consistent, and probably testable, explanation for the reported unusual
animal behavior prior to large earthquakes than the alleged sensitivity of animals
64 F. T. Freund
to minute vibrations or very low-intensity EM emission, which seem not to be
supported by other evidence.

 

Acoustic emission and earthquake precursors
G. Paparo(1), G.P. Gregori (1), I. Marson (2), M. Poscolieri (1)
(1) Ist. di Acustica CNR, Roma, Italy, (2) INOGS, Trieste, Italy,
gabriele.paparo@idac.rm.cnr.it

It makes nonsense investigating floods with no adequate monitoring and investigation
of the hydrologic regime of the river. It makes nonsense investigating seismic
and other extreme geodynamic occurrences, with no adequate monitoring and investigation
of the crustal stress propagation.

 We show how such rationale can be effectively implemented and practically exploited
by means of acoustic emissions (AE; ultrasound range). We find observational information
about earthquakes with several months advance and occurring within some wide area.
We stress, however, that we perform crustal diagnosis, rather than earthquake prediction.
We deal with different tectonic settings, i.e. the Italian peninsula, the Cephallonia Island (Greece),
and South America (Argentina, Provinces of San Juan and Mendoza). Some regular patterns and
time variations of AE are observed, in addition to the aforementioned phenomena associated
with a forthcoming seismic event. The potential ought to be stressed upon
investigating crustal stress on a continental or planetary scale, even by means of a
sparse set of simultaneously operated AE recording sites.

Earthquake precursors in magnetically conjugated ionosphere regions

Yu. Ya. Ruzhin, V. I. Larkina and A. Kh. Depueva

IZMIRAN, Troitsk-town, Moscow Region, 142092, Russia

Available online 19 June 1998.

Abstract

It is known that the highly conducting geomagnetic field flux tube and quasitrapped energetic particles provide strong electrodynamic coupling between magnetically conjugated regions in the ionosphere. Such a connection is likely to cause ionospheric disturbances before a pending earthquake not only over the epicentral zone but also in the opposite hemisphere. Data from the INTERCOSMOS-18 and ALOUETTE satellites were analyzed from this point of view. Earthquake precursors in magnetically-conjugated ionosphere regions were found in very low frequency (VLF) emission and F2-peak parameters. F2-precursors appear some days before the earthquake, and manifest themselves as an Appleton-type anomaly if the epicenter of the future earthquake is situated near the magnetic equator. Estimation of the electric field magnitude necessary to generate the observed anomaly was made. It was shown that an electric field of less than one mV/m must be generated in the ionosphere.

VLF precursors appearance some hours before the earthquake is localized close to the magnetic shell corresponding to the future earthquake epicenter and have a belt-like structure (longitudealigned for more than some tens thousands kilometers) in both hemispheres. VLF precursors are followed by energetic particle (electrons energy W>40 keV) precipitation.

 

Detecting Earthquake Precursor:
A Thermal Remote Sensing Approach


Arun K. Saraf*, Swapnamita Choudhury1, Vineeta Rawat*, Priyanka Banerjee2, Sudipta Dasgupta3 and J.D. Das**
*Department of Earth Sciences
**Department of Earthquake Engineering
Indian Institute of Technology Roorkee
Roorkee – 247667, India
saraffes@iitr.ernet.in

 
Abstract

Seismologists have not yet found a reliable earthquake precursor. In order to look for a reliable earthquake precursor a remote sensing based thermal technique has been employed recently based on the concept that stress accumulated in rocks in tectonically active regions may be manifested as temperature variation through a process of energy transformation. Rise in land surface temperature (LST) before an impending earthquake has been detected for 23 earthquakes using pre- and post-earthquake datasets of thermal sensors such as NOAA-AVHRR, Terra/Aqua-MODIS and SSM/I for the different parts of India and the world. Analysis revealed a transient short-term thermal rise in LST ranging from 2-12° C around epicentral areas of different earthquakes. In the studied earthquakes, pre-earthquake thermal anomalies generally started developing about 1-10 days prior to the main event depending upon the magnitude and focal depth, and then disappeared just before the main shock. This temperature increment prior to an impending earthquake can also be attributed to degassing from rocks under stress and/or to p-hole activation in stressed rock volume and their further recombination at the rock-air interface. A precise correlation of LST maps of Bam (26 Dec 2003, Iran) and Zarand (22 Feb 2005, Iran) with InSAR generated deformation maps; published by other workers (Stramondo et al., 2005a, 2005b, Parsons, 2005); also provides evidence that the thermal anomaly is ground related phenomena, not an atmospheric one.

1.0 Introduction
Satellite thermal remote sensing has emerged as a potential tool in detection of pre-earthquake thermal infrared (TIR) anomaly in land surface temperature in and around epicentral regions. Thermal data in seismic studies was first put to application in Russia in 1985 (Tronin, 2000) and results were published (Gorny et al., 1988).This paper presents observations made by post-earthquake analysis of 23 past earthquakes occurred in different parts of the world and using thermal data from different satellites (Choudhury, 2005, Saraf and Choudhury, 2005a, 2005b, 2005c and Panda et al. 2007). Bam earthquake, Iran (26 Dec 2003) and Bhuj earthquake, India (26 Jan 2001) are two classic examples and are discussed here in detail. The Bam earthquake (table 1) brought devastation over the SE Iran, killing about 31,000 people, injuring 75,600 and rendering 75,600 homeless (Choudhury et al., 2006). In case of Bhuj earthquake (table 1) the death toll was estimated about 20,083 according to Gujarat Government figures and was accompanied by the wide scale damage to property and economy (Saraf and Choudhury, 2005a).

2.0 Concept of Thermal Infrared Anomaly
Before an imminent earthquake earth’s crust passes through an earthquake preparatory phase. Accumulation of stresses and resultant pressure built up leads to the rise in LST. The enhanced TIR emission from earth’s surface retrieved by satellites prior to earthquakes is known as “thermal anomaly” (Freund et al., 2005).

 
Table1: Details of the studied earthquakes.

 
S.N. Earthquake Origin Location Mag. (USGS) Mw Focal Depth (km)
Date Time (UTC) Lat. (N) Long. (E)
1. Persian Gulf, Iran 28 Jun 06 21:02 26.82° 55.90° 5.8 10
2. Darb-e-Astaneh, Iran 31 Mar 06 01:17 33.58° 48.79 6.1 07
3. Fin, Iran 25 Mar 06 07:28 27.57° 55.69° 5.9 18
4. Faryab, Iran 28 Feb 06 07:31 28.12° 56.87° 6.0 18
5. Qeshm, Iran 27 Nov 05 10:22 26.77° 55.86° 6.0 10
6. POK, India 08 Oct 05 03:51 34.43° 73.53 7.6 10
7. Zarand, Iran 22 Feb 05 02.25 30.75° 56.82° 6.4 14
8. Banda-Aceh, Sumatra 26 Dec 04 00:58 3.31° 95.95 9.3 30
9. Firozabad-Kajoor, Iran 28 May 04 12:38 36.29° 51.59° 6.3 28
10. Bam, Iran* 25 Dec 03 01:56 29.00° 58.34° 6.6 10
11. Kerman, Iran 21 Aug 03 04:02 29.05° 59.77° 5.9 20
12. Jahron, Iran 10 Jul 03 17:40 28.35° 54.17° 5.8 10
13. Boumerdes, Algeria 21 May 03 18:44 36.90° 03.71 6.8 10
14. S. Xinxiang, China 24 Feb 03 02:03 39.61° 77.24° 6.4 11
15. Changureh-Avaj, Iran 22 Jun 02 02:58 35.63° 49.05° 6.5 10
16. Hindukush, Afghanistan 25 Mar 02 14:56 35.97° 69.18° 6.1 33
17. Double Earthquakes of Hindukush, Afghanistan 03 Mar 02 12:08 36.44° 70.45° 6.2 195
18. 03 Mar 02 12.08 36.54° 70.42° 7.4 256
19. Bhuj, India* 26 Jan 01 03:16 23.23° 70.18° 7.9 16
20. Izmit, Turkey 17 Aug 99 00:01 40.74° 29.86° 7.6 17
21. Zhangbei, China 10 Jan 98 03:50 41.08° 114.50 6.2 30
22. Jabalpur, India 22 May 97 04:21 23.08° 80.06° 6.0 35
23. Kalat, Pakistan 04 Mar 90 19:46 28.92 ° 66.33° 6.1 10

*Earthquakes here discussed in detail.


There may be various physical explanations for thermal anomalies appearing before an impending earthquake. However, authors prefer to explain the physical mechanism of thermal anomaly mainly by two leading theories: Earth-degassing Theory (Zu-ji et al., 1999) and p-hole activation Theory (Freund 2000, 2002 and 2003). During the preparatory phase of an earthquake, in high tectonic stress, pore spaces in the rocks are reduced due to subsurface pressure releasing the trapped gases. These gases on reaching earth’s surface are already at an elevated temperature with respect to the air temperature due to subsurface geothermal heat and secondly, it also creates a local greenhouse effect on the land surface thus serving as the source of outgoing anomalous radiation. Positive thermal anomalies at a regional scale were observed with the examination of around 9000 thermal images for the Middle Asian earthquake. These anomalies were attributed to the green house effect that was caused before the earthquake by an increase in gases like CO2, CH4 etc (Tronin, 1996). A new theory of charge generation in rocks prior to earthquakes is given by Freund (2000, 2002 and 2003). This theory keeps parity with laboratory experiments (Qiang et al. 1997, Ozounov and Freund 2004) and also provides explanation for other observed geophysical precursors. Electronic charge carriers can be free electrons or sites of electron-deficiency in the crystal structure, the latter known as p-holes (Freund, 2000). These p-holes in their dormant state are in form of Positive Hole Pairs (PHPs). A PHP is devoid of any charge, as outer shells of both of the two O- atoms are filled.

Introduction of PHPs in minerals during rock-genesis and alteration has been explained explicitly by Freund (2002). Interestingly, O- -O- bond-distance (1.5A°) is almost half of O2–-O2- bond distance (2.8–3.0A°). It implies that the peroxy-bound O- has a small partial molar volume, thus having a tendency to be favored by high pressure. Being dissociated under high-accumulated stress a PHP introduces two holes (charge deficiencies) into the valence bond to attain electrostatic stability, causing the insulator to become a p-type semiconductor. Positive holes propagate though an oxide or silicate structure by electron hopping and immediately rush towards the Earth’s surface with high speed. Immediately after reaching a non-solid medium like the atmosphere or water, a positive hole acquires an electron to become O2- from O- . The recombination leads to vibrationally highly excited O-O bonds. These O-O bonds can de-excite either by emitting specific mid-IR photons in the 11-12 µm region or by channeling the excess energy into neighbouring Si-O and Al-O bonds, which in turn will emit in the 8-10 µm region (Freund et al., 2005, 2007). The observation and the spectral signature of the emitted radiation provides strong evidence that the underlying effect is a kind of mid IR luminescence arising from the recombination of p-holes at the rock surface. The energy therefore emitted by this electron acquirement increases the surface temperature. Several Remote Sensing Rock Mechanics (RSRM) (Wu et al, 2000) experiments and other studies (Pellet et al., 2007) conducted to predict material behavior on stress application and time-space forecast of rock failure have validated the infrared emission from the rocks under stress. However, it should be noticed that this heat is starkly different from the frictional heat that develops at the fault surfaces during rupture. This frictional heat, that takes a lot more time to come up to the surface, actually develops at the time of the earthquake itself, and hence has no contribution towards the pre-earthquake thermal anomaly (Banerjee, 2007).

3.0 Data and Methodology
Day and nighttime High Resolution Picture Transmission (HRPT), Global Area Coverage (GAC), Local Area coverage (LAC); AVHRR datasets available from Indian Institute of Technology Roorkee-Satellite Earth Station, India (IITR-SES) and from National Environmental Satellite Data and Information Service (NESDIS) website; were mainly used in this study. Data from other satellites thermal sensors like Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra/Aqua satellites, Special Sensor Microwave Imager (SSM/I) on board DMSP satellites (SSM/I-DMSP) were also used. NOAA-HRPT has a spatial resolution of 1.1 km and a temperature resolution of 0.5° C. LST time series maps were prepared for day and night data sets.

Data calibration and temperature calculation is based on the Becker and Li (1990) split window algorithm, which uses the differential absorption effect in channel 4 and 5 of NOAA-AVHRR. A user specified temperature range for daytime as well as nighttime scenes was specified in each case. Temperature outside this range was masked. Cloud covers were delineated and avoided for any temperature calculation. Image co-registration and correction for different satellite view angles were not done since any image overlay analysis was not intended. Wherever, available air temperature data for previous and subsequent years along with the year of earthquake occurrence was collected and Temperature Variation Curves (TVCs) and Normal Curves were generated. A comparison has shown that the difference of LST and air temperature is within +/- 1° C. In the present study air temperature data has been used to support LST analysis. This means that just based on air temperature data such study of pre-earthquake thermal anomaly would not be possible logically from spatial point of view (Choudhury et al., 2006).

4.0 Observations
Both nighttime and daytime datasets were utilized for the study and an almost consistent pattern in the temporal and spatial development of thermal IR anomaly was observed. The two case studies are here discussed in detail.

4.1 Bam earthquake, Iran (26 Dec 2003, Mw 6.6)
Analyses of the daytime and nighttime datasets for the Bam showed a distinct anomaly in LST which appeared before the main shock (table 2 and figure 1). The temperature increase was about 5-7° C than the usual temperature of the region. At some places, the temperature was about 6-10° C higher than the normal temperature of the region in that period of the year.

In the nighttime maps, it was seen that on 18 Dec 2003 there was a complete normal temperature regime in the region. The appearance of an intense thermal anomaly was seen around the earthquake epicenter near Bam on 21 Dec 2003 (data was not available on 19 and 20 Dec 2003 for the same time of acquisition).

If this rise had started earlier and then reached its peak on 21 Dec, it could not be ascertained through the gap in the data. To bring continuity to this gap, data acquired on 18 and 19 Dec 2003 night at 16:30-18:15 UTC hours, with a difference of four hours, was analyzed. A thermally anomalous region had developed near the epicenter on 19 Dec 2003, which rose to the maximum amplitude of thermal anomaly on 21 Dec 2003. The temperature was around 7-13° C on 21 Dec 2003 (around 7-8° C higher than normal temperature in the region around that period of the year). On 22 Dec 2003, the temperature was less than this boost. The temperature was normal again on 23 Dec 2003. The earthquake came on 26 Dec 2003, five days after the boost in LST as observed through nighttime data.

Daytime LST time series maps show that the rise in temperature started on 22 Dec 2003. The anomaly stayed on till 24 Dec 2003 (just two days before the earthquake). The normal temperature was around 22-25° C on 21 Dec 2003. On 24 Dec 2003, the temperature was around 29-32° C (about 7-10° C higher than usual temperature of the region around that period of the year).

Analysis and similar processing of nighttime NOAA-AVHRR data of the year 2004 acquired at around the same time and on the same days as the 2003 data, showed that there was no such abnormal behavior of the LST on those days in that year. The anomalous region in the nighttime data on 21 Dec 2003 occupied an area of about 308,000 km2 and in the daytime LST map of 24 Dec 2003 covered an area of about 328,200 km2. The InSAR co-seismic surface deformation studied by Stramondo et al. (2005a, 2005b) of the same Bam earthquake, revealed the existence of an unmapped fault 4 km west of Bam fault, the only active tectonic feature in the epicentral area (figure 2). Analysis of the hypocenter distribution of aftershocks recorded by a temporal seismic network installed 6 weeks after the main shock revealed source of earthquake as the “Arg-e-Bam” fault (Nakamura et al., 2004). The unwrapped interferometric phase onto the SAR amplitude image displays the surface displacement pattern, extends to around 30 km in N-S and to around 20 km in the E-W direction (figure 2b). It was observed that co-seismic InSAR surface deformation filed is in concurrence with the pre-seismic thermal anomalous area developed to the east of the epicenter on the 26 Dec 2003. This analysis indicates a plausible relationship between the two and suggests that LST pre-earthquake phenomena are land related.

4.2 Bhuj Earthquake, India (26 Jan 2001, Mw 7.9)
On 14 Jan 2001, S-W Gujarat, with respect to the surrounding region started to pick up an increase in temperature. Within a few days there was a spread and increase of the thermal anomaly. The area experiencing this anomaly had spread in a NW-SE direction. On 23 Jan 2001 there developed an anomaly, which had an area of around 300 km2, SE of earthquake epicenter (table 3 and figure 3).

 

Figure 1: Daytime NOAA-AVHRR LST time series map of Iran before and after the earthquake in Bam, Iran on 26 Dec 2003 (Choudhury, 2005). Red colour filled star symbol indicates the day of the earthquake, whereas without colour filled star signifies location of the impending earthquake epicenter of 26 Dec 2003.


 

Figure 2: A comparison of thermal anomaly field of 26 Dec 2003 derived from nighttime NOAA-AVHRR image (a) and surface displacement field derived by unwrapping interferometric phase onto the SAR amplitude image (b). The black line in (b) indicates the SAR inferred earthquake fault location. Red star indicates the location of Bam earthquake epicenter.


Succeeding this boost, the increase started to wane out, after 23 Jan 2001, just three days before the earthquake of 26 Jan 2001. On 23 Jan, temperature was maximum between 28ºC - 31ºC in the anomalous region (about 5º C - 7º C higher than the normal temperature). On 28 Jan 2001, the anomaly disappeared and the region showed a normal temperature of around 24º C. In order to ascertain that the thermal anomaly discussed above was induced by the 26 Jan 2001 earthquake, same period satellite data from NOAA-AVHRR for the year 2003 of the same area were analyzed. Further, it was observed that during the January 2003 the thermal regime of the study area shows completely normal pattern.

Available weekly average temperature data (between 1951-1980) of epicentral area and surrounding region clearly shows that in the third week, the temperature trend touches the lowest temperature as compared to second and fourth weeks of the above years (figure 4). Whereas, in the year 2001, a peak was observed in the third week (instead in the second or fourth weeks). This is in contrary to the observed trend for the past thirty years (1951-1980). This ground observation is in accord with the observed satellite thermal detection during January 2001.

Table 2: Details of nighttime and daytime NOAA-AVHRR data (acquired by IITR-SES) of the year 2003 used to prepare LST time series maps as shown in figure 1 to study pre-earthquake thermal anomaly of 26 Dec 2003 Bam Earthquake.

 
S. N. Date Time of Acquisition (UTC)
    Nighttime Daytime
  21:30-23:00 hours 16:30-18:15 hours 09:00-10:00 hours
1 18 Dec 03 21:40 18:00 -
2 19 Dec 03 - 17:36 -
3 20 Dec 03 - 17:03 -
4 21 Dec 03 22:50 - -
5 22 Dec 03 22:38 - 10:01
6 23 Dec 03 10:25 - 09:46
7 24 Dec 03 22:14 - 09:34
8 25 Dec 03 22:02 - 09:23
9 26 Dec 03 21:50 - 09:11
10 27 Dec 03 21:39 - -


5.0 Discussion and Conclusions
The post-event investigation of pre-earthquake thermal anomalies by analyzing pre- and post-earthquake LST images reveals valuable information about the changes in the thermal IR regime of the affected area. The understanding of the buildup of pre-earthquake TIR anomalies and their satellite detection provides possibilities of a reliable and potential precursor.

On the basis of this study we can conclude that earthquakes with moderate to strong magnitudes may be preceded by the detectable TIR anomaly. So far studied 23 recent past earthquakes have showed temporal, short-term TIR anomalies before the event. The analyses of time series LST maps for these earthquakes showed a 2-13°C rise in LST 1-10 days before the earthquake struck (table 4). It was also noticed that magnitude and focal depth play a vital role in intensity and spatial extent of the thermal anomaly. Higher earthquake magnitude and shallower focal depth are favorable conditions for the appearance of intense thermal anomaly with larger spatial extent and vice versa.

 

Figure 3: Daytime NOAA-AVHRR LST time series map of Gujarat, India before and after the earthquake in Bhuj, India on 26 Jan 2001 (Saraf and Choudhury, 2005a).


 

Figure 4: Weekly average temperature variations in the second, third and fourth weeks (Julian weeks) of the years from 1951 to 1980 (data source: IMD 1991) (Saraf and Choudhury, 2005a).

 
A difference in the time of appearance of the thermal IR anomalies has been observed in daytime and nighttime data of the Bam, Bhuj and other studied earthquakes. This is probably due to the typical meteorological phenomena or the appearance of clouds during the day or night.

Table 3: Details of nighttime and daytime NOAA-AVHRR data (acquired by IITR-SES) of the year 2001 used to prepare LST time series maps as shown in figure 3 to study pre-earthquake thermal anomaly of 26 Jan 2001 Bhuj Earthquake.

 
S. No. Date Time of Acquisition (IST)
Daytime (Hours) Nighttime (Hours)
Scene 1 12 Jan 01 17:39 23:01
Scene 2 14 Jan 01 17:15 22:15
Scene 3 21 Jan 01 17:32 22:59
Scene 4 23 Jan 01 17:09 22:13
Scene 5 27 Jan 01 18:02 22:23
Scene 6 28 Jan 01 17:50 22:00
Scene 7 29 Jan 01 17:38 23:20


In case of cloud-free weather conditions this difference may be attributed to differential heating contrary to differential cooling. As such, the diurnal thermal regime may not be conformable at daytime or nighttime. But it was a regular observation that the anomalies appear first in nighttime thermal images. A prominent observation regarding the earthquakes of moderate magnitude, i.e. less than magnitude 6, is the appearance of a dual TIR peak in surface temperatures instead of single rise observed previously. In case of Kerman earthquake, Iran (21 Aug 2003) the first peak appears about ten days before the earthquake and second temperature peak relatively closer to the main shock i.e. 2-6 days (table 4). This may lead us to understand that perhaps the energy accumulated in the stressed rocks have been released sporadically in form of electromagnetic emission, apparent temperature increment or any other geophysical earthquake precursor, which in turn might reduce the magnitude of the main shock. Study of other earthquakes with prominent aftershock events reveals that the occurrence of aftershocks prevents the re-establishment of normal conditions even after the main event is over.

Due to prevailing residual stresses, the epicenter and adjoining areas experience aftershocks and thus the disappearance of the IR anomaly takes longer time. The observation of common field for TIR anomaly and surface displacement of Bam is a significant contribution. The concurrences of both the fields provide enough reason to believe plausible relationship between surface deformations and appearance of the thermal IR anomaly. The observation substantiates the explanations for pre-earthquake TIR anomalies given above. This is an important contribution to the knowledge of pre-earthquake thermal IR anomaly phenomena and for developing it as reliable and potential earthquake precursors.

Table 4: Characteristics of the short-term pre-earthquake thermal anomalies associated with the Iran earthquakes studied through NOAA-AVHRR and other thermal data sets.

 
S. N. Earthquake Mag. (Mw) Focal Depth (km) Pre-earthquake thermal anomaly (before the earthquake) Intensity of thermal rise Spatial Extent of Thermal Anomaly (km2)
Rise started Maximum rise observed
1. Persian Gulf, Iran 5.8 10 11 days 8 days 2-4 °C 770,457
2. Darb-e-Astaneh, Iran 6.1 07 7 days 4 days 5-10 °C 80,000
3. Fin, Iran 5.9 18 13 days 10 days& 2 days** 5-7 °C 963,072
4. Faryab, Iran 6.0 18 5 days 2 days 4-8 °C 729,344
5. Qeshm, Iran 6.0 10 7 days 3-4 days 2-3 °C 382,074
6. POK, India 7.6 10 8 days 6 days 6-8 °C 45,000
7. Zarand, Iran 6.4 14 5 days 1 day 10-12 °C 75,600
8. Banda-Aceh, Sumatra 9.3 30 >15 days 1day 6-12 °C -
9. Firozabad-Kajoor, Iran 6.3 28 5 days 1 day 4-6 °C 71,22,452
10. Bam, Iran* 6.6 10 4 days 2 days 7-10 °C 328,200
11. Kerman, Iran 5.9 20 11 days 10 days& 6 days** 5-10 °C 949,780
12. Jahron, Iran 5.8 10 9 days 7 days 5-7 °C 13,82,543
13. Boumerdes, Algeria 6.8 10 7 days Few hours 5-10 °C 91,100
14. S. Xinxiang, China 6.4 11 < one week 4-6 °C -
15. Changureh-Avaj, Iran 6.5 10 7 days 2 days 5-10 °C 163,243
16. Hindukush, Afghanistan 6.1 33 2 weeks 6-10°C -
17. Double Earthquakes of Hindukush, Afghanistan 6.2 195 Few days to a week 4-10 °C -
18. 7.4 256 -
19. Bhuj, India* 7.9 16 12 days 3 days 5-10 °C 179,150
20. Izmit, Turkey 7.6 17 1 week 6-10°C -
21. Zhangbei, China 6.2 30 3 weeks 4-8 °C -
22. Jabalpur, India 6 35 11 days 8 days 5-10 °C 154,072
23. Kalat, Pakistan 6.1 10 2 weeks 2-10 °C -


Acknowledgement:
We are greatly indebted to Department of Science and Technology (Seismological Divison)

(now Ministry of Earth Sciences), New Delhi for financial assistance, National Institute of Oceanography (NIO), Goa, and India Meteorological Department (IMD), Ahmedabad for providing valuable data for this study.

References
 
 

Earthquake precursor is reported

October 9, 1996

 

MIT scientists have discovered strong evidence that a 1994 earthquake in the mid-Atlantic off South America was initiated by a 100-second episode of slow slip along the fault, and that such precursors may be typical for large earthquakes located on the ocean-ridge transform system.

The discovery, which was reported in the October 4 issue of Science, gives clues about how some large earthquakes start and could shed light on the long-standing question of which types of earthquakes, if any, are predictable.

"What we learn about quakes in the oceanic ridge-transform system may shed light on earthquakes in general," said Professor Thomas H. Jordan, head of the Department of Earth, Atmospheric and Planetary Sciences. Professor Jordan, Jeffrey J. McGuire, a graduate student in EAPS, and Pierre F. Ihmle, a former MIT student now at the Institut de Physique du Globe in Paris, are co-authors of the study.

The scientists used a battery of sophisticated analytical techniques to look at data collected from a worldwide network of seismic stations in different ways. They specifically looked at low- frequency energy releases that seismologists typically don't study. The slow precursors are not easy to spot on seismometers because there is no obvious isolated burst of energy, but rather anomalies in the size and arrival times of low-frequency waves.

"Earthquakes in the middle of the ocean have a slow, smooth deformation that is difficult to detect. These earthquakes preferentially radiate energy at low frequencies," said Dr. Jordan, the Robert R. Shrock Professor of Earth and Planetary Sciences.

In 1993, Drs. Ihmle and Jordan and another former student, Paolo Harabaglia, said that the 1989 Macquarie Ridge earthquake, which occurred in an oceanic fault south of New Zealand, was initiated by a precursor that released energy slowly and smoothly beginning about six minutes before the main rupture. That large earthquake had a magnitude of 8.2.

The new paper in Science reports that the 1994 earthquake near the equator in the Romanche Transform in the mid-Atlantic ridge-transform system off South America was preceded by a similar episode of slow and smooth energy release that started 100 seconds before the main earthquake. That earthquake had a magnitude of 7.0. Earthquakes of this magnitude can cause major damage if they occur on land.

By exploring the spectral characteristics of the earthquake, or the energy of the seismic waves emitted from it at various frequencies, the scientists discovered that the slow precursor of the Romanche earthquake grew for at least 100 seconds before it triggered a fast, main rupture. The total amount of fault slippage that occurred during the precursor was about 15 percent of the main shock's slippage.

"Due to the low noise levels of the modern global network of seismic stations, this event marks the first time where a slow precursor's energy release can be seen as the first energy arriving at seismometers, providing more convincing evidence for the precursor's existance than in the past," Mr. McGuire said. He added that studying slow precursors may shed light on how big earthquakes get started.

This type of precursor does not occur in California, a highly active earthquake region. One reason may be that the fracture characteristics of continental rocks differ from those of the oceanic crust.

"We hope to learn something fundamental about earthquakes. You need to look at the full complexity of the earthquakes on the planet to understand earthquakes as a whole," Professor Jordan said.

To aid in further study of the precursor phenomenon, he and colleagues have started a global survey that uses very low-frequency waves to study earthquakes of magnitude 6.0 or greater that have occurred since 1995.

Funding for the research came from the National Science Foundation and NASA.
 

A version of this article appeared in MIT Tech Talk on October 9, 1996.

_________________________________________________________________________________________________________________________________________

Algae and earthquake precursors

Category: earthquakesgeohazardsgeology
Posted on: May 16, 2007 7:45 AM, by Chris Rowan

Concentrations of the natural pigment chlorophyll in coastal waters have been shown to rise prior to earthquakes.

These chlorophyll increases are due to blooms of plankton, which use the pigment to convert solar energy to chemical energy via photosynthesis.

This is based on an article by Singh et al. [doi] in Advances in Space Research, part of a special issue devoted to the use of satellite remote sensing for studying and predicting natural hazards such as earthquakes. The authors claim that you can detect a rise in sea surface temperature just before large coastal earthquakes. The blooms observed in this study are, they say, result from an increased flow of heat energy from the ocean to the atmosphere, enhancing the upwelling of cold, nutrient rich water and fuelling a boom in the growth of photosynthetic algae.

This is all very interesting, but what is unclear in the paper is the reason so much heat energy is being released prior to the earthquake, which is itself releasing accumulated strain energy. And, looking at the paper, I'm not really sure the relationship is quite as clear-cut as is claimed.
 

The figure below is from their paper, which shows the amount of upwelling measured at three localities in the month in which an earthquake was recorded. The peaks can be linked to increased concentrations of chlorophyll, as measured by satellite. The dark vertical lines mark the date of the earthquake.

algaeeqsb.jpg

As you can see, there are a number of peaks in each of these plots, not just the one associated with an earthquake. In all but perhaps the last case the peak preceding the earthquake doesn't seem to have a significantly different magnitude to all the others, and the timing of the earthquake seems quite variable - in one case it coincides with the peak, in another it is ten days afterwards.

You may beg to differ, but I'm rather unconvinced by this. It seems on a par with other proposed earthquake precursors like low frequency EM fluctuations (on which there are also some papers in this issue of Space Science Research) - not only is it not clear how they work, but the supposed precursor signals are not obviously distinguishable from background. Sadly, I don't think we'll be using satellite data to reliably forecast earthquakes just yet.

 

__________________________________________________________________________________________________________________________________________________________
In China

Through the Natural Hazards and Emergency Response Program, ESS (primarily the Geological Survey of Canada) is participating in the joint Investigation of Earthquake Precursors project with the China Earthquake Administration (CEA). This project is defined under the ESS-CEA Memorandum of Understanding for cooperation in the field of earthquake studies. As a part of this project, ESS and CEA scientists will:

ESS scientist (right) interviews earthquake survivors in Haicheng.

ESS scientist (right) interviews earthquake survivors in Haicheng.

The CEA has been running a multidisciplinary monitoring program to capture precursors of earthquakes for prediction purposes for over 40 years and has accumulated a wealth of knowledge. ESS can benefit by learning from the successes and failures of China's program, and will be able to apply new knowledge and techniques to its own efforts to understand how and when earthquakes recur in order to mitigate earthquake hazards in Canada.

 

 

EARTHQUAKE PRECURSORS’ DISPLAY FEATURES IN THE TERRITORY OF ARMENIA Dr. Hrachya M. PETROSYAN

Armenian National Survey for Seismic Protection (NSSP),
Head of Unified Center for Seismic Hazard Assessment,
375054, Davidashen-4, Yerevan, Armenia,
Phone: (374-10) 282811, 361562; Fax: (374-10) 366280,
E-mail: hrachikp@yahoo.com petrosyan@nssp-gov.am

Based on retrospective analysis of multidisciplinary observations in research
seismoactive zone for current seismic hazard assessment, practically all seismic events
which had occurred from 1983 till 2002 in the territory of Armenia and adjacent territories are
systematically tested. The precursory anomalies of 11 strong regional ( 6.0) and 7
perceptible local earthquakes ( 3.7) are classified.
Practice shows that small-but-perceptible local earthquakes of previous years
invariably caused some panic among the local population, and damage of buildings and
structures. For this reason, our study includes the local small-but-perceptible earthquakes
with 3.7< M< 5.0 as well.
Time-series providing many years of data on 15 different types of parameters from
more than 80 stations of the National observation network of Armenian NSSP are processed
(some stations having been in operation for more than 24 years). The medium-term, shortterm
and operative probably-seismogenic anomalies as well as co-seismic and post-seismic
effects anticipating the tested seismic events are distinguished. The vertical scales of timeseries
are defined so that 2 and more are included in 1 mm, where is an error of
measurements.

The classification is performed and the parameters of precursory anomalies are
determined. The graphic images of each anomaly are allocated and the probability of seismic
realization of anomalies P is determined as the relation m/n, where m is quantity of those
anomalies, after which the tested seismic events with ³3.7 occurred, and n is the total of
such anomalies at the given station and on given parameter available in the Armenian NSSP
database. The Catalogue of precursory anomalies of the tested earthquakes is completed. It
includes 167 anomalies-tests displaying the features of 656 anomalies. The precursory
anomalies were grouped by earthquake. The value of P, the epicentral distance of the
monitoring station, the moment of the seismic event and its magnitude were each shown.
The overwhelming majority of precursors was classified and represented in the Catalogue for
the first time. There are many examples, in particular, the imposing of precursory anomalies
of different order, as well as the presence of pre-, co-, post-seismic periods in observed
anomalies (fig.1). At the end of the Catalogue, the summary table is given (table 1), where
the data on parameters of the tested earthquakes and precursory anomalies are provided.
The tests were made within several years, with inclusion of new seismic events
occurred and in process of new data acquisition from observation stations. The Catalogue of
precursors was supplemented few times, and the anomalies included were critically
overestimated and sometimes rejected. Naturally, the Catalogue will be supplemented with
the tests of future strong regional and perceptible local earthquakes.
The Catalogue is in daily use at the Unified Center for Seismic Hazard Assessment of
Armenian NSSP for current seismic hazard assessment.
- 2 -
C, imp/min 1
C, ml/l 2
H,m 3
Fig.1. The pre-seismic anomaly of subsoil Radon concentration at the Noemberyan station
before the Baku (2000, M=6.5) earthquake (1); pre-, co-, and post-seismic anomalies of
Helium concentration at the Kajaran station, Spitak (1988, M=7.0) earthquake (2); pre-, co-,
and post-seismic anomalies of underground water level at the Noemberyan station,
Noemberyan (1997, M=4.4) earthquake (3).
The moments of seismic events are marked by vertical arrows, and anomalies - by horizontal
and inclined arrows.
M=6.5
D=445
M=7.0
D=257
- 3 -
The summary table of earthquakes’ test-precursors (the 1-st page)
Table 1
N
Earthquake
Parameter
Station
D
(km)
Duration
of timeseries
DT
(days)
Amplitude
of
anomaly
Type of
anomaly
Probability
P = m/n
1 2 3 4 5 6 7 8 9 10
1 Narman, Turkey
30.10.1983
M = 6.8
Geomagn.-
SD
He
He
He
Gyulag.- Jermuk
Ararat
Surenavan
Kajaran
262
264
267
397
since 80 /
80
since 07.78
since 12.79
since 01.83
510
236
180
142
22
195
360
1130
m/t
sh/t
sh/t
sh/t
1/1
2/4
1/1
1/4
2 Spitak, Armenia
07.12.1988
M = 7.0
Geomagn.-
SD
HGD
HGD
HGD
HGD
HGD
Cl
HCO3
pH
Cl
pH
He
HGD
HGD
Radon s/s
HGD
pH
Gyulag.- Jermuk
Metsamor
Goris
Sevan
Noemberyan
Dzorakhbyur
Surenavan
Kajaran
Kajaran
Kajaran
Ararat
Kajaran
Dzorakhbyur
Goris
Kumayri
Noemberyan
Surenavan
19
85
235
71
59
87
133
257
257
257
125
257
87
235
32
59
133
since 80 /
80
since 01.86
since 01.87
since 01.87
since 01.86
since 01.88
since12.79
since 01.83
since 01.83
since 01.83
since 07.78
since 01.83
since 01.88
since 01.87
since 01.86
since 01.86
since 12.79
1020
840
457
380
247
234
220
165 and
reaction
157 and
reaction
127 and
reaction
117
68 and
reaction
57
39
34
14
reaction
32
52
48
198
220
70
0.8
1.22
2.44
*
0.34
0.72
0.8
0.85
480
1280
13
9
18
247
0.43
m/t
m/t
m/t
sh/t
sh/t
sh/t
sh/t
sh/t
---
sh/t
---
sh/t
---
sh/t
sh/t
---
sh/t
sh/t
sh/t
sh/t
---
1/1
1/2
1/2
1/1
2/5
1/1
1/5
1/5
---
1/8
---
1/4
---
1/2
5/9
---
1/4
1/1
1/13
1/1
---
- 4 -
Legend:
HGD – hydrogeodynamics - underground water level (cm)
Hydrogeochemical regional network:
He – concentration of Helium in mineral water (ml/l ·10-6)
Cl - concentration of chlorine in mineral water (ml/l ·10-6)
pH - acidity of mineral water
HCO3 – concentration of HCO3 in mineral water (ml/l ·10-6)
Hydrogeochemical READINESS network:
Conduct. – electrical conductivity of mineral water (msim/cm)
Flow – flow of mineral water (l/sec)
Water temp. – temperature of mineral water (°C)
Radon in water – Radon concentration in mineral water (imp/min)
pHr – acidity of mineral water
Radon s/s – subsoil Radon concentration (imp/min)
INP – electromagnetic monitoring by “Irreversibility of non-stationary processes” techniques (W/m3)
Geomagn.- SD – geomagnetic observations – sinchronical difference (nTl)
Geomagn.- geomagnetic observations (nTl)
ULRW – ultra long radio waves (ionospheric observations) (mcs)
Fishes – activity of aquarium fishes (imp/hour)
Atm. press. – atmospheric pressure (mm)
* - non-amplitude anomaly
Type of precursor:
m/t – middle-term
sh/t – short-term
oper- operative
reaction – co-seismic and post-seismic effects
- 5 -
The earthquake precursors are characterized by the following features:
1. The preparation of earthquakes with various magnitudes (M³3.7) was accompanied
by the appearance of precursors in various physical fields (fig. 2 –5).
H, m
Fig. 2. The example of variations of underground waters level at the Dzorakhbyur station.
Spitak earthquake (Armenia, 07.12.1988, M=7.0).
2. The precursors and co-seismic effects were apparent at distances up to 1300 kms
from the epicenter of the seismic event (M=7.4) (fig. 6, 7).
H, m
Fig. 6. Variations of underground waters level at the Ashotsk station. Izmit earthquake
(Turkey, 17.08.1999, M=7.4).
M, F,
Fig. 7. Short-term anomalies and co-seismic effects of conductivity (a) and flow (b) of
mineralized waters in the Kajaran borehole. Izmit earthquake (Turkey, 17.08.1999, M=7.4).
M=7.0
D=87
P=1/1
short-term
short-term
co-seismic
effects
M=7.4
D=1.186
P=1/2
- 6 -
Active faults:
------------____well manifested
/////////////// buried
- - - - - - probable
º º º º º extension zones
epicentre of earthquake
• cities: Yerevan, Stepanakert
Methods and parameters of
monitoring:
HGD
subsoil Radon
geomagnetic
Geochemia of mineralized water:
pH
HCO3
Cl
He
Type of test-anomaly:
middle-term
short-term
reaction
Fig.3. An example of strong regional earthquakes (11 maps overall): sensitivity
of the observation stations in the territory of Armenia during the preparation and
realization of Spitak (Armenia, 07.12.1988, M=7.0) earthquake.
43.50 44.00 44.50 45.00 45.50 46.00 46.50 47.00 47.50
38.50
39.00
39.50
40.00
40.50
41.00
41.50
Spitak,
1988, M=7.0
Yerevan
Stepanakert
Goris
Dzorakhbyur
Noemberyan
Stepanavan
Metsamor
Sevan
Surenavan
Ararat
Kajaran
Gyumri
- 7 -
43.50 44.00 44.50 45.00 45.50 46.00 46.50 47.00 47.50
38.50
39.00
39.50
40.00
40.50
41.00
41.50
Active faults:
------------____well manifested
/////////////// buried
- - - - - - probable
º º º º º extension zones
epicentre of earthquake
• cities: Yerevan, Stepanakert
Methods and parameters of
monitoring:
HGD
INP
ULRW (ionospheric)
Type of test-anomaly:
short-term
operative
Fig. 4. An example of perceptible local earthquakes (7 maps overall): sensitivity
of the observation stations in the territory of Armenia during the preparation
and realization of Martuni (Armenia,10.12.1992, M=5.0) earthquake.
Stepanakert
Yerevan
Dzorakhbyur
Davidashen
Garni
Goris
Jermuk
Karchakhbyur
Martuni
1992, M=5.0
8
43.00 44.00 45.00 46.00 47.00
39.00
40.00
41.00
Noem
Ijev
Sevn
Kamo
Karc
Ekhg
Gors
Arta
Dzor
Mets
Azat Kuch
Shir
Amas
Asho
Active faults:
_______ well manifested
/ / / / / / / / buried
º º º º º extension zones
- - - - - probable
Type of test-anomaly:
middle-term
short-term
operative
reaction
Fig. 5. An example of sensitivity maps (15 maps overall): sensitivity of HGD-monitoring
stations in the territory of Armenia during the preparation and realization of earthquakes.
- 9 -
3. The quantity of precursors and the forms of their display are increased at the
last stage of the seismic event preparation (fig. 8, tables 2, 3).
The late short-term precursors (the time of display of precursor DT, that is a time
interval between the occurrence of a precursor and the seismic event, - more than 3 days
and less than 3 months) in the variations of 14 investigated parameters from 15 (except for
the flow) are showed. They reflected the preparation of 17 tested earthquakes out of 18.
More than half of 656 precursority anomalies and 167 test-anomalies make the late shortterm
precursors, and the forms of their display are diversified and cover almost all forms of
the anomal variations. It is obvious, that at the last stage of the preparation of seismic events
the quantity of precursors, anomal parameters and forms of display of precursority anomalies
are increased.
Fig. 8. The basic forms of display of the variations of underground waters level at the
NSSP stations: 1 and 2 - single local anomalies; 3 and 4 - infringements of the seasonal
course; 5 and 6 - double local anomalies; 7 and 8 - single positive anomalies at the
background of underground waters level falling; 9 and 10 - series of negative and positive
anomalies; 11 and 12 - sign-variable anomalies; 13 – step-form anomaly; non-amplitude
anomalies as an increase (14) and a full loss (15) of background dispersion. The vertical
arrows mark the moments of earthquakes.
1
3
2
4
5 6 7 8
9
10
11 12
13 14 15
10
- 10 -
Table 2
N Type of precursors
General
quantity of
anomalies
Quantity
of testanomalies
Pmin
Pmax
Pav
1 Long-term 0 0 --- --- ---
2 Middle-term 7 5 1/2 1/1 0.80
3 Short-term, including: 562 150 1/22 1/1 0.38
- HGD 142 43 1/15 1/1 0.51
- Geochemical (regional network): 140 28 1/22 1/1 0.29
- Geochemical (READINESS netw.): 60 11 1/12 1/1 0.30
- Subsoil Radon 128 31 1/13 1/1 0.40
- INP method 46 20 1/8 1/1 0.40
- Geomagnetic 15 5 1/5 1/2 0.37
- ULRW (ionospheric) 19 10 1/5 3/5 0.44
- Fishes activity 5 1 2/5 2/5 0.40
- Atmospheric pressure 7 1 2/7 2/7 0.29
4 Operative 87 12 1/15 1/1 0.25
TOTAL 656 167 1/22 1/1 0.48
Table 3
Regional
earthquakes
M 6.0
Local
earthquakes
M 3.7
N
Method (parameter)
Quantity
of tests Pav Quantity
of tests Pav
Pav
for all
earthquakes
1 HGD 33 0.52 15 0.47 0.50
Geochemical (READINESS network): 10 0.27 2 0.29 0.27
- flow of water 3 0.17 --- --- 0.17
2
- temperature 2 0.67 --- --- 0.67
- conductivity 2 0.33 2 0.29 0.31
- Radon 2 0.10 --- --- 0.10
- pH 1 0.08 --- --- 0.08
3 Geochemical (regional network): 20 0.31 8 0.27 0.30
- pH 3 0.26 1 0.33 0.28
- HCO3 2 0.27 2 0.35 0.31
- Cl 2 0.23 1 0.25 0.24
- He 13 0.48 4 0.15 0.39
4 Subsoil Radon 21 0.40 12 0.35 0.38
5 INP method 10 0.30 11 0.47 0.39
6 Geomagnetic 7 0.45 2 0.50 0.46
7 ULRW (ionospheric) 10 0.43 2 0.32 0.41
8 Fishes activity --- --- 2 0.40 0.40
9 Atmospheric pressure 1 0.29 --- --- 0.29
- 11 -
4. The sequence of earthquake preparation stages and the development of rock
destruction process are reflected in the anomalies of different orders available,
which are superimposed on each other with the approach of the moment of the
earthquake’s main fault formation (fig. 9 -11).
DT, nTl
Fig. 9.Variations of a synchronous difference of intensity DT of geomagnetic field at the
Stepanavan station (in relation to the Jermuk station) during the preparation of Narman
(Turkey, 30.10.1983, M=6.8) and Spitak (Armenia, 07.12.1988, M=7.0) earthquakes.
C, ml/l
Fig.10 Variations of helium concentration at the Karchakhbyur station. Baku earthquake
(Azerbaijan, 25.11.2000, M=6.5).
Dj, mcs
Fig. 11. The diagram of the difference of phases Dj for the indignant ionosphere on line
Reunyon-Yerevan. Ardebil earthquake (Iran, 28.02.97, M=6.7).
M=6.5
D=368
M=7.0
D=19
P=1/1
P=2/6
M=6.7
D=400
- 12 -
5. The forms of precursory occurrences are rather diverse (fig. 12 - 15):
5a. The middle-term precursors (DT more than 1 year) have the form of multimonthly
trend or multi-monthly sharp infringement of the seasonal trend
(fig. 12);
H, ò
Fig. 12. Variations of underground waters level at the Goris station. Spitak earthquake.
5b. The short-term precursors (DT - more than 3 days and less than 1 year) have
more often the form of short (from several weeks up to several months)
infringements of the seasonal trend, positive or negative smooth anomalies or
series of anomalies, steps, decrease or increase of a background dispersion
(fig. 13);
C, imp/min
Fig. 13. Variations of subsoil Radon contents at the Ekhegnadzor station. Bavra earthquake
(Armenia, 18.05.1994, M=4.2).
M=7.0
D=235
P=1/2
M=4.2
D=157
P=1/2
- 13 -
5c. The operative precursors (DT - 3 days and less) have usually the form of
spasmodic single or sign-variable pulses from several hours up to 1-2 days (fig.
14, 15).
Dj, mcs
Fig. 14. Variations of the difference of phases on the Reunyon-Yerevan line (at the
Davidashen station). Martuni earthquake (Armenia, 10.12.1992, M=5.0).
F, l/sec
Fig. 15. Variations of mineral waters flow in the Akhurik borehole. Bolvadin earthquake
(Turkey, 03.02.2002, M=6.2).
M=5.0
D=49
P=1/7
M=6.2
D=1.078
P=1/10
- 14 -
6. The dependence of the precursors distribution areas on magnitude of
earthquake (fig. 16, 17).
Fig. 16. The dependence of the precursors distribution areas on M:
a - early short-term, b - late short-term, c - operative precursors.
Fig. 17. The dependence of the precursors distribution areas on M: a - underground waters
level, b - helium, c - subsoil radon, d - irreversibility of non-stationary processes.
The areas of distribution of some precursority types are shown in fig.16, and in fig. 17
are the areas of distribution of some precursority parameters (methods), observed on the
territory of Armenia. The areas of different types and parameters (methods) are similar: they
are the zones limited to straight lines like y=kx+b. The distance of precursors display is
directly proportional with the magnitude of the earthquake.
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
à)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
á)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
â)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
ã)
a) b)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
â)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
á)
3 4 5 6 7 8
0
1
2
3
4
lg
M
,êì
1000
100
10
àa)) b) c)
c) d)
- 15 -
7. The dependence of the precursor amplitude on magnitude of earthquake and on
epicentral distance (fig. 18).
The dependence of A on Ì/lg with the use of all HGD precursors is shown as an
example on fig. 18a. Similar figures reflecting the absence of the dependence of A on M, on
and on Ì/lg , are characteristic for the other four parameters, too. For the exception of the
influence of different geological conditions of the positions of observation sites (the
neutralization of site- effects) on the investigated regularity, we have tried to receive the
amplitude dependences for each station separately in five investigated parameters. Certain
dependences of A on Ì/lgD have been received only for eight observation sites and for three
parameters (fig. 18 b, c).
Fig.18. The dependence of the HGD precursor's amplitude A on M (adjusted for epicentral
distance ) for all HGD precursors (a) and for separate HGD stations (b, c):
1 - Karchakhbyur, 2 - Ijevan, 3 - Noyemberyan, 4 - Ashotsk, 5 - Dzorakhbyur, 6 - Metsamor.
The amplitude of anomalies elaborated for these stations, grows with the increase of
Ì/lg . The successful results gained practically only on variations of the underground waters
level are probably explained by the fact, that the changes of the HGD-field is a direct
deformometer, reflecting the changes of intense-deformed conditions of the earth's crust. The
direct dependence of an amplitude of precursority anomalies on earthquake's magnitude at
certain stations may be explained by high sensitivity of these observation sites towards the
processes of the preparation of seismic events, which reflects the high activity of those parts
of active faults, within the range of which these monitoring sites are located.
1 2 3 4
0
20
40
60
80
M
â)
A, cì
lg
6
5
4
1 2 3 4
0
50
100
150
200
250
M
à)
A, cì
lg 1 2 3 4
0
50
100
150
200
250
M
á)
A, cì
lg
1
3
2
a) b) c)
- 16 -
Conclusion
It is desirable further:
1. to continue the testing of the strong regional and perceptible local seismic events
and, after the careful check and the analysis of data, to add the Catalogue of earthquakes
precursors by the tests of the occurrence of future seismic events;
2. to assess the reliability of precursors not only on monoparameter observations but
in their complexity as well, using the formulae of full probability and the Bayes’ formula, or the
theory of recognition of images.
References
1. Petrosyan H.M. et al. (1999) Testing of strong and perceptible earthquakes for current
seismic hazard assessment in the territory of Armenia. Gyumri, pp. 238-244 (in Russian).
2. Petrosyan H.M. (2004) Testing of seismic events in the territory of Armenia and adjacent
countries for current seismic hazard assessment. The thesis of Doctor of Geology. Yerevan,
181 p. (in Russian).
3. Petrosyan H. (2004) Testing of Izmit large earthquake in the territory of Armenia.
Proceedings of 5-th International Symposium on Eastern Mediterranean Geology.
Thessaloniki, Greece, 14 to 20 April 2004, pp. 635-638.
4. Sidorin A.Ya. The dependence of the times of occurrence of earthquakes’ precursors on
the epicentral distance. DAN SSSR, 1979, vol.245, N4, pp.825-828 (in Russian).
5. Sobolev G.A. The principles of the prediction of earthquakes. M., Nauka, 1993, 313 p. (in
Russian).
6. Zubkov S.I. The times of occurrence of earthquakes precursors. Izv. AN SSSR. The
Physics of the Earth. 1987, vol.5, pp. 87-91 (in Russian).


 

ELECTROMAGNETIC EARTHQUAKE PRECURSORS
MONITORING METHODOLOGY

R. Berkman, K. Hattori, V. Korepanow, A. Svenson
1. INTRODUCTION
Electromagnetic precursors of earthquakes seem to be the most often occurred and
the best studied phenomena. There are many practical evidences of the electric and
magnetic specific fluctuations occuring before the earthquakes. The most promising
frequency band of their occurrence, according to the information published, is ULF-ELF
range for short-term precursors and periods from hours to months for medium-term
precursors. The common feature of all these fluctuations is their small amplitude what
makes it difficult to extract them from the background noise – fluctuations, caused by the
sources of other nature. Also this imposes special requirements on the sensitivity
threshold and stability of the measuring devices - electrometers and magnetometers. By
this, taking into account that the corresponding measurements are executed mostly in
field conditions, both thermal and temporal stabilities of the measuring devices have to be
as high as possible.


The aim of this work is to review the existing methods of the measurement of
electric and magnetic fields oriented to the solution of the problem of short-term
precursors extraction.
2. MAGNETIC MEASUREMENTS
The analysis of the existing data gives approximative values of ELF-ULF magnetic
signals apparently connected with the seismic activity in the range from few nanotesla to
some tens of picotesla. It is evident that all wide amount of types of magnetometers is
restricted for this application to only four ones: flux-gate magnetometers (FGM), torsion
magnetometers (TM), search-coil magnetometers (SCM) and SQUID magnetometers.

The last magnetometer type, because of high price and practical impossibility to use such
magnetometer for constant monitoring, is not taken into account here. Further the noise
spectra of the three selected types of magnetometers – FGM, TM and SCM – will be
studied both theoretically and experimentally.

The most important feature of the magnetic measurements is the huge dynamic
band of the studied signal - more than 200 dB. Other peculiarity is the apparent nonsense
of the magnetic measurements, when the device with the reduced relative error about
0.1% is used for the measurements of the values about 0.0001% from its full scale. It
becomes clear that in such conditions the selection of parameters characterizing
measurement precision has to be specially made.

3. ELECTRIC MEASUREMENTS
The noise level of magnetometers, as it was stated above, is mostly determined by
the sensor parameters and construction. The electric measurements peculiarity is the
direct contact of the sensor – electrode – with the medium in which the measurements are
executed. It immediately causes the well known phenomenon – formation of double layer
along the electrode contact surface. That is why even ideally manufactured electrode will
have potential drop across this layer, highly dependent on many electrophysical
parameters. That is why not only electrodes manufacturing methodology but also the
methodology of their use during the experimental study of the electric fields in the
conductive media plays an important role for the decreasing of the electrometer noise.
The experimental study conducted at LCISR clearly stated that the lowest noise -
within the adopted frequency domain for the measurements in the ground - have copper
electrodes (Cu-Cu SO4), if the methodology of their application is properly followed.
This study showed that the potential difference of properly selected electrodes pair in
laboratory conditions can change within ±20÷30 microvolts within one mouth of
observations. Taking into account that the obtained in practice values of electric
precursory signals (from hundreds microvolts to some millivolts) are considerably
greater, such electrodes seem to solve the problem of their registration.

But unfortunately in field practice many other sources of noise and error appear.
Accepting that most of investigators make proper installation of the electrodes in the soil,
we shall analyse further only the methodological errors which give different ways of
electrodes mutual placing and of their matching with measuring amplifiers.

CONCLUSION
The presented study allows to recommend an optimal set of measuring
instrumentation and a corresponding methodology of its use for the monitoring of ELFULF
electric and magnetic fields, oriented at the selection of seismic activity
electromagnetic precursors.

For lower part of the frequency domain of interest FLUX GATE
MAGNETOMETER appear to be the best choice in order to get minimally possible noise
value.
But already starting from about 0.01 Hz and higher the SEARCH COIL
MAGNETOMETER overcome any other possible type of the magnetometers as to the
noise level. At higher frequency, say about 10 kHz, their noise can be as low as few
femtotesla (10-15T).

The Cu-CuCO4 electrodes are preferable for the measurements of electric fields in
the ground. The double 2-electrodes array is recommended for electrodes installation,
especially if long-term study of long-period electric field fluctuation is executed.

§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§

In this study we show that ELF-VLF emissions, occurred during the fracturing of rocks in laboratory experiments and during rock extraction in a cave, have been recorded also before earthquakes by 3 ELF-VLF stations installed in the central Apennines. The electric field emission is composed by two types of signals which are related to the evolution of micro-fractures within the rocks. During two years of monitoring, from August 2003 to September 2005, we recorded three signals, which all have frequency band and characteristics similar to the signal recorded during the fracturing of rocks in the laboratory and during the extraction of rocks. All these signals occurred a few days before the occurrence of M ≥ 4.5 earthquakes in the surrounding region (100-300 km). May be EM precursors of the earthquakes ?

Algae And Earthquake Precursors
 
Scienceblogs.com | May 27, 2007 9:44AM CDT
Concentrations of the natural pigment chlorophyll in coastal waters have been shown to rise prior to earthquakes. These chlorophyll increases are due to blooms of plankton, which use the pigment to convert solar energy to chemical energy via photosynthesis. This is based on an article by Singh et al. [doi] in Advances in Space Research, part of a special issue devoted to the use of satellite remote sensing for studying and predicting natural hazards such as earthquakes. The authors claim that you can detect a rise in sea surface temperature just before large coastal earthquakes. The blooms observed in this study are, they say, result from an increased flow of heat energy from the ocean to the atmosphere, enhancing the upwelling of cold, nutrient rich water and fuelling a boom in the growth of photosynthetic algae...
 

The Geomagnetic Quakes as Reliable Earthquake
Precursor- Beijing, Lanzhou regions, 2003
Mavrodiev S. Cht.
Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences
Tzarigradsko shose 72, 1784 Sofia, Bulgaria
mavrodi@inrne.bas.bg
China geophysicists, seismic earthquake data
China geophysicists, vector geomagnetic data, Beijing and Lanzhou geomagnetic observatories
The impressive development of the Earth sciences on the basis of new precise Crust condition
parameters measurements permits to estimates the probabilities for earthquakes risk. But the
prediction the time, epicentre and Magnitude of incoming earthquake is not a solved problem.
Many scientists are state that this is not solvable. Such pessimism is right because of very
scare time and space set of Crust movement parameters monitoring and the uncertainties of
our today knowledge for the Earth and its magnetic field. The local "when" Earthquake
prediction is based on the connection between geomagnetic "quakes" and the next incoming
minimum or maximum of tidal gravitational potential. The probability time window for the
predicted earthquake is +/-1 day for the minimum and +/-2 days for the maximum. The
preliminary statistic estimation on the basis of distribution of the time difference between
predicted and occurred earthquakes for the period 2002- 2003 for Sofia region and for Beijing
and Lanzhou regions are given. The solving of earthquake's prediction problem and creating
its theory need the efforts of wide interdisciplinary science group from physicists,
geophysicists, seismologists, Earth geomagnetism theory, Atmosphere and near space
physics, biologists, the application of temporary almost real time GIS for data acquisition,
visualization, archiving and analysis, the new possibilities for solving step by step the
nonlinear inverse problems for testing the adequateness of physical models and the reliability
of predictions. The monitoring should include standard geodetic data, seismic hazard map
developments, electromagnetic field monitoring under (electrical signals in VAN method and
its Thanassoulas's variant), on (electropotential distribution, geomagnetic variations) and over
(VLF and ULF, vertical electropotential distribution) Earth surface, atmosphere effects
(earthquake's clouds, electrical charge distribution), the behaviour of Earth radiation belts,
biological precursors. The statistical estimation for reliability of time, epicentre and
magnitude prediction is obligatory. The Beijing- Lanzhou regions are proposed as polygon for
testing the possibilities for creating research and short time earthquakes prediction
NETWORK. The important advantage of the proposal is that the geophysical seismic,
geomagnetic, atmosphere and near space monitoring exists and the research needs more
software than hardware for testing the approach and applying it in the practice. The main
trouble is not only the scepticism of the part of science community but also the governments
and municipalities authority's problems with practical using of the predictions for preventing
of the hazards consequences and the restoration of the environment. The creating of wide
international science group for every earthquake region can accelerate the solving of this
complex from science point of view and with high social meaning problem.
1. Introduction- short history of the developments of earthquake predictions results
1.1. Earthquake precursors- seismic, electromagnetic under on and over Earth surface,
atmospheric, Earth radiation belts anomalies, biological
2. Geomagnetic quake as reliable earthquake precursor
etc... etc...

 

Earthquake Precursors Pattern

 

The pattern of earthquake precursors could be described as “pre-earthquake signature” since every reoccurrence maintains common characteristics.

Observing the recorded signal, one realizes that, throughout the spectrum 1-40Hz, there are lines almost parallel to each other at any given time. Any spectral line differs 1-3Hz from the previous one.

The duration of the “pre-earthquake signature” varies from 15 minutes up to 3 hours, while for the specific spectral impressions, the intensity of the magnetic field varies from 2 PT to ~20 PT.

It is most likely that the common characteristics of the signal will still appear in frequencies higher than 40 Hz.  Further observation is required in order to verify the following:

  1. What is the co-relation of the magnetic field intensity and the duration of the “pre-earthquake signature” with the data of the earthquake to follow?
  2. What is the highest frequency in which the “pre-earthquake signature” is observed, as well as the co-relation of the frequency with the earthquake magnitude or the distance from the epicenter to the reception device?

Jedi Simon Cymatics Research

          Jedi Meta Search - Brainstorming and Integrate Scientific Procedures
                                                             QUAKE Notes 

page 1  jedi simon brainstorming

page 2  jedi simon brainstorming

page 3  jedi simon brainstorming

    page 4  jedi simon brainstorming

page 5  jedi simon brainstorming

 page 6  jedi simon brainstorming

page 8  jedi simon brainstorming

page 9  jedi simon brainstorming

page 9  jedi simon brainstorming                                                                                    

                                                                                                         use your heart ...  serve life