Monday, February 9, 2015

Breakdown of effects of GHCN adjustments.

In a previous post, I showed how running TempLS with unadjusted or adjusted GHCN station data made a small difference to the final global average. That is useful to know, because some say that AGW is a creature of GHCN adjustments. This is certainly not true.

That analysis can be extended. I can trace the contribution to this small adjustment effect from various categories of stations - by continent, rural status, or airport status. And as before, it shows the difference by year, or effect on trend (to present time) by year.

A useful addition to the present post will be some trends calculated for land only in a NOAA paper by Lawrimore et al, 2011, Table 4 (h/t Victor Venema):

v3 Uncorr v3 Corr
1880–20100.61°C/Century0.79°C/Century
1901–20100.70°C/Century0.91°C/Century
1951–20100.16°C/Decade0.18°C/Decade
1981–20100.27°C/Decade0.27°C/Decade

Note the change of units (in the original). I got an approx land only difference of 0.17°C/Cen for 1900-2014 by dividing the global value by the land area fraction. That was pretty rough; the value in this post is more accurate, and agrees well.

In this post, I'll start with the usual active plot - see here for guidance. I've added a new "Norm" button. This allows you to toggle to a mode where the average adjustment for each category is calculated. The initial mode gives the contribution of that category to the global average. I think that is most useful, as it is area weighted. The category average can show high average adjustments for relatively rare stations. I'm using TempLS mesh exclusively; the grid style allows too much mixing between categories.

As usual with the active graph, you can select the set to show, move and zoom the graph, and show the numerical data. Then I'll show some ordinary graphs comparing the categories.

Here is the active plot. You can click on Trendback to togle a mode that shows trends to present, and on Data to pop up the numbers you see on screen.

Global with urban/rural breakdown

Here is the plot. All have a slight smoothing over 3 years. Note that all these plots can be generated on the zoomable active plotter, along with the associated trend graphs.


So you can see that the total global contribution is negative to about 1970, then positive. This has the effect that the contribution to trends from about 1960 onward is actually cooling. These are the trends that are vigorously debated as possible evidence of AGW. They are not enhanced by adjustment - on the contrary.

Before 1970 the adjustments "cool the past", by up to 0.05°C. However, on a land basis, that is up to 0.2°C. That is representative of the maximum average adjustment.

A virtue of the left style is that it is additive. So urban adjustments make up the larger part, but not by all that much.

There is an interesting sharp rise to about 1975, and then a decade plateau. We'll see more of it. It seems here associated with rural. My theory is that it is an airport effect, and may be associated with the move of some rural classified stations to airports.

There is also a WW2 effect, principally associated with Mixed (between rural and urban). Again I think it may have been connected with the moving of civilian stations to military airports, and maybe back again.

Breakdown by airport status




Here both of those features are markedly associated with airport status (hence my theory). Despite what you sometimes hear, adjustments to airports are a relatively minor net contributor. Note that the GHCN inventory ratings that I use refer to current conditions. Many "airports" have century plus records; they obviously weren't always airports. Some adjustments probably relate to the time when the station was originally transferred to the airport.

Breakdown by continent





Well, as you might expect from recent fusses, South America has the most negative adjustments. Africa is mostly positive. Less adjusted are Europe and Oceania (BoM and NIWA take a bow!).

Special cases


Finally, I'll show CONUS (continental US) and the Arctic. There has been a lot of fuss in recent times. The Arctic adjustments have been passionately excoriated by Booker of the Telegraph.



Well, as you'll see, the Arctic has made very little contribution to the global average, and it is quite flat. On the right, the average adjustment is also quite low (and again flat in trend).

CONUS, on the other hand, has the highest average adjustments (also noted here). That is mainly TOBS. Because CONUS is a smallish fraction of total land, the effect overall is modest, and also not so strong in trend. As previously noted, the adjustments reach a minimum in the 1930's.








21 comments:

  1. OT - Nick - NCEP has not updated on your page since January 26.

    ReplyDelete
  2. Do you think TempLS is about right?

    ReplyDelete
    Replies
    1. JCH,
      No, the missing January data is still missing. I just re-ran hoping that Monday would have yielded something, but no.

      Delete
    2. I've fixed NCEP and a related problem in the daily TempLS mesh. Everything is there now except China. January is very similar to December.

      Delete
  3. Replies
    1. Thanks. Just a little warmer than Dec (0.72).

      NCEP has shot up in the last few days.

      Delete
    2. Yes. The way it is going it looks to me like the 12 months ending in two weeks could reach .70C.

      Delete
  4. Nick, meant to comment on this sooner. This is another great data research tool.

    Any chance you can add the 1900-now trends in the legend?

    ReplyDelete
    Replies
    1. To make more clear what I meant, I was suggesting, for the new active plot "Effect of GHCN adjustments on global average T", adding text each item in the length (e.g., "All")e to include the °C/decade trend for that series

      This would quantify the magnitude of the effect on estimate of trend, which I see as the most important long-term number coming out of the surface temperature data.

      Delete
  5. A "download time series" button would be pretty cool tool. ;-)

    ReplyDelete
    Replies
    1. Carrick,
      Thanks. Yes, I could add that trend to the legend. Or to the data pop-up. If you show the trendback mode, that number is at the top of the display.

      The data button is meant to display whatever the numbers are for the curves on the screen, including time series if that is what you are showing. I had improved it recently to make it more consistent in doing that.

      Delete
    2. ps I see there is a bug with zooming the time axis. It's meant to work if you drag below the x-axis. Vertical zoom (y axis) works. In progress.

      Delete
    3. Nick: If you show the trendback mode, that number is at the top of the display.

      I'm guessing it's a display error with my browser then, 'cause it doesn't show up on my computer.

      I didn't notice the Data button. Very cool! Thanks. (Offered at cost, I'd suggest moving it to a different location because it's not associated with the delay, and might be easier to spot.)

      Delete
    4. Is there a "reset display" button? Sorry if I'm being blind.

      Delete
    5. Sorry to keep spamming this, but I'd suggest increasing the amount of information on the "y" axis. E.g.,

      "effect of adjustment on temperature"

      and

      "effect of adjustment on trend".

      Finally, if you really feel like it, a button that toggles "total" versus "residual"—residual" being what we're seeing now, total meaning the unadjusted temperature series with the checked corrections added to them.

      Delete
    6. Carrick,
      "'cause it doesn't show up on my computer. "
      Sorry, I meant at the top of the window that pops up when you press "Data".

      "Is there a "reset display" button?"
      There is Redraw on the left. It shouldn't be needed, though. Any drag etc resets the display.

      Re the Y axis, I was trying to convey that info in the heading. I'll experiment.

      On the drag bug, I see that it is browser specific. Firefox is OK; Chrome is OK in vitro, but something goes wrong when in the blog environment.

      Delete
    7. I was referring to having a Reset button that "resets to default values".

      Again, just a suggestion, but having text fields that wire in the axis boundaries is a good backup for when the browser is refusing to play nicely.

      Delete
    8. Carrick,
      Sorry I misunderstood there. It's not always easy to think of a natural default. There is a quick way to go to a chosen state. You may notice as the mouse wanders, faint gold numbers vary on left and right. There is a faint diagonal separating two triangles. In the lower, it shows years. Mouse near origin gives short ranges near present. Move right (E) near the axis and the start gets earlier. Move N and the end gets earlier. If you click with Shift pressed, the x axis moves to that range. The upper tri does mirror things with the y axis.

      Delete
    9. ps I see the gold numbers have become pink.

      Delete
  6. Nick: It's not always easy to think of a natural default.

    On my drawing software, for reset, I set the view boundary to the maximum range (minimum non-clipping rectangle). Typically I wire this to CONTROL-R (or COMMAND-R on a Mac).

    ReplyDelete