=Paper=
{{Paper
|id=Vol-2608/paper66
|storemode=property
|title=Neural network forecasting of earth globe seismic activity level
|pdfUrl=https://ceur-ws.org/Vol-2608/paper66.pdf
|volume=Vol-2608
|authors=Vadym Tiutiunyk,Leonid Chernogor,Olha Tiutiunyk,Tural Agazade
|dblpUrl=https://dblp.org/rec/conf/cmis/TiutiunykCTA20
}}
==Neural network forecasting of earth globe seismic activity level==
Neural Network Forecasting
of Earth Globe Seismic Activity Level
Vadym Tiutiunyk1[0000-0001-5394-6367], Leonid Chernogor2[0000-0001-5777-2392],
Olha Tiutiunyk3[0000-0002-3330-8920], Tural Agazade1[0000-0002-0165-1118]
1
National University of Civil Defence of Ukraine,
Chernyshevska Str., 94, Kharkiv, 61023, Ukraine
tutunik_v@ukr.net, agazade.tural.2019@gmail.com
2
V.N. Karazin Kharkiv National University,
Svobody Sq., 4, Kharkiv, 61022, Ukraine
leonid.f.chernogor@univer.kharkov.ua
3
Simon Kuznets Kharkiv National University of Economics,
Science Ave., 9-А, Kharkiv, 61166, Ukraine
tutunik_o@ukr.net
Abstract. In order to develop the scientific and technical foundations of creat-
ing an artificial intelligence system for monitoring tectonic origin emergencies,
the paper in the results of modeling and forecasting conducted on the basis of
neural network technologies, total number of occurrences, with discretion in
one month, earthquakes on Earth as one of the main characteristic indicators of
the seismic activity of the Earth’s globe are presented.
Keywords. emergency, earthquake, tectonic source emergencies monitoring,
artificial intelligence system, artificial neural network, neural network simula-
tion, seismic activity level prediction
1 Introduction
The tendency to an abrupt increase in the number and destructive power of natural
disasters over the past few decades of the life of society leads to a deterioration of
socio-economic and environmental consequences. It indicates the need to develop
effective measures to prevent and eliminate emergencies (ES) of various nature on the
Globe [1–3].
A promising direction for solving this problem is the development of an effective
hazard detection system at the stage of their inception. Also, the causes will be estab-
lishing of the occurrence these factors manifestations and effects on them in order to
prevent the occurrence of emergencies. This has been implemented on the basis of the
classical control loop presented in Fig. 1 [4–7].
This article is part of a planned set of scientific studies aimed at developing a
safety system. This eliminates or minimizes losses as much as possible under condi-
tions of manifestation of an emergency. The work is focused on studying the proc-
Copyright © 2020 for this paper by its authors. Use permitted under Creative
Commons License Attribution 4.0 International (CC BY 4.0).
esses of emergence and spreading of emergencies of lithospheric origin, which repre-
sent or may pose a serious danger to the life of society [8–12].
DECISION SUPPORT
SYSTEM The person making
the decision
SYSTEM
Determination of PERFORMANCE
CHOICE
parameter SOLUTIONS
Management decisions
GEOGRAPHIC INFORMATION Modeling and Formalization of
SYSTEM
for the prevention and
prediction of management
elimination of
development of decisions
emergency emergencies
Analysis and
DataBase of Bringing a decision
knowledge
systematization to the performers
of information
CLASSICAL MONITORING Data processing
EMERGENCY SYSTEM
ES
Hazard
Influence on
registration
the hazard
LOCAL TERRITORY
Fig. 1. Diagram of the emergencies monitoring structure as a means of control
(Source: author's own elaboration)
When solving the problem of creating an artificial intelligence system for monitor-
ing emergencies of the tectonic nature, there is a need to studying and simulation the
processes of occurrence and propagation of the seismic activity level of the local terri-
tory under the conditions of seismic activity of the Earth’s globe, as an element of the
system of nonlinear energy interactions Sun – Earth – Moon.
2 Formal problem statement
Development of scientific and technical bases for creation of artificial intelligence
system for monitoring of tectonic origin of the emergency, which is realized by con-
structing neural network models of forecasting the level of seismic danger of the
Earth’s territory by the amount and destructive energy of the emergency of tectonic
origin.
The construction of neural network models for predicting the level of seismic dan-
ger in the Earth’s globe was performed using the STATISTICA statistical package
[13, 14].
3 Literature review
Today, the basis of predictive observations of Earth’s seismic activity is knowledge
of the physical laws of the earthquake mechanism and control of the physical fields in
the seismic zone [15–19]. The basis of these studies is the idea of the absence of ex-
ternal factors affecting the study area that arise in the system of nonlinear energy
interactions Sun – Earth – Moon.
The mechanistic approach, which has been developing in various countries for
quite some time, has provided answers to many questions about the preparation of the
earthquake. Currently, more than a dozen earthquake models have been developed,
the most famous of which is following papers [20–22]:
a model of avalanche-unstable crack formation (AUC), which consists in the rapid
growth of the number of cracks, their interaction with each other and eventually the
occurrence of a main or main rupture, the shift of which instantly drops the accumu-
lated elastic with the formation elastic waves;
dilatant diffusion model (DD), in which crack formation also plays a decisive role,
but unlike the AUC model, fractures occur and the presence of water in the rocks of
the epicentral region is essential;
consolidation model of earthquake preparation (Dobrovolskiy’s model) describes
earthquake preparation as a process of emergence and evolution of rigid inhomoge-
neities in a continuous environment, as well as some other models that emerged as a
generalization of the study of stress and fracture processes of solids samples in the
laboratory and their subsequent transfer to natural seismic events.
At the same time, the following methods are now used to obtain a comparative as-
sessment of the level of life-threat in the conditions of manifestation of a emergency
[23–30]:
statistical, based on the analysis of emergency statistics on local territories over
several years to determine the risk indicators;
probable, based on the application of mathematical models, which connects the
prerequisites for the appearance of emergency with the possibility of their manifesta-
tion;
expert based on expert judgment in combination with fuzzy set theory.
The advantage of the statistical method is objectivity. Probable and expert methods
allow to take into account the sources of potential danger, which are rare in the form
of emergencies, but the consequences of which are catastrophic. However, the prob-
able method is extremely cumbersome and time-consuming, requiring a large number
of outputs, which results in low accuracy of the results obtained. In the absence of
tried and tested mathematical models and reliable enough initial data for them, it is
advisable to carry out an expert evaluation of the impact on the possibility of imple-
menting large-scale emergencies of a large number of difficultly formalized initial
data.
The prospects of using a statistically probabilistic approach to predict earthquakes
in a separate Earth globe, without taking into account the external factors (as an ele-
ment of the Sun – Earth – Moon nonlinear energy interaction system), are given
in [31]. This paper presents the results of a study of earthquake prediction in the
northwestern area of Vietnam using neural network technologies.
Thus, the use of neural network technologies is one of the promising directions for the
development of approaches to predicting earthquakes across the globe into systems of
nonlinear Sun – Earth – Moon interactions. This determines the direction of our scientific
research in the field of monitoring of the tectonic origin emergencies [8–10, 32–34].
4 Particularities of the processes occurring in the Sun – Earth –
Moon system and affecting the level of seismic hazard
of the Earth’s globe functioning
The dynamics of the physical processes of the Sun – Earth – Moon system affect-
ing the seismic hazard level of the local territory functioning can be schematically
represented according to Fig. 2 – 4. This can be characterized by the following spatial
constructions within the solar galaxy.
1. The axis of rotation of the Earth in the celestial sphere describes a complex
wave-like trajectory. The points of the axis of rotation are at an angular distance of
about from the pole ecliptic (Fig. 2). The vertex of the cone coincides with the Earth
center. The points of equinoxes and solstices move along the ecliptic towards the sun.
Moments of gravitational forces influence on the equatorial bulges and vary depend-
ing on the positions of the Moon and the Sun relative to the Earth. When the Moon
and the Sun are in the plane of the Earth’s equator the moments of forces disappear. If
tilts of Moon and Sun are the maximum, then the magnitude of the torque will be
greatest. The nutations, owing to fluctuations in the moments of the forces of the axis
of rotation of the Earth have been observed by consist of a series of small periodic
oscillations. The main nutations have a period of 18.6 years - the time of the orbital
nodes of the Moon. Movement with this period occurs on an ellipse. The major axis
of the ellipse is perpendicular to the direction of the precessional motion and is equal
to; small - parallel to it and equal. Next in magnitude of the amplitude are the compo-
nents with a period of 0.5 year, 13.7 days., 9.3 years, 1 year, 27.6 days. etc., therefore
the trajectory has the form of “thin laces” (shown on the enlarged fragment in the left
part of Fig. 2) [35–38].
2. The pressure from the solid inner core and the surrounding melt (outer core)
onto the mantle arises as a result of the eccentric revolution of the Earth’s shell
around the displaced inner core, which squeezes the shell from the inside. The forces
compressing the shell of the sphere (planet) and drawing it inward to the core arise in
other parts of the planet. This process has two components: impact at the expense of
the annual displacement of the center nucleus relative to the center of the globe
(Fig. 2 – 4); impact at the expense of eccentric circulation of the core relative to the
lower mantle, when due to the difference in angular rotation velocity of the core and
lower mantle ( 1 – angular velocity of rotation of the mantle; 2 – angular velocity
of rotation of the outer core; 3 – angular velocity of rotation of the inner core;
2 1 – angular velocity of rotation of the outer core relative to mantle
(“western drift”)), therefore, there are zones of high pressure and vacuum
( P1 P2 where P1 and P2 are indicators of pressure of the inner core of the globe on
its surface), affecting the level of seismic activity of the surface of the Earth (Fig. 3).
As long as there is a difference in the angular velocity of rotation and displacement of
the nucleus, the appearance of such zones will be maintained [39, 40].
18,4 PRECESSION
n
NUTATION
P
PRECESSION AND NUTATION Axis of rotation
23,50
Axis tilt of orbital curve
1 ; 2 ; 3
Perpendicular to the orbit
Axis rotation
SUN
EARTH N
Orbit direction
Planetary
equator
MOON Moon center
Ecliptic The motion of the center of the
Earth's outer Earth (inner core) around the center
core of mass
S The center of mass in the Earth-Moon
Earth's inner
core Center of the Earth
Fig. 2. Motion diagram of the inner core of the Earth in the Sun–Earth–Moon sys-
tem (Source: author's own elaboration)
P1
Inner core
Outer core
3
ES
2 1
P2
Fig. 3. Influence diagram of internal core oscillations on seismic activity (Source:
author's own elaboration)
3. Internal elastic stresses arise in the process of moving lithospheric plates
(Fig. 4), which are energy sources of earthquakes. The occurrence depth of elastic
stresses depends on the nature of a movement plates. The relative motion of litho-
spheric plates leads to the emergence of shallow (not deeper than 20–25 km) earth-
quake sources and dipping of lithospheric plates into the mantle initiates the appear-
ance of sources of deep (exceeding 70 km) earthquakes. The probability of elastic
stresses - sources of earthquakes decreases with increasing distance from the interface
of the lithospheric separation plates.
4. Surface and bulk seismic waves are the propagation factors of earthquake haz-
ards Z0 , that can cause secondary earthquakes [41, 42].
5. The probability of mutual amplification or weakening of bulk seismic waves in-
creases in the process of spatial-vibrational movement of the Earth’s internal core and
its effect on the external core. Consequently, the possibility of secondary earthquakes
Z increases also [43].
7. The territorial-temporal changes in the intensity of the natural electromagnetic
field pulses of the Earth initiating anomalous processes in the atmosphere occur due
to the movement of the Earth’s inner core has been established [44–46].
Directions of Earthquake
movement probability
lithospheric
plates (arrows)
Earthquake
Elastic stresses Surface seismic
waves
x
Z0
Z’ Z’’
Plate І Plate ІІ
1
Lithospheric
plate section
Bulk seismic Mantle 2
waves 3
Outer core Inner core
Fig. 4. Process diagram of earthquake and the spread of seismic activity (Source:
author's own elaboration)
Thus, combining the analysis results of the impact dynamics and energy of the in-
ternal physicochemical processes of the Earth on the origin generating tectonic proc-
esses allowed to formulate an approach to studying the nature of seismic phenomen. It
is an important tool for analyzing the results of civil defense research on the devel-
opment of models for the development of ES tectonic nature.
5 Experiments and results of neural network prediction of the
level of seismic activity of the Earth
According to the purpose of the research, the solution of the scientific problem in
the work is ensured by constructing an artificial neural network (ANN) model – a
model for predicting the level of seismic activity of the planet from the conditions of
the system functioning of nonlinear energy interactions Sun – Earth – Moon. This is a
time series model, where the initial indicator, according to Fig. 5, is the total number
of earthquakes in the world N t (provided that M 5 the magnitude of the earth-
quake), depending on the current time of analysis ( t ), as well as the distance of the
Earth’s inner core from the center of the planet t and the change in the length of
r0
day LODt St 86400s (where St 86400s – the length of day,
r t
r0 7,292115 10 5 rad s – constant (average) angular velocity of the Earth’s own
rotation).
For the development presented in Fig. 5 models have a multilayer perception
(MLP) selected. The input parameters of this model are the results of the analysis of
the monthly dynamics of indicators, N t , t and LODt for 2009 – 2018,
which are presented in [8–10, 32–34] and Fig. 6.
Year
January
February
March
April
Artificial neural network
May
t June
July N t
August
September
October
November
December
* t
LOD* t
Fig. 5. Scheme of the ANN time series model for the monthly forecasting of the
Earth’s number of earthquakes, depending on the indices of the remoteness of the
Earth’s inner core from the center of the planet and changes in the length of day
(Source: author's own elaboration)
Fig. 6. Monthly dynamics of the variations in the rate of variation of the Earth’s axial
rotation ( LODt ), the distance of the inner core from the center of the Earth ( t )
and the total seismic activity ( N t ) with magnitude during 2009 – 2018 are sum-
marized on the Globe (Source: author's own elaboration)
When developing a mathematical model for predicting the level of seismic activity
of the Earth ball, based on the decision parameter N t on the complexity of the
MLP architecture, it was based on five analyzes of the learning results of networks,
which accidentally included five hundred neural networks. The criterion for choosing
the optimal network was the relationship between the forecast error and the complex-
ity of the architecture. The results of the analysis are presented in Table 1.
Table 1. Statistical characteristics of the three-layer perceptions suggested by the
“decision-maker” as the best for neural network time series analysis of the forecast of
the monthly dynamics of the Earth’s seismic activity ( N t )
Study Control Test Activation
Activation
№ Architecture produ- produ- produ- hidden
of yield
ctivity ctivity ctivity layer
1. MLP 15-199-1 0,944 0,859 0,976 Exp. Logistic
2. MLP 15-233-1 0,935 0,634 0,964 Identity Exp.
3. MLP 15-243-1 0,939 0,859 0,986 Exp. Logistic
4. MLP 15-187-1 0,938 0,884 0,985 Exp. Logistic
5. MLP 15-12-1 0,967 0,727 0,668 Tanh Identity
For neural network training, all observations were divided into three samples. By
default, random sampling was performed between samples to avoid retraining the
network and to ensure quality generalization (forecasting). The first sample (Educa-
tional – 50% of observations) was used to train the network; the second (Control –
25% of observations) – for cross-validation of the learning algorithm during its opera-
tion; third (Test – 25% of observations) – for the final independent testing of a trained
neural network. The training was done with speed 0,01 .
Presented in Table. 1 networks are characterized by a more effective balance be-
tween modeling error and architecture complexity for time series analysis of the fore-
cast of the monthly dynamics of Earth’s seismic activity ( N t ). This was the basis
for further construction of the three-layer MLP 15-12-1 network, which has fifteen
inputs (Fig. 5), 12 elements in the hidden layer and one logical output function of
thirteen inputs, presented in Fig. 7.
Fig. 7. MLP 15-12-1 Three-Layer Perceptron Architecture with Logical Signal
Transmission for Time Series Prediction of Globe Seismic Monthly Dynamics
( N t ) (Source: author's own elaboration)
In this case, the use of logical activation functions, with scaling parameters, was
based on a given fraction of span of a logical function equal 0,9 (respectively range
0,05; 0,95 ) to the range of training of the neural network. The function of activating
the hidden layer of the perceptron MLP 13-10-1 is hyperbolic tangent. This allows for
a slight extrapolation of the data. In addition, the use of logical functions stabilizes
learning.
The results of checking the adequate prognostic performance of MLP 15-12-1 net-
work are presented in Fig. 8, where the dependence observed ( N* t ) statistic the
values of the Earth’s seismic activity indicator via predicted ( N t ) values. The
coefficient of correlation between these indicators on the results of the training "Edu-
cational Choice" is equal rN2 * t N t 0,967 .
Fig. 8. Dependence of observed ( N* t ) statistic values of the Earth’s seismic activ-
ity indicator by via predicted values ( N t ) by MLP Network 15-12-1 (Source:
author's own elaboration)
Fig. 9. Graph of the time dependence (of observation number) where the depend-
ence observed ( N* t ) statistic the values of the Earth’s seismic activity indicator via
predicted ( N t ) values by MLP network 15-12-1 (Source: author's own elabora-
tion)
The graph of time dependence (from observation number), with the appropriate
level of prognostic adequacy, where the dependence observed ( N* t ) statistic the
values of the Earth’s seismic activity indicator via predicted ( N t ) values by the
MLP 15-12-1 network is presented in Fig. 9.
According to the data analysis Table 2 it is necessary to state that obtained within
the limits of the ideas about the dynamics of physical processes that occur in the sys-
tem Sun – Earth – Moon and affect the level of seismic danger of functioning of the
local territory of the planet (see Fig. 2 – 4), the MLP neural network 15-12-1 and the
results of its forecast allow us to ascertain the adequacy, in accordance with the data
of Figs. 8 and the rN2 * t N t 0,967 parameter figure presented in Fig. 5 models for
predicting the level of seismic activity of the globe.
Table 2. The result of the prediction by the MLP network 15-12-1 level of seismic
activity-those of the Globe
2019
t
01 02 03 04 05 06 07 08 09 10 11 12
t
*
3,00 4,002,50 -1,00 -2,00 -3,00 -4,50 -5,00 -3,00 -0,50 0,50 2,00
LOD t
*
0,81 1,00 0,96 0,98 0,70 0,19 0,01 0,05 0,17 0,27 0,37 0,24
N* t 120 104 118 113 93 103 117 127 131 137 138 123
N t 116 113 132 120 83 97 125 141 144 124 147 138
N t ,% 3 9 11 6 10 6 7 11 9 9 7 12
6 Conclusion
The paper presents a neural network model of the time series for the monthly fore-
casting of earthquakes with magnitude М 5 . The obtained neural network allows us
to predict the level of seismic activity of the globe with adequacy at the
level rN2 * t N t 0,967 .
The results obtained are the basis for the development of scientific and technical
bases for the creation of a system of artificial intelligence for the performance of tasks
of monitoring of tectonic origin.
References
1. First Report of the Chairman of the International Committee on Problems of Global
Changes in the Geological Environment „GEOCHANGE”. http://www.ru.geochnge-
report.org
2. Barishpolets, V.A.: System analysis of disasters occurring in the world. Radio Electronics.
Nanosystems. Information Technology, 2, 1 – 2, 162 – 176 (2010).
3. National report on the state of technogenic and natural security in Ukraine.
http://www.dsns.gov.ua/
4. Code of Civil Protection of Ukraine. In: Voice of Ukraine, 220 (5470), 4 – 20 (2012).
5. Resolution of the Cabinet of Ministers of Ukraine of January 9, 2014, No. 11 “On Ap-
proval of the Regulation on the Unified State Civil Protection System”.
http://zakon5.rada.gov.ua/laws/show/11-2014-%D0%BF
6. Kalugin, V.D., Tiutiunyk, V.V., Chernogor, L.F., Shevchenko, R.I.: Development of scien-
tific and technical basis for establishment of monitoring, prevention and liquidation of
emergency situations of natural and man-made nature, and also ensuring of environmental
of ecological security. Information processing systems, 9(116), 204 – 216 (2013).
7. Andronov V.A., Divizinyuk, M.M., Kalugin, V.D., Tiutiunyk V.V.: Scientific and design
bases of complex creation system of emergency monitoring situations in Ukraine. National
university of civil protection of Ukraine (2016).
8. Tiutiunyk, V., Chernogor, L. V., Kalugin, Agazade, T.: Development of the bases of
geoinformational the systems of emergency monitoring situations of tectonic origin. Ap-
plied Radio Electronics, 18, 1,2, 52 – 65 (2019).
9. Tiutiunyk, V., Kalugin, V., Pysklakova, O., Yaschenko, O., Agazade, T.: Hierarchical
clustering of seismic activity local territories Globe. EUREKA: Physics and Engineering,
4: 41 – 53 (2019).
10. Tiutiunyk, V., Chernogor, L. V., Kalugin, Agazade, T.: Functional description of the zon-
ing of local territories of the Globe by quantity and destructive energy of tectonic extreme
origin situations. Municipal economy of cities, 1, 154: 272 – 287 (2020).
11. Chernogor, L.F.: Geomagnetic Disturbances Accompanying the Great Japanese Earth-
quake of March 11, 2011. Geomagnetism and Aeronomy, 59, 1, 62 – 75 (2019).
12. Chernogor, L.F.: Possible Generation of Quasi-Periodic Magnetic Precursors of Earth-
quakes. Geomagnetism and Aeronomy, 59, 3, 374 – 382 (2019).
13. Kim, J.O., Muller, C.U., Klekka, U.R.: Factor, discriminant and cluster analysis. Finance
and Statistics (1989).
14. Khalafyan, A.A.: STATISTICA 6. Statistical Data Analysis. LLC Binom-Press (2007).
15. Sobolev, G.A., Ponomarev, A.V.: Earthquake Physics and Harbingers. Science (2003).
16. Guglielmi, A.V.: Foreshocks and aftershocks of strong earthquakes in the light of disaster
theory, Advances in Physical Sciences, 185, 4, 415 – 429 (2015).
17. B.V. Levin, M.V. Rodkin, I.N. Tikhonov, Great Japanese earthquake, Nature, 10: 14 – 22,
2011.
18. Tertyshnikov, A.V., Platonov, V.V.: Prospects for monitoring seismic conditions from
space. Electronic scientific journal "Investigated in Russia.
http://www.zhurnal.ape.relarn.ru/artecles/2007/031.pdf
19. Smirnov, V.M., Smirnova, E.V.: Study of the possibility of using satellite navigation sys-
tems for monitoring seismic phenomena. Electromechanical issues, 105, 94 – 104 (2008).
20. Romashov, A.N., Gypsy, S.S.: In search of a generalized geotectonic concept. Geotecton-
ics, 4, 3 – 12 (1996).
21. Gufeld, I.L.: Seismic process. Physicochemical aspects. TsNIIMash (2007).
22. Mishin, S.V., Panfilov, A.A., Hasanov, I.M.: Earth's gravity is the cause of earthquakes.
Geophysical Journal, 41, 6, 213 – 222 (2019).
23. Kotovenko, O.A., Sobolevska, L.I., Miroshnichenko, O.Yu.: Stochastic modeling in the
study of processes under the influence of environmental management in the region. East
European Journal of Advanced Technology, 2/14, 37 – 41 (2012).
24. Akimov V.A., Radiev, N.N.: Determination of the relative hazard of the territories. Emer-
gency Security Issues, 6, 129 – 140 (2000).
25. Lepikhin, A.M.: Integrated Territory Security Indicators. Emergency Security Issues, 5,
93 – 98 (2008).
26. Vladimirov, V.A., Kulba, V.V., Malinetsky, G.G., Makhutov, H.A.: Risk management.
Science (2000).
27. Reinschke, K., Ushakov, I.A.: Graph reliability assessment using graphs. Radio and com-
munication (1988).
28. Ryabinin, I.A.: Reliability and safety of structurally complex systems. Polytechnic (2000).
29. Ryabinin, I.A.: The concept of logical probability theory of security. Devices and control
systems, 10, 6 – 9 (1993).
30. Y.D. Vishnyakov, N.N. Radaev. General risk theory, Publishing Center "Academy", 2008.
31. Pupkov, K.A., Chong, Khao Dinh: The use of neural network technologies in earthquake
prediction problems (for example, the north-western region of Vietnam). Bulletin of the
N.E. Bauman Moscow State Technical University, Series "Instrument Making", 2, 70 – 78
(2012).
32. Tiutiunyk, V.V., Chernogor, L.F., Kalugin, V.D., Agazade, T.: Assessment of influence
power effects in system the Sun-Earth-Moon on the level of seismic activities the territory
of the Globe. Management, navigation and communication systems, 6(46), 238 – 246
(2017).
33. Tiutiunyk, V.V., Chernogor, L.F., Kalugin, V.D., Agazade, T.: Assessment dynamic and
power effects on the Earth and their influences on ratios between seismic activity levels of
Globe hemispheres. Scientific Bulletin: Civil protection and fire safety, 2(4), 101 – 117
(2017).
34. Tiutiunyk, V.V., Chernogor, L.F., Kalugin, V.D., Agazade, T.: Influence estimation of ro-
tation speed variation of the Earth on level of seismic activity of local Globe territory.
GEOINFORMATIKA, 3(67), 36 – 48 (2018).
35. Sidorin, A.Ya.: Influence of the Sun on seismicity and seismic noise. Seismic instruments,
40, 71 – 80 (2004).
36. Levin, B.V., Sasorova, E.V., Domansky, A.V.: Properties of “critical latitudes”, rotation
variations and seismicity of the Earth. Bulletin of the Far Eastern Branch of the Russian
Academy of Sciences, 3, 3 – 8 (2013).
37. Wiemer, S., Wyss, M.: Mapping spatial variability of the frequency-magnitude distribution
of earthquakes: An overview. Advances in Geophysics, 45, 259 – 302 (2002).
38. Atef, A.H., Liu, K.H., Gao, S.S.: Apparent weekly and daily earthquake periodicities in the
Western United States. Bull. Seismol. Soc. Amer., 99, 4, 2273 – 2279 (2009).
39. Klimenko, A.V.: Global properties of Earth's seismic activity and their relationship with its
rotation. The dissertation ... of the candidate Phys.-Math. of sciences (2005).
40. Malyshkov, Yu.P., Malyshkov, S.Yu.: Periodic variations of geophysical fields and seis-
micity, their possible connection with the motion of the Earth's core. Geology and geo-
physics, 2, 152 – 172 (2009).
41. Berezniakov, A.I., Niemets, К.А.: Earth Physics. V.N. Karazin Kharkiv National Univer-
sity (2010).
42. Chernogor, L.F.: Physics and ecology of disasters. V.N. Karazin Kharkiv National Univer-
sity (2012).
43. Malyshkov, Yu.P., Dzhumabaev, K.B., Malyshkov, S.Yu.: Earthquake prediction method.
Patent of the Russian Federation, No. 2238575, МPК G01V3/00 (2004).
44. Tertyshnikov, A.V.: Harbingers of Severe Earthquake in the Ozonosphere, Heliogeophysi-
cal research, 2, 54 – 59 (2012).
45. Garcia, R., Crespon, F., Ducic, V., Lognonne, P.: Three-dimensional ionospheric tomogra-
phy of post-seismic perturbations produced by the Denali earthquake from GPS data. Geo-
phys. J. Int., 163, 1049 – 1064 (2005).
46. Heki, K., Ping, J.: Directivity and apparent velocity of the coseismic traveling ionospheric
disturbances observed with a dense GPS array. Earth Planet. Sci. Lett., 236, 845 – 855
(2005).