生物环境物理学特论EnvironmentalBiophysics-森林水文学.ppt

生物環境物理学特論Environmental Biophysics小杉緑子 第５回：群落H2O交換＆CO2交換 ―生態系フラックスデータの特性と使い道 今日の参考書 　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 「フラックス観測マニュアル」 フラックス観測マニュアル編集委員会編、2009　　　 http://www2.ffpri.affrc.go.jp/labs/flux/manual_j.html 「植物と微気象-群落大気の乱れとフラックス」 文字信貴、大阪公立大学共同出版会、2003 「陸域生態系における二酸化炭素等のフラックス観測の実際」　Asia Flux運営委員会編、2003 http://www-cger.nies.go.jp/publication/M015/M015.pdf 乱流変動法（渦相関法）：乱流フラックスを求める方法として、用いる仮定が少なくもっとも直接的な方法であると考えられている。CO2などの微量気体の鉛直フラックスは次式のように表される。 cは物質濃度、wは風速の鉛直成分でバーは時間平均、プライム(‘)は平均からの偏差を表している。平坦で一様な地表面でwの平均wが0と仮定できる場合、鉛直フラックスはwとcの共分散の平均　　で表される。 乱流フラックスとは 生物体と環境の間の運動量?エネルギー?物質の交換過程は概ね乱流輸送であり、潜熱(latent heat flux)?顕熱(sensible heat flux)?CO2などは乱流(eddy)によって輸送される。 乱流フラックスとは 風にのって運ばれるH2O分子 鉛直風速の向き 鉛直風速が下向きのときH2O濃度が高く、 鉛直風速が上向きのときH2O濃度が低ければ、 HO2は下向きに輸送される（H2O吸収＝凝結） 鉛直風速が下向きのときH2O濃度が低く、 鉛直風速が上向きのときH2O濃度が高ければ、 H2Oは上向きに輸送される（ H2O放出＝蒸発散） H:顕熱フラックス(W m-2)、Cp：空気の定圧比熱(J K-1 kg-1) lE:潜熱フラックス(W m-2)、 l：気化潜熱(J kg-1)、E ：蒸発散量（H2Oフラックス）(kg m-2　s-1)、r:空気の密度(kg m-3) 、w：風速(m s-1)、q:比湿(kg kg-1), rv:水蒸気密度(kg m-3) 乱流フラックスとは CO2もH2Oや顕熱と同様、乱流によって輸送される。 乱流変動法での群落CO2Flux測定＝NEEが直接測定可能。 風にのって運ばれるCO2分子 鉛直風速の向き 鉛直風速が下向きのときCO2濃度が高く、 鉛直風速が上向きのときCO2濃度が低ければ、 CO2は下向きに輸送される（CO2吸収） 鉛直風速が下向きのときCO2濃度が低く、 鉛直風速が上向きのときCO2濃度が高ければ、CO2は上向きに輸送される（CO2放出） Fc:CO2フラックス(吸収が正）(mmol m-2 s-1) rmol:空気の密度(mol m-3)、w：風速(m s-1)、c:CO2濃度(mol mol-1), rc:CO2密度（mol m-３） 生態系フラックスが算定されるまで （１）観測 まず樹冠上のある高度で、応答のよい測器を使い0.05～0.2秒間隔程度でx,y,z方向の風速、温度、CO2濃度、H2O濃度を連続測定する。 x,y,z方向の風速と温度の変動 を測定する超音波風速温度計 CO2,H2O濃度の変動を測定する オープンパス式ガスアナライザー （２）データの分割 フラックスを計算するためには、まず延々と続く生データをある長さで分割する必要がある。この時間のことを「平均化時間」という。 平均化時間は、短すぎると乱流フラックスに寄与する波の全てを捕らえることができず、長すぎると日変化のような大きな変動までひらってしまうので、 30分から一時間程度がよいとされる。 生態系フラックスが算定されるまで 2003/6/29 9:00-9:30の赤穂常緑広葉樹林のデータ （３）生データのチェックや異常値の除外 生態系フラックスが算定されるまで （４）座標回転 w’c’の平均からフラックスを算出するためには、水平一様な地形上でwの平均が0でなければならない。しかし実際には地形の影響などからｚ風速の平均は0にならない。これを補正するために座標を回転し、計算にはz風速の値ではなく、座標回転後のw風速を用いる。 Double rotation (wを毎回ゼロに） Planer fit (wが長期平均的にゼロになる一平面を仮定） 生態系フラックスが算定されるまで （５）タイムラグを見積もる 風速?温度とCO2密度?H2O密度は別々の測器で測