Aattelin että vois olla asiaa aloittaa uusi ketju tälle...
Tuossa toisessa ketjussa kirjoittelin että metsästin pikkasen tuota NAO-dataa eli North Atlantic Oscillation, siitä ei sen enmpää, vielä.
Samalla kävi ilmi että se liikkuu myös nimellä "Teleconnection incidies" sekä että myös "500 mb Geopotential height anomaly (m)" (GPA) on yksi arvo joka kuuluu joukoon. NOAA-servulta jota GFS:ää haetaan löytyy tuo GPA ja se haetaan myös tästä lähtien (arvot gpa500mb ja gpa1000mb) :)
Tuo GPA ja NAO näyttää kävelevän aika käsi kädessä kun tätä siva katselin: http://www.nc-climate.ncsu.edu/climate/patterns/NAO.html
Tein testikartan kyseisellä GPA:lla joka liitteenä.
Seuaraava koeajo,
NOAA:lta napattua kuukausittaista NAO-dataa, liitteenä käppyrä 01/1990-12/2010. Jotain on tapahtunut vuoden 2009 lopulla joka jatkuu edelleen koska samanalaista monttua ei näy koko käppyrässä vaikka ottais 1950 asti :o
Toisaalta, siinä näkyy myös 2008/2009 lämmin talvi todella hyvin.
Seuraava testikäppyrä ;D
Kuuluisasta UAH-datasta tehty käppyrä pohjoisen pallonpuoliskon lämpötilapoikkeamasta. Huomioitavaa on että heillä jo on se FMI:kin tulossa oleva keskiarvo 1981-2010, otettu käyttöön viimeisen päivityksen yhteydessä. Heidän reamden mukaan ero 1971-2000 ja 1981-2010 väillä on yleisesti ottaen noin 0.1°.
Dataa löytyy 1979-.
Tuota liitteen lätkämailaa ei edes FMI saa käännettyä ylöspäin, sen verran rajusti ollaan alas tultu ;D
UAH = University of Alabama/National Space Science and Technology Center.
Ja viimeisenä, muttei pienimpänä käpyyränä,
HADCRUT3-käyrä, eli se suruisan kuuluisa mitä IPCC:kin käyttää ilmstahössötykseensä. Jep, kyseinen data löytyy kun löytyykin netistä UK MetOfficen/Hadley Centren toimesta.
Mutta löysin varsin mielenkiintoista tietoa siitä samaiselta sivulta josta datakin (http://www.cru.uea.ac.uk/cru/data/temperature/);
Data perustuu vain 3000 sääaseman tietoihin joka mielstäni on todella vähän. Ja mitä parasta/pahinta, näitä asemia ei edes ole joka puolella maailmaa vaan eniten länsimaissa kuten USA, Etelä-Canada, Eurooppa ja Japani. Sekä verkko on kasvanut 3000:een asemaan vasta vuosina 1951-90. Ja tätä maailman päättäjät/FMI & co uksovat kun raamattua, jep, jep...
Dataa kyllä löytyy vuodesta 1850 asti (mutta javaskriptin timestamppi ei osanut kun vuodesta 1902) ja keskiarvona käytetty 1961-90.
Ero UAH ja HADCRUT:in välillä näyttää olevan se että UAH käyttää NOAA:n sateliittidataa ja siten kattaa koko maapallon kun HADCRUT on noiden asemien mittaamaa.
Tämä nyt meni pikkasen ohi aiheen kun oli kyse NAO:sta mutta ei se mitään ;)
Löysin sen kolmannen data-lähteen jota ylseisesti käytetään maapallon lämpötiloissa nimittäin RSS. Ongelma näissä on se että kaikki neljä käyttää eri aikajaksoa keskiarvona (yllätys,yllätys) joten eivät edes ole suoraan vertailukelpoisia sekä se että UAH ja RSS-dataa ei ole saatailla kun vuodesta 1979. Mutta, tässä ne tulevat, tein kaikista neljästä käppyrät 1980-2010.
UAH
NOAA-sateliiteilla mitattu, k-arvo = 1981-2010
(http://www.nordicweather.net/cache/uah1980.png)
Data: University of Alabama/National Space Science and Technology Center
RSS
NOAA-sateliiteilla mitattu, k-arvo = 1979-1998
(http://www.nordicweather.net/cache/rss1980.png)
data: Remote Sensing Systems/NASA
HadCRUT3
Sääasemilla mitattu, asemia noin 3000, lähinnä länsimaissa, k-arvo = 1961-90
Tämä on se mitä IPCC & co käyttävät
(http://www.nordicweather.net/cache/hadcrut1980.png)
data: UK MetOffice
NCDC
K-arvo = 1901-2000
Maalueiden data: Global Historical Climate Network (GHCN)
Merialuiden data: Extended reconstructed sea surface temperature (ERSST).
(http://www.nordicweather.net/cache/ncdc1980.png)
data: NOAA/NCDC
Tuota UK MetOfficen HadCRUT3-dataa saa käyttää ilman sen kummempia lupia ""For personal study and research" ;D
EDIT: Löysin NOAA:lta datasetin jossa käytetty keskiarvona vuosia 1901-2000
Tein löydön liittyen käytettyyn keskiarvoon. :)
NOAA:n NCDC:n sivuilta löytyi dataa jossa käyetty keskiarvona 1901-2000 eikä jotain 20 tai 30 vuoden jaksoa.
Dataa näyttää löytyvän vuodesta 1880.
Ei kun pureskelemaan sitäkin käppyräksi ;D
EDIT Pureskeltu ja käppyrä lisätty edellisen viestiin muitten kaveriksi ;D
Dodiih, sivuille tehty 2 uutta sivua:
NAO-sivu, jhon on vielä tulossa lisää tavaraa:
http://www.nordicweather.net/index.php?sivu=70&lang=fi
Pohjoisen maapallopuoliskon lämpötilapoikkeamat, jossa kaikki 4 yllä olevaa käppyrää:
http://www.nordicweather.net/index.php?sivu=71&lang=fi
;D
Niin, ja kyseistä skriptiä?
Kunhan saan kasattua pikkasen sitä niin käytetty php-pätkä tulee jakoon ;)
Erinomaista työtä jälleen kerran! Mielenkiintoisia käppyröitä kaikki ja tulen ehdottomasti jossain vaiheessa ottamaan nämä käyttöön omallakin sivulla, mutta nyt vähän huilin noiden tilastopirulaisten jälkeen :D. Mielenkiintoista nähdä mihin suuntaan kukin käppyrä vuosien varrella kehittyy.
Lainaus käyttäjältä: weatherc - sunnuntai, 09.01.2011, 01:58
Mutta löysin varsin mielenkiintoista tietoa siitä samaiselta sivulta josta datakin (http://www.cru.uea.ac.uk/cru/data/temperature/);
Data perustuu vain 3000 sääaseman tietoihin joka mielstäni on todella vähän. Ja mitä parasta/pahinta, näitä asemia ei edes ole joka puolella maailmaa vaan eniten länsimaissa kuten USA, Etelä-Canada, Eurooppa ja Japani. Sekä verkko on kasvanut 3000:een asemaan vasta vuosina 1951-90. Ja tätä maailman päättäjät/FMI & co uksovat kun raamattua, jep, jep...
Juuri tuo on mielestäni se ongelma näissä mittauksissa. Pohjoisella pallonpuoliskolla on huomattavasti enemmän noita havaintoasemia, sillä eteläinen pallonpuolisko on pääosin meren vallassa. Jos nämä pohjoisen pallonpuoliskon asemat vielä sijaitsevat pääosin länsimaissa, ei täysin tarkkoja tilastoja voida tietenkään saada. Suuntaa antavia kylläkin, mutta virhelukema koko maapallon osalta saattaa pahimmillaan olla merkittävä.
Noniin tässä tulee skriptit :)
Ohje:- Seuraavat js:ät täytyy olla ladattuna sivulla jossa haluat näyttää käyriä:
<!--[if IE]><script type="text/javascript" src="http://www.nordicweather.net/jquery/excanvas.pack.js"></script><![endif]-->
<script type="text/javascript" src="http://ajax.googleapis.com/ajax/libs/jquery/1.4/jquery.min.js"></script>
<script type="text/javascript" src="http://www.nordicweather.net/jquery/jquery.flot.comb.js"></script>nordicweatherissä olevat skriptit voit ladata omalle serverillesi jos haluat samoin kuten skriptien lopussa oleva
<script type="text/javascript" src="http://www.nordicweather.net/jquery/jquery.uahgraph.own.js"></script>- Skripteissä ei ole css:ää definoituna joten jos se löytää sopivat tiedot omasta css:täsi se käyttää niitä
- Tiedostot ovat suorat kopiot omista joten pätkä joka pitää lisätä sivulle alkaa # SCRIPT START ja loppuu # SCRIPT END
- Molemmat hakevat/tarkistavat sekä NAO että lämpötiladatan tuoreuden ja hakee tuoreen jos tarvetta on, cacheaika on 1 vrk (86400 sekuntia)
- Ovat tehty siten että NAO-skripti näytetään eri sivulla kun lämpötila-skriptit eikä jquery.uahgraph.own.js-toimi ilman pientä muokkausta.
Kannattaa pitää mielessä että HadCRUT3-käyrän data on UK Metofficen jota saa käyttää "For personal study and research", muut ovat täysin "vapaata riistaa" kun ovat peräisin Jenkkilästä "julkiselta domainilta". ;)
LainaaJuuri tuo on mielestäni se ongelma näissä mittauksissa. Pohjoisella pallonpuoliskolla on huomattavasti enemmän noita havaintoasemia, sillä eteläinen pallonpuolisko on pääosin meren vallassa. Jos nämä pohjoisen pallonpuoliskon asemat vielä sijaitsevat pääosin länsimaissa, ei täysin tarkkoja tilastoja voida tietenkään saada. Suuntaa antavia kylläkin, mutta virhelukema koko maapallon osalta saattaa pahimmillaan olla merkittävä.
Osuit naulan kantaan tässä hommassa ;D
Sitten tähän väliin purkkakoodin mestarinäyte ;D 8)
Oppisin äsken että se on php:n strtotime/unix timestamp joka esti sen etten saanut parsittua vuodet ennen 1902 noihin käyriin.
Alkuperäinen data kun on muodossa pp-kk-vvvv eikä se toimi timestamppina vaan se pitää saada muotoon 1234567890 eli sekuntia sitten 1.1.1970, joka voi olla myös negatiivinen.
Kiertotie vuosille ennen 1902:
1. Lisättiin vuosilukuihin 100 jotta esim. vuodesta 1890 tuli 1990
2. Tehtiin timestamp strtotime:llä siitä
3. Poistettiin 3136320000 sekuntia siitä (= 100 vuotta)
(4. kerrattiin 1000:lla jotta siitä tuli javascript-timestamp [= millisekuntia] mutta se tehtiin kaikille)
Ja violá, meillä on timestampit myös vuosille 1850-1902 ;D ;D ;D
Oman sivun käyrissä näkyy nyt kaikki vuodet, liitteet eivät ole vielä päivitetty.
Seuraava askel: NAO-ennuste ;D
Löysin yhden sivun jossa NAO-ennustedataa on tehty GFS-ennusteesta. Poisulkemismetodilla, koska en löytänyt mistään infoa että mitä arvoa siihen käyetään, niistä arvoista jota GFS:ssä on saatavilla eli:
- mittayksikkö on 'm' eli metri
- kyseessä on keskiarvon poikkeama
- sivullani olevan Wikipediapätäkän mukaan siihen käytetään 500 mb arvoa
niin tulin siihen tulokseen että "Geopotential height anomaly" sopii enemmän kun hyvin siihen haarukkaan. Myös gpa500mb:n antamat lukemat ovat samoissa lukemissa (eli vaihtelee parisataa pakkaselta plussalle) ja jossain oli sivu jossa gpa:ta verrattiin NAO:on niin hep :)
Kaivoin esiin "GFS to text"-skriptini jossa GFS-data pureskellaan tekstimuotoon ja hain yhden ennusteen gpa500mb:stä ja nyt odotellaan että dedin OpenGraADS sais purekeltua koko pohjoisen maapallopuoliskon datat numeeriseen muotoon ;)
Enpä tehnytkään koko pohjoista maapallonpuoliskoa koska php:ssa loppus muisti kesken ;)
Tein suunnillen Pohjoismaiden kokoisen ruudun koordinaatit josta tuli "vaivalloiset" 184275 datapistettä, laskin ne yhteen ja jaoin sen määrällä jotta sain keskiarvon ja kappas kumma, lukemaksi tuli 57.2, joka täsmää pelottavan hyvin tuon löytämäni NAO-ennusteen saman ajan kanssa kun siinä oli "NAO East"-ennuste :)
Dodiih, erittäin experminetal NAO ennuste Pohjoismaille GFS-datasta ;D
- Ruudukkoa jota käytetty on 5°W - 40°E ja 45°N - 80°N
- Arvona 500 mb Geopotenitial height anomaly (GPA)
Miten se on laskettu:
Jokaisesta ennusteesta laskettu yhteen kaikki ruudukon GPA-arvot ja jaettu ne arvojen lukumäärällä (6462 kpl) jotta saatiin keskiarvo, yhteensä 78 lukemaa seuraavalle 14 päivälle.
(http://www.nordicweather.net/cache/naoennuste.png)
Miten päädyin että täytyy olla GPA? (Saa korjata jos olen väärässä ;))
Poissulkemismetodilla
- Jossain NAO:a koskekevassa jutussa oli myös GPA-kartta
- Wikipedian jutussa oli maininta 500 mb
- GFS-ennusteessa jota löysin oli lukeman yksikkönä m eli metriä, ja ainut data GFS-datasta joka sopii näihin on juuri GPA
- Uhkaavan samantapainen tulos yllä olevassa käyrässä jota oli löydetyn GFS-ennusteen NAO East-ennusteessa ottaen huomioon se että kyseisen ennuste kattaa koko itäisen osan pohjoisesta puoliskosta
EDITTIÄ
Kokeilin ajaa koko pohjhoista maapallon puoliskoa ja nyt toimis ok, oli vissiin jotain häikkää tuossa aikasemmin gs-filussa koska nyt tuli datat todella nopeasti, tuotos hivenen erilainen kun tuo "Pohjoismaat only" ;)
Nyt käytetyt koordinaatit: 90°W - 90°E ja 0°N - 88°N (viimeisistä 2 asteteesta Pohjoisnavalla ei kunnollista dataa).
Nyt on arvoja per ennuste hitusen vajaat 64000 ekalla viilolla (GFS 0.5°) ja 16000 toisella viikolla (GFS 1.0°) :P
(http://www.nordicweather.net/cache/naoennuste2.png)