Növényi produkció mérése mikrometeorológiai módszerekkel. Ökotoxikológus MSc, április 21

Növényi produkció mérése mikrometeorológiai módszerekkel

Ökotoxikológus MSc, 2015. április 21.

Felszín – légkör kölcsönhatások 

  

A legalapvetőbb kölcsönhatás a felszín és a légkör között: a szél, és annak súrlódása a felszínnel. vektor mennyiség u, v, w komponensek 3-as természete van:   



átlagos szélsebesség (advekció) periodikus hullámok (főleg éjszaka) véletlenszerű fluktuációk (turbulencia !!!  kicserélődési folyamatok)

A légkörben végbemenő folyamatok tanulmányozása komplex tudomány.  

rendezettek: portölcsérek, tornádók, hurrikánok kaotikusak: turbulens áramok

Felszín – légkör kölcsönhatások  

lamináris áramlás turbulens áramlás

Copyright(C)1980 Y. IRITANI, N. KASAGI and M. HIRATA, All rights reserved.

Felszín – légkör kölcsönhatások

Lamináris határréteg

Turbulens határréteg

Felszín – légkör kölcsönhatások

  



Súrlódás  mechanikai turbulencia Felhajtó erő  konvektív turbulencia A légkörben együtt vannak jelen és vesznek részt a felszín és a légkör közötti kicserélődési folyamatokban. 

A turbulencia:  

 





nem-lineáris 3D-s diffúz  átkeveredés disszipatív több mérettartományt felölel 3km-től 10-3 m-ig

Felszín – légkör kölcsönhatások

A hőmérsékleti rétegződés szerepe: magasság

labilis légrétegződés

stabilis légrétegződés

hőmérséklet

Teljesen kormányzott konvekció (neutrális eset) esetén az örvények kör alakúak. Átmérőjük (l) megegyezik a szabad úthosszal (kz). A horizontális és vertikális sebesség fluktuációk megegyeznek a súrlódási sebességgel. Instabil esetben: a vertikális sebesség fluktuációk nagyobb lesznek mint a horizontálisak, mivel az örvénynek horizontális kiterjedése nagyobb lesz, mint a szabad úthossz. Stabil esetben: az örvények horizontális irányban nyúlnak meg

Felszín – légkör kölcsönhatások

Logaritmikus szélprofil 



tapasztalat: a felszín közelében kisebb a szélsebesség. Könnyen belátható: u  a z

z

Felszín – légkör kölcsönhatások



Dimenzió analízis segítségével levezethető: a=[m s-1] u a k

k: von Kármán konstans (=0.4) u*: súrlódási sebesség, ~ momentum átvitel intenzitásával

u u   z kz



u*  z  u z   ln   k  z0 

z0: érdességi magasság: az a magasság, ahol u=0.



Magas vegetáció esetén a szélprofil megemelkedik

u*  z  d   u z   ln k  z0  z0= érdességi magasság d: 0-pont eltolódási szint d=0.6*h z0=0.1*h

Felszín – légkör kölcsönhatások

z0

d

Felszín – légkör kölcsönhatások

A növényzet szerepe: d, z0: változik DE u* is változik, hiszen a szélnyírás is változik

példa: erdő 30 m magas, fű 0,5 m magas

Példa: Hogyan alakul az u* erdő illetve gyep felett, ha feltesszük, hogy z=40 m magasan ugyanakkora a szélsebesség mindkét esetben. erdő

gyep

h

30 m

0,5 m

d

18 m

0,3 m

z0

3m

0.05 m

u z gyep u z erdő



u*gyep u*erdő

  z  d gyep     ln     z0, gyep   u*gyep     u*erdő   ln  z  d erdő     z0,erdő    

  40  0,3    ln  0,05   u    *gyep 3.39     40  18   u*erdő  ln  3    

Vagyis, ha ugyep = u erdő 40 méter magasan, akkor a hányadosuk 1, tehát:

1

u*gyep u*erdő

3.39

u*erdő  u*gyep 3.39

Fluxus: 1ségnyi idő alatt 1ségnyi felületen átáramló anyag mennyisége

Gradiens: adott mennyiség vmilyen irányú megváltozása Kiiindulás: Fick diffúziós törvénye: F  K

c z K: örvényes diffúzivitás [m2s-1]

u z

Momentum áram:

  a K m

Szenzibilis hőáram:

H  a c p K h

Látens hőáram: Szén-dioxid áram:

E 

Fco

2

 z

a c p e Kv  z

c  a K c z

A gradienseket meg tudjuk mérni, az örvényes diffúzivitást viszont becsülni kell Feltevés: Km=Kh=Kc=Kv

(Reynolds-féle hasonlósági elmélet = azonos forrás)

3 módszer ismeretes az örvényes diffúzivitás becslésére:

aerodinamikai módszer energia mérleg módszer „direkt” módszer Aerodinamikai módszer KM 

τ u ρ z

τ  ρu*2 u*  kz

u z

 u  k z   τ u*2  z   k 2 z 2 u KM    u u u z ρ  z z z 2

2 2

Szenzibilis hőáram: Látens hőáram:

H  a c p k 2 z 2

u  z z

ac p 2 2 u e E  k z  z z

Felszín – légkör kölcsönhatások

Szenzibilis hőáram:

Látens hőáram: (párolgás)

u  u T  T  ρc k u H  ρc k E    2

2

p

2

1

1

z d  ln z  d    

2 2

2

1

Ezeket tudjuk: : a sűrűség (1,2 kg m-3), cp: a levegő hőkapacitása (1005 J kg-1 K-1), k: von Kármán féle állandó (0,4), : pszichrometrikus állandó (0,65 mbar/°C) Ezeket megbecsüljük: z1, z2:a két szint felszín feletti magassága, d: kiszorítási rétegvastagság, (növényzet magasságának 60%-a)

p



2

 u1 e1  e2 

  z 2  d  ln z  d     1

Ezeket mérjük: u1, u2: a szélsebesség T1, T2: a hőmérséklet, e1, e2: gőznyomás

2

Energia mérleg módszer K 

Rn  G  e  a  c p    z  z  

Bowen arány módszer:

H  E

Ae  H  E Ennél a módszernél elég 2 szinten mérni a hőmérsékletet és a nedvességet, valamint egy szinten a sugárzási egyenleget. 

Ebből az áramok:

Ae E  1 

Visszahelyettesítve:

Ae H 1 

T H z   T   e E c e K  z

c K p

H

p

V

Feltesszük, hogy KH=KV

1

A rendelkezésre álló energia becsülhető, mint a sugárzási egyenleg 90%-a. 

Örvény-kovariancia módszer  Direkt árammérési technika a növényzet és a légkör közötti kölcsönhatás mérésére  Előnye: 24 órás mérés az év minden napján.

Mérés 10 Hz-es felbontásban  u,v,w  hőmérséklet  CO2  vízgőz  ózon  VOC  CH4

Szónikus anemométer:  3D szélmezőt méri

 Nincsenek mozgó alkatrészei  gyors válaszidejű  3 forrás / 3 detektor (ultrahang)  A kettő közti út megtételéhez szükséges időt méri  u, v, w, hangsebesség  Ts

Infravörös gázanalizátor (IRGA) A kibocsátott infravörös hullámok abszorpcióját méri.

Alul: forrás, fent: detektor Egyéb gyors válaszidejű szenzorok: O3, CH4, VOC

Adatgyűjtő berendezések és szoftverek

Reynolds átlagolás:

x  x  x

A jellemző pillanatnyi értéke (A) felírható, egy átlag (A), és az attól való eltérés (A’), az ún. fluktuáció összegeként.

 A  A  A  A  A  A  A

fluktuáció

átlag

A  0

A  B  A  B

x  x  x

cA  cA

A B   A B

Fluxus: a koncentráció és a vertikális sebesség szorzata







Fx   x w

Fx   x   w  w   x w   ' x w'   ' x w   x w' ' x

'

- az anyagmegmaradás elve miatt:

w0

- átlagolási szabályok alapján:

w  0

Fx   x w  Fx   x   x w  w 



Momentumáram:

  u *2   u ' w' v' w'

Szenzibilis hőáram:

H  c p w ' Ts '

Látens hőáram:

E  L w ' q'

CO2 áram:

FCO 2   w 'CO 2 '



1/ 2

Összetevőkre bontás Nettó ökoszisztéma széncsere (NEE): w – c + NEE: a növényzetből a légkör felé  légzés - NEE: a légkörből a növényzet felé  fotoszintézis dominál

NEE = -GPP + Reco éjjel: GPP=0  Reco

Fc  Rref e

1  1  E0     56, 02 t  46, 02 

t: a hőmérséklet °C-ban Rref: a referencia légzés 10°C-on E0: az aktivációs energia

NEE: Net Ecosystem Exchange