Jump to content

VAN method: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
usage of only summarized criticism
balancing by including the corresponding replies to the criticism (1996)
Line 39: Line 39:
Historically, the usefulness of the VAN method for prediction of earthquakes had been a matter of debate. Both positive and negative criticism on the early conception of the VAN method is summarized in the 1996 book "A Critical Review of VAN", edited by Sir James Lighthill<ref>{{Harvnb|Lighthill|1996}}.</ref> and in the debate issue of Geophysical Research Letters in 1996<ref name=":2" />. This criticism predates the time series analysis methods introduced by the VAN group in 2001. The main points of the criticism were:
Historically, the usefulness of the VAN method for prediction of earthquakes had been a matter of debate. Both positive and negative criticism on the early conception of the VAN method is summarized in the 1996 book "A Critical Review of VAN", edited by Sir James Lighthill<ref>{{Harvnb|Lighthill|1996}}.</ref> and in the debate issue of Geophysical Research Letters in 1996<ref name=":2" />. This criticism predates the time series analysis methods introduced by the VAN group in 2001. The main points of the criticism were:


* Predictive success: The first article in 1981 by VAN reporting the observation of the correlation between SESs and earthquakes was later criticized by claiming that the list of earthquakes used for the correlation was inaccurate.<ref>{{Cite journal|last=Wyss|first=M.|date=1996-05-27|title=Inaccuracies in seismicity and magnitude data used by Varotsos and Coworkers|url=http://onlinelibrary.wiley.com/doi/10.1029/95GL03545/abstract|journal=Geophysical Research Letters|language=en|volume=23|issue=11|pages=1299–1302|doi=10.1029/95GL03545|issn=1944-8007}}</ref> VAN refuted this criticism by explaining that the data came from the preliminary bulletin of the [[National Observatory of Athens]] (NOA), which is expected to differ somewhat from the final NOA bulletin and does not affect the observed correlation between SES and earthquakes.<ref>{{Cite journal|last=Varotsos|first=P.|last2=Eftaxias|first2=K.|last3=Lazaridou|first3=M.|last4=Dologlou|first4=E.|last5=Hadjicontis|first5=V.|date=1996-05-27|title=Reply to “Inaccuracies in seismicity and magnitude data used by Varotsos and Co-workers,” by M. Wyss|url=http://onlinelibrary.wiley.com/doi/10.1029/96GL01415/abstract|journal=Geophysical Research Letters|language=en|volume=23|issue=11|pages=1303–1306|doi=10.1029/96GL01415|issn=1944-8007}}</ref> An additional criticism was raised when examining the probability of chance correlation of 22 claims of successful predictions by VAN from January 1, 1987 through November 30, 1989.<ref>{{Cite journal|last=Wyss|first=M.|last2=Allmann|first2=A.|date=1996-05-27|title=Probability of chance correlations of earthquakes with predictions in areas of heterogeneous seismicity rate: The VAN Case|url=http://onlinelibrary.wiley.com/doi/10.1029/95GL03548/abstract|journal=Geophysical Research Letters|language=en|volume=23|issue=11|pages=1307–1310|doi=10.1029/95GL03548|issn=1944-8007}}</ref> This criticism was also refuted by VAN showing that the methodology on which the criticism was based cannot be correct since it rejects even an ideal prediction method (i.e., all earthquakes above a threshold are successfully predicted, no false alarm).<ref>{{Cite journal|last=Varotsos|first=P.|last2=Eftaxias|first2=K.|last3=Lazaridou|first3=M.|last4=Dologlou|first4=E.|last5=Hadjicontis|first5=V.|date=1996-05-27|title=Reply to “Probability of chance correlations of earthquakes with predictions in areas of heterogeneous seismicity rate: The VAN Case,” by M. Wyss and A. Allmann|url=http://onlinelibrary.wiley.com/doi/10.1029/96GL00908/abstract|journal=Geophysical Research Letters|language=en|volume=23|issue=11|pages=1311–1314|doi=10.1029/96GL00908|issn=1944-8007}}</ref>
* Predictive success: VAN has claimed to have observed at a recording station in Athens a perfect record of a one-to-one correlation between SESs and earthquake of magnitude ≥ 2.9 which occurred 7 hours later in all of Greece.<ref>{{Harvnb|Varotsos|Alexopoulos|Nomicos|1981b}}.</ref> However, it was later shown that the list of earthquake used for the correlation was false. Although VAN stated in their article that the list of earthquakes was that of the Bulletin of the [[National Observatory of Athens]] (NOA), it was found that 37% of the earthquakes actually listed in the bulletin, including the largest one, were not in the list used by VAN for issuing their claim. In addition, 40% of the earthquake which VAN claimed had occurred were not in the NOA bulletin.<ref>{{Harvnb|Wyss|1996a}}.</ref> Examining the probability of chance correlation of 22 claims of successful predictions by VAN of M > 4.0 from January 1, 1987 through November 30, 1989 <ref>{{Harvnb|Varotsos|Lazaridou|1991}}.</ref> it was found that 74% were false, 9% correlated by chance, and for 14% the correlation was uncertain.<ref>{{Harvnb|Wyss|Allmann|1996}}.</ref> No single event correlated at a probability greater than 85%, whereas the level required in statistics for accepting a hypothesis test as positive would more commonly be 95%.
* Proposed SES propagation mechanism: An analysis of the propagation properties of SES in the Earth’s crust showed that it is impossible that signals with the amplitude reported by VAN could have been generated by small earthquakes and transmitted over the several hundred kilometers between the epicenter and the receiving station.<ref>{{Harvnb|Bernard|1992}}; {{Harvnb|Bernard|LeMouel|1996}}.</ref> In effect, if the mechanism is based on piezoelectricity or electrical charging of crystal deformations with the signal traveling along faults, then none of the earthquakes which VAN claimed were preceded by SES generated an SES themselves. There is also some doubt that the phenomena of rock physics seen in certain laboratory conditions can be assumed to take place in the Earth’s seismogenic crust.
* Proposed SES propagation mechanism: An analysis of the propagation properties of SES in the Earth’s crust showed that it is impossible that signals with the amplitude reported by VAN could have been generated by small earthquakes and transmitted over the several hundred kilometers between the epicenter and the receiving station.<ref>{{Harvnb|Bernard|1992}}; {{Harvnb|Bernard|LeMouel|1996}}.</ref> In effect, if the mechanism is based on piezoelectricity or electrical charging of crystal deformations with the signal traveling along faults, then none of the earthquakes which VAN claimed were preceded by SES generated an SES themselves. There is also some doubt that the phenomena of rock physics seen in certain laboratory conditions can be assumed to take place in the Earth’s seismogenic crust.
* [[Electromagnetic compatibility]] issues: VAN’s publications are further weakened by failure to address the problem of eliminating the many and strong sources of change in the magneto-electric field measured by them, such as [[telluric current]]s from weather, and [[electromagnetic interference]] (EMI) from man-made signals.{{citation needed|date=July 2012}} One critical paper clearly correlates an SES used by the VAN group with digital radio transmissions made from a military base.<ref>{{Harvnb|Pham|Boyer|Chouliaras|Le Mouël|1998}}.</ref>
* [[Electromagnetic compatibility]] issues: VAN’s publications are further weakened by failure to address the problem of eliminating the many and strong sources of change in the magneto-electric field measured by them, such as [[telluric current]]s from weather, and [[electromagnetic interference]] (EMI) from man-made signals.{{citation needed|date=July 2012}} One critical paper clearly correlates an SES used by the VAN group with digital radio transmissions made from a military base.<ref>{{Harvnb|Pham|Boyer|Chouliaras|Le Mouël|1998}}.</ref>

Revision as of 11:48, 6 May 2016

The VAN method – named after P. Varotsos, K. Alexopoulos and K. Nomicos, authors of the 1981 papers describing it[1][2] – measures low frequency electric signals, termed "seismic electric signals" (SES), by which Varotsos and several colleagues claimed to have successfully predicted earthquakes in Greece.[3][4] Both the method itself and the manner by which successful predictions were claimed have been severely criticized[5][6][7] and defended[8].[9] The VAN method was advanced in 2001 by improving the accuracy of the estimation of the time window of the forthcoming earthquake through a new analysis, termed "natural time", and this new method has been independently reproduced.[10] Research related to the VAN method is currently carried out at the Solid Earth Physics Institute, University of Athens, Greece.

Description of the VAN method

Prediction of earthquakes with this method is based on the detection, recording and evaluation of seismic electric signals or SES. These electrical signals have a fundamental frequency component of 1 Hz or less and an amplitude the logarithm of which scales with the magnitude of the earthquake.[11] According to VAN proponents, SES are emitted by rocks under stresses caused by plate-tectonic forces. There are three types of reported electric signal:[12]

  • Electric signals that occur shortly before as a major earthquake. Signals of this type were recorded 6.5 hours before the 1995 Kobe earthquake in Japan, for example.[13]
  • Electric signals that occur some time before a major earthquake.
  • A gradual variation in the Earth's electric field some time before an earthquake.

Several hypotheses have been proposed to explain SES:

  • Stress-related phenomena: Seismic electric signals are perhaps attributed to the piezoelectric behaviour of some minerals, especially quartz, or to effects related to the behavior of crystallographic defects under stress or strain. Series of SES, termed SES activities (which are recorded before major earthquakes), may appear a few weeks to a few months before an earthquake, when the mechanical stress reaches a critical value.[14][2] The generation of electric signals by minerals under high stress leading to fracture has been confirmed with laboratory experiments.[15]
  • Thermoelectric phenomena: Alternately, Chinese researchers proposed a mechanism which relies on the thermoelectric effect in magnetite.[16]
  • Groundwater phenomena: Three mechanisms have been proposed relying on the presence of groundwater in generating SES. The electrokinetic effect is associated with the motion of groundwater during a change in pore pressure.[17] The seismic dynamo effect is associated with the motion of ions in groundwater relative to the geomagnetic field as a seismic wave creates displacement. Circular polarization would be characteristic of the seismic dynamo effect, and this has been observed both for artificial and natural seismic events.[18] A radon ionization effect, caused by radon release and then subsequent ionization of material in groundwater, may also be active. The main isotope of radon is radioactive with a half-life of 3.9 days, and the nuclear decay of radon is known to have an ionizing effect on air. Many publications have reported increased radon concentration in the vicinity of some active tectonic faults a few weeks prior to strong seismic events.[19] However, a strong correlation between radon anomalies and seismic events has not been demonstrated.[20]

While the electrokinetic effect may be consistent with signal detection tens or hundreds of kilometers away, the other mechanisms require a second mechanism to account for propagation:

  • Signal transmission along faults: In one model, seismic electric signals propagate with relatively low attenuation along tectonic faults, due to the increased electrical conductivity caused either by the intrusion of ground water into the fault zone(s) or by the ionic characteristics of the minerals.[21]
  • Rock circuit: In the defect model, the presence of charge carriers and holes can be modeled as making an extensive circuit.[22]

Seismic electric signals are detected at stations which consist of pairs of electrodes (oriented NS and EW) inserted into the ground, with amplifiers and filters. The signals are then transmitted to the VAN scientists in Athens where they are recorded and evaluated. Currently the VAN team operates 9 stations, while in the past (1989) they could afford up to 17.[23]

The VAN team claimed that they were able to predict earthquakes of magnitude larger than 5, with an uncertainty of 0.7 units of magnitude, within a radius of 100 km, and in time window ranging from several hours to a few weeks. Several papers confirmed this success rate, leading to statistically significant conclusion.[24] For example, there were eight M ≥ 5.5 earthquakes in Greece from January 1, 1984 through September 10, 1995, and the VAN network forecast six of these.[25]

The VAN method has also been used in Japan,[26] but in early attempts success comparable to that achieved in Greece was "difficult" to attain.[27] A preliminary investigation of seismic electric signals in France led to encouraging results.[28]

Earthquake prediction using "natural time" analysis

Since 2001 the VAN team has tried to improve the accuracy of the estimation of the time of the forthcoming earthquake. They introduced the concept of natural time, a time series analysis technique which puts weight on a process based on the ordering of events.[29] They define a variance term κ in the time series analysis. Their current method detects SES as valid when κ = 0.070. That is the first part of their protocol. Once the SES are deemed valid, a second analysis is started. All the subsequent seismic (not electric) events are now noted, and the region is divided up as a Venn diagram with at least two seismic events per overlapping rectangle. When the distribution of κ for the rectangular regions has its maximum at κ = 0.070, the critical seismic event is imminent, i.e. it will occur in a few days to one week or so, and a prediction is issued.[30]

Here is a brief introduction to the analysis. Two terms, χ and Q, characterize each event. χ is defined as k/N, where k is an integer (the k-th event) and N is the total number of events in the time sequence of data. Q is the energy released for each event. A related term, p, is the ratio Q/Q(total), and describes the fractional energy released. They introduce a critical term κ, defined as κ = [ Σ p(k) (k/N)^2 ] − [ Σ (k/N) p(k) ]^2, where κ is the variance in natural time, k is the event index, N is the number of events, Σ signifies summation, for k = 1 to N, and p(k) is the fractional quantity of energy for the k-th event. This variance puts extra weight on the energy term, p(k).

The VAN team published to have successfully predicted twenty five of the 28 major earthquakes from 2001 through 2010 in the region of latitude N 36° to N 41° and longitude E 19° to E 27° with this new analysis.[31] Predictions are part of papers housed in arXiv, and new predictions are currently being made and uploaded there.[32] A description of the updated VAN method was collected in a book published by Springer in 2011, titled "Natural Time Analysis: The New View of Time."[33]

Results

Natural time analysis revealed an additional fact showing the physical connection of SES activities with earthquakes as follows: Taking the view that the earthquake occurrence is a phase-change (critical phenomenon) where the new phase is the mainshock occurrence, the above mentioned variance term κ is the corresponding order parameter.[33] The κ value calculated for a window comprising a number of seismic events comparable to the average number of earthquakes occurring within a few months, fluctuates when the window is sliding through a seismic catalogue. It was found that these κ fluctuations exhibit a minimum a few months before a mainshock occurrence and in addition this minimum occurs simultaneously with the initiation of the corresponding SES activity. Such a simultaneous appearance of two precursory phenomena in independent datasets of different geophysical observables (electrical measurements, seismicity) has been observed for the first time in the literature.[34] Furthermore the natural time analysis of the seismic catalogue of Japan during the period from January 1, 1984 until the occurrence of the magnitude 9.0 Tohoku earthquake on March 11, 2011, revealed that such clear minima of the κ fluctuations appeared before all major earthquakes with magnitude 7.6 or larger. The deepest of these minima was found to occur on January 5, 2011, i.e., almost two months before the Tohoku earthquake occurrence.[35] Finally, by dividing the Japanese region into small areas, it has been found that some small areas show minimum of the κ fluctuations almost simultaneously with the large area covering the whole Japan and such small areas clustered within a few hundred kilometers from the actual epicenter of the impending major earthquake.[36][37]

Criticisms of VAN

Historically, the usefulness of the VAN method for prediction of earthquakes had been a matter of debate. Both positive and negative criticism on the early conception of the VAN method is summarized in the 1996 book "A Critical Review of VAN", edited by Sir James Lighthill[38] and in the debate issue of Geophysical Research Letters in 1996[9]. This criticism predates the time series analysis methods introduced by the VAN group in 2001. The main points of the criticism were:

  • Predictive success: The first article in 1981 by VAN reporting the observation of the correlation between SESs and earthquakes was later criticized by claiming that the list of earthquakes used for the correlation was inaccurate.[39] VAN refuted this criticism by explaining that the data came from the preliminary bulletin of the National Observatory of Athens (NOA), which is expected to differ somewhat from the final NOA bulletin and does not affect the observed correlation between SES and earthquakes.[40] An additional criticism was raised when examining the probability of chance correlation of 22 claims of successful predictions by VAN from January 1, 1987 through November 30, 1989.[41] This criticism was also refuted by VAN showing that the methodology on which the criticism was based cannot be correct since it rejects even an ideal prediction method (i.e., all earthquakes above a threshold are successfully predicted, no false alarm).[42]
  • Proposed SES propagation mechanism: An analysis of the propagation properties of SES in the Earth’s crust showed that it is impossible that signals with the amplitude reported by VAN could have been generated by small earthquakes and transmitted over the several hundred kilometers between the epicenter and the receiving station.[43] In effect, if the mechanism is based on piezoelectricity or electrical charging of crystal deformations with the signal traveling along faults, then none of the earthquakes which VAN claimed were preceded by SES generated an SES themselves. There is also some doubt that the phenomena of rock physics seen in certain laboratory conditions can be assumed to take place in the Earth’s seismogenic crust.
  • Electromagnetic compatibility issues: VAN’s publications are further weakened by failure to address the problem of eliminating the many and strong sources of change in the magneto-electric field measured by them, such as telluric currents from weather, and electromagnetic interference (EMI) from man-made signals.[citation needed] One critical paper clearly correlates an SES used by the VAN group with digital radio transmissions made from a military base.[44]
  • Scientific reporting: Most importantly, the method is hindered by a lack of statistical testing of the validity of the VAN hypothesis because the researchers keep changing the parameters of the hypothesis (the moving goalpost technique).[45]
  • Public policy: Finally, one requirement for any earthquake prediction method is that, in order for any prediction to be useful, it must predict a forthcoming earthquake within a reasonable time-frame, epicenter and magnitude. If the prediction is too vague, no feasible decision (such as to evacuate the population of a certain area for a given period of time) can be made. In practice, the VAN group issued a series of telegrams in the 1980s, warning of impending earthquakes that did not occur, or did not occur within the parameters listed in the telegrams. During the same time frame, the technique also missed major earthquakes.[46] These inaccurate predictions from the early VAN method led to public criticism and the cost associated with false alarms generated ill will.[47] Major opponents of VAN were the Greek seismologists Vassilis Papazachos and G. Stavrakakis. The debate between Papazachos and the VAN team has repeatedly caused public attention in their home country Greece and has been extensively discussed in the Greek media.[citation needed]

See also

Notes

  1. ^ Varotsos, Alexopoulos & Nomicos 1981a, 1981b
  2. ^ a b Varotsos & Alexopoulos 1984
  3. ^ Varotsos & Kuhlanek 1993 (preface to a special edition about VAN)
  4. ^ Varotsos, Alexopoulos & Lazaridou 1993
  5. ^ Mulargia & Gasperini 1992
  6. ^ Geller 1997, §4.5
  7. ^ ICEF 2011, p. 335
  8. ^ Lighthill 1996 (proceedings of a conference that reviewed VAN)
  9. ^ a b twenty articles in a special issue of Geophysical Research Letters (table of contents)
  10. ^ Uyeda, Kamogawa & Tanaka 2009
  11. ^ Varotsos, Alexopoulos & Nomicos 1981a; Varotsos et al. 1981; Varotsos, Alexopoulos & Nomicos 1982.
  12. ^ Varotsos, Alexopoulos & Lazaridou 1993.
  13. ^ Matsumoto, Ikeya & Yamanaka 1998.
  14. ^ Varotsos et al. 1986, p. 120.
  15. ^ Hadjicontis et al. 2007
  16. ^ Shen et al. 2011.
  17. ^ Gershenzon, Gokhberg & Yunga 1993.
  18. ^ Honkura et al. 2009.
  19. ^ Pulinets 2007.
  20. ^ ICEF 2011, p. 334.
  21. ^ Varotsos et al. 1998.
  22. ^ Freund 1998.
  23. ^ Varotsos & Lazaridou 1991
  24. ^ Hamada 1993
  25. ^ Uyeda 1996
  26. ^ Uyeda, Kamogawa & Tanaka 2009
  27. ^ Utada 1993, p. 153
  28. ^ Maron et al. 1993
  29. ^ Varotsos, Sarlis & Skordas 2002; Varotsos 2006.
  30. ^ Varotsos, Sarlis & Skordas 2011, Chapter 7.
  31. ^ Varotsos, Sarlis & Skordas 2011, p. 326
  32. ^ http://arxiv.org/, search for "Varotsos"
  33. ^ a b Varotsos, Sarlis & Skordas 2011.
  34. ^ Varotsos, P. A.; Sarlis, N. V.; Skordas, E. S.; Lazaridou, M. S. (2013-03-18). "Seismic Electric Signals: An additional fact showing their physical interconnection with seismicity". Tectonophysics. 589: 116–125. doi:10.1016/j.tecto.2012.12.020.
  35. ^ Sarlis, Nicholas V.; Skordas, Efthimios S.; Varotsos, Panayiotis A.; Nagao, Toshiyasu; Kamogawa, Masashi; Tanaka, Haruo; Uyeda, Seiya (2013-08-20). "Minimum of the order parameter fluctuations of seismicity before major earthquakes in Japan". Proceedings of the National Academy of Sciences. 110 (34): 13734–13738. doi:10.1073/pnas.1312740110. ISSN 0027-8424. PMC 3752201. PMID 23918353.
  36. ^ Sarlis, Nicholas V.; Skordas, Efthimios S.; Varotsos, Panayiotis A.; Nagao, Toshiyasu; Kamogawa, Masashi; Uyeda, Seiya (2015-01-27). "Spatiotemporal variations of seismicity before major earthquakes in the Japanese area and their relation with the epicentral locations". Proceedings of the National Academy of Sciences. 112 (4): 986–989. doi:10.1073/pnas.1422893112. ISSN 0027-8424. PMC 4313817. PMID 25548194.
  37. ^ Huang, Qinghua (2015-01-27). "Forecasting the epicenter of a future major earthquake". Proceedings of the National Academy of Sciences. 112 (4): 944–945. doi:10.1073/pnas.1423684112. ISSN 0027-8424. PMC 4313830. PMID 25583499.
  38. ^ Lighthill 1996.
  39. ^ Wyss, M. (1996-05-27). "Inaccuracies in seismicity and magnitude data used by Varotsos and Coworkers". Geophysical Research Letters. 23 (11): 1299–1302. doi:10.1029/95GL03545. ISSN 1944-8007.
  40. ^ Varotsos, P.; Eftaxias, K.; Lazaridou, M.; Dologlou, E.; Hadjicontis, V. (1996-05-27). "Reply to "Inaccuracies in seismicity and magnitude data used by Varotsos and Co-workers," by M. Wyss". Geophysical Research Letters. 23 (11): 1303–1306. doi:10.1029/96GL01415. ISSN 1944-8007.
  41. ^ Wyss, M.; Allmann, A. (1996-05-27). "Probability of chance correlations of earthquakes with predictions in areas of heterogeneous seismicity rate: The VAN Case". Geophysical Research Letters. 23 (11): 1307–1310. doi:10.1029/95GL03548. ISSN 1944-8007.
  42. ^ Varotsos, P.; Eftaxias, K.; Lazaridou, M.; Dologlou, E.; Hadjicontis, V. (1996-05-27). "Reply to "Probability of chance correlations of earthquakes with predictions in areas of heterogeneous seismicity rate: The VAN Case," by M. Wyss and A. Allmann". Geophysical Research Letters. 23 (11): 1311–1314. doi:10.1029/96GL00908. ISSN 1944-8007.
  43. ^ Bernard 1992; Bernard & LeMouel 1996.
  44. ^ Pham et al. 1998.
  45. ^ Mulargia & Gasperini 1992; Mulargia & Gasperini 1996; Wyss 1996b.
  46. ^ Hamada 1993.
  47. ^ Mulargia & Geller 2003, p. 318.

References