On the Breaking of the Milankovitch Cycles Triggered by Temperature Increase: The Stochastic Resonance Response

: Recent decades have registered the hottest temperature variation in instrumentally recorded data history. The registered temperature rise is particularly signiﬁcant in the so-called hot spot or sentinel regions, characterized by higher temperature increases in respect to the planet average value and by more marked connected effects. In this framework, in the present work, following the climate stochastic resonance model, the effects, due to a temperature increase independently from a speciﬁc trend, connected to the 10 5 year Milankovitch cycle were tested. As a result, a breaking scenario induced by global warming is forecasted. More speciﬁcally, a wavelet analysis, innovatively performed with different sampling times, allowed us, besides to fully characterize the cycles periodicities, to quantitatively determine the stochastic resonance conditions by optimizing the noise level. Starting from these system resonance conditions, numerical simulations for increasing planet temperatures have been performed. The obtained results show that an increase of the Earth temperature boosts a transition towards a chaotic regime where the Milankovitch cycle effects disappear. These results put into evidence the so-called threshold effect, namely the fact that also a small temperature increase can give rise to great effects above a given threshold, furnish a perspective point of view of a possible future climate scenario, and provide an account of the ongoing registered intensity increase of extreme meteorological events. Incomplete and Unequally Spaced Time-Series.


Introduction
The history of climate science is quite old. The Greek term "klima" designated the inclination of the Sun's rays in respect to the Earth's surface, so testifying the understanding, already in those ancient times, of the correlation between the solar energy flux and the daily and seasonal temperature variations. Starting from the end of the 1800s, the progress in physics sciences has allowed researchers to interpret the fundamental processes involved in the air and ocean movements, as well as to evaluate the planet radiative balance controlled by the energy exchanges between the Earth and the external space. This latter depends in part on the Earth's surface temperature and in part on the greenhouse effect, firstly formulated by Joseph Fourier in 1824. The greenhouse effect has inspired investigations into the role of the CO 2 emitted by human activities on the planet climate and has led to a rough estimation of the planet temperature evolution by the Sweden chemist Svante Arrhenius.
In parallel, the work of naturalists, geologists, and glaciologists showed that the Earth's climate has deeply changed during this time. At the beginning of the twentieth century, the Serb mathematician Milutin Milankovitch formulated the hypothesis that the slow modifications in the Earth's orbit around the Sun were the cause of the great glaciations. This theory was confirmed at the end of the twentieth century, thanks to the In the literature, different methods to investigate signal periodicities exist, with specific reference to the Fourier (FT) and to the wavelet (WT) transforms.
The calculated Fourier power spectrum of the not equispaced data record, performed by means of the Lomb-Scargle method [27][28][29][30], is shown in Figure 2. Such an analysis method allows researchers to detect periodic features in a series characterized by a not constant sampling. In the literature, different methods to investigate signal periodicities exist, with specific reference to the Fourier (FT) and to the wavelet (WT) transforms.
The calculated Fourier power spectrum of the not equispaced data record, performed by means of the Lomb-Scargle method [27][28][29][30], is shown in Figure 2. Such an analysis method allows researchers to detect periodic features in a series characterized by a not constant sampling. In order to extract more information on the frequency components of the temperature variation as a function of time during the last 5.5 M years, in the following, we applied the wavelet analysis. It is well known that, when it is important to single out the signal frequency localizations, one can take a huge advantage from WT analysis [31][32][33][34][35][36]. In fact, differently from FT, which shows all the signal frequencies without a time localization, WT shows the component frequencies and at what time they occur [37][38][39].
In Figure 3, the obtained WT scalogram of the temperature data, obtained by using an analytical Morlet wavelet with parameter 6 [40][41][42] is shown. The performed WT analysis allowed us to draw the following conclusions: i.
Variance increases, characterized by high ice volumes and by low temperature values, are present with a period of 100 K years; these maxima find a correspondence with the Earth's orbit eccentricity period (first Milankovitch cycle); ii.
Variance increases are present in the temperature sequence with a period of 41 K years; these maxima correspond to the period of the Earth's axis inclination (second Milankovitch cycle); In order to extract more information on the frequency components of the temperature variation as a function of time during the last 5.5 M years, in the following, we applied the wavelet analysis. It is well known that, when it is important to single out the signal frequency localizations, one can take a huge advantage from WT analysis [31][32][33][34][35][36]. In fact, differently from FT, which shows all the signal frequencies without a time localization, WT shows the component frequencies and at what time they occur [37][38][39].
In Figure 3, the obtained WT scalogram of the temperature data, obtained by using an analytical Morlet wavelet with parameter 6 [40][41][42] is shown. In order to extract more information on the frequency components of the temperature variation as a function of time during the last 5.5 M years, in the following, we applied the wavelet analysis. It is well known that, when it is important to single out the signal frequency localizations, one can take a huge advantage from WT analysis [31][32][33][34][35][36]. In fact, differently from FT, which shows all the signal frequencies without a time localization, WT shows the component frequencies and at what time they occur [37][38][39].
In Figure 3, the obtained WT scalogram of the temperature data, obtained by using an analytical Morlet wavelet with parameter 6 [40][41][42] is shown. The performed WT analysis allowed us to draw the following conclusions: i.
Variance increases, characterized by high ice volumes and by low temperature values, are present with a period of 100 K years; these maxima find a correspondence with the Earth's orbit eccentricity period (first Milankovitch cycle); ii.
Variance increases are present in the temperature sequence with a period of 41 K years; these maxima correspond to the period of the Earth's axis inclination (second Milankovitch cycle); The performed WT analysis allowed us to draw the following conclusions: i. Variance increases, characterized by high ice volumes and by low temperature values, are present with a period of 100 K years; these maxima find a correspondence with the Earth's orbit eccentricity period (first Milankovitch cycle); ii.
Variance increases are present in the temperature sequence with a period of 41 K years; these maxima correspond to the period of the Earth's axis inclination (second Milankovitch cycle); iii.
Lower intensity variance increases with a period of about 23 K years, corresponding to the variation of the Earth's axis during a double-conic motion, i.e., to solar precession, are present.
Furthermore, the obtained WT scalogram features are in good agreement with the prominent features reported by Crucifix [43]. Figure 4 shows a sketch of the three Milankovitch cycles. In particular, Figure 4a shows the Earth's eccentricity variation characterized by a period of 100 K years. More specifically, the Earth's orbit eccentricity has varied from nearly 0 to 0.07, and today its value is about 0.017. The Earth receives 20-30% more radiation when it is in the perihelion than when it is at aphelion. In particular, the average distance Sun-Earth is 1.49 × 10 11 m, while the Sun's average diameter is about 110 times that of the Earth (12.750 × 10 3 m). Figure 4b displays the cycle of the Earth's axis inclination, ranging from a minimum of 21 • 55 to a maximum of 24 • 20 , with a period of 41 K years. The inclination of the Earth's axis relative to the orbit plane generates seasons. Slight changes in the tilt changes the amount of solar radiation falling on certain planet locations. The Earth's axis is nowadays tilted of~23.5 degrees, with a decreasing trend. Finally, Figure 4c shows the variation of the Earth's axis following a double-cone motion, i.e., a precession motion, with a 23 K years period. This precession motion does not change the tilt of Earth's axis, but its orientation. During the past millennia, the Earth's axis pointed more or less north, towards Polaris, known as the Pole Star. However, about 5000 years ago the Earth was more turned towards another star, called Thubin. Moreover, in about 12,000 years, the axis will point to Vega. After a cycle of precession, the orientation of the planet is changed with respect to the aphelion and the perihelion: if a hemisphere is pointed towards the sun during the perihelion, it will be opposite during the aphelion, the vice versa being true for the other hemisphere. As a result, the hemisphere that is pointed towards the sun during perihelion and away during aphelion experiences more extreme seasonal contrasts than the other hemisphere. At present, winter in the southern hemisphere occurs near the aphelion and summer near the perihelion, which means that the southern hemisphere experiences more extreme seasons than the northern hemisphere [44][45][46][47][48][49].
As matter of fact, as shown by the cornerstone work of Benzi et al. [50], the variations in the distribution of solar radiation during Milankovitch cycles, alone, cannot explain the wide decrease in temperature that marks the transition from an interglacial age to a glacial age. Therefore, it clearly emerges that to enhance the effect of these causes, which foretell a temperature variation of only 0.3 K, a positive feedback should be taken into account. On that score, to figure out this puzzle, a stochastic resonance effect has been postulated. This effect concerns with nonlinear systems energetically characterized by a threshold (i.e., by an activation barrier), and consists of a system output (response) to a weak input signal (cause), which undergoes a resonance-like behavior, i.e., to a performance maximum, as a function of the noise level, differently to what occurs in linear systems where the output is maximized in the absence of noise [51][52][53]. In fact, the complex climate system contains several interacting components that give rise to a highly nonlinear behavior through its many feedback mechanisms. iii. Lower intensity variance increases with a period of about 23 K years, correspondin to the variation of the Earth's axis during a double-conic motion, i.e., to solar preces sion, are present. Furthermore, the obtained WT scalogram features are in good agreement with th prominent features reported by Crucifix [43]. Figure 4 shows a sketch of the thre Milankovitch cycles. In particular, Figure 4a shows the Earth's eccentricity variation char acterized by a period of 100 K years. More specifically, the Earth's orbit eccentricity ha varied from nearly 0 to 0.07, and today its value is about 0.017. The Earth receives 20-30% more radiation when it is in the perihelion than when it is at aphelion. In particular, th average distance Sun-Earth is 1.49 10 11 m, while the Sun's average diameter is about 11 times that of the Earth (12.750 10 3 m). Figure 4b displays the cycle of the Earth's axi inclination, ranging from a minimum of 21°55' to a maximum of 24°20′, with a period o 41 K years. The inclination of the Earth's axis relative to the orbit plane generates seasons Slight changes in the tilt changes the amount of solar radiation falling on certain plane locations. The Earth's axis is nowadays tilted of ~23.5 degrees, with a decreasing trend Finally, Figure 4c shows the variation of the Earth's axis following a double-cone motion i.e., a precession motion, with a 23 K years period. This precession motion does not chang the tilt of Earth's axis, but its orientation. During the past millennia, the Earth's axi pointed more or less north, towards Polaris, known as the Pole Star. However, about 500 years ago the Earth was more turned towards another star, called Thubin. Moreover, i about 12,000 years, the axis will point to Vega. After a cycle of precession, the orientatio of the planet is changed with respect to the aphelion and the perihelion: if a hemispher is pointed towards the sun during the perihelion, it will be opposite during the aphelion the vice versa being true for the other hemisphere. As a result, the hemisphere that i pointed towards the sun during perihelion and away during aphelion experiences mor extreme seasonal contrasts than the other hemisphere. At present, winter in the souther hemisphere occurs near the aphelion and summer near the perihelion, which means tha the southern hemisphere experiences more extreme seasons than the northern hemispher [44][45][46][47][48][49].
(a) As matter of fact, as shown by the cornerstone work of Benzi et al. [50], the variations in the distribution of solar radiation during Milankovitch cycles, alone, cannot explain the wide decrease in temperature that marks the transition from an interglacial age to a glacial age. Therefore, it clearly emerges that to enhance the effect of these causes, which foretell a temperature variation of only 0.3 K, a positive feedback should be taken into account. On that score, to figure out this puzzle, a stochastic resonance effect has been postulated. This effect concerns with nonlinear systems energetically characterized by a threshold (i.e., by an activation barrier), and consists of a system output (response) to a weak input signal (cause), which undergoes a resonance-like behavior, i.e., to a performance maximum, as a function of the noise level, differently to what occurs in linear systems where the output is maximized in the absence of noise [51][52][53]. In fact, the complex climate system contains several interacting components that give rise to a highly nonlinear behavior through its many feedback mechanisms.
Climatic driving mechanisms acting on human evolution constitute an open problem of interdisciplinary scientific interest that triggers an intense debate concerning possible evolutionary scenarios. The concentration increase of CO2 and other greenhouse gases in atmosphere influences the climate system significantly. The major role in determining the warming rate and its spatial distribution is played by ocean water, due to its heat capacity. On one hand, the climate response to temperature variation is delayed by the ocean; on the other hand, ocean water delays the warming effect on climate. Global surface temperatures show short and long time scales in the atmosphere-ocean general circulation, although fast and slow responses in the ocean are different from those in the atmosphere.
In this paper we addressed the issue of climate variability caused by internal processes, by external forcing, and by their interaction, and discuss mechanisms of a global Earth temperature increase. In particular, we explore the stochastic resonance scenario under the condition of a temperature increase, independently from peculiar, specific temperature trends. The results of the performed numerical simulations show that, starting from a resonance condition, an increase in temperature gives rise to a transition towards a chaotic regime.

Methods Section
It is well known that the Fourier approach is based on a definite set of sinusoidal functions, which provides a unique representation in terms of a standard orthogonal basis. Its weaknesses include the assumptions of linearity and stationarity and the a priori Climatic driving mechanisms acting on human evolution constitute an open problem of interdisciplinary scientific interest that triggers an intense debate concerning possible evolutionary scenarios. The concentration increase of CO 2 and other greenhouse gases in atmosphere influences the climate system significantly. The major role in determining the warming rate and its spatial distribution is played by ocean water, due to its heat capacity. On one hand, the climate response to temperature variation is delayed by the ocean; on the other hand, ocean water delays the warming effect on climate. Global surface temperatures show short and long time scales in the atmosphere-ocean general circulation, although fast and slow responses in the ocean are different from those in the atmosphere.
In this paper we addressed the issue of climate variability caused by internal processes, by external forcing, and by their interaction, and discuss mechanisms of a global Earth temperature increase. In particular, we explore the stochastic resonance scenario under the condition of a temperature increase, independently from peculiar, specific temperature trends. The results of the performed numerical simulations show that, starting from a resonance condition, an increase in temperature gives rise to a transition towards a chaotic regime.

Methods Section
It is well known that the Fourier approach is based on a definite set of sinusoidal functions, which provides a unique representation in terms of a standard orthogonal basis. Its weaknesses include the assumptions of linearity and stationarity and the a priori choice of basis functions, which do not change in time. This means that the projection of the basis results in a 'global' frequency analysis that may not provide a good match with the time scales of the underlying processes.
Some techniques have been designed to address the issue of nonstationarity, so as to provide a simultaneous time-frequency analysis. Two of such techniques are windowed Fourier analysis and wavelets. In windowed Fourier analysis, a Fourier analysis is applied to a sliding window, whose length is shorter than the given record; then, a projection is applied to each windowed subset of the data.
Wavelet analysis is based on a set of orthogonal functions with compact support or approximately compact support, vanishing outside a finite interval. The interval is rescaled to highlight different time scales and the center of interval can be shifted to address nonstationarity. The algorithm determining the decomposition coefficients is still a projection, but the resulting coefficients depend on the shape, the width, and location of interval of support. A set of orthonormal functions can be constructed from a "mother" wavelet via dilation and translation. More specifically, from a formal point of view, a Fourier series decomposes a function f (t) defined for 0 ≤ t ≤ 2π and having a period T = 2π into a sum of sinusoidal functions, i.e., of cosine and/or sine functions: a 0 , a j , and b j being the so-called Fourier coefficients, evaluated through the integrals: where a 0 furnishes the average value of the function f (t) over a period T = 2π and where the Fourier coefficients, a j and b j , furnish the weight of each cosine and sine contributions; for odd functions all the a j coefficients are zero whereas for even functions all the b j coefficients are zero. Fourier transform (FT) constitutes an extension of the Fourier series, which furnishes huge advantages for functions f (t) that are nonperiodic, and therefore when T → ∞ , FT decomposes the function f (t) into sets of sine and/or cosine functions, or of exponential functions e 2πiωt , each of them with weightf (ω), with continuously variable frequencies ω: As a matter of fact, FT suffers from the limitation that it does not allow a connection between the frequency spectrum and the signal evolution in time and, therefore, to overcome this difficulty, wavelet (WT) analysis has been formulated [54][55][56].
WT results are more advantageous in respect to Fourier transform; in fact, the WT approach allows users to decompose a signal into its wavelets components, by means of mother wavelet ψ: where the parameter a > 0 denotes the scale and its value corresponds to the inverse of the frequency, the parameter τ represents the shift of the scaled wavelet along the time axis, ψ * represents the complex conjugation of the mother wavelet, while the scaled and shifted mother wavelet ψ a,τ (t) can be expressed by: Differently from FT, which shows only which signal frequencies are present, WT, in addition, also shows where, or at what scale, they are. Furthermore, while FT allows researchers to decompose the signal only in cosine and sine component functions, the WT takes into account several wavelet mother functions.
More precisely, different types of mother wavelet can be used to perform a wavelet analysis. As a rule, mother wavelets are characterized by properties such as orthogonality, compact support, and symmetry. Although different mother wavelet functions satisfy the previous conditions, when applied to the same signal produce different results. For these reasons, the choice of the best mother wavelet is crucial. In the present study, a comparison among the different correlation degree between the signal and different mother wavelet functions has been performed [57]. As a result, in this work the wavelet Morlet has been chosen as mother wavelet: where ω 0 represents the center pseudo-frequency and σ provides the wavelet bandwidth [58,59]. It should be noticed that in the special case in which the wavelet mother is ψ(t) = e −2πit , the WT transform reduces to the FT. For this study, Matlab environment has been used.

The Physical Approach to Climate Change through the Stochastic Resonance Modelization
The stochastic resonance phenomenon takes place when the introduction of a noise in a nonlinear system characterized by a form of threshold amplifies a weak signal. In particular, it occurs that, by increasing the noise amplitude, the system performance parameters, e.g., the output signal-to-noise ratio (SNR), show a resonance. Therefore, such a behavior is different in respect to linear systems, where the SNR is maximized in the absence of noise.
The time evolution of the climate system configuration can be described by means of a configuration variable x(t), by means of a potential U(x) and of a random force f r (t) through the Langevin equation: where ε sin (ω 0 t) is a weak forcing term corresponding, within the climate model, to the 100,000-year Milankovitch cycle; f r is a white noise, of intensity I, which corresponds to the fast dynamics and whose autocorrelation function is < f r (t) f r (0) ≥ 2Iδ(t). In the following, we will assume that this latter hypothesis holds for the climate system. A simple example of potential U(x) is furnished by the symmetric bistable potential introduced by McNamara et al. [60]: which has two stable equilibrium configurations at x ± = ±c = ± a b , and an unstable equilibrium configuration at x = 0; the two stable minima correspond to the glacial and interglacial temperature states. In the absence of the forcing term, i.e., for ε = 0, for assigned system parameters, the configuration variable x(t) fluctuates in time between the two minima with a variance that is proportional to I, whereas in the presence of a periodic forcing term and of specific noise values, the transition probability of escaping from each well varies regularly in time. In particular, what emerges from the stochastic resonance modelling is that, for specific noise values, jumps between the two minima occur in synchrony with the frequency of the weak periodic signal.

Numerical Simulation of Climate Changes and Results
The climate system has a chaotic nature, where a very small perturbation can be amplified. This chaotic character is taken into account both explicitly, by introducing a noise term, or indirectly, by introducing small variations in the initial conditions or in some physical parameters; in this latter case, one usually copes with a huge number of runs for the same prevision and performs a probabilistic analysis [61][62][63][64][65][66]. In order to modelize the climate cycle behavior and the effects connected with the global temperature trend, in this section, we adapted the stochastic differential Langevin's equation, which originally took into account only a systematic and a random contribution, to a dynamic case in which, besides a friction and a stochastic term, an external contribution and a double well potential were present [67][68][69][70][71]: In other terms, the climate system is modelled as a particle of mass m, with a given average energy and configuration coordinate x(t), which experiences a force potential V 1 (x), and on which both a friction force, f f riction = −ζ dx(t) dt , i.e., a Stokes force, and a random force, f random = ξ(t), connected with the system "heat bath", are exerted [72][73][74][75][76][77][78][79].
In the "high friction limit", i.e., when the system's acceleration is negligible, by including the constant ζ into the potential function V 2 (x) = V 1 (x) ζ and into the random force ξ(t) = ξ 0 (t) ζ , Equation (11) becomes: Let us now introduce an explicit form for a generalized potential V(x, t), conceived in order to include a simple double-well potential and a periodic contribution: Under this hypothesis, Equation (12) can be written as: where the quantities α and η(t) allow researchers to describe with a greater detail the stochastic noise ξ(t). Such a mathematical model allows users to describe the planet temperature behaviour when one identifies x(t) the Earth's temperature T(t) [80,81].
V(x, t) is the Earth's climate potential, which takes into account the energy balance between the incoming R in and outcoming R out radiation and the two system levels describing the glacial (T 1 ∼ = 278.6 K) and interglacial (T 2 ∼ = 283.3 K) states [82][83][84]. In particular, it is: R in = Q(1 + Acos(ωt)), where Q is the solar constant and A is the amplitude of the modulation, estimated, by Milankovitch, to be equal to 5 × 10 −4 , with ω = 2π T = 2π 100000 years . Furthermore, R out = µ(T)R in + E(T), where µ(T) is the albedo, i.e., the reflected solar radiation which, for simplicity, has a linear decreasing trend with temperature, and E(T) = 4πR 2 σT 4 is the energy radiated by the Earth according to its temperature T. Then, it is c dT dt = − ∂V ∂T = R in − R out = Q(1 + A cos(ωt))(1 − µ(T)) − E(T) (c being the Earth's thermal capacity). Finally, the stochastic noise represents the short time scale processes connected with atmospheric dynamics and circulation. In short, in terms of temperatures value, we have: As shown by Benzi et al. [85], depending on the α values, in spite of the weakness of the periodic contribution ax cos (ωt), which in our modelization corresponded to the first Milankovitch cycle, it is possible to get an amplification of the signal, so justifying the remarkable gap of almost 10 • C between glacial and interglacial periods.
In the present work, with the aim to single out the effects of temperature increase on the above described SR scenario, at first we performed numerical simulations on a system characterized by a double-well potential, V = T 4 − T 2 , which allowed us to describe the glacial and interglacial states. In the model drawing up, we have taken into account a small periodic disturbance, ax cos(ωt), which allowed us to describe the Earth's orbit eccentricity cycle, and a noise level of amplitude √ α fixed on the SR condition. In this framework, starting from the obtained resonance conditions, we have performed new numerical simulations, taking into account increasing global Earth temperature values. This temperature increase has been taken into account by downshifting the double-well potential of a constant term V 0 with respect to the fixed total energy of the system, i.e., V = T 4 − T 2 − V 0 ; in such a case, in fact, although the driving force keeps constant, see Equation (15), the downshift of the potential well, keeping constant the total energy, corresponds to an increase of the system kinetic energy, and hence of temperature. In other words, the simulations have been performed at first by varying the intensity of the noise while keeping the periodic signal constant in order to single out the noise amplitude value, for which the jump of the system from one hole to the other one is synchronized with the external system periodic disturbance; in such a way, the resonance condition was obtained [86,87]. Successively, starting from this SR condition, we have simulated the modelled system answer for different values of the global planet temperature. Figure 5 reports the solution of Equation (14), i.e., the x(t) = T(t), by using the following parameters: a = 0.46 and α = 0.23. As it can be seen, under such a situation, the system jumped from the left well to the right well, giving rise to a periodic oscillation between the two wells corresponding to the glacial (at temperature T 1 ) and interglacial (at temperature T 2 ) states [88][89][90].
formed new numerical simulations, taking into account increasing global Earth temperature values. This temperature increase has been taken into account by downshifting the double-well potential of a constant term with respect to the fixed total energy of the system, i.e., = − − ; in such a case, in fact, although the driving force keeps constant, see Equation (15), the downshift of the potential well, keeping constant the total energy, corresponds to an increase of the system kinetic energy, and hence of temperature. In other words, the simulations have been performed at first by varying the intensity of the noise while keeping the periodic signal constant in order to single out the noise amplitude value, for which the jump of the system from one hole to the other one is synchronized with the external system periodic disturbance; in such a way, the resonance condition was obtained [86,87]. Successively, starting from this SR condition, we have simulated the modelled system answer for different values of the global planet temperature. Figure  5 reports the solution of Equation (14), i.e., the x(t)=T(t), by using the following parameters: = 0.46 and = 0.23. As it can be seen, under such a situation, the system jumped from the left well to the right well, giving rise to a periodic oscillation between the two wells corresponding to the glacial (at temperature T1) and interglacial (at temperature T2) states [88][89][90]. Figure 5. Solution x(t) = T(t) of the generalized Langevin's equation for the following parameters: a = 0.46 and α = 0.23. As it can be seen, under such conditions, the system jumped from the left well to the right well, giving rise to a periodic oscillation between the two wells, with the same period of the external disturbance, which correspond to the glacial, at temperature T1, and interglacial, at temperature T2, states.
To spotlight the SR condition and to quantitatively determine the signal-to-noise ratio, SNR, a WT scalogram of the obtained solution, i.e., x(t) = T(t), has been evaluated for the fixed value of a = 0.46 and for different values of α, i.e. α = 0.05 (left side), 0.23 (center), and 0.70 (right side), as reported in Figure 6. . Solution x(t) = T(t) of the generalized Langevin's equation for the following parameters: a = 0.46 and α = 0.23. As it can be seen, under such conditions, the system jumped from the left well to the right well, giving rise to a periodic oscillation between the two wells, with the same period of the external disturbance, which correspond to the glacial, at temperature T 1 , and interglacial, at temperature T 2 , states.
To spotlight the SR condition and to quantitatively determine the signal-to-noise ratio, SNR, a WT scalogram of the obtained solution, i.e., x(t) = T(t), has been evaluated for the fixed value of a = 0.46 and for different values of α, i.e., α = 0.05 (left side), 0.23 (center), and 0.70 (right side), as reported in Figure 6.  Such a WT approach allowed us to evaluate the SNR value, for different noise levels [91][92][93][94][95]. The result of such a calculation for 12 values of noise, in the range 0.00~0.83, is shown in Figure 7. These results put in the foreground that, starting from low noise levels, at first, the SNR value surged, attaining a maximum of a = 0.46 and α = 0.23 and then it dropped for higher noise values. Therefore, for a = 0.46, the value of α = 0.23 established a cutting-edge between increasing and decreasing SNR conditions [96][97][98][99][100][101]. Such a WT approach allowed us to evaluate the SNR value, for different noise levels [91][92][93][94][95]. The result of such a calculation for 12 values of noise, in the range 0.00∼0.83, is shown in Figure 7. These results put in the foreground that, starting from low noise levels, at first, the SNR value surged, attaining a maximum of a = 0.46 and α = 0.23 and then it dropped for higher noise values. Therefore, for a = 0.46, the value of α = 0.23 established a cutting-edge between increasing and decreasing SNR conditions [96][97][98][99][100][101].
Such a WT approach allowed us to evaluate the SNR value, for different noise levels [91][92][93][94][95]. The result of such a calculation for 12 values of noise, in the range 0.00~0.83, is shown in Figure 7. These results put in the foreground that, starting from low noise levels, at first, the SNR value surged, attaining a maximum of a = 0.46 and α = 0.23 and then it dropped for higher noise values. Therefore, for a = 0.46, the value of α = 0.23 established a cutting-edge between increasing and decreasing SNR conditions [96][97][98][99][100][101]. During the last few decades, there has been considerable progress in understanding climate variability by large observational efforts, model simulations, and theory development. In the next section, following the proposed SR climate model, we present the results of new numerical simulations that will show how, starting from a resonance condition, an increase of temperature gives rise to transition to a chaotic regime. During the last few decades, there has been considerable progress in understanding climate variability by large observational efforts, model simulations, and theory development. In the next section, following the proposed SR climate model, we present the results of new numerical simulations that will show how, starting from a resonance condition, an increase of temperature gives rise to transition to a chaotic regime.

Effects of the Global Climate Warming on the Basis of the SR Climate Model
A climate change is considered anomalous when it is not compatible with only its internal variability within a period of time longer than about 30 years. The evidence of a global climate change is based on the analysis of the long-term variability that affects the different components of the climate system, which are the atmosphere (i.e., the shell of gases and vapors around the Earth), the hydrosphere (i.e., seas, rivers, water basins), the cryosphere (i.e., sea ices, snows, glaciers, iced soils), the biosphere (i.e., the living organisms), and the lithosphere (i.e., Earth's crust). These elements continuously exchange energy and matter (i.e., water, mineral and/or organic substances) between each other; furthermore, their interactions are modulated by planet radiative balance fluctuations, which can be amplified or stabilized by positive or negative feedback physical processes (positive and negative retroactions). Atmosphere constitutes the rapid component of the climate system, since the air masses mix, at the planetary level, within some months, whereas the seas interact with the atmosphere from both short time scales (i.e., daily) through to several decades. Sea currents are due to the Earth's rotation, winds, and seas density, which in turn depends on sea temperature and salinity. Climate changes depend on external forcing, i.e., perturbations external to the climate systems (e.g., orbital parameters, variations of Sun's activity), or by internal forcing (e.g., El Nino). Furthermore, anthropic activities produce greenhouse gases, such as carbon dioxide (CO 2 ), methane (CH 4 ), nitrogen protoxyde (N 2 O), etc.; aerosols (small solid or liquid particles); and changes in the soil usage (e.g., deforestation, agriculture).
The climate trend can be inferred from the analysis of long-term datasets, which are essentially based on paleoclimate archives of different climate-dependent chemicalphysical quantities, often integrated through numerical and interpretative models. In fact, climate changes leave indirect information in the physical, chemical, and biological compositions of natural materials, i.e., in proxy data. In particular, the observations and the available measurements of the atmosphere, of the oceans, and of the Earth's surface parameters indicate unequivocally that since the nineteenth century, the planet has started to undergo to a warming process. In recent decades, soil and ocean temperatures, on average, have increased, even if not homogeneously, but with greater warming near polar latitudes, especially in the Arctic region.
In the meantime, the ice and snow coverage has decreased in different areas of the planet. The atmospheric water vapor content has increased everywhere due to the higher temperature of the atmosphere, while sea level has also risen. Furthermore, changes in several other climatic indicators have also been detected; these include the seasonal extension in many places of the planet; an increase in the tendency of extreme meteorological events, such as heat waves or heavy rainfall; and a general decrease in extreme cold events. As far as the cause of the registered warming is concerned, if one takes into account only the natural forcing, it is not possible to justify these climatic changes. Most of the warming effects observed globally over the last 50 years can only be explained if one takes into account the effects produced by human activities, and in particular both the emissions from fossil fuels and the deforestation processes [102][103][104][105].
The series of global temperatures registered over the period 1840-2017 show a warming trend interrupted only for short periods. In particular, the decade between 1900 and 1910 and that between 1942 and 1950 represent the only time intervals in which global temperatures have undergone a reversal trend. Moreover, starting from the late 1970s, the slope of the regression line has undergone a substantial increase, which has brought the heating values between 1970 and 2017 to about 0.17 • C/decade [106][107][108][109][110][111].
In this framework, with the aim to single out the effects of temperature increase on the above described SR scenario, starting from the resonance condition obtained at a = 0.46 for α = 0.23, we next focused our attention on the SR model outcomes for increasing values of the Earth's temperature. More specifically, to cope with the effects of global warming within the SR scenario, we performed new numerical simulations for a sequence of downshifted double-well potentials (see Figure 8).
Climate 2021, 9, x FOR PEER REVIEW 13 of 17 Figure 8. Downshifting of the double well potential, which, starting from the resonance condition, keeping constant the total energy, corresponded to an increase of the system kinetic energy and hence of temperature.
In fact, a temperature increase corresponded to a downshift of the system potential energy, furnishing a dummy run of the system behaviour. As an example, Figure 9 shows three WT scalograms of the generalized Langevin's equation solution for = 0, = 1, and = 2.5. As it can be seen, the Earth's temperature increase gave rise to a breaking of the stochastic resonance condition. Figure 8. Downshifting of the double well potential, which, starting from the resonance condition, keeping constant the total energy, corresponded to an increase of the system kinetic energy and hence of temperature.
In fact, a temperature increase corresponded to a downshift of the system potential energy, furnishing a dummy run of the system behaviour. As an example, Figure 9 shows three WT scalograms of the generalized Langevin's equation solution for V 0 = 0, V 0 = 1, and V 0 = 2.5. As it can be seen, the Earth's temperature increase gave rise to a breaking of the stochastic resonance condition.
In fact, a temperature increase corresponded to a downshift of the system potential energy, furnishing a dummy run of the system behaviour. As an example, Figure 9 shows three WT scalograms of the generalized Langevin's equation solution for = 0, = 1, and = 2.5. As it can be seen, the Earth's temperature increase gave rise to a breaking of the stochastic resonance condition. Figure 9. WT scalograms of the generalized Langevin's equation solution for = 0, = 1, and = 2.5. As it can be seen, the temperature increase gave rise to a breaking of the stochastic resonance condition.
In summary, in this study we have tested the climate stochastic resonance model under the condition of a temperature increase. More specifically, within this framework model, we have determined the resonance stochastic condition by varying at first the noise level and by applying the wavelet analysis; then, once the resonance condition has been reached, by keeping constant the noise level, new simulations were performed by varying the thermal potential. As a result, the simulation results indicate a transition towards a chaotic regime, where the Milankovitch cycle effects disappear.

Conclusions
The last few decades have registered the hottest temperature variation in recorded history. In this framework, the effects provided by a global temperature increase on climate behavior, as inferred by both applying the stochastic resonance interpretative model for climate and WT data analysis, are discussed.
In particular, a Milankovitch-cycle-breaking scenario induced by a triggering temperature increase, independently from a peculiar temperature trend, is forecasted by the climate stochastic resonance model. Numerical simulations for increasing planet temperatures have been performed. The simulations' outcomes beget the following scenario: Figure 9. WT scalograms of the generalized Langevin's equation solution for V 0 = 0, V 0 = 1, and V 0 = 2.5. As it can be seen, the temperature increase gave rise to a breaking of the stochastic resonance condition.
In summary, in this study we have tested the climate stochastic resonance model under the condition of a temperature increase. More specifically, within this framework model, we have determined the resonance stochastic condition by varying at first the noise level and by applying the wavelet analysis; then, once the resonance condition has been reached, by keeping constant the noise level, new simulations were performed by varying the thermal potential. As a result, the simulation results indicate a transition towards a chaotic regime, where the Milankovitch cycle effects disappear.

Conclusions
The last few decades have registered the hottest temperature variation in recorded history. In this framework, the effects provided by a global temperature increase on climate behavior, as inferred by both applying the stochastic resonance interpretative model for climate and WT data analysis, are discussed.
In particular, a Milankovitch-cycle-breaking scenario induced by a triggering temperature increase, independently from a peculiar temperature trend, is forecasted by the climate stochastic resonance model. Numerical simulations for increasing planet temperatures have been performed. The simulations' outcomes beget the following scenario: starting from the resonance condition, the Earth's temperature increase boosts a transition towards a chaotic regime, where the Milankovitch cycles disappear. It should be noticed that, although meteorology and climate drastically differ for the involved time scales, they constitute the two faces of the same coin. In this framework, the obtained results justify, on short time scales, the registered intensity increase of extreme meteorological events. If the greenhouse gas emissions continue to increase at the present rate during the XXI century, a temperature increase of 5 ± 1 • C is expected by the year 2100 in respect to the value corresponding to the beginning of the industrial era; this would also imply a rise of the sea level of 75 ± 20 cm. For these reasons, as well as in consideration of the climate system high inertia connected to the residence time of CO 2 in the atmosphere, it is important to stabilize the global temperature within +2 • C by 2100 compared to the preindustrial value, in agreement to the international treatises.
Author Contributions: M.T.C. and S.M. contributed equally; in particular M.T.C. and S.M. performed theoretical and numerical analysis, contributed analytic tools and interpreted the data, wrote the main draft of the paper, and participated in revising the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflict of interest.