Sunday, December 15, 2013

The insignificant effect of CO2 on climate is quantified at  with a near-perfect explanation of annual average global temperatures since before 1900; R^2 = 0.97+

Calculated Mean Global Temperatures 1610-2012

This monograph considers only average global temperature (AGT). It does not discuss weather, which is a complex study of energy moving about the planet. It does not even address local climate, which includes precipitation. It does, however, consider the issue of Global Warming and the tenaciously held but mistaken perception that human activity has a significant influence on it.

The word ‘trend’ is used here for temperatures in two different contexts. To differentiate, α-trend applies to averaging-out the uncertainties in reported average global temperature measurements to produce the average global temperature oscillation resulting from the net ocean surface oscillation. The term β-trend applies to the slower average temperature change of the planet which is associated with change to the average temperature of the bulk volume of the material (mostly water) involved.

The first paper to suggest the hypothesis that the sunspot number time-integral is a proxy for a substantial driver of average global temperature change was made public 6/1/2009. The discovery started with application of the first law of thermodynamics, conservation of energy, and the hypothesis that the energy acquired, above or below break-even (appropriately accounting for energy radiated from the planet), is proportional to the time-integral of sunspot numbers. The derived equation revealed a rapid and sustained global energy rise starting in about 1941. The true average global temperature anomaly change β-trend is proportional to global energy change.

Subsequent analysis revealed that the significant factor in calculating the β-trend is the sunspot number anomaly time-integral. The sunspot number anomaly is defined as the difference between the sunspot number in a specific year and an average sunspot number for many years.

Measured surface temperature anomaly α-trends oscillate above and below the temperature anomaly β-trend calculated using only the sunspot number anomaly time-integral.

The existence of multiple ocean oscillations led to the perception that there must be an effective net surface temperature oscillation for the planet with all named and unnamed ocean oscillations as participants. Plots of measured average global temperatures indicate that the net surface temperature oscillation has a period of 64 years with the most recent maximum in 2005.

Combination of the effects results in the effect of the ocean surface temperature oscillation (α-trend) decline 1941-1973 being slightly stronger than the effect of the rapid rise from sunspots (β-trend) resulting in a slight decline of the trend of reported average global temperatures. The steep rise 1973-2005 occurred because the effects added. A high coefficient of determination, R2, demonstrates that the hypothesis is true.

Over the years, several refinements to this work (often resulting from other's comments which may or may not have been corroborative) slightly improved the accuracy and led to the equations and figures in this paper.

Prior work
The law of conservation of energy is applied effectively the same as described in Reference 2 in the development of a very similar equation which calculates temperature anomalies. The difference is that the variation in energy ‘OUT’ has been found to be adequately accounted for by variation of the sunspot number anomalies. Thus the influence of the factor [T(i)/Tavg]4 is eliminated.

Change to the level of atmospheric carbon dioxide has no significant effect on average global temperature. This was demonstrated in 2008 at Reference 6 and is corroborated at Reference 2 and again here.  

As determined in Reference 3, reported average global temperature anomaly measurements have a random uncertainty with equivalent standard deviation ≈ 0.09 K. A substantial contributor to this variation appears to be the apparently random variation in magnitude and period of el Nino.

Global Warming ended more than a decade ago as shown here, and in Reference 4 and also Reference 2.

Average global temperature is very sensitive to cloud change as shown in Reference 5. An increase of approximately 186 meters in average cloud altitude or a decrease of average albedo from 0.3 to the very slightly reduced value of 0.2928 would account for all of the 20th century increase in AGT of 0.74 °C.

The value used for average sunspot number was 43.97 (average 1850-1940) in Ref. 1. It is set at 34 (average 1610-1940) in this paper. The procession of values for average sunspot number from 43.97 to 34 produces slight but steady improvement in R2 for the period of measured temperatures, and, progressively greater credibility of average global temperature estimates for the period prior to accurate, world wide, direct measurements becoming available.

Initial work is presented in several papers made public at

Ocean surface temperature oscillation
The ocean oscillation does not significantly add or remove planet energy. In the decades immediately prior to 1941 the amplitude range of the trends was not significantly influenced by any candidate internal forcing effect; so the observed amplitude of the effect on AGT of the net ocean surface temperature trend anomaly then, must be approximately the same as the amplitude of the part of the AGT trend anomaly due to ocean oscillations since then. This part is approximately 0.36 K total range with a period of approximately 64 years (verified below).

The AGT trajectory (Figure 1) suggests that the least-biased simple wave form of the ocean surface temperature oscillation is approximately saw-toothed. Ignoring the offset for the moment, the sea surface temperature anomaly oscillation can be described as varying linearly from 0.0 K in 1909 to approximately 0.36 K in 1941 and linearly back to the 1909 value in 1973. This cycle repeats before and after with a period of 64 years. This is consistent with the 50-70 year period previously observed by others.

Because the actual magnitude of the ocean oscillation in any year is needed, the expression to account for the contribution of the ocean oscillation to AGT is given by the following:

ΔTosc = (A,y)             K (degrees)                 (i)

where the contribution of the net of ocean oscillations to AGT change is the magnitude of the effect on AGT of the surface temperature anomaly trend of the oscillation in year y, and A is the maximum magnitude of the effect on AGT of the ocean surface temperature trend oscillation. 

Equation (i) is graphed in Figure 0.5.
Figure 0.5: Ocean surface temperature oscillations do not affect the bulk energy of the planet.

Although the peak-to-peak amplitude of the effect on AGT of the ocean oscillation is approximately 0.36 K, by definition, ocean oscillation is symmetrical with respect to zero. Therefore the contribution of Sea Surface Temperature (SST) trend oscillation to AGT since before 1900 has been approximately ±0.18 K.

The sunspot number anomaly time-integral drives the temperature anomaly β-trend
By definition, energy change divided by effective thermal capacitance is temperature change.

Table 1 in reference 2 shows the influence of atmospheric carbon dioxide (CO2) to be insignificant (tiny change in R2 if considering CO2 or not; also shown later in this paper to be insignificant) so it can be removed from the equation by setting coefficient ‘C’ to zero. With ‘C’ set to zero, Equation 1 in Reference 2 calculates average global temperature anomalies (AGT) since 1895 with 89.82% accuracy (R2 = 0.898220).

The current analysis determined that 34, the approximate average of sunspot numbers from 1610-1940, provides a slightly better fit (in fact, the best fit) to the measured temperature data. The influence, of Stephan-Boltzmann radiation change due to AGT change, on energy change is adequately accounted for by the sunspot number anomaly time-integral. With these refinements to Equation (1) in Reference 2 the coefficients become A = 0.3588, B = 0.003461 and D = ‑ 0.4485.  R2 increases slightly to 0.904906 and the calculated anomaly in 2005 is 0.5045 K. Also with these refinements the equation calculates lower early temperature anomalies and projects a slightly higher (0.3175 K vs. 0.269 K in 2020) future anomaly trend. The resulting equation for calculating the AGT anomaly for any year, 1895 or later, is then:

Anom(y) = (0.3588,y) + 0.003461/17 Σyi=1895 (s(i) – 34) – 0.4485    (ii)

            Anom(y) = calculated temperature anomaly in year y, K
            (0.3588,y) = approximate contribution of ocean cycle effect to AGT in year y
            s(i) = average daily Brussels International sunspot number in year i
            17 = effective thermal capacitance, W yr m-2 K-1

Measured temperature anomalies are from Figure 2 of Reference 3. The excellent match of the up and down trends since before 1900 of calculated and measured temperature anomalies, shown here in Figure 1, and, for 5-year moving average temperature anomaly measurements, in Figure 1.1, demonstrate the usefulness and validity of the calculations. Figure 1.2 uses the same data as Figure 1 (which is prior to the corrupting ‘changes’ by the reporting agencies) with reported SSN and ‘best guess’ estimate of AGT anomalies for the years 2013-2015.

Projections until 2020 use the expected sunspot number trend for the remainder of solar cycle 24 as provided 11 by NASA. After 2020 the limiting cases are either assuming sunspots like from 1925 to 1941 or for the case of no sunspots which is similar to the Maunder Minimum.

Some noteworthy volcanoes and the year they occurred are also shown on Figure 1. No consistent AGT response is observed to be associated with these. Any global temperature perturbation that might have been caused by volcanoes of this size is lost in the temperature measurement uncertainty. Note: The uncertainty is not in the method, or the measuring instruments themselves, but results from the effectively roiling (at this scale) of the object of the measurements.

Much larger volcanoes can cause significant temporary global cooling from the added reflectivity of aerosols and airborne particulates. The Tambora eruption, which started on April 10, 1815 and continued to erupt for at least 6 months, was approximately ten times the magnitude of the next largest in recorded history and led to 1816 which has been referred to as ‘the year without a summer’. The cooling effect of that volcano exacerbated the already cool temperatures associated with the Dalton Minimum.
 Figure 1: Measured average global temperature anomalies with calculated prior and future trends using 34 as the average daily sunspot number. (Last update 12/15/13)

Figure 1.1: Same as Figure 1 but with 5-year running average of measured temperatures. R2 = 0.973. (Added 5/26/15)

Figure 1.2: Same as Figure 1 but with 5-year running average of measured temperatures and extended through 2015. R2 = 0.973. (Added 1/6/16)

As discussed in Reference 2, ocean oscillations produce oscillations of the ocean surface temperature with no significant change to the average temperature of the bulk volume of water involved. The effect on reported AGT of the full range of surface temperature oscillation is given by the coefficient ‘A’. (A, B, C, and D are the coefficients in Equation 1 of Reference 2)

The influence of ocean surface temperature oscillations can be removed from the equation by setting ‘A’ to zero. To use all regularly recorded sunspot numbers, the integration starts in 1610. The offset, ‘D’ must be changed to -0.1993 to account for the different integration start point and setting ‘A’ to zero. Setting ‘A’ to zero requires that the anomaly in 2005 be 0.5045 - 0.3588/2 = 0.3251 K. The result, Equation (1) here, then calculates the trend 1610-2012 resulting from just the sunspot number anomaly time-integral.

Trend3anom(y) = 0.003461/17 * Σyi = 1610 [s(i)-34] – 0.1993                (1)

Trend3anom(y) = calculated temperature anomaly β-trend in year y, K degrees.
0.003461 = the proxy factor, B, W yr m-2.
17 = effective thermal capacitance of the planet, W Yr m-2 K-1
s(i) = average daily Brussels International sunspot number in year i
34 ≈ average sunspot number for 1610-1940.
-0.1993 is merely an offset that shifts the calculated trajectory vertically on the graph, without changing its shape, so that the calculated temperature anomaly in 2005 is 0.3251 K which is the calculated anomaly for 2005 if the ocean oscillation is not included.

Sunspot numbers back to 1610 are shown in Figure 2 of Reference 1.

Applying Equation (1) to the sunspot numbers of Figure 2 of Reference 1 produces the trace shown in Figure 2 below.

Figure 2: Anomaly trend (β-trend) from just the sunspot number anomaly time-integral using Equation (1).

Average global temperatures were not directly measured in 1610 (thermometers had not been invented yet). Recent estimates, using proxies, are few. The temperature anomaly trend that Equation (1) calculates for that time is roughly consistent with other estimates. The decline in the trace 1610-1700 on Figure 2 results from the low sunspot numbers for that period as shown on Figure 2 of Reference 1. 

How this phenomenon could take place
Although the connection between AGT and the sunspot number anomaly time-integral is demonstrated, the mechanism by which this takes place remains somewhat speculative.

Various papers have been written that indicate how the solar magnetic field associated with sunspots can influence climate on earth. These papers posit that decreased sunspots are associated with decreased solar magnetic field which decreases the deflection of and therefore increases the flow of galactic cosmic rays on earth.

Henrik Svensmark, a Danish physicist, found that decreased galactic cosmic rays caused decreased low level (<3 km) clouds and planet warming. An abstract of his 2000 paper is at Reference 13. Marsden and Lingenfelter also report this in the summary of their 2003 paper 14 where they make the statement “…solar activity increases…providing more shielding…less low-level cloud cover… increase surface air temperature.”  These findings have been further corroborated by the cloud nucleation experiments 15 at CERN.

These papers associated the increased low-level clouds with increased albedo leading to lower temperatures. Increased low clouds would also result in lower average cloud altitude and therefore higher average cloud temperature. Although clouds are commonly acknowledged to increase albedo, they also radiate energy to space so increasing their temperature increases radiation to space which would cause the planet to cool. Increased albedo reduces the energy received by the planet and increased radiation to space reduces the energy of the planet. Thus the two effects work together to change the AGT of the planet.

Simple analyses 5 indicate that either an increase of approximately 186 meters in average cloud altitude or a decrease of average albedo from 0.3 to the very slightly reduced value of 0.2928 would account for all of the 20th century increase in AGT of 0.74 °C. Because the cloud effects work together and part of the temperature change is due to ocean oscillation, substantially less cloud change would suffice.

Hind Cast Estimate of Combined Sunspot Effect and Ocean Oscillation Effect
As a possibility, the period and amplitude of oscillations attributed to ocean cycles demonstrated to be valid after 1895 are assumed to maintain back to 1610. Equation (1) is modified as shown in Equation (2) to account for including the effects of ocean oscillations. Since the expression for the oscillations calculates values from zero to the full range but oscillations must be centered on zero, it must be reduced by half the oscillation range.

Trend4anom(y) = (0.3588,y) – 0.1794 + 0.003461/17 * Σyi = 1610 [s(i)-34] – 0.1993   (2)

The ocean oscillation factor, (0.3588,y) – 0.1794, is applied to the period prior to the start of direct temperature measurements as a possibility. The effective sea surface temperature anomaly, (A,y), is defined above and in Reference 2.

Applying Equation (2) to the sunspot numbers from Figure 2 of Reference 1 produces the trend shown in Figure 3 next below. Available measured average global temperatures from Figure 2 in Reference 3 are superimposed on the calculated values.
 Figure 3: Calculated temperature anomalies from the sunspot number anomaly time-integral plus ocean oscillation using Equation (2) with superimposed available measured data from Reference 3 and range estimates determined by Loehle.

Figure 3 shows that temperature anomalies calculated using Equation (2) estimate possible trends since 1610 and actual trends of reported temperatures since they have been accurately measured world wide.

The match 1895-2012 has R2 = 0.9049 which means that 90.49% of average global temperature anomaly measurements are explained. All factors not explicitly considered (such as the 0.09 K s.d. random uncertainty in reported annual measured temperature anomalies, aerosols, CO2, other non-condensing ghg, volcanoes, ice change, etc.) must find room in that unexplained 9.51%. Note that a coefficient of determination, R2 = 0.9049 means a correlation coefficient of 0.95.

A survey 12 of non-tree-ring global temperature estimates was conducted by Loehle including some for a period after 1610. A simplification of the 95% limits found by Loehle are also shown on Figure 3. The spread between the upper and lower 95% limits are fixed, but, since the anomaly reference temperatures might be different, the limits are adjusted vertically to approximately bracket the values calculated using the equations. The fit appears reasonable considering the uncertainty in all values.

Smoothing of the measured temperatures using 5-year moving average, as shown in Figure 1.1, achieved R2 = 0.973. This accounts for most of the random uncertainty in reported annual measured temperature anomalies with only 2.7% left unexplained.

Calculated temperature anomalies look reasonable back to 1700 but indicate higher temperatures prior to that than most proxy estimates. They are, however, consistent with the low  sunspot numbers in that period. They qualitatively agree with Vostok, Antarctica ice core data but decidedly differ from Sargasso Sea estimates during that time (see the graph for the last 1000 years in Reference 6). Worldwide assessments of average global temperature that far back are sparse and speculative. Ocean oscillations might also have been different from assumed.

Possible lower values for average sunspot number
Possible lower assumed values for average sunspot number, with coefficients adjusted to maximize R2, result in noticeably lower estimates of early (prior to direct measurement) temperatures with only a tiny decrease in R2. Calculated temperature anomalies resulting from using an average sunspot number value of 26 are shown in Figure 4. The projected temperature anomaly trend decline is very slightly less steep (0.018 K warmer in 2020) than was shown in Figure 1.

Figure 4: Calculated temperature anomalies from the sunspot number anomaly time-integral plus ocean oscillation using 26 as the average sunspot number with superimposed available measured data from Reference 3 and range estimates determined by Loehle.

Carbon dioxide change has no significant influence
The influence that CO2 has on AGT can be calculated by including ‘C’ in Equation (1) of Reference 2 as a coefficient to be determined. The tiny increase in R2 demonstrates that consideration of change to the CO2 level has no significant influence on AGT. The coefficients and resulting R2 are given in Table 1.

Table 1: A, B, C, D, refer to coefficients in Equation 1 in Reference 2
Average daily SSN
ocean oscillation A
sunspots B
Coefficient of determination R2
% cause of 1909-2005 AGT change (0.9465 K)
Ocean oscillation
CO2 change
*Measured temperatures smoothed using 5-year moving average. (Added 5/26/15, rev 12/13/15)
**Equation calibrated through 1990
***Svalgaard 5 5 16 SSN, Average daily SSN set to 63

††Extended through 2015

The coefficient ‘B’ is effectively a combination of proxy factor and influence coefficient.

Assessment using Svalgaard SSN (rev 5/6/16)
Figure 5: SSN used in assessment results in Table 1 for Svalgaard 5 5 16. R2 = 0.968900 (1895-2012). (Pre 1700 SSN are Brussels International SSN/0.6).
Figure 6: Anomaly trend (β-trend) from just the sunspot number anomaly time-integral using SSN from Figure 5 (Svalgaard)

Possible explanation of why CO2 change has no significant effect on climate. (Added 10/2/14, revised 11/21/14, 11/29/14, 3/15/15)
1) Firmly acknowledge the established fact that gas molecules can absorb/emit photons only at specific discreet wavelengths (which might be broadened from pressure, etc.). This fact makes spectroscopy possible. Full spectrum (Plank’s law) Stephan-Boltzmann (S-B) radiation applies to liquids and solids, not to gases.
2) From gas kinetics, the time between atmospheric molecule collisions is extremely short (The Hyperphysics calculator calculates approximately 0.0001 microsecond at sea level pressure and temperature).
3) The elapsed time between absorption and emission of a photon by a CO2 gas molecule is perhaps shorter at higher temperature but must be greater than zero or there would be no evidence that absorption-emission had occurred.
4) At sea level conditions, some or all of the photon energy that is absorbed by a ghg molecule is immediately transferred to other molecules by collision. The process of absorbing a photon and transferring (thermal conduction in the gas) the added energy to other molecules is thermalization. A common observation of thermalization by way of water vapor is that cloudless nights cool faster when absolute water vapor content is lower.
5) The reduced radiation flux on both sides of the 15 micron CO2 absorption line, as observed in most Top Of Atmosphere (TOA) measurements 18 results because some of the EMR energy absorbed by CO2 has been thermalized.
6) Terrestrial radiation is nearly all in the wavelength range 6-100 microns. Thermalized energy carries no identity of the molecule that absorbed it.
7) Jostling between the molecules sometimes causes reverse-thermalization. At low to medium altitudes, EMR emission stimulated by reverse-thermalization is mostly by way of water vapor. The TOA spike at 15 microns results from reverse-thermalization to CO2  molecules at very high altitude.
8) The thermalized radiation warms the air, reducing its density, causing updrafts which are exploited by soaring birds, sailplanes, and occasionally hail. Updrafts are matched by downdrafts elsewhere, usually spread out but sometimes recognized by pilots and passengers as ‘air pockets’ and micro bursts.
9) The population gradient of ghg molecules, (especially water vapor above about 3 km, declining with increasing altitude) favors radiation to space. Ghg molecules that emit a photon are ‘recharged’ by reverse-thermalization. (or by absorbing a photon of appropriate wavelength)
10) Clouds (average emissivity about 0.5) consist of solid and/or liquid water particles (each particle containing millions of molecules) that radiate according to S-B law. Low amount of water vapor above clouds and widening molecule spacing allows substantial radiation directly to space.
11) The tiny increase in ghg from increased CO2 causes absorption/thermalization to occur at slightly lower altitude which very slightly increases the convection rate.

12) The increase in absorbing molecules near the surface is apparently compensated for by an equal increase in emitting molecules high in the atmosphere radiating energy from the planet.

Further discussion of ocean cycles (Added 6/23/14)
Ocean cycles are perceived to contribute to AGT in two ways: The first is the direct measurement of sea surface temperature (SST). The second is warmer SST increases atmospheric water vapor which acts as a forcing and therefore has a time-integral effect on temperature

The combined effect of temperature contribution to AGT of ocean cycles is approximated by a function that has a saw-tooth trajectory profile. It is represented earlier in this paper and in Equation (1) of Reference 2 by (A,y) where A is the total amplitude and y is the year. The uptrends and down trends are each determined to be 32 years long for a total period of 64 years. The total amplitude resulting from ocean oscillations was found here to be 0.3588 K for unsmoothed and 0.3482 for smoothed (cases highlighted in Table 1).

Thus, for an ocean cycle surface temperature uptrend, the contribution of ocean oscillations to AGT is approximated by adding (to the value calculated from the sunspot number anomaly time-integral) 0.3588 multiplied by the fraction of the 32 year period that has elapsed since a low. For an ocean cycle surface temperature down trend, the contribution is calculated by adding 0.3588 minus 0.3588 multiplied by the fraction of the 32 year period that has elapsed since a high. The calculated lows were found to be in 1909 and 1973 and the highs in 1941 and 2005. The resulting trajectory is shown in Figure 0.5, and, offset by half the amplitude, is shown as ‘approximation’ in Figure 5.

Temperature data is available for three named cycles: PDO index, ENSO 3.4 index and AMO index. Successful accounting for oscillations is achieved for PDO and ENSO when considering these as forcings (with appropriate proxy factors) instead of direct measurements. As forcings, their influence accumulates with time. The proxy factors must be determined separately for each forcing. The measurements are available since 1900 for PDO 16 and ENSO3.4 17. This PDO data set has the PDO temperature measurements reduced by the average SST measurements for the planet.

The contribution of PDO and ENSO3.4 to AGT is calculated by:
PDO_NINO = Σyi=1900 (0.017*PDO(i) + 0.009 * ENSO34(i))        (3)

            PDO(i) = PDO index 16 in year i
            ENSO34(i) = ENSO 3.4 index 17 in year i

How this calculation compares to the idealized approximation used in Equation (2) is shown in Figure 5. The high coefficient of determination in Table 1 and the comparison in Figure 5 corroborate the assumption that the saw-tooth profile provides an adequate approximation of the influence of all named and unnamed ocean cycles in the calculated AGT anomalies.
Figure 5: Comparison of idealized approximation of ocean cycle effect and the calculated effect from PDO and ENSO.

(added 9/21/2015)
The AMO index 22 is formed from area-weighted and detrended SST data. It is shown with two different amounts of smoothing in Figure 6 along with the saw-tooth approximation for the entire planet.
Figure 6: Comparison of idealized approximation of ocean cycle effect and the AMO index.

The high coefficient of determination in Table 1 and the comparisons in Figure 5 and 6 corroborate the assumption that the saw-tooth profile provides adequate approximation of the net effect of all named and unnamed ocean cycles in the calculated AGT anomalies.

Others that have looked at only amplitude or only duration factors for solar cycles got poor correlations with average global temperature. The good correlation comes by combining the two, which is what the time-integral of sunspot number anomalies does. As shown in Figure 2, the temperature anomaly trend, determined using the sunspot number anomaly time-integral, exhibits substantial change over the recorded period. Prediction of future sunspot numbers more than a decade or so into the future has not yet been confidently done although assessments using planetary synodic periods appear to be relevant 7,8.

As displayed in Figure 2, the time-integral of sunspot number anomalies alone appears to show the estimated true average global temperature trend (the net average global energy trend) during the planet warm up from the depths of the Little Ice Age.

The net effect of ocean oscillations is to cause the surface temperature α-trend to oscillate above and below the β-trend calculated using only the sunspot number anomaly time-integral. Equation (2) accounts for both and also, because it matches measurements so well, shows that rational change to the level of atmospheric carbon dioxide can have no significant influence.

Figure 1.1 shows the near perfect match with calculated temperatures which occurs when random fluctuation in reported measured temperatures is smoothed out with 5-year moving average.

Long term prediction of average global temperatures depends primarily on long term prediction of sunspot numbers.

9. Deleted.
10. Deleted.
11. Graphical sunspot number prediction for the remainder of solar cycle 24
13. Svensmark paper, Phys. Rev. Lett. 85, 5004–5007 (2000)
14. Marsden & Lingenfelter 2003, Journal of the Atmospheric Sciences 60: 626-636
18. Barrett, ‘Greenhouse molecules, their spectra and function in the atmosphere’, Energy & Environment, Vol. 16, No. 6, 2005.
19. Deleted
20. Deleted
21 Deleted.


  1. Your Reference 10 note is unhelpful to one without context. I suggest a link. And while I came to this article via a comment thread, the idea of it lying fallow without added commentary and updates seems a waste. Do you have a Disqus link to add ?

  2. The Reference list with links is at the end of the article.

    I am unfamiliar with 'Disqus link'. Please clue me in through email link accessed through 'View my complete profile' in upper right corner of article. The email link is displayed on the left side just below where my picture isn't. Thanks.

  3. Thanks Dan!
    This is the kind of information that we all need to read. I am totally frustrated with the Global Warming religion and their zealotry... I have resorted to cynicism on all my blogs. Thanks for your post! Keep up the good work!

    1. Thanks for the comments.

      The warmers lose cred as the CO2 continues to go up and average global temperature doesn't.

  4. 1816? The year without a summer? Man, that must have sucked.

  5. Replies
    1. Certainly a detailed piece. I just cannot bring myself to believe in man-made global warming. Wrote a satirical post some time back you might enjoy Paleo Indians and Global Warming. Almost forgot, thanks for your recent visit to one of my blogs (Right Wing Humor)>

  6. I have only started reading this after reading your comment on WUWT so I haven't much to add except a little criticism about the axis titles. I guess for blogs the units should be (%deg;C) or divide by K (and not K degrees as you have in one graph). I know its difficult to choose the right option for the common reader but mixing things up is not a good idea.

  7. The recommended symbol for degrees Kelvin is K with no degree symbol. Since a Celcius degree is the same size as a Kelvin degree that could have been used but requires also the degree sign. Practice in science work is to use K for either the temperature or the temperature difference which is what anomalies are. All of the graphs have the ordinate (anomaly) in degrees Kelvin using the symbol K. I added 'degrees' on the first graph in case someone might think that the K stood for 1000.

  8. Dan - your method produces results very similar to mine which is based simply on the working hypothesis that the recent peak in temperatures is a result of a synchronous peak in the 60 year and 1000 year quasi- periodicities in the temperature data and the use of the neutron count (and 10Be) as the most useful proxy for solar activity over the long term. For several posts on this subject and for forecasts of the timing and amount of the coming cooling see
    Here are the conclusions of the last post on the link
    "In earlier posts on this site at 4/02/13 and 1/22/13
    I have combined the PDO, ,Millennial cycle and neutron trends to estimate the timing and extent of the coming cooling in both the Northern Hemisphere and Globally.
    Here are the conclusions of those posts.
    1/22/13 (NH)
    1) The millennial peak is sharp - perhaps 18 years +/-. We have now had 16 years since 1997 with no net warming - and so might expect a sharp drop in a year or two - 2014/16 -with a net cooling by 2035 of about 0.35.Within that time frame however there could well be some exceptional years with NH temperatures +/- 0.25 degrees colder than that.
    2) The cooling gradient might be fairly steep down to the Oort minimum equivalent which would occur about 2100. (about 1100 on Fig 5) ( Fig 3 here) with a total cooling in 2100 from the present estimated at about 1.2 +/-
    3) From 2100 on through the Wolf and Sporer minima equivalents with intervening highs to the Maunder Minimum equivalent which could occur from about 2600 - 2700 a further net cooling of about 0.7 degrees could occur for a total drop of 1.9 +/- degrees
    4)The time frame for the significant cooling in 2014 - 16 is strengthened by recent developments already seen in solar activity. With a time lag of about 12 years between the solar driver proxy and climate we should see the effects of the sharp drop in the Ap Index which took place in 2004/5 in 2016-17.
    4/02/13 ( Global)
    1 Significant temperature drop at about 2016-17
    2 Possible unusual cold snap 2021-22
    3 Built in cooling trend until at least 2024
    4 Temperature Hadsst3 moving average anomaly 2035 - 0.15
    5 Temperature Hadsst3 moving average anomaly 2100 - 0.5
    6 General Conclusion - by 2100 all the 20th century temperature rise will have been reversed,
    7 By 2650 earth could possibly be back to the depths of the little ice age.
    8 The effect of increasing CO2 emissions will be minor but beneficial - they may slightly ameliorate the forecast cooling and help maintain crop yields .
    9 Warning !! There are some signs in the Livingston and Penn Solar data that a sudden drop to the Maunder Minimum Little Ice Age temperatures could be imminent - with a much more rapid and economically disruptive cooling than that forecast above which may turn out to be a best case scenario.

    How confident should one be in these above predictions? The pattern method doesn't lend itself easily to statistical measures. However statistical calculations only provide an apparent rigor for the uninitiated and in relation to the IPCC climate models are entirely misleading because they make no allowance for the structural uncertainties in the model set up.This is where scientific judgment comes in - some people are better at pattern recognition and meaningful correlation than others. A past record of successful forecasting such as indicated above is a useful but not infallible measure. In this case I am reasonably sure - say 65/35 for about 20 years ahead. Beyond that certainty drops rapidly. I am sure, however, that it will prove closer to reality than anything put out by the IPCC, Met Office or the NASA group. In any case this is a Bayesian type forecast- in that it can easily be amended on an ongoing basis as the Temperature and Solar data accumulate. If there is not a 0.15 - 0.20. drop in Global SSTs by 2018 -20 I would need to re-evaluate.

  9. Dr. Page,
    My assessment is an extrapolation of the past. I demonstrated that CO2 change had no significant effect in a paper made public in 2008 at The determination of the two climate drivers appears solid. Prediction beyond 2037 depends on prediction of sunspot numbers and ocean oscillations, neither of which have been confidently predicted very far into the future.

    I don't like long moving averages because they tend to obscure trend peaks. Long moving averages can be misleading, particularly in noisy data.

    Although we both predict future declines, it looks like my decline prediction for 2035 is about twice what you came up with. I haven't looked beyond 2037 but suspect sunspots may be related to planetary synodic periods and ocean oscillations may be influenced by some beat frequency with the lunar cycles.

    I believe that future-cast must also work as hind-cast and simple is better than complex. Predictions must be physics based.

    I agree that what the IPCC is doing has little to do with climate.

  10. Dan, I don't know if you do the Facebook thing or not, but here's an interesting thing about Bill Nye and Global Warming:

    Facebook Bill Nye Global Warming

  11. JB - Thanks for the link. I added a comment.

  12. I saw that. There was a guy named Christopher Robert that said this to you:

    Quote" I'll give you a hint: if it ends in .com and it's not a real journal (,, etc.), it's not a source, and no real scientist will ever even click your links. Richard Feynman would roll in his grave if he knew crazies like you were trying to cite him. Citing a blogspot to try to argue science....sheesh."Quote

  13. Wow, some very straight forward info for a change. Instead of the blah blah blah global warming because blah blah blah.

    1. Never knew the point on c02 either, as temps don't change. All $$$$

    2. GlobalWarmingTimes A Global Warming News Site offring each and every latest news and information about Global Warming I just follow this blog and read it daily basis it will keep you updats.. Global Warming News

    3. Mahesh - It might keep you misinformed.

  14. The CO2 data that was used (and change to which was shown to have no significant effect) is from Law Dome Antarctica 1180-1958 and Mauna Loa Hawaii after 1958. A graph is on page 11 at
    Several other early (most superseded) papers are at

  15. Thanks Dan,

    Best article I've read in weeks. Nice to see real mathematics.

  16. Thanks for the feedback.

    Here is more. The physics of thermalization explains why non-condensing ghg change has no significant effect on climate.

  17. A really good article, thank you very much.

  18. Your valuable blog aware us from the side effect of global warming.

    Thanks for the sharing a awaking blog on the global warming times.

    Global Warming Times

  19. Dan Pangburn: "The warmers lose cred as the CO2 continues to go up and average global temperature doesn't." Precisely. The counter by claiming the current hiatus is too short of a time period. All of the hullabaloo for the past year is an attempt to downplay and discredit the idea (data) and supporters of the hiatus (supporting real data is an contradiction in terms itself).

    I looked at this on WUWT in two related posts in which I look at the simple cross plot of NASA GISS mean monthly temperature vs. Mauna Loa Keeling CO2.

    Key points:
    Taking the entire data set (back to 1959, limited by the Mauna Loa data) the X-Y plot clearly breaks into three distinct periods a) 1959 - mid 1970s flat or no correlation; b) mid 1970s to late 1990s global mean temperature monotonically increase suggesting apparent correlation with with CO2; c) late 1990s to date - current hiatus with zero correlation between CO2 and global mean temperature. Re the first period in fact from the mid 1940s to mid 1970s not only was there was no relation between temperature and CO2 but even cooling during the period. The point in the second WUWT article is if manmade CO2 is the predominate factor "causing" global warming, this conclusion is not supported by the data. Therefore it is not a conclusion it is an unsupported premise. The data during the last period late 1990s to date is remarkable in that the X-Y data are a complete scattergram/shotgun with zero correlation.

    My conclusion is: "The accepted (by the consensus) hypothesis that global mean temperature (the dependent variable) can be explained by or is due to “mainly” a single variable, CO2 is patently false during the 18+ year hiatus...(d)id CO2 sensitivity go to sleep? ...(a)re other variables exactly canceling out the CO2 effect? It is also important to recognize that the Mauna Loa data includes manmade and non-manmade CO2. The integrated assessment climate models (IAMs) are deterministic physical models of the climate with built in predetermined physical model structures (i.e., defining of cause and effect). From the data, we can say they (models and their proponents) are wrong based on their (in)ability to explain the facts (data) during this hiatus. The zero “fit” / complete lack of correlation does allow us to conclude that during the hiatus, the assumption that CO2 is the major thing driving global mean temperature is not just a lousy hypothesis, it’s wrong and unsupported by the data (fact). The data present strong evidence that all of the variability (scatter) in the data during the hiatus is due to “not CO2.”

  20. You really share a valuable blog on the effect of global warming.
    The data and fact that you have been presented for the better understand the global warming effect is correct and easy to get.
    Global Warming News

  21. Pity it doesn't match reality.
    How are polar temperatures going ? Classic GHG phenomenon of warmer when little or no sun or at night.
    Record heat worldwide for last few years.
    First direct observation of carbon dioxide's increasing greenhouse effect at Earth's surface

    February 25, 2015
    Lawrence Berkeley National Laboratory
    Scientists have observed an increase in carbon dioxide's greenhouse effect at Earth's surface for the first time. They measured atmospheric carbon dioxide's increasing capacity to absorb thermal radiation emitted from Earth's surface over an 11-year period at two locations in North America. They attributed this upward trend to rising carbon dioxide levels from fossil fuel emissions.
    0.2 Watts per Sq meter over 11 years during which time CO2 rose 22ppm, since late 1700's CO 2 has risen 140 approx ppm largely due to man along with feedbacks and the starting to melt methane Clathrates. Just mans activities has increased CO2 back radiation by approx 1.4 Watts per Sq Meter,

    A maunder Minimum reduces incident energy by approx 0.5 Watts per Sq Meter.
    Do your own math

    1. Perhaps you did not read the analysis. The match is 97% with measured average global temperatures.

      LBNL got it wrong. If CO2 had an effect on climate, it would have appeared in the 542 million years of the Phanerozoic eon.

      You need to take the blinders off.

    2. This comment has been removed by the author.

    3. If existing albedo is 30% (of which two-thirds is due to cloud cover) then an increase to just 31% lowers the temperature by about 0.9 degree which is comparable with the difference between the Maunder minimum temperatures and the mean temperatures over the last ~1000 year natural cycle. The radiation from CO2 (like that from water vapor) works against the gravitationally-induced temperature gradient, thus leading to lower surface temperatures because, in case you don't know, the so-called wet adiabatic lapse rate is less in magnitude than the dry rate. How could you possibly maintain radiative balance with the Sun if the thermal profile did not rotate about a pivoting altitude and become lower at the surface end? You can't raise the surface temperature with more water vapor and, at the same time, have a less steep temperature gradient, because radiative balance would be way out and soon correct itself by lowering the whole thermal plot whilst retaining the same gradient.

      This is the epitome of self-contradiction in the greenhouse radiative forcing conjecture. A higher surface temperature due to water vapor and yet a less steep temperature gradient !?!?!

    4. Determination of Standard Atmosphere, first done in the 1950s, is described here and agrees with some of what you say. The population gradient of ghg contributes substantially to the radiation balance. Energy balance is calculated at

  22. Dan, great work, you need to get this published in a refereed journal. It needs an abstract before the intro, and the intro may benefit from a little wordsmithing, otherwise it's excellent.

    I had difficulty understanding why the value of R squared was different in Table 1 between the two rows in yellow? Looks all the parameters are the same but the correlation coefficient is different, why is that?

    1. Measured values are 5-year smoothed at second hi-lite.

  23. Great to see this post found really understanding about global warming, I get another blog related to Global Weirding News and Global Warming Newswhich have almost every news links for global warming.

  24. I agree that climate change is all 100% natural, with possible causes relating to planetary orbits (see and their effects which could pertain to variations in solar activity, sunspots, cosmic ray intensity, cloud formation and solar intensity varying because of variation in the distance from the Sun. I’ve said all that on my first website years ago.

    What I don’t agree with is the radiative forcing conjecture that radiation from the cold atmosphere helps the Sun to raise the surface temperature to observed values. You cannot add the two and use the total in Stefan-Boltzmann calculations, but Hansen and Pierrehumbert thought they could get away with doing so: after all, politicians are ignorant about such matters and so are brain-washed climatologists.

    Only maximum entropy production can do that in the process wherein new thermal energy absorbed in the upper troposphere makes its way to the surface by natural convective heat transfer which is restoring thermodynamic equilibrium or, in other words, maximizing entropy. On what point do you have difficulty in your understanding of this standard physics based on the Second Law?

    1. I don't know what Hansen and Pierrehumbert did but "add the two" is certainly wrong. Some other mistakes by the 'consensus' are discussed at

      It appears you could use a refresher in radiation heat transfer to understand what is going on. Your last paragraph is nonsense. Heated air rises.

  25. Why we need the heat creep process:
    (copied from the comment thread on my blog)

    We know that, apart from direct solar radiation, there must be additional thermal energy entering a planets surface each morning and early afternoon because otherwise we dont have any warming at those times, and we dont have overall energy balance at the surface for each complete rotation, such as each 24 hours for Earth. In typical Earth energy diagrams like the one on this page, they show back radiation of about 324W/m^2 into the surface, but the electro-magnetic energy in that radiation is not converted to thermal energy in the warmer surface, and so there is no heat transfer from cooler to warmer regions. Such heat transfers can never happen by radiation. Instead, a similar amount of thermal energy is actually delivered into the surface by conduction across the interface with the atmosphere. This was energy originally absorbed each morning from new solar radiation that is strong enough to raise the cold temperatures found mostly in the upper troposphere and above. This new energy must be what makes its way to the surface, but it does not do so by radiation. Only non-radiative processes (natural convective heat transfers) can transfer thermal energy up the sloping thermal plane, provided that the process is increasing entropy. This is because such processes depend upon energy being transferred by molecular collisions, and gravity acts on molecules in flight between collisions increasing their kinetic energy as their gravitational potential energy decreases. But, by definition, no energy is transferred across internal boundaries in a perfect state of thermodynamic equilibrium. Energy transfers only happen when new energy creates a state that is not in thermodynamic equilibrium, so that entropy can increase. If the new energy is at the top of a column of the troposphere, then a new state of thermodynamic equilibrium (with its associated temperature gradient) can only be attained if some of that new thermal energy at the top is dispersed downwards to warmer regions. That is what I have been first to explain does in fact happen, so that we know this heat creep process is possible, and it is indeed the correct physics which explains the missing energy that raises the surface temperature each morning, even under thick cloud cover.

    1. You have really talked yourself into a mess here. The clear atmosphere is essentially transparent to radiation from the sun. Clouds absorb some, reflect some, and transmit some. Start there or just look at