On the Possible Contribution of Natural Climatic Fluctuations to the Global Warming of the Last 135 Years

A number of numerical experiments with artificial random signals (the second order autoregressive processes), which have important statistical properties similar to that of the observed instrumental temperature (1850-2015), were carried out. The results show that in frame of the selected mathematical model the return period of climatic events, analogous to the current global warming (linear increase of temperature for 0.95 ̊C during the last 135 years) is 2849-5180 years (one event per 2849-5180 years). This means that global warming (GW) of the last 135 years can unlikely be fully explained by inherent oscillations of the climatic system. It was found however, that natural fluctuations of climate may appreciably contribute to the GW. The return period of climatic episodes with 0.5 ̊C warming during the 135 years (half of the observed GW) was less than 500 years. The result testifies that the role of external factors (emission of greenhouse gases, solar activity etc.) in the GW could be less than often presumed.


Introduction
It is generally well-known that the mean surface temperature of the globe has been rising during more than last century [1] [2].It is widely accepted also that this GW is caused primarily by anthropogenic increase of greenhouse gases concentration [1] [2].However debates on this question still continue.Some experts maintain that current warming does not exceed the natural fluctuations of climate.Evidence of appreciable contribution to GW of non-greenhouse factors has been obtained by many authors.E.g.Soon et al. [3] noted that if the urbanization effect is properly taken into account, one can conclude that solar variability is the dominant factor of Northern Hemisphere long-term temperature changes since at least 1881.Zhao and Feng [4] reported that variations in solar activity play an important role in changes of climate over global scale during the last more than 100 years.According to Harde [5] the Sun is the main contributor to global warming of the last century.Lüdeсke et al. [6] showed that variations of Central European temperature after 1757 were likely governed by periodic oscillations resulted from intrinsic climatic dynamics.Scafetta [7] [8] claimed that the global climate oscillations from 1950 to 2011 were appreciably influenced by astronomic planetary cycles, particularly by motion of Jupiter and Saturn.Swanson and Tsonis [9] noted that in the period 1900-2000 Northern Hemisphere climate variability might be partly explained by chaotic dynamics.
Privalsky and Fortus [10] modeled variations of global temperature during 1850-2009 as an autoregressive process of the fourth order.They arrived at a conclusion that GW of the last 150 years could be fully explained by natural climatic variability without any external forcing.In the present work we continue research of possible contribution of intrinsic variations of climatic system to the GW.

Modeling of the Changes of the Global Temperature during the Last 135 Years
Privalsky and Fortus [10] modeled variations of the global temperature as a random autoregressive process: ( ) ( ) where ε j is a zero mean white noise and values of p = 4 and autoregressive coefficients j φ were determined in order to obtain Fourier spectrum of the modeled signal coincides with the spectrum of the actual instrumental series (1850-2009).
It should be noted, however, that the length of experimental temperature record is too short for accurate estimation of spectral power over the low-frequency domain (periods longer than 40 -50 years).Here we follow the approach of [10] but with substantial changes: 1) We used data on instrumental temperature anomaly ΔT i over 1880-2015 (http://cdiac.ornl.gov/ftp/trends/temp/hansen/gl_land_ocean.txt, Figure 1(a)).We did not use the data on global temperature during 1850-1880 because of their large uncertainty (see e.g.[11] [12]).The observed temperature shows monotonic increase, determined by its linear trend, for ca 0.95˚C over 1880-2015 (see Figure 1(a)).
2) We modeled the global temperature with the less intricate autoregressive process of the second order (AR(2)).
3) We adjusted our model so that the simulated signal reproduced the following statistical properties of the actual temperature record:

Results and Discussion
We performed statistical experiments with three different AR(2) processes of the 1000 years length.Their properties are detailed in Table 1.We looked for time intervals of 135 years length during which the modeled temperature shows linear increase for 0.95˚C or more-the natural global warming (NGW) events.We used Table 1.Results of statistical examination of frequency occurrence of NGW events.2000 simulations for each kind of AR( 2) process (Table 1).
In Table 1 1 that statistical properties of the simulated temperatures were actually very similar to that of the observed data set.N 1 is the number of cases in which one NGW episode took place during 1000 years.N 2 and N 3 are numbers of cases with two and three such events.The case of the simulated temperature which has two NGW episodes during 1000 years is shown in Figure 3. Statistical experiments performed using three autoregressive processes, which have statistical characteristics close to that Then we estimated probability of the part of the GW that was induced by intrinsic climatic variability.We repeated the statistical experiment but with the amplitude of temperature rise of 0.5˚C during 135 years-the Half Global Warming (HGW) episode.We found that the frequency of HGW events is 2.13 -2.34 per 1000 years and their return period is 426 -469 years.As to events analyzed by Privalsky and Fortus [10] (rise of temperature for 0.5˚C during 65 years) we estimated their return period as ca 280 years that is less than determined in [10] (489 years).These results testify that inherent variations of climatic system could contribute to the GW appreciably.If the GW is at least partly a result of natural variability of the climatic system this warm fluctuation might continue in the next few decades.Actually, we reveal that the frequency of warm episodes with the amplitude of temperature rise of 0.6˚C during 165 years is 1.5 -1.6 per 1000 years.

Conclusion
Numerical experiments with artificial signals (second order random processes), which have important statistical properties similar to that of the observed time series of temperature (1850-2015), indicate that the return period of climatic events analogous to the current GW (linear increase of temperature for 0.95˚C during 135 years) is 2849-5180 years.Correspondingly, the probability that the current global terrestrial temperature is going through the NGW episode is 0.026 -0.047.This shows that the GW unlikely could be fully explained by natural variability of temperature.Therefore it is reasonable to regard the GW as a phenomenon exceptional from the point of view of intrinsic climatic oscillations, which need an additional external forcing factor for explanation.On the other hand, the statistical experiments showed that an appreciable part of the GW might be a result of natural fluctuations of climatic system.It was found that that the return period of climatic events with 0.5˚C warming during 135 years is less than 500 years.Thus, the probability that the half of the contemporary GW is produced by natural random oscillations of temperature is more than 0.25.Emission of greenhouse gases (carbon dioxide, methane, nitrous oxide) most likely is another possible contributor to the GW of the last centuries [1] [2] [13].Changes in the solar radiation at the Earth's surface (global brightening) might be important source of the warming of the last decades [14].However, our results show that the contribution of these external factors (including greenhouse effect) to the GW could be less than is often believed [1] [2].

Figure 1 . 1 φ , 2 φ and variance of white noise 2 ε σ were determined so that the mean parameters 1 mσ
Figure 1.(a) The observed global temperature anomaly (thin line), its 30-year average (thick line), linear trend (dash-dot line); (b) High-frequency fluctuations of the observed global temperature (30-year average is subtracted).hf σ is the stan- dard deviation of high-frequency temperature variability.ΔTt rend -is the linear increase of instrumental temperature during 1880-2015.