Stockwell asks: Is the Atmosphere Still Warming?

Guest post by Dr. David Stockwell

I suspect that the only really convincing evidence against global warming is a sustained period of no global warming or cooling — climate sensitivity and feedbacks are too esoteric.

I have followed the recent global temperature with some excitement, and started to prepare a follow up to a previous article I wrote on the failure of global temperature to meet AGW expectations.

The Nature publication “Recent Climate Observations Compared to Projections” by Rahmstorf, Hansen and others in 2007 claimed an up-tick in a graph showed that “global temperatures were increasing faster than expected”, and consequently climate change would be worse than expected. In “Recent Climate Observations: Disagreement with Projections”, using their methodology and two additional year’s data, the up-tick was shown to be an artefact of inadequate smoothing of the effects of a strong El Nino. Perhaps this rebuttal played some part in subsequent revisions of Rahmstorf’s graph with longer smoothing, which had the unfortunate effect (for him) of removing the up-tick, so they could no longer claim, “global temperatures were increasing faster than expected”.

Can we answer the question “Is the Atmosphere still warming” in a reasonable way?

From the field of econometrics comes empirical fluctuation processes (EFP), available to programmers in an R package called strucchange – developed to analyse such things as changes in exchange rates by the brilliant Achim Zeileis. The idea is to find a test of the null hypothesis that the slope parameter m for a section of a series has not changed over time:

H0: m1 = m2 versus the alternative H1: m1 not equal to m2

The idea is to move a window of constant width over the whole sample period, and compare local trends with the overall distribution of trends. The resulting process should not fluctuate (deviate from zero) too much under the null hypothesis and—as the asymptotic distributions of these processes are well-known—boundaries can be computed, which are only crossed with certain probability. If, on the other hand, the empirical process shows large fluctuations and crosses the boundary, there is evidence that the data contains a structural change in the parameter. The peaks can be dated and segmented regression lines fit between the breaks in slope.

I applied the strucchange function EFP to the five official global temperature data sets (CRU, GISS, NOAA, UAH and RSS) from 1978 using the latest values in 2011, and to mean global sea level. The results for the global temperature are below:

Click to enlarge - Figure 1. Fluctuation process, structural change model and information measures determining the number of structural breaks for the five global temperature data-sets (CRU, GISS, NOAA, UAH and RSS).

The fluctuation process (top panel) crosses the upper significance boundary a number of times, indicating that the trend parameter is unstable. For example, it crosses in 1998, coincident with the strong El Nino, and then relaxes. Most recently, three of the five data sets are at the lower boundary, indicating that at least the CRU, NOAA and RSS datasets have shifted away from the overall warming trend since 1978.

The middle panel shows the structural break model for the CRU data, with the optimal number of breaks given by the minimum of the Bayesian Information Criterion (BIC) (bottom panel). The locations of the breaks are coincident (with a lag) with major events: the ultra-Plinian (stratosphere reaching) eruptions of Mt Chichon and Mt Pinatubo, the Super El Nino and the Pacific Decadal Oscillation (PDO) phase change in 2005.

Sometimes these types of models are sensitive to the start and end point, so I re-ran the analysis with data from 1950. Figure 2 is the resulting structural break model for CRU. While the fluctuation process did not show the same degree of recent downtrend, the structural break model is similar to the shorter series in Figure 1, except the temperatures since 1998 are fit with a single flat segment.

The temperature is plotted over random multiple AR(1) simulations, showing the temperature has ranged between the extremes of an AR(1) model over the period.

click to enlarge Figure 2. Linear vs. segmented regressions for the global temperature dataset CRU, with the timing of significant climatic events.

Another indication of global temperature is the mean global sea level, both barometric and non-barometric adjusted. Global sea levels tell the same story as atmospheric temperature, with a significant deceleration in sea level rise around the PDO shift in 2005.

click to enlarge - Figure 3. The fluctuation process, structural break model and information measures for global mean sea level, both barometric and non-barometric adjusted.

By these objective criteria, there does appear to be a structural change away from the medium-term warming trend. Does this mean global warming has stopped?

What are the arguments that warming continues unabated?

Easterling and Wehner in their article “Is the climate warming or cooling?” lambasted “Numerous websites, blogs and articles in the media [that] have claimed that the climate is no longer warming, and is now cooling” for “cherry picking” the recent data. They examined the distribution of 10 year slopes of both the realized and modelled global temperature. They argued that because there were a small number of periods of flat 10 year temperatures that the long-term warming trend is intact.

Both E&W and EFP agree that there is a small chance of flat temperatures for 10 years (EFP says around 5%) during a longer-term warming trend. What E&W’s are saying is that given a small chance at one time, the chance of flat temperatures at any time, over the last 50 years say, is much higher. This doesn’t alter the fact that to an observer during any of those decades when temperature was flat (as now) there would still be a 5% chance of a break in the long-term trend.

Breusch and Vahid (2008 updated in 2011) chimed in with “Global Temperature Trends”, stating “there is no significant evidence for a break in trend in the late 1990s”, and “There is nothing to suggest that anything remarkable has happened around 1998.” As hard as I looked I could not find any estimates of significance to back up their claim of significant evidence.

The statement is even more puzzling as the last 15% at the ends of the series are typically not tested for breaks due to low power of the test on the diminishing numbers of data. The 1990’s fall in the outside 15%. Breaks the size of the break in 1976 would not have been detected on their data.

Of course, there are a variety of other observations of the Earth’s radiative balance and ocean heat content, supporting of the “no warming” claim, by top researchers such as Douglass and Loehle. There does not appear to be any credible empirical evidence from the AGW camp that the atmosphere is still warming.

I suspect that as in “Recent Climate Observations” where climate scientists were fooled into thinking that “climate change will be worse than expected” by the steep up-tick in global temperatures during a strong El Nino, they have also been fooled by a steep but longer-term up-tick in global temperatures associated with a positive phase of the PDO.

0 0 votes
Article Rating

Discover more from Watts Up With That?

Subscribe to get the latest posts sent to your email.

98 Comments
Inline Feedbacks
View all comments
April 13, 2011 6:57 pm

Elmer: The sea level in in mm. There is a replacement graph with the correct axis label coming.
BarvoZulu; Not sure who you are addressing, but the idea is to approach the question in a principled way, not an opinionated way. The EFP indicates that we are getting very close to 95% confidence that the atmosphere isn’t warming.

P. Solar
April 13, 2011 6:58 pm

It is unclear from your text and the graphs how you get from the EFP plot to your linear breakpoints. It would appear from the paper on strucchange that dates where the EFP breaks the boundary indicate a change in structure.
Now looking at fig. 1 , the most obvious dates are 1988 and 1997.5 , you chose to break at 1999 and 1988 falls exactly half way between two of your break points.
It appears from what is presented here that the EFP test and BIC has been used to suggest 4 break points but that the break points that were chosen were not derived from the results of the EFP plot.
Maybe you need to explain how you chose the break-points and how (of indeed if) these relate to the efp results.

P. Solar
April 13, 2011 7:10 pm

David Stockwell says:
April 13, 2011 at 6:08 pm

P Solar, the gray scribbles are: “The temperature is plotted over random multiple AR(1) simulations, showing the temperature has ranged between the extremes of an AR(1) model over the period.”

What AR model is that? How was it chosen , what does it mean? Why is it relevant?
Just showing a scribble plot “model” that has about the same amplitude does not prove/show anything unless you explain it.

The uncertainty of the breaks are indicated by the red bars over the x axis of the middle panels.

No, I meant the experimental uncertainty of the data. It’s an old-fashioned concept that used to exist in the hard sciences.
What you show with the red bars is the variable time lag between your apparently arbitrary breakpoints and the “events” you highlight.
Thanks for the reply.

April 13, 2011 7:11 pm

Ric Werme says:
April 13, 2011 at 5:58 pm
http://www.facebook.com/event.php?eid=148669818520369
According to their Facebook page, 3,116 FB members will be attending, 2,426 may be attending. I clicked “See all” and it only showed me a few. Trying to friend them all would take a while.
=========================================
If you find one that wants to play, let the rest of us WUWT FB members know, wouldja?

Katherine
April 13, 2011 7:21 pm

Sometimes these types of models are sensitive to the start and end point, so I re-ran the analysis with data from 1950.
If you re-run the analysis with data starting from 1930 or 1920 would the warming be similar?

April 13, 2011 7:32 pm

I like this cross-diciplinary approach / analysis.
Key observations :
1) The current segment is way flatter than any of the other segments
2) It quantifies what is visually obvious
Very intriguing !

April 13, 2011 7:36 pm

P Solar: The Zeileis approach, that I am trying to emulate, uses the EFP to determine if the slope parameter is unstable. If it is, by exceeding the 95%CL, then you proceed with the exploratory segmentation technique to get a model of how the trends might have been changing.
The EFP is not used for the dating of the breaks, which is just a brute force search based on the change in residuals for break vs no break.
EFP says it is not correct to fit a linear regression to the temperature — because it appears to be a process with structural changes, much like an exchange rate might change due to some government intervention.

April 13, 2011 7:42 pm

Alan+Katherine on time frame: “I’m sorry but could we stop “analyzing” the temperatures within such a narrow time frame window?…”
Its OK when you ask a specific question. The change in slope from 2000 to 2010 is a significant departure from the steep slope from 1978 to 2010. The change in slope from 2000 to 2010 is not significant relative to the lower slope from 1950, or from 1900 for that matter.
That is why I say that there appears to be a significant departure from the 1978 slope, and of course, this is the fist that is supposed to be mostly caused by CO2 according to the IPCC. But CO2 is still increasing, but temperature is not, why?

Luther Wu
April 13, 2011 7:59 pm
tommoriarty
April 13, 2011 8:00 pm

“…subsequent revisions of Rahmstorf’s graph with longer smoothing, which had the unfortunate effect (for him) of removing the up-tick…”
Rahmstorf must be getting used to the idea of updated data taking all the fun out of his predictions. See, for example…
http://climatesanity.wordpress.com/2010/11/17/rahmstorf-2009-part-9-applying-three-corrections/

April 13, 2011 8:06 pm

David Stockwell says:
“…I wrote on the failure of global temperature to meet AGW expectations.”
As Prof Richard Feynman explains regarding any new hypothesis: “We compare it directly with observation to see if it works.” If the hypothesis doesn’t agree with experiment, if it doesn’t agree with observation, it’s wrong.
The AGW hypothesis [and especially the Catastrophic AGW hypothesis] do not agree with experiment or observation. Therefore, they are both wrong.
The scientific method is ignored by the purveyors of the failed CAGW Hypothesis Conjecture. Their model is wrong, because it does not agree with observation. There is no runaway global warming, nor even warming. Despite a hefty ≈40% increase in CO2, temperatures are declining. Thus, the CO2=CAGW [and CO2=AGW] hypotheses are falsified; they do not agree with observation.
Who are you gonna believe? Michael Mann? Or Richard Feynman?

April 13, 2011 8:10 pm

Smokey: That’s my view. If you haven’t had your favourite theories demolished a few dozen times, you are not a scientist.

mike g
April 13, 2011 8:13 pm


Well, Alan, then I’ve got just the videos for you. They’re of Bob Carter doing exactly what you are asking for. And, doing it very convincingly. No wonder media down under has him in the cross hairs.

And you can get the other three parts from Lucy’s fine site:
http://www.greenworldtrust.org.uk/Science/Curious.htm

April 13, 2011 8:16 pm

P Solar on AR(1): This takes some explaining and not essential, but the literature on such tests (read B&V linked in the article) compares the realized temps to a presumed model of random variation. If you take the AR(1) model of the detrended temps, (ie Y_t= a Y_t-1) + e) then the a coefficient is about 0.7 and the SD of e is about 0.1C. Simulate that tons of times to get a feel for the spread of results due to pure, autocorrelated randomness. Thats the gray area in fig 5.
Most of the literature hovers around the 95%CL for tests of the realized temperature against random AR models (of various flavours).

April 13, 2011 8:29 pm

Thanks, Dr. Stockwell. I appreciate your honesty.
Have you seen this knocking of heads between two environmental scientists on global warming? – http://college-ethics.blogspot.com/2011/04/two-environmentalists-knock-heads.html
I’d be interested in you thoughts on this.

April 13, 2011 8:40 pm

Alice: I am more concrete about that sort of thing. The most basic assumption of linear regression is that the slope parameter m in y=mx+c is constant. It’s not. So the model is wrong. I fit a model with segments for each of the different slopes. Temperature is not increasing anymore. Either it increased from 1978 because of CO2 but has ‘maxed out’, or it was caused by the sun all along.

Lady Life Grows
April 13, 2011 8:45 pm

David Stockwell says:
April 13, 2011 at 8:10 pm
“If you haven’t had your favourite theories demolished a few dozen times, you are not a scientist.”
I had a favorite T-shirt I got from a roommate that I called my “science T-shirt.” It showed a honey dripper and said “oh, Lord, make my words as sweet as honey, for tomorrow I may have to eat them.”
But I don’t feel sweet as honey when I consider that practical consequences of the AGW hysteria. Aside from fighting photosynthesis and risking extinctions, the econazis are burning so much corn as ethanol that world food stocks have declined and the Arabs countries are rioting. Tens to hundreds of thousands have been killed so far, and it may touch off larger wars that could kill millions or even billions of people.
That program is so inefficient that it would not reduce warming even if the hysteria were right. So maybe we can get the US congress to have mercy on the world’s poor and stop ethanol subsidies. We need to cut plenty anyway, and that is one good place.

P. Solar
April 13, 2011 9:14 pm

David Stockwell says:
April 13, 2011 at 7:36 pm

P Solar: The Zeileis approach, that I am trying to emulate, uses the EFP to determine if the slope parameter is unstable. If it is, by exceeding the 95%CL, then you proceed with the exploratory segmentation technique to get a model of how the trends might have been changing.
The EFP is not used for the dating of the breaks, which is just a brute force search based on the change in residuals for break vs no break.
EFP says it is not correct to fit a linear regression to the temperature — because it appears to be a process with structural changes, much like an exchange rate might change due to some government intervention.

EFP says it is not correct to fit a linear regression to the temperature … but then you do. Including those areas that were found to be the most unstable , right across 1988 for example 😕
In fact it’s hard to see what use you make of efp once you’ve plotted it.
So the key work you do here is the “brute force” optimisation of segments, a process you don’t describe.
What is the rationale for the discontinuities you introduce ? What happens if you go to 8 or 10 segments, your optimisation does not seem to look beyond 5. That may well just be a local minimum, and a pretty shallow one at that.
I agree with the idea that post-2003 has been pretty flat but you don’t need a complex analysis to see that.
The rest of it just does not ring true for me. I don’t think this kind of method will be accepted as proving anything , apart from by some of the less critical minds here.

Bob K.
April 13, 2011 9:43 pm

I am sure this: “H0: m1 = m2 VERSES the alternative H1: m1 not equal to m2” was ment to read: “H0: m1 = m2 VERSUS the alternative H1: m1 not equal to m2”.

April 13, 2011 9:49 pm

Bob: “In looking at you Figure 1, the PDO switch that you’ve highlighted is the only event that lags the structural break.”
I wouldn’t take the location of the PDO shift as exact. When I say PDO I mean all those processes that PDO represents, that seem to explain the multi-decadal temperature fluctuations so well.

April 13, 2011 10:05 pm

P Solar: “EFP says it is not correct to fit a linear regression to the temperature … but then you do. ”
EFP says it is not correct to fit a SINGLE linear regression to the temperature. So you need a model with more than one slope. Which is what the method does.
The EFP justifies the use of more than one slope. Its a critical step. The next step is dating of the breaks and is more exploratory — there are other criteria than BIC. But BIC should be a global optimum — it would always be greater with more parameters I think, unless it is a really weird type of series (I could imagine one I suppose if it was highly periodic).
The rational for a break, in the case of one test of breaks, is the supremum ration of the sum of squares for a single line vs. a broken line, and is an F statistic: F = RSS1/RSS2. However, this is not the only statistic either.
You really want to go to the literature if you want the details as its technical.

RoHa
April 13, 2011 10:22 pm

“Stockwell asks: Is the Atmosphere Still Warming?”
And is his answer “No”?

Katherine
April 13, 2011 10:47 pm

David Stockwell says:
That is why I say that there appears to be a significant departure from the 1978 slope, and of course, this is the fist that is supposed to be mostly caused by CO2 according to the IPCC. But CO2 is still increasing, but temperature is not, why?
Thanks for the clarification, Dr. Stockwell. I think I understand your article better now.

April 13, 2011 10:53 pm

“I agree with the idea that post-2003 has been pretty flat but you don’t need a complex analysis to see that. The rest of it just does not ring true for me. I don’t think this kind of method will be accepted as proving anything , apart from by some of the less critical minds here.”
When I hear the words “pretty flat” or “pretty close” I immediately think “He must be a climate scientist”. As to whether it rings true for you, or me — who cares? The case for a segmented fit is this:
1. The model y=mx+c is WRONG as m is not a constant. Therefore you must use a model that allows m to fluctuate (and this goes for AR(1) with a constant drift term too).
2. The segmentation method “discovers” exogenous causes such as volcanic eruptions, without any knowledge of the timing of these events. So the method is interesting as it tells you something that you didn’t assume.
Now it could be that including TSI as a dependent variable explains enough of the variation that the EFP no longer crosses the boundary, then the case for a segmented model would not be so strong — I will have to look into that.
3. The method provides a way to evaluate the trends at the ends of a series without the problems that moving averages have. This is an important issue as we are most interested in what is happening at the ends, as that is the most recent data.

rogerthesurf
April 13, 2011 11:35 pm

Does anyone know, according to the IPCC, how soon will it be before we are fully effected by the 7 meter sea level rise caused by the melting of the Greenland ice cap then?
Cheers
Roger
http://www.rogerfromnewzealand.wordpress.com