Guest Post by Willis Eschenbach
[see Update at the end of this post]
I got to thinking about the (non) adjustment of the GISS temperature data for the Urban Heat Island effect, and it reminded me that I had once looked briefly at Anchorage, Alaska in that regard. So I thought I’d take a fresh look. I used the GISS (NASA) temperature data available here.
Given my experience with the Darwin, Australia records, I looked at the “homogenization adjustment”. According to GISS:
The goal of the homogenization effort is to avoid any impact (warming or cooling) of the changing environment that some stations experienced by changing the long term trend of any non-rural station to match the long term trend of their rural neighbors, while retaining the short term monthly and annual variations.
Here’s how the Anchorage data has been homogenized. Figure 1 shows the difference between the Anchorage data before and after homogenization:
Figure 1. Homogenization adjustments made by GISS to the Anchorage, Alaska urban temperature record (red stepped line, left scale) and Anchorage population (orange curve, right scale)
Now, I suppose that this is vaguely reasonable. At least it is in the right direction, reducing the apparent warming. I say “vaguely reasonable” because this adjustment is supposed to take care of “UHI”, the Urban Heat Island effect. As most everyone has experienced driving into any city, the city is usually warmer than the surrounding countryside. UHI is the result of increasing population, with the accompanying changes around the temperature station. More buildings, more roads, more cars, more parking lots, all of these raise the temperature, forming a heat “island” around the city. The larger the population of the city, the greater the UHI.
But here’s the problem. As Fig. 1 shows, until World War II, Anchorage was a very sleepy village of a few thousand. Since then the population has skyrocketed. But the homogeneity adjustment does not match this in any sense. The homogeneity adjustment is a straight line (albeit one with steps …why steps? … but I digress). The adjustment starts way back in 1926 … why would the 1926 Anchorage temperature need any adjustment at all? And how does this adjust for UHI?
Intrigued by this oddity, I looked at the nearest rural station, which is Matanuska. It is only about 35 miles (60 km) from Anchorage, as shown in Figure 2.
Figure 2. Anchorage (urban) and Matanuska (rural) temperature stations.
Matanuska is clearly in the same climatological zone as Anchorage. This is verified by the correlation between the two records, which is about 0.9. So it would be one of the nearby rural stations used to homogenize Anchorage.
Now, according to GISS the homogeneity adjustments are designed to adjust the urban stations like Anchorage so that they more closely match the rural stations like Matanuska. Imagine my surprise when I calculated the homogeneity adjustment to Matanuska, shown in Figure 3.
Figure 3. Homogenization adjustments made by GISS to the Matanuska, Alaska rural temperature record.
Say what? What could possibly justify that kind of adjustment, seven tenths of a degree? The early part of the record is adjusted to show less warming. Then from 1973 to 1989, Matanuska is adjusted to warm at a feverish rate of 4.4 degrees per century … but Matanuska is a RURAL station. Since GISS says that the homogenization effort is designed to change the “long term trend of any non-rural station to match the long term trend of their rural neighbors”, why is Matanuska being adjusted at all?
Not sure what I can say about that, except that I don’t understand it in the slightest. My guess is that what has happened is that a faulty computer program has been applied to fudge the record of every temperature station on the planet. The results have then been used without the slightest attempt at quality control.
Yes, I know it’s a big job to look at thousands of stations to see what the computer program has done to each and every one of them … but if you are not willing to make sure that your hotrod whizbang computer program actually works for each and every station, you should not be in charge of homogenizing milk, much less temperatures.
The justification that is always given for these adjustments is that they must be right because the global average of the GISS adjusted dataset (roughly) matches the GHCN adjusted dataset, which (roughly) matches the CRU adjusted dataset.
Sorry, I don’t find that convincing in the slightest. All three have been shown to have errors. All that shows is that their errors roughly match, which is meaningless. We need to throw all of these “adjusted datasets” in the trash can and start over.
As the Romans used to say “falsus in unum, falsus in omnibus”, which means “false in one thing, false in everything”. Do we know that everything is false? Absolutely not … but given egregious oddities like this one, we have absolutely no reason to believe that they are true either.
Since people are asking us to bet billions on this dataset, we need more than a “well, it’s kinda like the other datasets that contain known errors” to justify their calculations. NASA is not doing the job we are paying them to do. Why should citizen scientists like myself have to dig out these oddities? The adjustments for each station should be published and graphed. Every single change in the data should be explained and justified. The computer code should be published and verified.
Until they get off their dead … … armchairs and do the work they are paid to do, we can place no credence in their claims of temperature changes. They may be right … but given their egregious errors, we have no reason to believe that, and certainly no reason to spend billions of dollars based on their claims.
[Update – Alaska Climate Research Center releases new figures]
I have mentioned the effect of the Pacific Decadal Oscillation (PDO) below. The Alaska Climate Research Center have just released their update to the Alaska data. Here’s that information:
Figure 4. Alaska Temperature Average from First Order Observing Stations
In the Alaska Climate Research Center data, you can clearly see the 1976 shift of the PDO from the cool to the warm phase, and the recent return to the cool phase. Unsurprisingly, the rise in the Alaska temperatures (typically shown with a continuously rising straight trend line through all the data) have been cited over and over as “proof” that the Arctic is warming. However, the reality is a fairly constant temperature from 1949-1975, a huge step change 1975-1976, and a fairly constant temperature from 1976 until the recent drop. Here’s how the IPCC Fourth Assessment Report interprets these numbers …
Figure 5. How the IPCC spins the data.
SOURCE: (IPCC FAR WG1 Chapter 9, p. 695)
As you can see, they have played fast and loose with the facts. They have averaged the information into decade long blocks 1955-1965, 1965-1975, 1975-1985 etc. This totally obsures the 1975-1976 jump. It also gives a false impression of the post-1980 situation, falsely showing purported continuing warming post 1980. Finally, they have used “adjusted data” (an oxymoron if there ever was one). As you can see from Fig. 4 above, this is merely global warming propaganda. People have asked why I say the Alaska data is “fudged” … that’s a good example of why.





Willis,
Can you confirm something for me?
I see the GISS homogenization adjustments for Anchorage in fig #1 as follows:
1) In the mid 1920s GISS added .9C to the raw temperature data to get their homogenized adjusted temp data
2) around year 2000 GISS added 0.0C to the raw temperature data to get their homogenized adjusted temp data
3) in the ~85 years between GISS linearly (in stepwise fashion) decreased what the degrees C that they added to the raw temperature data to get their homogenized adjusted temp data
4) From late 1940s to present there was a significantly amount of urbanization as evidenced by the population curve.
Do I have the items 1), 2) and 3) above right?
Is so then I have some further observations to make, but first making sure I got the above right.
John
Willis Eschenbach (00:35:49) :
‘Robert may or may not provide an answer, but it gives lurkers and posters alike the chance to ponder the issues and answer the questions for themselves. And that to me is a worthwhile thing.’
Willis. The patience of Job and the wisdom of Solomon make you a rear beast indeed. Chuck in your dry any pithy humour and you might have the makings of a ‘universal man’!
But compliments aside, your comments have been enlightening for me and for that many thanks
Regards
Doug
I think I’ve figured out what is going on here. The relevant parts of the code are in directory STEP2 of the GISTEMP source. The first thing to note is the file v2.inv in ./input_files. It lists the station data, and for the two Alaska stations:
42570273000 ANCHORAGE/INT__________________ 61.17 -150.02__ 40____8U__173FLxxCO 1A 5WATER__________ C__ 53
42570274001 MANTANUSKA AES__________________61.57 -149.27__ 46__225R__ -9FLxxCO30x-9TUNDRA__________C__ 18
The second last number gives the GHCN brightness rating – C is highest. So both will be adjusted.
The adjustment is done in the subroutine adj() of padjust.f. Temperatures are held as integers, to tenths of a degree – 121 is 12.1C. The relevant code fragment is:
____ do iy=iy1,iy2
________sl=sl1
________if(iy.gt.knee) sl=sl2
________iya=iy
________if(iy.lt.iy1a) iya=iy1a
________if(iy.gt.iy2a) iya=iy2a
________iadj=nint( (iya-knee)*sl-(iy2a-knee)*sl2 )
iy is the year. You’ll see that there is provision for a “knee” and two slopes (switching at the knee). There are also limits beyond which the adjustment will be held stationary.
So the effect of the “nint” is that the adjustment is piecewise linear, and forced to the nearest integer value – ie 0.1C.
So now you can see where those plots come from. Anchorage has no knee, but just a steady slope adjustment, made stepwise by the nint. Mantanuska has a knee, and a steep slope following the knee.
The rationale seems to be that a broad slope correction is made to match the city to have the trend of its rural surrounds. The knee is presumably to allow for a transition in the past from rural to urban.
Sould be rare beast! sorry Willis!
The common factor in the errors is the expectation of those manipulating the Data: that means if the expectations are wrong the conclusions are wrong.
True scientists alter the theory to fit the facts not the facts to fit the theory.
‘ Steinar Midtskogen (22:36:00) :
“Falsus in unum, falsus in omnibus”
Actually, the expression should rather be “falsus in uno, falsus in omnibus”. If you can’t get the Latin right, can you then be right about the rest? ;)’
Yes, you can. As in “coito, ergo sum”.
Raw data is good for you and may be digested later as food for thought.
I was in the military and spent a large proportion of my sdervice in norfolk(uk) home of the uea. We used to call the locals carrot crunchers. So my question is, carrot eater, is that you philj or briffa
?
John Whitman (01:39:04)
John, I don’t understand how or why they constructed the stepwise change in temperature. Your explanation seems as good as any. However, why on earth would they want to change the mid 1920s temperature by adding almost a full degree to it?
In any case, make your observations, always welcome …
Re: Nick Stokes (Feb 22 02:00),
More on the algorithm used. In PApars.f they compute, as they stated, for each urban station (in the do 200 loop) the AVG() of yearly distance-weighted average values of nearby rural stations. They then calculate the difference between that and the urban values URB(). This information is passed to getfit() in t2fit.f (via the common block FITCOM).
In getfit(), the lines-with-knee fit to AVG-URB is laboriously computed by trying every possible knee (within 5 yrs of the end of range), and choosing the value with least residual SS (in trend2(), tr2.f). The resulting slopes and knee are what ends up in the adj() routine set out in my previous post, and with the steppedness of integer conversion become the plots shown in the head post here.
Willis Eschenbach (23:19:45) :
That’s the whole point of the GISS adjustment. The stations labeled as urban do not get to impact the long term trends in the final result. They’re left in there to contribute some short-term wiggles to the local averages.
So yes, they’re essentially tossing out the urban stations. If you don’t like that, then follow the GHCN adjustment set instead. And if you don’t like adjustments at all, you might note that the unadjusted data give about the same global mean trend, anyway.
As for your continual wonderment at what happened: just go and see the neighboring dark rural stations. Then apply the two-legged adjustment described in the papers, if you really want to recreate the math. Do you want to simply wonder, or do you want to learn the answers? There is nothing stopping you from the latter.
Willis Eschenbach (23:19:45) :
Incidentally, you’re using the word ‘homogenize’ incorrectly in this context. In this field, inhomogeneities refer to the things that mess with the record at any given station: station moves, instrument changes, etc.
Nick Stokes (03:15:17) :
The ‘knee’ used to be at a fixed point in time, I think 1950. Allowing it to float makes for a better calculation, but adds to the computation as you can see.
Rural it is then. Don’t quite understand how it scores 18. Moonlight reflected in they eyes of Wolves perhaps?
Wouldn’t it be great if the warmists and the sceptics could agree a number of rural and always-been-rural sites, and look at the raw data without any fiddle factors, and see just how much warming there really is…
Why elevate the erly temperatures for an urban site?
Well, the net effect is to reduce the apparent UHI effect.
That means that later uban data can be adopted in the later years with a lesser correction meaning that the average leter years temperatures is elevated.
This isn’t about reducing he apparent warming on this site but getting the later temperatures into the data set with no signinificant reduction – i.e. don’t pull down the later temperatures for the expanded city, bring up the earlier temeratures.
This works if the earlierr temperatures do not significantly increase the average global temp but the later urban dat, if included with minimal correction, does significantly increase the mean global temp.
AT least, that’s one possibility. unless I have misinterpretted this somwehere.
There are 29 rural (by light at night) stations in Alaska within 1000km of Anchorage. Here they are listed by increasing distance in km.
DistanceStation
109SKWENTNA
126TALKEETNA
177PUNTILLA
245SUMMIT/WSO AIRPORT
252FAREWELL FAA AP
256CORDOVA/MILE
265GULKANA/INTL.
281MIDDLETON ISLAND/AUTO
289MCKINLEY PARK
311ILIAMNA FAA AP
323MINCHUMINA
334CAPE SAINT ELIAS ALASKA, U
427YAKATAGA/AIRPORT
457TANANA
466NORTHWAY FAA AP
508ANIAK/AIRPORT
526HOLY CROSS
598YAKUTAT
600EAGLE
614ALLAKAKET
629UNALAKLEET
643BETTLES
648FT YUKON UNITED
722MOSES POINT
855GUSTAVUS/2 SW
880KOTZEBUE, RAL
917UMIAT
944ANNEX CREEK
994COLD BAY
It appears that the Weather Underground mean temperature for January in Anchorage is a full 1 degree plus C cooler than the GISS reading (16 F or -8.9 C, vs. -7.4 C). Given that 2009 readings at both places were similar, this seems a little strange.
I’ve also noticed this in the Vladivostok readings, which registered a full 3 degree C difference (-12.3 C vs. 4 F or -15.6 C)
It appears that this difference has existed for some time (different locations?), but that the gap between readings has been increasing over the past decade:
WU GISS
1997: 5 F (-15C); -13.7 C
1998: 4 (-15.6); -14.3
1999: 11 (-11.7); -10.3
2000: 3 (-16.1); -15.1
2001: 1 (-17.22); -16.2
2002: 12 (-11.1); -9.4
2003: 7 (-13.9); -12.1
2004: 7 (-13.9); -11.2
2005: 7 (-13.9); -11.8
2006: 6 (-14.4); -12.6
2007: 15 (-9.4); -7.1
2008: 1 (-17.2); –12.8
2009: 10 (12.2); – 10.8
2010: 4 (-15.6); -12.3
What gives?
Nick Stokes (05:31:29) :
I think they only use those within 500 km if there are enough stations within that range, which appears to be the case.
So now that we have a list of the stations we need to look at, to answer all of the questions Eschenbach raises…
Thank you Willis for providing that extra GISS link on index brightness.
It was interesting to look at the entries for my country, Australia, and other countries which interest me such as Estonia, Latvia and Lithuania.
I have some nitpiky comments.
The information on which the Brightness Index is calculated could be out of date. Hobart Australia has an index of 34 and Canberra an index of 30. Canberra’s population has overtaken Hobart’s by at least 100,000. They both however end up as C in the metadata so it doesn’t matter.
There is a similar disparity in the index for Vilnius and Kaunas (10 and 17) in Lithuania.
To find the cities for Estonia, knowledge of Baltic, Hanseatic and some modern history helps.
Tallinn: The capital and the second n is dropped in both sets of data.
Baltischport :This is the Russian name for the town of Paldisksi about 50 km west of Tallinn. It is marked as USSR.
Pjarnu: This seems to be the Russian spelling for Parnu.
Dorpat: This is the German name for Tartu, Estonia’s second largest city. It is also stated as being in the USSR.
Fellin: This is the German name for Viljandi. Again another USSR listing.
The Latvian data has the same oddities of German names and USSR locations.
Sure the co-ordinates (which I didn’t check for accuracy) are the critical elements for location and these are just nitpiks. But it is a tad sloppy.
pat (18:44:46)
AGW has been declared a religion by a UK court in 2009. An employee took his company to court for preventing him or sacking him because of his religious belief i AGW. Fact
Veronica (England) (05:20:11) :
This was done by a north american school boy last year and somewhere on the climate realist sites there is a video of his explanations and chart. Try a google search or perhaps one of the readers here can remember where it was.
Can anybody (or hopefully somebody) provide the world, or at least me, with the confidence that:
1) all these immense efforts and all this time and money spent by good as well as bad scientists on collecting and interpreting (or deliberately misinterpreting) hundreds of thousands of temperature data can ever lead to a believed average world temperature graph?
2) that there is any use in attempting to do this, BEFORE it has been established with certainty that (a possibly dangerous) upward or downward trend is caused by human activity rather than natural forces?
3) if neither, that the efforts mentioned under 1) are still worthwhile?
Thinking outside the box …
Steve Mosher advises that: “UHI is, according the theory, a positive bias that is created by changes to the physical (geometry) and material changes in a site over time, and waste heat which is tied to population and human activity.”
Since human activity produces UHI warming, and in light of the enormous growth in global population centers resulting from industrialization, perhaps the IPCC religionists have been right all along about AGW … but for the wrong reason.
Wanted to ask this again……along the lines of what the Carrot Eater mentioned above…
Question – can we compare global (or N Hemi) averages of raw data vs.
“adjusted data” to see the differences imposed by the cumulative adjustments? just curious.
Also along the same lines, what ever happened to the comparisons between the ‘good’ US stations vs. the combined ‘good plus bad’ US stations as determined by SurfaceStations project? Did this show a significant difference so that we can get a handle on the overall impact of the siting issues (?)
Thanks