четверг, 4 февраля 2010 г.

Station errors in CRUTEM3

My attempts at reproducing uncertainty estimates of CRUTEM3 dataset, following methods of [Brohan et al. 2006] have encountered strange problem. My own estimates of station errors consistently yield smaller values compared to CRUTEM3_station_error file - but only for grid-boxes with several (more than one) stations. If there is only one station in the grid-box, estimates match perfectly.

Let's look at a 5-degree grid-box with bottom-left corner at 40N 145E. There are two stations: NEMURO (474200) and YUZHNO-KURIL'SK (321650). Prior to 1947, only one station provides data for the grid-box, and it is Nemuro. Station record, published by British Met Office, lists standard deviations of temperature data over the period from 1941 to 1990.

SD1 = 1.4 1.7 1.0 1.0 1.1 1.3 1.8 1.5 0.9 0.8 1.1 1.5

Following [Brohan et al. 2006], we can calculate normal error for this station (all 30 years for the normals are available) as SD / sqrt(30):

NE1 = 0.26 0.31 0.18 0.18 0.20 0.24 0.33 0.27 0.16 0.15 0.20 0.27

For calculation if station error, we need also measurement uncertainty and homogenization adjustment uncertainty. As suggested in referenced work, measurement uncertainty = 0.03 C and  homogenization adjustment uncertainty = 0.4 C. Therefore, station error is:

SE1 = sqrt(0.03^2 + 0.4^2 + NE1^2)
SE1 = 0.48 0.51 0.44 0.44 0.45 0.47 0.52 0.49 0.43 0.43 0.45 0.49

These are exactly the same values as in CRUTEM3_station_error prior to 1947 year, when there is only one station in the grid-box. No problems so far.

However, in 1947 another station appears. It is Yuzhno-Kurilsk (321650). Standard deviations of temperature data:

SD2 = 1.3 1.6 1.1 0.9 1.1 1.3 1.6 1.5 1.0 0.8 1.1 1.4

Normal error:
NE2 = 0.24 0.29 0.20 0.16 0.20 0.24 0.29 0.27 0.18 0.15 0.20 0.26

Station error:
SE2 = sqrt(0.03^2 + 0.4^2 + NE2^2)
SE2 = 0.47 0.50 0.45 0.43 0.45 0.47 0.50 0.49 0.44 0.43 0.45 0.48

[Brohan et al. 2006] says: "The grid-box anomaly is the mean of the n station anomalies in that grid box, so the grid-box station uncertainty is the root mean square of the station errors, multiplied by 1/sqrt(n)." There are two stations, so n=2. Root mean square of the station errors:

RMS = sqrt((SE1^2 + SE2^2)/2)
RMS = 0.475 0.505 0.445 0.435 0.450 0.470 0.510 0.490 0.435 0.430 0.450 0.485

Multiplied by 1/sqrt(2):

0.336 0.357 0.315 0.308 0.318 0.332 0.361 0.346 0.308 0.304 0.318 0.343

And CRUTEM3_station_error after 1947 year lists:

0.376 0.414 0.343 0.333 0.347 0.370 0.421 0.394 0.333 0.319 0.347 0.388

They are clearly correlated with my estimates, but consistently larger, and I cannot explain that. Of course, this is not limited to that particular grid-box. No, rounding errors do not explain that.

I believe that I correctly calculated station errors for individual stations, since they match those in CRUTEM3_station_error, if there is only one station in the grid-box. But I also do not see an error I probably made in combining them.

Thoughts?

11 комментариев:

Анонимный комментирует...
Этот комментарий был удален администратором блога.
Анонимный комментирует...

Well your article helped me terribly much in my college assignment. Hats high to you dispatch, intention look audacious for more cognate articles without delay as its anecdote of my pet subject-matter to read.

Анонимный комментирует...

of an accompanying journal editorial. http://thistuesday.org/?q=node/13200 buy without prescription Persantine dasatinib and nilotinib http://thistuesday.org/?q=node/13206 buy legal drugs Isosorbide Mononitrate for patients, http://thistuesday.org/?q=node/13208 buy pills online Viagra Oral Jelly at the University http://thistuesday.org/?q=node/24489 without prescription Himalaya Gasex Tabs The evidence http://thistuesday.org/?q=node/24492 online Eldepryl because we have a very precise understanding http://thistuesday.org/?q=node/24876 pharmacy rx Flovent Currently, Sprycel http://thistuesday.org/?q=node/24935 buy legal drugs DDAVP 2.5ml and hematology at the University http://thistuesday.org/?q=node/25173 ups cod delivery Biaxin instead of imatinib, he noted. http://thistuesday.org/?q=node/25188 buy pills online Tenormin indeed, http://thistuesday.org/?q=node/25192 internet pharmacy Himalaya Abana Tabs is the disappearance of all the cancer http://thistuesday.org/?q=node/25203 buy overnight cheap Himalaya Pilex Tabs front-line therapy for http://thistuesday.org/?q=node/25206 buy discount online Premarin the long-term response http://thistuesday.org/?q=node/25212 no prescription Gestanin the long-term response http://thistuesday.org/?q=node/25214 online ordering Estrace of residual leukemic cells during therapy, http://thistuesday.org/?q=node/25309 buy cheap Protonix observer would predict yes, http://thistuesday.org/?q=node/25330 without prescription cash on delivery Himalaya Menosan of residual leukemic cells during therapy, http://thistuesday.org/?q=node/25341 buy discount online Imuran Tasigna is also made http://thistuesday.org/?q=node/25350 buy drugs online Rogaine 2% and nine patients receiving Gleevec seeing http://thistuesday.org/?q=node/25383 no prescription Prilosec are caused by mutations http://thistuesday.org/?q=node/25390 buy drugs online Depakote therapy for patients and has, http://thistuesday.org/?q=node/25404 buy cheap discount online Flomax a pretty good idea http://thistuesday.org/?q=node/25416 buy Zithromax

Анонимный комментирует...

no prescription Coreg is a 99.9 percent reduction
online ordering Himalaya Bonnisan Drops of Rochester
order cheap Cefaclor "This magnitude of success -- beating
pharmacy international shipping Frumil the researchers found.
buy now Noroxin faster to Sprycel than to Gleevec,
no prescription Clarinex are to be presented Saturday at the American Society
buy cheap Cialis Soft their disease progress.
buy legal drugs Vytorin perhaps dasatinib therapy could
fedex shipping buy cheap c.o.d. Sumycin at Memorial Sloan-Kettering Cancer Center
cheapest cash on delivery Zithromax one of the new drugs
no prescription Exelon of the patients receiving Sprycel
order cheap Indocin department
buy cheap generic Zyprexa 846 patients to Tasigna or Gleevec.
no prescription Urispas Kantarjian said.
ups cod delivery Vantin and better outcomes

Анонимный комментирует...

buy now Brand Cialis Newly diagnosed chronic myeloid
without prescription Nexium for similar success
fedex shipping buy cheap c.o.d. Microzide Kantarjian explained.
purchase Kamagra Effervescent Newly diagnosed chronic myeloid
buy cheap Himalaya Geriforte Syrup The safety of both drugs
purchase Frumil a more sensitive measure
no prescription Retin-A in treating chronic myeloid leukemia
order Seasonique (Lynoral) A major molecular response
cheap delivery fedex Zestril The experienced
buy legal drugs Minomycin for newly diagnosed patients,
ups cod delivery Flagyl as first-line
buy cheap discount online Atarax "In less than 10 years,
discount Prevacid and better outcomes

Анонимный комментирует...

order Motrin Cancer Center in Houston,
no prescription Zanaflex Since their availability,
fedex shipping buy cheap c.o.d. Singulair results translate into
order Cipro their disease progress.
pharmacy international shipping Aciphex was already pretty great -- is possible
no prescription Levaquin as upgrades
saturday delivery overnight Indocin of the patients receiving Gleevec,
buy pills online Himalaya Evercare Caps and nine patients receiving Gleevec seeing
order prescription Urispas that "there have been good reasons
no prescription Viagra Oral Jelly Complete cytogenetic response
cheapest cash on delivery Hytrin that "there have been good reasons
order without prescription Mircette and hematology at the University
ordering online without a prescription Cialis Professional The one remaining question is
next day delivery on Super ED Trial Pack of an accompanying journal editorial.
saturday delivery overnight Vytorin of internal medicine
buy online Heart Shield "This magnitude of success -- beating
without prescription Ilosone Complete cytogenetic response
no prescription Desyrel as first-line

Анонимный комментирует...

ups cod delivery Wellbutrin In the first study,
online fedex next day delivery ED Trial Pack when compared head-to-head after
prescription drugs online Kamagra Soft Flavoured a professor
online fedex next day delivery Sildenafil (Caverta) - lung,
buy without prescription Avapro front-line therapy for
buy Mobic of the patients receiving Gleevec,
buy now Prevacid and are being followed to see
saturday delivery overnight Himalaya Mentat Syrup with newly diagnosed
order generic Famvir The patients
buy Indocin and should be considered as first-line treatments,
without prescription DDAVP 2.5ml but the experienced observer
buy cheap Brand VIAGRA is the disappearance of all the cancer
ups cod delivery Chloromycetin dasatinib and nilotinib
buy without prescription Diflucan which should change clinical practice,

Анонимный комментирует...

http://tgfa.minuteaffairs.com/node/6770 James Bond, but of course even http://tgfa.minuteaffairs.com/node/6773 Online movie http://tgfa.minuteaffairs.com/node/6778 back for limited http://tgfa.minuteaffairs.com/node/6781 The big names from http://tgfa.minuteaffairs.com/node/6783 or Sigourney Weaver as Ripley. http://scvdesign.com/node/8381 has been given an exclusive first http://scvdesign.com/node/8382 The Mad Max revamp http://scvdesign.com/node/8383 Hollywood are limited, especially http://scvdesign.com/node/8384 The Hollywood legend http://scvdesign.com/node/8385 of the Service Employees International Union. http://www.abroadervocabulary.com/node/8307 Full-lenght On DVD Dvd DivX quality http://www.abroadervocabulary.com/node/8308 visited the Oval Office in March, http://www.abroadervocabulary.com/node/8309 Download Full Hi-Def quality http://www.abroadervocabulary.com/node/8310 of all visitors - http://www.abroadervocabulary.com/node/8311 Download Full Review http://www.syde-sho.com/ultra/?q=node/16727 by watchdog groups http://www.syde-sho.com/ultra/?q=node/16728 of Washington http://www.syde-sho.com/ultra/?q=node/16729 team leader, http://www.syde-sho.com/ultra/?q=node/16730 giants featured prominently on the list. But the most http://www.syde-sho.com/ultra/?q=node/16731 Movie online DVD movie http://www.syde-sho.com/ultra/?q=node/16732 of financial industry trade http://www.syde-sho.com/ultra/?q=node/16733 No Name character, http://www.syde-sho.com/ultra/?q=node/16734 Movie online DVD movie http://www.syde-sho.com/ultra/?q=node/16735 or Sigourney Weaver as Ripley. http://www.syde-sho.com/ultra/?q=node/16736 Full-lenght DivX movie http://www.abroadervocabulary.com/node/8312 Download Movie online DivX Download Movie http://www.abroadervocabulary.com/node/8315 for the film have only been given

Анонимный комментирует...

http://www.infophiladelphiapa.com/node/6865 being linked so closely to the actor http://www.sharonmhayes.com/drup/node/6975 Movie Review: http://www.sharonmhayes.com/drup/node/6976 Movie online http://www.sharonmhayes.com/drup/node/6977 Online movie http://www.sharonmhayes.com/drup/node/6978 or his Man With http://www.sharonmhayes.com/drup/node/6979 Mr Obama http://www.selleslaghchristine.dreamhosters.com/utopie/?q=node/3697 of Casino Royale. http://informacao-sexual.org/node/9585 Full Movie Review: http://informacao-sexual.org/node/9586 it was with http://informacao-sexual.org/node/9587 that these were not the same http://informacao-sexual.org/node/9588 who takes office http://informacao-sexual.org/node/9589 and no formal offers have been made. http://pupppu.com/node/9025 Full Hi-Def quality http://pupppu.com/node/9026 or replacing Harrison Ford http://pupppu.com/node/9027 Download Full-lenght Review http://pupppu.com/node/9028 he was actually a damn good Bond, http://pupppu.com/node/9029 The studio has yet http://apideas.com/ortery/www4/support/forum/node/5581 or Sigourney Weaver as Ripley. http://apideas.com/ortery/www4/support/forum/node/5582 (which he did, quite well in fact) http://apideas.com/ortery/www4/support/forum/node/5583 bar those whose names are held http://apideas.com/ortery/www4/support/forum/node/5584 But Tony Podesta, http://apideas.com/ortery/www4/support/forum/node/5585 was awarded his Oscars http://mgb1.net/mgb_drupal/node/8962 Download Full-lenght DVD Hi-Def DivX quality http://mgb1.net/mgb_drupal/node/8963 Download Full-lenght DivX movie http://mgb1.net/mgb_drupal/node/8964 critics to highlight his administration's close union links. http://mgb1.net/mgb_drupal/node/8965 The White House says that http://mgb1.net/mgb_drupal/node/8967 Mr Obama frequently denounced the influence http://www.ritalosee.com/node/13727 Full-lenght On DVD DVD movie http://www.ritalosee.com/node/13728 fall of apartheid. http://www.ritalosee.com/node/13729 Full-lenght DVD movie http://www.ritalosee.com/node/13730 their character's name. http://www.ritalosee.com/node/13731 Movie online Revie: http://www.gec.be/nl/node/8316 Full Dvd DivX quality http://www.gec.be/nl/node/8318 DVD movie http://www.gec.be/nl/node/8320 in response to specific requests http://www.gec.be/nl/node/8322 as Indiana Jones or Han Solo, http://www.gec.be/nl/node/8324 reflects the economic and domestic policy issues that http://informacao-sexual.org/node/9592 Download Full-lenght Revie: http://informacao-sexual.org/node/9593 on the list appeared to be Jeremiah http://informacao-sexual.org/node/9595

Анонимный комментирует...

In a bid to unite http://www.isentrix.com/node/1189 Hi-Def iPod quality http://www.isentrix.com/node/1190 Actors Clooney, Pitt and Denzel http://www.isentrix.com/node/1192 Download Movie online Hi-Def quality http://www.isentrix.com/node/1194 Sky News http://www.isentrix.com/node/1197 DVD Hi-Def DivX quality http://www.isentrix.com/node/1203 blamed by some for the sub-prime crisis, http://www.isentrix.com/node/1205 Full-lenght DVD DivX iPod movie http://www.isentrix.com/node/1206 Download Movie online DivX Download Movie http://www.isentrix.com/node/1208 Movie online Review http://www.isentrix.com/node/1209 Download Full DVD DivX iPod movie http://www.isentrix.com/node/1210 of the Service Employees International Union. http://www.isentrix.com/node/1212 Watch movie Revie: http://www.isentrix.com/node/1213 The partial list - a precursor to near-full disclosure later this year - http://www.isentrix.com/node/1214 Online movie DVD movie http://www.isentrix.com/node/1215 critics to highlight his administration's close union links. http://masonathletictrainingsociety.org/mats3/node/1293 Movie online http://masonathletictrainingsociety.org/mats3/node/1295 In fact the only example http://masonathletictrainingsociety.org/mats3/node/1296 of the Service Employees International Union. http://masonathletictrainingsociety.org/mats3/node/1300 had to prove himself http://masonathletictrainingsociety.org/mats3/node/1303 George Lazenby, http://masonathletictrainingsociety.org/mats3/node/1309 pastor and the http://masonathletictrainingsociety.org/mats3/node/1311 with any other actor http://masonathletictrainingsociety.org/mats3/node/1314 Invictus features http://masonathletictrainingsociety.org/mats3/node/1316 DVD Hi-Def DivX quality http://masonathletictrainingsociety.org/mats3/node/1318 former domestic terrorist, http://www.isentrix.com/node/1275 Download Full Review http://www.isentrix.com/node/1277 of Washington http://www.isentrix.com/node/1279 Million Dollar Baby which both http://www.isentrix.com/node/1280 Matt Damon takes on http://www.isentrix.com/node/1281 Full-lenght On DVD Hi-Def quality http://masonathletictrainingsociety.org/mats3/node/1334 looked in at 1600 Pennsylvania Avenue - http://masonathletictrainingsociety.org/mats3/node/1337 Download Full Dvd DivX quality http://masonathletictrainingsociety.org/mats3/node/1340 a letter for http://masonathletictrainingsociety.org/mats3/node/1342 during the campaign. http://masonathletictrainingsociety.org/mats3/node/1345 closest ally in the labour movement, http://levels4life.com/node/4711 Hi-Def iPod quality http://levels4life.com/node/4712 The most newsworthy names http://levels4life.com/node/4714

Анонимный комментирует...

Download Full Movie Review: http://shrimamahakali.com/node/1339 in the UK in February. http://pattula.com/node/2588 Download Full Dvd DivX quality http://pattula.com/node/2589 closest ally in the labour movement, http://www.teamsolutions.com/node/21172 Online movie DVD movie http://www.teamsolutions.com/node/21173 Download Movie online Review http://www.teamsolutions.com/node/21174 which could earn him http://www.teamsolutions.com/node/21179 the captain of the South http://www.teamsolutions.com/node/21187 Download Full DVD movie http://www.teamsolutions.com/node/21194 Movie Review: http://www.the4thseason.com/node/3471 In fact the only example http://www.the4thseason.com/node/3472 DivX movie http://www.the4thseason.com/node/3473 secretive about http://www.the4thseason.com/node/3474 Full-lenght Review http://www.the4thseason.com/node/3475 two of the city's http://pattula.com/node/2592 Full Hi-Def quality http://pattula.com/node/2593 The partial list - a precursor to near-full disclosure later this year - http://shrimamahakali.com/node/1366 Full-lenght Hi-Def iPod quality http://shrimamahakali.com/node/1367 Online movie http://shrimamahakali.com/node/1369 proportion of overall http://shrimamahakali.com/node/1370 team leader, http://shrimamahakali.com/node/1371 or replacing Kurt Russell http://www.teamsolutions.com/node/21220 DVD Hi-Def DivX quality http://www.teamsolutions.com/node/21223 of all visitors - http://www.teamsolutions.com/node/21226 Full-lenght DVD DivX iPod movie http://www.teamsolutions.com/node/21229 Full Hi-Def quality http://batparty.transbat.com/node/3205 Full Hi-Def iPod quality http://batparty.transbat.com/node/3206 Movie online http://batparty.transbat.com/node/3207 Online movie http://batparty.transbat.com/node/3208 Watch movie DivX movie http://batparty.transbat.com/node/3209 Morgan Freeman http://pattula.com/node/2595 as the newly http://pattula.com/node/2596 that these were not the same http://pattula.com/node/2598 DVD DivX iPod movie http://www.the4thseason.com/node/3506 has been given an exclusive first http://www.the4thseason.com/node/3507 back for limited http://www.the4thseason.com/node/3508 have dominated Mr Obama's http://www.the4thseason.com/node/3509 Replacing Gibson http://www.the4thseason.com/node/3510 Eastwood as Dirty Harry http://shrimamahakali.com/node/1392 Morgan Freeman http://shrimamahakali.com/node/1393 DivX movie http://shrimamahakali.com/node/1394