_{1}

We present the confirmed validity of the prediction (made in unpublished article (Simonenko, 2019b) presented on December 9, 2019 to the journal Energy Research) of the thermohydrogravidynamic theory (Simonenko, 2014) concerning the evaluated probabilities of the forthcoming intensification of the global seismotectonic, volcanic, climatic and magnetic processes of the Earth during the evaluated beginning (from December 7, 2019 to April 18, 2020) of the established in 2012 AD (Simonenko, 2012a) and finally confirmed in 2019 AD (Simonenko, 2019a) “first forthcoming subrange 2020 ÷ 2026 AD of the increased intensification of the Earth” determined by the combined non-stationary cosmic energy gravitational influences on the Earth of the planets (Mercury, Venus, Mars and Jupiter) and the Sun (owing to the gravitational interaction of the Sun with Jupiter, Saturn, Uranus and Neptune).

The problem of the long-term predictions of the strong earthquakes (Richter, 1958) is the significant problem of the modern geophysics (Sgrigna & Conti, 2012). It is well known (Sgrigna & Conti, 2012) that “the deterministic prediction of the time of origin, hypocentral (or epicentral) location, and magnitude of an impending earthquake is an open scientific problem”. It was pointed out (Sgrigna & Conti, 2012) with some regret that the modern “study of the physical conditions that give rise to an earthquake and the processes that precede a seismic rupture of an ordinary event is at a very preliminary stage and, consequently, the techniques of prediction of time of origin, epicentre, and magnitude of an impending earthquake now available are below standard”. In the special issue (Console, Yamaoka, & Zhuang, 2012) of the International Journal of Geophysics, Rodolfo Console, Koshun Yamaoka, and Jiancang Zhuang assessed the status of the art of earthquake forecasts and their applicability. It was conjectured (Console, Yamaoka, & Zhuang, 2012) that the recent destructive earthquakes occurred in Sichuan (China, 2008), Italy (2009), Haiti (2010), Chile (2010), New Zealand (2010), and Japan (2011) “have shown that, in present state, scientific researchers have achieved little or almost nothing in the implementation of short- and medium-term earthquake prediction, which would be useful for disaster mitigation measures”.

In this article, following the global prediction thermohydrogravidynamic principle (11) used for the real planetary configurations of the Earth and the planets of the Solar System, the author presents the confirmed validity of the prediction (of the unpublished article (Simonenko, 2019b) presented on December 9, 2019 to the journal Energy Research) of the thermohydrogravidynamic theory concerning the evaluated probabilities of the strongest earthquakes during the established beginning (from December 7, 2019 to April 18, 2020) of the predicted (Simonenko, 2012a, 2014) and confirmed (Simonenko, 2019a) “first forthcoming subrange 2020 ÷ 2026 AD of the increased intensification of the Earth” determined by the combined non-stationary cosmic energy gravitational influences on the Earth of the system Sun-Moon, Mercury, Venus, Mars, Jupiter and the Sun (owing to the gravitational interaction of the Sun with Jupiter, Saturn, Uranus and Neptune).

To do this, in Section 2.1 we present the established (Simonenko, 2006, 2007a) generalized differential formulation (7) of the first law of thermodynamics (for moving rotating deforming compressible heat-conducting stratified individual finite continuum region τ subjected to the non-stationary Newtonian gravitational field) extending the classical identical formulation (Gibbs, 1873; Landau & Lifshitz, 1976) by taking into account (along with the classical infinitesimal change of heat δ Q and the classical infinitesimal change of the internal energy dU τ ) the established (Simonenko, 2006, 2007a) differential (infinitesimal) energy gravitational influence dG (as the result of the Newtonian non-stationary cosmic and terrestrial gravitation) on the individual finite continuum region τ during the infinitesimal time interval dt . Based on the established (Simonenko, 2006, 2007a) differential energy gravitational influence dG , in Section 2.2 we present the established (Simonenko, 2012a, 2014) global prediction thermohydrogravidynamic principles (11) and (12) determining the maximal temporal intensifications of the global and regional natural (seismotectonic, volcanic, climatic and magnetic) processes of the Earth.

Using the global prediction thermohydrogravidynamic principle (11), we present in Section 3.1 the prediction of the thermohydrogravidynamic theory concerning the evaluated probabilities (14), (15), (16), (17), (18), (19), (20), (21), (22), (23), (24), (25), (26), (27), (28), (29), (30), (31), (32), (33), (34) and (35) of the strongest earthquakes during the beginning (from September 24, 2019 to April 18, 2020) of the established (Simonenko, 2012a, 2014, 2019a) first forthcoming subrange 2020 ÷ 2026 AD of the increased intensification of the Earth.

Taking into account the considered (on December 6, 2019 based on

Based on the corrected probabilities Pr cor,i (given by (37) for i from i = 9 to i = 22) of the possible strongest (near the calculated numerical time moment (13) from December 7, 2019 to April 18, 2020) intensifications of the seismotectonic, volcanic, climatic and magnetic processes of the Earth (Simonenko, 2019b), we present in Section 3.3 the calculation of the mean date t ¯ ( 2019 , 2020 ) of the probable most strongest earthquake of the Earth during the considered (Simonenko, 2019b) range December 7, 2019 ÷ April 18, 2020.

In Section 4 we present the main results and conclusions.

We shall assume that τ is an individual finite continuum region bounded by the closed continual boundary surface ∂ τ considered in the three-dimensional Euclidean space with respect to a Cartesian coordinate system K. We shall consider the individual finite continuum region τ in a Galilean frame of reference with respect to a Cartesian coordinate system K centred at the origin O and determined by the axes X 1 , X 2 , X 3 (see

The position-vector r c of the mass center C of the individual finite continuum region τ in the K-coordinate system is given by the following expression

r c = 1 m τ ∭ τ r ρ dV , (1)

where m τ is the mass of the individual finite continuum region τ , dV is the mathematical differential of the physical volume of the individual finite continuum region τ , ρ ≡ ρ( r ,t) is the local macroscopic density of mass distribution, r is the position-vector of the continuum volume dV . The speed of the mass centre C of the individual finite continuum region τ is defined by the

expression V c = d r с dt .

We shall use the differential formulation of the first law of thermodynamics (De Groot & Mazur, 1962) for the specific volume ϑ = 1 / ρ of the one-component deformed continuum with no chemical reactions:

du dt = dq dt − p d ϑ dt − ϑ П : Grad v , (2)

where u is the specific (per unit mass) internal thermal energy, p is the thermodynamic pressure, П is the viscous-stress tensor, v is the hydrodynamic

velocity of the continuum macrodifferential element (De Groot & Mazur, 1962), the operator d/dt = ∂ / ∂ t + v ⋅ ∇ denotes the total derivative (De Groot & Mazur, 1962; Batchelor, 1967; Gyarmati, 1970; Landau & Lifshitz, 1988) following the continuum substance, dq is the differential change of heat across the boundary of the continuum region (of unit mass) related with the thermal molecular conductivity described by the heat equation (De Groot & Mazur, 1962):

ρ dq dt = − div J q , (3)

where J q is the heat flux. The viscous-stress tensor П is taken from the decomposition Ρ = p δ + Π of the pressure tensor P (De Groot & Mazur, 1962), where δ is the Kronecker delta-tensor. The macroscopic local mass density ρ of mass distribution and the local hydrodynamic velocity v of the macroscopic velocity field is determined by the classical hydrodynamic continuity equation (De Groot & Mazur, 1962; Batchelor, 1967; Landau & Lifshitz, 1988):

∂ ρ ∂ t + div(ρ v ) = 0 (4)

under the absence of distributed space-time sources of mass output.

Using the differential formulation (2) for the specific (per unit mass) internal thermal energy u of the one-component deformed continuum with no chemical reactions, the heat equation (3), the decomposition Ρ = p δ + Π of the pressure tensor P , the hydrodynamic continuity equation (4), the classical equation (Batchelor, 1967; Gyarmati, 1970) for each variable f (such as v 2 / 2 , u and ψ ):

d dt ∭ τ f ρd V = ∭ τ d dt f ρd V, (5)

and the general equation of continuum movement (Gyarmati, 1970):

d v dt = 1 ρ div T + g (6)

for the deformed continuum characterized by the symmetric stress tensor Т = − P (Gyarmati, 1970) of a general form and taking into account the time variations of the potential ψ of the non-stationary gravitational field (characterized by the local gravity acceleration vector g = − ∇ ψ ) inside of the individual finite continuum region τ , we derived (Simonenko, 2006, 2007a) the generalized differential formulation (for the Galilean frame of reference) of the first law of thermodynamics (for moving rotating deforming heat-conducting stratified individual finite one-component continuum region τ subjected to the non-stationary Newtonian gravitational field and to non-potential stress forces characterized by the symmetric stress tensor Т ):

dU τ + dK τ + dπ τ = δQ + δA np, ∂ τ + dG (7)

where δA np, ∂ τ is the established generalized (Simonenko, 2006, 2007a) differential work done during the infinitesimal time interval dt by non-potential stress forces (pressure, compressible and viscous forces for Newtonian continuum) acting on the boundary surface ∂ τ of the individual finite continuum region τ ; δQ is the classical (Gibbs, 1873; Landau and Lifshitz, 1976) differential change of heat of the individual finite continuum region τ related with the thermal molecular conductivity of heat across the boundary ∂ τ of the individual finite continuum region τ ;

dG = dt ∭ τ ∂ ψ ∂ t ρd V (8)

is the established (Simonenko, 2007a, 2007b) differential (infinitesimal) energy gravitational influence (as the result of the Newtonian non-stationary gravitation) on the individual finite continuum region τ (during the infinitesimal time interval dt ); π τ is the established (Simonenko, 2006, 2007a) macroscopic potential energy (of the individual finite continuum region τ ) related with the non-stationary potential ψ of the gravitational field; U τ is the classical (Gibbs, 1873; Landau & Lifshitz, 1976) microscopic internal thermal energy of the individual finite continuum region τ ; K τ is the established (Simonenko, 2006, 2007a) instantaneous macroscopic kinetic energy of the individual finite continuum region τ .

Taking into account the general relation (8) for the infinitesimal energy gravitational influence dG on the individual finite continuum region τ , we obtained (Simonenko, 2014, 2019) the following relation for the combined differential cosmic non-stationary energy gravitational influence dG(τ c,r ) (during the infinitesimal time dt ) of the Solar System (the planets, the Moon and the Sun owing to the gravitational interaction of the Sun with Jupiter, Saturn, Uranus and Neptune) on the internal rigid core τ c,r of the Earth:

dG(τ c,r ) = dt ∭ τ c,r ∂ ψ comb ∂ t ρ c,r dV, (9)

where ρ c,r is the mass density of the internal rigid core τ c,r , ∂ ψ comb / ∂ t ≡ ∂ ψ comb (τ c,r ,t)/ ∂ t is the partial derivative (of the combined cosmic gravitational potential ψ comb ≡ ψ comb (τ c,r ,t) in the internal rigid core τ c,r of the Earth) approximated as follows (Simonenko, 2014, 2019)

∂ ψ comb (τ c,r ,t) ∂ t = ∑ i = 1,i ≠ 3 9 ∂ ψ 3i (C 3 ,t) ∂ t + ∑ j = 5 8 ∂ ψ 3j S (C 3 ,t) ∂ t + ∂ ψ 3MOON (C 3 ,t) ∂ t , (10)

where ∂ ψ 3i (C 3 ,t)/ ∂ t is the partial derivative (Simonenko, 2009, 2012a, 2013) of the gravitational potential ψ 3i (C 3 ,t) created by the planet τ i (of the Solar System) at the mass center C 3 of the Earth; ∂ ψ 3j S (C 3 ,t)/ ∂ t is the established (Simonenko, 2012a, 2019a) partial derivative of the gravitational potential ψ 3j S (C 3 ,t) created by the Sun (due to the gravitational interaction of the Sun with the outer large planet τ j , j = 5, 6, 7, 8) at the mass center C 3 of the Earth; ∂ ψ 3MOON (C 3 ,t)/ ∂ t is the partial derivative (Simonenko, 2009, 2012a, 2013) of the gravitational potential ψ 3MOON (C 3 ,t) created by the Moon at the mass center C 3 of the Earth.

Based on the generalization (7) of the first law of thermodynamics (used for the internal rigid core τ c,r of the Earth), we formulated (Simonenko, 2012a, 2014) the global prediction thermohydrogravidynamic principles determining the maximal temporal intensifications of the established thermohygrogravidynamic processes (in the internal rigid core τ c,r and in the boundary region τ rf between the internal rigid core τ c,r and the fluid core τ c,f of the Earth considered as a whole) subjected to the combined cosmic energy gravitational influence of the planets of the Solar System, the Sun (owing to the gravitational interaction of the Sun with the outer large planets) and the Moon. We concluded (Simonenko, 2012a, 2014) (based on the generalization (7) of the first law of thermodynamics used for the internal rigid core τ c,r of the Earth) that the maximal intensifications of the established thermohygrogravidynamic processes are related with the corresponding maximal intensifications of the global and regional natural (seismotectonic, volcanic, climatic and magnetic) processes of the Earth.

The rigorous global prediction thermohydrogravidynamic principles (determining the maximal temporal intensifications near the time moments t = t ∗ (τ c,r ) and t = t (τ c,r ) ∗ , respectively, of the thermohydrogravidynamic processes in the internal rigid core τ c,r and in the boundary region τ rf between the internal rigid core τ c,r and the fluid core τ c,f of the Earth) are formulated as follows (Simonenko, 2012a, 2014):

ΔG(τ c, r ,t ∗ (τ c,r )) = max t ∫ t 0 t d t ′ ∭ τ c,r ∂ ψ comb ∂ t ′ ρ c,r dV − local maximum for time moment t ∗ (τ c,r ), (11)

and

ΔG(τ c, r ,t (τ c,r ) ∗ ) = min t ∫ t 0 t d t ′ ∭ τ c,r ∂ ψ comb ∂ t ′ ρ c,r dV − local minimum for time moment t (τ c,r ) ∗ . (12)

The global prediction thermohydrogravidynamic principles (11) and (12) define the maximal and minimal combined cosmic integral energy gravitational influences ((11) and (12), respectively, for the time moments t = t ∗ (τ c,r ) and t = t (τ c,r ) ∗ ) on the considered internal rigid core τ c,r (of the Earth) subjected to the combined cosmic integral energy gravitational influence of the planets of the Solar System, the Sun (owing to the gravitational interaction of the Sun with the outer large planets) and the Moon.

To determine the forthcoming date t ∗ (τ c,r , 2020) corresponding to the first maximal combined planetary and solar integral energy gravitational influence on the internal rigid core τ c,r of the Earth (and on the Earth as a whole) during the established first forthcoming subrange (2023 ± 3)AD (Simonenko, 2012a, 2014, 2019a) of the increased intensification of the global seismotectonic, volcanic, climatic and magnetic activity of the Earth, we use the global prediction thermohydrogravidynamic principle (11) determining the maximal temporal intensifications of the global natural (seismotectonic, volcanic, climatic and magnetic) processes of the Earth. Based on the global prediction thermohydrogravidynamic principle (11) and considering the real planetary configurations of the Earth and the planets of the Solar System, we obtained on December 4, 2019 (Simonenko, 2019b) the forthcoming date (corresponding approximately to January 6, 2020):

t ∗ (τ c,r ,2020) = 2020 .016666667AD (13)

related with the first maximal (during the range (2020 ÷ 2 026 )AD ) combined planetary and solar integral energy gravitational influence (11) on the internal rigid core τ c,r (and on the Earth as a whole).

Considering the range (1960 ÷ 2018) AD and analyzing the previous strongest earthquakes (occurred near the calculated dates t ∗ (τ c,r , (1960 + m)) corresponding for m= 0, 1, …, 58 to the maximal combined planetary and solar integral energy gravitational influences (11) on the internal rigid core τ c,r of the Earth), we calculated on December 4, 2019 (Simonenko, 2019b) the following probabilities of the possible (near the calculated numerical time moment (13) from September 24, 2019 to April 18, 2020) strongest earthquakes and related (Simonenko, 2007b, 2009, 2012a, 2013, 2014, 2016, 2019a) strongest volcanic, climatic and magnetic processes) during the corresponding ranges (indicated in the round brackets):

Pr 1 {t e,max,2020 ∈ ( September24 ÷ September27,2019)} = 0 , (14)

Pr 2 {t e,max,2020 ∈ ( September 28 ÷ October 7, 2019)} = 5 55 = 0 .090909, (15)

Pr 3 {t e,max,2020 ∈ ( October8 ÷ October17,2019)} = 2 55 = 0 .036363, (16)

Pr 4 {t e,max,2020 ∈ ( October18 ÷ October27,2019)} = 1 55 = 0 .018181, (17)

Pr 5 {t e,max,2020 ∈ ( October 28 ÷ November 6, 2019)} = 2 55 = 0 .036363, (18)

Pr 6 {t e,max,2020 ∈ ( November7 ÷ November16,2019)} = 4 55 = 0 .072727, (19)

Pr 7 {t e,max,2020 ∈ ( November17 ÷ November26,2019)} = 2 55 = 0 .036363, (20)

Pr 8 {t e,max,2020 ∈ ( November 27 ÷ December 6, 2019)} = 4 55 = 0 .072727, (21)

Pr 9 {t e,max,2020 ∈ ( December 7 ÷ December 16, 2019)} = 3 55 = 0 .054545, (22)

Pr 10 {t e,max,2020 ∈ ( December17 ÷ December26,2019)} = 3 55 = 0 .054545, (23)

Pr 11 {t e,max,2020 ∈ ( December27,2019 ÷ January5,2020)} = 3 55 = 0 .054545, (24)

Pr 12 {t e,max,2020 ∈ ( January6 ÷ January15,2020)} = 5 55 = 0 .090909, (25)

Pr 13 {t e,max,2020 ∈ ( January16 ÷ January25,2020)} = 2 55 = 0 .036363, (26)

Pr 14 {t e,max,2020 ∈ ( January26 ÷ February4,2020)} = 1 55 = 0 .018181, (27)

Pr 15 {t e,max,2020 ∈ ( February5 ÷ February14,2020)} = 2 55 = 0 .036363, (28)

Pr 16 {t e,max,2020 ∈ ( February 15 ÷ February 24, 2020)} = 3 55 = 0 .054545, (29)

Pr 17 {t e,max,2020 ∈ ( February25 ÷ March5,2020)} = 4 55 = 0 .072727, (30)

Pr 18 {t e,max,2020 ∈ ( March6 ÷ March15,2020)} = 3 55 = 0 .054545, (31)

Pr 19 {t e,max,2020 ∈ ( March16 ÷ March25,2020)} = 3 55 = 0 .054545, (32)

Pr 20 {t e,max,2020 ∈ ( March26 ÷ April4,2020)} = 1 55 = 0 .018181, (33)

Pr 21 {t e,max,2020 ∈ ( April5 ÷ April14,2020)} = 1 55 = 0 .018181, (34)

Pr 22 {t e,max,2020 ∈ ( April15 ÷ April18,2020)} = 1 55 = 0 .018181 . (35)

Magnitude M max,i of strongest earthquake in range i, Region of the Earth | Date t e,i of strongest earthquake in range i, in yr | The time range i | The probability Pr i | Δ i = | t e,i − t ∗ (τ c,r , 2 020) | in days |
---|---|---|---|---|

M max,8 = 6.0 , 41 km NW of Platanos, Greece | November 27, 2019 = 2019.907255 | November 27, 2019 ÷ December 6, 2019 | 4 55 | 39.96 days before the date t ∗ (2020) |

M max,7 = 6.4 , 16 km WSW of Mamurras, Albania | November 26, 2019 = 2019.904517 | November 17, 2019 ÷ November 26, 2019 | 2 55 | 40.96 days before the date t ∗ (2020) |

M max,6 = 7.1 , 138 km E of Bitung, Indonesia | November 14, 2019 = 2019.871663 | November 7, 2019 ÷ November 16, 2019 | 4 55 | 52.96 days before the date t ∗ (2020) |

M max,5 = 6.6, 136 km W of Neiafu, Tonga | November 4, 2019 = 2019.844284 | October 28, 2019 ÷ November 6, 2019 | 2 55 | 62.96 days before the date t ∗ (2020) |

M max,5 = 6.6 , 14 km E of Bual, Philippines | October 29, 2019 = 2019.829397 | October 28, 2019 ÷ November 6, 2019 | 2 55 | 68.4 days before the date t ∗ (2020) |

M max,4 = 6.4 , 63 km NNE of Isangel, Vanuatu | October 21, 2019 = 2019.807494 | October 18, 2019 ÷ October 27, 2019 | 1 55 | 76.4 days before the date t ∗ (2020) |

M max,3 = 6.4 , 7 km ENE of Columbio, Philippines | October 16, 2019 = 2019.793805 | October 8, 2019 ÷ October 17, 2019 | 2 55 | 81.4 days before the date t ∗ (2020) |

M max,2 = 6.7 , 69 km WSW of Constitucion, Chile | September 29, 2019 = 2019.746064 | September 28, 2019 ÷ October 7, 2019 | 5 55 | 98.83 days before the date t ∗ (2020) |

M max,1 = 6.5 , 11 km S of Kairatu, Indonesia | September 25, 2019 = 2019.735113 | September 24, 2019 ÷ September 27, 2019 | 0 | 102.83 days before the date t ∗ (2020) |

8) and the corresponding dates t e,i of strongest earthquake in the time ranges (from i =1 to i = 8) occurred in different regions of the Earth.

We saw on December 6, 2019 (Simonenko, 2019b) that the considered (in

Pr cor {t e,max,2020 ∈ ( December7,2019 ÷ April18,2020)} = ∑ i = 9 22 Pr cor,i = 1 (36)

of realization of the strongest (from September 24, 2019 to April 18, 2020) earthquake in the range from December 7, 2019 (the beginning of the time range i = 9) to April 18, 2020 (the end of the range i = 22).

Assuming the proportionalities Pr cor,i ~ Pr i for the time ranges i from i = 9 to i = 22, we calculated on December 6, 2019 (Simonenko, 2019b) the corrected probabilities (instead of the probabilities (22)-(35)) for i from i = 9 to i = 22

Pr cor,i = 55 55 − 20 Pr i = 55 35 Pr i (37)

of realization of the strongest (from September 24, 2019 to April 18, 2020) earthquake during the ranges i from i = 9 to i = 22.

Taking into account the corrected (on December 6, 2019) probabilities (37) (for i from i = 9 to i = 22) of the possible (near the calculated numerical time moment (13) from December 7, 2019 to April 18, 2020) strongest earthquakes (and related strongest intensification of volcanic, climatic and magnetic processes of the Earth from September 24, 2019 to April 18, 2020) during the corresponding ranges (indicated in the round brackets), we obtained on December 7, 2019 (Simonenko, 2019b) the maximal probability (for i = 12):

Pr cor,12 {t e,max,2020 ∈ ( January 6 ÷ January 15, 2020)} = 55 35 ⋅ 5 55 = 5 35 = 1 7 (38)

of realization of the strongest earthquake and related (Simonenko, 2007b, 2009, 2012a, 2013, 2014, 2019a) strongest volcanic, climatic and magnetic processes of the Earth from September 24, 2019 to April 18, 2020) during the corresponding range ( January6 ÷ January15,2020) numbered by i = 12. We saw (on December 6, 2019) that the most probable (strongest from September 24, 2019 to April 18, 2020) intensification of the seismotectonic, volcanic, climatic and magnetic processes of the Earth may occur during the evaluated range ( January6 ÷ January15,2020) adjoining to the time moment t ∗ (τ c,r ,2020) = 2020 .016666667AD (given by (13) and corresponding approximately to January 6, 2020) related with the first maximal (during the subrange (2020 ÷ 2 026 )AD ) combined planetary and solar integral energy gravitational influence (11) on the internal rigid core τ c,r of the Earth. This result is in agreement with the physical essence (Simonenko, 2012a, 2014, 2019a) of the global prediction thermohydrogravidynamic principle (11) determining the maximal temporal intensifications near the time moment t = t ∗ (τ c,r ) of the thermohydrogravidynamic processes in the internal rigid core τ c,r and in the boundary region τ rf between the internal rigid core τ rf and the fluid core τ c,f of the Earth.

Using the probabilities Pr cor,i (given by (37) for i from i = 9 to i = 22), we calculate (on October 6, 2020, especially for the presentation on the 10^{th} International Conference on Geology and Geophysics) the mean date t ¯ ( 2019 , 2020 ) (in yr) of the probable most strongest earthquake during the considered (Simonenko, 2019b) range (from December 7, 2019 to April 18, 2020):

t ¯ ( 2019 , 2020 ) = Pr cor,9 t e,max,9 + Pr cor,10 t e,max,10 + ... + Pr cor,22 t e,max,22 = ∑ i = 9 22 Pr cor,i t e,max,i , (39)

where t e,max,i is the mean date of the corresponding ranges in the probabilities (22)-(35). Using the formula (39), we calculate the mean date t ¯ ( 2019 , 2020 ) (in yr) of the probable most strongest earthquake during the considered (Simonenko, 2019b) range (from December 7, 2019 to April 18, 2020):

t ¯ ( 2019 , 2020 ) = ∑ i = 9 22 Pr cor,i t e,max,i = 2020 .0995139 . (40)

which corresponds to the theoretical date

t ¯ ( 2019 , 2020 ) = February5 , 2 0 2 0 AD (41)

According to the U.S. Geological Survey, the most strongest (during the considered range from December 7, 2019 to April 18, 2020) earthquake (characterized by the maximal 7.7-magnitude) of the Earth occurred 123 km NNW of Lucea, Jamaica on

January 28, 2020 AD. (42)

We obtain the difference of 8 days between the calculated theoretical date t ¯ ( 2019 , 2020 ) (given by (41)) and the real date (42).

We have presented in Section 2.1 the established (Simonenko, 2006, 2007a, 2009, 2012a, 2013) generalized differential formulation (7) (given for the Galilean frame of reference) of the first law of thermodynamics deduced rigorously based on the postulates of the non-equilibrium thermodynamics (De Groot & Mazur, 1962; Gyarmati, 1970) and hydrodynamics (Batchelor, 1967; Landau & Lifshitz, 1988). The generalized differential formulation (7) is valid for moving rotating deforming heat-conducting stratified individual finite one-component continuum region τ (characterized by the symmetric stress tensor T ) subjected to the non-stationary gravitational field. The generalized differential formulation (7) of the first law of thermodynamics extends the classical (Gibbs, 1873; Landau & Lifshitz, 1976) identical formulation by taking into account (along with the classical infinitesimal change of heat δQ and the classical infinitesimal change of the internal energy dU τ ) the established (Simonenko, 2004, 2006, 2007a) infinitesimal increment dK τ of the macroscopic kinetic energy K τ , the established (Simonenko, 2006, 2007a) infinitesimal increment dπ τ of the gravitational potential energy π τ , the established generalized (Simonenko, 2006, 2007a) infinitesimal work δA np, ∂ τ done on the individual finite continuum region τ by the surroundings of τ , and the established (Simonenko, 2007a) differential energy gravitational influence dG on the individual finite continuum region τ during the infinitesimal time interval dt due to the non-stationary gravitational field.

We have presented in Section 2.2 the established (Simonenko, 2012a, 2014, 2019a) global prediction thermohydrogravidynamic principles (11) and (12) determining the maximal temporal intensifications of the global and regional natural (seismotectonic, volcanic, climatic and magnetic) processes of the Earth.

Based on the global prediction thermohydrogravidynamic principle (11) and considering the real planetary configurations of the Earth and the planets of the Solar System, we have presented in Section 3.1 the obtained on December 7, 2019 (Simonenko, 2019b) the forthcoming date t ∗ (τ c,r ,2020) = 2020 .016666667AD (corresponding approximately to January 6, 2020) related with the first maximal (during the subrange (2020 ÷ 2 026 )AD ) combined planetary and solar integral energy gravitational influence (11) on the internal rigid core τ c,r and on the Earth as a whole. Based on the global prediction thermohydrogravidynamic principle (11) and analyzing the previous strongest earthquakes (occurred near the calculated dates t ∗ (τ c,r ,(1960 + m)) corresponding for m= 0, 1, …, 58 to the maximal combined planetary and solar integral energy gravitational influences (11) on the internal rigid core τ c,r of the Earth) during the considered range (1960 ÷ 2018) , we calculated on December 7, 2019 (Simonenko, 2019b) the probabilities (14), (15), (15), (16), (17), (18), (19), (20), (21), (22), (23), (24), (25), (26), (27), (28), (29), (30), (31), (32), (33), (34) and (35) of the possible strongest (near the calculated numerical time moment (13) from September 24, 2019 to April 18, 2020) intensifications of the seismotectonic, volcanic, climatic and magnetic processes (during the corresponding ranges indicated in the round brackets of the probabilities (14), (15), …, and (35)) determined by the combined integral energy gravitational influences on the Earth of the planets (Mercury, Venus, Mars and Jupiter) and the Sun due to the gravitational interactions of the Sun with Jupiter, Saturn, Uranus and Neptune.

Taking into account that the considered (on December 6, 2019 (Simonenko, 2019b) based on

We have presented in Section 3.3 the calculated (on October 6, 2020, especially for the presentation on the 10^{th} International Conference on Geology and Geophysics) the mean date t ¯ ( 2019 , 2020 ) (given by (41) of the probable most strongest earthquake of the Earth during the considered range December 7, 2019 ÷ April 18, 2020. We have obtained the reasonably small difference of 8 days between the calculated theoretical date t ¯ ( 2019 , 2020 ) (given by (41)) and the real date January 28, 2020 AD (given by (42)) of the most strongest earthquake (characterized by the maximal 7.7-magnitude) of the Earth occurred 123 km NNW of Lucea, Jamaica during the established (Simonenko, 2019b) range from December 7, 2019 to April 18, 2020. We can conclude that the reasonably small difference of 8 days is based on the used simple approximation of the circular orbits (Simonenko, 2007b, 2012a, 2013) of the planets around the Sun and by disregarding the established partial derivative ∂ ψ 3MOON (C 3 ,t)/ ∂ t (Simonenko, 2009, 2012a, 2013) of the gravitational potential ψ 3MOON (C 3 ,t) created by the Moon at the mass center C 3 of the Earth.

Taking into account the prediction (Simonenko, 2012a, 2014, 2019a) of the thermohydrogravidynamic theory (Simonenko, 2006, 2007b, 2009, 2012a, 2013, 2014, 2016, 2019a) concerning the first forthcoming subrange 2020 ÷ 2026 AD of the increased intensification of the Earth and taking into account especially that the calculated (Simonenko, 2019a) date

2021.1 AD (43)

is related with maximal (during the all first subrange 2020 ÷ 2026 AD) energy gravitational influence on the Earth of the planets (Mercury, Venus, Mars and Jupiter) and the Sun (owing to the gravitational interaction of the Sun with Jupiter, Saturn, Uranus and Neptune), we have concluded (on October 6, 2020 as the main conclusion for the presentation on the 10^{th} International Conference on Geology and Geophysics) that it is very important to make in the near future the urgent identical calculations (of the probabilities of the possible strongest intensifications of the seismotectonic, volcanic, climatic and magnetic processes for the time range containing the calculated date 2021.1 AD) for Japan (Simonenko, 2018), Italy (Simonenko, 2018), Greece (Simonenko, 2018), China (Simonenko, 2018) and Chile (Simonenko, 2018), for which we have established (Simonenko, 2018) the confirmed cosmic energy gravitational genesis of the strongest earthquakes.

The validity of this main conclusion (of the author’s presentation on the 10^{th} International Conference on Geology and Geophysics) is supported also by the fact that the previous prediction (Simonenko, 2017) of the thermohydrogravidynamic theory (Simonenko, 2006, 2007b, 2009, 2012a, 2013, 2014, 2016, 2019a) concerning the strongest intensifications of the seismotectonic and climatic processes in California (since 9 August, 2017 and before 3 March, 2018) was confirmed really (Simonenko, 2019c).

Avoid the stilted expression, S. V. S. thanks the Editor, Prof. Rolin Zhang and Organizing Committee for publication of the article in this conference journal. The author thanks also the Editor with gratitude for the editorial comments and corrections improving the final text of the article.

The author thanks also Dr. Ma Deyi, Academician of CAS and Dr., Prof. Victor A. Akulichev, Academician of RAS for support (on the 2^{nd} Russia-China Symposium on Marine Science occurred on 10-13 October, 2012 in Vladivostok) of the Fundamentals of the thermohydrogravidynamic theory of the global seismotectonic, volcanic and climatic variability of the Earth (Simonenko, 2012b) resulted to this article.

The author declares no conflicts of interest regarding the publication of this paper.

Snko, S. V. (2020). The Confirmed Validity of the Thermohydrogravidynamic Theory Concerning the Forthcoming Intensification of the Global Natural Processes from December 7, 2019 to April 18, 2020 AD. Journal of Geoscience and Environment Protection, 8, 351-367. https://doi.org/10.4236/gep.2020.810022