Spencer: Direct Evidence that Most U.S. Warming Since 1973 Could Be Spurious

by Roy W. Spencer, Ph. D.

Where the lights are - the CONUS has population that roughly tracks with brightness

INTRODUCTION

My last few posts have described a new method for quantifying the average Urban Heat Island (UHI) warming effect as a function of population density, using thousands of pairs of temperature measuring stations within 150 km of each other. The results supported previous work which had shown that UHI warming increases logarithmically with population, with the greatest rate of warming occurring at the lowest population densities as population density increases.

But how does this help us determine whether global warming trends have been spuriously inflated by such effects remaining in the leading surface temperature datasets, like those produced by Phil Jones (CRU) and Jim Hansen (NASA/GISS)?

While my quantifying the UHI effect is an interesting exercise, the existence of such an effect spatially (with distance between stations) does not necessarily prove that there has been a spurious warming in the thermometer measurements at those stations over time. The reason why it doesn’t is that, to the extent that the population density of each thermometer site does not change over time, then various levels of UHI contamination at different thermometer sites would probably have little influence on long-term temperature trends. Urbanized locations would indeed be warmer on average, but “global warming” would affect them in about the same way as the more rural locations.

This hypothetical situation seems unlikely, though, since population does indeed increase over time. If we had sufficient truly-rural stations to rely on, we could just throw all the other UHI-contaminated data away. Unfortunately, there are very few long-term records from thermometers that have not experienced some sort of change in their exposure…usually the addition of manmade structures and surfaces that lead to spurious warming.

Thus, we are forced to use data from sites with at least some level of UHI contamination. So the question becomes, how does one adjust for such effects?

As the provider of the officially-blessed GHCN temperature dataset that both Hansen and Jones depend upon, NOAA has chosen a rather painstaking approach where the long-term temperature records from individual thermometer sites have undergone homogeneity “corrections” to their data, mainly based upon (presumably spurious) abrupt temperature changes over time. The coming and going of some stations over the years further complicates the construction of temperature records back 100 years or more.

All of these problems (among others) have led to a hodgepodge of complex adjustments.

A SIMPLER TECHNIQUE TO LOOK FOR SPURIOUS WARMING

I like simplicity of analysis — whenever possible, anyway. Complexity in data analysis should only be added when it is required to elucidate something that is not obvious from a simpler analysis. And it turns out that a simple analysis of publicly available raw (not adjusted) temperature data from NOAA/NESDIS NOAA/NCDC, combined with high-resolution population density data for those temperature monitoring sites, shows clear evidence of UHI warming contaminating the GHCN data for the United States.

I will restrict the analysis to 1973 and later since (1) this is the primary period of warming allegedly due to anthropogenic greenhouse gas emissions; (2) the period having the largest number of monitoring sites has been since 1973; and (3) a relatively short 37-year record maximizes the number of continuously operating stations, avoiding the need to handle transitions as older stations stop operating and newer ones are added.

Similar to my previous posts, for each U.S. station I average together four temperature measurements per day (00, 06, 12, and 18 UTC) to get a daily average temperature (GHCN uses daily max/min data). There must be at least 20 days of such data for a monthly average to be computed. I then include only those stations having at least 90% complete monthly data from 1973 through 2009. Annual cycles in temperature and anomalies are computed from each station separately.

I then compute multi-station average anomalies in 5×5 deg. latitude/longitude boxes, and then compare the temperature trends for the represented regions to those in the CRUTem3 (Phil Jones’) dataset for the same regions. But to determine whether the CRUTem3 dataset has any spurious trends, I further divide my averages into 4 population density classes: 0 to 25; 25 to 100; 100 to 400; and greater than 400 persons per sq. km. The population density data is at a nominal 1 km resolution, available for 1990 and 2000…I use the 2000 data.

All of these restrictions then result in thirteen 24 to 26 5-deg grid boxes over the U.S. having all population classes represented over the 37-year period of record. In comparison, the entire U.S. covers about 31 40 grid boxes in the CRUTem3 dataset. While the following results are therefore for a regional subset (at least 60%) of the U.S., we will see that the CRUTem3 temperature variations for the entire U.S. do not change substantially when all 31 40 grids are included in the CRUTem3 averaging.

EVIDENCE OF A LARGE SPURIOUS WARMING TREND IN THE U.S. GHCN DATA

The following chart shows yearly area-averaged temperature anomalies from 1973 through 2009 for the 13 24 to 26 5-deg. grid squares over the U.S. having all four population classes represented (as well as a CRUTem3 average temperature measurement). All anomalies have been recomputed relative to the 30-year period, 1973-2002.

The heavy red line is from the CRUTem3 dataset, and so might be considered one of the “official” estimates. The heavy blue curve is the lowest population class. (The other 3 population classes clutter the figure too much to show, but we will soon see those results in a more useful form.)

Significantly, the warming trend in the lowest population class is only 47% of the CRUTem3 trend, a factor of two difference.

Also interesting is that in the CRUTem3 data, 1998 and 2006 would be the two warmest years during this period of record. But in the lowest population class data, the two warmest years are 1987 and 1990. When the CRUTem3 data for the whole U.S. are analyzed (the lighter red line) the two warmest years are swapped, 2006 is 1st and then 1998 2nd.

From looking at the warmest years in the CRUTem3 data, one gets the impression that each new high-temperature year supersedes the previous one in intensity. But the low-population stations show just the opposite: the intensity of the warmest years is actually decreasing over time.

To get a better idea of how the calculated warming trend depends upon population density for all 4 classes, the following graph shows – just like the spatial UHI effect on temperatures I have previously reported on – that the warming trend goes down nonlinearly as population density of the stations decrease. In fact, extrapolation of these results to zero population density might produce little warming at all!

This is a very significant result. It suggests the possibility that there has been essentially no warming in the U.S. since the 1970s.

Also, note that the highest population class actually exhibits slightly more warming than that seen in the CRUTem3 dataset. This provides additional confidence that the effects demonstrated here are real.

Finally, the next graph shows the difference between the lowest population density class results seen in the first graph above. This provides a better idea of which years contribute to the large difference in warming trends.

Taken together, I believe these results provide powerful and direct evidence that the GHCN data still has a substantial spurious warming component, at least for the period (since 1973) and region (U.S.) addressed here.

There is a clear need for new, independent analyses of the global temperature data…the raw data, that is. As I have mentioned before, we need independent groups doing new and independent global temperature analyses — not international committees of Nobel laureates passing down opinions on tablets of stone.

But, as always, the analysis presented above is meant more for stimulating thought and discussion, and does not equal a peer-reviewed paper. Caveat emptor.

The climate data they don't want you to find — free, to your inbox.
Join readers who get 5–8 new articles daily — no algorithms, no shadow bans.
0 0 votes
Article Rating
162 Comments
Inline Feedbacks
View all comments
JJ
March 16, 2010 2:55 pm

Dr. Spencer,
How raw is the raw data that you used? Assume no TOB … what about instrumentation changes, etc?
What do your satellite data for CONUS show over the ’73-09 period?
JJ

harrywr2
March 16, 2010 2:56 pm

Sean Peake (13:01:03) :
“To me that seems kind of low and may reflect all of the US, not the continental US—Alaska is the likely reason because it’s half the size and has just over 1 person/sq m.”
Pop density in the US not including Alaska is about 36/km2.
OT but related.
In 1960 there where a total of 108 billion airline passenger-kilometers flown. In 2000 three were 3 trillion passenger-kilometers.
We’ve had a 30 fold increase in air travel…the vast majority of the worlds thermometers are at the airport. But the world population has only doubled.
Business’s like to locate near transportation hubs, people who need jobs locate near business’s.

NickB.
March 16, 2010 3:06 pm

R. Gates (12:26:11) :
I am wondering though, if the UHI is so widespread, at what point is it no longer just an effect that causes skewing of the data, but part of the environment? What I’m getting at is, all the extra heat that urban areas generate (and certainly this does skew the data as urbanization continues), but that extra heat in those urban area is still extra heat added to the total heat balance of the planet. Taken to an extreme, if somehow the entire planet were a giant concrete parking lot, you’d have a giant “heat island”, or would you simply have a different climate? These urban heat islands are, on some level, part of the environment, growing rapidly, and averaged over the whole globe. will add to warming.
Wow, I’m having a civil (agreeable?) conversation with you and Wren on a thread? Haha, will wonders never cease? 🙂
I stand firm that I believe these are not measurement “errors” per se, but from a surface temperature record standpoint should be, more or less, red dots instead of giant orange blobs. They should also be classified and represented as land use changes. They are, as much as I hate the term, a form of AGW but it’s not CO2-driven AGW and we call it AGW… everyone’s still going to say “CO2!!!! CO2!!!!!” So what do we call it?
We’re getting back to the “It’s more complicated than the RC folks think” – an idea that has always run through my head when I look at the “consensus” storyline of the climate.
What really gets me is that the mainstream surface temperature reconstructions (GISS in particular) allegedly have already “over-corrected” (according to Hansen or Schmidt, not sure which one) for UHI and still trend higher than Dr. Spencer’s raw reconstruction. They’re already off by a significant margin before you even start to consider UHI – which makes it almost farcical.
Then you also have to consider that per the IPCC, “land use” changes are a negative forcing (~.25 W/m2). If they are really a net positive forcing due to UHI then it means something on the positive forcing side has been exaggerated above and beyond the implications of the lower average global trend that – according to Pat Michaels – would probably be .13 instead of .17 (so ~30% overstated?).

Alba
March 16, 2010 3:06 pm

Paul Z. (08:03:52) :
“This is a very significant result. It suggests the possibility that there has been essentially no warming in the U.S. since the 1970s.”
Tell this to the pope, will ya:
http://news.bbc.co.uk/2/hi/uk_news/scotland/glasgow_and_west/8570294.stm
Scottish Secretary Jim Murphy, who is leading UK government preparations for the tour, said it was “an historic visit at an important time”.
“The Papal visit represents an unprecedented opportunity to strengthen ties between the UK and the Holy See on action to tackle poverty and climate change as well as the important role of faith in creating strong and cohesive communities,” he said.
===
Is the Vatican invested heavily in carbon emissions trading?
Well, there’s one born every minute, isn’t there. Here’s a fellow who actually thinks he can believe what politicians say!!! (LOL) Even worse, he thinks that he can believe whatever Jim Murphy says. (Even more LOL!!)

NickB.
March 16, 2010 3:27 pm

Richard M (14:34:43)
That is a very good point. While I thought I had read that “urban” areas accounted for 1.5% of the land surface (maybe the less than 1% is total earth surface), that only accounts for less than 50% of the world’s population. The other ~3.9 billion folks live out there in the “non-urban” areas – whatever that means.
I went digging and found this: http://www.demographia.com/db-intlua-area2000.htm
So that means for every country in the world “urban” areas are at least 950 people/km2 – which is already almost off Dr. Spencer’s chart. When I went looking I was thinking “oh this is going to just be noise”… but when you look at Dr. Spencer’s chart and then start to think about, for example, what % of the world’s land surface is 25 people/km2 or greater, and then it’s not quite so insignificant.
Interesting!

wayne
March 16, 2010 4:01 pm

If we had sufficient truly-rural stations to rely on, we could just throw all the other UHI-contaminated data away.
There you go, correct.
Way to go Spence! Well done.
In my opinion it might also be shown that during the average twelve hours of nighttime, light winds that tend to carry the UHI effect (energy) over rural areas affect the temperatures to a greater degree than in the nighttime to daytime transition. Warm or even hot surfaces transfer the heat energy during early evening hours as they cool but there is no such effect in the night to day transition, nothing to ever counteract this UHI effect. Follow the energy flow in time as it flexes every diurnal period! There is excess energy in the UHI areas but there is no such thing as a “Rural Cool Island” effect to carry ‘negative’ energy to the UHI areas in the night to day transition. That’s one funky way to look at it, but should be correct in the diurnal sense!
That should account for what we see in the slopes of minimum and maximum temperature graphs in rural stations. The nighttime slopes are always greater and the maximum show little effect across ~115 years except in rare cases and these probably have other effects occurring to account for this discrepancy.
My comment above should not take away from Dr. Spencer’s presentation in any way, with max and min separate graphs not showing this effect, you just never see it but instead it is inferred to me by his concept.
Another thing comes to mind, and this is in an area I have little expertise in, can the lights at nighttime that are making it impossible to find a truly night sky if you want to take a telescope and look at the nighttime sky cause a tiny UHI effect during the nighttime? Can this SW radiation from mere lighting reflect and absorb enough energy in the gases of the atmosphere to cause some of that increase of the slope of minimum temperature readings? Gases can hold little energy compared to solid or liquid matter, in other words air is easily heated (low specific heat) and is hard to radiate it away (low emissivity). That I don’t know.
But in high population dense areas, many people do not to have lights on at night, it’s already bright enough. But in isolated homes and businesses in rural areas, they are in near total darkness and most would naturally add outside lighting for safety reasons. That is why we no longer have dark skies. That single factor should take on an exponential curve also against population density but may prove too tiny to bother with, but maybe it doesn’t.
Roy, an additional analysis with day and night being separated may prove fruitful with proper physics principles backing it.
Also, while looking up NOAA/NCDC stations, the population for Boulder stood out (being over a million) and I don’t think that is correct. 🙂 How many other populations are not correct?
Remember, this was only my opinion, being a transfer of possible knowledge or an alternate viewpoint. If you are an AGW proponent and want it, take it, otherwise, just ignore it. These comments are just that, comments, not scientific papers, they come later.

March 16, 2010 4:01 pm

pat michaels (11:06:21) :
Roy–
This is exactly the result that McKitrick and I got globally (J. Geophys 08). We found 50% of the land warming in the NASA record was from “non-climatic” factors.
You can’t have a UHI over the ocean, so these findings apply to about 30% of the surface. After you factor in the relative warmings of each hemisphere separately, I believe the decadal trend will drop by .04deg–not an inconsiderable amount.
The drop will be from 0.17 to 0.13; isn’t it odd that the 0.13 is almost exactly what the MSU shows?

This would make things very neat and tidy – if it weren’t for the fact that MSU shows a US warming trend of 0.22 deg per decade, i.e. about the same as the GHCN trend.

Ron1
March 16, 2010 4:04 pm

I don’t know if this point has been raised, but what does a flat temperature line during the period in question say about sea surface temperatures as the oceans warm or cool the land masses.
This would seem to suggest that they too must have had near zero anomalies. While not tested, I would not be surprised that this pattern repeated in other regions of the world which if true would also point to stable SSTs during this period.

Capn Jack.
March 16, 2010 4:12 pm

I dont have the math anymore to argue the statistics, but I like the approach of testing the UHI contamination to the thermometer record (raw) and then using that logically to filter the UHI component out off the temperature record as an attempt to measure real national warming over a significant land mass. The UHI is a testable.
There is nothing wrong with the arbitrary starting point as it meets basic criteria, maximum data integrity and the time line of policy approaches to warming attributed to man made CO2.

R. Gates
March 16, 2010 4:18 pm

Nick B.,
I won’t tell anyone, if you won’t, that we’re having a somewhat civil conversation. 🙂
You really do raise some valid points. I know that this thread is supposed to be about the unreliability of the temperature data because of UHI, but it certainly does make sense that the very thing (i.e. urbanization) that is “skewing” the data, is also creating real heat that spreads into the global environment because so far at least there are not any bubbles built over the large maga-cities like NY, Tokyo, London, Singapore, Syndey, LA, etc. But the heat from these cities spreads globally, and certainly, by some fraction of degree, they certainly must raise the overall global temperature.
Interestingly enough, the energy to produce this heat comes from many sources, and not just fossil fuels, with the primary source being probably sunlight. Sunlight striking brick, concrete, asphalt, metal etc. will turn that energy into infrared heat energy, versus the sunlight that falls on grass and trees, which is mainly converted into chemical energy. In effect then, the cities become turn into little (or mega) UHI generating islands. The final twist on this is of course the fact that the megacities both increase the overall heating of the earth, but also are net producers of CO2, whereas the plants or ground cover in existence before there were cities probably used more CO2, or at least were carbon neutral. To the extent then that rising CO2 levels really are a global climate issue, the heat and CO2 produced by cities becomes a multi-threat to the global climate status quo.
Where will this all end?
1) Humans will continue to expand and use up all resources and we all die.
2) Humans will perfect fusion technology and with an infinitie amount of energy we can make the earth into a garden of eden and we then expand into the Solar System and beyond, until we contact an intelligent civilization that either greets us with open arms, enslaves us, or finds us tasty and serves us as snacks at their alien parties.
3) Humans will muddle along in their human like way until the next glacial period comes and everything we’ve ever built is covered with 2 miles of ice for the next 100,000 years.
4) A coronal mass ejection wipes out 99.9% of life on earth and the whole thing starts again.
No matter what happens…you can be sure that the AGW believers and the AGW skeptics will be at each others throat until the bitter end…

Geoff Sherrington
March 16, 2010 5:07 pm

Lon Hocker (11:01:45) :
“In short, I feel that imagining a zero value of temperature rise for the last 30 years in the US is pushing your idea quite a bit too far”.
There are many places in the Southern hemisphere that are truly rural and which need no UHI correction. Among them are many that show zero or negative temperature change over the last 40 years, but there is a lot of noise. Trend fitting is questionable in individual cases but might have significance when you combine many cases.
The logical first assumption is that there has been no global temperature change. If you then observe that there might have been a change, you seek an explanation. That’s all Dr Spencer has done.
If you believe that it is unlikely that parts of the world have avoided a temperature change, then you still have to explain the many rural sites where there has been no temperature change.
Look again at Macquarie Island, out in the middle of nowhere, influenced almost not at all by the hand of man. There is very minor infilling of missing data. Data from Bureau of Meteorology, Australia.
http://i260.photobucket.com/albums/ii14/sherro_2008/2010/mACQUARIE.jpg?t=1268784258
Why does it look immune from GHG?

John from MN
March 16, 2010 5:17 pm

Thanks so very much for your diligent work Roy.
Your work aligns perfectly of my observsions in my rural are of S.MN……Temprature over the last 120 years is just very slightly positive. So if all would be contributed from increased level of Co2 in the atmosphere, we could surmize that yes Co2 has a small positive effect on surface tempratures. But as I suspected it is so very small the slight warming is definitely a postive effect on mankind. To bad the scientists can not un-lock themselves from the politics and see the light. Yes the increase from 350 ppm to 390 ppm warm the surface ever so much. But the empirical data shows as the co2 increases it become a smaller and smaller positive and surface tempratures concur with that……..John….

March 16, 2010 5:18 pm

Roy
The thing that bothers me is that I’m not sure there is a completely satisfactory explanation for your results. Let’s consider a hypothetical case. In 1973 the population density of Smallsville is 10 per sq km; in 2009 it is 20 per sq km. Although its population has doubled it’s still in the low population class. Metropolis (can’t think where I got these names from), on the other hand, had a pop density 500 per sq km in 1973 which is the same as in 2009, i.e. it’s unchanged. If population density were a factor then we might expect the temperature trend in Smallsville to be more affected than the trend in Metropolis. What I’m saying is that there doesn’t necessarily seem to be an obvious reason why a low density population station should be any less ‘contaminated’ than a high density population station. In fact, in one of your previous posts, you implied that there was a log relationship and that a population rise in a rural location would, therefore, have a disproportionately greater effect than a similar rise in an urban location.
Now I’m not familiar with US demographics so it’s quite possible someone is going to tell me that population growth in the past 30 odd years has been restricted almost entirely to the major cities. If so – fair enough – that answers my point. But failing that – are we quite sure there isn’t a regional effect going on, i.e. is it possible that a number of the low pop stations are clustered in a region of the country which just happens to have warmed at a lower rate. Unlikely, perhaps, but I still think it needs to be ruled out.
Interesting post, anyway.

TERRY BEAR
March 16, 2010 5:52 pm

GREAT WORK DOC!
This research plus what Dr. Jones admitted in his most recent BBC interview:
NO statistically significant warming from 1995-2010 while C02 increased by 30+%. In my book that is enough to reject the NULL: Mean Global Temperature (>) greater than (i.e. AGW) than historic mean and supports the ALTERNATIVE HYPOTHESIS Mean Global Temperature (=<) (Skeptic) equal to or less than historic mean temperature. Plus Jones admitted that from 2002 to 2010 there has been slight cooling. He further admitted that the MWP was as hot if not hotter than today. This is a critical to getting some sanity back in the discussion of science in our society today.
THANKS DOC!
tjl

David44
March 16, 2010 5:53 pm

John Finn (17:18:26) “Now I’m not familiar with US demographics so it’s quite possible someone is going to tell me that population growth in the past 30 odd years has been restricted almost entirely to the major cities.”
Yes, cities and suburbs. I don’t have numbers, but with the mechanization of agriculture, the loss of small factories in small towns, and the attraction of urban amenities, there has been a strong rural to urban migration in the U.S. since at least the 1950’s. (Not counting migration of retirees to Sun Cities and such, which I think are more suburban than rural.)

rbateman
March 16, 2010 6:00 pm

Keep at it Roy, you’re doing fine.
Is there an upper limit to UHI?

E.M.Smith
Editor
March 16, 2010 6:02 pm

Lon Hocker (11:01:45) : In short, I feel that imagining a zero value of temperature rise for the last 30 years in the US is pushing your idea quite a bit too far.
I think it’s fine.
I made a variation on First Differences that I call dT/dt. The variation is that I start in the past going forward to make the anomalies, but start in the present and take my anomalies running total going backwards in time (on the assumption that “now” is a better record from better instruments than 200 years ago… an assumption I’m increasingly doubting 😉 and on a “gap” in the data FD will “reset” to zero while I just keep waiting for a valid data item and put all the change in when one shows up. (Thus, if there were a 0.1 C / year trend and a 5 year gap I’d expect to put +0.5 C in when that valid datum shows up…). This preserves trends better especially with holey data (and the data has a LOT of holes in it…) while not having any “imagined data” used at all. (No in-filling or any other estimating… just pure data, straight.) Oh, and I calculate the anomaly by “month of the year” comparing only January to January and December to December. This avoids all the fallacies of averaging temperatures (an intensive variable…) and all the issues around data dropouts being different by season et.al.
In exchange for all that, it loses the “feature” of automatically adjusting to things like equipment changes. (That ‘reset on change’ behaviour…). To me, this is an added feature as I want to see the impact of the changes. That’s what I’m hunting for: “What Changed” that could bugger things?
OK, enough preamble. For North America the graph is:
http://chiefio.files.wordpress.com/2010/03/north_america_sh.png
This is ALL DATA, so those early years are ONE thermometer (that’s why it’s so volatile in the early years…). About 1816 you can clearly see the Year Without A Summer. By about 1825 we’re up to around 25 thermometers and the graph settles down at the zero line. We have a cold period, then a very hot 1934 era, cooling into the GIStemp baseline, and warming back to the 1998 peak. Ending back at zero today. (“Now” is defined as zero as I run time from now into the past…) So, with all that preamble, notice that about 30 years ago we are again at zero? Oh, and this is the GHCN “unadjusted” data and has, in theory, no UHI adjustments done to it at all.
Basically, this was as close to the “raw” data as I could figure out how to get and as close to “pure anomaly” as I could create. All “unadorned”.
To make it easier to see exactly, here is the text version of the data for The USA which is most of the data used to make that chart. You want the “dT” line. I’ve cut this table at 1800 when we have all of 4 thermometers, as before that I don’t see much real usability for averaging reports. You will notice that 1979 and 1980 are both -0.38 while 1978 is +0.34 which to me says “pretty much zero +/- annual jitter”.

[chiefio@Hummer DTemps]$ cat Temps.rM425.yrs.dT
Produced from input file: ./DTemps/Temps.rM425
Thermometer Records, Average of Monthly dT/dt, Yearly running total
by Year Across Month, with a count of thermometer records in that year
-----------------------------------------------------------------------------------
YEAR     dT dT/yr  Count JAN  FEB  MAR  APR  MAY  JUN JULY  AUG SEPT  OCT  NOV  DEC
-----------------------------------------------------------------------------------
2010  -0.02  0.02  134   0.2  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
2009  -0.08  0.06  136  -0.7  1.1  0.5  0.0  0.8 -0.3 -0.7 -0.1  0.0 -0.8  1.6 -0.7
2008   0.66 -0.73  136  -0.7  1.4 -2.2  0.2 -1.4  0.0 -0.1 -1.5 -0.8 -2.2 -0.6 -0.9
2007   0.89 -0.23  134  -3.3 -2.0  1.6 -2.1  0.4  0.0 -0.9  0.8  1.9  2.4 -0.1 -1.5
2006   0.68  0.21 1177   3.6 -1.5  0.3  0.1  0.1  0.0  0.0  0.0 -0.2 -0.1  0.0  0.2
2005   0.48  0.20 1214   1.5  2.1 -2.5 -0.1 -1.6  0.9  1.0  1.6  0.8  0.0  0.1 -1.4
2004   0.45  0.03 1382  -1.0  0.4  1.7  0.3  1.0  0.1 -1.0 -1.9  0.5  0.0  0.5 -0.2
2003   0.68 -0.23 1411  -2.0 -1.9  1.6 -0.8  0.7 -1.4 -0.4  0.6 -1.2  1.7  0.4  0.0
2002   0.78 -0.10 1421   1.9  1.0 -0.4  0.0 -1.6  0.8  0.8 -0.4  1.1 -0.9 -2.6 -0.9
2001   0.51  0.27 1434  -0.7 -2.7 -2.8  0.7 -0.2  0.0  0.3  0.2 -0.1 -0.6  4.4  4.7
2000   0.93 -0.43 1429  -0.2  0.2  2.2  0.0  1.2  0.1 -0.8  0.3  0.3  0.5 -4.5 -4.4
1999   1.27 -0.33 1447  -1.0 -0.1  0.0  0.3 -1.3  0.0  0.0 -0.4 -2.2 -0.6  1.4 -0.1
1998   0.08  1.18 1428   2.8  1.4 -1.3  1.6  2.3  0.1  0.9  1.1  1.4  0.7  2.2  1.0
1997  -0.17  0.26 1431   0.2  0.9  2.9 -1.1 -1.0 -0.6  0.1 -0.3  1.2  0.1  0.7  0.0
1996   0.34 -0.52 1464  -1.9 -0.5 -2.6  0.2  0.7  0.6 -0.5 -1.2 -0.3 -0.5 -0.9  0.7
1995   0.38 -0.04 1495   2.5  2.2  0.0 -1.6 -0.6 -1.3  0.3  1.4 -0.3  0.2 -1.4 -1.9
1994  -0.26  0.64 1519  -1.4  0.4  1.5  1.5 -0.3  1.6  0.1 -0.4  0.9  0.7  1.9  1.2
1993   0.28 -0.53 1529  -1.4 -4.5 -1.6 -1.0  0.2  0.2  0.8  1.4 -0.5 -0.4 -0.4  0.8
1992   0.80 -0.53 1536   2.0  0.0  0.0 -0.7 -1.2 -1.2 -1.3 -1.7 -0.3 -0.5  0.4 -1.8
1991   0.94 -0.14 1549  -3.9  1.2 -0.4  0.5  2.1 -0.1  0.3  0.3 -1.0  0.2 -2.8  1.9
1990  -0.10  1.04 1572   1.4  3.8  1.3  0.0 -0.6  0.8 -0.2  0.5  1.4 -0.1  1.3  2.9
1989   0.33 -0.43 1597   3.4 -2.0 -0.1  0.2 -0.6 -1.0 -0.4 -1.3 -0.2  1.0 -0.5 -3.7
1988   0.79 -0.46 1598  -1.8 -1.9 -0.3 -0.8 -0.9  0.0  0.5  1.0 -0.2  0.1 -0.4 -0.8
1987   0.72  0.07 1589  -1.1  0.9 -1.1  0.1  0.9  0.1  0.0  0.4  0.3 -1.2  1.4  0.2
1986  -0.20  0.92 1590   3.8  2.1  0.8 -0.6 -0.3  1.3  0.0  0.1  0.5 -0.1  0.3  3.1
1985   0.16 -0.36 1594  -1.1 -3.4  2.2  2.3  1.2 -0.7  0.4 -1.2  0.1 -0.1 -0.4 -3.6
1984   0.13  0.03 1592  -2.2  0.5 -1.6  1.0  0.7  0.7 -0.7 -1.0 -1.2  0.0 -0.8  5.0
1983  -0.14  0.27 1594   3.8  2.1  0.3 -0.5 -1.7  0.6  0.5  1.9  0.7  0.6  0.6 -5.7
1982   0.60 -0.74 1605  -3.1 -2.2 -0.4 -3.3  1.2 -1.9 -0.2 -0.1 -0.2  0.6 -1.2  1.9
1981   0.23  0.37 1623   0.2  2.5  1.7  1.8 -0.7  0.8 -0.8 -0.7 -0.9 -0.1  1.1 -0.5
1980  -0.38  0.62 1650   4.2  2.6 -1.6  0.6  0.5  0.2  1.3  1.0  0.2 -1.2  0.5 -0.9
1979  -0.38 -0.00 1657  -1.4 -0.2  1.0 -0.6 -0.1 -0.5 -0.3 -0.4 -0.1  0.5 -0.4  2.5
1978   0.34 -0.73 1660   1.5 -4.4 -1.8 -1.3 -1.2 -0.4 -0.4  0.1  0.1  0.2 -0.3 -0.8
1977  -0.19  0.53 1660  -3.6 -1.9  0.6  1.1  1.9  0.8  0.9  0.5  0.9  1.7  2.4  1.1
1976  -0.02 -0.18 1669  -1.7  3.3  2.2  2.6 -1.3  0.1 -0.3 -0.6  0.7 -2.6 -2.7 -1.8
1975   0.25 -0.27 1670   0.5 -0.6 -2.8 -2.5  0.6  0.0 -0.2  0.6  0.1  0.7  0.4  0.0
1974   0.41 -0.16 1679   0.7  0.4 -0.6  1.2  0.5 -0.7  0.3 -1.0 -1.2 -1.4 -0.2  0.1
1973  -0.22  0.63 1685   0.0  0.2  1.3 -0.3 -0.8  0.8  0.5  0.4  0.0  2.1  1.9  1.5
1972   0.03 -0.26 1689   0.8 -0.2  1.9  0.0  1.2 -1.0  0.0  0.0 -0.2 -2.1 -1.2 -2.3
1971   0.02  0.02 1693   1.1 -0.8  0.1  0.0 -1.7  0.5 -0.8 -0.6  0.0  1.6  0.0  0.8
1970  -0.11  0.13 1797  -1.3  1.0  1.1 -1.2  0.1  0.7  0.0  0.2  0.0  0.6  0.1  0.2
1969  -0.12  0.02 1813   0.0  0.4 -3.7  0.7  1.6 -0.6  0.5  0.7  0.5 -1.4  0.0  1.5
1968   0.05 -0.17 1822  -2.4  0.0  0.2  0.0  0.4  0.2  0.4  0.1  0.2  0.4  0.0 -1.6
1967  -0.09  0.14 1823   3.6  0.1  0.2  0.8 -1.1  0.0 -1.3  0.0 -0.2  0.5 -1.1  0.2
1966   0.05 -0.14 1830  -2.4 -0.2  3.4 -1.0 -1.1  0.6  1.3 -0.1  0.7 -0.8 -0.5 -1.6
1965   0.08 -0.03 1835  -0.4  0.1 -1.6  0.2  0.0 -0.7 -1.0  0.3 -0.5  0.5  0.7  2.0
1964   0.21 -0.12 1841   3.6 -0.4 -2.3 -0.3  0.4 -0.2  0.4 -0.6 -0.9 -3.1 -0.5  2.4
1963   0.13  0.07 1850  -1.4 -1.1  2.8  0.2 -0.8  0.3  0.6 -0.2  0.9  1.7  0.6 -2.7
1962   0.18 -0.05 1801  -1.0 -1.2 -2.5  1.5  2.0 -0.5 -0.6 -0.2 -0.4  0.9  0.7  0.7
1961  -0.02  0.20 1761  -0.4  2.9  4.3 -2.2 -0.4  0.1 -0.1  0.1 -0.9 -0.3 -0.9  0.2
1960   0.41 -0.42 1763   0.5 -0.8 -3.2  0.5 -1.0 -0.5  0.0 -0.6  0.5  0.5  1.9 -2.9
1959   0.21  0.20 1767  -1.1  0.9  1.3  0.4 -0.2  1.0  0.2  0.2  0.0 -0.4 -2.2  2.3
1958   0.45 -0.24 1769   2.1 -3.4 -1.6 -0.2  0.9 -0.7 -0.6  0.7  0.4  1.3  1.0 -2.8
1957   0.41  0.04 1765  -1.6  2.6  0.5  1.0 -0.6 -0.3  0.6 -0.1  0.0 -2.1  0.3  0.2
1956   0.20  0.21 1757   0.3  0.5  0.0 -1.9  0.0  1.8 -1.0 -1.0 -0.6  0.6  1.1  2.7
1955   0.90 -0.70 1755  -0.4 -4.6  0.4 -0.4  1.6 -1.4  0.0  0.9 -0.3 -0.1 -2.9 -1.2
1954   1.07 -0.17 1828  -2.7  2.1 -2.0  2.5 -1.0 -0.6  0.5  0.0  0.2 -0.6  0.0 -0.4
1953   0.51  0.56 1818   1.8  0.2  2.6 -1.2  0.0 -0.2 -0.1 -0.1  0.1  1.9  1.5  0.2
1952  -0.07  0.58 1804   0.9  0.7 -0.2  0.8 -0.3  2.1  0.5  0.3  0.8 -1.0  1.4  1.0
1951  -0.02 -0.06 1791  -0.7 -0.1  0.0  0.7  0.5 -0.5  1.1  0.8  0.2 -1.4 -1.1 -0.2
1950   0.50 -0.52 1763   1.9  0.9 -1.2 -1.7 -1.0 -0.7 -1.6 -1.0  0.0  1.2 -2.2 -0.8
1949   0.18  0.33 1754   0.3  0.5  0.8 -0.7  0.8  0.3  0.5  0.1 -1.1  0.9  1.1  0.4
1948   0.38 -0.20 1624  -2.3  0.2  0.3  1.0  0.0  0.8  0.3 -1.2 -0.2 -2.9  1.7 -0.1
1947   0.92 -0.54 1507   0.1 -2.1 -4.5 -1.4  0.4 -0.7 -0.5  2.0  0.9  2.5 -1.9 -1.3
1946   0.30  0.62 1485   0.9 -0.1  0.4  1.5  0.4  1.0  0.4 -0.7 -0.3  0.1  0.2  3.6
1945   0.45 -0.15 1482  -1.4  0.1  4.1  1.0 -2.5 -1.5 -0.1 -0.1  0.0 -0.6  0.0 -0.8
1944   0.38  0.07 1435   1.9 -0.5  0.0 -1.5  1.5 -0.3 -0.6 -0.5  0.7  0.7  0.5 -1.0
1943   0.40 -0.03 1427  -0.5  2.5 -1.7 -0.9  0.0  0.7  0.1  0.8  0.0 -0.6 -1.0  0.3
1942   0.84 -0.44 1430  -0.8 -1.0  1.6  0.2 -1.4  0.0  0.0 -0.2 -0.7 -0.5 -0.1 -2.4
1941   0.18  0.67 1423   4.6 -0.3 -0.9  1.7  1.1 -0.5  0.0  0.2  0.0  0.0  2.0  0.1
1940   1.00 -0.83 1415  -5.8  1.3 -0.9 -0.4 -1.1  0.1 -0.2 -0.4 -1.2  0.5 -1.3 -0.5
1939   1.05 -0.05 1416   1.3 -2.7 -1.6 -0.6  1.4  0.2  0.3 -0.7  0.6 -0.9  0.5  1.6
1938   0.17  0.88 1415   2.2  2.3  3.5  1.1 -0.9 -0.2 -0.2 -0.3  0.5  1.6  0.1  0.9
1937   0.41 -0.24 1414   0.1  3.8 -2.6  0.1 -0.9 -0.7 -1.2  0.0 -0.7 -0.1  0.5 -1.2
1936   0.26  0.15 1412  -2.1 -5.8 -0.1  0.1  3.3  1.5  0.7  1.0  1.0  0.0 -0.2  2.4
1935   1.27 -1.01 1407  -2.0  1.6  1.2 -1.9 -3.7 -2.0 -0.5 -0.1  0.2 -1.2 -2.6 -1.1
1934   0.86  0.41 1406  -0.2  1.0  0.2  1.5  2.1 -0.2  0.9  0.7 -1.8  0.8  1.7 -1.8
1933   0.33  0.53 1408   1.2 -3.4  2.2 -0.5  0.0  1.2  0.3 -0.4  1.5  0.8  0.7  2.7
1932   1.28 -0.95 1404  -0.1 -0.5 -1.1  0.0  0.6 -0.6 -0.9  0.3 -2.0 -1.9 -2.2 -3.0
1931   0.39  0.89 1398   4.7 -0.8 -0.6 -1.3 -0.1  1.1  0.4 -0.4  1.2  2.5  1.6  2.4
1930  -0.25  0.64 1391  -1.1  6.5 -1.8  1.2  0.5  0.5  0.8  0.4  1.2 -0.9  1.3 -0.9
1929   0.27 -0.52 1386  -2.9 -3.7  0.8  1.8 -1.1  0.9  0.0  0.0  0.4 -0.8 -1.6  0.0
1928   0.38 -0.11 1383   0.4 -1.9  0.0 -1.7  0.8 -0.6  0.4  1.6 -1.3 -0.3 -1.1  2.4
1927   0.22  0.16 1376   0.2  0.5  1.7  0.8 -0.8 -0.4 -0.4 -1.9  0.5  0.8  2.0 -1.1
1926   0.51 -0.29 1370   0.9 -0.4 -2.5 -2.6  0.9 -1.3 -0.1  0.5 -1.5  3.1  0.0 -0.5
1925  -0.50  1.01 1367   1.1  1.9  3.3  2.3  1.0  1.0  1.1 -0.1  2.6 -3.4 -0.9  2.2
1924   0.18 -0.68 1360  -3.9  2.5 -0.4  0.1 -0.9 -0.1 -1.1  0.3 -1.5  1.8 -0.2 -4.8
1923   0.43 -0.25 1360   3.4 -1.6 -1.4 -0.4 -1.2 -0.9  0.4 -0.6 -1.0 -1.9  0.2  2.0
1922   1.27 -0.83 1353  -3.4 -2.2 -3.1 -0.4  0.6 -0.1 -1.1  0.3  0.0  0.1  0.0 -0.7
1921  -0.17  1.43 1350   2.9  2.0  3.3  2.5  0.7  1.5  1.3  0.5  0.8 -0.2  1.2  0.7
1920   0.18 -0.35 1342  -1.9  0.0 -0.5 -2.1 -0.4 -1.0 -1.0 -0.6  0.0  0.8  0.0  2.5
1919   0.18  0.00 1334   4.8 -0.2 -1.9  1.1 -0.8 -0.3  1.1 -0.6  2.2 -1.2 -0.7 -3.5
1918  -0.90  1.08 1326  -3.0  2.0  3.3  0.0  3.3  1.8 -0.9  1.1 -1.1  3.4 -0.7  3.8
1917  -0.32 -0.58 1313  -0.1 -1.3 -0.9 -0.7 -2.4  0.4 -0.2 -0.7  0.0 -1.5  1.0 -0.6
1916   0.03 -0.35 1297   0.1 -2.1  2.0 -2.5  0.7 -0.1  1.8  1.3 -0.7 -1.5 -1.3 -1.9
1915   0.17 -0.13 1276  -2.9  3.8 -1.7  2.1 -1.6 -1.9 -1.5 -1.1  0.3  0.0  0.0  2.9
1914   0.24 -0.07 1264   1.5 -0.4  0.3 -0.3  0.7  0.5  0.3 -0.8  0.2  1.8 -0.9 -3.8
1913  -0.54  0.78 1252   4.2 -0.3  1.8  0.1 -0.4  1.0  0.4  1.7  0.0 -0.8  1.5  0.2
1912   0.42 -0.96 1234  -4.7 -2.0 -3.9  0.6 -0.5 -2.1 -0.2 -0.4 -1.3  0.3  2.5  0.2
1911   0.32  0.10 1222   1.5  2.4 -3.1 -1.7  1.8  1.6 -0.3  0.0  0.4 -1.4 -1.5  1.5
1910  -0.05  0.37 1191  -1.1 -3.0  4.8  2.2  0.1 -0.5  0.8 -0.9  0.6  1.8 -2.7  2.3
1909   0.35 -0.40 1179  -0.1  1.5 -1.9 -1.7 -0.6  0.8 -0.5  1.0 -0.8  0.0  1.4 -3.9
1908  -0.12  0.47 1150   0.4 -0.5 -0.7  3.0  1.8  0.6  0.2  0.0  0.7 -0.2  0.9 -0.5
1907   0.17 -0.29 1127  -1.1  0.3  4.9 -3.3 -2.1 -0.9  0.3 -0.8 -1.2  0.1  0.0  0.3
1906  -0.26  0.42 1097   3.9  3.4 -4.7  1.2  0.0 -0.2  0.0  0.2  0.5  0.2 -0.5  1.1
1905  -0.45  0.19 1067  -0.2 -1.9  1.9  1.1 -0.2  0.5  0.3  0.8  0.5 -0.8 -0.1  0.4
1904  -0.37 -0.07 1054  -2.0 -0.2 -1.2 -1.0  0.0  0.9 -0.4 -0.1  0.9 -0.1  1.4  0.9
1903   0.02 -0.39 1016   0.3  0.1  0.7  0.0 -0.9 -1.1 -0.2 -0.2  0.2 -0.3 -2.7 -0.6
1902   0.06 -0.04  983  -0.9  0.5  1.0  0.5  0.9 -0.8 -2.0 -1.0 -0.6 -0.3  2.3 -0.1
1901   0.56 -0.50  966  -0.8 -0.4  0.6 -1.3 -0.6 -0.1  1.7 -0.5 -1.1 -1.1 -0.7 -1.7
1900  -0.16  0.72  952   1.7  1.9  1.0  0.5  0.5  0.1  0.1  0.9  0.9  1.3 -1.5  1.2
1899   0.11 -0.27  927  -0.9 -4.6 -2.7  0.5  0.2 -0.1 -0.2 -0.1 -1.0  1.5  3.1  1.1
1898   0.08  0.03  909   1.8  0.4  1.4 -0.6  0.0  0.7 -0.2  0.8 -0.2 -1.8 -1.1 -0.9
1897   0.21 -0.12  886  -1.3 -0.2  1.1 -1.1 -1.6 -0.5  0.1 -0.7  1.9  2.1  0.5 -1.8
1896  -0.40  0.61  851   2.0  3.8 -1.1  0.3  1.5  0.0  1.0  0.2 -1.8  0.6 -0.1  0.9
1895   0.25 -0.65  822  -1.8 -1.5 -1.7  0.3 -0.1  0.0 -1.0  0.0  0.5 -1.8  0.0 -0.7
1894  -0.48  0.73  765   2.5 -0.5  2.4  1.0  1.0  0.0  0.0  0.5  0.6  0.5  0.2  0.6
1893  -0.30 -0.18  728  -1.1 -2.2  0.0  0.0  0.1  0.1  0.5 -0.3 -0.2 -0.3  0.0  1.2
1892  -0.01 -0.29  662  -2.0  1.0  0.6 -1.1 -0.4  0.1  0.7  0.2 -0.6  0.6  0.0 -2.6
1891   0.16 -0.17  595   0.2 -0.9 -0.6  0.1 -0.2 -0.7 -1.3  0.4  1.4  0.0 -1.5  1.1
1890   0.12  0.04  536   0.5  2.4 -1.8 -0.4 -0.2  1.1  0.2 -0.3  0.0  0.4  1.4 -2.8
1889  -0.37  0.48  511   2.9 -1.3  2.3  0.3  0.5 -0.5 -0.2 -0.2  0.0  0.0 -0.4  2.4
1888  -0.17 -0.20  441  -1.2  0.3 -1.3  0.5 -1.9 -0.1 -0.6  0.2 -0.1 -0.1  0.5  1.4
1887  -0.30  0.13  386   0.9  0.1  0.4 -0.4  0.8  0.4  0.7 -0.6 -0.6 -1.1  0.4  0.6
1886  -0.52  0.22  356  -0.7  1.8  0.6  0.6  1.0  0.0 -0.2  0.8  0.6  1.0 -1.1 -1.8
1885  -0.35 -0.17  335   0.3 -2.0 -0.8  0.7 -0.3  0.0  0.8  0.0 -0.8 -1.8  0.4  1.5
1884  -0.47  0.12  316   0.0  0.6  0.5 -0.6  0.9 -0.6 -0.6  0.0  1.2  1.3 -0.2 -1.1
1883   0.03 -0.50  277  -2.3 -2.3 -1.6  0.0  0.3  0.4  0.5 -0.6 -0.6 -1.2  0.7  0.7
1882   0.20 -0.17  259   2.3  2.1  1.1  0.2 -2.6  0.0 -1.0 -0.6 -1.1  0.1 -0.1 -2.4
1881  -0.09  0.29  212  -6.2 -1.5  0.0 -0.6  0.0 -0.5  0.6  0.9  1.9  1.5  3.0  4.4
1880   0.04 -0.13  203   5.2  1.7 -1.8  0.2  1.0  0.6 -0.6  0.0  0.2 -2.6 -3.5 -2.0
1879   0.43 -0.39  201  -2.0 -2.2 -1.5 -1.6  0.9  0.2 -0.4 -0.9 -0.6  2.1 -0.5  1.8
1878   0.20  0.23  192   1.7 -0.4  3.3  1.9  0.1 -0.5  0.5  0.4 -0.4 -0.2  0.5 -4.1
1877  -0.39  0.59  179  -3.1  1.5  1.0  0.0 -0.4 -0.6 -0.2  0.0  1.1  1.4  0.3  6.1
1876  -1.01  0.62  177   5.6  4.2  0.1  1.1 -0.3  0.7  0.5  0.7  0.0 -0.6  0.8 -5.4
1875  -0.11 -0.90  176  -4.2 -3.1 -1.4  0.9 -0.2 -0.8 -0.4 -0.4 -1.2 -0.6 -0.9  1.5
1874  -0.65  0.54  172   2.2  0.6  0.2 -1.5  1.0  0.1  0.3  0.0  1.0  1.2  1.1  0.3
1873  -0.64 -0.01  171  -0.7 -0.6  1.5 -0.9 -0.9  0.0 -0.2 -0.5 -0.4 -0.8  0.7  2.7
1872   0.02 -0.67  163  -1.1 -0.8 -3.9 -0.8  0.0 -0.1  0.5  0.0  0.9 -0.8 -0.5 -1.4
1871   0.10 -0.07  146  -0.6  0.2  3.0  0.4 -0.2 -0.1 -0.9  0.6 -1.1  0.2 -1.9 -0.5
1870  -0.64  0.74  125   0.2 -0.4 -0.2  1.0  1.5  1.5  1.2  0.0  0.8  2.5  1.6 -0.8
1869  -0.85  0.21  103   4.0  2.2 -2.8  0.9 -0.2 -0.5 -2.1  0.5  1.1 -1.4 -0.9  1.7
1868  -0.53 -0.32   97  -0.3 -3.0  4.2 -1.2  1.3 -0.7  2.2 -0.7 -1.7 -1.3 -1.3 -1.3
1867  -0.46 -0.07   88  -1.4  2.1 -1.8 -1.3 -1.3  0.7 -1.1  1.7  1.1  0.0  0.2  0.2
1866  -0.08 -0.38   82   0.7 -0.9 -1.8  0.7 -1.1 -1.1  1.5 -1.0 -2.8  0.9  0.2  0.2
1865  -0.38  0.29   75  -1.5 -0.7  1.3  1.4 -0.4  1.0 -1.4 -0.9  2.7  0.7  0.9  0.4
1864  -0.48  0.10   77  -1.6  0.9  0.7 -0.6  0.2  1.2  0.8  0.3  0.3  0.0 -0.4 -0.6
1863  -0.57  0.09   67   2.7  1.4 -0.6  0.6  0.8 -0.4  0.0 -0.1 -1.1 -1.6  0.6 -1.2
1862  -0.21 -0.36   64  -0.1 -3.0 -0.7 -1.2  0.9 -1.1  0.3  0.7  0.6  0.0 -0.5 -0.2
1861  -0.14 -0.07   71  -1.4  1.1 -1.5  0.3 -1.8  0.1 -0.5 -0.4  0.3  0.5  0.0  2.5
1860  -0.48  0.33   59   0.3 -0.6  0.0  1.4  0.3  0.7  0.2  0.3  0.0  1.2 -0.6  0.8
1859  -0.36 -0.12   77  -2.3  3.2  1.6 -0.9  1.6 -1.2 -0.3 -0.1 -0.6 -2.0  2.4 -2.8
1858  -0.93  0.57   81   6.8 -4.1  1.1  2.4  0.2  1.2  0.1  0.1 -0.1  1.3 -1.1 -1.0
1857  -1.10  0.17   82  -0.7  4.7  1.1 -3.3 -0.7 -1.6 -0.8  0.3  0.2 -0.5 -0.3  3.6
1856  -0.48 -0.62   74  -4.3 -0.2 -1.5  0.0 -0.7  1.9  0.6 -0.2 -0.8  0.3 -1.0 -1.6
1855  -0.06 -0.42   64   1.3 -2.3 -1.6  1.0  0.0 -1.0 -1.0 -1.0 -0.2 -1.7  0.9  0.6
1854  -0.12  0.06   59  -0.5 -0.1  0.3 -0.6  0.1 -0.5  1.1  0.3  0.4  1.4 -0.6 -0.6
1853  -0.41  0.29   46   2.4 -0.9  0.2  1.9 -0.4  0.8 -0.4  0.5  0.7 -1.6  1.9 -1.6
1852  -0.28 -0.12   56  -3.0 -0.7 -1.2 -1.4  0.4  0.4  0.3  0.0 -1.0  0.5  0.5  3.7
1851  -0.32  0.03   54  -0.7  0.8  1.7  1.1  1.3 -0.6 -0.5 -0.6  0.4  0.6 -2.1 -1.0
1850  -0.50  0.18   48   2.8  3.0 -1.2 -0.6 -0.6 -0.2  0.8  0.1  0.2  0.5 -1.6 -1.0
1849  -0.26 -0.24   48  -2.8 -2.5  1.2 -0.9 -1.8  0.3  0.3  0.0  1.0 -0.1  4.5 -2.1
1848  -0.47  0.21   42   2.3  0.3  1.1  0.5  1.5  0.8 -1.2  0.3 -1.5  0.6 -3.2  1.0
1847   0.32 -0.79   43  -1.6  1.6 -2.8 -1.7 -1.6  0.0  0.4 -1.4 -2.2 -0.2 -0.5  0.5
1846  -0.13  0.45   43  -0.5 -2.8 -0.6  0.2  1.9 -0.5 -0.1  0.4  2.5 -0.4  1.7  3.6
1845  -0.34  0.22   47   4.3  0.6  1.1 -2.0 -1.7  0.2  0.4  1.0  0.0  0.8  0.6 -2.7
1844  -0.88  0.53   50  -4.7  3.6  4.6  2.9  1.4  0.1  0.0 -0.7 -1.0  0.5  0.5 -0.8
1843  -0.10 -0.77   48   0.6 -5.5 -7.0 -1.4  0.6  0.8 -0.1  0.7  1.3 -1.4  0.5  1.6
1842  -0.42  0.32   44   0.7  3.3  3.4  2.7 -0.4 -2.3 -0.2 -0.8 -1.2  1.5 -1.4 -1.5
1841  -0.29 -0.12   43   4.0 -3.4 -1.2 -2.6 -1.5  0.9 -0.1  0.0  1.7 -1.5  0.2  2.0
1840  -0.24 -0.05   42  -3.0  1.6  1.0 -0.1  0.6  1.2 -0.3  0.6 -0.4 -1.7  0.8 -0.9
1839  -1.07  0.82   39  -1.4  5.5 -0.9  3.6  1.6 -2.5 -0.8 -0.7 -0.5  2.9  0.9  2.2
1838  -1.09  0.03   38   3.8 -4.6  2.4 -0.1 -0.4  1.8  2.1  1.3  0.9 -1.0 -3.4 -2.5
1837  -1.47  0.37   41  -1.5  2.6  0.9 -1.0 -1.7  0.1 -0.3  0.7 -0.4  2.1  1.9  1.1
1836  -0.81 -0.66   38  -0.8 -0.5 -1.9 -0.3  0.2 -0.8  0.0 -0.9  1.1 -3.7 -0.8  0.5
1835   0.27 -1.08   34   2.1 -6.7 -1.8 -2.0  0.6  0.1 -1.8 -1.2 -1.9  1.7 -0.1 -2.0
1834   0.22  0.06   36  -3.6  3.9  1.5 -0.4 -1.9  0.8  1.1  1.0 -0.3 -0.1  0.1 -1.4
1833  -0.08  0.29   34   2.0  0.2 -1.2  2.2  2.2 -0.7  0.5 -0.3  0.6 -1.2 -0.9  0.1
1832  -0.23  0.16   36   2.2  0.9 -1.3 -1.1 -1.3 -1.9 -0.8 -0.9 -0.2 -0.2  0.4  6.1
1831   0.93 -1.17   31  -1.8 -1.9  0.5 -1.3  0.5  1.6 -1.1  0.0  0.2 -0.5 -3.5 -6.7
1830   0.03  0.90   30  -0.5  2.2  2.5  2.1 -1.4 -0.7  2.0  0.2  1.1  0.8  4.3 -1.8
1829   1.28 -1.25   29  -2.1 -6.1 -2.7  0.8  0.6 -1.2 -0.7 -0.9 -1.4  0.3 -1.9  0.3
1828   0.41  0.87   24   3.9  2.1  0.5 -2.6  0.2  1.8 -0.1  1.0 -0.3 -0.3  2.3  2.0
1827   0.53 -0.12   21  -1.7  1.1  0.3  2.4 -2.3 -0.7  0.1  0.0  0.0  0.0 -1.2  0.5
1826   0.67 -0.14   17  -0.5 -0.1 -1.1 -2.3  2.6 -0.4 -1.3  0.0  0.5 -0.1  0.4  0.6
1825  -0.01  0.68   18  -1.1  1.9  2.8  1.1  0.8  1.8  1.3  0.8 -0.3  0.7  0.0 -1.6
1824  -0.58  0.57   17   1.7  2.5 -0.2 -0.6 -0.2 -0.2 -0.1 -0.7  0.8  0.5  1.4  1.9
1823  -0.03 -0.55   14   2.3 -2.9 -1.6  1.1 -1.4 -0.2  0.0  0.1 -2.0 -0.7 -2.7  1.4
1822  -0.63  0.60   13   1.0 -0.7  2.1  1.4  1.1 -0.3  1.3 -0.7  0.3  0.3  1.4  0.0
1821  -0.41 -0.22   10  -0.4 -1.2 -0.3 -3.3  0.2  0.5 -1.6  1.5  0.1  0.9  0.8  0.2
1820  -0.09 -0.32    9  -1.6  0.0  0.7  0.3 -0.4 -0.1  0.6 -0.2 -0.2  0.0 -1.6 -1.3
1819  -0.90  0.81    4   2.0  4.2 -1.6  1.1  0.0  0.0 -0.5  1.1  2.2 -0.3 -0.2  1.7
1818  -0.93  0.03    3   0.8  0.6  0.9 -3.2  0.3  2.4  1.6 -0.1 -1.1  0.6  0.5 -2.9
1817  -0.82 -0.12    4  -0.8 -3.6  0.0  0.3  0.6  0.4  1.0  0.7  1.8 -1.4 -0.6  0.2
1816  -0.39 -0.42    3   1.2  1.7 -1.0 -0.9  0.0 -2.4 -3.2 -0.6 -2.4  0.8  0.7  1.0
1815   0.27 -0.67    2  -0.3 -3.3  1.5 -1.7 -4.4  0.1  1.4 -1.2 -0.6 -0.6  0.7  0.4
1814   0.05  0.22    3   0.2  1.6  1.1  0.2  3.4 -0.4 -0.4 -1.3 -1.9  0.7 -0.2 -0.3
1813  -0.83  0.88    4   0.0  0.0  0.2  0.9  1.0  0.7  0.3  1.9  2.7  0.2  1.3  1.4
1812   0.46 -1.29    3  -1.8  0.0 -3.8 -0.3 -1.9 -1.4 -0.8 -0.5 -1.8 -2.0 -0.7 -0.5
1811   0.30  0.16    2   0.1 -3.0  3.3 -1.8 -1.3 -0.1  0.8 -0.4  0.0  2.5  1.2  0.6
1810  -0.38  0.68    2   2.4  4.4  0.2  1.2  0.5  0.8  0.9  0.8  2.5 -3.9  1.8 -3.4
1809   0.29 -0.68    2  -1.8 -3.7 -2.1 -0.5  1.2 -0.8 -2.3  0.1 -1.6  4.9 -3.0  1.5
1808  -0.24  0.53    2   1.2  2.9  3.4  1.6  0.0  1.9 -1.3 -1.7  0.3 -1.6  1.8 -2.1
1807  -0.07 -0.17    2  -1.8 -3.8  0.5  1.2 -1.9 -1.8  2.5  1.3 -1.3  0.8 -1.3  3.5
1806   1.23 -1.30    2   1.8  1.4 -5.3 -4.0  0.3 -0.1 -2.2 -2.5 -1.8  0.9  0.6 -4.7
1805   0.26  0.98    2  -0.2  0.9  3.4  2.6 -0.9 -0.4  0.9  1.0  0.4 -0.2 -1.6  5.8
1804   0.68 -0.43    3  -2.1 -1.9 -1.4 -1.6  3.2 -0.6 -0.8  0.0  3.3 -1.4  2.3 -4.1
1803   0.77 -0.09    4  -4.0  0.5 -0.3  0.0 -0.6  1.0  1.1  0.2 -1.4  0.4 -0.8  2.8
1802   0.67  0.11    4   3.7 -0.2 -0.8  2.0 -2.0 -0.3 -0.6  0.9 -0.5 -0.3  0.4 -1.0
1801  -0.00  0.67    4   1.2  2.8  2.5 -2.5  1.3 -0.1 -0.3 -0.5  1.3  1.6  0.8 -0.1
1800  -0.41  0.41    4  -1.2  0.0  2.5  1.8  0.1  0.1  0.6 -0.1  0.2  0.6 -1.5  1.8
For Country Code 425

For folks wanting to see more about North America, the link to that particular article is:
http://chiefio.wordpress.com/2010/03/16/north-america-flat-canada-not-so-much/
So at the end of the day, it looks to me like Dr. Spencer has it “just right”.

suricat
March 16, 2010 6:19 pm

Now I see what you are doing Roy.
Well done. You’ve given credibility to the data by using a 1 km grid square for population density. The temp data is just as chaotic for sample area, but the UHI effect associated with population density is proved for a station’s report of temp.
Best regards, suricat.

HumanityRules
March 16, 2010 6:41 pm

Interesting work. As with all data on climate change it tends to raise more question than answers in me. Whether it confirms or contradicts my own personnal outlook. I’m hoping that’s the correct position to take as a true sceptic.
So I have some questions which either the author or somebody in the know might answer for me.
1) Has anybody in the literature taken this sort of approach to UHI? What were their conclusions?
2) Can you do the same analysis using the satellite data? Would that back-up the findings?
3) Can you analyse temporal changes with this method?
4) What are the limitations/weaknesses of this method compared with other attempts to calculate UHI? (sorry bit of a hard one)

James Sexton
March 16, 2010 6:46 pm

Dr. Spencer,
What a serendipitous occasion to come back from a vacation and have such a clear analysis to read. Thanks!! However, your earlier article did show a slight variance from the U.S. and the rest of the world especially in the lower pop density. So, I think a “rest of the world” analysis would be in order to be able to apply your findings globally. Given the clear methods you’ve given us, I’d do it myself, but then I’d be stealing your thunder. (lol, not really, I’m lazy and am more comfortable critiquing and cheering from the sidelines.) Also, (and I’m not sure if you’re planning to go this way or not) a time study should probably be done going further back, comparing UHI per points in history. Obviously, there are significant differences in UHI components in a town of 50,000 in 1920 compared to one presently.
I see you’ve brought some detractors to the site. GREAT WORK!!! You know you’re onto something when that happens.
pat michaels (11:06:21) :……I’d say that his would probably show a more significant divergence. I’d also say, that your right about the oceans, but that’s not what was measured to come to the warming conclusion(not significantly) by the Jones’ and Hansen’s ect. and then I’d probably commence about the oscillation of currents and winds ect., but then, I’ll probably wait for the final chapter of Dr. Spencer’s analysis.
Dr. Spencer, thanks again.

HumanityRules
March 16, 2010 6:50 pm

Few more questions I forgot about.
5) There is usually a call for confidences and uncertainty measurements with this sort of data. Any available?
6)From your update “If the population densities for a pair of stations are exactly the same, we do not include that pair in the averaging.” Would these not act as good ‘negative controls’ for this experiment?

Pete H
March 16, 2010 6:51 pm

Roy,
Living in Shanghai, A thought occurred to me as I read your article. There is a huge migrant worker population, living and working in the major Chinese cities. Every Chinese New year the vast majority of these workers leave the UHI cities and travel home to the provinces for at least 2 weeks.
Indeed, Shanghai became a ghost town just a few weeks ago. Traffic levels drop to probably 10/20% of normal levels. House heating is turned off (Mostly reversible air conditioner units). Factories shut down etc.
I was thinking along the line that there must be a change in the temperature of the major conurbations taking into account the UHI effect. Has anyone noticed this or recorded any difference?

Dave N
March 16, 2010 7:18 pm

David44 (12:43:50) :
It’s apparent from the video that the selection criteria were sites over 150,000 in population where there is a relatively close rural comparison site, and both sites have full temp history over a century. As far as I can tell, they included all sites that fit that criteria.

Ray R.
March 16, 2010 7:23 pm

A VS is putting on a great show…..link to comments. Eli and Tamino roughed up a bit.
http://ourchangingclimate.wordpress.com/2010/03/01/global-average-temperature-increase-giss-hadcru-and-ncdc-compared/#comment-1626

Paul K2
March 16, 2010 7:24 pm

Hmmm… spurious warming for 2% of the surface of the planet, IF you believe this rather labored and complicated reasoning is statistically significant.
OTOH we see record breaking month after month high global temperature anomalies for the UAH record. Getting back to reality Dr. Spencer, it is looking increasingly likely that the trailing 13 month UAH moving average will set an all-time record when the July results come in. January was the highest in the database (before adjustment), February was just a whisker off Feb ’98, and March is looking like this month will blow out the Mar ’98 record. With El NIno lag effects expected to continue to show up in the atmosphere through June, the message is becoming extremely clear…
… the UAH record is confirming the GISS record and near-record highs set in the last five years. The planet’s atmosphere is heating up at a high rate. The oceans are absorbing massive amounts of heat that swamp the heat buildup in the atmosphere.
Isn’t this the important message?