International Journal of Geosciences
Vol.09 No.12(2018), Article ID:89580,19 pages

Study of the Possibility of Predicting Earthquakes

Strachimir Cht Mavrodiev1, Lazo Pekevski2, Emil Botev3, Ali Pinar4, Giorgi Kikuashvili5, Alexader Vol6, Arie Gilat7

1INRNE-BAS, Sofia, Bulgaria

2Skopje Seismological Observatory, Skopje, Macedonia

3NIGGG BAS, Sofia, Bulgaria

4Department of Earthquake Engineering, Bogaziçi University, Istanbul, Turkey

5Tbilisi, Giorgia

6OSHADI Drug Administration Ltd., Rehovot, Israel

7Geological Survey of Israel, Jerusalem (ret.), Israel

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: October 30, 2018; Accepted: December 26, 2018; Published: December 29, 2018


It is already well known that the “when, where and how strong” earthquake prediction problem cannot be solved by only analyzing the database from former earthquakes. A possible solution to this problem is proposed herein based on the analysis of the physicochemical processes as participants in earthquake preparation and on the characteristic rate of reflection of these processes on the Earth’s surface. The proposed procedure includes monitoring of correlation of electromagnetic fields variations with tidal waves. This solution provides a way of selecting a complex of reliable earthquake precursors using the Inverse Problem Method for earthquakes which will occur in the region around the monitoring point (radial distance ≈ 700 km) in the next seven-day period [1] .


Earthquakes Prediction, Earth Degassing, Physical Chemistry, Explosions, Electromagnetic Field Monitoring

1. Introduction

In this article, tectonic earthquake and eruption processes are presented and discussed as a series of chemical explosions caused by physicochemical processes, which are partly reflected on Earth’s surface. Energy accumulation in the hypocenter includes the latent energy of chemically active (explosive) substances which induce exothermic syntheses (explosions) depending on local changes of Pressure-Temperature (PT). Also discussed are chemical explosions which are presented here as branched chain reactions that form concentrations of highly explosive clusters whose explosions generate earthquake hypocenters. Those intermediate hypocenters (foreshocks - major shock - aftershocks) ascend along faults and fractures from a depth of about 700 km to the Earth’s surface, becoming the complex of an earthquake.

Cold nuclear synthesis (fusion) and natural fission reactions are also presented as the major internal energy sources.

Whole assemblage of hypocenter preparation processes is accompanied by the generation of electromagnetic fields, which in contrast to other processes, are instantly reflected on the Earth’s surface.

Presented are principles and results of forecasting of regional imminent seismic activity based on the analysis of one minute of INTERMAGNET geomagnetic field data and NASA codes for Sun-Moon tides―Geomagnetic Quake (GQ). Examples of prediction of the period, magnitude, depth, and coordinates of the hypocenter of an impending earthquake are based on the Inverse Problem Method for the analysis of monitoring data variations of geoelectromagnetic fields. The necessary and sufficient conditions for the existence of a solvable inverse problem are formulated based on the Dubna method for discovering the hidden dependences. The accuracy of prediction depends on the values of depth, coordinates, time, magnitude of the impending earthquake, number of monitoring points, geology of the region, and on the ill-posed quality of the received overdetermined non-linear algebraic system. Monitoring of additional data provides a way of selecting a complex of reliable earthquake precursors.

2. The Reason for Tectonic Earthquakes Is a Series of Chemical Explosions

The classical Reid’s model [2] of “earthquake as a result of rock displacement under accumulated elastic stress” cannot explain the observed earthquake-cycle, foreshocks - main shock - aftershocks, of the strong earthquake, and from where comes its monstrous energy [3] [4] . There are clear demonstrations that the reasons for most of the “tectonic earthquakes” are underground natural explosions which cause the acceleration in the vertical direction exceeding the acceleration of gravity [5] and are shown in Table 1. Anomalous flow of the Earth’s core lower mantle helium enriched by its light isotope 3He, and hydrogen, which accompanies earthquakes and volcanic eruptions, provides evidence that volcanic eruption is a variety of “tectonic earthquake” wherein its hypocenter rises to the Earth’s surface. Both events are a release of high density energy from the deep-rooted mantle-plumes through physicochemical reactions.

“There is evidence, which should be regarded less as established fact than as working hypothesis, that in the neighborhood of the epicenter the vertical component of motion is larger relative to the horizontal components than elsewhere” [6] . The elastic energy of “semi-solid mantle and lithosphere” breaks cannot

Table 1. Database of ground motion components in the neighborhoods of the epicenters.

cause such motion and energy release in principle [7] . The energy release and the dominant vertical component of Earth’s surface motion [6] [8] [9] support the claim that explosion is a basic mechanism of tectonic earthquakes.

Flow of H and He from the core and the lower mantle where pressure-temperature (PT) exceeds 1,000,000 atm. and 5000 K is undeniably a major energy source [3] [4] . This energy is:

1) Quasi-constantly released during billions of years of the Earth’s existence and practically limitless;

2) Can be quickly concentrated and focused;

3) Is of very high density;

4) Offers very high velocities of energy release;

5) Has small losses during transportation over long distances.

Authors have proposed [3] [4] a conceptual system of hypotheses, which explains that during Earth’s accretion, primordial helium and hydrogen were trapped and stored in the planet’s interior as He- and H-interstitial solutions and compounds, stable only under ultrahigh PT-conditions, which were discovered in recent experiments. The endothermic reactions of their generation provided effective cooling of the planet and prevented its evaporation, where the end-products of those reactions were more compact than the initial gases. Since the stabilization of our planet, exothermic processes of H and He degassing became dominant, releasing the energy invested in their generation. The specific energy of the core-mantle H and He was calculated with 3He serving as a unique measuring transformer correlative to the internal heat flow. Multiplying its flow from the lower mantle by the highest coefficient of correlation results in 5.12 × 1020 J/year, an amount of energy five-fold greater than the entire energy loss involved in earthquake and volcanic activity [3] [4] .

In distinction to other main sources of the Earth’s internal energy (cold fusion and fission nuclear reactions, radioactive decomposition of U, Th and 40K, gravitational differentiation in the Earth’s liquid core and the energy of lunar tides) the chemical energy can be concentrated and focused in the mantle-plumes supplied hypocenters generating great earthquakes and volcanic eruptions [3] [4] .

In summary: H- and He-sublimation from the solid and convection in the liquid core with flux-melting the solid mantle, generates gas-liquid scavenging plumes. H- and He-release are accompanied by an intense release of their stored specific (latent) energy, and their ionization and incorporation into different chemical compounds was followed by decomposition due to local and gradual PT changes. Ongoing compressions-decompressions (foreshocks - major shock - aftershocks) within upward moving hypocenters (magma chambers), accompanied by additional releases of energy, cause release of elemental H, O, C, S, Cl, F etc., and induce synthesis of H2O, SO2, H2SO4, CO2, H2S, HCL, HF and other compounds in accordance with local changes of PT and energy, in accordance with Principle Le Chatelier-Braun (“The Equilibrium Law”: “whenever a system in equilibrium is disturbed the system will adjust itself in such a way that the effect of the change will be nullified”).

3. Physics and Chemistry of the Hypocenter Preparation and Earthquake Prediction

The earthquake preparation includes a number of processes having different rates. Only the electromagnetic waves approach Earth’s surface immediately. This difference in the rates allows us to use them as a reliable earthquake precursor. A possible solution of this problem is proposed herein based on the analysis of the physicochemical processes as participants in earthquake preparation and on the characteristic rate of reflection of these processes on the Earth’s surface. This solution provides a way of selecting a complex of reliable earthquake precursors using the Inverse Problem Method for earthquakes which will occur in the region around the monitoring point (radial distance ≈ 700 km) in the next seven day period.

N.N. Semenov declared already in 1956 in his Nobel Prize lecture titled “Some problems relating to chain reactions and to the theory of combustion” that the trains of chemical explosions are chemical branched chain-reactions [10] . This declaration is supported by a comparison of seismograms from earthquakes and nuclear explosions where the complexity of natural events (earthquakes) is higher than that of artificial events (explosions). Micro- or macro-foreshocks are forerunners of the major shock. Natural earthquakes are more complex than nuclear explosions and at teleseismic distances the difference between them is clear. This difference is observed very clearly in the relationship of solids to surface-wave amplitudes [11] . The nuclear weapons test is just a point explosion, whereas earthquake is the superposition of the totality of explosions which are distributed in space and time.

Prerequisites for the chemical explosions are the critical concentrations of reactants and their ratio which depends on PT conditions [10] [12] [13] . The critical concentration of the reactants and critical size of the explosive substances cluster is the first necessary condition of the local explosion [14] .

The possibility of explosion propagation (or detonation) to other clusters depends on the distance between clusters or on the cluster volume concentration. The critical or more than critical concentration of ready to detonate volume of explosive substances clusters is the second condition of earthquake. Too large a distance between clusters limits the propagation of detonation possibility due to the possibility of absorption of the local explosion’s energy by surrounding matter. This absorption causes local heating of matter and the formation of the chemically active substances [15] .

Relatively small concentrations of explosive clusters before an earthquake produce foreshocks, which is preparation for an earthquake’s major shock. “Combustion” of most of the clusters during the earthquake process decreases their concentration and generates aftershocks, (which take part in the rising hypocenter) and cause the relaxation of the surrounding matter.

Formation and accumulation of the explosive substances cluster, and preparation of the earthquake, is a totality of process. A hypocenter is an open thermodynamic system which uses all the possible degrees of freedom. This system is non-linear due to a principally different rate of separate processes including diffusion and filtration of molten matter through porous rock and cracks, heating and cooling, and stress and strain flow. An earthquake may be described as a bifurcation which returns part of the mantle-lithospere system to their main trajectory of development and corresponds to minimal internal free energy of the system and maximal rate of entropy production in the macrosystem. The following possible reactants participate in an earthquake explosion: hydrogen - oxygen; hydrogen - halogens; hydrogen - sulfur; alkanes (methane, etc.) - oxygen; alkanes - halogens; alkanes - nitrates, etc.

Explosive substances are produced and accumulated due to the energy which is released in the Earth’s core, mantle, and lithosphere, by five main sources:

1) Cold nuclear synthesis (fusion reactions), which are accompanied by a generation and release of energy, 3He, 4He, 3H and Earth neutrinos [16] - [23] .

2) The natural fission nuclear reactors with fast neutrons on the boundary of Earth’s solid/liquid core, and possibly, liquid core/mantle [24] [25] [26] . The capacity of those reactors depends on the Sun-Axions flow-intensity [27] .

3) Tidal waves cause dissipation of energy in the Earth’s core, mantle, and lithosphere; tidal waves are the main source of energy on the Jupiter and Saturn moons [28] .

4) Gravitational differentiation promotes the solid core formation and plume activity [29] .

5) Earth’s degassing of hydrogen and helium [3] [4] [29] [30] [31] , which generates anomalies of these gases in vicinities of active faults and forms a halo of hydrogen and helium surrounding our planet and comprising its exosphere.

Results of the analyses of the earthquake and volcanic eruption related gases clearly show that the mechanism of quake and eruption is a release of chemical energy. Physicochemical processes are the most effective way of transformation, transportation, and accumulation of energy into the hypocenter chamber.

Total heat losses of the Earth are 46 ± 3 TW [32] including about 20 TW of the released radiogenic heat. According to neutrino flow measurements radiogenic component in heat losses would be near 19.9 ± 9 TW [16] [26] . Heat from the core is about 8 TW ( [32] table 11). Heat flow from the convective mantle including gravitational differentiation and degassing is about 39 TW [3] [4] [32] .

The average annual energy of earthquakes is about 0.44 TW and energy of tidal wave dissipation is about 0.1 TW ( [32] table 13). Heat and radiation produced by the nuclear reactions cause dissociation and an increase of reactants temperature, which create thermal currents due to the Seebeck effect. Partial or complete melting of solid matter generates electrochemical processes on phase boundaries and corresponding galvanic currents. Local heating increases local pressure and accelerates the movement of semi-melt or melt matter through porous rock and cracks creating the streaming potential, sedimentation potential, and corresponding currents. Temperature gradient of the heat flow causes electric gradient or electrode potentials on phase boundaries, which create electrophoresis and electroosmosis. Hydrogen, oxygen, sulfur, halogens, carbon, nitrogen, aluminum, and alkaline metals are possible participants in energy transport due to their chemical activity. The presence of helium increases flexibility and rate of stream of solids and melts [30] [33] [34] . Moreover, hydrogen, lithium, and boron may take part in nuclear fusion reactions [16] - [21] [35] [36] . Movement of the matter creates triboelectricity [37] due to friction between boundaries of the mantle-fragments and of gas-liquid plume and mantle, and also promotes electrokinetic phenomena (electrophoresis and electroosmosis).

4. Sun-Moon Tides as Earthquake Triggers

Tidal waves cause dissipation of energy in the mantle and lithosphere [28] [38] [39] [40] and periodic stress-strain waves create a peristaltic effect and increase the rate of rising of plume matter. Tidal waves have a very high velocity (more than 460 m/s in the lithosphere and more than 240 m/s on the liquid core boundary in the equatorial zone). This velocity is higher than the critical rate of propagation of brittle cracks. A high rate of loading promotes brittle behavior of the semi-solid viscose mantle and generates cracks. Cracks and cavities are filled by melt, steam, gas, suspension, etc. Coexistence of the liquid and solid phase provides “adiabatic” heat transport with maximal efficacy:

K = n V λ C V 3 N A (1)

where: K is thermal conductivity; n V is number of particles per unit volume, which is close for liquid and solid phase near melting point; λ is mean concentration of free particles (electrons, for example); CV is molar heat capacity; NA is Avogadro’s number.

It is well known, that the specific heat capacity during phase transition is two orders of magnitude higher than for the mixed liquid and solid phase [41] [42] [43] [44] due to latent heat of the heat transition. It means that solid-liquid state is thermodynamically preferable for the mantle matter and for the earthquake-hypocenter heat-transport. Multiphase structure is preferable also for the matter-transport due to significantly higher flexibility and super-elasticity during phase transition [43] [44] .

Thermal-current loops at the phase boundaries can enhance the local magnetic field created by the matter flows and by physicochemical processes.

B = 1 n μ 0 4 π r n 2 d l n r ^ n (2)

where: B is magnetic field intensity (a vector); dln is differential element of the current in the direction of the corresponding current (vector); r is distance from the current to the point of magnetic field measurement; r ^ n is the unit vector from the current element to the point of magnetic field measurement.

Tidal waves are a natural trigger of the totality of physicochemical processes, which comprise earthquake-preparation [45] [46] [47] . Most of these processes are accompanied by electromagnetic phenomena. The rate of the magnetic field propagation is 300,000 km/sec. This means that the geomagnetic signal approaches the Earth’s surface without any delay. However, the time taken for the relaxation processes, for creating electrical currents, and for changing the local geomagnetic field is much longer than that of the magnetic field propagation. The rate of detonation at atmospheric pressure varies from 3 to 11 km/s (more than the velocity of sound) and the rate of the longitudinal and transverse waves in the solid mantle changes from 8 to 13.5 km/s for P-waves and from 4.5 to 7 km/s for S-waves. The time propagation of shock wave from hypocenter to epicenter is a function of depth and may vary between 1 to 90 seconds. The rate of all other processes may be much smaller. For example, the rate of plume matter movement, diffusion, or filtration through fractured or porous rock may be very low also.

So, processes of earthquake-hypocenter preparation comprise a multi-parametric non-linear system, which compensates differences in times of response or relaxation of different processes by bifurcation.

The monitoring system has to use parameters with a characteristic time of response equal or shorter than the duration of hypocenter matter relaxation. Moreover, the time it takes for the measurement of these parameters has to be shorter than the time of earthquake preparation. Time and rate of the processes involved are variable and may accelerate toward earthquake or bifurcation. It means that relatively short-time reliable predictions may be based only on monitoring the changes of the electromagnetic fields and viscous-elastic waves as a response to tides only.

For the longer time prognosis, other reliable precursors have to be included (see below).

5. Development of the “Geomagnetic Quake” Earthquake Prediction Method

・ The beginning of this investigation dates back to the period of the International Scientific Group work in the framework of the “Clean and Peaceful Black Sea” program (1986-1991)―Crimea 1927 earthquakes. In 1989 INRNE-BAS started an investigation of the physics of the sea flames, observed at the time of the Crimea 1927 earthquake with the help of the magnetometer, with a relative accuracy of less than or equal to 1 nT, 2.4 samples per second (Boris Vasiliev, JINR, Dubna, private communication).

・ The second stage of the research was carried out within the framework of the program “1998-2004 INRNE-BAS monitoring” and consisted of monitoring the geomagnetic field-Boris Vasiliev, JINR, Dubna, with a one component magnetometer.

These measurements were made with an interval of one minute and the obtained data allowed suggesting the following variables: the geomagnetic quake (GQ) of the field in the vicinity of a tidal wave maximum, the average daily magnitude of the magnetic field strength, and its hourly and minute average deviation [1] . This was the first statistical evidence that after a geomagnetic quake there is an increase of the seismic activity around a monitoring point estimated by the number of earthquakes with different magnitudes.

Further progress (2004-2010) was reached with the use of a three-component Danish fluxgate magnetometer, operated by the Skopje Seismological Observatory with the analysis of one minute INTERMAGNET data [48] [49] [50] .

Introduction of the diurnal sum of earthquake energy as a numerical estimation of regional seismic activity.

The distributions for variable Day Diff of 874532 world earthquakes with magnitude M = 3.5, 5, 6, 7 which occurred in the period of 1981-2017 (International Seismological Center: support the old hypothesis that the Sun-Moon tides are the earthquakes trigger. DayDiff is as follows:

DayDiff = EqTime − NearTideExtreme,

where EqTime is the time of earthquake and NearTideExtreme is the time of nearest extreme of the diurnal mean of the Sun-Moon tide values and is presented in Figure 1.

As is seen from the graphs, 92% of all analyzed earthquakes occurred in the time period of +/− 2.23 days around the time of a tide extreme in the locality of their epicenter.

The value at Day Diff can be interpreted as a count of aftershocks.

Table 2 illustrates the cases of magnitude (M), the values of sigma (half wide), total number of earthquakes, number of earthquakes outside of sigma intervals (+, −), and percent of earthquakes that occurred in the sigma interval.

It is important to note that the daily estimation of regional (around monitoring point) seismic activity can be achieved by summing up the variable reflecting

Figure 1. The Gauss fit of DayDiff distributions in case of different M.

Table 2. Distribution of earthquakes relative to tides maximum.

density of surface energy:

S Cht = 10 1.5 M + 4.8 / ( D + Depth + Dist Mp ) 2 [ J / km 2 ] (3),

where D = 40 km for the Balkan region, for every earthquake that occurred. Probably this value is determined by average lithosphere thickness in the Balkans. The explicit form of function S Cht was obtained using the Dubna inverse problem method [51] - [56] for estimation of the correlations between GmPrecSig and regional seismic activity in the tide’s extreme period.

The inverse problem solution allowed selection of the following variables based on monitoring data:

The estimation of the variability of one experimental series of measurements -

T = { T i , i = 1 , , n } includes the calculation of the mean value mT:

m T = i = 1 n T i n (4)

and the value of its dimensionless standard deviation SdT:

Sd T = i = 1 n ( 1 T i m T ) 2 n (5)

Obviously, the biggest standard deviation value means the biggest variability of the series.

Therefore, the procedure of monitoring includes numerical comparison of two sequential series of measurements.

The data used:

• Japan one minute INTERMAGNET geomagnetic stations MMB (Memambetsu, Lat 43.907˚N, Lon 144.193˚E, Altitude = 42 m), KAK (Kakioka, Lat 36.232˚N, Lon 140.186˚E, Altitude = 36 m) KNY (Kanoya, Lat 31.42˚N, Lon 130.88˚E, Altitude = 107 m) minute data (,

• Software for calculation of the daily and minute Earth tide behaviour (Dennis Milbert, NASA, Solid Earth Tide Software updated 07 June, 2018,,

• Earth tide extremes (daily average maximum, minimum and inflexed point) as a trigger of earthquakes,

• Data for World A-indices (National Oceanic and Atmospheric Administration Space Weather Prediction Center (NOAA SWPC)

The geomagnetic signal was calculated as a simple function of the relative standard deviations of the geomagnetic vector components. The precursor signal is the difference between today and yesterday’s geomagnetic signal corrected by the A- indices values. An increase of the precursor signal means an increase of the geomagnetic field variability. Such a positive leap was named a geomagnetic quake in analogy with an earthquake. An analysis of the correlation between the earthquakes occurrences, and the time of Sun-Moon Earth tide extremes on the basis of the variable earthquake’s surface energy density SChtM permits the forecast of an imminent regional seismic activity. Calculation of the day differences (DayDiff) between the time of the earthquakes occurrence and the time of the nearest tide extreme permits the building of the curve of DayDiff and its Gaussian fit. The comparison of Gauss widths for all the earthquake occurrences and those with the biggest SChtM permits the formulation of the hypothesis for “predictable” earthquakes.

The typical data and results of its analysis are shown in Figure 2

Figure 2. Kakioka diurnal geomagnetic and earthquakes monitoring in the time period around the Fukushima earthquake with geomagnetic field data on March 11, 2011.

The Geomagnetic field components Northm, Eastm, Downm, m = 1440, are the minute averaged values of the geomagnetic vector F, and the variables SdNorthh, SdEasth, dSDownh are their standard deviation, calculated for 1 hour, h = 1, … , 24:

SdNorth h = m = 1 60 ( 1 North m m North h ) 2 60 n (6)


m North h = m = 1 60 North m 60 (7)

The geomagnetic signal GeomHourSigh is the geometrical sum of hour standard deviation normed by the module of hour geomagnetic vector:

GmSig h = SdNorth h + SdEast h + SdDown h m North h + m East h + m Down h (8)

The A indices are the Low, Medium and High indices, calculated by the NOAA, Space Weather Prediction Center: In this paper we use ALow.

The variable GmSighDay is the diurnal mean value of GmSigh:

GmSig h Day = m = 1 60 GmSig h 60 (9)

and PrecSigDay is

PrecGmSig Day = 2 PrecGmSig Day PrecGmSig YesterDay A Day + A YesterDay (10)

The indices of earthquake’s magnitude value are the distance of one hundred km between the epicenter and the monitoring point.

The variable Periodic SChtSum [J/km2] is the sum of the variables SCht for all earthquakes that occurred in the time period +/− 2.7 days before and after the tide extreme in the 700 km region around the monitoring point. Obviously, its value can serve as an estimation of the regional seismic activities for the time period around the tide’s extreme.

The variable Diurnal SChtMSum [J/km2 per day] is the sum of the variable SCht, calculated for all earthquakes that occurred during the day in the 700 km region around the monitoring point. This variable can serve as a quantitative measure of diurnal regional seismicity.

It is worth noting that the explicit form of the variable SCht was established in the framework of the inverse problem [48] - [54] in an effort to get a clearer correlation between the variables PrecSigDay and PeriodicSChtSum.

The variable Tide Minute [cm] is the module of tide vector calculated every 15 minutes.

The variable Tide Day [cm] is the diurnal mean value in the time calculated in analogy of the mass center formulae

Time TideDay = m = 1 360 m TideDay m m = 1 360 TideDay m (11)

Note: For the seconds and number of samples per second, the generalization has been to calculate geomagnetic field characteristics for every minute, and correspondingly the values of GmSigday have to be the average for 1440 minutes.

The positive value of the variable PrecursorSigday means that the geomagnetic field variability, which is calculated via standard deviations of geomagnetic field components, is increasing (formulae 1, 2).

In an analogy with an earthquake we call such an increase a geomagnetic quake.

As one can see from Figure 2, after the appearance of the geomagnetic quake, in nine of twelve cases (75%), the regional seismic activity is increasing (the bigger value of the Periodic SChtMSum variable) in the time period around the following tide extreme. Therefore, the geomagnetic quake approach described herein can serve as a forecast of imminent regional seismic activity.

In Figure 3 the values of the variable Periodic SChtMSum are calculated not

Figure 3. Memambetsu diurnal geomagnetic and earthquakes monitoring for the period Jul. 1, 2014 ? Jan. 1, 2015.

only for the time periods around the extremes, but also for the time period between them. We can see that the values in almost every extreme period are higher.

The use of the above described analysis for a longer time period with the calculation of distribution of day difference between the “predicted” earthquakes (earthquakes with the highest value SChtM) demonstrates the reliability of this approach for forecasting imminent regional seismic activity for regions with seismic risk (Figure 4).

The correspondence between GmPrecSig and the predicted increase of regional seismic activity permits the formulation of the hypothesis that the “predicted” earthquake is defined as the quake with the bigger value of function SCht. The GQ analysis of one minute INTERMAGNET data of Kakioka, Memambetsu and Kanoya, Japan stations demonstrated that such earthquakes are “predicted” in the same time from two stations, which permits the formulation of the overdetermined algebraic system for the explicit form of function

GmPrecSig (Mi, Depti, Lati, Longi), i = 1, …, number of “predicted” earthquakes (12)

Using (12), with the number of monitoring points, number of reliable precursors one can formulate the conditions for existence of overdetermined algebraic

Figure 4. The distribution and its Gauss fit of DayDiff for all earthquakes occurred in Memambetsu (700 km radial distance) region [56] .

system, solution of which are the values of Mi, Depti, Lati, Longi), i = 1, …, number of predicted earthquakes in the time period on the next tide extreme.

At the end is a brief presentation of the real time data acquisition system for researching the reliability of earthquakes precursors, formulation and solution of the overdetermined inverse problem for the magnitude, depth, and coordinates of epicentre of impending earthquake.

The other precursor of earthquakes is the bore-hall water level. The preliminary results of such findings were positive (Georgia, Armenia, BlackSeaHazNet proceedings), but for its reliability, more statistical analysis is needed.

6. Conclusions

The approach proposed for solving the problem of the “how, where and when” of earthquake prediction does not decrease the importance of commonly accepted investigations based on seismology, geology, geoelectromagnetism and JPS data.

Presented herein is a new approach for forecasting of the regional weekly seismic and volcanic activity based on the analysis of the INTERMAGNET geomagnetic field and NASA code for Sun-Moon tides data and is a Geomagnetic Quake approach. The results are based on the inverse problem method for analysis of the geomagnetic field which instantaneously reflects terrestrial currents in hypocenter. The necessary and sufficient conditions for the existence of a solvable inverse problem are formulated on the basis of the existence of reliable precursors.

A real time data acquisition system for researching the reliability of earthquake precursors, formulation, and solution of the overdetermined inverse problem for the magnitude, depth, and coordinates of the epicentre of impending earthquake must include the following additional types of monitoring:

1) Increase of helium 3, 4 and radon surface concentration;

2) Increase of the local heat flows through the Earth’s surface;

3) The water-level in boreholes and steam release into atmosphere;

4) Change of the ionosphere height, conductivity for the radio waves and meteorological conditions;

5) Appearance of low frequency waves in the atmosphere and the Earth’s crusts which are not observed by nuclear test monitoring;

6) Variations in Schumann resonance.

An inverse problem could be solved for prediction time (±2.7 days), magnitude, depth and epicenter coordinates of an upcoming earthquake―with 4 parameters, which means that at least 4 monitoring points in a region of a radial distance of 700 km are necessary for formulation of the solvable overdetermined algebraic system.

The combination of the geomagnetic measurements and the above-listed additional reliable precursors is bound to make feasible an overdetermined algebraic system.

The solution of such a multi-parametric system will provide the possibility for estimating an earthquake’s magnitude and epifocal coordinates’ prediction accuracy and will be very useful in the further research of the nature of the tectonic processes.


One of us (SChtM) would like to thank my co-authors from BlackSeaHazNet FP7 IRSES (2011-2014) for cooperation and discussions, Boris Vasilev, JINR, Dubna for the one-component magnetometer, the mathematician and programmer Lubomir Aleksandrov for his code REGN and discussions on the inverse problem methods, Oleg Molchanov and Alexey Sissakyan for many years of fruitful polemics, as well as my late Prof. Vladimir Kadyshevsky, who taught me how to try to gradually understand the physics with analyzing the dynamics and kinemathics of the processes in different time, space and energy scales.


This paper is supported by the Project MON, FSI DN 14/1, from 11.12.2017, “Exploration of changes in some geophysical fields preceding the occurrence of earthquakes in the Balkans”.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Mavrodiev, S.C., Pekevski, L., Botev, E., Pinar, A., Kikuashvili, G., Vol, A. and Gilat, A. (2018) Study of the Possibility of Predicting Earthquakes. International Journal of Geosciences, 9, 688-706.


  1. 1. Mavrodiev, S.Cht. and Thanassoulas, C. (2001) Possible Correlation between Electromagnetic Earth Fields and Future Earthquakes. Seminar Proceedings, INRNE-BAS, Sofia, Bulgaria, 23-27 July 2001.

  2. 2. Reid, H.P. (1910) The California Earthquake of April 18, 1906; the Mechanics of the Earthquake (Carnegie Institute, Washington D.C., 1910).

  3. 3. Gilat, A. and Vol., A. (2005) Primordial Hydrogen-Helium Degassing, an Overlooked Major Energy Source for Internal Terrestrial Processes. HAIT Journal of Science and Engineering B, 2, 125-167.

  4. 4. Gilat, A. and Vol, A. (2012) Degassing of Primordial Hydrogen and Helium as the Major Energy Source for Internal Terrestrial Processes. Geoscience Frontiers, 3, 911-921.

  5. 5. Bolt, B.A. (1978) Earthquakes: A Primer. W.H. Freeman and Company, San Francisco, 112-115.

  6. 6. Richter, C.F. (1958) Elementary Seismology. W.F. Freeman and Company, San Francisco, and Bailey Bros. & Swinfen Ltd., London

  7. 7. Bürgmann, R. and Dresen, G. (2008) Rheology of the Lower Crust and Upper Mantle: Evidence from Rock Mechanics, Geodesy, and Field Observations. Annual Review of Earth and Planetary Sciences, 36, 531-567.

  8. 8. Shrestha, B. (2009) Vertical Ground Motions and Its Effect on Engineering Structures: A State-of-the-Art Review.

  9. 9. Bozorgnia, Y. Vertical Ground Motion.

  10. 10. Semenov, N.N. (1956) Nobel Price Lecture: “Some Problems Relating to Chain Reactions and to the Theory of Combustion. Nobel Lecture.

  11. 11. Dahy, S.A. and Hassib, G.H. (2009) Discriminating Nuclear Explosions from Earthquakes at Teleseismic Distances. European Journal of Applied Sciences, 1, 47-52.

  12. 12. Gieras, M., et al. (2006) Determination of Explosion Parameters of Methane-Air Mixtures in the Chamber of 40 dm3 at Normal and Elevated Temperature. Journal of Loss Prevention in the Process Industries, 19, 263-270.

  13. 13. Zhang, Q., et al. (2012) Effect of Spark Duration on Explosion Parameters of Methane/Air Mixtures in Closed Vessels. Safety Science, 50, 1715-1721.

  14. 14. Maranda, A. and Szymański, R. (2013) Tests on Critical Diameter and Detonation Velocity of Mixtures of Ammonium Nitrate (V) and Selected Organic Substances. Chemik, 67, 13-18.

  15. 15. Antillon, E. and Strachan, A. (2015) Mesoscale Simulations of Shockwave Energy Dissipation via Chemical Reactions.

  16. 16. Fukuhara, M. (2016) Possible Generation of Heat from Nuclear Fusion in Earth’s Inner Core. Scientific Reports, 6, Article No. 37740.

  17. 17. Ratis, Y.L. (2009) Neutrino Catalysis of Nuclear Synthesis Reactions in Cold Hydrogen. Concepts of Physics, 4, 525-541.

  18. 18. Rofman, V.M. (2012) To the Question of Genesis of Tritium in the Water of River Chagan (Semipalatinsk Anomaly, Kazachstan,). (In Russian)

  19. 19. Shomer, I. (2010) Is Deuterium Fusion Catalyzed by Antineutrinos?

  20. 20. Prados-Estévez, F., Subashiev, A. and Nee, H. (2017) Nuclear Fusion by Lattice Confinement. Journal of the Physical Society of Japan, 86, Article ID: 074201.

  21. 21. Prados-Estévez, F., Subashiev, A. and Nee, H. (2017) Strong Screening by Lattice Confinement and Resultant Fusion Reaction Rates in Fcc Metals. Nuclear Instruments and Methods in Physics Research Section B, 407, 6772.

  22. 22. Terez, E. and Terez, I. (2015) Fusion Reactions as the Main Source of the Earth’s Internal Energy. Herald of the Russian Academy of Sciences, 85, 163-169.

  23. 23. Storms, E. (2012) An Explanation of Low-Energy Nuclear Reactions (Cold Fusion). The Journal of Condensed Matter Nuclear Science Electronic, 9, 86-107.

  24. 24. Meijer de, R.J. and van Westrenen, W. (2008) Assessing the Feasibility and Consequences of Nuclear Georeactors in Earth Core-Mantle Boundary Region. South African Journal of Science, 104, 111-118.

  25. 25. Araki, T., et al. (2005) Experimental Investigation of Geologically Produced Antineutrinos with KamLAND. Nature, 436, 499-503.

  26. 26. Gando, A. (2011) The KamLAND Collaboration. Partial Radiogenic Heat Model for Earth Revealed by Geoneutrino Measurements. Nature Geoscience, 4, 647-651.

  27. 27. Rusov, V.D., et al. (2011) Solar Dynamo as Host Power Pacemaker of Earth Global Climate, EU FP7 IRSES 2011 Project. Complex Research of Earthquake’s Forecasting Possibilities, Seismicity and Climate Change Correlations.

  28. 28. Běhounková, M., et al. (2009) Coupling Mantle Convection and Tidal Dissipation: Applications to Enceladus and Earth-Like Planets. Journal of Geophysical Research, 115, e09011.

  29. 29. Toulhoat, H., et al. (2012) Chemical Differentiation of Planets: A Core Issue.

  30. 30. Recanati, A., Kurz, M.D., Warren, J.M. and Curtice, J. (2012) Helium Distribution in a Mantle Shear Zone from the Josephine Peridotite. Earth and Planetary Science Letters, 359-360, 162.

  31. 31. Subashiev, A. and Nee, H. (2017) Hydrogen Trapping at Divacancies and Impurity-Vacancy Complexes in Nickel: First Principles Study. Journal of Nuclear Materials, 487, 135-142.

  32. 32. Jaupart, C., Labrosse, S. and Mareschal, J.S. (2015) Temperatures, Heat and Energy in the Mantle of the Earth.

  33. 33. Gufeld, I.L. and Matveeva, M.I. (2011) Barrier Effect of Degassing and Destruction of the Earth’s Crust. Doklady Earth Sciences, 438, 677-680.

  34. 34. Vol, A. (2014) Thermo-Electrochemical Processes of the Earth’s Degassing Creating Geomagnetic Field and Changing Its Value and Direction (Thermodynamic Approach). International Journal of Geosciences, 5, 1219-1230.

  35. 35. Cowen, R. (2013) Two-Laser Boron Fusion Lights the Way to Radiation-Free Energy. Nature.

  36. 36. Hora, H., Eliezer, S., Kirchhoff, G.J., Nissim, N., Wang, J.X., Lalousis, P., Xu, Y.X., Miley, G.H. and Martinez-Val, J.M. (2017) Road Map to Clean Energy Using Laser Beam Ignition of Boron-Hydrogen Fusion. Laser and Particle Beams, 35, 730-740.

  37. 37. Galembeck, F., et al. (2014) Friction, Tribochemistry and Triboelectricity: Recent Progress and Perpectives.

  38. 38. Ray, R.D. and Egbert, G.D. (2012) Fortnightly Earth Rotation, Ocean Tides and Mantle Anelasticity. Geophysical Journal International, 189, 400-413.

  39. 39. Venedikov, A. and Arnoso, R. (2001) Program VAV/2000 for Tidal Analysis of Unevenly Spaced Data with Irregular Drift and Colored Noise. Journal of the Geodetic Society of Japan, 47, 281-286.

  40. 40. Venedikov, A.P., Arnoso, R. and Vieira, R. (2003) A Program for Tidal Data Processing. Computer & Geosciences, 29, 487-502.

  41. 41. Losada-Pérez, P., Tripathi, C.S.P. and Glorieux, C. (2011) Measurements of Heat Capacity and Enthalpy of Phase Change Materials by Adiabatic Scanning Calorimetry. International Journal of Thermophysics, 32, 913-924.

  42. 42. Yi, L., Jin, C. and Liu, J. (2013) Liquid-Solid Phase Transition Alloy as Reversible and Rapid Molding Bone Cement.

  43. 43. Hayashi, Y., et al. (2007) Phase Transformation Behaviour and Bending Properties of Hybrid Nickel-Titanium Rotary Endodontic Instruments. International Endodontic Journal, 40, 247-253.

  44. 44. Kirby, S.H. (1987) Localized Polymorphic Phase Transformations in High-Pressure Faults and Applications to the Physical Mechanism of Deep Earthquakes. Journal of Geophysical Research. Solid Earth, 92, 13789-13800.

  45. 45. Knopoff, L. (1964) Earth Tidevs as a Triggering Mechanism for Earthquakes. Bulletin of the Seismological Society of America, 54, 1865-1870.

  46. 46. Tamrazyan, G.P. (1967) Tide-Forming Forces and Earthquakes. ICARUS, 7, 59-65.

  47. 47. Tamrazyan, G.P. (1968) Principal Ragularities in the Distribution of Major Earthquakes Relative to Solar and Lubar tades and Other Cosmic Forces. ICARUS, 9, 574-592.

  48. 48. Mavrodiev, S.C. (2004) On the Reliability of the Geomagnetic Quake as a Short Time Earthquake’s Precursor for the Sofia Region. Natural Hazards and Earth System Sciences, 4, 433-447.

  49. 49. Mavrodiev, S.C. and Pekevski, L. (2008) Complex Regional Network for Earthquake Researching and Imminent Prediction. In: Singh, B., Ed., Electromagnetic Phenomena Related to Earthquakes and Volcanoes, Narosa Pub. House, New Delhi, 135-146.

  50. 50. Mavrodiev, S.C., Pekevski, L. and Jimseladze, T. (2008) Geomagnetic-Quake as Imminent Reliable Earthquake’s Precursor: Starring Point for Future Complex Regional Network. In: Singh, B., Ed., Electromagnetic Phenomena Related to Earthquakes and Volcanoes, Narosa Pub. House, New Delhi, 116-134.

  51. 51. Alexandrov, L. (1970) Auto-Regularized Gauss-Newton-Kantorovich Iteration Process. Comm. JINR, P5-5515, Dubna.

  52. 52. Alexandrov, L. (1971) Regularized Computational Process of Newton-Kantorovich Type. Computational Mathematics and Mathematical Physics, 1, 36-43.

  53. 53. Alexandrov, L. (1973) The Program REGN (Regularized Gauss-Newton Iteration Method) for Solving Nonlinear Systems of Equations. Comm. JINR P5-7259, Dubna.

  54. 54. Alexandrov, L. (1983) Program Code REGN, PSR 165 RSIK ORNL, Oak Ridge, Tennessee, USA.

  55. 55. Andrault, D., Monteux, J., Le Bars, M. and Samuel, H. (2016) The Deep Earth May Not Be Cooling Down. Earth and Planetary Science Letters, 443, 195-203.

  56. 56. Mavrodiev, S.C. (2016) Imminent Earthquake Forecasting on the Basis of Japan Intermagnet Stations, NEIC, NOAA and Tide Code Data Analysis. Open Journal of Earthquake Research, 5, 62-78.