The mass extinction events occurred regularly through the Phanerozoic, see for detail review in [1]. In particular, in the Jurassic period, to which this study is devoted, the Triassic-Jurassic (201.64 Ma) and Toarcian (182.60 Ma) mass extinctions happened.

The mass extinctions were often related to kill mechanisms such as marine anoxia, global warming, ocean acidification coupled with changes in atmospheric greenhouse gases, toxic metal poisoning, meteorite impact and cosmic gamma rays. Now it is increasingly widely thought that large igneous province (LIP) eruptions might be the driver of many of the purported proximal kill mechanisms [1]. However, let’s recall that mass extinction theories have developed from the simple “death-by-sea-level-change” hypothesis, which was proposed almost fifty years ago by Newell [2]. Hallam and Wignall confirmed Newell’s regression hypothesis for at least some major and minor extinction events [3] - [7]. In particular, in these studies it was proposed sea-level rose during the period from the latest Permian to the earliest Triassic, and that the oceanic anoxia caused by the continuing sea level raised that triggered the mass extinction. In [8], it was pointed out that the spread of anoxic bottom waters associated with marine transgression, sometimes, but not always, preceded by a major regression. It was also a potent extinction mechanism, presumably because of the severe reduction of viable habitat area. Authors took attention on the fact that the ultimate cause of the sea-level changes is generally unclear due to a glacioeustatic driving mechanism that can be convincingly demonstrated only for the end Ordovician and end Devonian events.

However, in [9] it was asserted that it is unlikely that sea-level fall played a significant role in the Triassic-Jurassic boundary extinctions in either a local or a global context. Detail discussion of a sea level role in mass extinction and discussions could be found in [9] - [13]. Thus, it was formed opinion that the Triassic-Jurassic mass extinction was related to a pronounced eustatic sea-level rise and partly to tectonic collapse anticipating the formation of ocean nearby, a phenomenon bound up with the creation of the proto Mediterranean and Atlantic oceans.

In particular in [14] [15], it was pointed that the cause of the end-Triassic mass extinction was probably linked to the contemporary activity of the Central Atlantic Magmatic Province (CAMP), which heralded the breakup of the super continent Pangaea. In that way, the possible kill mechanisms associated with magmatic activity include sea-level changes, marina anoxia, climatic changes, release of toxic compounds, and acidification of seawater.

As Triassic-Jurassic boundary mass extinction, the Toarcian mass extinction that happened in the Early Jurassic was an object for many studies and it was excellently described in the literature [1]. The Toarcian extinctions happened in the several parts of the world, such as Northwest Europe, South America and North America, Tibet and Japan, which proved that these mass extinctions have revealed the global nature of the crisis. According to [1] [15] [16], the Toarcian mass extinction has relation to the Karoo and Ferrar Traps in southern Gondwana.

In a number of studies the relations between mercury and mass extinction were investigated, see e.g. [16] - [27]. The Hg enrichment in sediments could be derived from massive volcanism and LIPs, from the combustion of coal deposits, from a meteoritic source, or from biomass burning due to wildfires and soil erosion [28] or from post-depositional processes [23]. In [25], it was shown that the Hg and paleontological evidences from the same archive indicate that significant biotic recovery did not begin until CAMP eruptions ceased.

However, these studies did not provide a physical explanation for the correlation between sedimentary Hg enrichments and massive volcanism. Indeed, what is the fundamental geophysical difference between the massive CAMP and Karoo-Ferrar LIPs from St. Helen or Pinatubo eruptions? We will notice that on the slopes of St. Helen and Pinatubo the mercury rivers did not flow. In this study, we answer on a question of what distinguishes between volcanic emissions during shallow and deep convection in the inner layers of the Earth.

Further, we should also mention the studies, in which the possibility of the galaxy influence on extinction processes was discussed.

Napier and Clube proposed the idea that gravitational disturbances caused by the Solar System crossing the plane of the Milky Way galaxy are enough to disturb comets in the Oort cloud surrounding the Solar System [29] [30] [31] [32]. The disturbance sends comets towards the inner Solar System. It raises the chance of an impact. According to the hypothesis, this results in the Earth experiencing large impact events about every 30 million years. Further, this hypothesis was evolved to the “Shiva” hypothesis (Shiva-Hindu God of Destruction) and has been investigated in the series of studies [33] - [39].

Also, note that periodicity of extinctions in the geologic past was investigated by Raup and Sepkoski Jr. [40]. In this study, a definition of the conception of “bottlenecking” effect of mass extinction was introduced. In study [41], it was written about periodic mass extinctions and the Sun’s oscillation around the galactic plane (Figure 1(a)).

The NASA image of the Milky Way Galaxy is presented in Figure 1(b). In this study a photo of the Milky Way Galaxy, NASA/JPL-Caltech/R. Hurt (SSC/Caltech), [42], which is available at the NASA/JPL-Caltech website, was used. The general structure of the Milky Way galactic arms and the parameters of the galaxy can be found in [43] [44]. The places, where the trajectory of the Sun intersects the galaxy arms are indicated by red arrows.

The Sagittarius, Scutum-Crux, Norma and the Perseus arm’s location and their relation with Neogene-Paleogene, Cretaceous-Jurassic, Persian-Carboniferous and Silurian-Ordovician stages are discussed in [45] - [53]. Since the primary source of changes in the sea level is the galaxy arms, we briefly highlight

the question of the periodicity of galactic processes. During discussion in [54] of a correlation between long-term cyclicities in Phanerozoic sea-level sedimentary record and their potential drivers, it was highlighted that the potential drivers, in addition to major plate tectonic motions, are galaxy cosmic rays and the motions of the Solar System in the Milky Way, see also e.g. [48] [55] [56] [57]. Note that the galaxy influence was also considered in [45] [46] [49] [53] and [58] and references therein. However, in all these works, there is no mention about the possible mechanism of galaxy influence. In our opinion, this effect is through the activation of the earth natural reactor.

The analysis of cyclical nature of geological processes and their connection with galactic frequencies can be found in [59] [60]. The hypothesis about comets as the cause of mass extinction was based on the study by Alvarez et al. [61], in which was reported about iridium increases and extraterrestrial cause for the Cretaceous-Tertiary extinction. The discussion about the Alvarez impact theory could be found in [62] - [68] and in numerous references in them. The cyclical frequency of comet, which shower from the Oort Cloud, was studied also, e.g. see in [69] [70] and [71]. By contrast of geological records, sometimes the spectral analysis studies of comet report little or no evidence of statistically significant cycles in impact of crater ages, see for example discussion in [38] [72].

As it was shown above, there are many mass extinction hypotheses, in which the different reasons of mass extinction are substantiated. The arguments in these hypotheses are proved, but these hypotheses form the knowledge mosaic. Is it possible to merge together all these hypotheses? In this study, we try to do it.

The relationship between the mass extinctions, Milky Way Galaxy, natural reactors, LIPs, MORBs, continental drifts and the variation of ⁸⁷Sr/⁸⁶Sr, outlined above, is summarized and presented in Figure 2.

In this study we will link together the galactic processes, the activation of nuclear processes inside the Earth and on/in the Sun, as well as the drastic changes in the habitat conditions of land species. Actually, we will renovate the oldest Newell’s hypothesis, but we will give explanations based on absolute other principle. The problem of mass extinction of species will be considered using the example of the species extinction in Hettangian and Sinemurian (201.6 - 190 Ma, Early Jurassic), while a statistical analysis of the Saurischia temporal distribution (unranked species) will be given for the whole Jurassic.

Geosciences study the processes occurring in the atmosphere, in the depths of the seas and oceans, and inside our planet. Three basic branches of physical

sciences, such as atmospheric physics, oceanology and geophysics accordingly are presented. The geophysicists explore the Earth’s core and mantle as well as the tectonic and seismic activity of the lithosphere. Moreover, geophysics in the wide sense includes many branches of knowledge: volcanology, seismology, geodesy, geochemistry, geomorphology, paleontology, stratigraphy, structural geology, engineering geology, and sedimentology. The degree of knowledge of various geosciences objects varies greatly. The most studied are the processes occurred in the atmosphere, hydrosphere and lithosphere.

The huge mass inside our planet remains poor or completely unexplored. In this situation, it is necessary to ask geophysicists the simplest questions such as: What is the physical principal of 4-layer Earth’s core-mantle models? Is the size of the planet constant? Why is the pressure under the lithosphere higher in the southern hemisphere than in the northern hemisphere? Why did the process of nucleating elements inside the planet stop precisely on iron and nickel? How did elements with an atomic number higher than iron form on our planet? Why continents suddenly break up and start to drift? Is it possible to record a ⁴⁰K, ²³⁵U layer inside the Earth by measuring a geoneutrino? Unfortunately, the geophysics has not answers on these simple questions.

Therefore this study will be interesting not only to specialists who are investigated in the fate of dinosaurs, but also will be interesting to geophysicists, volcanologists, seismologists, astrophysics, paleontologists and other specialists, including climatologists, who could not be successful to reconstruct the paleoclimate. In this sense, this study is fundamental.

Due to our sea level data in this study is significantly above the values obtained by other researchers groups, in this section we presented and discussed the results of previous studies.

In the Cretaceous-Tertiary, the higher sea-levels above the present-day value are: 225 ± 42 m at 82 Ma after EXXON Petroleum Company [73] ; 361 m at 84 Ma after [4] ; 266 m at 91 Ma in [74] [75] ; 242 m at 86 Ma after [5] ; 79 m at 53 Ma by Miller et al., 2005) [76] and in the range of 85 to 270 m in the Cretaceous period (~145 to 65 Ma) after Müller et al. [77]. Although there is consensus concerning on the crude shape of the curve with two maxima in the Cretaceous-Tertiary and the Ordovician-Silurian, the magnitude of the fluctuation is controversial (see various models and references in [78] [79]. Later the Cretaceous eustasy sea levels maximum was updated in 2014 by Haq [80]. The average sea levels throughout the Cretaceous remained higher 75 - 250 m above than the present day mean sea level. Sea level reached two maximum, the first was in early Barremian (~160 - 170 m) and the second (~240 - 250 m), the highest peak of the Cretaceous, was in the earliest Turonian.

As it is well-known, Vail et al. [73] divided sea-level depositional sequences temporally into six orders ranging from tens hundreds of millions years (first- and second-order) to tens of thousands years (sixth order). First- and second-order sea level sequences were ascribed to tectono-eustatic changes in the global ocean volume, while from fourth-order to sixth-order sea level sequences were attributed to climate change within the Milankovitch frequency band. However, the third-order sea level sequences, assigned to time intervals of ~ 0.5 to 3 Ma, were interpreted as the result of climate or tectonic forcing.

According to Miller et al. [76] the eustasy changes in the global sea level happened due to changes in the water volume in the ocean or due to changes in the volume of ocean basins. Thus the water-volume changes are dominated by growth and decay of continental ice sheets, producing high amplitude, rapid eustatic changes up to 200 m. Other processes that affect water volume occurred at high rates and low amplitudes (5 - 10 m): desiccation and inundation of marginal seas, thermal expansion and contraction of seawater, and variations in groundwater and lake storage. Changes in ocean basin volume are dominated by slow variations in sea-floor spreading rates or ocean ridge lengths in 100 to 300 m amplitude. Variations in sedimentation cause moderate amplitude up to 60 m. Thus, an increase in the sea level of more than 300 m is not expected, see also [75] [81] [82].

The exception is Carter [78] and Watts [83], in which it was noted that on a scale of 5 - 100 Ma, the Phanerozoic sea-level cycle associated with 2nd-order sea-level fluctuations could reach 5000 m as a result of thermo-tectonic subsidence on the selected sites of the terrestrial surface. In [84], it is paid the attention to possible influence of Middle Ocean Ridge Basalts (MORBs) and Large Igneous Provinces (LIPs) on the sea level and it was cautiously suggested of 500 - 1000 m range of sea level values.

Thus, without having physical mechanism of a great sea level lifting, the researchers rather carefully express opinions about more than 300 m sea level.

As it is well-known, the crust of our planet consists of light elements, since heavy elements sank down in the melt of the magma. Thus, the presence of heavy elements in the center of the planet, such as Th, U, Pu, is entirely acceptable. Note that an idea about an existence of georeactor was discussed after Kuroda 1960 [85]. Also, the presence of Th and U heat layers in the planet center, natural nuclear georeactor and about a thermal convection in the outer core were widely discussed in the serial studies by Herndon and colleagues [86] - [91]. The possibility of natural reactor presence and possible nuclear reactions in Mars were discussed by Brandenburg, in [92].

Herndon in [86] it was demonstrated the feasibility of planetary-scale nuclear fission reactors as energy sources for the giant outer planets, three of which radiate approximately twice as much energy as they each receive from the Sun. In [87] it was written about the feasibility of a planetary-scale nuclear fission reactor in the center of the Earth as the principal energy source for the geomagnetic field and as a contribution of energy source for other geodynamic processes, such as plate movement. Herndon, in [88] suggested that an U driven georeactor with thermal power < 30 TW presents in the Earth’s core and it is confined in its central part within the radius of about 4 km. In [90] it was pointed out that the georeactor numerical simulation results and the observed high ³He/⁴He ratios measured in Icelandic and Hawaiian oceanic basalts indicate that the demise of the georeactor is approaching. Herndon has proposed that a large drop of uranium has been collected at the center of the Earth, forming a natural 3 - 6 TW breeder reactor, so in this case, nuclear fission should provide the energy source for terrestrial magnetic field, a contribution to missing heat, and the source of the anomalous ³He/⁴He flow from the Earth. The results of numerical simulations of a deep-Earth nuclear fission reactor demonstrated that ³He and ⁴He could be produced by the georeactor [89].

Other aspects of natural nuclear georeactor were investigated in several studies by different research groups [93] - [106]. According to the published studies, a natural georeactor probably exists at the different deep-earth locations, including the center of the core [87] [89] [107] [108] ; on the inner core boundary [100] [101] and on the core-mantle boundary [102] [104].

Also it is well-known the electrons antineutrinos that would be emitted from such hypothetical georeactor have energies above the end-point of geoneutrinos from “standard” natural radioactive decays.

The main reaction of geoneutrino (antineutrino, υ ⌣ e ) registration from natural sources is the inverse beta decay reaction:

υ ⌣ e + p → e + + n ,       Q = 1.806   MeV (1)

Using the registration of geoneutrinos, it is possible to determine a part of the terrestrial heat flux from the radioactive elements (²³²Th and ²³⁸U). It will permit to obtain the vertical distribution of these radioactive elements inside of the Earth and, accordingly, to answer the question about presence and power of a natural nuclear reactor in the center of the Earth. Details of two liquid-scintillator neutrino experiments, such as KamLAND in Japan and Borexino in Italy, in which the geoneutrino signals were measured, could be found in [100] [103] [109] - [122] and in many other publications.

However, there are some difficulties in an implementation of geoneutrino measurements, namely, a small number of inverse beta decay events were recorded per year; usually it was recorded only several hundred events per year. Therefore, the Borexino and KamLand geoneutrino experiments should continue for at least 10 years. In addition, geoneutrino measurements are hampered by a strong background from nuclear reactors and nuclear test sites, by the effect of neutrino-antineutrino oscillations in processes occurring on the Sun, as well as by the burst of supernova stars. Note, that at the small statistics it is difficult to determine stream directions of geoneutrinos from various sources.

Besides the experimental difficulties of detecting the ²³²Th and ²³⁸U antineutrino spectra, for calculation of the natural georeactor power it is necessary to known information about crust thickness (upper, middle, lower) and information about the vertical distribution of ²³²Th and ²³⁸U radionuclides inside the Earth. Usually to calculate neutrino fluxes from the Earth’s interior the Bulk Silicate Earth (BSE) model or the Preliminary Reference Earth Model (PREM) of the Earth’s structure were used. However, the vertical distribution of radionuclides and the model of the Earth structure also cause the numerous discussions, which significantly complicate the processing of the obtained spectra.

Thus, for a long time, it was believed that the antineutrino detection provides a sensitive tool to test the natural georeactor hypothesis; however, the difficulties with recording and processing spectra were appeared as underestimated. Also we will remind that according to geophysical data, the heat flux from the depths of our planet is equal to 44 TW by Pollack et al. [123], 46 ± 3 TW by Jaupart and Mareschal [124] and 47 ± 2 TW by Davies and Davies [125]. These values are much higher than the values obtained in geoneutrino experiments.

Further in [109] and in [110], it was pointed that a detailed analysis excludes a natural reactor producing more than about 20 TW, however, based on the same KamLAND and Borexino Experiments data. In 2015, the researchers of Borexino Collaboration [121] informed that the model-independent analysis yields a radiogenic heat interval that is equal to 11 - 52 TW (69% C.L.) for U and Th radionuclide decay, which be compared with the global terrestrial power output of 47 TW. Later in 2016 Borexino Collaboration restrict the radiogenic heat production for U and Th between 23 and 36 TW, but they set an upper limit for a 3.4 TW georeactor at 90% C.L. or 4.2 TW at 95% C.L. These estimations were based on the statement that geo-neutrinos are produced by the decay of radioactive isotopes present in the crust and the mantle of our planet.

However, early Rusov and colleges in [100] [103] already shown a presence of slow nuclear burning on the boundary of the liquid and solid phases of the Earth’s core with georeactor of 30 TW. Therefore, Rusov and colleagues partially confirmed the theory of the author, presented in this study.

The significant difference consists of the following: ⁴⁰K and ²³⁵U fuel layers cannot be determined by using inverse beta decay reactions. The lower threshold of inverse beta decay reaction is equal to 1.806 MeV, while the upper boundaries of ⁴⁰K and ²³⁵U geoneutrino spectra are below this value. Thus, the ⁴⁰K yield is equal to 1.311 MeV, see Equation 2:

K 40 → C 40 a + e − + υ ⌣ e ,       Q = 1.311   MeV (2)

Except to neglect the decay chains of ⁴⁰K and ²³⁵U isotopes, the natural reactor power calculation method also neglects the decays of ⁸⁷Rb, ¹³⁸La, ¹⁷⁶Lu, ²³⁹Pu and ²⁴¹Pu.

Discussion sometimes takes forms far beyond the limits scientific knowledge. So the conflict between Herndon, the pioneer of geo reactor studying, and his NSF opponents turned into an open troublesome conflict [126]. In this study, we note that neither Herndon nor his opponents were right. The possibility of registration only minor fuel elements (²³²Th and ²³⁸U) casts doubt on advisability of carrying out the long and expensive experiments such as the KamLAND and Borexino Experiments.

Moreover, the near-surface powerful heat layer ⁴⁰K as well as Equation (2) is the basis for the revision of geophysical theories of subduction and of continents drift, as well as introduces changes in the Darwin’s evolutionary theory.

The theoretical basis of this study is the newly developed theory of the internal structure of the Earth by author in [105]. This theory is based on the buoyancy of the chemical elements, which are a part of the Earth. We will shortly remind an essence of our theory.

As it well-known, the Earth’s crust is composed mainly of oxides of light elements such as SiO₂ (59.7%), Al₂O₃ (15.4%), CaO (4.9%), MgO (4.36%), Na₂O (3.55%), K₂O (2.8%), H₂O (1.52%). The exception makes only oxides of Fe (z = 26) and Ti (z = 22), which content in the crust is equal to FeO (3.5%), Fe₂O₃ (2.6%) and TiO₂ (0.6%), respectively [127] [128]. Heavy elements are absent in the Earth crust in a significant amount. Thus the geoscientists draw a conclusion that the chemical elements were separated in the melted magma and core.

From the point of view of the nuclear science, it is quite natural to assume that the gravitational field of the Earth will separate not only the chemical elements, but also will disunite terrestrial isotopes. Let’s remind that the weighting of the elements inside of the Earth occurs due to the capture of neutrons. Below the equations of slow neutron capture and of chemical element transformation F 56 e → N 63 i are presented:

F 56 e + n → F 57 e + n → F 58 e + n → F 59 e → − β C 59 o + n → C 60 o → − β N 60 i + n → N 61 i + n → N 62 i + n → N 63 i → − β ⋯ (3)

The alternative reactions with involving of the ⁶⁰Fe, ⁶¹Co, ⁶⁴Ni and in the s-process (slow) neutron capture equations are not represented in Equation (3). Due to the s-process, it is possible to explain formation of all elements up to Z = 83. Nuclei with Z, greater than 84, do not have stable isotopes and are radioactive. Therefore, the isotope ²³²Th is formed from the ²³²Pb nucleus as a result of eight consecutive β decays. The initial ²³²Pb nucleus formed in the r-process (rapid) and it has 24 neutrons more than the stable ²⁰⁸Pb isotope. Thus, as a result of slow and rapid processes of neutron capture, there is a formation of heavy elements inside the Earth that are slowly deposited deep into the planet. More details about nucleosynthesis of the chemical elements could be found in several reviews [129] - [133].

The linear distribution of the chemical elements inside Earth at the non-perturbed state of natural terrestrial reactor (“cold” planet), according to elemental buoyancy theory, which was developed by author [105], was presented in Figure 3.

The decay products, which were formed in the Th-U layers, will rise up. With the decay of Th and U elements, an inhomogeneous distribution of decay products is formed, with a predominance of light and heavy decay products. Conventionally, we can assume that the peaks of the decay products correspond to such elements as Sr and Cs. According to the model, the Cs level will correspond to the level of 2860 km, which constitutes the boundary between the outer core and lower mantle; this boundary also called as Gutenberg discontinuity.

Thus, inside the Earth, both the processes of lifting up of light decay isotopes and the sink down of elements, which capture neutrons, will occur. The presence of the ⁴⁰K nuclear fuel layer defines the boundary between upper and lower mantle; the presence of heat-generating Th-U isotopes corresponds to the boundary between the inner and outer core (Figure 3). Therefore the Earth is a system of layers consisting of isotopes of chemical elements, which were vertically selected by the gravitational field of the Earth.

Further, it was noted that due to the existence of the hot ⁴⁰K fuel level and shallow convection in the upper mantle, the theory of subduction and continental

drift should be revised. At the period then planet was hot, the viscosity of magma in the upper mantle was low and the probability of continent drift was more than at present. Now abstract theories of subduction and continental drift will get in the framework of the Elemental Buoyancy Theory, described above, another completely convincing physical meaning.

Thus, we have led the reader to understanding of the fact that the standard geological table is a log of the registration of galactic events, which have passed during the existence of the Earth. Such treatment is completely new; no one has previously suggested reconstruction of the structure of the galactic arms by using terrestrial geological sedimentary.

Based on the model investigated in [105], at the period, when natural reactor was operated in a cool, unperturbed mode, the small convective processes dominated inside the magma and core, at which the molten masses do not leave the localization zones. The presence of nuclear fuel or product of nuclear decay at the boundaries leads to an increase in temperature and forms a thermocline or thermopause. In particular in the upper mantle in the convective process, only light chemical elements with atomic numbers up to potassium will be involved. This explains why silicon and sulfur compounds currently dominate in volcanic plumes.

A similar thermocline will be formed at the boundary of the inner and outer core, where the Th-U layers are located. As it is known, the U decay products mainly include chemical elements, which distribution has two maxima, separated by a minimum at the level of elements with atomic numbers, equal to ~50. Decay products, warmer and lighter than U, will rise approximately to ¹³⁷Cs level. Since ¹³⁷Cs also will form a thermocline, it will prevent the rise up of light decay products, such as ⁸⁵Sr-⁹⁰Sr. In [105], it was noted that on seismograms the Sr signal is absent. Thus, the Sr element and its isotopes can be used as a marker of perturbation of the terrestrial reactor (“hot” planet) and as a marker of sea level changes. Also in [105], it was noted that the isotopes ⁸⁵Sr-⁹⁰Sr had to be formed, but this level of discontinuity has not been presented in the seismic records. This means that the Cs layer forms a thermopause, which prevents the light decay elements of the Th-U decay lift up above the Cs layer. The shallow convection process places inside outer core. However, any external disturbance can lead to the destruction of this unstable equilibrium and light decay isotopes such as Sr can lift up in the lower mantle. Further, these isotopes will get into the Middle Ocean Ridge Basalts (MORBs) and into the Large Igneous Provinces (LIPs).

Also the sharp changes in δ¹³С content during the period of mass extinctions were considered in many studies, including for the Early Jurassic mass extinctions, see e.g. [134] - [138]. However, perturbation in the carbon-isotope record recovered very quickly; therefore it cannot be directly connected with the process of extinction of biological species and with process of Jurassic oceanic anoxic events. It is more likely that δ¹³С changes are determined by tectonic processes and astronomical processes. It is necessary to specify that changes in the isotopic composition are a common process at any operating mode reactor changes.

Therefore, it is possible to make an assumption that changes happened in operating mode of CNO nuclear cycle in the upper mantle at sharp activation of the ⁴⁰K layer during period of galactic gravitational compression. Thus, sharp changes in δ¹³С values should correlate with formation LIP, SO₂ volcanic emissions and lift up the heavy chemical elements from depths to upper mantle. The relations between mass extinction, δ¹³С, ⁸⁷Sr/⁸⁶Sr and δ³⁴S were presented in the reviews [139] [140]. Any above isotopes, it is possible to use for the characteristic of a condition of the terrestrial reactor, but we have chosen the ⁸⁷Sr/⁸⁶Sr ratio as a marker of a hot planet. The deep convection, taking place inside the hot planet, involves layers of upper and lower mantle.

The ⁸⁷Sr/⁸⁶Sr ratio during the Phanerozoic, based on [141] [142], were shown in Figure 4 as blue and black lines, respectively. Additionally, the bold color lines indicate the stages corresponding to Triassic and Jurassic. The correlation between the reduced values of ⁸⁷Sr/⁸⁶Sr, corresponding to periods of raised tectonic activity, and the Scutum-Crux and Norma Arms is observed. The more information about the revised ⁸⁷Sr/⁸⁶Sr ratio during the Jurassic could be found in [143], also see details in [144] [145] [146].

As it is known, at gravitational compression the nuclear reactions of hydrogen

and helium burning are possible. The nuclear reactions of burning hydrogen are the next:

H 1 + H 1 → H 2 + β + + ν H 2 + H 1 → H 3 e + γ H 3 e + H 3 e → H 4 e + 2 H 1 S u m m a r y : 4 H 1 → H 4 e + 2 β + + 2 ν + 26.7   M e V (4)

As a result, there is a full combustion of hydrogen and its transformation into helium [130]. A feature of helium combustion reactions is that after main reaction, when two 4He nuclei merge, the second reaction occurs with the formation of an unstable ⁸Be nucleus. However, due to the high density of ⁸Be nuclei (usually at M > 0.7 Msun), before it again will break up on two α-particles, it has time to interact with another ⁸Be nucleus. The process, so-called as the “triple” α-process, occurs with the formation of an excited ¹²C isotope:

H 4 e + H 4 e + H 4 e → B 8 e + H 4 e → C 12 + γ ,       Q = 7.16   M e V (5)

The process of fusion of two nuclei is schematically presented on plate in Figure 5.

With dense spherical packaging approach, the size of chemical elements (D) is associated with atomic numbers (A) as:

D A = r 0 A 1 / 3 (6)

Assuming that in the depths of the Earth’s interior the element packing is spherical, the ratio of atomic sizes before and after nuclear fusion is equal to D_A/D_2A = 0.63, where D_A and D_2A are the diameters of the elements with atomic numbers A and 2A. Also, the dependencies of the degree of compression on the

element number at the fusion of ¹H and ³He/⁴He are shown in Figure 5.

Thus, unlike metals, liquids or gases, at heating the nuclear substance will be compressed and, conversely, at cooling the nuclear substance will be expanded. Note, that applied to the Sun, the effect of expansion is well-known [147]. The Sun in the future will expand and absorb the planets located not far from the star. Fundamentally new is the application of this principle to the Earth’s natural reactor.

Thus, in this study in some sense, we return to the ideas of Roberto Mantovani and Samuel Warren Carey [148] - [152] in which early the possibility of Earth expansion was considered.

Initially soon after the formation the Earth was smaller in size, was hotter and was almost completely covered by the waters of the world ocean. Therefore, it is not surprising that the primary forms of life have arisen in the depths of the ocean. Secondly, as the Earth cools slowly, there is its gradual expansion. Accordingly, the area of the world ocean has decreased and the land area has increased. Hence, the statement that in the process of evolution the reptiles come out from ocean and clime up to land should be considered inconsistent with reality. Those reptiles, which did not have the opportunity to migrate into the depths of the ocean, were forced to adapt to life on land. Note that the process of changing the level of the ocean due to the expansion of the Earth is not monotonous and have its inner (reactor), planetary and galactic episodes. As the question of the Earth expansion as a result of its cooling is obvious, in this study we will focus on the Earth compression and on the worldwide floods resulting from galaxy compression.

First we will show in Jurassic the correlation between mass extinction of land dinosaurs, blossoming of sea fauna and terrestrial nuclear reactor activation. The temporal distribution for Saurachia and Ammonitida, obtained by the PaleoDB, was presented in Figure 7(a) & Figure 7(b). The reducing segment in the ⁸⁷Sr/⁸⁶Sr ratio in Figure 7 is related to the period of terrestrial nuclear reactor activation and to the perturbations in the mantle-core structure (“hot planet”).

As one would expect the amount of Ammonitida fossils is in an antiphase with extinction of land dinosaurs. In the small warm seas which occupy the big areas during the period of “hot planet”, the Ammonitida have quickly reproduced at III, IV, VIII and IX of Jurassic stages. Note that the sharp dinosaurs’ extinction happened at I - III stage (Hettangian-Pliensbachian), which is corresponding to the first Jurassic flood. Below for simplicity of consideration at the altitude analysis, we will be limited to this period, 201.6 - 190 Ma.

It is of interest to investigate in detail the result of statistical analysis of the altitude distribution during the first Jurassic flood, at I - II stages (Hettangian and Sinemurian). The fossil coordinates have been projected on the ETOPO1 map. The received values of Ammonitida, Pectinida, Nautilida and Saurischia altitudes are presented in Figure 8.

The distribution shows a bimodal character in Figure 8(d), which means either the presence of two dinosaurs’ populations or the altitudinal migrations of dinosaurs. Below it will be shown that such a distribution is a result of dinosaurs’ migration during the global floods. As it is possible to see from this distribution (Figure 8), the significant part of dinosaurs escaped from the flood at altitudes of 1600 - 1800 m. Also from marine distributions it can be seen that sea mollusks at the first Jurassic worldwide flood lived at altitudes of 400 - 600, 1400 - 1600 and 2200 - 2300 m, Figures 8(a)-(c). Thus, the structure at the first Jurassic worldwide flood period was heterogeneous and there were at least three flood waves with heights of 400, 1600 and 2200 m. Note that heights of ~1600 m presented in both Ammonitida and Pectinida distributions (Figure 8(a) & Figure 8(b)).

Since the maximum for shallow-water Pectinida practically coincides with the maximum of Saurischia, it is possible to assert that most dinosaurs preferred to live at the sea beaches. At the same time, of course, it is impossible to exclude possibility that all these dinosaurs died at the sea shore due to the giant tsunami caused by tectonic changes. The sea sediments and mudflow, raised by the tsunami, could create favorable conditions for the conservation of dinosaurs’ fossils and marine inhabitant shells.

Also, as it can be seen from the Ammonitida distribution, there is still a maximum at 2200 m, while the corresponding maximum in the dinosaurs’ distribution is practically degenerate. At these heights, a minor maximum of double-clam mollusks was observed. Therefore, the dinosaurs either could not climb to heights of 2100 - 2200 m, or they did not have time to do it, while the well-floating Ammonitida could easily be pushed up at such altitudes by water streams. Further, to show that the given estimations of sea level changes and accordingly of reduction of planet radius are not artifacts it is of interest to investigate the spatial distribution of dinosaurs’ fossils and to define the circumstances of their death.

The Saurischia migration path near the east coast of South Africa and fossil locations are shown in Figure 9(a) & Figure 9(b). In both cases, the dinosaurs moved in the canyon, first to the west, then turned to the south-east and climbed up on the plateau. According to Figure 9(a) the maximum height, in which fossils were found, is equal to 1323 m. At high altitudes, the remains of Thecodontosaurus and Massospondylus have been identified. The approximate dating of dinosaurs’ death is ~195.7 ± 5.4 Ma.

Further, according to Figure 9(b), in second area in South Africa the maximum altitude, where the fossils were found, is equal to 2032 m. At high altitudes, the remains of the Massospondylus, Melanorosaurus and Syntarsus have been identified. The approximate dating of their death is ~195.8 ± 4.9 Ma. In spite of the fact that both areas separated from each other on considerable distance the general behaviors of dinosaurs are similar. It is expressed in aspiration to rise above 1400 - 1800 m.

The analysis showed that dinosaurs repeatedly tried to escape in isolated lagoon systems that were located in the middle of the world-wide oceans. An example of such a lagoon system is the mountain ranges nearby the west coast of North America. The distributions of Saurischia fossils in this lagoon at the first (201 - 190 Ma) and at the second Jurassic flood (168 - 160 Ma) are shown in Figure 10(a)

and Figure 10(b), respectively.

During the first Jurassic flood, dinosaurs climbed up into canyon to the high– mountainous lagoon, which is the ideal place to preserve the population, Figure 10(a). The storms of the Jurassic World Ocean did not reach these places; the waves were breaking on the ridges surrounding the lagoon. The lagoon was connected to the World Ocean through two narrow canyons, which mitigated changes in sea level. At high altitudes, the remains of Kayentapus, Otozoum and Dilophosaurus were identified. The maximum height for finding tracks and fossils is equal to 2001 m. The approximate death date is ~193.7 ± 6.2 Ma. The Saurischia tracks in lagoons are clearly visible, Figure 10(a).

A similar picture of the Saurischia distribution in the lagoon at second Jurassic flood, 168 - 160 Ma is shown in Figure 10(b). The difference between these two analyzed cases is that some of the remains are found outside the lagoon. In the second case, the Saurischia tracks in lagoons are not clearly recognized, Figure 10(b). For obvious flooding, the spatial distribution of Ammonitida and Pectinida is additionally shown. These mollusks formed a characteristic reef on the north-western side of the lagoon. Ammonitida fossils were found at heights of 2160 m (164.8 Ma). The Pectinida fossils were located at altitudes of 2394 m (167.2 Ma), 2040 m (164.5 Ma) and 1440 m (160.4 Ma). The maximum height of the Saurischia remains is equal to 1952 m. At high altitudes, the remains of the Therangospodus, Megalosauripus, Carmelopodus and Brontopodus have been identified. Approximate dating death is ~161 Ma. Thus, the sea level during the second Jurassic flood was significantly higher than that in first Jurassic flood (1200 - 1800 m).

Therefore, the dinosaurs of North America successfully escaped in mounted lagoon during the all Jurassic period, in other words this lagoon is a kind of the “Noah Ark” preserving the variety of terrestrial species in North America during the worldwide floods. Contrary to successful escape events the unsuccessful attempt of dinosaurs escape occurred on the east coast of the North America (Figure 11).

In this case, the dinosaurs’ trails are well-known due to their prints left on the ground. The height of these hills is approximately 600 - 800 m. According latitude-longitude projection on ETOPO1, the maximum height for footprints on trail was recorded at altitude ~144 m. The soil composition at the excavation places was mixing of sandstone, mudstone and “shalle” (pebble). The approximate dating of the tracks and remains is equal to ~197.3 ± 4.0 Ma, which corresponds to the initial period of the first Jurassic worldwide flood, which had 400 m of altitude. As the hill heights are insignificant, it is not surprising that dinosaur fossils, dating at subsequent Jurassic periods, were not found in this region.

Summarizing all, we can conclude about the transience of the flood processes, after which the system slowly returns to its original state, that is, we are dealing with a pronounced sawtooth-shaped process. The maximum length of migration routes is equal to ~700 km. Assuming that dinosaurs moved at a speed of about

1 - 3 km/h, we get that the first Jurassic flood fill huge area in a record short period in range of 10 - 30 days. Note that the Sun and our planet Earth in the future will entered into the zones of gravitational compression, therefore it is necessary to plan in advance creation of galactic rescue bases. North America and Tibet plateaus, rising above sea level to altitude more than 2500 m, are an ideal place to create such a galactic rescue bases.