Next Article in Journal
Mathematics and Poetry • Unification, Unity, Union
Next Article in Special Issue
How to Achieve Compliance with GDPR Article 17 in a Hybrid Cloud Environment
Previous Article in Journal
Changing Sci from Post-Publication Peer-Review to Single-Blind Peer-Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Atmospheric Temperature and CO2: Hen-Or-Egg Causality?

by
Demetris Koutsoyiannis
1,* and
Zbigniew W. Kundzewicz
2
1
Department of Water Resources and Environmental Engineering, School of Civil Engineering, National Technical University of Athens, 157 80 Athens, Greece
2
Institute for Agricultural and Forest Environment, Polish Academy of Sciences, 60-809 Poznań, Poland
*
Author to whom correspondence should be addressed.
Submission received: 7 September 2020 / Accepted: 16 November 2020 / Published: 25 November 2020
(This article belongs to the Special Issue Feature Papers 2020 Editors' Collection)

Abstract

:
It is common knowledge that increasing CO2 concentration plays a major role in enhancement of the greenhouse effect and contributes to global warming. The purpose of this study is to complement the conventional and established theory, that increased CO2 concentration due to human emissions causes an increase in temperature, by considering the reverse causality. Since increased temperature causes an increase in CO2 concentration, the relationship of atmospheric CO2 and temperature may qualify as belonging to the category of “hen-or-egg” problems, where it is not always clear which of two interrelated events is the cause and which the effect. We examine the relationship of global temperature and atmospheric carbon dioxide concentration in monthly time steps, covering the time interval 1980–2019 during which reliable instrumental measurements are available. While both causality directions exist, the results of our study support the hypothesis that the dominant direction is T → CO2. Changes in CO2 follow changes in T by about six months on a monthly scale, or about one year on an annual scale. We attempt to interpret this mechanism by involving biochemical reactions as at higher temperatures, soil respiration and, hence, CO2 emissions, are increasing.

Graphical Abstract

Πότερον ἡ ὄρνις πρότερον ἢ τὸ ᾠὸν ἐγένετο (Which of the two came first, the hen or the egg?).
Πλούταρχος, Hθικά, Συμποσιακὰ Β, Πρόβλημα Γ (Plutarch, Moralia, Quaestiones convivales, B, Question III).

1. Introduction

The phrase “hen-or-egg” is a metaphor describing situations where it is not clear which of two interrelated events or processes is the cause and which the effect. Plutarch was the first to pose this type of causality as a philosophical problem using the example of the hen and the egg, as indicated in the motto above. We note that in the original Greek text, “ἡ ὄρνις” is feminine (article and noun), meaning the hen rather than the chicken. Therefore, here, we preferred the form “hen-or-egg” over “chicken-or-egg”, which is more common in English (Very often, in online Greek texts, e.g., [1,2], “ἡ ὄρνις” appears as “ἡ ἄρνις”. After extended searching, we contend that this must be an error, either an old one in copying of manuscripts, e.g., by monks in monasteries, or a modern one, e.g., in OCR. We are confident that the correct word is “ὄρνις”).
The objective of the paper is to demonstrate that the relationship of atmospheric CO2 and temperature may qualify as belonging to the category of “hen-or-egg” problems. First, we discuss the relationship between temperature and CO2 concentration by revisiting intriguing results from proxy data-based palaeoclimatic studies, where the change in temperature leads and the change in CO2 concentration follows. Next, we discuss the databases of modern (instrumental) measurements related to global temperature and atmospheric CO2 concentration and introduce a methodology to analyse them. We develop a stochastic framework, introducing useful notions of time irreversibility and system causality while we discuss the logical and technical complications in identifying causality, which prompt us to seek just necessary, rather than sufficient, causality conditions. In the Results section, we examine the relationship of these two quantities using the modern data, available at the monthly time step. We juxtapose time series of global temperature and atmospheric CO2 concentration from several sources, covering the common time interval 1980–2019. In our methodology, it is the timing rather than the magnitude of changes that is important, being the determinant of causality. While logical, physically based arguments support the “hen-or-egg” hypothesis, indicating that both causality directions exist, interpretation of cross-correlations of time series of global temperature and atmospheric CO2 suggests that the dominant direction is T → CO2, i.e., the change in temperature leads and the change in CO2 concentration follows. We attempt to interpret this latter mechanism by noting the positive feedback loop—higher temperatures increase soil respiration and, hence, CO2 emissions.
The analysis reported in this paper was prompted by observation of an unexpected (and unfortunate) real-world experiment: during the COVID-19 lockdown in 2020, despite the unprecedented decrease in carbon emissions (Figure 1), there was an increase in atmospheric CO2 concentration, which followed a pattern similar to previous years (Figure 2). Indeed, according to the International Energy Agency (IEA) [3], global CO2 emissions were over 5% lower in the first quarter of 2020 than in that of 2019, mainly due to an 8% decline in emissions from coal, 4.5% from oil, and 2.3% from natural gas. According to other estimates [4], the decrease is even higher: the daily global CO2 emissions decreased by 17% by early April 2020 compared with the mean 2019 levels, while for the whole 2020, a decrease between 4% and 7% is predicted. Despite that, as seen in Figure 2, the normal pattern of atmospheric CO2 concentration (increase until May and decrease in June and July) did not change. Similar was the behaviour after the 2008–2009 financial crisis, but the most recent situation is more characteristic because the COVID-19 decline in 2020 is the most severe ever, even when considering the periods corresponding to World Wars. It is also noteworthy that, as shown in Figure 1, there are three consecutive years in the 2010s where there are no major increases, in emissions while there was an increase in CO2 concentration. (At first glance, this does not sound reasonable and we have therefore cross-checked the data with other sources, namely the Global Carbon Atlas [5], and the database of Our World In Data [6]; we found only slight differences.) Interestingly, Figure 1 also shows a rapid growth in emissions after the 2008–2009 global financial crisis, which is in agreement with Peters et al. [7].

2. Temperature and Carbon Dioxide—From Arrhenius and Palaeo-Proxies to Instrumental Data

Does the relationship of atmospheric carbon dioxide (CO2) and temperature classify as a “hen-or-egg”-type causality? If we look at the first steps of studying the link between the two, the reply is clearly negative. Arrhenius (1896, [8]), the scientist most renowned for studying the causal relationship between two quantities, regarded the changes in atmospheric carbon dioxide concentration as the cause and the changes in temperature as the effect. Specifically, he stated:
Conversations with my friend and colleague Professor Högbom together with the discussions above referred to, led me to make a preliminary estimate of the probable effect of a variation of the atmospheric carbonic acid [meaning CO2] on the temperature of the earth. As this estimation led to the belief that one might in this way probably find an explanation for temperature variations of 5–10 °C, I worked out the calculation more in detail and lay it now before the public and the critics.
Furthermore, following the Italian meteorologist De Marchi (1895, [9]), whom he cited, he rejected what we call today Milanković cycles as possible causes of the glacial periods. In addition, he substantially overestimated the role of CO2 in the greenhouse effect of the Earth’s atmosphere. He calculated the relative weights of absorption of CO2 and water vapour as 1.5 and 0.88, respectively, or a ratio of 1:0.6.
Arrhenius [8] also stated that “if the quantity of carbonic acid increases in geometric progression, the augmentation of the temperature will increase nearly in arithmetic progression”. This Arrhenius’s “rule” (which is still in use today) is mathematically expressed as
T T 0 = α ln ( [ CO 2 ] [ CO 2 ] 0 )
where T and [CO2] denote temperature and CO2 concentration, respectively, T0 and [CO2]0 represent reference states, and α is a constant.
Here, it is useful to note that Arrhenius’s studies were not the first on the subject. Arrhenius [8] cites several other authors, among whom Tyndall (1865, [10]) for pointing out the enormous importance of atmospheric absorption of radiation and for having the opinion that water vapour has the greatest influence on this. Interestingly, it appears that the first experiments on the ability of water vapour and carbon dioxide to absorb heat were undertaken even earlier by Eunice Newton Foote (1856, [11]), even though she did not recognize the underlying mechanisms or even the distinction of short- and long-wave radiation [12,13,14]).
While the fact that the two variables are tightly connected is beyond doubt, the direction of the simple causal relationship needs to be studied further. Today, additional knowledge has been accumulated, particularly from palaeoclimatic studies, which allow us to examine Arrhenius’s hypotheses on a sounder basis. In brief, we can state the following:
  • Indeed, CO2 plays a substantial role as a greenhouse gas. However, modern estimates of the contribution of CO2 to the greenhouse effect differ largely from Arrhenius’s results, attributing 19% of the long-wave radiation absorption to CO2 against 75% of water vapour and clouds (Schmidt et al. [15]), or a ratio of 1:4.
  • During the Phanerozoic Eon, Earth’s temperature varied by even more than 5–10 °C, which was postulated by Arrhenius—see Figure 3. Even though the link of temperature and CO2 is beyond doubt, this is not clear in Figure 3, where it is seen that the CO2 concentration has varied by about two orders of magnitude and does not always synchronize with the temperature variation. Other factors may become more important at such huge time scales. Thus, an alternative hypothesis of the galactic cosmic ray flux as a climate driver via solar wind modulation has been suggested [16,17], which has triggered discussion or dispute [14,18,19,20,21,22,23]. The T–CO2 relationship becomes more legible and rather indisputable in proxy data of the Quaternary (see Figure 4). It has been demonstrated in a persuasive manner by Roe [24] that in the Quaternary, it is the effect of Milanković cycles (variations in eccentricity, axial tilt, and precession of Earth’s orbit), rather than of atmospheric CO2 concentration, that explains the glaciation process. Specifically (quoting Roe [24]),
    variations in atmospheric CO2 appear to lag the rate of change of global ice volume. This implies only a secondary role for CO2—variations in which produce a weaker radiative forcing than the orbitally-induced changes in summertime insolation—in driving changes in global ice volume.
Despite falsification of some of Arrhenius’s hypotheses, his line of thought remained dominant. Yet, there have been some important studies, based on palaeoclimatological reconstructions (mostly the Vostok ice cores [25,26]), which have pointed to the opposite direction of causality, i.e., the change in temperature as the cause and that of CO2 concentration as the effect. Such claims have been based on the fact that temperature change leads and CO2 concentration change follows. In agreement with Roe [24], several papers have found the time lag to be positive, with estimates varying from 50 to 1000 years or more, depending on the time period and the particular study [27,28,29,30,31,32]. Claims that CO2 concentration leads (i.e., a negative lag) have not generally been made in these studies. At most, a synchrony claim has been sought on the basis that the estimated positive lags are often within the 95% uncertainty range [31], while in one publication [29], it has been asserted that a “short lead of CO2 over temperature cannot be excluded”. With respect to the last deglacial warming, Liu et al. [32], using breakpoint lead–lag analysis, again find positive lags and conclude that the CO2 is an internal feedback in Earth’s climate system rather than an initial trigger.
Figure 3. Proxy-based reconstructions of global mean temperature and CO2 concentration during the Phanerozoic. The temperature reconstruction by Scotese [33] was mainly based on proxies from [21,34,35,36], while the CO2 concentration proxies have been taken from Davis [37], Berner [38], and Ekart et al. [39]; all original figures were digitized in this study. The chronologies of geologic eras shown in the bottom of the figure have been taken from the International Commission on Stratigraphy [40].
Figure 3. Proxy-based reconstructions of global mean temperature and CO2 concentration during the Phanerozoic. The temperature reconstruction by Scotese [33] was mainly based on proxies from [21,34,35,36], while the CO2 concentration proxies have been taken from Davis [37], Berner [38], and Ekart et al. [39]; all original figures were digitized in this study. The chronologies of geologic eras shown in the bottom of the figure have been taken from the International Commission on Stratigraphy [40].
Sci 02 00083 g003
Figure 4. (upper) Time series of temperature and CO2 concentration from the Vostok ice core, covering part of the Quaternary (420,000 years) with time step of 1000 years. (lower) Auto- and cross-correlograms of the two time series. The maximum value of the cross-correlation coefficient, marked as ◆, is 0.88 and appears at lag 1 (thousand years) (adapted from Koutsoyiannis [30]).
Figure 4. (upper) Time series of temperature and CO2 concentration from the Vostok ice core, covering part of the Quaternary (420,000 years) with time step of 1000 years. (lower) Auto- and cross-correlograms of the two time series. The maximum value of the cross-correlation coefficient, marked as ◆, is 0.88 and appears at lag 1 (thousand years) (adapted from Koutsoyiannis [30]).
Sci 02 00083 g004
Since palaeoclimatic data suggest a direction opposite to that assumed by Arrhenius, Koutsoyiannis [30], using palaeoclimatic data from the Vostok ice cores at a time resolution of 1000 years and a stochastic framework similar to that of the present study (see Section 4.1), concluded that a change in temperature precedes that of CO2 by one time step (1000 years), as illustrated in Figure 4. He also noted that this causality condition holds for a wide range of time lags, up to 26,000 years, and, hence, the time lag is positive and most likely real. He asserted that the problem is obviously more complex than that of exclusive roles of cause and effect, classifying it as a “hen-or-egg” causality problem. Obviously, however, the proxy character of these data and the overly large time step of the time series reduce the reliability and accuracy of the results.
Studies exploring the rich body of modern datasets have also been published. Most of the studies have been based on the so-called “Granger causality test” (see Section 4.2). To mention a few, Kodra et al. [41], after testing several combinations and lags within the Granger framework, did not find any statistically significant results at the usual 5% significance level (they only found 2 cases at the 10% significance level; see their Tables 2 and 3). Stern and Kaufmann [42] studied, again within the Granger framework, the causality between radiative forcing and temperature, and found that both natural and anthropogenic forcings cause temperature change, and also that the inverse is true, i.e., temperature causes greenhouse gas concentration changes. They concluded that their results
show that properly specified tests of Ganger [sic] causality validate the consensus that human activity is partially responsible for the observed rise in global temperature and that this rise in temperature also has an effect on the global carbon cycle.
By contrast, Stips et al. [43] used a different method [44] to investigate the causal structure and concluded that their
study unambiguously shows one-way causality between the total Greenhouse Gases and GMTA [global mean surface temperature anomalies]. Specifically, it is confirmed that the former, especially CO2, are the main causal drivers of the recent warming.
Here, we use a different path to study the causal relation between temperature and CO2 concentration with the emphasis given on the exploratory and explanatory aspect of our analyses. While we occasionally use the Granger statistical test, this is not central in our approach. Rather, we place the emphasis on time directionality in the relationship, which we try to identify in the simplest possible manner, i.e., by finding the lag, positive or negative, which maximizes the cross-correlation between the two processes (see Section 4.1). We visualize our results by plots, so as to be simple, transparent, intuitive, readily understandable by the reader, and hopefully persuading. For the algorithmic-friendly reader, we also provide statistical testing results which just confirm what is directly seen in the graphs.
Another difference of our study, from most of the earlier ones, is our focus on changes, rather than current states, in the processes we investigate. This puts the technique of process differencing in central place in our analyses. This technique is quite natural and also powerful for studying time directionality [30]. We note that differencing has also been used in a study by Humlum et al. [45], which has several similarities with our study, even though it is not posed in a formal causality context, as well as in the study by Kodra et al. [41]. However, differencing has been criticized for potentially eliminating long-run effects and, hence, providing information only on short-run effects [42,46]. Even if this speculation were valid, it would not invalidate the differencing technique for the following reasons:
  • The short-term effects deserve to be studied, as well as the long-term ones.
  • The modern instrumental records are short themselves and only allow the short-term effects to be studied.
  • For the long-term effects, the palaeo-proxies provide better indications, as already discussed above.

3. Data

Our investigation of the relationship of temperature with concentration of carbon dioxide in the atmosphere is based on two time series of the former process and four of the latter. Specifically, the temperature data are of two origins, satellite and ground-based. The satellite dataset, developed at the University of Alabama in Huntsville (UAH), infers the temperature, T, of three broad levels of the atmosphere from satellite measurements of the oxygen radiance in the microwave band using advanced (passive) microwave sounding units on NOAA and NASA satellites [47,48]. The data are publicly available on the monthly scale in the forms of time series of “anomalies” (defined as differences from long-term means) for several parts of earth as well as in maps. Here, we use only the global average on monthly scale for the lowest level, referred to as the lower troposphere. The ground-based data series we use is the CRUTEM.4.6.0.0 global T2m land temperature [49]. This originates from a gridded dataset of historical near-surface air temperature anomalies over land. Data are available for each month from January 1850 to the present. The dataset is a collaborative product of the Met Office Hadley Centre and the Climatic Research Unit at the University of East Anglia. We note that both sources of information, UAH and CRUTEM, provide time series over the globe, land, and oceans; here, we deliberately use one source for the globe and one for the land.
The two temperature series used in the study are depicted in Figure 5. They are consistent with each other (and correlated, r = 0.8), though the CRUTEM4 series shows a larger increasing trend than the UAH series. The differences can be explained by three reasons: (a) the UAH series includes both land and sea, while the chosen CRUTEM4 series is for land only, in which the increasing trend is substantially higher than in sea; (b) the UAH series refers to some high altitude in the troposphere (see details in Koutsoyiannis [50]), while the CRUTEM4 series refers to the ground level; and (c) the ground-based CRUTEM4 series might be affected by urbanization (many ground stations are located in urban areas). In any case, the difference in the increasing trends is irrelevant for the current study, as the timing, rather than the magnitude, of changes is the determinant of causality. This will manifest in our results.
The most famous CO2 dataset is that of Mauna Loa Observatory [51]. The Observatory, located on the north flank of Mauna Loa Volcano on the Big Island of Hawaii, USA, at an elevation of 3397 m above sea level, is a premier atmospheric research facility that has been continuously monitoring and collecting data related to the atmosphere since the 1950s. The NOAA also has other stations that systematically measure atmospheric CO2 concentration, namely at Barrow, Alaska, USA and at South Pole. The NOAA’s Global Monitoring Laboratory Carbon Cycle Group also computes global mean surface values of CO2 concentration using measurements of weekly air samples from the Cooperative Global Air Sampling Network. The global estimate is based on measurements from a subset of network sites. Only sites where samples are predominantly of well-mixed marine boundary layer air, representative of a large volume of the atmosphere, are considered (typically at remote marine sea level locations with prevailing onshore winds). Measurements from sites at high altitude (such as Mauna Loa) and from sites close to anthropogenic and natural sources and sinks are excluded from the global estimate. (Details about this dataset are provided in [52]).
The period of data coverage varies, but all series cover the common 40-year period 1980–2019 which, hence, constitute the time reference of all our analyses. As a slight exception, the Barrow (Alaska) and South Pole measurements have not yet been available in final form for 2019 and, thus, this year was not included in our analyses of these two time series. The data of the latter two stations are given in irregular-step time series, which was regularized (by interpolation) to monthly in this study. All other data series have already been available on a monthly scale. While some of the earlier studies refer to a longer time span (e.g., [42,43] which start from 1850s), here, we avoid using non-systematic data earlier than 1980 due to their low reliability and bypass the raised controversies explained in Appendix A.1.
All four CO2 time series used in the study are depicted in Figure 6. They show a superposition of increasing trends and annual cycles whose amplitudes increase as we head from the South to the North Pole. The South Pole series has opposite phase of oscillation compared to the other three.
The annual cycle is better seen in Figure 7, where we have removed the trend with standardization, namely by dividing each monthly value by the geometric average of the preceding 5-year period. The reason why we used division rather than subtraction and geometric rather than arithmetic average (being thus equivalent to subtracting or averaging the logarithms of CO2 concentration) will become evident in Section 5. In the right panel of Figure 7, which depicts monthly statistics of the time series of the left panel, it is seen that in all sites but the South Pole, the annual maximum occurs in May; that of the South Pole occurs in September.

4. Methods

4.1. Stochastic Framework

A recent study [30] investigated time irreversibility in hydrometeorological processes and developed a theoretical framework in stochastic terms. It also studied necessary conditions for causality, which is tightly linked to time irreversibility. A simple definition of time reversibility within stochastics is the following, where underlined symbols denote stochastic (random) variables and non-underlined ones denote values thereof or regular variables.
A stochastic process x _ ( t ) at continuous time t, with nth order distribution function:
F ( x 1 , x 2 , , x n ; t 1 , t 2 , , t n ) P { x _ ( t 1 ) x 1 , x _ ( t 2 ) x 2 , , x _ ( t n ) x n }
is time-symmetric or time-reversible if its joint distribution does not change after reflection of time about the origin, i.e., if for any n, t 1 , t 2 , , t n ,
F ( x 1 , x 2 , , x n ; t 1 , t 2 , , t n ) = F ( x 1 , x 2 , , x n ; t 1 , t 2 , , t n ) .
If times t i   are   equidistant ,   i . e . ,   t i t i 1 = D , the definition can be also written by reflecting the order of points in time, i.e.,
F ( x 1 , x 2 , , x n 1 , x n ; t 1 , t 2 , , t n 1 , t n ) = F ( x 1 , x 2 , , x n 1 , x n ; t n , t n 1 , , t 2 , t 1 ) .
A process that is not time-reversible is called time-asymmetric, time-irreversible, or time-directional. Important results related to time (ir)reversibility are the following:
  • A time reversible process is also stationary (Lawrance [53]).
  • If a scalar process x _ ( t ) is Gaussian (i.e., all its finite dimensional distributions are multivariate normal) then it is reversible (Weiss [54]). The consequences are (a) a directional process cannot be Gaussian; (b) a discrete-time ARMA process (and a continuous-time Markov process) is reversible if and only if it is Gaussian.
  • However, a vector (multivariate) process can be Gaussian and irreversible at the same time. A multivariate Gaussian linear process is reversible if and only if its autocovariance matrices are all symmetric (Tong and Zhang [55]).
Time asymmetry of a process can be studied more conveniently (or even exclusively in a scalar process) through the differenced process, i.e.,
x ˜ _ τ , ν x _ τ + ν x _ τ
for an appropriate time-step ν of differencing. The differenced process represents change of the original process within a time period of length ν. We further define the cumulative process of x _ τ for discrete time κ as
X _ κ x _ 1 + x _ 2 + + x _ κ .
Through this, we find that the time average of the original process x _ τ for discrete time scale κ is
x _ τ ( κ ) x _ ( τ 1 ) κ + 1 + x _ ( τ 1 ) κ + 2 + + x _ τ κ κ = X _ τ κ X _ ( τ 1 ) κ κ .
Similar equations for the cumulative and averaged processes for the differenced process x ˜ _ τ , ν are given in Appendix A.2.
The variance of the process x _ τ ( κ ) is a function of the time scale κ, which is termed the climacogram of the process:
γ κ var [ x _ τ ( κ ) ] .
The autocovariance function for time lag η is derived from the climacogram through the relationship [56]
c η = ( η + 1 ) 2 γ | η + 1 | + ( η 1 ) 2 γ | η 1 | 2 η 2 γ | η | .
For sufficiently large κ (theoretically as κ → ∞), we may approximate the climacogram as
γ κ κ 2 H 2
where H is termed the Hurst parameter. The theoretical validity of such (power-type) behaviour of a process was implied by Kolmogorov (1940 [57]). The quantity 2H—2 is visualized as the slope of the double logarithmic plot of the climacogram for large time scales. In a random process, H = 1/2, while in most natural processes, 1/2 ≤ H ≤ 1, as first observed by Hurst (1951 [58]). This natural behaviour is known as (long-term) persistence or Hurst–Kolmogorov (HK) dynamics. A high value of H (approaching 1) indicates enhanced presence of patterns, enhanced change, and enhanced uncertainty (e.g., in future predictions). A low value of H (approaching 0) indicates enhanced fluctuation or antipersistence (sometimes misnamed as quasi-periodicity as the period is not constant).
For a stationary stochastic process x _ τ ,   the   differenced   process   x ˜ _ τ has mean zero and variance:
γ ˜ ν , 1 var [ x ˜ _ τ , ν ] = var [ x _ κ + ν ] + var [ x _ τ ] 2 cov [ x _ τ + ν , x _ τ ] = 2 ( γ 1 c ν )
where γ 1   and   c ν are the variance and lag ν autocovariance, respectively, of x _ τ . Furthermore, it has been demonstrated [30] that the Hurst coefficient of the differenced process x ˜ _ τ precisely equals zero, which means that x ˜ _ τ is completely antipersistent, irrespective of γ κ .
To study irreversibility in vector processes, we can use second-order moments and, in particular, cross-covariances among the different components of the vector. In particular (adapting and simplifying the analyses and results in Koutsoyiannis, [30]), given two processes x _ τ   and   y _ τ , we could study the cross-correlations:
r x ˜ y ˜ [ ν , η ] corr [ x ˜ _ τ , ν , y ˜ _ τ + η , ν ] .
Time (ir)reversibility could then be characterized by studying the properties of symmetry or asymmetry of r x ˜ y ˜ ( ν , η ) as a function of the time lag η. In a symmetric bivariate process, r x ˜ y ˜ [ ν , η ] = r x ˜ y ˜ [ ν , η ] , and if the two components are positively correlated, the maximum of r x ˜ y ˜ [ ν , η ]   will   appear   at   lag   η = 0 . If the bivariate process is irreversible, this maximum will appear at a lag η 1 0 and its value will be r x ˜ y ˜ [ ν , η 1 ] .
Time asymmetry is closely related to causality, which presupposes irreversibility. Thus, “no causal process (i.e., such that of two consecutive phases, one is always the cause of the other) can be reversible” (Heller, [59]; see also [60]). In probabilistic definitions of causality, time asymmetry is determinant. Thus, Suppes [61] defines causation thus: “An event Bt [occurring at time t′] is a prima facie cause of the event At [occurring at time t] if and only if (i) t < t ,   ( ii )   P { B t } > 0 , (iii) P ( A t | B t ) > P ( A t ) ″. In addition, Granger’s [62] first axiom in defining causality reads, “The past and present may cause the future, but the future cannot”.
Consequently, in simple causal systems, in which the process component x _ τ is the cause of y _ τ (like in the clear case of rainfall and runoff, respectively), it is reasonable to expect r x ˜ y ˜ [ ν , η ] 0 for any η 0 , while r x ˜ y ˜ [ ν , η ] = 0 for any η = 0 . However, in “hen-or-egg” causal systems, this will not be the case, and we reasonably expect r x ˜ y ˜ [ ν , η ] 0 for any η . Yet, we can define a dominant direction of causality based on the time lag η 1 maximizing cross-correlation. Formally, η 1 is defined for a specified ν as
η 1 arg max η | r x ˜ y ˜ ( ν , η ) | .
We can thus distinguish the following three cases:
  • If η1 = 0, then there is no dominant direction.
  • If η 1 > 0 , then the dominant direction is x _ τ y _ τ .
  • If η 1 < 0 , then the dominant direction is y _ τ x _ τ .
Justification and further explanations of these conditions are provided in Appendix A.3.

4.2. Complications in Seeking Causality

It must be stressed that the above conditions are considered as necessary and not sufficient conditions for a causative relationship between the processes x _ τ   and   y _ τ . Following Koutsoyiannis [30] (where additional necessary conditions are discussed), we avoid seeking sufficient conditions, a task that would be too difficult or impossible due to its deep philosophical complications as well as the logical and technical ones.
Specifically, it is widely known that correlation is not causation. As Granger [62] puts it,
when discussing the interpretation of a correlation coefficient or a regression, most textbooks warn that an observed relationship does not allow one to say anything about causation between the variables.
Perhaps that is the reason why Suppes [61] uses the term “prima facie cause” in his definition given above which, however, he does not explain, apart for attributing “prima facie” to Jaakko Hintikka. Furthermore, Suppes discusses spurious causes and eventually defines the genuine cause as a “prima facie cause that is not spurious”; he also discusses the very existence of genuine causes which under certain conditions (e.g., in a Laplacean universe) seems doubtful.
Granger himself also uses the term “prima facie cause”, while Granger and Newbold [63] note that a cause satisfying a causality test still remains prima facie because it is always possible that, if a different information set were used, it would fail the new test. Despite the caution issued by its pioneers, including Granger, through the years, the term “Granger causality” has become popular (particularly in the so-called “Granger causality test”, e.g., [64]). Probably because of that misleading term, the technique is sometimes thought of as one that establishes causality, thus resolving or overcoming the “correlation is not causation” problem. In general, it has rarely been understood that identifying genuine causality is not a problem of choosing the best algorithm to establish a statistical relationship (including its directionality) between two variables. As an example of misrepresentation of the actual problems, see [65], which contains the statement:
Determining true causality requires not only the establishment of a relationship between two variables, but also the far more difficult task of determining a direction of causality.
In essence, the “Granger causality test” studies the improvement of prediction of a process y _ τ by considering the influence of a “causing” process x _ τ through the Granger regression model:
y _ τ = j = 1 η a j y _ τ j + j = 1 η b j x _ τ j + ε τ
where a j and b j are the regression coefficients and ε τ is an error term. The test is based on the null hypothesis that the process x _ τ is not actually causing y _ τ , formally expressed as
H 0 : b 1 = b 2 = = b η = 0
Algorithmic details of the test are given in [64], among others. The rejection of the null hypothesis is commonly interpreted in the literature with a statement that x _ τ “Granger-causes” y _ τ .
This is clearly a misstatement and, in fact, the entire test is based on correlation matrices. Thus, it again reflects correlation rather than causation. The rejection of the null hypothesis signifies improvement of prediction and this does not mean causation. To make this clearer, let us consider the following example: people sweat when the atmospheric temperature is high and also wear light clothes. Thus, it is reasonably expected that in the prediction of sweat quantity, temperature matters. In the absence of temperature measurements (e.g., when we have only visual information, like when watching a video), the weight of the clothes algorithmically improves the prediction of sweat quantity. However, we could not say that the decrease in clothes weight causes an increase in sweat (the opposite is more reasonable and becomes evident in a three-variable regression, temperature–clothes weight–sweat, as further detailed in Appendix A.4).
Cohen [66] suggested replacing the term “Granger causality” with “Granger prediction” after correctly pointing out this:
Results from Granger causality analyses neither establish nor require causality. Granger causality results do not reveal causal interactions, although they can provide evidence in support of a hypothesis about causal interactions.
To avoid such philosophical and logical complications, here, we replace the “prima facie” or “Granger” characterization of a cause and, as we already explained, we abandon seeking for genuine causes by using the notion of necessary conditions for causality. One could say that if two processes satisfy the necessary conditions, then they define a prima facie causality, but we avoid stressing that as we deem it unnecessary. Furthermore, we drop “causality” from “Granger causality test”, thus hereinafter calling it “Granger test”.
Some have thought they can approach genuine causes and get rid of the caution “correlation is not causation” by replacing the correlation with other statistics in the mathematical description of causality. For example, Liang [44] uses the concept of information (or entropy) flow (or transfer) between two processes; this method has been called “Liang causality” in the already cited work he co-authors [43]. The usefulness of such endeavours is not questioned, yet their vanity to determine genuine causality is easy to infer: It suffices to consider the case where the two processes, for which causality is studied, are jointly Gaussian. It is well known that in any multivariate Gaussian process, the covariance matrix (or the correlation matrix along with the variances) fully determines all properties of the multivariate distribution of any order. For example, the mutual information in a bivariate Gaussian process is (Papoulis, [67])
H [ y _ | x ] = ln σ y 2 π e ( 1 r 2 )
where σ and r denote standard deviation and correlation, respectively. Thus, using any quantity related to entropy (equivalently, information), is virtually identical to using correlation. Furthermore, in Gaussian processes, whatever statistic is used in describing causality, it is readily reduced to correlation. This is evident even in Liang [44], where, e.g., in his Equation (102), the information flow turns out to be the correlation coefficient multiplied by a constant. In other words, the big philosophical problem of causality cannot be resolved by technical tricks.
From what was exposed above (Section 4.1), the time irreversibility (or directionality) is most important in seeking causality. In this respect, we certainly embrace Suppes’s condition (i) and Granger’s first axiom, as stated above. Furthermore, we believe there is no meaning in refusing that axiom and continuing to speak about causality. We note though that there have been recent attempts to show that
coupled chaotic dynamical systems violate the first principle of Granger causality that the cause precedes the effect.
[68]
Apparently, however, the particular simulation experiment performed in the latter work which, notably, is not even accompanied by any attempt for deduction based on stochastics, cannot show any violation. In our view, such a violation, if it indeed happened, would be violation of logic and perhaps of common sense.
Additional notes for other procedures detecting causality, which are not included in the focus of our study, are given in Appendix A.4.

4.3. Additional Clarifications of Our Approach

After the above theoretical and methodological discourse, we can clarify our methodological approach by emphasizing the following points.
  • To make our assertions and, in particular, to use the “hen-or-egg” metaphor, we do not rely on merely statistical arguments. If we did that, based on our results presented in the next section, we would conclude that only the causality direction T → [CO2] exists. However, one may perform a thought experiment of instantly adding a big quantity of CO2 to the atmosphere. Would the temperature not increase? We believe it would, as CO2 is known to be a greenhouse gas. The causation in the opposite direction is also valid, as will be discussed in Section 6, “Physical Interpretation”. Therefore, we assert that both causality directions exist, and we are looking for the dominant one under the current climate conditions (those manifest in the datasets we use) instead of trying to make assertions of an exclusive causality direction.
  • While we occasionally use statistical tests (namely, the Granger test, Equations (14) and (15)), we opt to use, as the central point of our analyses, Equation (13) (and the conditions below it) because it is more intuitive and robust, fully reflects the basic causality axiom of time precedence, and is more straightforward, transparent (free of algorithmic manipulations), and easily reproducible (without the need for specialized software).
  • For simplicity, we do not use any statistic other than correlation here. We stress that the system we are examining is indeed classified as Gaussian and, thus, it is totally unnecessary to examine any statistic in addition to correlation. The evidence of Gaussianity is provided by Figure A1 and Figure A2 in Appendix A.5, in terms of marginal distributions of the processes examined and in terms of their relationship. In particular, Figure A2 suggests a typical linear relationship for the bivariate process. We note that the linearity here is not a simplifying assumption or a coincidence as there are theoretical reasons implying it, which are related to the principle of maximum entropy [67,69].
  • All in all, we adhere to simplicity and transparency and, in this respect, we illustrate our results graphically, so they are easily understandable, intuitive, and persuasive. Indeed, our findings are easily verifiable even from simple synchronous plots of time series, yet we also include plots of autocorrelations and lagged cross-correlation, which are also most informative in terms of time directionality.

5. Results

5.1. Original Time Series

Here, we examine the relationship of atmospheric temperature and carbon dioxide concentration using the available modern data (observations rather than proxies) in monthly time steps, as described in Section 3. To apply our stochastic framework, we must first make the two time series linearly compatible. Specifically, based on Arrhenius’s rule (Equation (1)), we take the logarithms of CO2 concentration while we keep T untransformed. Such a transformation has also been performed in previous studies, which consider the logarithm of CO2 concentration as a proxy of total radiative forcing (e.g., [41]). However, by calling this quantity “forcing”, we indirectly give it, a priori (i.e., before investigating causation), the role of being the cause. Therefore, here, we avoid such interpretations; we simply call this variable the logarithm of carbon dioxide concentration and denote it as ln [ CO 2 ] .
A synchronous plot of the two processes (specifically, UAH temperature and ln [ CO 2 ] at Mauna Loa) is depicted in Figure 8. Very little can be inferred from this figure alone. Both processes show increasing trends and thus appear as positively correlated. On the other hand, the two processes appear to have different behaviours. Temperature shows an erratic behaviour while ln [ CO 2 ] has a smooth evolution marked by the annual periodicity. It looks impossible to infer causality from that graph alone.
Somewhat more informative is Figure 9, which depicts lagged cross-correlations of the two processes, based on the methodology in Section 4.1 but without differencing the processes. Specifically, Figure 9 shows the cross-correlogram between UAH temperature and Mauna Loa ln [ CO 2 ] at monthly and annual scales; the autocorrelograms of the two processes are also plotted for comparison. In both time scales, the cross-correlogram shows high correlations at all lags, with the maximum attained at lag zero. This does not hint at a direction. However, the cross-correlations for negative lags are slightly greater than those in the positive lags. Notice that to make this clearer, we have also plotted the differences r j r j in the graph. This behaviour could be interpreted as supporting the causality direction [CO2] → T. However, we deem that the entire picture is spurious as it is heavily affected by the fact that the autocorrelations are very high and, in particular, those of ln [ CO 2 ] are very close to 1 for all lags shown in the figure.
In our investigation, we also applied the Granger test on these two time series in both time directions. To calculate the p-value of the Granger test, we used free software (namely the function GRANGER_TEST [70,71]). It appears that in the causality direction [CO2] → T, the null hypothesis is rejected at all usual significance levels. The attained p-value of the test is 1.8 × 10−7 for one regression lag (η = 1), 1.8 × 10−4 for η = 2, and remains below 0.01 for subsequent η. By contrast, in the direction T → [CO2], the null hypothesis is not rejected at all usual significance levels. The attained p-value of the test is 0.25 for η = 1, 0.22 for η = 2, and remains above 0.1 for subsequent η.
Therefore, one could directly interpret these results as unambiguously showing one-way causality between the total greenhouse gases and temperature and, hence, validating the consensus view that human activity is responsible for the observed rise in global temperature. However, these results are certainly not unambiguous and, most probably, they are spurious. To demonstrate that they are not unambiguous, we have plotted, as shown in the upper panels of Figure 10, the p-values of the Granger test for moving windows with a size of 10 years for number of lags η = 1 and 2. The values for the entire length of time series, as given above, are also shown as dashed lines. Now the picture is quite different: each of the two directions appear dominating (meaning that the attained significance level is lower in one over the other) in about equal portions of the time. For example, for η = 2, the T → [CO2] dominates over [CO2] → T for 58% of the time. The attained p-value for direction T → [CO2] is lower than 1% for 1.4% of the time, much higher than in the opposite direction (0.3% of the time). All of these observations favour the T → [CO2] direction.
To show that the results are spurious and, in particular, affected by the very high autocorrelations of ln [ CO 2 ] and, more importantly, by its annual cyclicity, we have “removed” the latter by averaging over the previous 12 months. We did that for both series and plotted the results in the lower panels of Figure 10. Here, the results are stunning. For both lags η = 1 and 2 and for the entire period (or almost), T → [CO2] dominates, attaining p-values as low as in the order of 10−33. However, we will avoid interpreting these results as unambiguous evidence that the consensus view (i.e., human activity is responsible for the observed warming) is wrong. Rather, what we want to stress is that it is inappropriate to draw conclusions from a methodology which is demonstrated to be so sensitive to the used time windows and data processing assumptions. In this respect, we have included this analyses in our study only (a) to show its weaknesses (which, for the reasons we explained in Section 4.2, we believe would not change if we used different statistics or different time series) and (b) to connect our study to earlier ones. For the sake of drawing conclusions, we contend that our full methodology in Section 4.1 and Section 4.3 is more appropriate. We apply this methodology in Section 5.2.

5.2. Differenced Time Series

We have already explained the advantages of investigating the differenced processes, which quantify changes from a mathematical and logical point of view. In our case, taking differences is also physically meaningful as both CO2 concentration and temperature (equivalent to thermal energy) represent “stocks”, i.e., stored quantities and, thus, the mass and energy fluxes are indeed represented by differences.
The chosen time step of differencing is equal to one year (ν = 12 for the monthly time step of the time series). For instance, from the value of January of a certain year, we subtract the value of January of the previous year and so forth. A first reason for this choice is that it almost eliminates the effect of the annual cycle (periodicity). A second reason is that the temperature data are given in terms of “anomalies”, i.e., differences from an average which varies from month to month. By taking ν = 12 , the varying means are eliminated, and “anomalies” are effectively replaced by the actual processes (as the differences in the actual values equal the differences of “anomalies”).
We perform all analyses on both monthly and annual time scales. Figure 11 shows the differenced time series for the UAH temperature and Mauna Loa CO2 concentration at monthly scale; the symbols Δ T and Δ ln [ CO 2 ] are used interchangeably with x ˜ _ τ , 12 and y ˜ _ τ , 12 , respectively.
Comparing Figure 8 (undifferenced series) with Figure 11 (differenced series), one can verify that the latter is much more informative in terms of the directionality of the relationship of the two processes. While Figure 8 did not provide any relevant hints, Figure 11 clearly shows that, most often, the temperature curve leads and that of CO2 follows. However, there are cases where the changes in the two processes synchronize in time or even become decoupled.
Figure 12 shows the same time series at the annual time scale, with the year being defined as July–June for Δ T and February–January for Δ ln [ CO 2 ] . The reason for this differentiation will be explained below. Here, it is more evident that, most of the time, the temperature change leads and that of CO2 follows.
It is of interest here that the variability of global mean annual temperature is significantly influenced by the rhythm of ocean–atmosphere oscillations, such as ENSO, AMO, and IPO [72]. This mechanism may be a complicating factor, in turn influencing the link between temperature and CO2 concentration. This is not examined here (except a short note in the end of the section) as, given the focus in examining just the connection of the latter two processes, it lies out of our present scope.
The climacograms of the differenced time series used (actually four of the six to avoid an overcrowded graph) are shown in Figure 13. It appears that the differenced temperature time series are consistent with the condition implied by stationarity, i.e., H = 0 for the differenced process. The same does not look to be the case for the CO2 time series, particularly for the Mauna Loa time series, in which the Hurst parameter appears to be close to 1/2. Based on this, one would exclude stationarity for the Mauna Loa CO2 time series. However, a simpler interpretation of the graph is that the data record is not long enough to reveal that H = 0 for the differenced process. Actually, all available data belong to a period in which [CO2] exhibits a monotonic increasing trend (as also verified by the fact that all values of Δ ln [ CO 2 ] in Figure 11 and Figure 12 are positive, while stationarity entails a zero mean of the differenced process). Had the available database been broader, both positive and negative trends could appear. Indeed, a broader view of the [CO2] process based on palaeoclimatic data (Figure 3 and Figure 4) would justify a stationarity assumption.
The preliminary qualitative observation from graphical inspection of Figure 11 and Figure 12 suggests that the temperature change very often precedes and the CO2 change follows in the same direction. We note, though, that temperature changes alternate in sign while CO2 changes are always positive.
A quantitative analysis based on the methodology in Section 4.1 requires the study of lagged cross-correlations of the two processes. Figure 14 shows the cross-correlogram between UAH temperature and Mauna Loa CO2 concentration; the autocorrelograms of the two processes are also plotted for comparison. The fact that the cross-correlogram does not have values consistently close to zero at any of the semi-axes eliminates the possibility of an exclusive (unidirectional) causality and suggests consistency with “hen-or-egg” causality.
The maximum cross-correlation of the monthly series is 0.47 and appears at a positive lag, η 1 = 5 months, thus suggesting T → [CO2], rather than [CO2] → T, as dominant causality direction. Similar are the graphs of the other combinations of temperature and CO2 datasets, which are shown in Appendix A.5 (Figure A3, Figure A4, Figure A5, Figure A6 and Figure A7). In all cases, η 1 is positive, ranging from 5 to 11 months.
To perform similar analyses on the annual scale, we fixed the specification of a year for temperature for the period July–June, as already mentioned, and then slid the initial month specifying the beginning of a year for CO2 concentration so as to find a specification that maximizes the cross-correlation at the annual scale. In Figure 14, maximization occurs when the year specification is February–January (of the next year), i.e., if the lag is 8 months. The maximum cross-correlation is 0.66. If we keep the specification of the year for CO2 concentration the same as in temperature (July–June), then maximization occurs at lag one year (12 months) and the maximum cross-correlation is 0.52. Table 1 summarizes the results for all combinations examined. The lags are always positive. They vary between 8 and 14 months for a sliding window specification and are 12 months for the fixed window specification. Most interestingly, the opposite phase in the annual cycle of CO2 concentration in the South Pole, with respect to the other three sites, does not produce any noteworthy differences in the shape of the cross-correlogram nor in the time lags maximizing the cross-correlations.
While, as explained in Section 4.2 and Section 5.1, the Granger test has weaknesses that may not help in drawing conclusions, for completeness and as a confirmation, we list its results here:
  • For the monthly scale and the causality direction [CO2] → T, the null hypothesis is not rejected at all usual significance levels for lag η = 1 and is rejected for significance level 1% for η = 2–8, with minimum attained p-value 1.8 × 10−4 for η = 6.
  • For the monthly scale and the causality direction T → [CO2], the null hypothesis is rejected at all usual significance levels for all lags η, with minimum attained p-value 2.1 × 10−8 for η = 7.
  • For the monthly scale, the attained p-values in the direction T → [CO2] are always smaller than in direction [CO2] → T by about 4 to 5 orders of magnitude, thus clearly supporting T → [CO2] as dominant direction.
  • For the annual scale with fixed year specification and the causality direction [CO2] → T, the null hypothesis is not rejected at all usual significance levels for any lag η, thus indicating that this causality direction does not exist.
  • For the annual scale with fixed year specification and the causality direction T → [CO2], the null hypothesis is not rejected at significance level 1% for all lags η = 1–6, with minimum attained p-value 5% for lag η = 2, thus supporting this causality direction at this significance level.
  • For the annual scale with fixed year specification, the attained p-values in the direction T → [CO2] are always smaller than in direction [CO2] → T, again clearly supporting T → [CO2] as the dominant direction.
We note that the test cannot be applied for the sliding time window case and, hence, we cannot provide results for this case.
We add a final remark, in view of a comment by Masters and Benestad [73] on the already mentioned study by Humlum et al. [45], in which they claim that “the inter-annual fluctuations in atmospheric CO2 produced by ENSO can lead to a misdiagnosis of the long-term cause of the recent atmospheric CO2 increase”. Inspired by this comment, we have made a preliminary three-variable investigation using differenced temperatures (UAH), logarithms of CO2 concentrations (Mauna Loa), and Equatorial South Oscillation index (SOI) characterizing ENSO. The investigation has been made on a monthly scale. Δ ln [ CO 2 ] has been linearly regressed with Δ T and the running average of SOI for the previous 12 months. At synchrony (without applying any time lag), the correlation of SOI with Δ T is 0.40, higher than that of Δ T and Δ ln [ CO 2 ] (0.24, as seen in Figure 14 at lag 0). The highest determination coefficient for the three regressed quantities is obtained when the time lag between Δ ln [ CO 2 ] and Δ T is again 5 months, as in the two-variable case (the optimal lag for SOI is 0, but the regression is virtually insensitive to the change of that lag). Its value is r 2 = 0.23 , corresponding to r = 0.48 , i.e., only slightly higher than the maximum cross-correlation coefficient of the two variable-case (which is 0.47 as seen in Table 1). In other words, by including ENSO in the modelling framework, the results do not change.
In brief, all above confirm the results of our methodology that the dominant direction of causality is T → [CO2].

6. Physical Interpretation

The omnipresence of positive lags on both monthly and annual time scales and the confirmation by Granger tests reduce the likelihood that our results are statistical artefacts. Still, our results require physical interpretation which we seek in the natural process of soil respiration.
Soil respiration, Rs, defined to be the flux of microbially and plant-respired CO2, clearly increases with temperature. It is known to have increased in the recent years [74,75]. Observational data of Rs (e.g., [76,77]; see also [78]) show that the process intensity increases with temperature. Rate of chemical reactions, metabolic rate, as well as microorganism activity, generally increase with temperature. This has been known for more than 70 years (Pomeroy and Bowlus [79]) and is routinely used in engineering design.
The Figure 6.1 of the latest report of the IPCC [75] provides a quantification of the mass balance of the carbon cycle in the atmosphere that is representative of recent years. The soil respiration, assumed to be the sum of respiration (plants) and decay (microbes), is 113.7 Gt C/year (IPCC gives a value of 118.7 including fire, which along with biomass burning, is estimated to be 5 Gt C/year by Green and Byrne [80]).
We can expect that sea respiration would also have increased. Moreover, outgassing from the oceans must also have increased as the solubility of CO2 in water decreases with increasing temperature [14,81]. In addition, photosynthesis must have increased, as in the 21st century the Earth has been greening, mostly due to CO2 fertilization effects [82] and human land-use management [83]. Specifically, satellite data show a net increase in leaf area of 2.3% per decade [83]. The sums of carbon outflows from the atmosphere (terrestrial and maritime photosynthesis as well as maritime absorption) amount to 203 Gt C/year. The carbon inflows to the atmosphere amount to 207.4 Gt C/year and include natural terrestrial processes (respiration, decay, fire, freshwater outgassing as well as volcanism and weathering), natural maritime processes (respiration) as well as anthropogenic processes. The latter comprise human CO2 emissions related to fossil fuels and cement production as well as land-use change, and amount to 7.7 and 1.1 Gt C/year, respectively. The change in carbon fluxes due to natural processes is likely to exceed the change due to anthropogenic CO2 emissions, even though the latter are generally regarded as responsible for the imbalance of carbon in the atmosphere.

7. Conclusions

The relationship between atmospheric concentration of carbon dioxide and the global temperature is widely recognized, and it is common knowledge that increasing CO2 concentration plays a major role in enhancement of the greenhouse effect and contributes to global warming.
While the fact that these two variables are tightly connected is beyond doubt, the direction of the causal relationship needs to be studied further. The purpose of this study is to complement the conventional and established theory, that increased CO2 concentration due to anthropogenic emissions causes an increase of temperature, by considering the concept of reverse causality. The problem is obviously more complex than that of exclusive roles of cause and effect, qualifying as a “hen-or-egg” (“ὄρνις ἢ ᾠὸν”) causality problem, where it is not always clear which of two interrelated events is the cause and which the effect. Increased temperature causes an increase in CO2 concentration and, hence, we propose the formulation of the entire process in terms of a “hen-or-egg” causality.
We examine the relationship of global temperature and atmospheric carbon dioxide concentration using the most reliable global data that are available—the data gathered from several sources, covering the common time interval 1980–2019, available at the monthly time step.
The results of the study support the hypothesis that both causality directions exist, with T → CO2 being the dominant, despite the fact that CO2T prevails in public, as well as in scientific, perception. Indeed, our results show that changes in CO2 follow changes in T by about six months on a monthly scale, or about one year on an annual scale.
The opposite causality direction opens a nurturing interpretation question. We attempted to interpret this mechanism by noting that the increase in soil respiration, reflecting the fact that the intensity of biochemical process increases with temperature, leads to increasing natural CO2 emission. Thus, the synchrony of rising temperature and CO2 creates a positive feedback loop. This poses challenging scientific questions of interpretation and modelling for further studies. In this respect, we welcome the review by Connolly [14], which already proposes interesting interpretations within a wider epistemological framework and in connection with a recent study [84]. In our opinion, scientists of the 21st century should have been familiar with unanswered scientific questions as well as with the idea that complex systems resist simplistic explanations.

Author Contributions

Conceptualization, D.K.; methodology, D.K.; software: D.K.; validation, Z.W.K.; formal analysis, D.K.; investigation, D.K. and Z.W.K.; data curation, D.K.; writing—original draft preparation, D.K. and Z.W.K.; writing—review and editing, D.K. and Z.W.K.; visualization, D.K. and Z.W.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding but was motivated by the scientific curiosity of the authors.

Acknowledgments

Some negative comments of two anonymous reviewers for an earlier submission of this manuscript in another journal (which we have posted online: manuscript at http://0-dx-doi-org.brum.beds.ac.uk/10.13140/RG.2.2.29154.15045/1, reviews at http://0-dx-doi-org.brum.beds.ac.uk/10.13140/RG.2.2.14524.87681) helped us improve the presentation and strengthen our arguments against their comments. We appreciate these reviewers’ suggestions of relevant published works, which we were unaware of. We are grateful to the reviewers of the submission of our study to Sci, Yog Aryal, Ronan Connolly, and Stavros Alexandris, whose constructive comments helped us to improve the paper, mostly by adding a lot of additional information and clarification in appendices. The discussion of DK with Antonis Christofides helped to substantially improve an initial version of Appendix A.4, which he informally reviewed. We thank the journal Sci for the interesting experience it offered us.

Conflicts of Interest

The authors declare no conflict of interest.

Data Availability

The two temperature time series and the Mauna Loa CO2 time series are readily available on monthly scale from http://climexp.knmi.nl. All NOAA CO2 data are available from https://www.esrl.noaa.gov/gmd/ccgg/trends/gl_trend.html. The CO2 data of Mauna Loa were retrieved from http://climexp.knmi.nl/data/imaunaloa_f.dat while the original measurements are in https://www.esrl.noaa.gov/gmd/dv/iadv/graph.php?code=MLO. The Barrow series is available (in irregular step) in https://www.esrl.noaa.gov/gmd/dv/iadv/graph.php?code=BRW, and the South Pole series in https://www.esrl.noaa.gov/gmd/dv/data/index.php?site=SPO. All these data were accessed (using the “Download data” link in the above sites) in June 2020. The global CO2 series is accessed at https://www.esrl.noaa.gov/gmd/ccgg/trends/gl_data.html, of which the “Globally averaged marine surface monthly mean data” are used here. The palaeoclimatic data of Vostok CO2 were retrieved from http://cdiac.ess-dive.lbl.gov/ftp/trends/co2/vostok.icecore.co2 (dated January 2003, accessed September 2018) and the temperature data from http://cdiac.ess-dive.lbl.gov/ftp/trends/temp/vostok/vostok.1999.temp.dat (dated January 2000, accessed September 2018).

Appendix A

Appendix A.1. On Early Non-Systematic Measurements of CO2

This Appendix (not contained in Version 1 of our paper) addresses comments by all three reviewers of Version 1, Yog Aryal [85], Ronan Connolly [14], and Stavros Alexandris [86], about the reasons why we delimit our analysis to the period 1980–2019. The two latter reviewers suggested using earlier data compiled by Beck (2007), who referred to old chemical analyses of atmospheric concentration of CO2.
We are sympathetic to the passion of the late Ernst-Georg Beck who, being a biology teacher, sacrificed a lot of time and effort to the exciting exercise of digging out old CO2 measurements. Indeed, it could be worthwhile to have a critical look at the historical data and to try to make order in them and utilize them. However, this would certainly warrant an individual paper with this particular aim.
Historically, it was not the first review paper of this sort. For instance, in his Table 1, Beck [87] refers to old works by Letts and Blake (~1900; [88]), who considered 252 papers with data (all in 19th century), and to Stepanova [89], who considered 229 papers with data (130 in 19th century and 99 in 20th century). Beck himself [87] considered 156 papers with data (82 in 19th and 74 in 20th century).
As usual, it is instructive to consider the paper by Beck [87] jointly with critical commentaries published later in the journal where the original paper appeared [90,91]. In particular, R.F. Keeling [90] opined that the old chemical measurements examined by Beck [87] “exhibit far too much geographic and short-term temporal variability to plausibly be representative of the background. The variability of these early measurements must therefore be attributed to ‘local or regional’ factors or poor measurement practice”. Keeling [90] also noted “basic accounting problems”. “Beck’s 11-year averages show large swings, including an increase from 310 to 420 ppm between 1920 and 1945 (Beck’s Figure 11)”. “To drive an increase of this magnitude globally requires the release of 233 billion metric tons of carbon to the atmosphere. The amount is equivalent to more than a third of all the carbon contained in land plants globally. […] To make a credible case, Beck needed to offer evidence for losses or gains of carbon of this magnitude from somewhere. He offered none.”
Meijer [91] expressed the opinion that Beck’s work “contains major flaws, such that the conclusions are wrong”. He also wrote: “The measurements presented in the paper are indeed useless for the purpose the author wants to use them, certainly in the way the author interprets them”. There is a lack of interpretation of diurnal and seasonal variability (effects called the “diurnal” and the “seasonal” rectifier in the literature) and consideration of atmospheric mixing or lack thereof. Meijer also criticized the lack of meta-data: “The necessary data to judge, namely measurement height, consecutive length of a record and especially temporal resolution, are lacking in [Beck’s] Table 2. In the light of the above, the whole ‘Discussion and Conclusion’ section is invalid, including [Beck’s] Figure 11, Figure 12, Figure 13 and Figure 14”. Indeed, the records mentioned in Beck’s Table 2 were local and short-lasting, with the longest periods being 1920–1926. Beck’s Figure 11 and Figure 13 show concatenated short segments of data at different places.
There are some other puzzling elements in Beck’s paper. For instance, in his Figure 5, referring to data from a meteorological station near Giessen, the variability of high amplitude seems suspicious and not physically realistic. In particular, from June to August 1940, the measured CO2 concentration increases from 340 to 550 ppm (much more than in Beck´s Figure 5 discussed by Keeling [90] and Meijer [91] as quoted above), with weird seasonal behaviour. Beck himself admitted that the results for Giessen “need to be adjusted downwards to take account of anthropogenic sources of CO2 from nearby city, an influence that has been estimated as lying between 10 and 70 ppm […] by different authors”.
The controversy and disputes among these authors extended beyond pure scientific issues. Thus, Beck [87] wrote “[t]he data accepted […] had to be sufficiently low to be consistent with the greenhouse hypothesis of climate change controlled by rising CO2 emissions from fossil fuel burning”. On the other hand, Meijer [91] wrote: “The author even accuses the pioneers Callendar and [Charles David] Keeling of selective data use, errors or even something close to data manipulation”. In addition, [R.F.] Keeling [90] noted: “Beck is […] wrong when he asserts that the earlier data have been discredited only because they don’t fit a preconceived hypothesis of CO2 and climate. […] Instead, the data have been ignored because they cannot be accepted as representative without violating our understanding of how fast the atmosphere mixes”.
In view of the above questions about the data reliability as well as the controversies and disputes, we decided to limit the period of our study in 1980–2019 in which the measurements are systematic and verifiable because they are made in several locations simultaneously.

Appendix A.2. Some Notes on the Averaged Differenced Process

The cumulative process of the differenced process x ˜ _ τ , ν will be
X ˜ _ κ , ν x ˜ _ 1 , ν + x ˜ _ 2 , ν + + x ˜ _ κ , ν = x _ 1 + ν x _ 1 + x _ 2 + η x _ 2 + + x _ κ + ν x _ κ = Χ _ κ + ν Χ _ ν Χ _ κ .
Note that for η = 1, this simplifies to
X ˜ _ κ , 1 = Χ _ κ + 1 Χ _ 1 Χ _ κ = x _ κ + 1 x _ 1 = x ˜ _ κ , 1 x ˜ _ κ .
Following Equation (7), the average differenced process at discrete time scale κ = η will be
x ˜ _ τ ( κ ) = X ˜ _ τ κ , κ X ˜ _ ( τ 1 ) κ , κ κ = ( Χ _ τ κ + κ Χ _ κ Χ _ τ κ ) ( Χ _ ( τ 1 ) κ + κ Χ _ κ Χ _ ( τ 1 ) κ ) κ
which, noting that in the rightmost part, the two terms Χ _ κ cancel each other and by virtue of (7), simplifies to
x ˜ _ τ ( κ ) = x _ τ + 1 ( κ ) x _ τ ( κ ) = x ˜ _ τ , 1 ( κ ) .
In other words, the average differenced process equals the differenced average process in the case that the differencing time step η has been chosen as equal to the averaging time scale κ. For κ = η = 1, we have x ˜ _ τ ( 1 ) x ˜ _ τ , 1 x ˜ _ τ .

Appendix A.3. Some Notes on Time Directionality of Causal Systems

In a unidirectional causal system in continuous time t, in which the process x _ ( t ) is the cause of y _ ( t ) , an equation of the form
y _ ( t ) = 0 α ( s ) x _ ( t s ) d s
should hold [67], where α ( t ) is the impulse response function. The causality condition is thus
α ( t ) = 0   for   t < 0 .
Here, we consider systems with positive dependence, in which α ( t ) 0 for t 0 , which are possibly also excited by another process v _ ( t ) , independent of x _ ( t ) . Working in discrete time, we write
y _ τ = j = 0 α j x _ τ j + v _ τ .
Assuming (without loss of generality) zero means for all processes, multiplying by x _ τ η , taking expected values, and denoting the cross-covariance function as c x y [ η ] E [ x _ τ η y _ τ ] and the autocovariance function as c x [ η ] E [ x _ τ η x _ τ ] , we find
c x y [ η ] = j = 0 α j c x [ η j ] .
For η > 0 , using the property that c x [ η ] is an even function ( c x [ η ] = c x [ η ] ), we get
c x y [ η ] = j = 0 α j c x [ j η ] = j = 0 η 1 α j c x [ η j ] + j = η α j c x [ j η ] ,
and for the negative part
c x y [ η ] = j = 0 α j c x [ j + η ] .
With intuitive reasoning, assuming that the autocovariance function is decreasing ( c x [ j ] < c x [ j ] for j > j ), as usually happens in natural processes, we may see that the rightmost term of Equations (A9) and (A10) should be decreasing functions of η (as for j > j it will be c x [ j η ] < c x [ j η ] and c x [ j + η ] < c x [ j + η ] ). However, the term j = 0 η 1 α j c x [ η j ] of Equation (A9) is not decreasing. Therefore, it should attain a maximum value at some positive lag η = η 1 . Thus, a positive maximizing lag, η = η 1 > 0 , is a necessary condition for causality direction from x _ τ   to   y _ τ . Conversely, the condition that the maximizing lag is negative is a sufficient condition to exclude the causality direction exclusively from x _ τ   to   y _ τ .
All above arguments remain valid if we standardize (divide) by the product of standard deviations of the processes x _ τ and y _ τ and, thus, we can replace cross-covariances c x y [ η ] with cross-correlations r x y [ η ] (or, in the case of differenced processes, r x ˜ y ˜ [ ν , η ] ).

Appendix A.4. Some Notes on the Alternative Procedures on Causality

Reviewer Yog Aryal [85] opined that we missed referring to the recent relevant works by Hannart et al. [92] and Verbitsky et al. [93]. In response to this comment, we include this Appendix (not contained in Version 1 of our paper) explaining, in brief, why we do not compare our results with the ones of those studies, also noting that only the latter study contains material that is prima facie comparable to ours. The former study, focusing on the so-called causal counterfactual theory, is more theoretical and also much more interesting. While we, too, are preparing a theoretical study, in which we will discuss some theories in detail, in this Appendix, we give some key elements of our theoretical disagreements and a counterexample that illustrates the disagreements.
We first note that in order to define causality, Hannart et al. [92] refer to the work on the 18th century philosopher David Hume and, in particular, his famous book Enquiry concerning Human Understanding [94] first published in 1748. From this book, we wish to quote the following important passage, which emphasizes the difficulties even in defining causality:
Our thoughts and enquiries are, therefore, every moment, employed about this relation: Yet so imperfect are the ideas which we form concerning it, that it is impossible to give any just definition of cause, except what is drawn from something extraneous and foreign to it.
Hannart et al. [92], while studying the probability of occurrence of an event Y, introduced the two-valued variable X f to indicate whether or not a forcing f is present, and continue as follows:
The probability p 1 = P ( Y = 1 | X f = 1 ) of the event occurring in the real world, with f present, is referred to as factual, while p 0 = P ( Y = 1 | X f = 0 ) is referred to as counterfactual. Both terms will become clear in the light of what immediately follows. The so-called fraction of attributable risk (FAR) is then defined as
FAR = 1 p 0 p 1
The FAR is interpreted as the fraction of the likelihood of an event that is attributable to the external forcing.
They also show that under some conditions, FAR is a probability which they denote PN and call probability of necessary causality. They stress that it “is important to distinguish between necessary and sufficient causality” and they associate PN (or FAR) “with the first facet of causality, that of necessity”. They claim to have “introduced its second facet, that of sufficiency, which is associated with the symmetric quantity 1 ( 1 p 1 ) / ( 1 p 0 ) ”; they denote it as PS, standing for probability of sufficient causality.
Central to the logical framework of Hannart et al. [92] is the notion of intervention of an experimenter, which is equivalent to experimentation with the ability to set the value of the assumed cause to a desired value. Clearly, this is feasible in laboratory experiments and infeasible in natural processes. The authors resort to the “so-called in silico experimentation” which, despite the impressive name chosen, is intervention in a mathematical model that represents the process. Hence, objectively, they examine the “causality” that is embedded in the model rather than the natural causality. One may argue that this it totally unnecessary. It would be better to inspect the model’s equations or code to investigate what causality has been embedded in the model instead of running simulations and calculating probabilities. In particular, if the models used are climate models as in [92], their inability to effectively describe (perform in “prime time”) the real-world processes [50,95,96,97,98,99,100] makes the entire endeavour futile. Another notion these authors use is exogeneity, which is related to the so-called causal graph, reflecting the assumed dependencies among the studied variables. Specifically, they state “a sufficient condition for X to be exogenous wrt any variable is to be a top node of a causal graph”.
Here, we will use the simple example of Section 4.2, temperature–clothes weight–sweat, to show that using the quantities FAR (or PN) and PS may give spurious results that do not correspond to necessary or sufficient conditions for causality, at least with their meaning in our paper.
We use the two-valued random variables x _ , y _ , z _ to model the states of temperature, clothes weight, and sweat, respectively. We designate the following states:
  • x = 1: being hot above a threshold;
  • y = 1: wearing clothes with weight below a threshold;
  • z = 1: sweat quantity above a threshold;
and the opposite states with x = 0 , y = 0 , z = 0 , respectively. We choose the threshold of temperature so that P { x _ = 0 } = P { x _ = 1 } = 0.5 and that of clothes weight so that P { y _ = 0 } = P { y _ = 1 } = 0.5 . We choose a small probability, 0.05, of wearing light clothes when cold, or heavy clothes when hot, i.e., P { y _ = 1 | x _ = 0 } = P { y _ = 0 | x _ = 1 } = 0.05 (generally, we avoid choosing zero probabilities; rather the minimum value we choose is 0.05).
Using the definition of conditional probability,
P { y _ = y | x _ = x } = P { y _ = y , x _ = x } P { x _ = x } ,
we find the probability matrix A with elements a i j = P { x _ = i , y _ = j } as follows:
A = [ 0.475 0.025 0.025 0.475 ] x = 0 x = 1 y = 0 y = 1 .
Now, we assign plausible values to the conditional probabilities of high sweat, P { z _ = 1 | x _ = x , y _ = y } , as follows:
  • Cold, heavy clothes:     P { z _ = 1 | x _ = 0 , y _ = 0 } = 0.2
  • Cold, light clothes:     P { z _ = 1 | x _ = 0 , y _ = 1 } = 0.1
  • Hot, heavy clothes:     P { z _ = 1 | x _ = 1 , y _ = 0 } = 0.95
  • Hot, light clothes:     P { z _ = 1 | x _ = 1 , y _ = 1 } = 0.80
Again, we have avoided setting any of the conditional probabilities to 0 (or 1), and we have used multiples of 0.05 for all of them.
Using the definition of conditional probability in the form
P { z _ = z | x _ = x , y _ = y } = P { z _ = z , y _ = y , x _ = x } P { y _ = y , x _ = x } ,
we find the joint probabilities for each of the triplets { x , y , z } that are shown in Table A1.
Table A1. Joint probabilities P { x _ = x , y _ = y , z _ = z } for all triplets { x , y , z }
Table A1. Joint probabilities P { x _ = x , y _ = y , z _ = z } for all triplets { x , y , z }
xyz = 0z = 1
000.380.095
010.02250.0025
100.001250.02375
110.0950.38
P { z _ = z } = 0.498750.50125
Now, assume that we let an “artificial intelligence entity” (AIE) decide on causality based on the probability rules of the Hannart et al. [92] framework. Our AIE has access to numerous videos of people and is “trained” to assign accurate values of y and z, referring to clothes and sweat, based on the images in videos. In the video images, no thermometers are shown and, thus, our AIE cannot assign values of x, nor can it be aware of the notion of temperature. Our AIE tries to construct a causal graph putting, say, y _ as a top node and z _ as an end node; hence, it assumes that y _ is exogenous. Based on the huge information it can access, our AIE can (a) claim that it has constructed a prediction model based on one part of the data (e.g., using the so-called deep-learning technique) and, hence, is able to perform “in silico experimentation” (even though this is not absolutely necessary) and (b) accurately estimate the joint and conditional probabilities related to { y , z } using either the model, the data, or both. Provided that the dataset is large enough, it will come up with the true values for the conditional probabilities, which are b i j = P { y _ = i , z _ = j } and c i j = P { z _ = j | y _ = i } , and form the matrices B and C, respectively, with values as follows:
B = [ 0.38125 0.11875 0.1175 0.3825 ] y = 0 y = 1 z = 0 z = 1 ,     C = [ 0.7625 0.2375 0.235 0.765 ] y = 0 y = 1 z = 0 z = 1 .
Here, the true values b i j have been determined from the values of Table A1 noting that
b i j = P { y _ = i , z _ = j } = P { z _ = j , y _ = i , x _ = 0 } + P { z _ = j , y _ = i , x _ = 1 }
and the true values c i j have been determined from the definition of conditional probability:
P { z _ = z | y _ = y } = P { z _ = z , y _ = y } P { y _ = y } .
Our AIE will then implement the causality conditions of sweat on clothes weight, assigning p 0 = P { z _ = 1 | y _ = 0 } = 0.2375 and p 1 = P { z _ = 1 | y _ = 1 } = 0.765 . It will further calculate the probability of necessary causality as PN = 0.690, and the probability of sufficient causality even higher, PS = 0.692. Hence, our AIE will inform us that there is all necessary and sufficient evidence that light clothes cause high sweat.
Now, coming to the study by Verbitsky et al. [93], we notice that it assumes that “each time series is a variable produced by its hypothetical low dimensional system of dynamical equations” and uses the technique of distances of multivariate vectors for reconstructing the system dynamics. As demonstrated in Koutsoyiannis [101], such assumptions and techniques are good for simple toy models but, when real-world systems are examined, low dimensionality appears as a statistical artifact because the reconstruction actually needs an incredibly high number of observations to work, which are hardly available. The fact that the sums of multivariate vectors of distances is a statistical estimator with huge uncertainty is often missed in studies of this type, which treat data as deterministic quantities to obtain unreliable results. We do not believe that the Earth system and Earth processes (including global temperature and CO2) are of low dimensionality, and we deem it unnecessary to discuss the issue further. We only note the fact that global temperature and CO2 virtually behave as Gaussian, which enables reliable estimation of standard correlations and dismiss the need to use the overly complex and uncertain correlation sums.

Appendix A.5. Additional Graphical Depictions

Figure A1. Normal probability plots of ΔT and Δln[CO2] where T is the UAH temperature and [CO2] is the CO2 concentration at Mauna Loa at monthly scale.
Figure A1. Normal probability plots of ΔT and Δln[CO2] where T is the UAH temperature and [CO2] is the CO2 concentration at Mauna Loa at monthly scale.
Sci 02 00083 g0a1
Figure A2. Scatter plot of Δ T and Δ ln [ CO 2 ] where T is the UAH temperature and [CO2] is the CO2 concentration at Mauna Loa at monthly scale; the two quantities are lagged in time using the optimal lag of 5 months (Table 1). The two linear regression lines are also shown in the figure.
Figure A2. Scatter plot of Δ T and Δ ln [ CO 2 ] where T is the UAH temperature and [CO2] is the CO2 concentration at Mauna Loa at monthly scale; the two quantities are lagged in time using the optimal lag of 5 months (Table 1). The two linear regression lines are also shown in the figure.
Sci 02 00083 g0a2
Figure A3. Auto- and cross-correlograms of the differenced time series of UAH temperature and Barrow CO2 concentration.
Figure A3. Auto- and cross-correlograms of the differenced time series of UAH temperature and Barrow CO2 concentration.
Sci 02 00083 g0a3
Figure A4. Auto- and cross-correlograms of the differenced time series of UAH temperature and South Pole CO2 concentration.
Figure A4. Auto- and cross-correlograms of the differenced time series of UAH temperature and South Pole CO2 concentration.
Sci 02 00083 g0a4
Figure A5. Auto- and cross-correlograms of the differenced time series of UAH temperature and global CO2 concentration.
Figure A5. Auto- and cross-correlograms of the differenced time series of UAH temperature and global CO2 concentration.
Sci 02 00083 g0a5
Figure A6. Auto- and cross-correlograms of the differenced time series of CRUTEM4 temperature and Mauna Loa CO2 concentration.
Figure A6. Auto- and cross-correlograms of the differenced time series of CRUTEM4 temperature and Mauna Loa CO2 concentration.
Sci 02 00083 g0a6
Figure A7. Auto- and cross-correlograms of the differenced time series of CRUTEM4 temperature and global CO2 concentration.
Figure A7. Auto- and cross-correlograms of the differenced time series of CRUTEM4 temperature and global CO2 concentration.
Sci 02 00083 g0a7

References

  1. Πλούταρχος, Συμποσιακά Β’ (Plutarch, Quaestiones Convivales B)—Βικιθήκη. Available online: https://el.wikisource.org/wiki/Συμποσιακά_Β’ (accessed on 15 October 2020).
  2. Plutarchus-Quaestiones Convivales (612c–748d). Available online: http://www.poesialatina.it/_ns/Greek/html/map-en.html (accessed on 15 October 2020).
  3. IEA (International Energy Agency). Global Energy Review 2020; IEA: Paris, French, 2020; Available online: https://www.iea.org/reports/global-energy-review-2020 (accessed on 1 July 2020).
  4. Le Quéré, C.; Jackson, R.B.; Jones, M.W.; Smith, A.J.; Abernethy, S.; Andrew, R.M.; De-Gol, A.J.; Willis, D.R.; Shan, Y.; Canadell, J.G.; et al. Temporary reduction in daily global CO2 emissions during the COVID-19 forced confinement. Nat. Clim. Chang. 2020, 10, 647–653. [Google Scholar] [CrossRef]
  5. Global Carbon Atlas. Available online: http://www.globalcarbonatlas.org/en/CO2-emissions (accessed on 1 July 2020).
  6. Our World in Data. Available online: https://ourworldindata.org/grapher/annual-co-emissions-by-region (accessed on 1 July 2020).
  7. Peters, G.P.; Marland, G.; Le Quéré, C.; Boden, T.; Canadell, J.G.; Raupach, M.R. Rapid growth in CO2 emissions after the 2008–2009 global financial crisis. Nat. Clim. Chang. 2012, 2, 2–4. [Google Scholar] [CrossRef]
  8. Arrhenius, S. On the influence of carbonic acid in the air upon the temperature of the ground. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1896, 41, 237–276. [Google Scholar] [CrossRef] [Green Version]
  9. De Marchi, L. Le Cause dell’Era Glaciale; Premiato dal R. Instituto Lombardo: Pavia, Italy, 1895. [Google Scholar]
  10. Tyndall, J. Heat a Mode of Motion, 2nd ed.; Longmans: London, UK, 1865; Available online: https://archive.org/details/heatamodemotion05tyndgoog/ (accessed on 1 September 2020).
  11. Foote, E. Circumstances affecting the heat of the sun’s rays. Am. J. Sci. Arts 1856, 22, 382–383. [Google Scholar]
  12. Jackson, R. Eunice Foote, John Tyndall and a question of priority. Notes Rec. 2020, 74, 105–118. [Google Scholar] [CrossRef]
  13. Ortiz, J.D.; Jackson, R. Understanding Eunice Foote’s 1856 experiments: Heat absorption by atmospheric gases. Notes Rec. 2020. [Google Scholar] [CrossRef]
  14. Connolly, R. Review of “Atmospheric temperature and CO2: Hen-or-egg causality?” (Version 1). Sci 2020. Available online: https://0-www-mdpi-com.brum.beds.ac.uk/2413-4155/2/4/77/review_report (accessed on 9 October 2020).
  15. Schmidt, G.A.; Ruedy, R.A.; Miller, R.L.; Lacis, A.A. Attribution of the present-day total greenhouse effect. J. Geophys. Res. 2010, 115, D20106. [Google Scholar] [CrossRef]
  16. Shaviv, N.J.; Veizer, J. Celestial driver of Phanerozoic climate? Geol. Soc. Am. Today 2003, 13, 4–10. [Google Scholar] [CrossRef]
  17. Shaviv, N.J. On climate response to changes in the cosmic ray flux and radiative budget. J. Geophys. Res. Space Phys. 2005, 110, A08105. [Google Scholar] [CrossRef] [Green Version]
  18. Rahmstorf, S.; Archer, D.; Ebel, D.S.; Eugster, O.; Jouzel, J.; Maraun, D.; Neu, U.; Schmidt, G.A.; Sever-Inghaus, J.; Weaver, A.J.; et al. Cosmic rays, carbon dioxide, and climate. Eos 2004, 85, 38–41. [Google Scholar] [CrossRef] [Green Version]
  19. Shaviv, N.J.; Veizer, J. Comment on Cosmic rays, carbon dioxide, and climate. Eos 2004, 85, 510. [Google Scholar] [CrossRef]
  20. Rahmstorf, S.; Archer, D.; Ebel, D.S.; Eugster, O.; Jouzel, J.; Maraun, D.; Neu, U.; Schmidt, G.A.; Sever-Inghaus, J.; Weaver, A.J.; et al. Reply to Comment [by Shaviv and Veizer]. Eos 2004, 85, 511. [Google Scholar] [CrossRef]
  21. Royer, D.L.; Berner, R.A.; Montañez, I.P.; Tabor, N.J.; Beerling, D.J. CO2 as a primary driver of phanerozoic climate. Geol. Soc. Am. Today 2004, 14, 4–10. [Google Scholar] [CrossRef]
  22. Shaviv, N.J.; Veizer, J. Comment on “CO2 as a primary driver of Phanerozoic climate” by Royer et al. Geol. Soc. Am. Today 2004, 14, 18. [Google Scholar]
  23. Royer, D.L.; Berner, R.A.; Montañez, I.P.; Tabor, N.J.; Beerling, D.J. Reply to Comment on “CO2 as a primary driver of Phanerozoic climate” by Royer et al. Geol. Soc. Am. Today 2004, 14, 18. [Google Scholar]
  24. Roe, G. In defense of Milankovitch. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef] [Green Version]
  25. Jouzel, J.; Lorius, C.; Petit, J.R.; Genthon, C.; Barkov, N.I.; Kotlyakov, V.M.; Petrov, V.M. Vostok ice core: A continuous isotope temperature record over the last climatic cycle (160,000 years). Nature 1987, 329, 403–408. [Google Scholar] [CrossRef]
  26. Petit, J.R.; Jouzel, J.; Raynaud, D.; Barkov, N.I.; Barnola, J.-M.; Basile, I.; Bender, M.; Chappellaz, J.; Davis, M.; Delayque, G.; et al. Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature 1999, 399, 429–436. [Google Scholar] [CrossRef] [Green Version]
  27. Caillon, N.; Severinghaus, J.P.; Jouzel, J.; Barnola, J.M.; Kang, J.; Lipenkov, V.Y. Timing of atmospheric CO2 and Antarctic temperature changes across Termination III. Science 2003, 299, 1728–1731. [Google Scholar] [CrossRef] [Green Version]
  28. Soon, W. Implications of the secondary role of carbon dioxide and methane forcing in climate change: Past, present, and future. Phys. Geogr. 2007, 28, 97–125. [Google Scholar] [CrossRef]
  29. Pedro, J.B.; Rasmussen, S.O.; van Ommen, T.D. Tightened constraints on the time-lag between Antarctic temperature and CO2 during the last deglaciation. Clim. Past 2012, 8, 1213–1221. [Google Scholar] [CrossRef] [Green Version]
  30. Koutsoyiannis, D. Time’s arrow in stochastic characterization and simulation of atmospheric and hydrological processes. Hydrol. Sci. J. 2019, 64, 1013–1037. [Google Scholar] [CrossRef]
  31. Beeman, J.C.; Gest, L.; Parrenin, F.; Raynaud, D.; Fudge, T.J.; Buizert, C.; Brook, E.J. Antarctic temperature and CO2: Near-synchrony yet variable phasing during the last deglaciation. Clim. Past 2019, 15, 913–926. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, Z.; Huang, S.; Jin, Z. Breakpoint lead-lag analysis of the last deglacial climate change and atmospheric CO2 concentration on global and hemispheric scales. Quat. Int. 2018, 490, 50–59. [Google Scholar] [CrossRef]
  33. Scotese, C.R. Phanerozoic Temperatures: Tropical Mean Annual Temperature (TMAT), Polar Mean Annual Temperature (PMAT), and Global Mean Annual Temperature (GMAT) for the last 540 Million Years. In Proceedings of the Earth’s Temperature History Research Workshop, Smithsonian National Museum of Natural History, Washington, DC, USA, 30–31 March 2018; Available online: https://www.researchgate.net/publication/324017003 (accessed on 8 September 2020).
  34. Grossman, E.L. Oxygen isotope stratigraphy. In The Geological Time Scale; Gradstein, F.M., Ogg, J.G., Schmitz, M.D., Ogg, G.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2012; pp. 181–206. [Google Scholar]
  35. Veizer, J.; Prokoph, A. Temperatures and oxygen isotopic composition of Phanerozoic oceans. Earth Sci. Rev. 2015, 146, 92–104. [Google Scholar] [CrossRef]
  36. Boucot, A.J.; Xu, C.; Scotese, C.; Morley, R.J. Phanerozoic Paleoclimate: An Atlas of Lithologic Indications of Climate; SEPM Concepts Sedimentol. Paleontol. No 11; SEPM Publication Editors: Tulsa, OK, USA, 2013; p. 30. ISBN 978-1-56576-282-4. [Google Scholar]
  37. Davis, W.J. The relationship between atmospheric carbon dioxide concentration and global temperature for the last 425 million years. Climate 2017, 5, 76. [Google Scholar] [CrossRef] [Green Version]
  38. Berner, R.A. Addendum to “Inclusion of the weathering of volcanic rocks in the GEOCARBSULF model” (R. A. Berner, 2006, v. 306, p. 295–302). Am. J. Sci. 2008, 308, 100–103. [Google Scholar] [CrossRef]
  39. Ekart, D.D.; Cerling, T.E.; Montanez, I.P.; Tabor, N.J. A 400 million year carbon isotope record of pedogenic carbonate: Implications for paleoatmospheric carbon dioxide. Am. J. Sci. 1999, 299, 805–827. [Google Scholar] [CrossRef]
  40. International Commission on Stratigraphy. Available online: https://stratigraphy.org/chart (accessed on 1 September 2020).
  41. Kodra, E.; Chatterjee, S.; Ganguly, A.R. Exploring Granger causality between global average observed time series of carbon dioxide and temperature. Theor. Appl. Climatol. 2011, 104, 325–335. [Google Scholar] [CrossRef] [Green Version]
  42. Stern, D.I.; Kaufmann, R.K. Anthropogenic and natural causes of climate change. Clim. Chang. 2014, 122, 257–269. [Google Scholar] [CrossRef]
  43. Stips, A.; Macias, D.; Coughlan, C.; Garcia-Gorriz, E.; Liang, X.S. On the causal structure between CO2 and global temperature. Sci. Rep. 2016, 6, 21691. [Google Scholar] [CrossRef] [Green Version]
  44. Liang, X.S. Information flow and causality as rigorous notions ab initio. Phys. Rev. E 2016, 94, 052201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Humlum, O.; Stordahl, K.; Solheim, J.E. The phase relation between atmospheric carbon dioxide and global temperature. Glob. Planet. Chang. 2013, 100, 51–69. [Google Scholar] [CrossRef]
  46. Richardson, M. Comment on “The phase relation between atmospheric carbon dioxide and global temperature” by Humlum, Stordahl and Solheim. Glob. Planet. Chang. 2013, 107, 226–228. [Google Scholar] [CrossRef]
  47. Spencer, R.W.; Christy, J.R. Precise monitoring of global temperature trends from satellites. Science 1990, 247, 1558–1562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Christy, J.R.; Norris, W.B.; Spencer, R.W.; Hnilo, J.J. Tropospheric temperature change since 1979 from tropical radiosonde and satellite measurements. J. Geophys. Res. 2007, 112, D06102. [Google Scholar] [CrossRef]
  49. Jones, P.D.; Lister, D.H.; Osborn, T.J.; Harpham, C.; Salmon, M.; Morice, C.P. Hemispheric and large-scale land surface air temperature variations: An extensive revision and an update to 2010. J. Geophys. Res. 2012, 117, D05127. [Google Scholar] [CrossRef] [Green Version]
  50. Koutsoyiannis, D. Revisiting global hydrological cycle: Is it intensifying? Hydrol. Earth Syst. Sci. 2020. [Google Scholar] [CrossRef]
  51. Keeling, C.D.; Bacastow, R.B.; Bainbridge, A.E.; Ekdahl, C.A.; Guenther, P.R.; Waterman, L.S. Atmospheric carbon dioxide variations at Mauna Loa observatory, Hawaii. Tellus 1976, 28, 538–551. [Google Scholar]
  52. Global Monitoring Laboratory—Carbon Cycle Greenhouse Gases. Available online: https://www.esrl.noaa.gov/gmd/ccgg/about/global_means.html (accessed on 1 September 2020).
  53. Lawrance, A.J. Directionality and reversibility in time series. Int. Stat. Rev. 1991, 59, 67–79. [Google Scholar] [CrossRef]
  54. Weiss, G. Time-reversibility of linear stochastic processes. J. Appl. Probab. 1975, 12, 831–836. [Google Scholar] [CrossRef]
  55. Tong, H.; Zhang, Z. On time-reversibility of multivariate linear processes. Stat. Sin. 2005, 15, 495–504. [Google Scholar]
  56. Koutsoyiannis, D. Generic and parsimonious stochastic modelling for hydrology and beyond. Hydrol. Sci. J. 2016, 61, 225–244. [Google Scholar] [CrossRef]
  57. Kolmogorov, A.N. Wienersche spiralen und einige andere interessante Kurven im Hilbertschen Raum. Dokl. Akad. Nauk SSSR 1940, 26, 115–118. [Google Scholar]
  58. Hurst, H.E. Long term storage capacities of reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 776–808. [Google Scholar]
  59. Heller, M. Time, causality, and the quantum theory. Rev. Metaphys. 1983, 37, 408–409. [Google Scholar]
  60. Kline, A.D. Are there cases of simultaneous causation? PSA Proc. Bienn. Meet. Philos. Sci. Assoc. 1980, 1980, 292–301, Philosophy of Science Association. [Google Scholar] [CrossRef]
  61. Suppes, P. A Probabilistic Theory of Causality; North-Holland Publishing: Amsterdam, The Netherlands, 1970. [Google Scholar]
  62. Granger, C.W. Testing for causality: A personal viewpoint. J. Econ. Dyn. Control 1980, 2, 329–352. [Google Scholar] [CrossRef]
  63. Granger, C.; Newbold, P. Forecasting Economic Time Series; Academic Press: San Diego, CA, USA, 1986. [Google Scholar]
  64. Gujarati, D.N.; Porter, D.C. Basic Econometrics, 5th ed.; McGraw Hill: Boston, MA, USA, 2009. [Google Scholar]
  65. McGraw, M.C.; Barnes, E.A. Memory matters: A case for Granger causality in climate variability studies. J. Clim. 2018, 31, 3289–3300. [Google Scholar] [CrossRef]
  66. Cohen, M.X. Analyzing Neural Time Series Data: Theory and Practice; MIT Press: Cambridge, MA, USA, 2014. [Google Scholar]
  67. Papoulis, A. Probability, Random Variables and Stochastic Processes, 3rd ed.; McGraw-Hill: New York, NY, USA, 1991. [Google Scholar]
  68. Paluš, M.; Krakovská, A.; Jakubík, J.; Chvosteková, M. Causality, dynamical systems and the arrow of time. Chaos Interdiscip. J. Nonlinear Sci. 2018, 28, 75307. [Google Scholar] [CrossRef] [PubMed]
  69. Koutsoyiannis, D.; Yao, H.; Georgakakos, A. Medium-range flow prediction for the Nile: A comparison of stochastic and deterministic methods. Hydrol. Sci. J. 2008, 53, 142–164. [Google Scholar] [CrossRef]
  70. Zaiontz, C. Real Statistics Using Excel. Available online: http://www.realstatistics.com/ (accessed on 1 September 2020).
  71. Zaiontz, C. Real Statistics Examples Workbooks. Available online: http://www.real-statistics.com/free-download/real-statistics-examples-workbook/ (accessed on 1 September 2020).
  72. Kundzewicz, Z.W.; Pińskwar, I.; Koutsoyiannis, D. Variability of global mean annual temperature is significantly influenced by the rhythm of ocean-atmosphere oscillations. Sci. Total Environ. 2020, 747, 141256. [Google Scholar] [CrossRef]
  73. Masters, T.; Benestad, R. Comment on “The phase relation between atmospheric carbon dioxide and global temperature”. Glob. Planet. Chang. 2013, 106, 141–142. [Google Scholar] [CrossRef]
  74. Bond-Lamberty, B.; Thomson, A. Temperature-associated increases in the global soil respiration record. Nature 2010, 464, 579–582. [Google Scholar] [CrossRef] [PubMed]
  75. IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. Available online: http://www.climatechange2013.org/report/ (accessed on 14 February 2020).
  76. Raich, J.W.; Schlesinger, W.H. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate. Tellus B Chem. Phys. Meteorol. 1992, 44, 81–99. [Google Scholar] [CrossRef] [Green Version]
  77. Makita, N.; Kosugi, Y.; Sakabe, A.; Kanazawa, A.; Ohkubo, S.; Tani, M. Seasonal and diurnal patterns of soil respiration in an evergreen coniferous forest: Evidence from six years of observation with automatic chambers. PLoS ONE 2018, 13, e0192622. [Google Scholar] [CrossRef] [Green Version]
  78. ÓhAiseadha, C.; Quinn, G.; Connolly, R.; Connolly, M.; Soon, W. Energy and Climate Policy—An Evaluation of Global Climate Change Expenditure 2011–2018. Energies 2020, 13, 4839. [Google Scholar] [CrossRef]
  79. Pomeroy, R.; Bowlus, F.D. Progress report on sulfide control research. Sew. Work. J. 1946, 18, 597–640. [Google Scholar]
  80. Green, C.; Byrne, K.A. Biomass: Impact on carbon cycle and greenhouse gas emissions. In Encyclopedia of Energy; Cleveland, C.J., Ed.; Elsevier: Amsterdam, The Netherlands, 2004; pp. 223–236. [Google Scholar] [CrossRef]
  81. Weiss, R. Carbon dioxide in water and seawater: The solubility of a non-ideal gas. Mar. Chem. 1974, 2, 203–215. [Google Scholar] [CrossRef]
  82. Zhu, Z.; Piao, S.; Myneni, R.B.; Huang, M.; Zeng, Z.; Canadell, J.G.; Ciais, P.; Sitch, S.; Friedlingstein, P.; Arneth, A.; et al. Greening of the earth and its drivers. Nat. Clim. Chang. 2016, 6, 791–795. [Google Scholar] [CrossRef]
  83. Chen, C.; Park, T.; Wang, X.; Piao, S.; Xu, B.; Chaturvedi, R.K.; Fuchs, R.; Brovkin, V.; Ciais, P.; Fensholt, R.; et al. China and India lead in greening of the world through land-use management. Nat. Sustain. 2019, 2, 122–129. [Google Scholar] [CrossRef] [PubMed]
  84. Connolly, R.; Connolly, M.; Carter, R.M.; Soon, W. How Much Human-Caused Global Warming Should We Expect with Business-As-Usual (BAU) Climate Policies? A Semi-Empirical Assessment. Energies 2020, 13, 1365. [Google Scholar] [CrossRef] [Green Version]
  85. Aryal, Y. Review of Atmospheric Temperature and CO2: Hen-or-Egg Causality? (Version 1). Sci 2020. Available online: https://0-www-mdpi-com.brum.beds.ac.uk/2413-4155/2/4/77/review_report (accessed on 25 September 2020).
  86. Alexandris, S. Review of Atmospheric Temperature and CO2: Hen-or-Egg Causality? (Version 1). Sci 2020. Available online: https://0-www-mdpi-com.brum.beds.ac.uk/2413-4155/2/4/77/review_report (accessed on 26 September 2020).
  87. Beck, E.-G. 180 Years of atmospheric CO2 gas analysis by chemical methods. Energy Environ. 2007, 18, 259–282. [Google Scholar] [CrossRef]
  88. Letts, E.A.; Blake, R.F. The carbonic anhydride of the atmosphere. Proc. R. Soc. Dublin 1900, 9, 107–270. [Google Scholar]
  89. Stepanova, N.A. A Selective Annotated Bibliography of Carbon Dioxide in the Atmosphere. Meteorol. Abstr. 1952, 3, 137–170. [Google Scholar]
  90. Keeling, R.F. Comment on “180 Years of Atmospheric CO2 Gas Analysis by Chemical Methods” by Ernst-Georg Beck. Energy Environ. 2007, 18, 637–639. [Google Scholar]
  91. Meijer, H.A.J. Comment on “180 Years of Atmospheric CO2 Gas Analysis by Chemical Methods” by Ernst-Georg Beck. Energy Environ. 2007, 18, 635–636. [Google Scholar]
  92. Hannart, A.; Pearl, J.; Otto, F.E.L.; Naveau, P.; Ghil, M. Causal counterfactual theory for the attribution of weather and climate-related events. Bull. Am. Meteorol. Soc. 2016, 97, 99–110. [Google Scholar] [CrossRef]
  93. Verbitsky, M.Y.; Mann, M.E.; Steinman, B.A.; Volobuev, D.M. Detecting causality signal in instrumental measurements and climate model simulations: Global warming case study. Geosci. Model Dev. 2019, 12, 4053–4060. [Google Scholar] [CrossRef] [Green Version]
  94. Hume, D. An Enquiry Concerning Human Understanding; Oxford World’s Classics; Oxford University Press: New York, NY, USA, 2008. [Google Scholar]
  95. Koutsoyiannis, D.; Efstratiadis, A.; Mamassis, N.; Christofides, A. On the credibility of climate predictions. Hydrol. Sci. J. 2008, 53, 671–684. [Google Scholar] [CrossRef]
  96. Anagnostopoulos, G.G.; Koutsoyiannis, D.; Christofides, A.; Efstratiadis, A.; Mamassis, N. A comparison of local and aggregated climate model outputs with observed data. Hydrol. Sci. J. 2010, 55, 1094–1110. [Google Scholar] [CrossRef]
  97. Kundzewicz, Z.W.; Stakhiv, E.Z. Are climate models “ready for prime time” in water resources management applications, or is more research needed? Hydrol. Sci. J. 2010, 55, 1085–1089. [Google Scholar] [CrossRef]
  98. Koutsoyiannis, D.; Christofides, A.; Efstratiadis, A.; Anagnostopoulos, G.G.; Mamassis, N. Scientific dialogue on climate: Is it giving black eyes or opening closed eyes? Reply to “A black eye for the Hydrological Sciences Journal” by D. Huard. Hydrol. Sci. J. 2011, 56, 1334–1339. [Google Scholar] [CrossRef] [Green Version]
  99. Tsaknias, D.; Bouziotas, D.; Koutsoyiannis, D. Statistical comparison of observed temperature and rainfall extremes with climate model outputs in the Mediterranean region. ResearchGate 2016. [Google Scholar] [CrossRef]
  100. Tyralis, H.; Koutsoyiannis, D. On the prediction of persistent processes using the output of deterministic models. Hydrol. Sci. J. 2017, 62, 2083–2102. [Google Scholar] [CrossRef]
  101. Koutsoyiannis, D. On the quest for chaotic attractors in hydrological processes. Hydrol. Sci. J. 2006, 51, 1065–1091. [Google Scholar] [CrossRef]
Figure 1. Annual change in global energy-related CO2 emissions (adapted from International Energy Agency (IEA) [3]).
Figure 1. Annual change in global energy-related CO2 emissions (adapted from International Energy Agency (IEA) [3]).
Sci 02 00083 g001
Figure 2. Atmospheric CO2 concentration measured in Mauna Loa, Hawaii, USA, in the last four years.
Figure 2. Atmospheric CO2 concentration measured in Mauna Loa, Hawaii, USA, in the last four years.
Sci 02 00083 g002
Figure 5. Plots of the data series of global temperature “anomalies” since 1980, as used in the study, from satellite measurements over the globe (UAH) and from ground measurements over land (CRUTEM4).
Figure 5. Plots of the data series of global temperature “anomalies” since 1980, as used in the study, from satellite measurements over the globe (UAH) and from ground measurements over land (CRUTEM4).
Sci 02 00083 g005
Figure 6. Plots of the data series of atmospheric CO2 concentration measured in Mauna Loa (Hawaii, USA), Barrow (Alaska, USA), and South Pole, and the global average.
Figure 6. Plots of the data series of atmospheric CO2 concentration measured in Mauna Loa (Hawaii, USA), Barrow (Alaska, USA), and South Pole, and the global average.
Sci 02 00083 g006
Figure 7. Plots of atmospheric CO2 concentration after standardization: (left) each monthly value is standardized by dividing with the geometric average of the 5-year period before it. (right) Monthly statistics of the values of the left panel; for each month, the average is shown in continuous line and the minimum and maximum in thin dashed lines of the same colour as the average.
Figure 7. Plots of atmospheric CO2 concentration after standardization: (left) each monthly value is standardized by dividing with the geometric average of the 5-year period before it. (right) Monthly statistics of the values of the left panel; for each month, the average is shown in continuous line and the minimum and maximum in thin dashed lines of the same colour as the average.
Sci 02 00083 g007
Figure 8. Synchronous plots of the time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa at monthly scale.
Figure 8. Synchronous plots of the time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa at monthly scale.
Sci 02 00083 g008
Figure 9. Auto- and cross-correlograms of the time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa.
Figure 9. Auto- and cross-correlograms of the time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa.
Sci 02 00083 g009
Figure 10. Plots of p-values of the Granger test for 10-year-long moving windows for the monthly time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa for number of lags (left) η = 1 and (right) η = 2. The time series used are (upper) the original and (lower) that obtained after “removing” the periodicity by averaging over the previous 12 months.
Figure 10. Plots of p-values of the Granger test for 10-year-long moving windows for the monthly time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa for number of lags (left) η = 1 and (right) η = 2. The time series used are (upper) the original and (lower) that obtained after “removing” the periodicity by averaging over the previous 12 months.
Sci 02 00083 g010
Figure 11. Differenced time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa at monthly scale. The graph in the upper panel was constructed in the manner described in the text. The graph in the lower panel is given for comparison and was constructed differently by taking differences of the values of each month with the previous month and then averaging over the previous 12 months (to remove periodicity); in addition, the lower graph includes the CRUTEM4 land temperature series.
Figure 11. Differenced time series of UAH temperature and logarithm of CO2 concentration at Mauna Loa at monthly scale. The graph in the upper panel was constructed in the manner described in the text. The graph in the lower panel is given for comparison and was constructed differently by taking differences of the values of each month with the previous month and then averaging over the previous 12 months (to remove periodicity); in addition, the lower graph includes the CRUTEM4 land temperature series.
Sci 02 00083 g011
Figure 12. Annually averaged time series of differenced temperatures (UAH) and logarithms of CO2 concentrations (Mauna Loa). Each dot represents the average of a one-year duration ending at the time of its abscissa.
Figure 12. Annually averaged time series of differenced temperatures (UAH) and logarithms of CO2 concentrations (Mauna Loa). Each dot represents the average of a one-year duration ending at the time of its abscissa.
Sci 02 00083 g012
Figure 13. Empirical climacograms of the indicated differenced time series; the characteristic slopes corresponding to values of the Hurst parameter H = 1/2 (large-scale randomness), 0 (full antipersistence) and 1 (full persistence) are also plotted (note, H = 1 + slope/2).
Figure 13. Empirical climacograms of the indicated differenced time series; the characteristic slopes corresponding to values of the Hurst parameter H = 1/2 (large-scale randomness), 0 (full antipersistence) and 1 (full persistence) are also plotted (note, H = 1 + slope/2).
Sci 02 00083 g013
Figure 14. Auto- and cross-correlograms of the differenced time series of UAH temperature and Mauna Loa CO2 concentration.
Figure 14. Auto- and cross-correlograms of the differenced time series of UAH temperature and Mauna Loa CO2 concentration.
Sci 02 00083 g014
Table 1. Maximum cross-correlation coefficient (MCCC) and corresponding time lag in months. The annual window for temperature is July–June, while for CO2, it is either different (sliding), determined so as to maximize MCCC, or the same (fixed).
Table 1. Maximum cross-correlation coefficient (MCCC) and corresponding time lag in months. The annual window for temperature is July–June, while for CO2, it is either different (sliding), determined so as to maximize MCCC, or the same (fixed).
Temperature—CO2 SeriesMonthly Time SeriesAnnual Time Series—Sliding Annual WindowAnnual Time Series—Fixed Annual Window
MCCCLagMCCCLagMCCCLag
UAH—Mauna Loa0.4750.6680.5212
UAH—Barrow0.31110.70140.5912
UAH—South Pole0.3760.54100.3812
UAH—Global0.4760.60110.6012
CRUTEM4—Mauna Loa0.3150.55100.5212
CRUTEM4—Global0.3390.55120.5512
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koutsoyiannis, D.; Kundzewicz, Z.W. Atmospheric Temperature and CO2: Hen-Or-Egg Causality? Sci 2020, 2, 83. https://0-doi-org.brum.beds.ac.uk/10.3390/sci2040083

AMA Style

Koutsoyiannis D, Kundzewicz ZW. Atmospheric Temperature and CO2: Hen-Or-Egg Causality? Sci. 2020; 2(4):83. https://0-doi-org.brum.beds.ac.uk/10.3390/sci2040083

Chicago/Turabian Style

Koutsoyiannis, Demetris, and Zbigniew W. Kundzewicz. 2020. "Atmospheric Temperature and CO2: Hen-Or-Egg Causality?" Sci 2, no. 4: 83. https://0-doi-org.brum.beds.ac.uk/10.3390/sci2040083

Article Metrics

Back to TopTop