In July 2018, 30 short-focus weak earthquakes (with the hypocenter depths of 1 - 15 km) occurred in the region of Tiberias (Israel). According to the world classification  , an earthquake with a magnitude of up to 5 is considered weak, a magnitude up to 7—average, and a magnitude above 7—strong or destructive. The leadership of the Homeland Service has concluded that a global cataclysm could occur in the country in the near future, which would lead to the largest civilian casualties. Currently developed systems, capable of predicting earthquakes over a fairly long period of time before an earthquake occurs, have not yet achieved significant success (see Appendix). A new multi-year plan to protect the state from earthquakes has been announced. This article describes a proposal to determine the epicenter of a future short-focus earthquake in tens of hours and reduce the power of an impending catastrophic earthquake. The presentation of the project is based on the proposed physical and physical-mathematical models of short-focus earthquakes, which allow us to understand and describe the occurrence of such large periods of time for the harbingers of a catastrophic short-focus earthquake. It is intended that the results of the work are to be protected by patent. The prototype of the patent is a patent  and an article  based on measuring the flow of mercury vapor from the Earth to the atmosphere, which is a precursor of earthquakes.
Section 2 provides information on the use of mercury measurements as earthquake precursors. Section 3 describes the physical concepts of earthquakes. The physical model of a short-focus earthquake, which formulates the conditions for the occurrence of an earthquake, is presented in Section 4. Section 5 presents the physical-mathematical model of an earthquake. Section 6 describes the proposed method for determining epicenter of earthquakes. Conclusions are presented in Section 7.
2. Information on the Use of Mercury Flux Measurements as Earthquake Precursors
An increase in mercury gas Hg was recorded in the Tashkent region after the 1966 earthquake, and it was suggested that Hg may be a precursor of an earthquake  . In Russia, equipment and methods for measuring mercury gas to determine information about a future earthquake have been developed since 1968  . In the paper  and in the patent  , the invention of the method for predicting earthquakes was described, and it forms the prototype for our further work. In the article  , the author cites the conclusion, made by Chinese experts, that mercury anomalies are the best precursor of earthquakes. The triumph of the program in China was the prediction of the earthquake of February 4, 1975 and the rescue of more than 20,000 people. The article  and the patent  describe the invention of a method for predicting earthquakes in seismically active areas based on recording the density of mercury vapors flux from the Earth to the atmosphere, which is a geochemical precursor of earthquakes. The essence of the invention: a periodic registration of a precursor of an earthquake (the density of mercury flux in the soil gas, rising from the Earth) is carried out, and then the possibility of an earthquake can be predicted based on the anomalous change in the value of the density of flux recorded over time. Empirical data accumulated during the monitoring at the observation point suggested that the content of mercury vapor in a closed cavity inside loose layers may change with an increased amplitude (10 - 1000 times) in a space of just a few hours, which correlates with subsequent earthquakes. The accumulated data shows that an earthquake is preceded by the appearance of a predecessor, evident in a short-term increase in the mercury content in the gas contained in the chamber, 4 - 40 hours before the earthquake. In addition, there are positive correlations between the amplitude of the predecessor, the power of the earthquake and the length of the period between the appearance of the precursor and the moment of the earthquake. Depending on the magnitude, such precursors can be recorded at a distance of 100 km or more from the sources of the earthquake. Figure 1 (Figure 1 of  ) presents the measurement of the density of mercury flux before the September 29, 1981 earthquake in Dushanbe, with a magnitude of 3 points and an epicenter located at a distance of 20 km from the measurement point, which showed a 90 times increase in the flow of mercury vapor 19 hours before the earthquake.
3. Conceptual Ideas about Earthquakes
In this section physical representations of earthquakes are considered, following the conceptual system of hypotheses described by Gilat A. and Vol A.    . This hypothesis system can explain all the actual facts about earthquakes and volcanic eruptions, which the conceptual system of hypotheses of Reid H.F.  could not explain with the hypothesis of an earthquake as a result of the displacement of rocks with accumulated elastic stress. The Reid model has problems in explaining the numerous geophysical and geochemical anomalies that precede and accompany earthquakes. The following phenomena are among these anomalies: 1) powerful gas emissions before and during an earthquake; 2) abnormally high infrared radiation, in other words, heat radiation, that precedes
earthquakes; 3) radio and acoustic noise before the earthquake; 4) diffused glow and ball lightning; 5) completely measureable rise of the ground, before an earthquake; 6) stress cycles associated with earthquakes and changes in groundwater regimes  .
During the Earth’s accretion period primordial hydrogen and helium, comprising 98% - 99% of the space matter, were trapped and stored in the Earth’s core and mantle through endothermic reactions as both solid and liquid solutions and chemical compounds. After the planet stabilized, the energy gained by the capture of H and He has been quasi-continuously released by the exothermic reactions of the degassing of the Earth. The resulting heat and continuous explosions produce the manifestations of magmatic activity in general, particularly volcanic eruptions and earthquakes. Analyses of gases in lavas from Kamchatka volcanoes show that primary explosive gases uncontaminated by meteoric water and air (H2, Cl2, CO, OH, F2, Br2, H2S, and CH4) comprise 10% - 70% of total volcanic gases. Saturated with energy gases break through fragile rocks and form passages for magma into magmatic chambers, as well as supplying with energy the earthquakes and volcanic eruptions. Tectonic earthquakes are presented and discussed as a series of chemical explosions caused by chemical chain branching reactions  . This declaration is supported by comparative analysis of seismograms from earthquakes and nuclear explosions, because the complexity of earthquakes is greater than that of artificial explosions. The accumulation of energy in the hypocenter includes the latent energy of chemically active (explosive) substances which induce exothermic explosions. Those gases ascend along faults and fractures from the depth of about 700 km to the Earth’s surface and form the complex of an earthquake. The following possible reactants participate in earthquakes: hydrogen—oxygen; hydrogen—halogens; hydrogen—sulfur; alkanes (methane, etc.)—oxygen; alkanes—halogens; alkanes—nitrates, etc. In the next section an explanation of the required conditions for the onset of an earthquake in Tiberias region is given.
4. Necessary Conditions for the Onset of an Earthquake
What conditions are necessary for an earthquake and what is the cause of more than 30 earthquakes in the region of Tiberias within one month of 2018? The Syrian-African Rift passes through Tiberias and opens up the possibility of the passage of explosive gases to the surface of the Earth. Statistics show that July is the month when earthquakes are most likely. Explosive gases on the way to the surface of the Earth get trapped, surrounded by poorly permeable rocks, which causes accumulation of gases in the traps. Explosiveness of combustible gases depends on their concentration and the concentration, of air, as described in  . The concentration of gas vapors in the air is calculated as a percentage of a given volume of air (% volume). A mixture of gas and air is explosive only if the content of the combustible gas in it is between the lower and upper limits of explosiveness. The concentration of a combustible gas in the air, below which the combustible mixture cannot explode, but above which it can burn as a result of the explosion, is called the lower explosive limit (LEL). The concentration of gas or vapor, below which an explosion occurs and above which the mixture is not capable of igniting, is called the upper explosion limit (UEL). The range between the LEL and the UEL is known as the explosion range for this gas. Table 1 (from  ) presents the explosive limits of certain gases of Kamchatka volcanoes (described in Section 3).
Thus, the necessary condition for the onset of an earthquake is that an explosive gas has to be in the explosiveness range for that gas. The natural gas consists mainly of methane (CH4)—from 70% to 98%. The composition of the natural gas may include heavier hydrocarbons—methane homologues: C2H6, C3H8, and C4H10. Natural gas also contains other substances that are not hydrocarbons: H2, H2S, N2, Hg, He, and other gases. Therefore, we assume that the main possible explosives associated with earthquakes in the Tiberias region are alkanes (methane, etc.)—oxygen. The next triggers that can cause an explosion can be: a series of branched chemical chain reactions, tidal interaction, and atmospheric pressure gradients. For deep-focus earthquakes, other compounds can play the role of an oxidizer, but in this article only small-focus earthquakes will be considered. Thus, the proposed earthquake model provides a physical explanation for the existence of earthquake precursors as follows. Let us consider the experimental measurements presented in Figure 1. The delta-shaped over time flow of Hg gas and an explosive gas in the area adjoining the fault fills earthquake traps with an explosive gas throughout the region, turning them into sources of future earthquakes. The density of an explosive gas after it fills the traps is significantly higher than the upper explosion limit (UEL) (see Table 1); therefore, an earthquake cannot immediately occur. The experimentally measured density of gas flow, as presented in Figure 1, is the forerunner of the earthquake. In this case, the delta-shaped gas flow should be observable throughout the region. The source of the earthquake, having been filled with gas, begins to diffuse the explosive gas in all directions. This process will continue until the gas concentration falls into the explosion range in a period of time T, which will be the time of an earthquake precursor, after which an earthquake will occur. Thus, this process allows us to explain the physical reason for the appearance of the precursor such a long period of time ahead of an earthquake. The origin of the 2 earthquakes in Figure 3 of  can be explained by the fact that there were two traps in the area of the fault, which, having filled with gas, exploded on December 27 and December 28, 1993.
Table 1. Explosive limits of various gases from the Earth.
Similar reasons (filling of all 30 traps by gas, coming from the Syrian-African fault) can explain the origin of 30 earthquakes in July 2018 in the region of Tiberias. The existence of traps that have the ability to become earthquake foci is determined by the requirement: the diffusion coefficient D (center) D (surrounding rocks). For example, for mercury D (samples of dry rocks) 0.1 (cm2/s) D (samples of natural moisture) 10−5 (cm2/s)  . The hydrothermal shell in the Tiberias region is saturated with water resources: first, it is the Kinneret; secondly, it is the Thermal springs of Khamei Tiberias located near the city of Tiberias. Therefore, the probability of the formation of earthquake foci in the region of Tiberius in the future is very high.
Earthquakes as a process are quite complicated and have different mechanisms. It can be assumed, that the following chemical reactions in the focus of a short-focus earthquake are an important mechanism influencing the onset of an earthquake.
In the conditions of lack of oxygen all hydrocarbons emit carbon dioxide and water. For example:
With an even smaller amount of oxygen, fine carbon (soot) is released:
For an earthquake to occur, oxygen is also needed in its focus. Therefore, the more time passes from the moment of filling the trap with gas, the less oxygen is present in the focus due to reactions (1), (2), which ultimately leads to an increase in the precursor time and, therefore, to a decrease in the power of the earthquake. It can be assumed that another important mechanism affecting the onset of an earthquake is a decrease in atmospheric pressure on the Earth’s surface (i.e., a cyclone burst). One of the signs of a cyclone is the prevailing area of calm weather in the center of the cyclone. Inside this area, a warmer temperature is formed than in the rest of the air flow and lower humidity is observed  .
The analysis of the observational data for the largest earthquakes with M ≥ 7.5 for the period 2000-2010  indicates a correlation of preceding negative atmospheric pressure gradients with the largest earthquakes of the 21st century. The atmospheric pressure gradients for 42 earthquakes are shown in Table 2 and in Figure 2.
The following explanation can be given for the trigger of atmospheric pressure gradients at the beginning of an earthquake, which follows from the presented
Table 2. Atmospheric pressure gradient distribution for earthquakes with M ≥ 7.5.
Figure 2. An example of atmospheric gradients pressure before an earthquake (depending on the number of days relative to an earthquake) with M ≥ 7.5, registered on the 2 closest to the epicenter of observation points. 11/16/2003. M = 7.8 Rat Islands (Aleutian Islands), Alaska. Figure 2 is given in  . Figure 2 shows the graphs for atmospheric pressure 30 days before and 5 days after the earthquake that occurred during a cyclone.
physical model of the onset of an earthquake. The formation of a cyclone (leading to a decrease in pressure) initiates the release of an explosive gas into the atmosphere; which in turn accelerates the decrease in the concentration of gas into the atmosphere; which in turn accelerates the decrease in the concentration of gas in the focus. This process will continue until the gas concentration falls within the explosion range, after which an earthquake will occur. Thus, the origin of thirty earthquakes in July 2018 may be associated precisely with the preceding negative gradients of atmospheric pressure in the Tiberias region. Note, that the overwhelming number of earthquakes in Israel occurred precisely in the summer months, when quiet weather with warm temperatures is usually observed  : July 11, 1927, August 3, 1993, June 2004, June 12, 2008, July 13, 2008, July 15, 2008, April 17, 2009, April 12, 2010, March 9, 2010, March 20, 2010, July 15, 2010, August 23, 2010.
5. Physical and Mathematical Models of a Short-Focus Earthquake
In this section, we consider the simplest physical and mathematical models of a short-focus earthquake associated with the diffusion of an explosive gas in all directions from the earthquake source, which will allow us to obtain a formula for estimating the time interval of an earthquake precursor.
The piezoelectric conductivity equation in the theory of filtration (describing the motion of a gas in a porous medium) is written as follows 
where piezoconductivity coefficient ; this range for was fixed for various environments. The accumulation of gas mass in the trap can be interpreted as the appearance of an instantaneous point source of mass m0, which begins to diffuse the gas in space. Then, the gas P(r, t) can be described by the formula (if the origin is placed in the earthquake source)  .
where , m0 is the mass of gas in the focus, , kB is Boltzmann constant, Tk is temperature.
The gas density u(r, t) can be described by the formula (if the origin is placed in the earthquake source) in the following form:
where U1 is the gas density at which an earthquake occurs, for example, (for CH4 and air) U1 = 0.08 (kg/m3). The earthquake energy E is equal to:
The formula linking magnitude M and energy of earthquakes was presented in  in such form:
where E is the earthquake energy, J is one joule. Then the formula for the time precursor of the earthquake T, using (4)-(8), can be written as following:
We calculate T = 19 hours, using (9) and , , , .
where a and b are free parameters. Thus, the derived formula (9) coincides in functional structure with the empirical formula (10). Moreover, formula (9) represents expressions for the free parameters a and b. Deviations from formulas (5) and (9) are observed in the processing of experimental measurements: earthquakes represented in a circle in Figure 3(a) in  with M < 5.8 represent a cloud of points that cannot be described by a straight line. We can offer the following explanation for these deviations from formulas (5) and (9) for earthquakes represented in a circle:
· the indicated earthquakes had a relatively small, different concentration of the oxidizer in the foci compared with the concentration of the oxidizer in the foci for earthquakes satisfying formulas (5) and (9);
· this scatter of small initial conditions of oxidizer concentration led to a significant increase in the time of the precursor T when the gas diffusion occurred due to the absence of the necessary explosion conditions and to the scatter of points;
· in turn, due to reactions (1) and (2) this led to an even greater increase in the time of the precursor T and, accordingly, to a significant decrease in the concentration of an explosive gas, when the necessary explosion condition became satisfied after the diffusion of the gas;
· as a result, the magnitude of earthquakes M decreased significantly and the points describing earthquakes, due to the initial spread of oxidizer concentrations in the source, fell into the area of the circle.
6. A Proposal of a System of Measures to Prepare for Future Earthquakes in the Tiberias Region
In the paper  and in the patent  the invention of the method for predicting earthquakes was described, which is the prototype of our further work. The summary of the prototype analysis: the examples of determining the hypocenter and epicenter were based on the empirical formula (2) in  , which uniquely describes only earthquakes with M > 5.5, the physical justification of the formula (2) in  is not presented.
The specific form of detection of mercury in rocks discovered in the course of research  under normal conditions corresponds to a state close to gaseous, in which atoms “impregnate” minerals throughout their entire volume and in the interstitial space of the lattices.
It is established that the reaction of rocks to deformations, accompanied by the release of Nd vapors from them, serves as an effective tool for studying the dynamic processes occurring in rocks, which is used in sensors.
1) A patent is proposed for a new device for recording information about earthquake precursors using mercury measurements as an extremely strong and sensitive forerunner occurring tens of hours before an earthquake. The proposed equipment for recording information should be able to determine the epicenter of a future earthquake.
2) It is proposed to create a network of stations in the north of Israel. At each station it is proposed to place a device for determining the epicenter of a future earthquake, as well as a device for determining the dependence of air pressure on time.
The prototype  gives a description of an invented method for predicting earthquakes in seismically active regions based on space-time variations in the density of mercury vapour flux from the Earth to the atmosphere, which is a geochemical precursor of earthquakes. The essence of the invention is to carry out periodic registration of the precursor of an earthquake in the set area of control. In order to do this, the magnitude of the density of mercury in the gas flux ascending from the Earth is registered at a depth of not less than 1 m. Then, according to an anomalous change of the density value over time, the possibility of an earthquake can be forecast. In the described patent, the mercury flux can come from any direction into one sensor; the information about the direction angle of where the flow came from is lost. For example, in Figure 1, after the trap is filled, the gas diffuses from the traps, but the information about the direction from which the gas flow came to the measuring device is not recorded. It is proposed in the prototype, that the directions were the flow from the trap came from to the measuring device can be found using the formula (2) from , which does not describe the situation with high accuracy.
In our proposed patent, the flow of mercury from any direction of a certain solid angle Ω0, comes into only one sensor. The number of sensors located in the hemisphere is N = 2π/Ω0. Thus, once the delta-shaped gas flow from the fault fills all the traps in the region, the gas begins to diffuse in all directions from all the traps in the region. Now, looking at the gas flows in the entire set of sensors in the hemisphere, we can immediately determine the sensor with the maximum gas flow from all traps in the region and find the solid angle from which the maximum gas flow came to every station from all earthquakes. Considering together the maximum gas flows from several stations located in the region, we can determine the epicentres and hypocentres of future earthquakes.
7. Prevention of Catastrophic Earthquakes by Means of Significantly Decreasing the Magnitude of a Future Earthquake
The method of finding the epicenter of a future earthquake is described in the Section 6. This section describes the proposed reduction in the magnitude of a future earthquake by slowing down the process of an explosive gas release from the earthquake source, which should prevent catastrophic earthquakes basing on the paper  . Studies of excited seismicity in the reservoir area. The dams of the Nurek hydroelectric station showed that with an increase in the level of vibrations during spillways, the average daily seismic energy normalized to an area of 1000 km2 decreases by 3 orders of magnitude and strong earthquakes disappear. Figure 3 shows that with an increase in the level of vibrations from 0.02 μm (the natural level of microseisms) to values of 0.4 μm, the average daily released seismic energy, normalized to an area of 1000 km2, decreases by 3 orders of magnitude. In separate periods of time, when the amplitudes of technogenic vibrations in the adjacent territories reached a maximum value of 0.5 μm during spillway, earthquakes completely disappeared.
The authors suggest the following set of activities:
1) Creation of a network of seismic observations in a seismically active fault
Figure 3. Schedule of the average daily released earthquake energy per unit area of 1000 km2 in the water area of the Nurek hydroelectric dam, depending on the level of microseisms at different distances and at different time periods for gradually receding zones within a radius of 25 km from the spillway from  .
zone with an energy level of representative recording of K ≥ 7 (M = 1.5). 2) Determination of the most probable places of expected strong earthquakes within a large seismic generating structure. 3) Creation of a complex of measurements of horizontal and vertical movements. Earth surface using GPS/GLONAS systems, including deformation measurements in adits. 4) Creation of a complex for regulating water injection into wells and vibration effects in areas of the likely occurrence of strong earthquakes. 5) Periodic injection of fluid into the wells in the phases of gravity lunar-solar tides of the earth’s crust, causing the expansion of tension environments. 6) Conducting vibrations of the earth’s crust in areas of likely occurrence of earthquakes with an intensity of up to 0.5 microns.
Based on the derivation of the formula for the time precursor of the earthquake T in Section 5, we write the time precursor of the earthquake T1 when the vibro-seismic device creates an oscillation with amplitude a and frequency ώ
, (c0= 0.625 µm−1). (12)
Figure 4 is constructed using formulas (9)-(10). Vibration exposure slows the release of gases and oxygen. As a result of a decrease in the emission of gases and oxygen, chemical reactions (1) and (2) will play a greater role, which should lead to an even greater violation of the necessary earthquake conditions and, consequently, to an even greater decrease in the magnitude of the earthquake. Thus, we propose to fix the epicenter of the earthquake with the help of a new sensor. Then we turn on vibrator in the epicenter of an earthquake.
Figure 4. The ratio of the time precursor of the earthquake T1 when the vibro-seismicdevice creates a wobble to time T precursor of the earthquake.
The article presents a physical model of the conditions necessary for an earthquake to start, described for a method based on the registration of flows of mercury vapors in a gas rising from the Earth. This model gives an explanation of why a precursor appears such a long period of time ahead of an earthquake (from several hours to hundreds of hours), using a physical-mathematical model. The characteristic times of an earthquake precursor for seismic methods are tens of seconds. The formula for the time of the precursor of the future earthquake, derived from the proposed physical and mathematical model of the earthquake, made it possible to explain and describe the increase in the time when the precursor appears, depending on the magnitude of the earthquake. The developed physical model of a short-focus earthquake is the basis for the proposed project to determine the epicenter of a future short-focus earthquake tens of hours ahead of its start, and to reduce the power of an impending catastrophic earthquake. The proposed system of measures should save the lives of citizens from future catastrophic short-focus earthquakes. This article is a breakthrough in the following directions: a physical model of the onset of an earthquake has been proposed; it was formulated a necessary condition for the onset of an earthquake, which allowed explaining and theoretically describing the anomalously large time of the appearance of an earthquake precursor; Based on the proposed chemical formulas, an explanation of the deviations of the experimental observations from the formulas in the prototype was given; The article contains a patent proposal for a sensor and methods for reducing the magnitude of future catastrophic short-focus earthquakes. The results of the work are supposed to be protected by a patent. To implement the author’s proposals, it is necessary to carry out the following program: development and testing of the proposed sensor, which allows determining the epicenters of future earthquakes. It is proposed to create a network of stations in the north of Israel. At each station it is proposed to place a device for determining the epicenter of a future earthquake, as well as a device for determining the dependence of air pressure on time. The method of reducing the magnitude of an impending catastrophic earthquake is based on the proposed physical model of the onset of an earthquake and is implemented by the action of a vibration source in the region of the detected earthquake epicenter. Theoretical modeling of chemical processes should be the next stage of theoretical studies.
The author is grateful to Lev Gilat for sending of some materials and for helpful discussions. The author also thanks Valery Cherkasky for their assistance with making of figures.
The Canadian company “Nanometrix” offered an earthquake early warning system; this technology will allow being informed about an earthquake with an interval of 10 - 30 seconds between the first strong earthquake and the destructive wave following it  . In addition, the Geological Survey of Israel (GSI) is developing a system called “Troy” (“Trumpet Sound”). According to information received from the head of this seismological project, when the alarm signals are heard, tremors will already be felt in Tiberias if the epicenter is in the Kinneret area  . Thus, attempts to find and implement a system capable of predicting earthquakes over a sufficiently long period of time before an earthquake have not yet achieved significant success.
The delta-shaped appearance of increased pressure in the earthquake source leads to the appearance of a filtering pressure wave propagating in all directions, and, accordingly, to a measurement point on the Earth’s surface (pressure filtration waves have long been known and, in particular, are used to study layer parameters [Ovchinnikov, GG Kushtanova, AG Gavrilov, MV Filtrational Pressure Waves as a Method of Reservoir Flow Parameters Investigation, Oil and Gas business (electronic journal http://ogbus.ru/files/ogbus/issues/6_2015/ogbus_6_2015_p124-161_OvchinnikovMN_en_en.pdf)].
For calculations, we use the theory of propagation of intense acoustic waves in a medium without dispersion [Vinogradova MB, Rudenko OV, Sukhorukov AP Wave Theory: Nauka, M., 1979, 383 pp. https://nashol.com/2015033183698/teoriya-voln-vinogradova-m-b-rudenko-o-v-suhorukov-a-p-1979.html]. The problem of changing the amplitude and duration of a single triangular disturbance (see §3, chapter 6, c) described by the following expression for the initial form of the disturbance
where u is the vibrational velocity, T is the pulse duration. Formulas describing the change in the peak value of the vibrational velocity u2(x) and the pulse duration T(x) have the form
where x is the distance, , , is the adiabatic index of the medium, with the speed of the acoustic wave. Based on formulas (2.1), we write the formulas for the flux density of mercury vapor Fs(t) in the sensor  .