The Eruption Over the IPCC AR5

Guest Post by Willis Eschenbach

In the leaked version of the upcoming United Nations Intergovernmental Panel on Climate Change (UN IPCC) Fifth Assessment Report (AR5) Chapter 1, we find the following claims regarding volcanoes.

The forcing from stratospheric volcanic aerosols can have a large impact on the climate for some years after volcanic eruptions. Several small eruptions have caused an RF for the years 2008−2011 of −0.10  [–0.13 to –0.07] W m–2, approximately double the 1999−2002 volcanic aerosol RF.

and

The observed reduction in warming trend over the period 1998–2012 as compared to the period 1951–2012, is due in roughly equal measure to a cooling contribution from internal variability and a reduced 2 trend in radiative forcing (medium confidence). The reduced trend in radiative forcing is primarily due 3 to volcanic eruptions and the downward phase of the current solar cycle.

Now, before I discuss these claims about volcanoes, let me remind folks that regarding the climate, I’m neither a skeptic nor am I a warmist.

I am a climate heretic. I say that the current climate paradigm, that forcing determines temperature, is incorrect. I hold that changes in forcing only marginally and briefly affect the temperature. Instead, I say that a host of emergent thermostatic phenomena act

quickly to cool the planet when it is too warm, and to warm it when it is too cool.

One of the corollaries of this position is that the effects of volcanic eruptions on global climate will be very, very small. Although I’ve demonstrated this before, Anthony recently pointed me to an updated volcanic forcing database, by Sato et al. Figure 1 shows the amount of forcing from the historical volcanoes.

volcanic forcing 1850 2012 Sato

Figure 1. Monthly changes in radiative forcing (downwelling radiation) resulting from historical volcanic eruptions. The two large recent spikes are from El Chichon (1983) and Pinatubo (1992) eruptions. You can see the average forcing of -0.1 W/m2 from 2008-2011 mentioned by the IPCC above. These are the equilibrium forcings Fe, and not the instantaneous forcing Fi.

Note that the forcings are negative, because the eruptions inject reflective aerosols into the stratosphere. These aerosols reflect the sunlight, and the forcing is reduced. So the question is … do these fairly large known volcanic forcings actually have any effect on the global surface air temperature, and if so how much?

To answer the question, we can use linear regression to calculate the actual effect of the changes in forcing on the temperature. Figure 2 shows the HadCRUT4 monthly global surface average air temperature.

hadCRUT4 1850-2012 and gaussianFigure 2. Monthly surface air temperatures anomalies, from the HadCRUT4 dataset. The purple line shows a centered Gaussian average with a full width at half maximum (FWHM) of 8 years.

One problem with doing this particular linear regression is that the volcanic forcing is approximately trendless, while the temperature has risen overall. We are interested in the short-term (within four years or so) changes in temperature due to the volcanoes. So what we can do to get rid of the long-term trend is to only consider the temperature variations around the average for that historical time. To do that, we subtract the Gaussian average from the actual data, leaving what are called the “residuals”:

residual hadcrut4 monthly anomaliesFigure 3. Residual anomalies, after subtracting out the centered 8-year FWHM gaussian average.

As you can see, these residuals still contain all of the short-term variations, including whatever the volcanoes might or might not have done to the temperature. And as you can also see, there is little sign of the claimed cooling from the eruptions. There is certainly no obvious sign of even the largest eruptions. To verify that, here is the same temperature data overlaid on the volcanic forcing. Note the different scales on the two sides.

residual hadcrut4 monthly anomalies plus forcingFigure 4. Volcanic forcing (red), with the HadCRUT4 temperature residual overlaid.

While some volcanoes line up with temperature changes, some show increases after the eruptions. In addition, the largest eruptions don’t seem correlated with proportionately large drops in temperatures.

So now we can start looking at how much the volcanic forcing is actually affecting the temperature. The raw linear regression yields the following results.

R^2 = 0.01 (a measure from zero to one of how much effect the volcanoes have on temperature)

"p" value of R^2 = 0.03 (a measure from zero to one how likely it is that the results occurred by chance) (adjusted for autocorrelation).

Trend = 0.04°C per W/m2, OR 0.13°C per doubling of CO2 (how much the temperature varies with the volcanic forcing)

"p" value of the TREND = 0.02 (a measure from zero to one how likely it is that the results occurred by chance) (adjusted for autocorrelation).

So … what does that mean? Well, it’s a most interesting and unusual result. It strongly confirms a very tiny effect. I don’t encounter that very often in climate science. It simultaneously says that yes, volcanoes do affect the temperature … and yet, the effect is vanishingly small—only about a tenth of a degree per doubling of CO2.

Can we improve on that result? Yes, although not a whole lot. As our estimate improves, we’d expect a better R^2 and a larger trend. To do this, we note that we wouldn’t expect to find an instantaneous effect from the eruptions. It takes time for the land and ocean to heat and cool. So we’d expect a lagged effect. To investigate that, we can calculate the R^2 for a variety of time lags. I usually include negative lags as well to make sure I’m looking at a real phenomenon. Here’s the result:

rsquared forcing and temperatureFigure 5. Analysis of the effects of lagging the results of the volcanic forcing. 

That’s a lovely result, sharply peaked. It shows that as expected, after a volcano, it takes about seven-eight months for the maximum effects to be felt.

Including the lag, of course, gives us new results for the linear regress, viz:

R^2 = 0.03 [previously 0.01]

"p" value of R^2 = 0.02 (adjusted for autocorrelation) [previously 0.03]

Trend  = 0.05°C per W/m2, OR 0.18 ± 0.02°C per doubling of CO2 [previously 0.13°C/doubling]

"p" value of the Trend = 0.001 (adjusted for autocorrelation). [previously 0.02]

As expected, both the R^2 and the trend have increased. In addition the p-values have improved, particularly for the trend. At the end of the day, what we have is a calculated climate sensitivity (change in temperature with forcing) which is only about two-tenths of a degree per doubling of CO2.

Here are the conclusions that I can draw from this analysis.

1) The effect of volcanic eruptions is far smaller than generally assumed. Even the largest volcanoes make only a small difference in the temperature. This agrees with my eight previous analyses (see list in the Notes). For those who have questions about this current analysis, let me suggest that you read through all of my previous analyses, as this is far from my only evidence that volcanoes have very little effect on temperature.

2) As Figure 5 shows, the delay in the effects of the temperature is on the order of seven or eight months from the eruption. This is verified by a complete lagged analysis (see the Notes below). That analysis also gives the same value for the climate sensitivity, about two tenths of a degree per doubling.

3) However, this is not the whole story. The reason that the temperature change after an eruption is so small is that the effect is quickly neutralized by the homeostatic nature of the climate.

Finally, to return to the question of the IPCC Fifth Assessment Report, it says:

There is very high confidence that models reproduce the more rapid warming in the second half of the 20th century, and the cooling immediately following large volcanic eruptions.

Since there is almost no cooling that follows large volcanic eruptions … whatever the models are doing, they’re doing it wrong. You can clearly see the volcanic eruptions in the model results … but you can’t see them at all in the actual data.

The amazing thing to me is that this urban legend about volcanoes having some big effect on the global average temperature is so hard to kill. I’ve analyzed it from a host of directions, and I can’t find any substance there at all … but it is widely believed.

I ascribe this to an oddity of the climate control system … it’s invisible. For example, I’ve shown that the time of onset of tropical clouds has a huge effect on incoming solar radiation, with a change of about ten minutes in onset time being enough to counteract a doubling of CO2. But no one would ever notice such a small change.

So we can see the cooling effect of the volcanoes where it is occurring … but what we can’t see is the response of the rest of the climate system to that cooling. And so, the myth of the volcanic fingerprints stays alive, despite lots of evidence that while they have large local effects, their global effect is trivially small.

Best to all,

w.

PS—The IPCC claims that the explanation for the “pause” in warming is half due to “natural variations”, a quarter is solar, and a quarter is from volcanoes. Here’s the truly bizarre part. In the last couple decades, using round numbers, the IPCC predicted about 0.4°C of warming … which hasn’t happened. So if a quarter of that (0.1°C) is volcanoes, and the recent volcanic forcing is (by their own numbers) about 0.1 W/m2, they’re saying that the climate sensitivity is 3.7° per doubling of CO2.

Of course, if that were the case we’d have seen a drop of about 3°C from Pinatubo … and I fear that I don’t see that in the records.

They just throw out these claims … but they don’t run the numbers, and they don’t think them through to the end.

Notes and Data

For the value of the forcing, I have not used the instantaneous value of the volcanic forcing, which is called “Fi“. Instead, I’ve used the effective forcing “Fe“, which is the value of the forcing after the system has completely adjusted to the changes. As you might expect, Fi is larger than Fe. See the spreadsheet containing the data for the details.

As a result, what I have calculated here is NOT the transient climate response (TCR). It is the equilibrium climate sensitivity (ECS).

For confirmation, the same result is obtained by first using the instantaneous forcing Fi to calculate the TCR, and then using the TCR to calculate the ECS.

Further confirmation comes from doing a full interative lagged analysis (not shown), using the formula for a lagged linear relationship, viz:

T2 = T1 + lambda (F2 – F1) (1 – exp(-1/tau)) + exp(-1/tau) (T1 – T0)

where T is temperature, F is forcing, lambda is the proportionality coefficient, and tau is the time constant.

That analysis gives the same result for the trend, 0.18°C/doubling of CO2. The time constant tau was also quite similar, with the best fit at 6.4 months lag between forcing and response.

<strong>

In this case it’s the Sato paper, which provides a dataset of optical thicknesses “tau”, and says:

The relation between the optical thickness and the forcings are roughly (See “Efficacy …” below):

instantaneous forcing Fi (W/m2) = -27 τ

adjusted forcing Fa (W/m2) = -25 τ

SST-fixed forcing Fs (W/m2) = -26 τ

effective forcing Fe (W/m2) = -23 τ

And “Efficacy” refers to

Hansen, J., M. Sato, R. Ruedy, L. Nazarenko, A. Lacis, G.A. Schmidt, G. Russell, et al. 2005. Efficacy of climate forcings. J. Geophys. Res., 110, D18104, doi:10.1029/2005/JD005776.

Forcing Data

For details on the volcanic forcings used, see the Sato paper, which provides a dataset of optical thicknesses “tau”, and says:

The relation between the optical thickness and the forcings are roughly (See “Efficacy …” below):

instantaneous forcing Fi (W/m2) = -27 τ

adjusted forcing Fa (W/m2) = -25 τ

SST-fixed forcing Fs (W/m2) = -26 τ

effective forcing Fe (W/m2) = -23 τ

And “Efficacy” refers to

Hansen, J., M. Sato, R. Ruedy, L. Nazarenko, A. Lacis, G.A. Schmidt, G. Russell, et al. 2005. Efficacy of climate forcings. J. Geophys. Res., 110, D18104, doi:10.1029/2005/JD005776.

(Again, remember I’m using their methods, but I’m not claiming that their methods are correct.)

Future Analyses

My next scheme is that I want to gin up some kind of prototype governing system that mimics what it seems the climate system is doing. The issue is that to keep a lagged system on course, you need to have “overshoot”. This means that when the temperature goes below average, it then goes above average, and then finally returns to the prior value. Will I ever do the analysis? Depends on whether something shinier shows up before I get to it … I would love to have about a dozen bright enthusiastic graduate students to hand out this kind of analysis to.

I also want to repeat my analysis using “stacking” of the volcanoes, but using this new data, along with some mathematical method to choose the starting points for the stacking … which turns out to be a bit more difficult than I expected.

Previous posts on the effects of the volcano.

Prediction is hard, especially of the future. 

Pinatubo and the Albedo Thermostat

Missing the Missing Summer

Volcanic Disruptions

Dronning Maud Meets the Little Ice Age 

New Data, Old Claims about Volcanoes 

Volcanoes: Active, Inactive and Interactive

Stacked Volcanoes Falsify Models

Get notified when a new post is published.
Subscribe today!
0 0 votes
Article Rating
463 Comments
Inline Feedbacks
View all comments
September 22, 2013 5:23 pm

Forrest M. Mims III (September 22, 2013 at 1:33 pm) “The aerosols reduced the photocurrent from the solar cell at noon by about 5 percent during this time. These observations were associated with a reduction in temperature of about 2 degrees F, which is similar to other reports elsewhere”
What happened at night? My guess is that nights became warmer.

September 22, 2013 5:26 pm

sunshinehours1 says:
September 22, 2013 at 5:13 pm
There were prior volcanic eruptions, too.
Not summer temperatures in 1816, among lowest ever recorded in the CET series.
To quote yourself:
HADCET: http://www.metoffice.gov.uk/hadobs/hadcet/mly_cet_mean_sort.txt
out of 353/354
1816
Coldest July ever
20th coldest August
26th coldest June
38th coldest May
50th coldest April
34th coldest September
238th coldest October – the outlier
23rd coldest November
11th coldest December
1817
3rd coldest May
15th coldest July
10th coldest August
2nd coldest October.

September 22, 2013 5:39 pm

Steven Mosher (September 22, 2013 at 4:19 pm) “Imagine if I tried to calculate the sensitivity from the drop in termperature when the sun goes down.”
Imagine if you tried calculating sensitivity without the correct daily rise in temperature over the ocean that creates the convection that cools the planet (i.e. one of the thermostats being referred to here). That’s what climate models do, see fig 3: http://echorock.cgd.ucar.edu/cas/adai/papers/DaiTrenberth_JC04.pdf

kuhnkat
September 22, 2013 5:39 pm

Willis,
try it again using RAW Data if you can find it. You may be a victim of the reduction of historic temps.

September 22, 2013 5:46 pm

Hmm, my guess might actually be correct: “it is shown that the dominant decreasing trend of mean maximum temperature and the dominant increasing trend of mean minimum temperature over periods 1992–1994 and 1985–1987 relative to that over the period 1988–1990 are consistent with the distribution of stratospheric volcanic aerosols and predictions from aerosol radiative forcing in the southeastern U.S.”
From http://www.sciencedirect.com/science/article/pii/S1352231097002446

September 22, 2013 5:49 pm

Willis, if you look at my links you’ll find references to Douglas’ paper. Hmm arguably one of the first words, but not the last word in trying to estimate sensitivity from volcano data.

September 22, 2013 6:06 pm

“dp says:
September 22, 2013 at 4:24 pm
The sun doesn’t go down – it goes somewhere else. It is not an impulse when the sun goes somewhere else.”
The sensitivity parameter is a measure of the system response over a given time period.
Lets take a simple example. The speed of your car is D/T
Now, lets apply a force to your rear tires by stomping on the pedal. there will be an instaneous response, a transient response as you gain speed, and an equlibrium response when you hit top speed.
The ECS is what we are mostly interested in, the full response after all feedbacks,
Estimating that from a tap on the breaks requires a bit more work than Willis presents.
It’s an active area of research. What I can say is that most of the models cant replicate the effect perfectly they tend to over estimate the cooling and overestimate the rebound. This has been known for a while.
At first I surmised it had to do with sensitivity. But it doesnt.

pat
September 22, 2013 6:17 pm

an ambiguous ‘central’ in the headline; the Nobel Prize gets a mention of course, plus some classic quotes:
23 Sept: BBC: Matt McGrath: Global warming pause ‘central’ to IPCC climate report
Scientists will underline, with greater certainty than ever, the role of human activities in rising temperatures.
But many governments are demanding a clearer explanation of the slowdown in temperature increases since 1998.
One participant told BBC News that this pause will be a “central piece” of the summary…
In the latest draft summary, seen by the BBC, the level of scientific certainty has increased…
This slowdown, or hiatus as the IPCC refers to it, has been leapt upon by climate sceptics to argue that the scientific belief that emitting carbon dioxide into the atmosphere increases the temperature of the planet, is wrong.
Scientists have attempted to explain the pause in a number of ways, with many arguing that the Earth has continued to warm but that the heat has gone into oceans…
But there is no certainty and little agreement among scientists on the mechanisms involved…
Prof Arthur Petersen is the chief scientist at the Netherlands Environmental Assessment Agency and part of the Dutch delegation that will review the IPCC report.
“Governments are demanding a clear explanation of what are the possible causes of this factor,” he told BBC News.
“I expect that this will be a central piece of the summary.”…
Any changes to the text will need to be approved by the scientists, who will want to make sure that they are consistent with the underlying reports. This could lead to some tense moments.
“I wouldn’t say there is a reluctance of the authors to take up such an issue as the pause, but they want to do it in a proper way,” said Prof Petersen.
“There will remain a tension between how much you can deliver based on the peer-reviewed science and what the governments would like to have.”***
***In the wake of that year’s report, a small number of embarrassing errors were detected in the underlying material. The organisation’s reputation was also questioned in the Climategate rumpus.
“Overall, the message is, in that sense more conservative I expect, for this IPCC report compared to previous ones,” said Prof Petersen.
“The language has become more complicated to understand, but it is more precise.
“It is a major feat that we have been able to produce such a document which is such an adequate assessment of the science. That being said, it is virtually unreadable!”***
http://www.bbc.co.uk/news/science-environment-24173504

Chuck Nolan
September 22, 2013 6:19 pm

Stephen Wilde says:
September 22, 2013 at 11:39 am
One has to ask what sets the governing mechanism and in my view it must be atmospheric pressure which is a consequence of mass alone held within a gravitational field:
http://www.newclimatemodel.com/the-setting-and-maintaining-of-earths-equilibrium-temperature/
—————————————————-
Excellent read Stephen but does that mean temperature/climate/weather are somehow affected or controlled by earth’s gravity?
I believe the gravity vector is different throughout the world but it is also constantly changing.
Does that mean as gravity shifts due to earth quakes and core movement we get “climate change”?
Or would these changes be too small to really affect atmospheric pressure?
cn

September 22, 2013 6:25 pm

Willis Eschenbach says:
September 22, 2013 at 5:51 pm
I showed the data for yield & production & linked to those figures & the price numbers. Don’t know how you could have missed them. Your own FAO chart shows exactly what I posted, ie that wheat production fell from 1991 & didn’t recover until 1998. It’s you who ignore the data.
Your claim was that every crop of every type did not fall after Pinatubo. That is simply false. Wheat is the number two crop in the world by tonnage & number one by area planted. Yield & production went down after Pinatubo, & price went up. Sorry if you don’t like inconvenient truths, but that’s a fact.
That you don’t like this fact hardly means that I’ve gone over to the dark side. Quite the opposite. I rely on facts & data even when I might wish them to be otherwise.
I knew you were wrong about wheat. I didn’t check other crops, but since you were wrong about wheat, could be you’re wrong on the others as well, although my recollection is that corn went up. In science it’s best to avoid categorical statements, especially without checking.

gbaikie
September 22, 2013 6:38 pm

“The amazing thing to me is that this urban legend about volcanoes having some big effect on the global average temperature is so hard to kill. I’ve analyzed it from a host of directions, and I can’t find any substance there at all … but it is widely believed.”
One problem seems that you using air temperature.
So I don’t Earth is warmed by changes in air temperature.
So what would expect is for the oceans to gain less energy.
And I think people think volcano affect global average temperature because
affects the weather- you get cooler days. As in it causes it to snow in the summer.
Such weather may not show up in temperature records as significant.
It also seems that dust in atmosphere will have less affect at noon compared to when
sun is lower in the sky. So it could have less to do with daytime highs or nighttime
lows.
So tropical temperature could not affected as much as higher latitude regions.
So to measure affects one needs to careful measure ocean temperature in top
100 meters- I think it affects entire ocean, but top 100 meters would easier to measure.
Or measure long term affects of slightly cooler ocean.
And I think matters how big the eruption is. Pinatubo ejected about 5 cubic kilometers,
which quite different than ejecting 100 cubic km.

September 22, 2013 6:39 pm

Tom McCord says:
September 22, 2013 at 9:58 am
Is it possible that the “Year Without a Summer” in 1816 was not really caused by the eruption of Mount Tambora in Indonesia after all?
———————————————–
Looking at Tony Brown,s reconstruction of the CET chart 1816 appears to be the heart of the Dalton Minimum. There are 3 years together that form the low point. The first of which is about -1.2, then -1.7C, the last shows -1.5+C.

thingadonta
September 22, 2013 6:54 pm

It’s so easy being in climate alarmism. You can get rid of the missing meat through volcanos, the deep sea, or maybe the hot water at the bottom of the ocean is being subducted at ocean trenches (funding for deep sea missions?), or going into deep sea volcanos. Maybe the soil is hotter now?. The possibilities are just endless.

Greg Goodman
September 22, 2013 6:59 pm

Joe Born says:
September 22, 2013 at 4:35 pm
You may have an interesting point. But how about you post a graph produced by you code so that we can see what the effect is without every man jack of us having to do it just to see your point?
It seems you take a deliberately extreme case of tau=2x sample interval (which is fine to make the point) and imply that tau=5x sample may also be corrupted. Paul_K agreed with the basic formula , it don’t recall him commenting on what you say here.
One added word of caution is needed on the use of iterative response formulae like the one in question. They can take a very long time to spin up, ie converge to an accurate result from the initial conditions.
I looked into the possibility of using iterative filters to low-pass temperature time series. Often to get within even 5% accuracy it look more that half the length of the dataset. (Depends upon specifics.)

Philip Bradley
September 22, 2013 7:53 pm

Willis said,
The binomial distribution says that if you flip a coin seven times, you get five or more heads A QUARTER OF THE TIME!!! And that’s exactly the odds of five of seven volcanoes occurring in half of the year.
So you look at that result, which has NO STATISTiCAL SIGNIFICANCE AT ALL, and you build a whole theory about how volcanoes operate out of it …
Sorry to make an example out of your foolishness, Philip, it’s not personal. I just can’t tell you how tired I am of people claiming significance where none exists.

I don’t take this stuff personally. Although, I don’t why you blew up at me. Fig 5 says to me that volcanic cooling is limited by one or more sub-annual processes, and monsoons are an obvious candidate. If so, then time of year of the eruption will affect the lag to maximum cooling.
I wasn’t constructing a theory from 7 data points. If anything I was making a prediction based on what I know of monsoon processes. I’ll see if I can find some data that supports my prediction or not.

Greg Goodman
September 22, 2013 8:00 pm

Willis: “The amazing thing to me is that this urban legend about volcanoes having some big effect on the global average temperature is so hard to kill. ”
This was part of the cold war government by fear strategy. Nuclear winter was based of the same idea and made reference to the already accepted idea about volcanoes. The principal has been embedded for at least two generations. It is not going to be displaced because Willis said so on WUWT. It takes more than that to change an orthodoxy.
However, we all seem agreed that there is some effect. It’s just that the IPCC are spinning it to be many times larger in order balance their GHG forcings being many time larger than life too.
The problem is Chichon and Mt P coincided with cooling periods that were already under way, but anyone expecting to see a volcano effect will get bias confirmation and latch onto a false attribution.

1 4 5 6 7 8 19