In this post I present a 5 panel trend charts which show the year-to-date anomaly trends for the 5 major global temperature anomaly series and a table that shows how 2010 YTD ranks over the entire record for each series. The source data and RClimate script file links are provided.
Click to Enlarge
The 2010 GISS YTD is the highest in the 131 year record. Both of the satellite anomaly series, RSS and UAH, and the Hadley and NOAA series had their 2nd highest YTD anomalies in 2010, 1998 had the highest YTD for these 4 series.
Please not that the plots show the climate agency monthly anomaly values with respect to each agency’s baseline period. For RSS and UAH, that is 1979-1998. For GISS it is 1951-1980 and for NOAA and Hadley, it is 1961-1990.
RClimate Script Details
Here are the data and RClimate Script links:
- Data Link (Use this link if you want to download CSV data file. See example script below if you’d like to read file directly into your R script.)
- RClimate Script link
I do not really know much about this stuff, but it appears to me that NOAA is saying January through October, 2010 is the warmest on record:
http://lwf.ncdc.noaa.gov/sotc/?report=global&year=2010&month=10#year-to-date
When I see long-term trends like this, I often wonder what this would look like with the measurement error overlaid for each point (maybe as a grayed and alpha-blended “ribbon”). Do you know if the measurement error for historical measures of the temperature anomaly is known?
Works.
1) In line 16: temps <- subset(temps, temps$yr<= cur_yr & temps$mo <= which_mo) the objects ‘cur_yr’ and ‘which_mo’ do not exist (yet).
To run, I replaced ‘cur_yr’ with ‘last_year’ and moved line 16 to after line 19: which_mo <- temps$mo[length(temps$mo)]
2) The function ‘addtable2plot()’ requires plotrix package.
Joel
Thanks for the quick double check of my script.
I’ve made the noted changes.
Can you retry to make sure I got them ?
Thanks