Overview of the Spectral Coherence between Planetary Resonances and Solar and Climate Oscillations

: The complex dynamics of solar activity appear to be characterized by a number of oscillations ranging from monthly to multimillennial timescales, the most well-known of which being the 11-year Schwabe sunspot cycle. Solar oscillations are important because they also characterize the oscillations observed in Earth’s climate and can thus be used to explain and forecast climate changes. Thus, it is important to investigate the physical origin of solar oscillations. There appear to be two possibilities: either the oscillations in solar activity are exclusively controlled by internal solar dynamo mechanisms, or the solar dynamo is partially synchronized to planetary frequencies by planetary forcings. The latter concept has recently gained support from a growing amount of evidence. In this work, we provide an overview of the many empirical facts that would support a planetary hypothesis of the variability of solar activity and emphasize their importance for climate research. We show that the frequencies produced by the complex interactions of all of the planets are coherent with the major solar activity and climate cycles, from monthly to multimillennial timescales, including the well-known Schwabe 11-year solar cycle. We provide some persuasive theoretical and empirical support for the planetary hypothesis of solar and climate variability.


Introduction
Throughout its 4.7-billion-year history, the Earth has seen significant climate changes. For example, throughout the Phanerozoic period, that is, the last 600,000 years, the global surface temperature oscillated between ice ages, with temperatures up to 5°C below 20th-century temperatures, and hot eras, with temperatures up to 14°C above current temperatures. Natural climate changes have occurred at all timescales. So far, comprehending the natural forcings that drive climate changes remains a contentious issue due to a lack of physical understanding of the various involved phenomena. Yet, most of this variability appears to be tied to a wide number of astronomical and solar factors, many of which exhibit various sorts of cyclical behavior. The lack of accurate representation of such forcings in climate models, particularly those responsible for the decadal to millennial oscillations, has a significant impact on climate science and projections for the future decades and centuries and, therefore, this issue must be thoroughly examined. Figure 1 shows a composite of temperature records from the Phanerozoic up to 2015 AD with a schematic spectrum of its climate variance at several timescales [1]. The figure divides the main known periodicities into three main timescales: (a) the galactic periodicities, such as those at about 150 [2,3] and 32 millions years [4], which appear to be due to the solar system's trajectory around the Galaxy, where it moves through its spiral arms and its vertical oscillation with respect to the galactic plane, respectively; (b) the Milankovitch cycles associated with the astronomical orbital forcings that are related to variations in the eccentricity of the Earth's orbit (100,000 and 400,000 years), of the Earth's axial tilt (41,000 years), and of the axial precession (21,000-26,000 years) [5]; and (c) the cycles ranging from decadal to multimillennial scales, which are likely related to solar activity variability cycles such as the 2100-2500-year Bray-Hallstatt cycle, the 1000-year Eddy cycle, the 170-240-year Jose and Suess-de Vries cycles, the secular Gleissberg cycles, and others [6,7]. The interannual timescale oscillations (e.g., the El Niño-Southern Oscillations) are typically associate with the so-called internal variability of the climate system, but they too may have a solar or astronomical origin because there are very similar interannual solar and orbital oscillations as well [8][9][10][11]. Moreover, climate data also present long decadal and multidecadal soli-lunar tidal cycles [12]. Climate change is once again associated to the intensity of cosmic ray fluxes on billion-year timeframes [13].  [14][15][16][17][18][19][20][21]). (B) Schematic spectrum of the climate variance at multiple timescales (adapted from Vinós [1]).
It is worth noting that the link between atmospheric CO 2 content and global temperature has been fairly poor throughout the last 600 million years [2,22,23]. CO 2 concentration, for example, has also lagged for centuries behind temperatures during deglaciation and glaciation periods, as occurred during the last 420,000 years, as shown in the Vostok ice core record [24]; although, if the data are processed in some way and some specific places are analyzed, the two variables appears much more tightly coupled [20,25]. Hence, carbon dioxide cannot have typically been the primary driver of climate changes for nearly all of Earth's history, but rather, it worked as one of the climate (positive) feedback mechanisms in response to solar, astronomical, orbital, and other natural forcings, although it was likely less important than water vapor and clouds. In fact, the atmospheric CO 2 concentration significantly depends on the surface temperature of the oceans and on the status of the biosphere, although it might also be suddenly altered by volcanic activity.
This review complements another recent review from the same authors [26] by focusing on the impact on climate change produced by solar activity variability cycles. These cycles account for most of the decadal to multimillennial climatic variability across the Holocene. We show that these variations appear to be linked to orbital planetary oscillations, thus making the solar activity and climate changes theoretically predictable.

Overview of the Planetary Hypothesis of the Origin of Solar Activity Cycles
The hypothesis of the planetary origin of the 11-year sunspot cycle was first put forward when Wolf [27] suggested that "the fluctuations of spot-frequency rely on the influences of Venus, Earth, Jupiter, and Saturn". A small number of scientists supported the planetary hypothesis of solar activity oscillations during the 20th century [28][29][30], but, during the last 15 years, this interest has grown significantly, leading to an increasing number of publications [31]. Recently, Scafetta and Bianchini [26] offered an extended review that was largely based on Scafetta's prior work. The planetary origin of solar activity oscillations is not restricted to the 11-year solar cycle but includes all major solar oscillations on timescales ranging from monthly to multimillennial. Although the physical mechanisms explaining how the Sun may be modulating climate change are still debated (e.g., there are proposals of additional solar-related forcings [2,3] as well as of complex dynamical responses of the climate system to insolation variability [32,33]), the fact that the Earth's temperature and the solar records display common oscillations throughout the Holocene has substantial implications for climate science [7,[34][35][36][37][38][39][40][41][42][43][44][45]. Indeed, numerous criteria show that solar variability is driven by planetary beats that also affect a variety of terrestrial variables, including 14 C and 10 Be production, Earth's rotation, ocean circulation, paleoclimate, geomagnetism, and so on [46].
At the moment, there are just two main hypotheses regarding the origin of solar oscillations. One idea is that solar activity changes are controlled solely by the solar dynamo mechanisms (let us call it the "reductionist" hypothesis). The second possibility is that the solar dynamo itself is partially synchronized by external harmonic planetary forcings, such as tides or other possible mechanisms (let us call it the "holistic" hypothesis). The debate is still ongoing since the physics of solar activity changes is still poorly understood. However, we think that the arguments for a planetary synchronized solar dynamo concept should be better investigated because they could be the key to understand both solar activity and climate oscillations.
In recent years, there has been increasing interest in the concept that planets can influence solar activity [31,47]. We believe that this approach could yield the development of methodologies useful for long-time-range solar activity forecasts, similarly to how ocean tides are currently predicted using harmonic constituent models. This research should also have an impact on space weather and climate change forecasting [12,37,48]. We also address several frequent critiques of the planetary hypothesis of solar and climatic variability.
One of the primary objectives of this research is to find repeating patterns and harmonics in solar activity records that could be linked to planetary harmonics. The frequency and phase of planetary harmonics is compared with solar activity records. However, the amplitude of the secular solar activity changes cannot be determined by the planetary harmonics. Actually, the scientific literature is still hotly debating the precise amplitude of the secular solar activity changes.
Numerous total solar irradiance (TSI) proxy models have been proposed [38,42,[49][50][51][52][53][54][55][56], which greatly differ from each other regarding the amplitude of the solar activity secular changes. For example, the TSI rise from the Maunder Minimum to the minimum of solar cycle 24 (2008-2019) varies from 0.75 W/m 2 to 6.3 W/m 2 according to the proposed TSI proxy model. Some models show a very limited range of secular variability [38,[50][51][52]56], while other models show a much larger secular variability [42,49,54,55]. Yeo et al. [57] claimed that the dimmest state of the Sun cannot be below 2 W/m 2 . However, Schmutz [58] objected that a too low TSI secular variability would poorly explain the good correlations found between paleoclimatic and solar records during the Holocene.
The fact that all TSI proxy models indicate simultaneous, albeit not identical, multidecadal periods of low and high solar activity is a key aspect. Secular/bisecular large oscillations, for example, have been observed. The grand minima of these solar activity cycles throughout the last millennium are named Oort, Wolf, Spörer, Maunder, and Dalton. There are other millennial oscillations known as Eddy cycles, as well as multimillennial oscillations, such as the Bray-Hallstatt cycles, which last roughly 2100-2500 years [37,38,53,59,60]. Similar oscillations are also found in climate records [38,39,41,42,[60][61][62]. The origin of solar activity oscillations is unknown, but the phenomenon could be due to complex patterns generated by the interference of several oscillations, which could be possibly linked to planetary harmonics. If this is so, a planetary model of solar activity changes could be empirically adopted to hindcast and forecast them.
There are also shorter multidecadal oscillations, such as quasi-45-, 60-, and 80-100-year solar activity oscillations [34,35]. Some of them might already be confirmed using the TSI satellite data, available since 1978. However, even here there are uncertainties.
If solar activity changes are partially controlled by planetary harmonics, for example, through tides, a planetary model could help in solving the TSI satellite composite controversy between ACRIM and PMOD. Solving this controversy is rather important for both solar and climate change studies. In fact, the current global climate models (GCMs) adopt a solar forcing derived from TSI proxy models with low secular variability, which match the PMOD TSI satellite composite decreasing trend from 1980-2020 [71]. If, on the other hand, TSI grew from 1980 to 2000, as suggested by the ACRIM TSI satellite composite, the solar contribution to climate change should properly evaluated [72,73].

Empirical Evidences for the Planetary Origin of Solar and Climate Variability Cycles
Several academics have recently investigated the idea that the planets of the solar system might have an impact on solar activity since, especially over the last 15 years, more and more scientific data have been discovered to support this hypothesis. Many of these findings were examined by Scafetta and Bianchini [26]. Let us now overview the most supporting evidences. Section 3.1 demonstrates that realistic tidal models predict a tidal cycle of 10-12 years, which corresponds to the observed Schwabe 11-year sunspot cycle. Section 3.2 provides additional data in support of a planetary synchronized solar dynamo hypothesis.

The Venus-Earth-Jupiter-Saturn Model for the Schwabe 11-Year Solar-Activity Cycle
Since the second half of the 19th century, it has been clear that complex combinations of several planets should have been examined in order to find a potential connection between solar activity and planetary motion. Admittedly, all the planets should influence the Sun concurrently. However, according to Wolf [27], the quartet Venus, Earth, Jupiter, and Saturn had to be the basis of any minimal model that could account for the 11-year solar cycle. In fact, single-planet models are ineffective. For example, due to its eccentric orbit, Jupiter would produce a tidal oscillation of 11.86 years, but this period is too long to suit the 11-year sunspot cycle, albeit fortuitous correlations may appear during certain time intervals.
Wolf's hypothesis was based on the observation that the planetary order of the relative tidal significance is Jupiter, Venus, Earth, Mercury, Saturn, Mars, Uranus, and Neptune [74]. Mercury, on the other hand, appeared to move too quickly, whereas Mars, Uranus, and Neptune seemed to be less important because their tides are much weaker than those of the other planets. Wolf's suggestion of the planets Venus, Earth, Jupiter, and Saturn as a potentially working quartet made sense. Some authors even focused on an even more constrained planetary collection that only consists of Venus, Earth, and Jupiter, which are the first three planets with the largest tidal influence, but, as we argue below, this triplet appears to be too reductive.
In any case, rudimentary versions of the Venus-Earth-Jupiter triple syzygies tidal alignment model had already been proposed over a century ago [28,29,75]. More recently, various complex variants of the same model were developed [74,[76][77][78][79][80][81][82][83]. The general conjecture on which these models are based is that the smaller the difference in days between the conjunctions/oppositions of Jupiter and Venus with Earth and the Sun, the greater the action they exerted on the Sun. The planetary functions created by these models appear to be tightly associated with the 11-year sunspot cycle from 1700 to the present; see Figure 3.
Scafetta [74] demonstrated that the principal recurrent period of the Venus, Earth, and Jupiter triple-syzygies tidal alignment model is given by the invariant inequality where P V = 224.701 days, P E = 365.256 days, and P J = 4332.589 days are the sidereal orbital periods of Venus, Earth, and Jupiter, respectively. Equation (1) can be also rewritten as which has significant mathematical properties, because P VEJ can be interpreted as the complex beat function between the third harmonic of the Venus-Earth synodic cycle (2P VE = (1/P V − 1/P E ) −1 = 1.59867 year, P VE = 0.79934 year) and the second harmonic of the Earth-Jupiter synodic cycle (2P EJ = 1/P E − 1/P J −1 = 1.09207 year). The factor 1/2 in Equation (1) transforms synodic cycles into tidal cycles. Therefore, Equation (2) can also be interpreted as the spring beat period between the Venus-Earth synodic period (2P VE ) and the beat between the spring tide between Venus and Earth (P VE = 0.799 year) and the spring tide between Earth and Jupiter (P EJ = 0.546 year). Table 1 reports the spring periods of all planet pairs. Equation (1) is also an important example of orbital invariant inequality, which is the orbital resonance characterized by the property of being invariant with respect to any rotating system. This property should be physically necessary for synchronizing the solar dynamo, which is generated by the differential rotation of the Sun [7]. The tidal beats of Equation (1) can be simulated by the function where t VE = 2002.8327 and t EJ = 2003.0887 are epochs of one Venus-Earth conjunction and one Earth-Jupiter conjunction, respectively. The cosine function is used to predict the spring tidal maxima during the conjunction epochs for each spring tide. As a result, Equation (3) is entirely based on astronomical data. Figure 3A compares Equation (3) with the annual sunspots number from 1700 to 2021. The graphic demonstrates that there is good timing between the sunspot cycle and the 11.07-year Venus-Earth-Jupiter alignment beats. In essence, it seems that more sunspots form when the three planets are more often closely aligned, which is when their tides can more forcefully twist the tachocline and affect the solar dynamo. Figure 3B,C demonstrates that different approaches of the Venus, Earth, and Jupiter triple-syzygies model are as effective at predicting both the timing and the frequency of the sunspot cycles across centuries [76,82].
We notice that without the 1/2 factor in Equation (1), the major period of the Venus, Earth, and Jupiter triple-syzygies alignments would be 22.14 years, which corresponds to the 22-year Hale solar magnetic cycle [83].
A realistic, although simplistic, tidal function can be written as where d P is the planet's average distance from the Sun, m P is its mass, T S = 27 days is the solar rotation period, and T P is the orbital period of the planet P. Moreover, d Pa and t Pa are the aphelion distance and one of its occurrence epochs, respectively, whereas α PJ, 2000 is the angular distance of the planet P from Jupiter on 1 January 2000 at 00:00. The planetary data are reported in Table 2.  [74] pointed out that what physically matters is the power dissipated by the tides inside the Sun, i.e., their work done in a time unit, which is proportional to the following function Equation (5) is proportional to the tidal power dissipated inside the Sun. Scafetta [74] hypothesized that the fraction of the dissipated power in the solar core would trigger an increase and a modulation of its luminosity production. Actually, Equation (4) approximates the full luminosity-tidal equation proposed by Scafetta [74], which was given as which is similar to Equations (4) and (5) but uses the actual ephemeris coordinates of all eight planets from Mercury to Neptune. Moreover, the tidal function is integrated from the solar core to the surface, and an amplification factor is introduced to take into account the H-burning rate increase due to the tidal gravitational power dissipated in the solar core, as suggested by Scafetta [74]. The original paper explains the various parameters and functions of Equation (6), which are R SP (t), the distance of a planet from the Sun; m P , the mass of the planet P; α P,t , the angles denoting the location of the planet P relative to the angular position φ on the Sun; Q −1 , the effective tidal dissipation factor; R S , the radius of the Sun; ∆t, the integration time interval, which is set to 1 day; K(χ), the function for converting gravitational power into TSI at 1 AU from the Sun; ρ(χ), the solar density function; χ = r/R S , the normalized distance from the Sun's center; and G, the universal gravitational constant. Figure 4A depicts the function f (t) (Equation (5)) with only Jupiter and Saturn as inputs. A quasi-11.86-year oscillation is observed, which is tied to Jupiter's orbital period, because its distance from the Sun varies between the aphelion and the perihelion. Yet, the P JS = 9.93 year neap-spring tidal cycle between Jupiter and Saturn also influences such oscillations. As a result, at the decadal timeframe, the tidal signal is clearly defined by two cycles with durations of 9.93 and 11.86 years, with an average length of close to 11 years. The figure also shows a clear beat with a period of about 60 years, whose maxima occurred in 1940-1950 and 2000-2010, which correspond to the maxima of a quasi-60-year oscillation observed in global surface temperature and other climatic records [37,85]. Figure 4B uses Wolf [27]'s planetary quartet to illustrate the function f (t) using Venus, Earth, Jupiter, and Saturn. The synthetic data clearly show the same 10-12-year oscillation. Once again, Figure 4C demonstrates that the patterns of quick tidal changes recur every 11 years. Figure 4D shows the hypothetical brightness signal produced by all planetary tides estimated by Scafetta [74] using Equation (6).   Figure 4B using Venus, Earth, Jupiter, and Saturn and of the tidally induced luminosity function (Equation (6)), using all the planets derived by Scafetta [74] and depicted in Figure 4D.  Figure 4B. (Bottom) The power spectrum of the tidal-luminosity function calculated in Scafetta [74] and shown in Figure 4D.
The two power spectra share a number of harmonics, particularly the planets' neapspring tides. The power spectrum clearly increases between 9 and 13 years because of the clustering of the 9.93-year neap-spring tidal cycle between Jupiter and Saturn and the 11.86-year tidal oscillation associated with the elliptic orbit of Jupiter. The neap-spring Jupiter-Saturn tidal oscillation is approximately three times smaller than Jupiter's orbital tidal cycle, which allows the two tidal cycles to produce complex beats. The measured Schwabe cycle spectral band fits the two tidal cycles well. As a result, using the real orbits of Jupiter and Saturn is critical for appropriately connecting the planetary tidal power spectrum to the 11-year sunspot cycle. The result clearly refutes the assertion by Okal and Anderson [86] and Nataf [87] that there is no proof for a periodicity close to 11 years in the tidal functions [84].
Scafetta [37,74] showed that the 11-year sunspot cycle is actually composed of three closely spaced frequencies. In fact, the presence of the two tidal frequencies with periods at about 10 and 12 years implies that the Schwabe cycle could be characterized by a primary spectral peak close to 11 years (which could be generated by the primary synchronization of the solar dynamo), surrounded by two minor spectral peaks at about 9.93 and 11.86 years, which correspond to P SJ and P J . This prediction is confirmed in Figure 6A, which shows several power spectra of different sunspot number records and suggests that the solar dynamo could be actually controlled by the neap-spring tide between Jupiter and Saturn and the tidal oscillation associated with Jupiter's eccentric orbit, because the two frequencies appear at the left and the right side of the main 11-year spectral peak. The same result was obtained also by other authors, such as by Tan and Cheng [88]. Scafetta [74] corroborated this conclusion by also showing that the distribution of the solar cycle lengths is bimodal instead of being a simple Gaussian centered at 11 years. In fact, of the 23 solar cycles that were examined between February 1755 and December 2008, 11 of them had a period between 9.0 and 10.5 years, with an average of 10 ± 0.51 years (which optimally corresponds to the period of the Jupiter-Saturn neap-spring tide), and 12 of them had a period between 11.25 and 13.6 years, with an average of 12 ± 0.75 years. No sunspot cycle length was recorded with a period larger than 10.5 years or smaller than 11.25 years. To date, solar cycle #24 from December 2008 to December 2019 is the only "Schwabe 11-year solar cycle" to have lasted exactly 11 years. This cycle should be included in the second list of the longer cycles, because the mean between P JS and P J is about 10.9 years. Thus, including this last cycle, the average of the 13 solar cycle lengths larger than 10.9 years is 11.9 ± 0.7 years, which agrees well with the tidal period associated with the orbit of Jupiter. This result supports the hypothesis that two physical attractors operating at the periodicities of the Jupiter-Saturn neap-spring tide (P SJ = 9.93 years) and of the Jupiter tide (P J = 11.86 years) could be responsible for the apparent dynamics of the sunspot cycle. Figure 6B shows the likely bimodality of the distribution of the solar cycle lengths, because two peaks appear on the left and right side of the 10.9-year period. The depicted distribution function was obtained using a smooth k-density function with a bandwidth of σ k = 0.4 years, estimated by averaging the default bandwidth (σ k = σ(4/3N) 0.2 , with σ being the standard deviation of N samples) for the two sets under the hypothesis that the distribution is bimodal. Table 3 reports the main characteristics of the 24 solar cycles from 1755 to 2019. Scafetta [37] provided even additional confirmation of this finding by showing that a three-frequency model based on the two tidal harmonics produced by Jupiter and Saturn hindcast a number of traits found in long solar data throughout the Holocene; see Section 3.2.
Using Wolf [27]'s planetary set, the five strongest tidal spectral peaks are P V J , P EJ , P VE , P SJ , and P J , as shown in Figure 5. The periodicities P SJ and P J fit the Schwabe 11-year sunspot cycle perfectly. Let us now discuss the key physical property associated with the Venus-Earth-Jupiter triple-syzygies tidal alignment paradigm, which is based on the other three spring tidal harmonics: P V J , P EJ , and P VE .
The high-frequency, which is connected with the spring tides between Venus and Jupiter (P V J = 0.3244 year), Venus and Earth (P VE = 0.7993 year), and Earth and Jupiter (P EJ = 0.5460 year), shows periodic patterns. The main patterns can be identified by looking for combinations of integers η 1 ,η 2 , and η 3 , such that and such that recurrence times are as close as possible to each other The three best combinations (η 1 , η 2 , η 3 ) are (32, 19, 13) = 10.38 ± 0.01 years, (34,20,14) = 11.05 ± 0.1 years, and (35,21,14) = 11.34 ± 0.1 years. The best centered combination between P JS and P J is (34,20,14). By averaging the latter combination of the three fast spring tides among Venus, Earth, and Jupiter, and the two tidal harmonics of Jupiter and Saturn, we obtain We also found that the recurrent pattern (η 1 , η 2 , η 3 ) = (32, 19, 13) = 10.38 ± 0.01 years optimally fulfills the more general condition η 1 · P V J ≈ η 2 · P EJ ≈ η 3 · P VE , when the integers are varied between 1 and 50. Interestingly, a quasi 10.4-year cycle is found in quite a number of climatic records [12,85] and even in meteorite fall records [90]. The beat between the 10.4-year recurrent tidal period and the spring tide between Jupiter and Saturn (P SJ = 9.93 years) produces a fast beat periodicity of P VEJS f ast = 2/(1/10.38 + 1/P JS ) = 10.15 years. If this cycle beats with the Jupiter orbital period, a new periodicity is produce at P central = 2/(1/10.15 + 1/P J ) = 10.95 year, which is very close to the central spectral peak of the the sunspot number record shown in Figure 6.
In conclusion, the five most important tides produced by Venus, Earth, Jupiter, and Saturn (P V J , P EJ , P VE , P JS , and P J ) fit the whole spectral band of the Schwabe 11-year sunspot cycle. The above considerations show that the Schwabe 11-year solar cycle is caused by the combined influence of Venus, Earth, Jupiter, and Saturn on the Sun, as initially proposed by Wolf [27] in the 19th century.

Empirical Evidences for Planetary Control of Solar Variability across Several Timescales
The spectral compatibility of the tidal planetary models with an 11-year solar cycle is just the first piece of evidence supporting the planetary synchronized solar dynamo theory. There is a substantial body of data showing that the primary planetary frequencies do indeed characterize changes in solar activity at all timescales, from the monthly to the multimillennial ones.
The main planetary harmonics that could have a possible physical relevance are the orbital periods, spring tidal periods, synodic conjunction periods, and their combinations, which were labeled orbital invariant inequalities [7]. Here, we give a brief summary of what we consider to be the most important empirical findings supporting the planetary theory of solar activity changes.

Monthly to Annual Timescales
Bigg [91] discovered a minor but constant periodicity in the relative daily sunspot number records at the sidereal period of the planet Mercury, modified by Venus, Earth, and Jupiter. This author came to the conclusion that extrasolar factors could have an impact on the creation of sunspots.
The planetary tidal signal, as seen in Figure 5, is made of a number of quick oscillations that are mostly related to the spring tidal periods between pairs of planets (see Table 1). The most important ones are P V J , P EJ , and P VE , which, as we saw in Section 2, could be connected to the Schwabe 11-year solar cycle. Therefore, the question is whether the solar records show any evidence of these rapid tidal oscillations and/or of their corresponding synodic planetary oscillations, which have periods twice as long. Scafetta and Willson [92,93] presented a study in which numerous daily total solar irradiance (TSI) satellite records were compared with the spring and synodic periods of the planets from the monthly to yearly timescales. Figure 7 summarizes the key findings. Figure 7A shows the ACRIM TSI satellite composite for the years 1992.5 through 2013.5, which collects the data from ACRIM II and ACRIM III high-precision measurements. Here, the TSI record is preferred to other solar observables (such as the number of sunspots), because it is the most accurate and complete indicator of changes in solar activity. The ACRIM TSI composite was chosen above other TSI datasets due to its duration, for the perception of being more sensitive to minute TSI changes, and for being more accurate than other proposed satellite records [65,93]. Similar results were also obtained by using the PMOD/VIRGO TSI composite and SORCE/TIM TSI satellite record [92,93].  Table 1 [92,93]. The period labeled " * " appears related to Mercury and the solar rotation. Figure 7B shows the periodogram of this TSI record. The colored vertical boxes in Figure 7B illustrate the predicted primary planetary frequencies; the black labels indicate the neap-spring tidal periods, while the blue labels indicate the planetary synodic periods. See Table 1 and References [92,93].
The figure shows that the TSI periodogram includes all of the primary spring and synodic planetary harmonics from monthly to the annual timescales. The prominent spectral peaks between 0.30 and 0.33 years refer to Venus's spring tidal cycles with Jupiter, Saturn, Uranus, and Neptune; see Table 1. Other TSI spectral peaks are seen in the frequency range of 0.065-0.105 years, which corresponds to the differential solar rotation: close to Mercury and Venus's spring tidal periods (P MV = 0.198 years), Earth and Jupiter (P. E.J. = 0546 years), and Venus and Earth (P VE = 0.799 years), as well as near the synodic cycles between Mercury and Venus (2P MV = 0.396 years) and Earth and Jupiter (2P EJ = 1.092 years). A 0.5-year oscillation is also observed; this harmonic is most likely caused by the Earth's orbit crossing the solar equatorial plane twice a year, as the solar luminosity depends on the Sun's latitude. In fact, sunspots almost never appear below 5 • or above 40 • , both in north and south latitudes. Figure 7B also shows a significant TSI spectral peak with a period P ≈ 20 days. This period may be related to Mercury and to the solar rotation. In fact, assuming an average solar rotation of T S = 27.5 days (at about 45 • latitude, [94]), its synodic cycle with Mercury (P M = 88 days) would be 1/(1/T S − 1/P M ) ≈ 40 days. Thus, the relative tide would have a period of about 20 days. The 1.09-year oscillation is outstanding in the TSI record, particularly from 1998 to 2004, during the maximum of solar cycle 23. More specifically, it was discovered that TSI rises during the conjunction of Earth and Jupiter, that is, when the two planets are at their closest angular separation relative to the Sun ( Figure 7A, top), as if the solar side facing Jupiter were slightly brighter so that the Earth can receive more light when it crosses the line connecting the Sun and Jupiter [92]. In fact, it has been discovered that stars with hot Jupiter-like planets exhibit a comparable hot spot brightening, albeit more intense. These brightenings are most likely caused by a direct magnetic contact between the planet and the stellar surface [95,96]. The finding is further supported by the existence of a near 13-month periodicity in the flux of cosmic rays, probably related to the Earth-Jupiter synodic cycle [97]. This result underlines the importance of looking into how the solar system's space weather and interplanetary fluxes of particles are structured in relation to the relative positions of the planets.

Multidecadal to Millennial Timescales
A large number of scientific data indicate that the Sun's activity is also characterized by a variety of unique oscillations that occur across timeframes ranging from decades to millennia. These periods seem to fall into distinct ranges. The most frequently documented solar oscillations cover the periods of 40-45 years, 55-65 years, 80-105 years (Gleissberg cycle), 115-150 years, 170-240 years (Jose and Suess-de Vries cycles), 800-1200 years (Eddy cycle), and 2000-2500 years (Bray-Hallstatt cycle) [35,59]. These oscillations fall into the major clusters created by the body of the planetary periods, as described by Scafetta and Bianchini [26] and their cited references. Let us now review the key findings.
(1) Scafetta [37] developed a multiscale harmonic solar and climatic model based on the Jupiter-Saturn neap-spring tidal oscillation, Jupiter's orbital tidal oscillation, and a solar dynamo cycle that was hypothesized to have a period of 10.87 years (see Figure 6). Four primary beats with periods of P S13 = 60.95 years, P S12 = 114.78 years, P S23 = 129.95 years, and P S123 = 983 years result from the combination of the three cycles. The model successfully predicts all the major secular solar activity maxima and minima over the last millennium, including the Oort, Wolf, Spörer, Maunder, and Dalton grand solar minima ( Figure 8A), which are approximately given by the quasi-115-130-year oscillation predicted by the model and are observed in all solar activity proxy records, including the 14 C and 10 Be radionucleotide [38,53]. The same model also predicts a quasi-1000-year oscillation, which has been well observed in detailed temperature reconstructions of the Northern Hemisphere over the last 2000 years ( Figure 8B). Finally, Figure 8C compares the same solar harmonic model against a global surface temperature record from 1850 to 2010 [37], where a similar 60-year modulation plus a secular upward trend is observed. Figure 8C, which was published in 2012 by Scafetta [37], predicted an incoming grand solar minimum from 2015 to 2045, induced by the interference of the 115-130-year oscillation and the 60-year oscillation produced by the planetary model following the Maunder, the Dalton, the 1900-1920, and the 1970s grand solar minima. An alternative planetary model based on the Venus, Earth, and Jupiter tides did the same [98,99]. Courtillot et al. [100], still using planetary resonances, predicted the same modern solar minimum. Scafetta [37] also predicted a new solar maximum between 2050 and 2070, which would also coincide with the millennial solar cycle maximum. Indeed, Scafetta [37] predicted the incoming grand solar minimum from 2015 to 2045, also from the similitude of solar cycles 1-6 with the solar cycles 20-24. Scafetta [7] and Scafetta and Bianchini [26] observed that planetary resonances are not distributed randomly but rather cluster around specific frequency ranges. This could explain why different models based on planetary resonances can produce comparable findings.  10 Be and 14 C cosmogenic isotopes [52,53]. (B) The identical model (red) overlaid over a reconstruction of the Northern Hemisphere's proxy temperature [101]. (C) The same planetary-solar model (red) is displayed against the HadCRUT4 global surface temperature record (black) [102], which was combined in 1850-1900 with Moberg et al. [103] proxy model (blue). The green curves illustrate the skewness of the quasi-millennial oscillation, which approximates the millennial temperature oscillation from 1700 to 2013 [37,48]. The black dots represent the maxima of a quasi-60-year oscillation produced by the model and also observed in the climate records.
(2) During the course of the last 10,000 years (the Holocene), the quasi-millennial three-frequency beat cycle properly predicts the observed solar and climate millennial cycles [37,48]. Figure 9A shows a comparison between the TSI proxy model by Steinhilber et al. [52] (red) and the millennial oscillation indicated by the three-frequency planetarysolar model (blue) [37,48]. Figure 9B shows a quasi-millennial oscillation in the summer temperatures of the European Alps throughout the Holocene [104]. The millennial oscillation predicted by the planetary model appears to be well synchronized for 12,000 years with the temperature record. An exception is observed in the 5500-5000 BC time interval, which, however, could be due to a poor coverage of data, since the two peaks are clearly seen in the GISP2 ice core air temperature record of the Greenland Ice Sheet [105], as shown Figure 9C. Another discrepancy is seen in 1500-1000 BC in the phase shift between the GISP2 temperature peak and that of the model, which, however, is better reproduced by the alpine record. These discrepancies could be explained by chaotic climatic dynamics and by possible dating and measurement uncertainties. (3) Scafetta [48] demonstrated that the same planetary-solar three-frequency model also generates two sets of oscillations with periods of 57, 61, and 65 years and 103, 115, 130, and 150 years. The latter four periods, which control grand solar minima and maxima, were also identified in records of solar activity from 340,000 to 320,000 years ago. Thus, Scafetta [48] demonstrated that Cauquoin et al. [106] were wrong when they asserted that there was no evidence for a planetary modulation of solar activity around 330,000 years ago. The finding refuted the main claim by Smythe and Eddy [107] that planetary tidal models were unable to account for phenomena, such as the Maunder sunspot grand minimum, and also refutes critiques such as that all solar records from the secular to millennial range, with the exception of the Schwabe cycle, could be some kind of random fluctuation [108,109].
(4) Scafetta et al. [6] and Scafetta [7] specified that variations in solar activity on timescales ranging from multidecades to millennia are spectrally coherent with the orbital invariant inequalities generated by Jupiter, Saturn, Uranus, and Neptune. The orbital invariant inequalities are the beats among the planetary synodic cycles. In the differentially rotating Sun, these periodicities were shown to have mathematical features that could favor the setting of synchronized dynamics.  [6,59]. The latter is accurately predicted by the (1, −3, 1, 1) invariant inequality produced by the beats between the synodic cycles of Saturn with Jupiter and Uranus and Neptune-the period being 2318 years. Figure 10 depicts the outcome. Figure 10A compares a solar (∆ 14 C) and a climate (δ 18 O) record from 9500 and 6000 years ago (adapted from (Figure 5 in [34])). Figure 10B compares the orbital invariant inequalities of the solar system (red bars) determined by Scafetta [7], with the cross-spectral analysis peaks between the two records shown in Figure 10A. Finally, Figure 10C directly compares the (1, −3, 1, 1) invariant inequality oscillation (blue), whose phase is calculated by astronomical data only [7,110], with the ∆ 14 C (‰) record ( [111], IntCal04.14c, black and red). In both cases, a perfect match between theoretical and experimental oscillations exists. These same cycles also emerge in the solar inertial motion [110,[112][113][114]. The Bray-Hallstatt cycle of 2000-2500 years is observed in several solar [59,115,116] and climate records [6,117,118] as well. For example, the Greek Dark Ages (between 1000-550 B.C.) and the Little Ice Age (between about 1350-1850 A.D.) are well-established and historically confirmed colder periods of the Holocene, characterized by numerous historical climatic crises which occurred during the solar activity lows related to the Bray-Hallstatt cycle [116], as Figure 10C shows.

Miscellaneous Evidences
There are additional empirical evidences that support the planetary hypothesis of solar activity variations.
For example, Hung [76] discovered that at least 25 of the 38 biggest solar flares began when one or more of the planets-Mercury, Venus, Earth, or Jupiter-were either immediately above or just on the opposite side of the Sun. Bertolucci et al. [119] discovered a statistically significant increase in solar flares and EUV emission from 1976 to 2015 and 1999 to 2015, when one or more planets had heliocentric longitudes mostly between 230 • and 300 • . They also found a link between total atmospheric electron content (TEC) and the orbital positions of the inner three planets from 1995 to 2012. Petrakou [120] discovered connections between solar flares and planetary harmonics as well.
Planetary harmonics were also discovered in aurora [121,122] and climate records [12,48]. Sharp [123] showed that several angular momentum perturbations on the Sun-mostly induced by Uranus and Neptune-coincide with the maxima observed in 14 C and 10 Be solar proxy records during the Holocene. Similarly, McCracken et al. [115] suggested that solar activity is modulated by planetary harmonics, because the incoming cosmic ray flux intensity was found to be low when Uranus and Neptune were in superior conjunction and higher when they were in inferior conjunction. Moreover, excellent agreement was found by Abreu et al. [124] between the long-term cycles in solar activity proxies and the periodicities in the planetary torque acting on a hypothesized nonspherical tachocline at the base of the Sun's convective zone. In particular, they found that some periodicities remained phase-locked for more than 9400 years. An alternative model of the Venus-Earth-Jupiter tidal torque was proposed by Wilson [83]. The Chandler wobble also exhibits harmonics of the Jovian planets [125].
In this overview, it is not possible to examine all of the empirical results suggesting that planetary motion affects solar activity. We acknowledge, however, that not all published results are equally credible: some may be dubious, and others may appear to be nothing more than coincidences and numerology. However, it should be acknowledged that there already is a substantial body of evidence in the scientific literature that strongly supports the planetary synchronized solar dynamo hypothesis. This mass of evidence must be seriously taken into consideration.

Evidences for Interannual and Multidecadal Planetary Periods in Global Surface Temperature Records
In the previous subsections, we have seen that planetary cycles characterize the variability observed in both solar and climate records (Figures 8-10). Additional evidence for a signature of multiple planetary cycles from the interannual to the multidecadal scales in the global surface temperature records from 1850 to 2020 has been discussed in several studies [11,12,48,110,126]. Figure 11A,B compare the time-frequency analyses between the HadCRUT3 global surface record [127] and the speed of the Sun relative to the center of mass of the solar system. As first observed by Scafetta [48], and later confirmed by sophisticated spectral coherence investigations [110,126], the observed changes in global surface temperatures mirror multiple astronomical cycles at the decadal and multidecadal scales. The Sun's speed has major spectral peaks at periods of 5. 93, 6.62, 7.42, 9.93, 11.86, 13.8, 20, and 60 years. The majority of them are connected to Jupiter and Saturn's orbits. The temperature record has major spectral peaks at periods of about 5. 93, 6.62, 7.42, 9.1, 10.4, 13.8, 20, and 60 years. Therefore, the majority of the climate oscillations appear to match the astronomical cycles, with the exception of the 9.1-and 10.4-year cycles. Scafetta [12,128] argued that the 9.1-year global surface temperature cycle is related to the lunar apsidal line rotation period of 8.85 years, the Saros eclipse cycle's first harmonic of 9.0 years, and the soli-lunar nodal cycle's first harmonic of 9.3 years [9,30,128,129]. The 10.4-year temperature cycle falls within the sunspot cycle frequency band and should represent the temperature signature of the 11-year solar cycle, because it is consistent with the variability of the solar cycle observed from 1900 (when it was longer) to 2000 (when it was shorter).
The aforementioned periodicities were essential in order to build a semiempirical global surface temperature model based on astronomical oscillations [12,85]. The same model was later upgraded with the addition of some interannual cycles by Scafetta [11]. To validate the model predictions made after 2014, the calibration was only based on the HadCRT4 data [102] from 1850 to that year. It was discovered that the model was accurate in predicting the two strong El Niño events that took place in 2015-2016 and 2020.
The model assumed two components: an astronomically based oscillating component superimposed on an anthropogenic plus volcano component. The astronomical component is based on a number of interannual periodicities, the 9.1-year solar-lunar cycle, and the astronomical-solar cycles at 10.5, 20, and 60, as well as the two longer cycles at 115 years and at about 1000 year. The millennial cycle was assumed to be asymmetric (because of additional multisecular and multimillennial cycles) with maxima in 1080 and 2060 and a minimum in 1700 (the Maunder Minimum). The millennial cycle was deduced from the Jupiter-Saturn model that is shown in Figure 8.
The equation of the harmonic model is which uses 13 harmonics whose frequency, amplitude, and phase of the adopted oscillations including the fast interannual ones, reported in Table 4.  [102,127,130]. (D) The same record compared against the solar astronomical harmonic climate model by Scafetta et al. [131], as updated in Scafetta [11].
The anthropogenic plus volcano component was estimated by dividing the mean of the ensemble of CMIP6 global climate model simulations of global surface temperature (GCMs). Scafetta [85]'s model made use of the CMIP5 GCMs. This halving was suggested by the fact that the adopted natural oscillations already accounted for at least 50% of the warming observed from 1970 to 2000, whereas the models used a set of radiative forcings that hindcasted that the observed warming was primarily anthropogenic during the same period [85,128]. According this estimate, natural variability has contributed to the global warming since 1900 substantially more than what climate models and some research suggest [132][133][134].
In fact, the specific solar forcing adopted by those models used the the low secular variability TSI record by Matthes et al. [71]. However, there exist different alternative reconstructions of the solar forcing that suggest a much larger secular variability [42,54,65], and those based on the ACRIM total solar irradiance satellite composite [42,65] show a significant solar luminosity increase between 1970 and 2000. We notice that an increase in solar activity from 1970 to 2000 is also predicted by the planetary-solar model shown in Figure 8. Table 3 also shows that solar cycle 22 (from September 1986 to August 1996) was very short, 9.92 years, and short cycles are usually linked to higher solar activity [135]. Figure 11C plots the ensemble average simulations of the CMIP6 global circulation models (GCMs) [136], utilizing historical forcings (1850-2014) that are extended from 2015 to 2100, according to three distinct shared socioeconomic pathway (SSP) scenarios compared with the HadCRUT3, HadCRUT4 and HadCRUT5 global surface temperature records [102,127,130]. The adopted emission scenarios are SSP2-4.5 (intermediate GHG emissions), SSP3-7.0 (high GHG emissions), and SSP5-8.5 (very high GHG emissions). Using the same SSP scenarios, Figure 11D compares the three temperature records to the proposed semiempirical astronomical harmonic model. The comparison of panels C and D shows that the semiempirical harmonic model performs better than the global climate models in hindcasting temperature changes from 1850 to 2020. The proposed model also forecasts mild warming for the ensuing decades [11,85].

Addressing the Critiques
The last 15 years have witnessed an increase in the number of publications advocating the planetary hypothesis of solar activity and climate changes. This is a classic example of science progressing from empirical relationship towards hypotheses and predictions and, hopefully, to a full theory in the future. It is clear that we are not there yet and, as a consequence, only a minority of scholars are today endorsing the planetary control hypothesis of solar and climate changes. Many critiques still exist. Some are physically founded and others empirically founded, while still others are blatantly unscientific. Let us go over them quickly.

How should the Sun's Wobbling be Understood?
Several authors found evidence supporting the planetary hypothesis for the origin of the solar activity oscillations by comparing solar activity records with the inertial motion of the Sun around the barycenter of the solar system [6,12,112,114,[137][138][139]. The Sun's wobbling mechanism is sometimes criticized as having no effect on solar activity, because the center of the Sun is in free-fall movement with no significant forces stressing the star. This critique ignores the fact that the Sun's wobbling function can be used as a proxy for conveniently evaluating most of the oscillations of the gravitational field within the solar system due to the orbiting planets. It is the observed frequencies that must be compared with the solar activity oscillations. This approach works because all the dynamical functions of the planetary orbits share a common set of frequencies related to the orbital periods, planetary conjunction cycles, and their various beats. For this reason, the Sun's wobbling can be used for partially testing the planetary hypothesis of solar activity variability from the spectral point of view, because it should share at least a set of frequencies with those of the actual physical mechanisms responsible for the observed oscillations.
The critique also ignores that the Sun is observed wobbling from any point outside the solar system. Similarly, the fluxes of anything coming from far away toward the Sun-such as cosmic rays, interstellar dusts, etc.-would fluctuate with the same frequencies of the solar wobbling. Thus, any forcings associated with these fluxes would also fluctuate in the same way.
Charvátová [112] first noted that solar wobbling presents alternating complicated ordered and disordered dynamics that are correlated with, for example, the Bray-Hallstat solar and climate oscillations. Scafetta et al. [6] noted that such ordered and disordered dynamics, being related to changes in the gravitational and magnetic properties of the heliosphere, could modulate the particle fluxes from outside and within the solar system. These could then affect both the solar activity and/or the Earth's climate.
Suggestions of probable physical mechanisms associated with the wobbling, which might affect the dynamo mechanism, the large scale structure of the heliosphere and the climate system, have been discussed in Scafetta et al. [6] and in Scafetta and Bianchini [26].
Cionco and Compagnucci [137] discussed solar acceleration radial impulses in the Sun-barycenter direction and angular momentum inversion characterizing the last prolonged solar minima and Cionco and Soon [138] proposed both a precise calculation and a physical model for exchange between two rotating flux-mass elements that conserve angular momentum, also predicting future solar activity and grand solar minimum events.
One incorrect interpretation of the Sun's wobbling is to assume that the focus of the Earth's orbit is placed in the barycenter of the solar system, and that as a result, the sunlight that the Earth receives is heavily modulated by the solar wobbling (see for example Zharkova et al. [140]). This is incorrect, because the Earth's orbit is not Keplerian around the solar system's barycenter, which is not the gravitational attraction center of the solar system and, in particular, of the Earth. As a matter of fact, the Earth does not clearly see the solar wobbling, because it wobbles nearly simultaneously together with the Sun.

The "Homeopathic" Tiny-Tidal Problem
It is currently believed that a solar dynamo operating in the convective solar zone is responsible for fluctuations in solar activity [141]. A common critique against the possibility that the planetary motions could somehow have synchronized the solar dynamo activity and the climate on Earth is that the amplitude of the planetary tidal acceleration on the Sun (and/or the Earth) is very small [142,143]. This is a fact that has been known since the 19th century, which is still presented as if it were a conclusive argument against the planetary hypothesis. Moreover, Charbonneau [144] noted that the planetary tidal forcings could only have what he labeled a "homeopathic" effect on the solar tachocline and came to the conclusion that they should not be able to synchronize the solar dynamo. However, he also rightly stressed that the issue was not insurmountable, because it could be solved under the condition of discovering a solar mechanism that greatly amplifies the impact of the tidal forces.
Scafetta et al. [131] answered Jager and Versteegh [142] and Callebaut et al. [143] by arguing that such a kind of reasoning is only based on a classical Newtonian understanding of the tidal phenomenon. In reality, the most obvious objection to such a critique could be that more mechanisms instead of just a basic Newtonian one could be simultaneously at play. Therefore, the critique could not serve as a persuasive refutation of the planetary hypothesis of solar variability.
In fact, the ability of the Sun to generate large amounts of energy through H-burning could make possible a significant amplification of the tidal harmonic perturbations [74,145]. Under specific circumstances, these periodic perturbations might become strong enough to synchronize the solar dynamo-which is an oscillator in and of itself-to external planetary frequencies [26]. Thus, in order to examine and evaluate the planetary hypothesis of solar activity variation, it should be required to look into this option.
Scafetta [74] specifically addressed the tiny-tidal problem and proposed that the Sun may act as a nuclear fusion amplifier of the planetary tidal forcing, as previously proposed by Wolff and Patrone [145]. More specifically, Scafetta [74] calculated that any tidal gravitational power dissipated in the solar core might be up to four million times enhanced by a modulation of the H-burning rate. Using this amplification factor, the cumulative impact of all planetary tides on the Sun was estimated to increase total solar irradiance (TSI) by about 0.3-0.8 W/m 2 . The frequencies of tidally induced oscillations could be preserved even if a signal is produced by a tidally induced variation of the Hburning rate, which would changes the local temperature of the plasma, as first proposed by Scafetta [74]. In fact, in the radiative central region of the Sun, tidal forces produce gand p-mode oscillations [146][147][148] that can quickly propagate with the speed of sound (a few hundred kilometers per second) up to the tachocline layer, which is the base of the solar dynamo. As a result, the perturbations caused by the arriving g-waves may be strong enough to synchronize the solar dynamo with the planetary tidal frequencies. As a matter of fact, the original 4 million amplification factor of the dissipated tidal energy in the solar core predicted by the tidal-to-luminosity model of Scafetta [74] could produce more energetic g-waves that are capable to synchronize the solar dynamo to the frequencies of the planetary tides [26,74].
Scafetta [74] did not assert that the physical issue has been fully resolved. In any case, as detailed in Scafetta and Bianchini [26], the scientific literature also suggests several alternative physical mechanisms via which the planets may be modulating solar activity. These include torques, electromagnetic forcings, and the modulation of solar wind and interplanetary particle fluxes that may also directly influence the Earth's climate by modulating the albedo by affecting cloud formation [6,90,96,119,124,149]. Consequently, it should at least be acknowledged that efforts are being made to solve the "homeopathic" tiny-tidal problem and/or to develop alternative solutions to explain the apparent spectral coherence between solar activity and climate change oscillations and the frequency set emerging from the complex interactions of all the planets. In general, the current magneto-hydrodynamic dynamo models are still unsatisfactory because, although the models predict a periodicity in solar activity changes, they can only roughly mimic the Schwabe 11-year cycle by choosing ad hoc values for their free parameters [150,151]. The timing of the maxima and minima of the 11-year solar cycle is also neither hindcasted nor predicted, but results from the initial conditions imposed on the model equations, which are also free parameters. Similarly, all the other known solar oscillations, from the monthly to the Bray-Hallstatt timescales, cannot be obtained by the solar dynamo models that are now available. Figures 3, 7, 8 and 10 demonstrate that harmonic models based on planetary ephemeris and fundamental astronomical data appear to be capable of accurately predicting the period and timing of all the major observed solar oscillations. In this way, the first challenge posed by Callebaut et al. [143]-that planetary influences should be able to reproduce at least the most basic periodicities emerging in solar records-is met. Theses results appear to give the planetary theory an empirical advantage over the reductionist theory that assumes that solar activity changes are due to internal solar mechanisms alone. Here, we support the idea that the planetary hypothesis aims not just to substitute but to complement the solar dynamo models, as also commented by Charbonneau [31]. In fact, the existence of a solar dynamo that can naturally oscillate could be a required condition for solar activity to be synchronized by weak external periodic forcings, as the synchronization of pendulums discovered by Huygens in the 17th century demonstrates [152].
However, it is possible to argue that at least part of the observed solar and/or climatic oscillations could be more directly linked to space weather being modified by the oscillations of the magnetic and gravitational fields of the planets [6,7,90]. Space weather variations are connected to particle flux variations interacting with the physical and chemical properties of the terrestrial atmosphere. For example, they could modulate the cloud formation and, therefore, the Earth's albedo. Finally, the planets induce small oscillations in the orbital motion of the Earth [9,10,125,153,154]. These sub-Milankovic orbital oscillations cause only slight changes in the incoming solar radiation, which may nonetheless affect the climate.

Rebuttals of the Claims That Planetary Harmonics are Incompatible with Solar Activity Cycles
In Section 2, we showed that planetary harmonics are compatible with solar activity cycles. Here, the main arguments against the planetary theory of solar activity changes are succinctly discussed.
(1) Okal and Anderson [86] evaluated the power spectrum of a planetary tidal function made up of Mercury, Venus, Earth, and Jupiter and compared it with that of the sunspot number record. These authors concluded that the sunspot cycle had a periodicity of 11 years, while the tidal function only displayed the 11.86-year cycle associated with Jupiter's orbit. The difference between the two periods was significant. In essence, this critique confirms that Jupiter alone cannot fully account for the Schwabe 11-year sunspot cycle. The problem, however, could be solved in two ways (see Section 3.1): (a) by the Venus-Earth-Jupiter triple-syzygies tidal alignment models (some of them having been already proposed a century ago [28,29,[74][75][76]82,83]) and (b) by the tidal model function produced by Jupiter and Saturn, which includes the tidal oscillation related to the neap-spring tide of Jupiter and Saturn (9.93 years) and the orbital related tide of Jupiter (11.86 year), as discussed above [37,74,84].
(2) Smythe and Eddy [107] analyzed the planetary conjunctions during the Maunder minimum (1645-1715) and concluded that "the failure of planetary conjunction patterns to reflect the drastic drop in sunspots during the Maunder Minimum casts doubt on the tidal theory of solar activity". This critique was rebutted by Scafetta [37], who showed that grand solar minima, such as Oort, Wolf, Spörer, Maunder, and Dalton, likely derive from the interference among the neap-spring tide of Jupiter and Saturn (9.93 years), the orbital tide of Jupiter (11.86 years), and a period close to 11 years that represents the mean solar dynamo cycle; see Section 3.2.2.
(3) Callebaut et al. [143] noted that the solar activity records show multiple periodicities: the Schwabe 11-year sunspot cycle, the Hale 22-year magnetic cycle, the 88-year Gleissberg cycle, the De Vries (Suess) period of 203-208 years, and the Hallstatt cycle of about 2300 years, and we can add also the Eddy quasi-1000-year cycle. They claimed that it "should be remarked . . . that virtually none of the papers on planetary influences on solar variability succeeded in identifying these . . . periodicities in the planetary attractions". However, Callebaut et al. [155] later retracted their claim in response to the rebuttal by Scafetta et al. [131] and admitted that "it is well-known that there are some periodicities that are common to solar activity and planetary motions". See the various planetary modelings discussed in Section 2 and Scafetta et al. [6] and Scafetta [7].
(4) Cauquoin et al. [106] investigated the spectral properties of high-resolution records of 10 Be in the EPICA Dome C (EDC) ice core from Antarctica during the Marine Interglacial Stage 9.3 (MIS 9.3), 325-336 k years ago and concluded that there are "No evidence for planetary influence on solar activity 330,000 years ago". However, their power spectra showed four significant spectral peaks at about 103, 115, 130, and 150 years. Scafetta [48] showed that these frequencies emerge from the three-frequency model based on the tides of Jupiter and Saturn beating with the quasi-11-year solar dynamo cycle [37], once some nonlinear coupling between planetary and solar harmonics is assumed. In fact, these frequencies are those that generate Grand Solar Minima such as the Maunder and Dalton ones.
(5) Nataf [87] "built synthetic tidal forcings to illustrate the lack of evidence for a" quasi 11 "years periodicity" and concluded that "the astrological quest for a link between solar cycles and planetary tides remains as unfounded as ever". He also stated that it was "therefore surprising that this idea has received a renewed attention". Nataf actually made an even more serious error than that made by Okal and Anderson [86], because he ignored the role of Saturn and, in addition, assumed the orbits to be circular. In fact, his hypothesis that the orbits are circular is evidently erroneous. Nataf [87]'s critique was already rebutted by numerous works Scafetta [37,74], Hung [76], Tattersall [82], Wilson [83]; see Section 3.1 and also the formal recent rebuttal by Scafetta [84]. In his reply, Nataf [156] erroneously claimed that Saturn could be neglected because its tidal oscillation is 20 times smaller than that of Jupiter. However, Scafetta [84] was talking about the oscillation related to the neap-spring tide of Jupiter and Saturn (9.93 years), whose amplitude is only three times smaller than that related to Jupiter's orbital movement (11.86 year), as shown by the power spectra in Figure 5.

A Brief Response to Some Unscientific Critiques
Investigating open issues is the way for science to develop. Suitable physical mechanisms are identified only after empirical proofs. In fact, science moves from merely observing a phenomenon and finding empirical correlations among the observables to a physical understanding of it. However, even when only empirical data are available, this process must go on through trial and error until a suitable model is found, which might take a long time. For example, in antiquity, it was well known that ocean tides track the motion of the Moon [157]. In the early 17th century, Kepler formulated the three laws of planetary motion without using any physical mechanism. Only in 1687, Newton [158] postulated his universal gravitational law as a potential physical process for both the tides and Kepler's laws.
The history of science provides numerous examples where a "tiny-tide"-like objection was used to reject hypotheses based on empirical data, which were later fully confirmed. For example, in 1863, Lord Kelvin (William Thomson) rejected the link between sunspot activity and geomagnetic activity, such as auroras, as mere coincidences [159], because he reasoned that a solar magnetic field could not alter the Earth's magnetosphere, since the Sun and the Earth are too far away. In his 1892 presidential address to the Royal Society, Lord Kelvin reiterated the same prejudice, despite some scientists having already made discoveries suggesting that the Sun was expelling matter that might have transported the solar magnetic field over long distances. Today, the existence of solar wind is confirmed. Similarly, according to modern skeptics of the planetary synchronized solar dynamo hypothesis, the planets are considered to be too far away to have any impact on the Sun.
Another case is when, shortly before World War I, the German scientist Alfred Wegener [160] put forth a theory of continental drift based on a large number of topographical, sedimentary, and fossil patterns across continents. Wegener's theory, however, received unreasonable criticism by the majority of conservative scientists, who were hostile to any novel geophysical theory for the simple reason that he could not propose a fully convincing scientific justification for the migration of the continents. At that time, from a geophysical point of view, it was frequently objected that the solar and lunar tidal torques on the Earth's crust were too weak to move the continents. Once again, a large body of clear empirical facts was ignored or exotically interpreted, and the most obvious interpretation of the evidences was rejected only because its full physical explanation was still missing. As we know today, Wegener's work is fully acknowledged. Similarly, the many evidences supporting the planetary hypothesis of solar activity variations are ignored only because the underlying physical processes are still controversial.
Finally, we wish to comment on the tendency of dismissing as "astrology" the possibility that the solar dynamo itself might be partially synchronized by planetary harmonic forcings. We point out that the modern understanding of how the Sun, Moon, and planets can influence processes occurring on Earth developed through a rigorous examination of many ancient astrological claims regarding the origin of the seasons, ocean tides, weather events, climate changes, seismic activity, etc. This investigation was mostly conducted in the 16th and 17th centuries by professional astrologers, such as Kepler, who had the patience and skills to distinguish between the facts and the superstitions that tainted earlier astrological works, including the Tetrabiblos by Ptolemy [157], who already tried to build up a scientific astrology in opposition to a mythological view of natural phenomena. Kepler himself was convinced that celestial configurations of the Sun, Moon, and planets had a real physical impact on people by modifying the Earth's climate and weather [161]. In any case, ancient astrologers did discover correlations between astronomical configurations and events, such as the rise and fall of civilizations, probably associated with cyclical climate changes [162]. For example, important references were made about the 60-year and quasi-millennial cycles, which were associated with the Trigon of Great Conjunctions of Jupiter and Saturn and its millennial rotation through the constellations, know as the Great Inequality between Jupiter and Saturn [162,163]. It is also known that traditional Chinese and Indian calendars adopted a sexagesimal system associated with the cycle of Brihaspati, that is, Jupiter, probably because a quasi-60-year cycle has been observed in the monsoon rainfall cycles since antiquity [164][165][166]. As a matter of fact, the quasi-60-and 1000-year cycles, together with other planetary cycles, are observed both in climate and solar records, as shown here in Figures 8 and 10 [35,37,38,48,52,60,85,89,166]. Now, if it were demonstrated that the planets can somehow synchronize solar activity changes and, directly or indirectly, the space weather conditions around the Earth, a number of ancient astrologers' claims might find a physical interpretation and answer. For this reason, it would be wise to approach these topics with an open mind, as Kepler himself did. Figure 1 depicts how the Earth's climate has substantially changed over time. The most well-known climatic periodicities are associated with galactic cycles [2][3][4], various astronomical orbital forcings associated with Milankovitch cycles [5], and numerous solar activity variability cycles, ranging from decadal to multimillennial timescales. Since there are interannual solar, orbital, and soli-lunar tidal oscillations, the interannual timescale oscillations (e.g., the El Niño-Southern Oscillations) may also have a solar or astronomical origin [9][10][11]. Herein, we examined the solar activity variability cycles, their possible planetary origin, and their implications for climate variations, which range from decadal to multimillennial timescales during the Holocene.

Conclusions
The scientific evidences presented in Section 2 suggest that the Sun is subject to planetary influence from monthly to millennial timescales. The most frequently documented solar cycles are those with periods of 0. Many of these oscillation are also found in climate records and could be used to forecast an important component of the natural climate variability [48,85]. On the contrary, the idea that changes in solar activity might be only governed by internal dynamo mechanisms appears to still fall short of fully explaining why all of these solar cycles coincide with the most significant planetary periods derived from the planets' orbits, spring tides, synodic cycles, and invariant inequalities.
We acknowledge that the exact physical mechanisms underlying a planetary synchronization of solar activity are still being investigated and are somehow still mysterious, but the empirical evidences are quite compelling and cannot be easily discounted. Gravitational tides may be able to explain some of the observed behavior, especially if internal solar mechanisms dramatically increase their impact on the solar luminosity production rate, as suggested by Wolff and Patrone [145], Scafetta [74] and Scafetta and Bianchini [26]. Additional mechanisms may be at work. For example, certain planetary orbital configura-tions and conjunctions, even those that include tidally weak planets such as Uranus and Neptune, may be able to modify the interplanetary space weather conditions. For example, they could modulate the large scale structure of the heliospheric electromagnetic field and the particle fluxes coming from interstellar space (such as cosmic rays or even dark matter), or solar wind and the fluxes of interplanetary dust [6,7]. In addition, even small variations in the orbital motion of the Earth and lunar tides could induce some other oscillations in the climate system. Even if the precise nature of the physical processes that could account for the findings is not yet fully understood, it might be in the future. Thus, we call for additional investigation into these issues given the empirical data now accessible and the fact that the specific physical mechanisms underlying the solar cycles are still unknown. Figure 12 proposes a schematic representation of the planetary model of solar activity oscillations and climate changes, together with various possible involved mechanisms. A planetary hypothesis of solar activity could actually shed light on the physical mechanisms that may have gone unnoticed or underestimated in the past, and we advise searching for planetary-induced synchronization mechanisms of the solar dynamo or of the space weather inside the heliosphere. A planetary model of solar activity cycles would be extremely valuable for predicting changes in solar activity as well as climate oscillations on several timescales.   [52]: https:// www.ncei.noaa.gov/access/paleo-search/study/8744 9400-year solar activity reconstruction [38]: https://www.ncei.noaa.gov/access/paleo-search/study/12894 The ∆ 14 C (‰) record ( [111], Int-Cal04.14c): https://www.radiocarbon.org/IntCal04.htm.