The Earth and the Sun. An Overview of Two Books
Joseph J. Smulsky
The paper provides an overview of two books "Future space problems and their solutions" and "New Astronomical Theory of the Ice Ages", which were published in 2018. Scientists' opinions about books and also the books’ contents and their main results are given. In the first book, the stability of the Solar system by the example of the evolution of the Mars orbit for 100 million years is shown. The optimal motion of the spacecraft to the vicinity of the Sun is considered. The results of an exact solution to the problem of the interaction of N bodies, which form a rotating structure, are presented. It is shown the evolution of two asteroids: Apophis and 1950DA, as well as ways to transform them into Earth satellites. The cause of the excess rotation of the Mercury perihelion is explained. The examples of the simulation of globular star clusters are given. The second book contains the foundations and issues of the Astronomical Theory of Climate Change. As a result of the interaction of the bodies of the Solar System, the parameters of the orbital and rotational motions of the Earth change. This change leads to a change in the distribution of solar heat over the surface of the Earth. In contrast to the works of previous authors, a new understanding of the evolution of the orbital and rotational motions of the Earth has been obtained. In particular, it has been established that the Earth's orbit and its axis of rotation precess relative to different directions in space. The paper compares the distribution of the amount of heat, i.e. insolation, on the surface of the Earth in different epochs. The insolation periods of climate change are shown. They coincide with the known changes of the paleoclimate. The changes in insolation for different time intervals, including for 20 million years ago, are given. Changes in insolation in the contemporary epoch and in the next million years are also shown. Sunrises and sunsets are also evolving. The paper presents the distribution by latitude of the duration of the polar days and nights in different epochs. Most of the phenomena discussed in the book are related to the Earth, and they are caused by the Sun. Therefore, the names of these two celestial bodies are rendered the title of the paper.
Keywords: interaction, orbits, Earth’s axis, evolution, insolation, climate, ice ages.
Author: Institute of Earth’s Cryosphere, Tyum. SC of SBRAS, Federal Research Center, Russia.
In 2018, two monographs  and  were published, based on the results of solving problems of interaction between bodies and the evolution of their movements. The problems were being solved within several decades. The movement of the bodies of the Solar system, including Earth, Mars, Mercury, asteroids and other bodies, was investigated. Also the movement of the bodies that make up the constellations was investigated. Tasks were considered, both translational and rotational. As a result, a new understanding of the processes occurring on the Earth and in the Solar System has been obtained. This paper has reviewed the results of these two books. The main result is long-term climate change on Earth. It is due to the interactions of the bodies of the Solar system. The climate changes occur on Earth, and the Sun is the main body that changes the Earth’s climate. Therefore, the name of these two bodies is used in the title of the paper.
- The book “Future Space Problems and Their Solutions”
2.1. The book cover
The front cover in fig. 1 shows 50 such planets like the Earth, which are located on a sphere around the Sun. They make a revolution in 1 year and receive as much solar heat as the Earth is now. This is the result of an exact solution to the problem of the interaction of bodies uniformly located on a sphere. It provides unlimited opportunities for the development of Mankind .
The author is presented on the back cover (fig. 1) during a presentation at the Madrid exhibition on May 14, 2011, about the Galactica system, which is intended for solving problems of gravitational interaction of bodies.
Fig. 1: Cover of the book “Future Space Problems and Their Solutions” .
Before the Content of the book are comments of scientists about it.
2.2 The Researcher’s Comments
2.2.1. Academician Vasily M. Fomin, one of the leaders of the Siberian Branch of the Russian Academy of Sciences
The manuscript of J.J. Smulsky's book "Future space problems and their solutions" appeared in my hands. Everyone, sitting at a fire, admired the nature surrounding us on the Earth and the starry sky of our Universe, its abundance of stars, planets, asteroids and comets. And always the question arose, how they move, what interactions between them occur. That is way this book is devoted to the calculation of bodies’ movements in various cases of interaction, based on gravitational forces. The developed system "Galactica" ,  –  has free access and is designed for numerically solving the problems of gravitational interaction of N-bodies. As a result of working with this system, one can study the optimal movements of spacecrafts, the evolution of the Solar system for 100 million years, the influence of the Sun on the Mercury’s perihelion, the evolution of the Earth's rotation axis and much more, which allows us to gain new knowledge about our world.
The system "Galactica" is well described and allows the user with higher education to work with it. All calculation results are displayed for convenience on the display screen. Using this system, the trajectories of asteroids "Apophis" and "1950 DA", their possible approaches to the Earth are studied and various ways of transforming these trajectories are indicated.
Thus, the reader will find a lot of interesting and will be able to explain some of the phenomena that occur in the long-term climate change on Earth, the formation of the future habitat of mankind. Maybe some insights will not coincide with existing conceptions at this time; however, they will serve as a basis for further research in this interesting field of knowledge.
2.2.2 Prof. Vladislav V. Cheshev,
The great philosopher from Tomsk, Tomsk State University, Russia
Today it is customary to speak of the Universe as a product of the Big Bang, which gave rise to the mystery of dark matter, expanding Universe and other tempting but badly comprehend realities. Joseph J. Smulsky gives the chance to feel the harmony of the world, to see it in the image of reasonable classical simplicity peculiar to Isaak Newton's mechanics. The elegance of simplicity and the complexity of the system appear in their mathematically comprehended unity.
2.2.3. Doctor of Technical Sciences Boris F. Boyarshinov,
Insitute of Thermophysics, SB RAS, Novosibirsk, Russia
The problems of gravitational interaction of bodies are extremely diverse. The author considers the situations most significant for the space exploration and to the preservation of life conditions on the Earth-planet. The developed program Galactica can be of interest for professionals and students of the relevant specialties, as for as it allows to apply the received data for the analysis of various problems of celestial mechanics.
2.2.4. Dr. Sc. in Physics and Mathematics Michael P. Anisimov,
Head of the Nanoaerosol Laboratory, SB RAS, Novosibirsk, Russia;
Research Professor, Clarkson University, NY, USA
The book is very impressive. Prof Joseph J. Smulsky considers grandiose projects: the optimal approach of the spacecraft to the Sun, the transformation of asteroids into satellites of the Earth, the creation of new planets like the Earth etc. Without doubt that the space agencies will use the solutions and soft, which are discussed and presented in the book. That will raise the computation culture in the area under consideraion.
2.2.5. Doctor of Technical Sciences Nikolay E. Shishkin,
Insitute of Thermophysics, SB RAS, Novosibirsk, Russia
Environment and human well-being are due to many factors, including the influence of movements of cosmic bodies both in the Solar system and beyond. The proposed monograph provides an excellent opportunity to predict the motion and rotation of planets, to optimize the flights of spaceships by means of a sufficiently simple to assimilate the program Galactica. A lot of worries are caused by possible collisions of comets and asteroids with the Earth. The examples of the most dangerous two asteroids show the capabilities of the Galactica program to calculate their evolution. These opportunities are wider than they are presented in the book. But even what is given impresses with its multifaceted tasks and their scale.
Qverview of Galactica Applications
1.2. Differential equations and method of their solution
1.3. The evolution of the Solar system for 100 million years
1.4. The compound model of the Earth's rotation
1.5. The compound model of the Sun’s rotation
1.6. The motion of the asteroid Apophis
1.7. Optimization of passive orbit with the use of gravity maneuver
1.8. Multilayer ring structures
1.9. The free-access Galactica system
Optimization Of The Mission To The Sun
2.2. Dependence of approaching the Sun on the launching velocity
2.3. Optimal start for approaching the Sun
2.4. Optimal gravitational maneuver for approaching the Sun
Multi-layer Rotating Structures
3.2. Statement of the problem
3.3. The forces acting between the bodies
3.4. Motion equations of the rotating structure
3.5. Solution of the equations
3.6. Examples of rotating structures
3.7. Stability problems for examined structures and their applications
Dangerous Asteroids And Their Possible Use
4.2. Problem statement
4.3. The evolution of the planets’ orbits
4.4. Preparation of initial data of asteroids
4.5. Apophis’ encounters
4.6. Apophis orbit evolution
4.7. Influence of initial conditions
4.8 Apophis’ trajectory in the vicinity of Earth
4.9. Possible use of asteroid Apophis
4.10. Asteroid 1950 DA approaches to the Earth
4.11. Evolution of the 1950 DA orbit
4.12. The trajectory of the 1950 DA in the encounter epoch
4.13. The asteroid 1950 DA and Earth-bound satellite
The Mercury’s Perihelion Precession
5.2. Observed motion of the vernal equinox
5.3. The relative velocity of the perihelion on the observations
5.4. Absolute velocity of the perihelion on observations
5.5. Computing Mercury’s perihelion precession
5.5.1. Gravitational tugs of the other planets
5.5.2. Influence of finite propagation speed of gravity
5.5.3. The precession of the perihelion under action of the Sun
5.5.4. The compound models of the Sun's oblateness
5.6. New sources of the Mercury’s perihelion precession
The Main Problems of Long-term Climate Change
6.2. The theory of the Earth insolation
6.3. The disadvantages of the previous theories
6.3. The solution for the orbital problem
6.4. The research into the rotational motion
6.5. The compound model of the rotational motion of the Earth
The Astronomical Theory Of Climate Change
7.2. Geometric characteristics of insolation
7.3. Insolation of the Earth
7.4. Evolution of orbital motion
7.5. Evolution of rotational motion of the Earth
7.6. Evolution of insolation
7.7. Checking reliability of the Earth rotation problem solution
New Insolation And Paleoclimate Of West Siberia
8.2. Evolution of obliquity and insolation at 65°N
8.3. Latitudinal change in insolation
8.4. Evolution of insolation in other latitudes
8.5. Insolation and the most recent glaciations in West Siberia
8.5.1. The most significant recent glaciations
8.5.2. The Last Glacial Maximum
8.5.3. Warm interglacial
8.5.4. Optimum insolation in the Holocene
8.6. Insolation periods of climate change
8.7. Insolation periods and paleoclimate classifications
Space Mechanics and Future Mankind Progress
9.2. Axisymmetric interactions of N-bodies
9.3. Multi-layer rotating structures
9.4. 3-D structures of N bodies
9.5. Simulation of globular star clusters
9.6. Formation of the future habitat for mankind
The Manual of the Galactica System
10.2. Essentials of Galactica
10.3. Differential equations of motion
10.4. Solution Method
10.5. Input file galacf26.dat
10.5.1. General description
10.5.2. Interaction region parameters
10.5.3. General problem solution parameters
10.5.4. Parameters for trajectories output
10.5.5. On-screen display parameters
10.6. Initial condition and integration data files
10.7. Trajectory files 1t, 2t… traekt.prn
10.8. Close-approach file Dice.dat
10.9. On-screen display
10.10. Examples of input files galacf26.dat
10.11. Initial condition files for several solved problems
10.12. Executable files
10.13. Generation of initial condition file
10.14. Problem solving procedure
Appendix 1. MathCad program for initial condition file
Appendix 2. Modes of display in the standard graphics of FORTRAN
Appendix 3. FORTRAN text of the Galactica program
Index of Authors.
After the Contents, there is a preface to the book, in form of the Foreword, written by Vitaly F. Novikov, a professor at Tyumen Industrial University. He is familiar with the author’s works because for several decades Novikov V.F. repeatedly was a reviewer of Smulsky’s monographs.
Over a period of fifty years, Professor Joseph J. Smulsky has been successively developing and bringing into fruition a novel method for reality cognition. He called his method a hypothesis-free method. The method was described in a book titled “The Theory of Interaction”, issued in 1999 . According to J.J. Smulsky, the action of one on another body consists in that the second body starts moving with acceleration. Mankind has developed mechanics, a science in which the interactions are measured as forces and the acceleration of a body is expressed in terms of the force and the relation between this body and some standard body called mass. By dividing the force by mass, we obtain the body’s acceleration. All forces and masses can be measured, and both the interactions among the bodies and the bodies’ motions can, therefore, be calculated.
In this approach, no assumptions or hypotheses about the nature of the interactions are necessary: it does not matter why the interactions exist, how they are transferred, and so on. That is why many theories such as the gravitational theory, the space-time theory, and the theory of relativity turn out to be unnecessary structures.
In the present book, the hypothesis-free method was employed to study the motion of space bodies. For calculating the interactions, J.J. Smulsky has developed software Galactica ,  – , which allows one to solve, with high accuracy, problems that could not be solved previously. One of such problems, for instance, is the evolution of the Solar system for 100 million years.
Also, J.J. Smulsky has analytically exactly solved some previously unsolved N-body problems, including the following ones: what motions are exercised by bodies on a ring ? Can a system of bodies involving several of rings rotate as an entity ? Those and other exact solutions have allowed J.J. Smulsky to solve several previously unsolved problems: about an optimum mission to the Sun , about the motion of Mercury’s perihelion , about the evolution of the Earth axis -, and so forth.
The solutions obtained provide new impressions of our world, explaining why the climate on the Earth suffers changes and what achievements are possible for mankind in the future.
In his new book, Smulsky explains the solution to these problems, provides detailed proofs, solution methods, computer programs and initial data. The whole toolkit allows young researchers to retrace the path of the pioneer of the problems presented, as well as to attempt the solution of unsolved other problems.
Vitaly F. Novikov,
Doctor of Physical-Mathematical Sciences,
Professor, Tyumen Industrial University, Russia.
2.5. Some book’s findings
Periodic changes are taking place on the Earth: day gives way night, and winter comes after summer. There are also more significant changes when the climate gets colder, and Ice Ages begin, and then followed by warm periods. We know that the change of day and night times on the Earth is caused by the rotation of the Earth around the Sun, and the change of seasons is due to the Earth's orbital motion. Apparently, the longer-term changes are also due to the changes in the motion of Solar-system bodies. After several years of working on that problem, the author has developed the following research plan (fig. 2). First, it is necessary to study the evolution of the orbits of the Solar system bodies and, then, the evolution of the axis of Earth’s rotation. At the third stage, it is required to investigate the evolution of insolation, and at the fourth, to compare the change in insolation with the change in natural processes. Then, at the fifth stage, we can finally tell how natural processes will develop on the Earth.
Fig. 2: The problems and relationships between them in the study of the evolution of the Earth’s Cryosphere as a result of the interaction of Solar system bodies .
Fig. 3: The evolution of the Mars orbit in the span of 100 million years: T is the time in millions of years into the past from the epoch 30.12.1949. wpm = 1687 "/century is the average angular velocity of the perihelion .
The graphs (fig. 3) show the changes of the four parameters of the Mars orbit: the eccentricity e, the angle of the orbit’s node ϕΩ, the inclination angle of the orbit plane i, the perihelion angle ϕp, and the perihelion’s angular velocity ωp over two time intervals, each interval lasting 50 million years. Such parameters are given in fig. 3: Te2 is the second period of eccentricity, TΩ and Ti are the shortest periods ascending node and inclination, respectively, in kyr (thousand years); em, im andϕΩm are the average values of appropriate parameters; Tp is a rotation perihelion period averaged over 50 Myr.
The parameters of the Mars orbit fluctuate within unchanged limits. The graphs showing the data over the two intervals of 50 million years are almost indistinguishable. Thus, the Mars orbit is regular and stable, and there are no signs for the violation of this stability. The same results were obtained for the other planets . These results, therefore, prove the stability of the entire Solar system.
Fig. 4. The angle between the planes of the orbits of Mars and its equator at different initial conditions in the solutions by Laskar et al. : the angle (Y-axis) is in degrees, and the time (X-axis) is in the past 250 million years.
Before the author’s work, the researchers came to the conclusion about the instability of the Solar system and about its possible destruction in one billion years. The graph (fig. 4) by Laskar and co-authors shows that the angle of inclination of the Mars orbit may change irregularly and, depending on the initial conditions, those changes may follow any of the four variation patterns.
Fig. 5. The trajectories and orbits of the spacecraft when launched on January 20, 2001, with different initial velocities v. The flight is passive. After the influence of Venus (when it crosses its orbit), the device goes into an elliptical orbit. 1 is the orbit of Venus .
These graphs (fig. 5) show the optimal trajectories of the spacecraft’s approach to the Sun. The position of the bodies, the Earth, Venus, Mercury, and the Sun, is shown by the moment of launch. The spacecraft is launched against the Earth’s orbital motion, and then it moves with the propulsor off. At the end of the propulsor’s operation, the coordinates and the velocity of the spacecraft are calculated so that it would have passed at a certain distance from Venus, letting the Venus to reduce its speed. Thus, here, a closer passage of the spacecraft to the Sun is ensured in comparison with those that can be achieved using other launching schemes. The time required for reaching the closest distance to the Sun is also shortened. For example, with a launch speed of 15 km/s, this time is equal to 0.25 years, whereas, with other launching schemes, this time turns to be equal to several year.
Fig. 6: Images of axisymmetric multilayer rotating structures as displayed on the computer screen. The structures were obtained by integration of differential motion equations performed using the Galactica program: N2 = 15; N3 = 30; period Prd =1 year; the mass of the central body is equal to the Sun mass. Given in the table are the relative radii of the rings and the masses of one body in the rings .
Here (fig. 6), results of an exact solution of the problem of N-bodies located in N2 layers around a central body are shown. Each of the layers contains N3 bodies, and the whole system rotates as a single entity. In this slide, we see three configurations of such a structure depending on the location of the bodies in adjacent layers. The structures have 15 layers with 30 bodies locating in each layer.
The graphs (fig. 7) show the distances between the asteroid Apophis and the Mercury, Venus, Moon, Earth and Mars for 100 forthcoming years (a), 100 past years (b) and 1000 future years (c) .
At the moment A, April 13, 2029, Apophis will pass at a distance of six Earth’s radii, and there will be no more such approaches.
Therefore, during the execution of this work in 2011, we proposed to transform the trajectory of Apophis into the orbit of an Earth satellite.
Fig. 7: Apophis’s encounters with celestial bodies during the time ΔT to a minimum distance Rmin, km: Mars (Ma), Earth (Ea), Moon (Mo), Venus (Ve) and Mercury (Me); a, b – ΔT = 1 year; c – ΔT = 10 years. T, cyr (1 cyr = 100 yr) is the time in Julian centuries from epoch JD0 (November 30.0, 2008). Calendar dates of approach in points: A – April 13, 2029; B – April 13, 2067; C – September 5, 2037; E – October 10, 2586.
Fig. 8: The trajectories of Earth (1) and 1950 DA (2) in the barycentric equatorial coordinate system xOy over 2.5 years in the encounter epoch of March 6, 2641 (point Ae): A0 and E0 are the starting points of the 1950 DA and Earth trajectories; Af and Ef are the endpoints of the 1950 DA and Earth trajectories; 3 - 1950 DA trajectory after the correction applied at the point Aa is shown approximately; the coordinates x and y are given in AU.
The same studies were performed for the asteroid 1950 DA, weighing 30 times more than Apophis. The asteroid 1950 DA will pass after 600 years at a distance from the Earth no closer than 2 million km. Fig. 8 shows a diagram of transformation of the trajectory of 1950 DA into the orbit of an Earth satellite. 1 is Earth's orbit; 2 is 1950 DA’s orbit. At the point Aa, the asteroid accelerates slightly, and it bypasses the Earth from the night side, point Ab, at a distance of six Earth’s radii. At the latter point, 1950 DA gets decelerated and becomes a satellite in a geostationary orbit. In this orbit, the asteroid will immovably hang over a certain equatorial point of the Earth.
Fig. 9: A comparison of the results of integrating differential equations of the Solar System bodies interaction in a span of 6,000 years with the form of point Sun (left) and with a compound model No. 15 (right) for the secular changes of the Mercury’s orbital elements. The interval between the points in centuries is two for left graphs and one for the right graphs; e is the eccentricity; i is the angle of inclination of the orbit plane to the equator plane of 1950.0; ϕΩ is the angle of the ascending node of the orbit; ϕp0 is the angle of the perihelion; Δa in m is the semimajor axis deviation from the average of 6,000 years, the values in meters; ΔP in centuries is the orbital period deviation from an average of 6,000 years. The angles are in radians, and the time T is in centuries in the past and the future from the epoch December 30.1949. 1 shows the results of the numerical solution of the Galactica program; 2 shows the secular changes according to S. Newcomb, 1895 ; and 3 shows the secular changes according to J.L. Simon et al., 1994 .
The graphs (fig. 9) show the results of calculations of the changes in the parameters of the Mercury orbit due to the interaction of the bodies in the Solar System for two cases: on the left, all bodies are treated as material points, and on the right, the oblate shape of the Sun is taken into account in the form of its 15th model.
Here, points 1 show the calculated parameters for six thousand years, i.e., for spans of ± 3 thousand years from the contemporary epoch; and lines 2 and 3, the results of approximation of the observational data.
The following parameters were considered: the eccentricity e, the inclination angle of the orbit plane i, the angle of the orbit’s ascending node ϕΩ, and the perihelion angle ϕp. All these parameters suffer changes. Since two other parameters, the major semi-axis a and the period P undergo no changes, the oscillations of the deviations of those parameters, Δa and Δp are shown here.
A noticeable difference between the calculated and observational data is only observed for ϕp0. As seen from the right graph, taking into account the oblate shape of the Sun leads to the coincidence between the calculated and observational data.
Thus, the Newton law of gravitation allows one to rather accurately predict the motion of all bodies in the Solar system.
The fig. 10 shows a structure formed by 99 bodies uniformly distributed on a sphere around the Sun. This structure was obtained as a result of the exact solution of the interaction problem of these bodies. The bodies are evenly distributed along the three-turn curve. On the right, these bodies are shown in motion. Each body moves in a circular orbit, and it receives as much heat from the Sun as the Earth presently does. As it is already noted, the results obtained while solving this problem offer unlimited possibilities for the development of mankind.
Fig. 10: A regular spatial structure with N = 100, radius of the sphere R = 1 AU (astronomical unit is equal to 1.496 million km), and period Prd = 1 year and with the central-body mass equal to the Sun mass: a – in 3D coordinate system, the lines at bodies 1 and 49 indicate the velocity vectors; b – image of the structure as displayed on the computer screen after 100 revolutions of a peripheral body around central body in the integration of the differential motion equations (1.3) with the Galactica system; a – projection onto the frontal plane, b – projection onto the horizontal plane.
Fig. 11: The M 53 (or NGC 5024) globular star cluster in the Coma Berenices constellation (https://ru.wikipedia.org/wiki/M_53.
Fig. 11 shows a globular star cluster in the Coma Berenices constellation. The globular star clusters are the oldest formations in the universe. Thousands of stars move in their orbits and never collide. How does this happen?
Fig. 12 shows three models of globular star clusters, which were derived from the spherically distributed solutions. This images are the models appearance after 100 revolutions of each body. After 50 revolutions, the configuration does not change. The models differ in the formation scenarios of the spherical structure. Structure 1 is shown in dynamics in the right bottom. The body 2 has an elongated orbit. In the dynamics, it becomes evident that the bodies move at differently directed velocities. At the beginning of the formation process, some of the bodies suffer collisions, and some of them merge into one body. The bodies become larger; as time passes, the collisions grow rare, and there comes a time when no further collisions occur. Each body acquires its own orbit close to an ellipse. The plane of this orbit intersects the planes of the orbits of other bodies, but the bodies pass the common points of their orbital planes at different places or at different times.
In contrast to the globular clusters shown in fig. 11, in these models (fig. 12) the bodies are concentrated at the periphery. Such globular clusters are also found.
Fig. 12: Spherically distributed structures with N ≈ 1000 after 100 revolutions of the peripheral bodies: 1, 2, and 3 — versions of the initial structures. Bottom right shows structure 1 in motion. The line bars at bodies 1 and 999 and so on indicate their velocity vectors v1 and v999. The distance L to the second body is been diminished in 2 times.
Fig. 13: Multilayer spheroidal model of a globular star cluster: shown are five layers involving 1488 peripheral bodies; a — at the time the structure was formed; b — after the 196th revolution of the first-layer body; the line bars at bodies 1 and 1488 indicate their velocity vectors, v1 and v1488.
The model of a spherical star cluster is developed based on a multilayer structure (Fig. 13). The figure on the right shows this structure after 196 revolutions of the bodies in the inner layer. The structure was formed after the 20th revolution, and during further 180 revolutions it remained almost unchanged.
- THE BOOK “NEW ASTRONOMICAL THEORY OF THE ICE AGES”
3.1. The book cover
Fig. 14: Cover of the book “New Astronomical Theory of the Ice Ages” .
The front side of the cover (fig. 14) shows how the Earth's illumination by the Sun is related to the orbital and rotational motion of the Earth. The Earth rotates around its axis, and therefore there is a change of day and night times on it. Rotation proceeds in a counterclockwise direction. The Earth moves around the Sun in its orbit also counterclockwise. The Sun is located not at the center of the orbit, but in one of the foci of the ellipse. Therefore, in perihelion (in winter, the Northern Hemisphere), it approaches the Sun to the smallest distance Rp, and in aphelion (in summer) it moves away from the Sun to the greatest distance Ra. The axis of the orbit is perpendicular to its plane. The axis of Earth’s rotation is inclined to the Earth axis at an angle of 23.4°. Therefore, on December 21, the Southern Hemisphere is better illuminated by the Sun, and the polar night stands behind the Arctic Circle. This is the winter solstice. On June 21, the Northern Hemisphere is better illuminated, and the polar day is above the Arctic Circle. At the position of March 20 and September 22, both hemispheres are equally illuminated. This is the day of the spring and autumn equinoxes, respectively. This is how the seasons change.
The parameters of the orbital and rotational motions of the Earth at the contemporary epoch create the climate on Earth that we know. However, these parameters undergo changes. The axis of the orbit rotates, or precesses, around the vector clockwise, while the Earth’s axis precesses also clockwise around another vector . Also, the both axes oscillate. Additionally, the orbit, shaped like an ellipse, rotates in its plane counterclockwise, and its shape changes: the ellipse either becomes more elongated or its turns into a circle. All this leads to climate change on the Earth, and it becomes the one we do not know yet. The back cover of the book contains the book’s annotation.
This book examines the factual aspects of the following problem: how our world works, why the observable phenomena occur, and why they evolve. The problems of the evolution of the orbital and rotational motions of the Earth, which underlie the new Astronomical theory of the ice ages, are among the most difficult problems in contemporary science. These problems were solved, and the book presents the obtained results, which can be comprehended by a wide circle of the readers. The comprehensibility of the book has been achieved due to 63 drawings and graphs, which also make the book compact and informative. The material contained in the illustrations far exceeds the text of the book, and it will serve as the basis for further research in various fields of science. For those who wish to understand the origins of the results obtained, to check and verify their validity, and maybe even obtain them themselves, the book is supplemented with the material given in Appendices.
The Evolution of the Orbital And Rotational Movements of The Earth
1.1. The movements of the Earth and their change
1.2. Evolution of the Earth's orbital motion
1.3. Evolution of Earth's rotational motion
1.4. Evolution of the Earth's orbital motion relative to its rotational motion
Insolation of the Earth, Its Structure And Evolution For 200 Ka (thousand Years Ago)
2.1. Evolution of the obliquity and insolation
2.2. Change in insolation along the latitude of the Earth
2.3. Evolution of insolation at other latitudes
2.4. Periods and grades of Earth's climate change
2.5. Insolation periods and paleoclimate
2.5.1. Insolation periods for 50 ka
2.5.2. Insolation periods and existing paleoclimatic classifications
Evolution of Insolation And Earth Movements For 5 Ma
3.1. Evolution of the obliquity and insolation of the Earth for 1 Ma
3.2. Evolution of the Earth's insolation and its movements for 5 Ma
Evolution Of Insolation For 20 Ma
4.1. Evolution of the obliquity and insolation for 20 Ma
4.2. Change in insolation along the latitude of the Earth in the warmest and coldest epochs for 20 Ma
Statistics of Insolation for 20 Ma
5.1. Analysis of insolation at equivalent latitudes for 20 Ma
5.2. Analysis of the summer insolation Qs65N for 20 Ma
Contemporary and Future Changes of Insolation
6.1. A contemporary changes of insolation
6.2. The evolution of insolation in future epochs
Comparisons, Additions, and Verifications
7.1. Comparison of insolation periods of paleoclimate with marine isotope stages
7.2. The duration of the Sun’s phenomena
7.2.1. Evolution of the Sun's phenomena
7.2.2. The annual and daily movement of the Sun in the middle latitudes
7.2.3. The annual and daily movement of the Sun in the polar latitudes
7.2.4. The daily movement of the Sun on the ground
7.2.5. Duration of a light day at the time of the solstice
7.2.6. Moments of the onset and end of the polar days and nights
7.3. Evidence of the reliability of the results
7.3.1. Statements of the problems and their differences
7.3.2. Reliability of the solution of the orbital problem
7.3.3. Reliability of the solution of the Earth’s rotation problem
7.3.4. Reliability of the solution of the third and fourth tasks
7.3.5. The physical cause of the difference between the new theory and the previous one
7.3.6. Final check
Appendix 1. Differential equations of motion of the Solar system bodies
Appendix 2. Evolution of the Earth's orbit over a time interval of 100 million years and the stability of the Solar system
Appendix 3. Differential equations of rotational motion of the Earth
Appendix 4. Evolution of the Earth's rotation axis for 20 Ma
Appendix 4.1. Precession of the axes of the Earth's orbit and rotation
Appendix 4.2. Oscillations of the Earth's rotation axis in the second inertial coordinate system
Appendix 4.3. The statistics of the oscillation of the obliquity ε for 20 Ma
3.4. Some book’s findings
Fig. 15: The Earth's position in its orbit in 2025 at the days of the spring equinox (March 20), of the summer solstice (June 21), of the autumn equinox (September 22) and of the winter solstice (December 21), and the time of its movement in days in spring (92.7 d), in summer (93.7 d ), in autumn (89.9 d) and in winter (89.0 d): is the Earth’s axis of its rotation, and is the vector, relative to which the axis precesses with a period of 25.74 thousand years; is the axis of the Earth's orbit, and is the vector relative to which the axis precesses with a period of 68.7 thousand years .
Consider the periods during which the parameters of the orbital and rotational motions chang (fig. 15). The oscillation periods of the orbit’s axis are: 97.35, 1164, and 2320 thousand years .
Eccentricity is e = (Ra-Rp)/(Ra+Rp). Oscillations periods of e are: 94.5, 413, 2310 thousand years.
The average period of perihelion’s rotation is 147 thousand years.
The orbit’s axis (fig. 15) precesses around the vector clockwise with a period of 68.7 thousand years, and the Earth’s axis also precesses clockwise around the vector with a period of 25.74 thousand years.
The periods of oscillation of the Earth’s axis are: 1 day, half a month, half a year, 18.6 years, and tens of thousands of years.
The amount of heat coming from the Sun to the Earth, i.e. insolation of the Earth, depends on the parameters of its orbit relative to the moving equatorial plane, and is determined by the following parameters: eccentricity of the orbit e, obliquity ε, which is the angle between the axes and (fig. 15), and the perihelion angle φpγ.
Fig. 16: φ-latitude distribution of summer Qs, winter Qw and total, i.e. annual QT insolations for three epochs: 1 is contemporary one; 2 is the warmest; 3 is the coldest epoch for 200 thousand years ago; Qs65N is insolation in GJ/m2 for the summer caloric half-year at the northern latitude of 65°; T, kyr is time in thousand years from December 30, 1949. .
The graphs (fig. 16) show the change in the amount of heat along the Earth's latitude φ from the South Pole (φ = -90°) to the North Pole (φ = 90°) in three epochs over the past 200 thousand years: 1 – contemporary epoch, 2 – the warmest epoch, occurred 31.28 thousand years ago, and 3 – the coldest epoch, which occurred 46.44 thousand years.
The dashed red lines show that from cold epoch 3 to warm epoch 2, the amount of heat per summer half-year (Qs) at the poles had changed two-fold. The range of the variation as we move toward to the equator decreases, and at the equator, the change in heat over the summer half-year is insignificant.
The blue lines show that, during cold epoch 3, the Earth receives a greater amount of heat in winter half-year (Qw) in comparison with warm epoch 2. Warmer winters in cold epochs occur over all terrestrial globe.
The orange lines show that from cold epoch 3 to warm epoch 2, the amount of heat per year (QT) at the poles, as well as in the summer half-year, changes two-fold. However, the range of change rapidly decreases as we move to low latitudes, and at the latitude of 45 , the annual amount of heat does not change. In the equatorial latitudes, the annual heat changes in the opposite way: in cold epoch 3 at the equator it is warmer than in the warm epochs.
So, first, the warming periods and the cooling periods such as the ice ages occur on the Earth at high latitudes. At low latitudes, in the glacial periods, becomes warmer.
Second, the winters during the glacial periods are warmer than those during the warm periods. These are the climates that exist on the Earth, and which we do not know.
Fig. 17: Insolation periods OI, 1I, 2I,…, 12I over the last 200 kyr and their bounds: 1 is mean insolation Qsm; 1t and 2t are first and second boundaries of warm levels; 1c and 2c are first and second boundaries of cold levels; m.w., w. and e.w. are moderately warm, warm and extremely warm levels; m.c., c., and e.c. are moderately cold, cold and extremely cold levels .
The graph (fig. 17) shows the changes in heat over the summer half-year at the latitude of 65 degrees in the Northern Hemisphere Qs65N for 200 thousand years ago. The parameter Qs65N is an indicator of climate change at high latitudes. All the extremes of Qs65N are numbered, the insolation periods 0I, 1I,…12I of climate change are introduced, and their boundaries are identified. Those periods correspond to: 0I – Holocene optimum, 1I – Sartan ice age in Western Siberia, 2I – Karginsky warming, 3I – Ermakovsky ice age. The Ermakovsky ice age is one of the strongest glaciations in the past 200 thousand years.
Fig. 18: Climate change levels for 1 Myr: 1 - average insolation Qsm; 1t and 2t - the first and second boundaries of warm levels; 1c and 2c are the first and second boundaries of the cold levels; m.w., w., and e.w. are warm levels; m.c., c., and e.c. are cold levels; tm.cl. is periods of temperate climate.
These periods also occur in Europe and North America. However, in those parts of the world, they are called differently  – .
Also, three levels of warm climate are introduced (fig. 17): a moderately warm climate, a warm climate, and an extremely warm climate. Similarly, three levels of cold climate are introduced; these are a moderately cold, a cold, and an extremely cold climate.
From the graph (fig. 18) it is seen that, over the past million years, there have been four extremely warm periods and six extremely cold periods. Those periods occur unevenly. During the first 300 thousand years, all strong cooling and warming periods had occurred. And then there were no such strong changes over the next 700 thousand years.
Fig. 19: The evolution of the orbit eccentricity e, the perihelion position angle φpγ, the obliquity ε, as well as the summer insolations Qs65N and I for 5 Ma (million years ago); Tpγt, thousands years, is a moving average perihelion rotation period for 20 thousand years; Tpγm = 21.7 thousand years is the average period of the perihelion rotation.
The graphs (fig. 19) show the changes in the parameters of the orbital and rotational motion of the Earth and insolation during the past 5 million years; those parameters are the eccentricity e, the perihelion angle ϕpγ reckoned from the moment of equinox (fig. 15), the angle ε between the axes and , the summer insolation Qs65N and insolation I in equivalent latitudes. The large and small values of Qs65N are distributed unevenly in time. The large values of insolation I at equivalent latitudes indicate that at the latitude of 65° during those epochs there was more heat than nowadays at the equator. In other words, it is really very strong warming.
Simultaneously, the low values of I suggest that at the latitude of 65° the heat received by the Earth during the summer half-year was lesser than nowadays at the pole. Indeed, these are really ice ages.
Fig. 20: The evolution of summer insolation Qs65N for 20 Ma in the form of 20 graphs for 1 million years (line 1); line 2 is summer insolation Qs65N by Laskar et al. ; Tl is the left reading of the time axis in million years, given in the left corner of each graph; 1t, 2t and 1c, 2c are the first and second boundaries of the warm and cold levels, respectively.
The graphs (fig. 20) show the changes in summer insolation Qs65N during 20 million years in the form of 20 graphs, each graph representing a period of one million years. In the first graphs of the left panel, Qs65N begins on the right for the first million years with the left countdown of time Tl = -1 and, then, it continues beyond the second graphs Tl = -2. The graphs end in the right bottom panel (-19) and (-20), i.e. with the 20th million of years.
The graphs show the unevenness of very strong warming and cooling of the climate. Here, one can find calm time intervals lasting one million years, when there were no strong climate fluctuations. This data is further evidence showing that, on the Earth, there were such climates that we do not know. Previously, all researchers were looking for the periodicity in the onset of glacial periods. Here, no such periodicity is observed. However, these changes occur with strict certainty and inevitability.
The lines 1 of the graphs (fig. 21) show the variation of the obliquity ε, of insolations Qs65N and I over the summer half-year for the forthcoming 200 thousand years. Considering insolation I, one can see that very strong warming will occur in 90 thousand years. The same graphs for the future one million years (fig. 22) show that there will be another strong warming that will occur in 900 thousand years. Thus, the next million years is characterized by smaller climate changes compared with the previous million years (fig. 18).
Fig. 21. The evolution of the obliquity ε, as well as summer insolations Qs65N and I for the span of 200 thousand years into the future. Comparison of new results 1 with the results of previous theories on the example of the works by Laskar et al.  2 and by Sharaf and Budnikova  3: obliquity ε is the angle of inclination in radians of the Earth's equator to the plane of its orbit; Qs65N is insolation in GJ/m2 for the summer caloric half-year at the northern latitude of 65°; I is insolation in equivalent latitudes over the summer caloric half year at north latitude 65°. T is the time in a thousand years from December 30, 1949 .
Fig. 22: The evolution of the obliquity ε, as well as summer insolations Qs65N and I for the span of 1 million years into the future. Comparison of new results 1 with the results of previous theories 2 on the example of Laskar et al . The interval between points is 200 years. The remaining designations see fig. 21.
Fig. 23: The duration of the polar night and the polar day in the Northern Hemisphere in three epochs: 2025, 15.9 and 31.3 thousand years ago. Latitudes from 90° to 66° are shown to the left of the picture, and from 68° to 58° are to the right of the picture.
On the sectorial map of the Arctic (fig. 23), the duration of polar nights is shown against the dark background, and the duration of polar days is shown against the yellow background in the three epochs: in the contemporary epoch 2025 year, 15.9 thousand years ago and 31.3 thousand years ago  – .
In the contemporary epoch (2025), the polar night at the pole lasts 175 days, and the polar day – 191 days. At the latitude of 66 degrees, there are no polar nights anymore, and the polar day lasts 18 days. In the Sartan Ice Age, which occurred 15.9 thousand years ago, at the latitude of 70°, there are no longer polar days and nights. Next, in the extremely warm Karginsky period, which occurred 31.3 thousand years ago, at the latitude of 58°, the polar day lasts 25 days. The 58° is almost the latitude of Tyumen. That is, in an extremely warm period, here in Tyumen, as today at the Arctic Circle, the polar day lasts almost a month, and a short distance to the north we have a polar night. These data are further evidence for the Earth's climates which we do not know.
Fig. 24: Projections of the precession trajectories of the Earth's orbit axis (a) and its axis of rotation (b) for 5 Myr on a plane perpendicular to the vectors and , respectively. The positions of the axes at time points and angles between them: 1S and 1N for T = 0 kyr, ε = 23.443°; 2S and 2N for T = -0.2326 Myr, ε = 30.778°; 3S and 3N for T = -2.6582 Myr, ε = 32.680° .
In the left graph (fig. 24), the red line shows the motion of the orbit’s axis over past 5 million years. The graph is the projection of the axis onto the plane normal to the precession vector . The blue line on the right shows the precession of the axis of Earth rotation projected onto the plane normal to the precession axis . On the same plane, the red line is the projection of the orbit’s axis . The angles ε between the axes and in the corresponding epochs vary in proportion to the distances 1N-1S, 2N-2S, 3N-3S.
The presented data on the variation of insolation are due to the large fluctuations of the angle ε. Those fluctuations seven to eight times exceed the range of oscillations of this angle in the works by previous authors. These exceedances are due to the structure of the precession of the Earth orbit’s axis and the axis of Earth’s rotation revealed in the present study.
Let us return to the problems posed at the beginning of the paper (fig. 2). Then, 20 years ago, I thought: this is how one has to proceed with the work on the five stages. In those times, I would feel happy if I could accomplish at least the first three stages. Today, however, I see that I have all the five stages completed. The cause of the long-term climate fluctuations on the Earth is the interactions between bodies in the Solar system. These interactions lead to the fluctuations of the solar heat falling onto the Earth. The minima and maxima of this heat were identified with high accuracy: over 200 thousand years, accurate to several minutes, and over 100 million years, accurate to several days. These fluctuations are manifested in marine and continental sediments, in the form of relief, and the changes observed in the plant and animal world, including humans. Today the researchers, which work in these areas, that is geologists, geographers, geocryologists, biologists, historians and, even philosophers, face the problem of identifying those manifestations and correlating them to the extremes of insolation. As a result, we will get accurate and consistent knowledge about the past of our planet.
The materials of this work were obtained as a result of research at the Institute of Earth's Cryosphere, Tyumen SC of SB RAS, Federal Research Center for three decades, and in recent years, the research has been carried out under the project IX.135.2.4.
The results are based on the solution of the problems about the orbital and rotational motions obtained on supercomputers of the share use center at the Siberian Supercomputer Center, Institute of Computational Mathematics and Mathematical Geophysics SB RAS, Novosibirsk, Russia.
- Smulsky J.J., 2018. Future Space Problems and Their Solutions. Nova Science Publishers, New York, 269 p. ISBN: 978-1-53613-739-2. http://www.ikz.ru/~smulski/Papers/InfFSPS.pdf
- Smulsky J.J., 2018. New Astronomical Theory of Ice Ages. “LAP LAMBERT Academic Publishing”, Riga, Latvia,. 132 p. ISBN 978-613-9-86853-7.http://www.ikz.ru/~smulski/Papers/InfNwATLPEn.pdf. (In Russian).
- Smulsky J.J., 2017. Advances in Mechanics and Outlook for Future Mankind Progress. International Journal of Modern Education and Computer Science (IJMECS), Vol. 9, No. 1, 15-25.http://www.mecs-press.org/ijmecs/ijmecs-v9-n1/IJMECS-V9-N1-2.pdf.
- Smulsky J.J., 1999. The Theory of Interaction. Novosibirsk: Publishing house of Novosibirsk University, Scientific Publishing Center of United Institute of Geology and Geophysics Siberian Branch of Russian Academy of Sciences. 294 p. http://www.ikz.ru/~smulski/ TVfulA5_2.pdf. (In Russian).
- Smulsky J.J., 2012. Galactica Software for Solving Gravitational Interaction Problems. Applied Physics Research, Vol. 4, No. 2, 110-123.http://dx.doi.org/10.5539/apr.v4n2p110.
- Smulsky J.J., 2012. The System of Free Access Galactica to Compute Interactions of N-Bodies. I. J. Modern Education and Computer Science, 11, 1-20. http://dx.doi.org/10.5815/ijmecs. 2012.11.01.
- Smulsky J.J., 2003. The Axisymmetrical Problem of Gravitational Interaction of N-bodies. Mathematical modeling, Vol. 15, No 5, 27-36.(In Russian).
- Smulsky J.J., 2015. Exact Solution to the Problem of N Bodies Forming a Multi-layer Rotating Structure. SpringerPlus, 4:361, 1-16, DOI:10. 1186/s40064-015-1141-1, URL: http:/ /www.springerplus.com/content/4/1/361.
- Smulsky J.J., 2008. Optimization of Passive Orbit with the Use of Gravity Maneuver. Cosmic Research, Vol. 46, No. 5, 456–464. http://www.ikz.ru/~smulski/Papers/COSR456.PDF.
- Smulsky J.J., 2011. New Components of the Mercury's Perihelion Precession. Natural Science, Vol. 3, No.4, 268-274.doi:10.4236/ ns. 2011.34034.http://www.scirp.org/journal/ns.
- Smulsky J.J., 2011. The Influence of the Planets, Sun and Moon on the Evolution of the Earth’s Axis. International Journal of Astronomy and Astrophysics, Vol. 1, Issue 3, 117-134.http://dx.doi.org/10.4236/ijaa.2011.13017.
- Smulsky J.J., 2016. Fundamental principles and results of a new astronomic theory of climate change. Advances in Astrophysics, Vol. 1, No. 1, 1–21.http://www.isaacpub.org,http:// www.isaa cpub.org/Journal/AdAp.
- Melnikov V.P. and Smulsky J.J., 2009. Astronomical theory of ice ages: New approximations. Solutions and challenges. Novosibirsk: Academic Publishing House. http://www.ikz.ru/~smulski/Papers/AsThAnE.pdf.
- Grebenikov E.A. and Smulsky J.J., 2007. Evolution of the Mars Orbit on Time Span in Hundred Millions Years. Reports on Applied Mathematics. Russian Academy of Sciences: A.A. Dorodnitsyn Computing Center. Moscow. 63 p. http://www.ikz.ru/~smulski/Papers/ EvMa100m4t2.pdf. (In Russian).
- Laskar J., Correia A. C. M., Gastineau M., Joutel F., Levrard B. and Robutel P. 2004. Long-Term evolution and chaotic diffusion of the insolation quantities of Mars. Icarus, Vol. 170, Issue 2, 343–364.
- Smulsky J. J. and Smulsky Ya. J., 2012. Dynamic Problems of the Planets and Asteroids, and Their Discussion, International Journal of Astronomy and Astrophysics, Vol. 2, No. 3, 129-155. http://dx.doi.org/10.4236/ ijaa.2012.23018.
- Newcomb, S., 1895. The Elements of the Four Inner Planets and the Fundamental Constants of Astronomy. Government Printing Office, Washington. 202 p.
- Simon, J.L., Bretagnon, P., Chapront, J. Chapront-Touze, G., Francou, Laskar, J., 1994. Numerical Expression for Precession Formulae and Mean Elements for the Moon and the Planets. Astron. Astrophys., 282, 663-683.
- Smulsky J. J., 2016. Fundamental principles and results of a new astronomic theory of climate change. Advances in Astrophysics., Vol. 1, No. 1, 1–21. http://www. Isaacpub. org/Journal/AdAp.
- Smulsky J.J., 2016. New results on the Earth insolation and their correlation with the Late Pleistocene paleoclimate of West Siberia. Russian Geology and Geophysics, No. 57, pp. 1099-1110.http://dx.doi.org/10.1016/j.rgg.2016.06.009.
- Smulsky, J.J., 2016. Evolution of the Earth’s axis and paleoclimate for 200,000 years. Saarbrucken, GR: LAP Lambert Academic Publishing. 228 p. ISBN 978-3-659-95633-1. http://www.ikz.ru/~smulski/Papers/InfEvEAPC02MEn.pdf. (In Russian).
- Smul’skii I.I., 2013. Analyzing the lessons of the development of the orbital theory of the paleoclimate. Her. Russ. Acad. Sci., No. 83 (1), 46–54.
- Laskar J., Robutel P., Joutel F., Gastineau M. Correia A. C. M. and Levrard B., 2004. A Long-term numerical solution for the insolation quantities of the Earth. Astronomy & Astrophysics, Vol. 428, No. 1, 261–285.
- Sharaf Sh. G. and Budnikova N. A. 1969. Secular changes in the orbital elements of Earth and the astronomical theory of climate change. In: Trans. Inst. Theor. Astron., Issue XIV. – L, Nauka, 48 – 109. [In Russian].
- Smulsky J. J., 2016. The Phenomena of the Sun in the Historical Perspective. Institute of Earth's Cryosphere SB RAS. Tyumen. Deposited in VINITI RAS 11.01.2016, 9-V2016, p. 66 [In Russian]. http://www.ikz.ru/~ smulski/Papers/SunPhnmen.pdf.
- Smulsky J.J., 2018. The Sun’s Movement in the Sky Now and in the Past. Open Access Library Journal, 5, e4250. http://dx.doi.org/10. 4236/ oalib.1104250.