Guest post by Mike Jonas
A few days ago, on Judith Curry’s excellent ClimateEtc blog, Vaughan Pratt wrote a post “Multidecadal climate to within a millikelvin” which provided the content and underlying spreadsheet calculations for a poster presentation at the AGU Fall Conference. I will refer to the work as “VPmK”.
VPmK was a stunningly unconvincing exercise in circular logic – a remarkably unscientific attempt to (presumably) provide support for the IPCC model[s] of climate – and should be retracted.
Background
The background to VPmK was outlined as “Global warming of some kind is clearly visible in HadCRUT3 [] for the three decades 1970-2000. However the three decades 1910-1940 show a similar rate of global warming. This can’t all be due to CO2 []“.
The aim of VPmK was to support the hypothesis that “multidecadal climate has only two significant components: the sawtooth, whatever its origins, and warming that can be accounted for 99.98% by the AHH law []“,
where
· the sawtooth is a collection of “all the so-called multidecadal ocean oscillations into one phenomenon“, and
· AHH law [Arrhenius-Hofmann-Hansen] is the logarithmic formula for CO2 radiative forcing with an oceanic heat sink delay.
The end result of VPmK was shown in the following graph
Fig.1 – VPmK end result.
where
· MUL is multidecadal climate (ie, global temperature),
· SAW is the sawtooth,
· AGW is the AHH law, and
· MRES is the residue MUL-SAW-AGW.
Millikelvins
As you can see, and as stated in VPmK’s title, the residue was just a few millikelvins over the whole of the period. The smoothness of the residue, but not its absolute value, was entirely due to three box filters being used to remove all of the “22-year and 11-year solar cycles and all faster phenomena“.
If the aim of VPmK is to provide support for the IPCC model of climate, naturally it would remove all of those things that the IPCC model cannot handle. Regardless, the astonishing level of claimed accuracy shows that the result is almost certainly worthless – it is, after all, about climate.
The process
What VPmK does is to take AGW as a given from the IPCC model – complete with the so-called “positive feedbacks” which for the purpose of VPmK are assumed to bear a simple linear relationship to the underlying formula for CO2 itself.
VPmK then takes the difference (the “sawtooth”) between MUL and AGW, and fits four sinewaves to it (there is provision in the spreadsheet for five, but only four were needed). Thanks to the box filters, a good fit was obtained.
Given that four parameters can fit an elephant (great link!), absolutely nothing has been achieved and it would be entirely reasonable to dismiss VPmK as completely worthless at this point. But, to be fair, we’ll look at the sawtooth (“The sinewaves”, below) and see if it could have a genuine climate meaning.
Note that in VPmK there is no attempt to find a climate meaning. The sawtooth which began life as “so-called multidecadal ocean oscillations” later becomes “whatever its origins“.
The sinewaves
The two main “sawtooth” sinewaves, SAW2 and SAW3, are:
Fig.2 – VPmK principal sawtooths.
(The y-axis is temperature). The other two sinewaves, SAW4 and SAW5 are much smaller, and just “mopping up” what divergence remains.
It is surely completely impossible to support the notion that the “multidecadal ocean oscillations” are reasonably represented to within a few millikelvins by these perfect sinewaves (even after the filtering). This is what the PDO and AMO really look like:
Fig.3 – PDO.
(link) There is apparently no PDO data before 1950, but some information here.
Fig.4 – AMO.
(link)
Both the PDO and AMO trended upwards from the 1970s until well into the 1990s. Neither sawtooth is even close. The sum of the sawtooths (SAW in Fig.1) flattens out over this period when it should mostly rise quite strongly. This shows that the sawtooths have been carefully manipulated to “reserve” the 1970-2000 temperature increase for AGW.
Fig.5 – How the sawtooth “reserved” the1980s and 90s warming for AGW.
Conclusion
VPmK aimed to show that “multidecadal climate has only two significant components”, AGW and something shaped like a sawtooth. But VPmK then simply assumed that AGW was a component, called the remainder the sawtooth, and had no clue as to what the sawtooth was but used some arbitrary sinewaves to represent it. VPmK then claimed to have shown that the climate was indeed made up of just these two components.
That is circular logic and appallingly unscientific. The poster presentation should be formally retracted.
[Blog commenter JCH claims that VPmK is described by AGU as “peer-reviewed”. If that is the case then retraction is important. VPmK should not be permitted to remain in any “peer-reviewed” literature.]
Footnotes:
1. Although VPmK was of so little value, nevertheless I would like to congratulate Vaughan Pratt for having the courage to provide all of the data and all of the calculations in a way that made it relatively easy to check them. If only this approach had been taken by other climate scientists from the start, virtually all of the heated and divisive climate debate could have been avoided.
2. I first approached Judith Curry, and asked her to give my analysis of Vaughan Pratt’s (“VP”) circular logic equal prominence to the original by accepting it as a ‘guest post’. She replied that it was sufficient for me to present it as a comment.
My feeling is that posts have much greater weight than comments, and that using only a comment would effectively let VP get away with a piece of absolute rubbish. Bear in mind that VPmK has been presented at the AGU Fall Conference, so it is already way ahead in public exposure anyway.
That is why this post now appears on WUWT instead of on ClimateEtc. (I have upgraded it a bit from the version sent to Judith Curry, but the essential argument is the same). There are many commenters on ClimateEtc who have been appalled by VPmK’s obvious errors. I do not claim that my effort here is in any way better than theirs, but my feeling is that someone has to get greater visibility for the errors and request retraction, and no-one else has yet done so.
Discover more from Watts Up With That?
Subscribe to get the latest posts sent to your email.
Stephen Rasey:
re your post at December 13, 2012 at 11:00 am.
Every now and then one comes across a pearl shining on the sand of WUWT comments. The pearls come in many forms.
Your post is a pearl. Its argument is clear, elegant and cogent. Thankyou.
Richard
Should it not be milliKelvin?
Steveta_uk says:
What VP has said repeatedly on JC’s site is basically that if you can provide a better fit, please do..
A “better fit” is not useful. In cases like this the correct model will not provide a better fit, because the correct model has deviations between theory and reality due to noise..
If I have a 100% normally distributed population and take 100 samples then the result will never be an exact normal distribution. I could model a “distribution” that better matched my samples, but it would most certainly not tell me anything useful. In fact it would lead me to believe my actual population was not normal.
This is why I am highly suspicious of any model that is trained on old data. It is basically an exercise in wiggle matching, not an exercise in getting the underlying physics correct. The best climate models will have pretty poor fit to old temperatures.
Substitute a 1300 year wave length sine with a max at MWP and min at LIA for the AGW function and you will get similar results.
lsvalgaard says:
December 13, 2012 at 11:21 am
vukcevic says:
December 13, 2012 at 11:11 am
Since you couldn’t fail my calculations
Of course, one cannot fail made-up ‘data’. What is wrong with your approach is to compute a new time series from two unrelated time series, and to call that ‘observed data’.
………………..
Wrong Doc.
Magnometer at the Tromso does it every single minute of the day and night.
http://www.vukcevic.talktalk.net/Tromso.htm
In red are Incoming variable solar magnetic field sitting on the top of the variable Earth’s magnetic field.
Rudolf Wolf started it with a compass needle, Gauss did it with a bit more sophisticated apparatus, and today numerous geomagnetic stations do it as you listed dozens in your paper on IDV.
So it is OK for Svalgaard of Stanford to derive IDV from changes of two combined magnetic fields, but is not for Vukcevic.
Reason plain and obvious, it would show that the SUN DOES IT !
Here is how geomagnetic field (Eart + solar) is measured and illustrated by our own Dr. Svalgaard
http://www.leif.org/research/Rudolf%20Wolf%20and%20the%20Sunspot%20Number.ppt#8
and he maintains they are not added together in his apparatus.
Can anyone spot 3 magnets?
Dr. S are you really serious to suggest that no changes in the Earth field are registered by your apparatus?
Case closed!
Steveta_uk says:
December 13, 2012 at 8:55 am
What VP has said repeatedly on JC’s site is basically that if you can provide a better fit, please do.
Nobody that I’ve seen has yet done so. Now I’m not in any way suggesting that what VP has done is in any way useful science. But still, can you not simply alter the 4 or 5 sines waves to show that you can provide just as good a fit without the AHH curve?
If you can, please present it here.
And if you cannot, then VP remains uncontested.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
HUH?
Let me get this correct.
I publish a paper showing that the rise and fall of women’s skirts plus a saw tooth pattern provides a good fit to the curve. Since no one can provide a better ‘fit’ than that the paper has to stand?
vukcevic says:
December 13, 2012 at 12:22 pm
So it is OK for Svalgaard of Stanford to derive IDV from changes of two combined magnetic fields, but is not for Vukcevic.
It is OK to derive two time series of the external field driven by the same source, but not to confuse and mix the external and internal fields that have different sources and don’t interact. Case is indeed closed, as you are incapable of learning.
As mentioned by Steveta_uk and others ….
Rather than engage in histrionics, the way to refute the Vaughan Pratt poster is to create a similar spreadsheet (or modify his spreadsheet) to show a nominal AGW signal.
As nearly as I can tell, Dr. Pratt has done everything responsible skeptics ask:
– Formulated a hypothesis
– Presented all of the supporting data
– Published in an accessible forum
– Asked for feedback.
I have not dropped in on the thread for a couple of days. However, I suspect that if someone has published a spreadsheet model that refutes Dr. Pratt’s, then it would have been mentioned here.
I do not care whether four parameters can fit an elephant. I would like to see someone mathematically refute Dr. Pratt’s model. I took a look and realized I do not have the time to reacquire the expertise to do it. (I had the expertise years ago and even have the optimization code I wrote for my AI class that could be adapted to this problem).
In my case, I strongly believe there are contradictory cases (it is just math and there are a lot of variables), but until someone devotes the mental sweat to create one (maybe Nick Scafetta has per Steven Mosher @ur momisugly 9:04 am), Dr. Pratt’s result stands as he has described it. He asks people to show the contradictions.
Finally re circularity. I agree that that post hoc curve fitting can be described as circular. All of the GCMs do it to reproduce historical temperature. What Dr. Pratt has done is simplify the curve fitting to a spreadsheet we can all use.
I am with Steveta, on this one. Unless someone comes up with better numerology, Professor Pratt’s numerology stands.
Steveta_uk says: “What VP has said repeatedly on JC’s site is basically that if you can provide a better fit, please do.”
Why would anyone want spend time searching for a “better fit” of an exaggerated exponential, bend down my a broken filter, plus a non physically attributable wiggle to ANYTHING?
Please explain the motivation and rewards of such an exercise.
RobertInAz: I have not dropped in on the thread for a couple of days. …. Dr. Pratt’s result stands as he has described it. He asks people to show the contradictions.
Then you ought to do so before commenting , no?
He asks for criticisms but it’s fake openness. He clearly has no intent of admitting even the most blatant errors in his pseudo-paper-poster.
Oops is not in the vocabulary of this great scientific authority.
We have a few people with logic problems here.
Steven Mosher:
1. “its not circular“. Of course it’s circular. Circular Logic is when you claim to have proved something that you assumed in the first place.
2. “its not a proof or support for models“. It certainly looks like an attempt to support the models. But VP’s motives are private, which is why I used “(presumably)“.
3. “you cant retract a poster“. As Richard Courtney has pointed out, of course you can retract a poster.
4 “This is basically the same approach that many here praise when scafetta does it“.. What people think about Scafetta or anyone else is irrelevant to the content of VP’s poster. However, if it illogically helps you to accept what I say now, here is what I once commented on a Loehle and Scafetta paper: “I’m calling BS on this paper.“. http://tinyurl.com/cxxo4lw
Steveta_uk:
You say “if you cannot [provide a better fit], then VP remains uncontested“. Utter tosh. I have just contested VP’s poster. Therefore VP’s poster has indeed been contested. You do not have to put up alternatives in order to contest something.
Matthew R Marler”
You say “That is over-wrought. Vaughan Pratt described exactly what he did and found, and published the data that he used and his result.“. Yes it’s true that VP published his data, methods and result. I congratulated him on publishing all the data and workings, and I truly meant it. If only all climate scientists did that then surely climate science couldn’t have got into its current mess. But the result was still obtained by circular logic.
jack hudson says:
December 13, 2012 at 11:07 am
……………….
On statistics –
I have notice that since computers and statistical packages became readily available in the 1980’s there has been a shift away from using a trained statistician to do-it-yourself statistics. ‘Six Sigma’ in industry is an example.
The statistical training I got from the ‘Six Sigma’ program at work was absolute crap. All they taught was how to use the computer program with not even a basic explanation of different types of distribution to go with it or even the warning to PLOT THE DATA so you could see the shape of the distribution. They did not even get into attributes vs variables!
It reminds me of the shift from the use of well trained secretaries who would clean up a technonut’s English and pry the need infor out of him to having everyone write their own reports. My plant manager in desperation insisted EVERYONE in the plant take night courses in English composition.
Too bad Universities do not insist that anyone using statistics must take at least three semesters of Stat.
lsvalgaard says:
December 13, 2012 at 12:43 pm
It is OK to derive two time series of the external field driven by the same source, but not to confuse and mix the external and internal fields that have different sources and don’t interact.
As this apparatus does:
http://www.leif.org/research/Rudolf%20Wolf%20and%20the%20Sunspot%20Number.ppt#8
records combined solar and Earth’s fields
or as Vukcevic does in here:
http://www.vukcevic.talktalk.net/EarthNV.htm
calculates combined solar and Earth’s fields
Do you suggest than the combined field curve that happens to match temperature change in the N. Hemisphere is coincidental, it just appeared by chance ?
fhhaynie says:
December 13, 2012 at 12:17 pm
Substitute a 1300 year wave length sine with a max at MWP and min at LIA for the AGW function and you will get similar results.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
That is pretty much what the Chinese, Liu Y, Cai Q F, Song H M, et al. did. They used 1324 years.
@Gail. anyone using statistics must take at least three semesters of Stat.
I agree. Three semesters of statistics would be more generally useful throughout adult life than three semesters of Calc. I’m 34 years out of my B.Sc, 31 from my Ph.D. The college texts I return to most often are my Stat books, Johnson & Leone 1977.
Not that calc isn’t useful. Not that it isn’t required for “Diff_E_Q”. But statistics is the one of the only courses that by design gives you training in uncertainty, to quantify what you don’t know.
vukcevic says:
December 13, 2012 at 1:22 pm
records combined solar and Earth’s fields
It records the external and internal fields superposed by Nature and thus existing in Nature
or as Vukcevic does in here:
calculates combined solar and Earth’s fields
no, you calculate a field that does not exist in Nature by combining two that are not physically related
Do you suggest than the combined field curve that happens to match temperature change in the N. Hemisphere is coincidental, it just appeared by chance ?
I say that the quantity you calculate does not exist in Nature and therefore that any correlation is spurious or worse. But I thought your case was closed. Keep it that way, instead of carpet bombing every thread on every blog with it.
Can’t you just take the fourier transform of the raw data to find the frequency components and phase shifts and whatever else is left over?
I do applaud the guy’s work with sines and exponentials. At least it isn’t the standard linear regression garbage!!!
Vaughan Pratt has commented on ClimateEtc that “The game therefore is to come up with two curves we can call NAT and HUM for nature and humans, such that HadCRUT3 = NAT + HUM + shorter-term variations, where NAT oscillates without trending too far in a few centuries, and some physically justifiable relation between HUM and the CDIAC CO2 data can be demonstrated.
To within MRES I’ve proposed SAW for NAT, AGW for HUM, […]“.
In order to show that VP’s version of NAT and HUM are unsupported, it is sufficient to show, as I have done, that they are unsupported. It is not necessary to propose a different version.
However, others have looked at NAT. eg, http://wattsupwiththat.com/2010/09/30/amopdo-temperature-variation-one-graph-says-it-all/
I haven’t investigated their workings, so I am not in a position to say whether their graph is worth anything, but at least it is using real data on the PDO and AMO. If they have got it right (NB. that’s an “If”), then they have nailed NAT, and HUM looks to be around a flat zero.
lsvalgaard says:
December 13, 2012 at 1:41 pm
………
I calculate combined effect of two variables, as you could calculate combined effect of wind and temperature on the evaporation, but in my case it happens that both variables are magnetic fields.
Are you happy now?
I post on other blogs so the readers also should be aware You don’t need to follow me around, if you think it is not worth your attention. Why are you so concerned ?
In a way I am pleasantly surprised that you are devoting all your attention to my ‘nonsense’ rather than the ‘brilliant’ work of your Stanford colleague, discussed above. Either you think Dr. Pratt’s work is of a superb quality or utter rubbish, in either case no comment of yours is required.
Good night.
vukcevic says:
December 13, 2012 at 3:01 pm
I calculate combined effect of two variables, as you could calculate combined effect of wind and temperature on the evaporation, but in my case it happens that both variables are magnetic fields.
Are you happy now?
Wind and temperature and evaporation are physically related. Your inputs are not. that they are both magnetic fields is irrelevant, it makes as much sense to combine them as it would the fields of the Sun and Sirius.
I post on other blogs so the readers also should be aware
I think the readers are ill served with nonsense.
You don’t need to follow me around, if you think it is not worth your attention. Why are you so concerned ?
Because scientists have an obligation to combat pseudo-science and provide the public with correct scientific information. Even though not all do that.
In a way I am pleasantly surprised that you are devoting all your attention to my ‘nonsense’ rather than the ‘brilliant’ work of your Stanford colleague, discussed above.
You should be ashamed of peddling your nonsense, not pleased when found out.
Either you think Dr. Pratt’s work is of a superb quality or utter rubbish
Curve fitting is what it is. If one believes in it has little bearing on the mathematical validity of the fitting procedure. I asked him to make an experiment for me and the result was that what he called the ‘solar curve’ was different in solar data and in CET and HadCRUT3 temperature data and between the latter two as well. This settled the matter for me at least.
Mike Jonas: But the result was still obtained by circular logic.
In filtering, there is a symmetry: if you know the signal, you can find a filter that will reveal it clearly; if you know the noise, you can design a filter to reveal the signal clearly. Pratt assumed a functional form for the signal (he said so at ClimateEtc), and worked until he had a filter that revealed it clearly.
The thought process becomes “circular” if you “complete the circle”, so to speak, and conclude that: since he found what he assumed, then it must be true. My only claim is that, given what he did, the result can be, and should be, tested on future data. I have written about the same regarding the modeling of Vukcevic and Scafetta. I would say the same regarding the curve-fitting of Liu et al cited by Gail Combs above. Elsewhere I have written the same of the modeling of Latif and Tsonis, and of the GCMs. I do not expect any extant model to survive the next 20 years’ worth of data collection, but I think that the data collected to date do not clearly rule out very much — though alarmist predictions made in 1988-1990 look less credible year by year.
Mike Jonas: However, others have looked at NAT. eg, http://wattsupwiththat.com/2010/09/30/amopdo-temperature-variation-one-graph-says-it-all/
I haven’t investigated their workings, so I am not in a position to say whether their graph is worth anything, but at least it is using real data on the PDO and AMO. If they have got it right (NB. that’s an “If”), then they have nailed NAT, and HUM looks to be around a flat zero.
Well said.
RobertInAz – “As nearly as I can tell, Dr. Pratt has done everything responsible skeptics ask:
– Formulated a hypothesis” – Yes, he did that.
“- Presented all of the supporting data” – Yes, he did that, and I congratulated him for it.
“- Published in an accessible forum” – Yes, he did that.
“- Asked for feedback.” – Yes, he did that.
And just in case you haven’t noticed, I have given him feedback. That feedback stated in no uncertain terms, and with full explanation, using the same data as Dr. Pratt, that he had used circular logic and that his findings were worthless.
Matthew R Marler – “The thought process becomes “circular” if you “complete the circle”, so to speak, and conclude that: since he found what he assumed, then it must be true. My only claim is that, given what he did, the result can be, and should be, tested on future data.“.
VP’s claimed results flowed from his initial assumptions. That’s what makes it circular.
And just in case you didn’t notice, I tested the key part of his result (the sawtooth) against existing data (the PDO and AMO) and found that it did not represent the “multidecadal ocean oscillations” as claimed.
So (a) the logic was circular, and (b) it was tested anyway and found wanting.
I submit that a simpler and better fit of the unfiltered data is 0.573-0.973sine(x/608+.96)+0.108sine(x/63+1.21)+0.038sine(x/20+1.46) where x=2PI*year. AGW may be covarient with that 608 year cycle and contributes a little bit to the magnitude of the coefficient -.973. Most of the residual looks like a three to five year cycle.