370 likes | 492 Views
LakeLoadResponse(LLR)-työkalun esittely 11.2.2008 Niina Kotamäki, Anita Pätynen, Olli Malve, Timo Huttula, Kai Rasmus. http://coherens/~tht_kr/lakestate/. LLR:n rakenne . 5 erilaista mahdollisuutta ajaa malleja Tavoitekuormien laskeminen fosforin Hyvän/tyydyttävän luokkarajan avulla
E N D
LakeLoadResponse(LLR)-työkalun esittely11.2.2008Niina Kotamäki, Anita Pätynen, Olli Malve, Timo Huttula, Kai Rasmus http://coherens/~tht_kr/lakestate/
LLR:n rakenne • 5 erilaista mahdollisuutta ajaa malleja • Tavoitekuormien laskeminen fosforin Hyvän/tyydyttävän luokkarajan avulla • Tavoitekuormien laskeminen typen Hyvän/tyydyttävän luokkarajan avulla • Tavoitekuormien laskeminen sekä typen että fosforin H/T luokkarajojen perusteella • Tavoitekuormien laskeminen eo. lisäksi myös klorofylli a:n H/T luokkarajan avulla • (Tavoitekuormien laskeminen biomassan H/T luokkarajan avulla )
Tavoitekuormien laskeminen • Hyvän ja tyydyttävän (H/T) vedenlaadun luokkarajalle on määritetty järvityyppikohtaiset pitoisuusarvot*) • Kokonaisfosforille (TotP, µg/l) • Kokonaistypelle (TotN, µg/l) • Klorofylli a:lle (Chla, µg/l) • Kasviplanktonbiomassalle (mg/l) *) Suomen ympäristökeskus, Riista- ja kalatalouden tutkimuslaitos. Pintavesien ekologisen luokittelun vertailuolot ja luokan määrittäminen. 22.1.2008
Tavoitekuormien laskeminen • LLR:llä voidaan laskea kokonaisfosforin ja -typen tavoitekuormat eri muuttujien H/T –luokkarajojen mukaan • Esim. voidaan laskea millaisella järveen tulevalla kokonaisfosforikuormalla saavutetaan järven hyvää tilaa vastaava järviveden kokonaisfosforipitoisuus
Esimerkkitapaus: Hiidenvesi • Tutustutaan LLR:n käyttöön esimerkkikohteen avulla: • Hiidenvesi on Länsi-Uudellamaalla, Karjaanjoen valuma-alueella sijaitseva keskisuuri järvi • Muodoltaan hyvin monimutkainen • Tarkastellaan vain syvänneosia
Hiidenvesi Olkkalanjoki Karjaanjoki (Vanjoki) Väänteenjoki
Tavoitekuormien laskeminen • Tarkastellaan esimerkkinä vaihtoehtoa 4: Laske tavoitekuorma Chl-a H/T luokkarajan perusteella • Vaihtoehdossa yhdistyvät kaikki edeltävät vaihtoehdot • Vaihtoehto 5 on vastaava, mutta perustuu kasviplanktonbiomassan H/T -luokkarajaan • Klikkaamalla painiketta siirrytään seuraavalle sivulle antamaan laskemiseen tarvittavia järvitietoja
Järven tiedot • Jos järvi löytyy järjestelmästä saadaan syöttötietokenttiin automaattisesti järven • Tilavuus • Keskisyvyys • Järvityyppi • Chla- ja ravinnehavainnot • väriluku • Jos järveä ei löydy, kaikki tiedot on syötettävä itse. • Chla- ja ravinnetietojen syöttö käydään läpi jäljempänä
Järven etsiminen järjestelmästä • Klikataan Etsi järvi järvitaulukosta • Avautuvaan kenttään kirjoitetaan järven nimi tai nimen osa (esim. ”Pyhäjärvi” tai ”Pyhä”) ja klikataan Submit • Hakusanalla ”järvi” saa näkyviin kaikki järjestelmässä olevat järvet • Näyttöön ilmestyvät taulukoituna hakua vastaavat järvet (tai ilmoitus ettei hakusanoja vastaavia järviä löydy) • Jos oikea järvi löytyy taulukosta, valitaan se klikkaamalla Valitse
Ennustetodennäköisyys (oletus) • Käyttäjä ei voi muuttaa, aina 50 • Ennustetodennäköisyyden ollessa 50, mallien antamilla tavoitekuormilla voidaan saavuttaa esim. järviveden H/T –luokkarajaa vastaava klorofylli a-pitoisuus, mutta 50 % todennäköisyydellä lopullinen pitoisuus voi olla tavoiteltua suurempikin
Ennustetodennäköisyys2 ja Ennustetodennäköisyys3 • Käyttäjä voi valita • Esim. antamalla arvon 80, klorofylli a –pitoisuuden ylittymisen todennäköisyys mallin antamilla tavoitekuormilla on enää 20% • Voit vertailla miten tavoitekuormat muuttuvat, kun ennustetta tarkennetaan
Havaintotaulukko • Seuraavaksi tarvitaan havaintotaulukko • Taulukkoon kootaan viipymäajalle laskettuina keskiarvoina: • Kokonaisfosforikuorma (kg/d) • Kokonaistyppikuorma (kg/d) • Kokonaisfosforipitoisuus (µg/l) • Kokonaistyppipitoisuus (µg/l) • pitoisuudet järven (keskeisimmästä) syvänteestä saatujen havaintojen (mielellään syvyyspainotettuina) keskiarvoina • Luusuan virtaama (m3/s)
Havaintotaulukko • Havainnot etsitään mahdollisimman pitkältä ajanjaksolta • Jos havaintoja ei ole, klikataan Minulla ei ole havaintoja :( • Tavoitekuormat voidaan laskea vaihtoehtoisesti kuormitusarvon ja virtaaman avulla (tästä lisää myöhemmin)
Havaintotaulukon kokoaminen • Järven viipymä (vuorokausissa) lasketaan kaavalla missä V = tilavuus Q = virtaama • Viipymä pyöristetään ylöspäin vuosiksi
Havaintotaulukon kokoaminen • Jos viipymä on esim. yksi vuosi, otetaan kaikki vuoden aikana saadut näytteenottotulokset ja lasketaan niiden keskiarvo = viipymäajan keskiarvo • Jos näytteenottoväli on pidempi kuin järven viipymä, voidaan taulukkoon merkitä suoraan näytteenottojen tulokset • tulosten edustavuutta arvioitava tarkoin !
Havaintotaulukon kokoaminen • Viipymäaikojen keskiarvot kootaan Excel-taulukkoon, joka tallennetaan csv-muodossa • Sarakkeet nimetään (huom. isot ja pienet kirjaimet!) : LN = kokonaistyppikuorma (kg/d) LP = kokonaisfosforikuorma (kg/d) TotN = kokonaistyppipitoisuus (µg/l) TotP = kokonaisfosforipitoisuus (µg/l) Q = luusuan virtaama (m3/s)
Havaintotaulukko • Kun taulukko on valmis ja tallennettu • Palataan LLR tietojensyöttösivulle • Klikataan havaintotaulukkorivin painiketta Browse • etsitään taulukkotiedosto hakemistosta ja klikataan Open
Chla- ja ravinnetiedot • Pitää taulukoida ja hakea, jos järvi ei löydy järjestelmästä • Taulukkoon kootaan eri vuosilta kaikki (keskeisimmän) syvänteen heinä-elokuun näytteenottokertojen määritystiedot kokonaisfosforista, kokonaistypestä ja klorofylli a:sta. (Ei keskiarvoja!)
Chla ja ravinnetiedot • Tiedot kootaan Exel-taulukkoon, joka tallennetaan csv-muodossa • Sarakkeet nimetään (huom. pienet kirjaimet!) : totp = kokonaisfosfori (µg/l) totn = kokonaistyppi (µg/l) chla = klorofylli a (µg/l)
Chla ja ravinnetiedot • Kun taulukko on valmis ja tallennettu • Palataan LLR tietojensyöttösivulle • Klikataan Chla ja –ravinnetiedot -rivin painiketta Browse • Haetaan taulukkotiedosto hakemistosta ja klikataan Open
Tavoitekuormien laskeminen • Kun kaikki tiedot ovat valmiina taulukossa, klikataan Laske tavoitekuormat • LLR tekee laskelmat ja siirtyy tulossivulle
Tulosten tulkinta • Ensimmäisenä on laskentatietotaulukko ja ilmoitus siitä, onko tiedot syötetty oikein. • Sitten kuvaajat • Lopussa taulukko
Kuvaajasta näet miten hyvin mallin antamat arvot vastaavat järvessä havaittuja pitoisuuksia • havaittu pitoisuus │ 90 %:n luottamusraja +―+ mallilla laskettu ennuste Kuvaajista näet miten järven kokonaisfosforipitoisuus eri mallien mukaan muuttuu, kun järveen tuleva kokonaisfosfo-rin pintakuorma muuttuu. Vaakaviivat näyttävät järven eri tilojen luokkarajat kokonaisfosforipi-toisuuksille. Tärkein on Hyvä/Tyydyttävä –raja, jonka alapuolelle jäävät pitoisuudet ovat vesipuitedirektiivin vaatimusten mukaisia. Ohjelma laskee tavoitekuormalle kolmen eri mallin antamat ennusteet. Mallit perustuvat erilaisiin järviaineistoihin. Lisätietoja malleista saat Kuvaajista näet miten Chapran järvikohtaiseen aineistoon perustuvan mallin ennustama tavoitekuorma muuttuu eri ennustetodennäköisyyk-sillä
Kuvaajasta näet miten hyvin mallin antamat arvot vastaavat järvessä havaittuja pitoisuuksia • havaittu pitoisuus │ 90 %:n luottamusraja +―+ mallilla laskettu ennuste Kuvaajista näet miten järven kokonaistyppipitoisuus eri mallien mukaan muuttuu, kun järveen tuleva kokonaistypen pintakuorma muuttuu. Vaakaviivat näyttävät järven eri tilojen luokkarajat kokonaistyppipitoi-suuksille. Tärkein on Hyvä/Tyydyttävä –raja, jonka alapuolelle jäävät pitoisuudet ovat vesipuitedirektiivin vaatimusten mukaisia. Ohjelma laskee tavoitekuormalle kolmen eri mallin antamat ennusteet. Mallit perustuvat erilaisiin järviaineistoihin. Lisätietoja malleista saat Kuvaajista näet miten Chapran järvikohtaiseen aineistoon perustuvan mallin ennustama tavoitekuorma muuttuu eri ennustetodennäköisyyk-sillä
Kuvaajista näet, millaisilla kokonaisfosforin ja kokonaistypen pintakuormilla järviveden klorofylli a -pitoisuus saadaan alle Hyvä/Tyydyttävä -luokkarajan Klorofylli a-ennuste on laskettu Chapran mallia, jonka tarkemman kuvauksen löydät Prosenttiluvut tarkoittavat eri ennustetodennäköisyyksiä. Voidaan siis tarkastella, miten erilaiset ennustetodennäköisyydet vaikuttavat tavoitekuormiin Nuolet näyttävät mikä on mallin mukaan järven klorofylli a -pitoisuus kuormituksen nykytasolla. (Mallin antaman pitoisuuden ylittymismahdollisuus on 50 %.)
Kuvaajista näet, millaisilla kokonaisfosforin ja kokonaistypen pintakuormilla järviveden klorofylli a -pitoisuus saadaan alle Hyvä/Tyydyttävä -luokkarajan Klorofylli a-ennuste on laskettu Dillon&Riglerin mallilla, jonka tarkemman kuvauksen löydät Prosenttiluvut tarkoittavat eri ennustetodennäköisyyksiä. Voidaan siis tarkastella, miten erilaiset ennustetodennäköisyydet vaikuttavat tavoitekuormiin Nuolet näyttävät mikä on mallin mukaan järven klorofylli a -pitoisuus kuormituksen nykytasolla. (Mallin antaman pitoisuuden ylittymismahdollisuus on 50 %.)
Tavoitekuorman laskeminen ilman järvihavaintoja • Kuormitusarvon saa tarvittaessa haettua VEPSistä • Tavoitekuormien laskeminen tulisi kuitenkin aina jos vain mahdollista, tehdä järvihavaintoihin perustuen!
Tavoitekuorman laskeminen ilman järvihavaintoja • Jos havaintoja ei ole, klikataan Minulla ei ole havaintoja :( • Tietojensyöttötaulukko muuttuu, ja havaintotaulukkokohdan tilalle saadaan syöttötietokentät kuormitusarvolle (kg/d) ja virtaamalle (m3/s)
Kuvaajasta näet miten järven kokonaisfosforipitoisuus eri mallien mukaan muuttuu, kun järveen tuleva kokonaisfosforin pintakuorma muuttuu Vaakaviivat näyttävät järven eri tilojen luokkarajat kokonaisfosforipi-toisuuksille. Tärkein on Hyvä/Tyydyttävä –raja, jonka alapuolelle jäävät pitoisuudet ovat vesipuitedirektiivin vaatimusten mukaisia. Ohjelma laskee tavoitekuormalle kahden eri mallin antamat ennusteet. Mallit perustuvat erilaisiin järviaineistoihin. Lisätietoja malleista saat
Kaaviokuva LLR- ympäristöstä LLR/internet • syöttötiedot Internet-työkalun näkymä ja toiminnallisuus Linux • tiedot vastaan • tallennus PERL, HTML, Java Kuvat, tiedot taulukoihin R-Linuxissa • tilastolliset analyysit • kuvien muodotaminen R/S-Plus Parametrien estimaatit pääasiassa täältä, mallien ajot LinBUGS • Bayes-estimointi • MCMC:llä OpenBUGS /LinBUGS
LLR:ään tulossa lähiaikoina (testauksessa) • Jos käyttäjällä ei ole havaintoja järvestä, voidaan typen ja fosforin kokonaiskuorma sekä virtaama arvioida • Biomassaraja-arvon avulla lasketaan tavoite • Internetiin • Herttaan ja VEPSiin kytkeminen • Järjestelmässä olevien järvien tarkemmat sijaintitiedot sekä klorofylli- ja ravinnehavainnot • Uusien havaintojen syöttö chla-malliin