GFDセミナー2011夏

接地境界層乱流とフラックス観測玉川一郎（岐阜大学流域圏科学研究センター）

自己紹介
玉川一郎

tama@green.gifu-u.ac.jp
http://tama.green.gifu-u.ac.jp/~tama/
岐阜大学流域圏科学研究センター 勤務
1964年大阪市生まれ，岐阜市在住
京都大学理学研究科地球物理学専攻修了（防災研究所暴風雤災害部門），博士（理），気象予報士
職歴：
気象庁東京航空気象台予報課 （１９８９－１９９１）運輸技官
名古屋大学大気水圏科学研究所（1995-1999)助手
岐阜大学工学部土木工学科(1999-2002)助教授
岐阜大学流域圏科学研究センター（2002-現在）助教授（准教授）第一部• 境界層，大気境界層，接地境界層• Monin-Obukhov の相似則
• Kolmogorov の-5/3乗則• 超音波風速計，Kansas 実験 (Businger & Dyer
EQ.)
• Kadar & Yaglom 1990

第2部
• 渦相関法，熱収支，水収支，炭素収支• センサー• 渦相関法の仮定• 現実のデータ• その処理と補正，問題点
サンプリング，平均，座標軸，移流• Gap filling第一部• 境界層，大気境界層，接地境界層• 乱流輸送• 基礎方程式• Monin-Obukhov の相似則
• Kolmogorov の-5/3乗則，マイクロスケール• 超音波風速計，Kansas 実験 (Businger &Dyer)• Kadar & Yaglom 1990

境界層
流体力学における境界層

図：Wikipedia より流れを（乱流）粘性の影響のある層（境界層）とその外にある非粘性流体に分けて考える近似面から離れる方向のスケールが非常に小さいという近似をする

境界層方程式
U, ℓ, δv = O(

ℓ≫δ

で

δU
)
ℓ

U2

u 運動方程式の左辺 O( ) で，粘性項の第一項＜＜第二項なので，
ℓ

δ=O
δ
ℓ

1

∝ R−2 を無視する近似で

νℓ
U

1
2

1
−
= O(R 2 ℓ)境界層方程式
∂u ∂v
+
=0
∂x ∂y
∂u
∂u
∂u
1 ∂p
∂2 u
+u +v
=−
+ν 2
∂t
∂x
∂y
ρ ∂x
∂y
1 ∂p
0=−
ρ ∂y

外縁で u，全体に −

1 ∂p
ρ ∂x

を与えて

non slip で解く

巽友正「流体力学」より

大気境界層小倉義光 一般気象学より

大気境界層一般気象学より煙ー大気中における振る舞と姿
横山長之 著

日中の境界層 （混合層）

煙 より

夜間の境界層 （安定層）

煙 より

夜間の境界層 （安定層）

煙 より

接地境界層• 大気境界層の下端部，フラックスの鉛直方向変化が無視できる範囲 （constant flux layer)• 水平一様，定常の下で，コリオリ力を無視し，
気圧傾度力を無視
粘性，分子拡散 を無視する範囲でOK• 運動量フラックス変化10%くらいを目安• 0.1H大気境界層乱流• 超音波風速温度計の開発以降，接地境界層乱流の直接観測が盛んに行われた。
• 乱流の統計力学としての中心的な理論である
Monin-Obukhov の相似則は，それ以前，第2次大戦中に提唱されていた。
• 壁面乱流の対数則は，もっと前
von Kármán, Th. (1930), “Mechanische Ähnlichkeit und Turbulenz”,
Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen,
Fachgruppe 1 (Mathematik) 5: 58–76 (also as: “Mechanical
Similitude and Turbulence”, Tech. Mem. NACA, no. 611, 1931).
From Wikipedia

対数則• 水平一様定常，コリオリ，粘性を無視
0=−

∂u′w′
∂z

レイノルズ応力が支配。
−u′ w ′ = u∗ 2

長さの次元に z を持ってきて無次元化
z ∂u 1=u∗ ∂z k

u=

u∗
𝑧
log
k
𝑧0

u′ ~ w ′ ~ℓ

∂u
∂z

,

相関係数分も含めて ℓ = kzとして−u′w′ = (kz)2

∂u ∂u
∂z ∂z

Monin-Obukhov の相似則
• Turbulence in atmosphere with a non-uniform
temperature, A.M. Obukhov, Boundary-Layer
Meteorology, 2, 7-29, 1971 (Orig. 1946)(http://www.springerlink.com/content/p21902m3p810
6311/)• Basic laws of turbulent mixing in the surface layer of
the atmosphere, A.S.Monin&A.M.Obukhov, Tr. Akad.
Nauk SSSR Geiphiz. Inst. 24(151):163-187,1954(http://mcnaughty.com/keith/papers/Monin_and_Obu
khov_1954.pdf ) 最初の実験的証拠
渦粘性係数 t  u * (kz) の対数則の拡張を考えた

Monin-Obukhovの相似則

3

u*
L
g
k   w' v ' v 

gkz  w' v '

z
v 

M-Oの安定度3Lu*u*u*  u' w'

スケーリングパラメタ

u*T*  w' 'uは主流方向，wは鉛直
Tv は仮温度，qは比湿 （CO2等でも同様）
フラックスは一定とみなせる
水平一様定常の仮定
境界層近似


  
 , ,zx y tu*Tv*  w' v 'u*q*  w' q'

Monin-Obukhovの相似則• フラックス―勾配関係
kz ∂u
z
= ϕm
u∗ ∂z
L

kz ∂θ
z
= ϕh
T∗ ∂z
L

• 変動量
w

z
 w  
u*
L

T

z
 T  
T*
L

kz ∂q
= ϕh
q ∗ ∂zz
L

この関係は，観測や（主に極限での）理論的考察によって確か
められ，広く気象モデルの下端
境界条件に使われている。

• 他 消散率，リチャードソン数，スペクトルも  kz / u*3

 z   
Ri    h2  L  m 

z/L は• MOの安定度(z/L)は，水平一様定常のもとでの，乱流運動エネルギー方程式の浮力消滅項／シアー生成項 に相当する。
gkz  w' v '
v 
z


3
L
u*

植物と微気象（文字2003）より

超音波風速温度計
Sonic Anemometer Thermometer•Send and receive sonic pulse between two proves•Count the time
•Calculate average sound speed and the difference ofthe speed•Calc. wind speed of sound path direction•Calc. temperature from average speed (sonic virtual
temperature)
Developped in1960s• Kaimal, J. C. and Businger, J. A.: 1963, A Continuous Wave Sonic
Anemometer-Thermometer, J. Appl. Meteorol. 2, 156–164.
• Mitsuta, Y.: 1966, Sonic Anemometer for General Use, J. Meteorol.
Soc. Jap. 44, 12–24.

今や トラ技に制作記事， 高校生が作って観測。
極寒冷地では，ヒーターをつける。降水くらいなら観測可→ATI→Kaijo

長期安定は難しい超音波風速計• 音波を往復させ，その音速の差から媒体である大気等の流速を計測する装置• 超音波パルスを使った実用システムは，光田＆海上電機（現 株式会社ソニック）の開発である。
• Mitsuta, Y.: 1966, Sonic Anemometer for
General Use, J. Meteorol. Soc. Jap. 44,
12–24.

First Commercial
Sonic Anemometer超音波風速計• 応答に優れることから，乱流観測に使用• 可動部がないので，メンテナンスフリーの風速計として• 原理的に直接風速をはかるので高精度• 近年，非常に安価なタイプが出現• 2次元タイプのものが，既存の風向風速計を置き換えつつある。
• 分解能と強風対応はトレードオフ

製品例

Vaisala のWWWページよりソニックのwwwページよりDavis のWWWページより

超音波風速計の使用例• トンネル，橋，タワーなどの風速監視• 室内，スポーツ競技場の風速• パイプの流速 （最近は電磁流速計やレーザー
流速計にとってかわられている）• 通常気象観測（プロペラや風杯を置き換えつつある）
大気乱流の観測• 地上／船舶上のポール，あるいはタワー• 航空機（低速）他の大気乱流観測用風速計• 熱線風速計 （主に実験室）• ピトー管 ５孔ピトー管なら3次元で，乱流変動まで計測可能。航空機に使われるが，地上用に販売していたこともある。RosemountAerospace（当時）（ピトー管と圧力センサーの間に水滴受けをつけるなどの工夫が必要）

SpaceAge Control社

http://www.spaceagecontrol.com/Main/Home より

Flux measurement
Sonic Anemometer Thermometer
Sonic temperature includes cross wind effect andwater vapor effect.

音仮温度と気象で使う仮温度はわずかに係数が違う

Kaimal and Finnigan 1994 “Atmospheric boundary layer flows:
their structure and measurement”

風速計自身の影響
Sonic Anemometer Thermometer
Also the effect of probes and other structure are also
the source of erros

In wind tunnel
伊藤，林，玉川“超音波風速TR-61A型プローブ特性試験
京都大学防災研究所 共同研究１２G-5
熱・水収支観測の高精度評価に関する研究
Hanafusa et.al.1982, Meteorol, Geophys.33,1-19
Kaimal and Gaynor 1983, J. Appl. Met.22, 863-880

In natural turbulent wind
( unpub.)

つづきAlso the effect of probes
and other structure arealso the source of erros

The amount of these effect isdifferent for each probe design,
so please find suitable articles
or make test by yourself.

=> Deep world of measurement

大規模な野外実験
Kansas 実験 (Haugen et. al. 1971, QJRMS)• Flux-Profile 関係 (Businger, Wyngaard, Izumi
and Bradley 1971, JAS)
32mタワー，18cmの wheat stubble, 2400mのフェッチ
風速の鉛直分布は，23, 36, 50, 70, 100cm, 2, 4, 5.66, 8, 11.3, 16,22.6, 32m
温度は，0.5, 1, 2, 4, 8, 16, 22.6, 32m

現在使われている主なMOS式はこのデータから，と言ってもよいFlux-Profile
Businger et. al. 1971 の解析では，無次元シアー関数が中立に近い範囲 (-2<z/L<2)
で，実験的に求められた。ただ，カルマン定数
が k=0.35 であるなどその後の議論を呼んだ

Businger et. al. 1971 よりMast の影響，風杯の回りすぎなどが懸念された。

Wyngaard, Businger, Kaimal and Larsen,
Comments on ‘A revaluation of the Kansas
近年は，k=0.4で整理された
無次元シアー関数を使うことが多い。 Mast Influence on Measurements of stress andcup anemometer overspeeding’, Boundary（水理学では，0.41が普通）Layer Meteorology, 22, 245-250, 1982 より

再掲

Monin-Obukhovの相似則• フラックス―勾配関係
kz 
z
 h  
T * z
L

kz U
z
 m  u * z
 L

• 変動量
w

z
 w  
u*
L

T

z
 T  
T*
 L

kz U
z
 h  
q * zL

この関係は，観測や（主に極限での）理論的考察によって確か
められ，広く気象モデルの下端
境界条件に使われている。

• 他 消散率，リチャードソン数，スペクトルも  kz / u *3

 z   h 
Ri    2 
 L   m 

文字２００３より
文字信貴 2003「植物と微気象」
群落大気の乱れとフラックス
大阪公立大学共同出版会

文字２００３より

文字２００３より

文字２００３より

文字２００３より

文字２００３より

このあたり、近年の教科書では
近藤裕昭「人間空間の気象学」２００１が詳しい

Kadar&Yaglom 1990

上図で ζ = (1/k) (z/L) カルマン定数を含まない，MO安定度
水平方向と鉛直方向を分ける次元解析を行い，３つの領域に分けた。
2
ちなみに [k] = (Lz/Lx)1/2 [z/L]=(Lz/Lx)
Kadar BA and Yaglom AM (1990). Mean fields and fluctuation moments in unstably
stratified turbulent boundary layers. J Fluid Mech 212: 637–662

変動量

文字２００３より

文字２００３より

文字２００３より

変動に関して(手元のデータより）• 不安定の極限で，u*によらない形で，無次元化された標準偏差がz/Lの関数になる。

Ｗ

不安定

安定

T

Tamagawa 1996
BLM 77, 1-20

変動 （U ｑ），怪しくなってくる

Ｑ

不安定
安定
U

Tamagawa 1996

あやしくなるのは，• 仮定の水平一様定常が崩れていく
砂漠上での観測 -> 水蒸気の鉛直微分が小さい＝相対的に水平方向の影響を受けやすい
U は，境界層全体に及ぶような循環の影響を受ける （Zi の影響を受けるという解析結果： だれだっけ？ 乱流として解析する周波数帯に依存）
メソスケール循環との重ね合わせとして理解しようとする方法もある。
もとに戻って，どういう変動が乱流変動か？ → スペクトルスペクトル （文字2003）

-5/3乗則 （文字2003より）

Power Spectrum
慣性小領域も見える

Kansas実験の結果
Kaimal and Finnigan 1994
教科書Atmospheric Boundary Layer Flowsより

文字2003

Co spectrum

Kaimal and Finnigan 1994

乱流の研究として ７０－８０‘ｓ• MO相似則は，確かに成立 ＝＞ 密度成層流体の乱流統計理論， でも，安定成層時は，今でも難しい。
(z less profile, wave, intermittent turbulence)• 水平一様定常仮定を極力満たすように観測条件を設定。観測場所のみ
ならず，時間も，たとえば，朝夕は観測せず日中，夜間のみ観測して統計するなどしていた。
• 他に，粗度高など地面の空気力学的条件• もちろん，レイノルズ応力や顕熱フラックスは計測された。熱収支の直接観測はまだ。潜熱は残差• 他に，乱流の中に組織的運動があることも，当時から指摘されていた。
立体的な理解はLES、ライダー観測などまで待つ必要がある。
• 慣性小領域でのコルモゴロフの－5／3乗則の観測，コルモゴロフ定数からLESのスマゴリンスキー定数。対数則から k-ε のモデル定数など，基礎的な定数を大気でも確認した （と思う）。風洞とはレイノルズ数が違う。

粗度高• 対数則（＝中立時のMO則）を積分（MO則でも同様）

u*u2  u1 log( z 2 / z1)
あてはめて，k u1 = 0 となる高度z1を z0 とする。

地面の特性量になるとされている。
• 温度に対しても同様にする。地面の温度を放射温度計で計測したもの等に置き換えて処理 （圧力による抵抗，ソースの位置の違い，表面値の定義の相違等でかなり違う値になる）• 森林のような高い粗度物体では z の原点を変える必要がある
z/z0  (z-d)/z0• キャノピーモデルを使う現在でも，キャノピー直上でこの値が使われている。
• しかし，z0, d は大気安定度の関数であるらしい(例えばZilitinkevich et. al. 2008)
Zilitinkevich et. al. , The Effect of Stratification on the Aerodynamic
Roughness Length and Displacement Height, Boundary-Layer Meteorol (2008) 129:179–190
DOI 10.1007/s10546-008-9307-9粗度高ほか上、人間環境の気象学 近藤裕昭著より
数値計算に使用したパラメタ
右、基礎気象学 浅井冨雄、新田尚、松野太郎著
より

第2部• 渦相関法，熱収支，水収支，炭素収支• センサー• 渦相関法の仮定• 現実のデータ• その処理と補正，問題点
サンプリング，平均，座標軸，移流• Gap filling地表面熱収支λλ

この内のH,LEが乱流観測から評価されるもちろん。このEは蒸発量なので水収支にも現れる。

また、炭素収支にも、乱流変動によるCO2輸送項がでる。

Observed data (daytime)
u

T

v

c

w

q

TKC での観測値の例

Observed data (nighttime)
u

T

v

c

w

q

TKC での観測値の例

渦相関法による乱流輸送量推定
Eddy covariance:

Then assume horizontal
homogeneity H  C p w'T 'E   w'q'

顕熱フラックス
1 R
d    dT
 
とすると

c p  cv  R

ちょっとこの辺で密度変動を無視して連続式を使って
る

より



T

 T
u j 'T '  J
cp 
uj

x j x j


 t


J: heating to unit mass

これを，両辺にρをかけて，地表におかれた単位面積のカラムで積分して，水平一様，定常，加熱は下端のみとすれば，H  C p w'T 'Holton, An introduction to Dynamic Meteorology 4th ed. 2004

Montgomery 1948, J.Met.5,265-274 に多分最初期の議論がある。

顕熱、潜熱フラックスの推定• 1980年代，水蒸気の乱流変動観測が可能になってきた。
• 細線熱電対乾湿球温度計が使われたが，その後，赤外線あるいは紫外線の吸収を利用した
Gas Analyzer が使われるようになった。
(Ohtaki and Matsui 1998 BLM 24, 109-119)
現在，LI-COR社の製品(LI-7000, LI-7500)が定番

結果，顕熱と潜熱のフラックスが測れるようになったので，熱収支観測が盛んに行われるようになった。その後，水蒸気，メタンなども可能になったフラックス測定 Field Projects• 多分， FIFE (The First ISLSCP Field
Experiment ) あたりが最初(1987-1989)。
http://daac.ornl.gov/FIFE/fife.shtml• 他に，TOGA-CORE, HEIFE 1991,92に集中観測(http://ssrs.dpri.kyoto-u.ac.jp /~heife/), BOREAS 1994,96に集中観測 (http://daac.ornl.gov/BOREAS/bhs/BOREAS_Home.html) など多数• GEWEX/GAMEでも，多数の地点で乱流観測が行われた。

START OF MY RESEARCH

HEIFE
http://ssrs.dpri.kyoto-u.ac.jp/~heife/(Heihe River Basin Field Experiment)
Observation of heat
and water exchangeover desert.

Mitsuta et.al. 1995

1991 Jun, Aug, Oct, Dec

I was a student!

HEIFE 1989-1993の乱流観測の様子

年間の熱収支、水収支• 放射、平均場、地中温度の観測は通年• 乱流観測は集中観測時のみ（有人）（記録計と電源の問題が大きかった）• Profile-flux関係式を現地で作成して年間に適用Mitsuta et. al. 1995
乾燥地なので水蒸気のバルク式はうまく作成できず
残差で推定

GAME/Tibet
1998 Intensive Observation

安多I belonged to Nagoya univ.

Tanaka et.al. 2001

観測装置
• 超音波風速温度計 （昔は，ほとんど Kaijo だった。最近は，Gillが多いが，いろいろなメーカーが使われる）• 赤外線Gas Analyzer (LI_COR LI-7500 オープンパ
ス，あるいは 同 LI-7000クローズドパス）• 昔は，ここに細線抵抗（あるいは熱電対）温度計を付けたが最近はあまり付けない。
• 必要に応じて，傾斜計，加速度センサー• 共分散のみが目的なら不要だが，スペクトル解析のためには，ローパスフィルターも必要。（電源ノイズ対策でもある）今時ならオーバーサンプリングして演算処理でも。

Flux measurement
Infrared H2O CO2 analyzer
Mirror Measure the relative light decay comparing theinfrared rays in absorption band and in anotherband.
Measure the amount of gas in the light path, inother words, density of gases.
Old generation sensors were unstable.
Light source & detecter
It also has deep world!

Open path 型は、２０００年頃、実用モデル
それ以前は、実用型としては、吸引して分析する
closed path型しかなかった。

石田 祐宣, 松島 大, 樋口 篤志, 檜山 哲哉, 戸田 求, 浅沼 順, 玉川 一郎, 宮崎 真, 田中 賢治, 杉田
倫明, 永井 秀幸, 田中 久則, 飯田 真一, 小林 菜花子: “2001年筑波大学陸域環境研究センター(TERC)における乱流計測機集中観測:機器比較と校正による誤差の解析”, 水文・水資源学会誌,
Vol. 17, No. 1, pp.43-60, (2004) .

Advanet E009
Kaijo AH300

２００１TERC比較観測の写真（石田＠弘前大）

LI-COR 社のWWWより

定番センサー

年１回程度の較正で安定して観測可CAMPBELL社のWWWより

Flux measurement
Infrared H2O CO2 analyzer
MirrorIt also has deep world! Ex. Heating by sensor itself

Light source & detecter

筺体の熱の影響の補正

記録に関して• 1990年代以前は，統計値のみ，あるいは 30
分前後の期間の乱流生データを保存。
• それ以降は，連続的に乱流データを保存するのが普通。
• サンプリングは，10Hz 程度• 例えば，90年代のFlux-PAM
は，10分毎に平均値と共分散を記録，GAME-Siberia では連続記録。
（太陽電池駆動）

耐久性に優れる専用ロガー

容量等性能が高いPC利用

Flux-PAM

NCARとGAME(筑波大学)共同開発。太陽電池で動作。
残念ながら当時とてもよく壊れた。記録制御部の問題

乱流フラックス計測• サンプリング間隔，平均時間，トレンド除去？
乱流変動のコスペクトル，計測器の安定性，応答性• Ｗの向き，座標系設定
更に現実的には• 密度変動による質量フラックスに対する補正• 鉛直風• 乱流以外の輸送（夜間の斜面流などが特に問題）

平均• コスペクトルを考慮して，適切な時間で平均処理。
• 乱流（特に，不安定時のプリューム)が十分平均されて統計的に安定する程度に長く，日変化等の影響を受けない程度に短く[ｗには影
響しないはずなので，大丈夫？]• 高周波側も十分取得する必要
Kaimal and Finnigan 1994

From GAME-Tibet 1998

Covariance w and T

Covariance w and mixing ratio ofwater vapor

30分以上なら，あまり問題はないが，水蒸気が雲活動の影響が見える

Unpub. Data by meトレンド除去• センサーのトレンドがあれば不可欠• 長周期の変動を取り除き、乱流だけにする機能• しかし、論理的には、トレンド間で共分散＝フラックスがある• １次、または２次のトレンド、High pass filter など使われたが、どちらかと言うと、ブロック平均が好まれる。
このあたり、Lee et. al. ed. 2004, Handbook of. Micrometeorology に詳しい

座標系
• 地表面を覆う面で積分するのが多分正しいが，現実には１点のタワー観測• ３つの選択1run毎に，ｚ軸周りに回転して v  0 & y’軸周りに回転して w  0
1runごとに w  0,
かつ v' w'  0v0長期間の平均風の乗る平面を求める(planar fit)
後述の平均鉛直風の問題や物理的理解のしやすさから好まれるが実際に実行できるとは限らない密度変動補正(WPL補正）
文字２００３より

密度変動補正(WPL補正）
文字２００３より

WPL補正• ｗの計測誤差以下の平均鉛直風を評価• 実際の計測で見える平均鉛直風は無視することが多い。
（単純に wT wq を計算すると非常に大きなフラックス値がでて明らかにおかしいことが多い。 大規模渦の片側のみの計測になるのでしないほうがまし。センサーやタワーのせいで出た鉛直風の可能性もある。）

鉛直風の補正 (Lee 1998)• 鉛直移流の効果 （planar fit によるw)

c p w Ta  T0



水平移流では温度変化なしと考えるなどの補正を加えると、熱収支がよく閉じるという報告があるが、年単位での解析例は多分ない
Lee, X.: 1998, ‘On Micrometeorological Observations of Surface-Air Exchange over Tall Vegetation’,
Agric. For. Meteorol. 91, 39–49.

精度良く計測するには• Ｗ の方向 （double rotation, planer fitting)• WPL 補正 （密度変動による乾燥空気の質量輸送の補償平均風を考慮） CO2等では輸送の正負まで変わることもある。
• 風速計自身の影響などセンサーの詳細な補正• 森林等では，計測高度以下の貯留や非一様性による鉛直風の考慮
それでも，森林の場合 (Rn-G)/(H+LE)は 0.7 – 0.9 程度。 たとえばFoken, T. 2008. The energy balance closure problem:
an overview. Ecol. Appl. 18, 1351-1367.

そもそも１点では。。。話も•神田 学，渡辺 力，マルコス オリバー レッツェル，ジークフリード ラッシュ：
ＬＥＳによるインバランス問題に対する検討（第１報）大気境界層スケールの対流構造の影響，水文・水資源学会誌, 15, 243-252, 2002.

点計測の時系列データに基礎を置く限り,測定点数を増やして厳密に空間平均をとっても,その乱流フラックスは真の水平平均乱流フ
ラックスを過小評価する(負のインバランスが生じる).これは,計測点に於けるローカル
な平均鉛直流の熱輸送寄与によるものである

境界層内の準定常的循環の影響

Kanda et. al. 2004

Kanda M, Inagaki A, Letzel MO, Raasch S (2004) LES study of the energy imbalance
problem witheddy covariance fluxes. Boundary-Layer Meteorol 110:381–404

データ処理の流れ• データを保存 （1 run = 30分程度)
この検討のために全データ保存• 不良データ，ノイズ除去• 座標回転，補正
逆にすることも多い• 平均，共分散等計算• 運動量，顕熱，潜熱，CO2交換量など算出
この後，目的に応じてデータを選別，欠けているところを実験式で補って，年間の熱収支，水収支，
炭素収支の解析へ最近の観測例として
岐阜大学高山常緑針葉樹林サイト

TKC Tower•Electric power supply•Information network

Information Network
TKY site
OfficeTKC site

Takayama AP
Gifu University
The Internet

Gifu Information Super Highway
by Gifu prefecture

Use PC at TKC from my laboratory
in University

PC in university

PC in TKC site

VNC server and viewer

WWW camera & PEN

Sensors on Tower (1/2)

Sensors on Tower (2/2)

Also flux measuring system• Sonic anemometerand thermometer(R3-50, Gill)• Infrared
H2O
and CO2 analyzer
LI-7500(LI-COR)

Observed data (daytime)
u

T

v

c

w

q

TKC での観測値の例

Observed data (nighttime)
u

T

v

c

w

q

TKC での観測値の例

乱流測定以外に平均場の計測も乱流測定以外に、平均場の鉛直分布を計測•
乱流輸送を決める平均場•
渦相関フラックス観測高度以下でのソース／シンク•
そもそも、森林の研究のために必要

CO2 profile measuring system

Heat budget (1 dim.)• Heat budget must be satisfied by energy
conservation (if horizontallyhomogeneous).

Radiation observation

Ground observation

So many sensors for profiles

Saitoh et. al. の紹介
•

•
•

•

Saitoh T.M., Tamagawa, I., Muraoka H., Lee N-Y., Yashiro Y., and Koizumi H. (2010)
Carbon dioxide exchange in a cool-temperate evergreen coniferous forest overcomplex topography in Japan during two years with contrasting climates. Journalof Plant Research doi:10.1007/s10265-009-0308-7.
斎藤琢・玉川一郎・村岡裕由・小泉博, 森林における光合成・呼吸活動に伴う貯熱量, 農業気象66巻4号,289-298,2010
Taku M. Saitoh, Ichiro Tamagawa, Hiroyuki Muraoka, and Hiroshi Koizumi, Energybalance closure over a cool temperate forest in steeply sloping topography duringsnowfall and snow-free periods, Journal of Agr. Met., 2011,
doi:10.2480/agrmet.67.3.4
斎藤琢・玉川一郎・村岡裕由, 渦相関法を用いた炭素収支評価におけるCO2貯留
変化量の影響 日本森林学会誌, 2011

CO2フラックス• 渦相関のよいデータを選ぶ（平均湿度計測の一致、風向の選択、定常性検査）• 夜間のデータの内、 u* の大きいものを使って温度と呼吸の関係式• 森林光合成と光合成有効放射の関係式を渦相
関データから作成• 林内の貯留• 年間を埋める推定値ー＞年間の収支

Saitoh T.M., Tamagawa, I., Muraoka H., Lee N-Y., Yashiro Y., and Koizumi H. (2010)
Carbon dioxide exchange in a cool-temperate evergreen coniferous forest overcomplex topography in Japan during two years with contrasting climates. Journalof Plant Research doi:10.1007/s10265-009-0308-7.

熱収支• 光合成、呼吸による熱の出入りは、問題になるかならないかぎりぎりくらい。林内の顕熱、潜熱の貯留の次に問題になる大きさ• 斜面なので、鉛直を図る放射観測と乱流フラックス観測の方向が違って、結果、面積が違う。
• 十分補正して考えると、平坦な森林と同程度の H+LE があるべき値の 70% 程度までとらえられている。
• Corrections of the storage term, the AOA error, and the inclined
surface improved EBR by 2.0–2.6%, 8.3–8.9% and 4.9–5.9%, andimproved CR by 5.3–5.8%, 5.6–7.9%, and 4.0–5.4%. Consequently,
EBRand CR reached 0.86 and 0.79 during the snow-free period. By
also considering the heat of fusion, EBR reached 0.79 during thesnowfall period. Our findings suggest that heat fluxes of eddycovariance measurements over a forest on steep slopes have asimilar accuracy as those over other topography.
Taku M. Saitoh, Ichiro Tamagawa, Hiroyuki Muraoka, and Hiroshi Koizumi, Energybalance closure over a cool temperate forest in steeply sloping topography duringsnowfall and snow-free periods, Journal of Agr. Met., 2011,
doi:10.2480/agrmet.67.3.4

近年の観測• 近年の乱流観測は，おもに大気と地面との間の交換量の観測を目的としている。
• 風速，気温，湿度，CO2濃度等の鉛直分布と，その交換に関わる，放射量，地温，水分量など多くの観測とともに行われる。
• 観測したということ自身が重要な未知の場所は，地球では多分もうない。（極域，チベット，ヒマラヤ，砂漠，
海洋上，森林，都市）• 多数地点の観測もすでに CASES-97, 99(http://www.cora.nwra.com/cases/CASES-99.html), CAPS2002, 03(http://tama.green.gifu-u.ac.jp/~tama/flux/ ), ADVEX などされている。

• FLUXNETによる陸面モニター （炭素が中心的話題）

Asia Flux FLUXNETのアジア版• Mizuguchi et. al. 2009, J. For
Res. 各サイトの計測開始時期
合計109サイト

森林サイト
森林以外のサイト
ただし、すでに論文化されているサイトは少ない!!
年積算NEE（またはNEP）の報告があるサイトは
森林サイトで 23
非森林サイトで 8

観測のマニュアル化• 乱流などの気象系観測の研究者が少ない• 森林の方は、観測系研究者が主流• 森林研究者のタワー観測への参加

データは自分で取らないとという文化

マニュアル化、データ処理の標準化（ブラック
ボックス化でもある）の進行
多地点比較のために品質をそろえる意味もあるタワーフラックス観測マニュアル内容
1. 観測計画とフラックス観測サイトの選定
1.1 観測サイトの選定
1.2 インフラの整備
1.3 観測項目の選定

2. 乱流系観測http://www.ffpri.affrc.go.jp/
labs/flux/manual/manual_in
dex_j.html
現在ver1.0 が公開中

環境省地球環境保全試験研究プロ
ジェクト 「アジア陸域炭素循環観測のための長期生態系モニタリングとデー
タのネットワーク化促進に関する研究」
森林総合研究所、農業環境技術
研究所、産業技術総合研究所、国立環境研究所の四機関合同

2.1 超音波風速温度計
2.2 オープンパス型CO2/H2O分析計
2.3 クローズドパス型CO2分析計
2.4 貯留変化量
2.5 簡易渦集積法
2.6 データロガー
2.7 ノイズのチェックと対策

3. 微気象観測
3.1 放射3.2 風向・風速
3.3 気温3.4 湿度3.5 地温・地中熱流量
3.6 土壌水分
3.7 降水量・積雪調査
3.8 水位・水温・灌漑・排水量
3.9 データロガー

PP166

内容例 （ノイズ対策）
ノイズの問題は対策を施そうとしても，一体どこをどのように対策してよいかが非常に難しい問題となる場合が多くある。ノイズに強い観測システムの設計のためには以下の点に注意すべきである。
1) 計測器の電源には，電力機器用の電源とは別系統の電源を用いる。ノイズカットトランスなどを用いて分離することが有効である。
2) トランスやフィルタを使用した時は，1次側と2次側のラインは接近しないように分離して配線する。
3) 電力ラインと信号ラインをなるべく離す。やむを得ず交叉する場合は直交させる。観測小屋や観測タワーの配線では違う取り入れ口や鉛直パイプを使い，両者の配線位置を明確に分離する。
4) 信号線はできるだけ短くし，余りの部分をループ状にしない。

5) 測器やデータロガー，信号線などを，ノイズを発生していると思われる機器からなるべく離す。
6) アースの取り方には充分注意する。計測器同士のアースをつないで，しっかりしたアース端子に接続する。
タワーフラックス観測マニュアルv1.0 pp88

内容例 （Tips!）
3.5 地温・地中熱流量
Tips!
温度センサと一体となった防水型のデータロガーの内部にシリカゲルを入れてしっ
かりと密閉すると，基盤の結露による不具合を防ぐことが出来る｡シリカゲルは適時，交換する｡
Tips!
スパイラルチューブあるいは樹脂製の保護管でケーブルを保護することにより，ネズミなどに噛み切られることを防ぐことが出来る｡
Tips 3.5-2
Tips!
１ヶ所で複数の深度の地温を測定する場合，センサを埋設した後には，どのケー
ブルがどの深度のセンサのものか解らなくなる｡そのため地表面付近やロガー結線部などに，測定深度を記したタグをつけておくと，その後のメンテナンス時に便利である｡
Tips 3.5-1

タワーフラックス観測マニュアルv1.0 pp133

近年の観測まとめ• 乱流と通常気象観測だけではなく，総合的観測が求められている• マニュアル化，標準化が進んで，多くの研究者が参入できるようになっている。
• 観測目的は，地面 （植生） の活動のモニター• モデル化および検証のためのデータ，衛星観測の地上検証データに使われる。
• 多地点を集めて，地球のモニター

終ありがとうございました

参考文献
•

•

•

•

•

•

Asanuma, J., Brutsaert, W., 1999: Turbulence variance characteristics oftemperature and humidity in the unstable atmospheric surface layer
above a variable pine forest, Water Resources Res., 35, 515-521.
Businger, J. A., 1988: A note on the Businger-Dyer profiles, Bound. Layer
Meteor., 42, 145-151.
Businger, J. A., Wyngaard, J. C., Izumi, Y., Bradley, E. F., 1971: Flux-profile
relationships in the atmospheric surface layer, J. Atmos. Sci., 28, 181-189.
Dyer, A.J., 1974: A review of flux-profile relationships, Bound. Layer
Meteor., 7, 363-372.
Foken, T., 2008: The energy balance closure problem: an overview, Ecol.
Appl., 18, 1351-1367.
Grelle, A., Burba, G., 2007: Fine-wire thermometer to correct CO2 fluxes
by open-path analyzers for artificial density fluctuations, Agric. For.
Meteorol., 147, 48-57.

•

•

•

•

•

•

Hanafusa, T., Fujitani, T., Kobori, Y., Mitsuta, Y., 1982: A new type sonicanemometer-thermometer for field operation, Papers in Meteorology andGeophysics, 33, 1-19.
Haugen, D. A., Kaimal, J. C., Bradley, E. F., 1971: An experimental study of
Reynolds stress and heat flux in the atmospheric surface layer, Quart. J.
Roy. Meteor. Soc., 97, 168-180.
石田祐宣, 松島大, 樋口篤志, 檜山哲哉, 戸田求, 浅沼順, 玉川一郎, 宮崎真, 田中賢治, 杉田倫明, 永井秀幸, 田中久則, 飯田真一, 小林菜花子,
2004: 2001 年筑波大学陸域環境研究センター (TERC) における乱流計測
機集中観測: 機器比較と校正による誤差の解析, 水文・水資源学会誌,
17, 43-60.
Kader, B. A., Yaglom, A. M., 1990: Mean fields and fluctuation moments inunstably stratified turbulent boundary layers, J. Fluid Mech., 212, 637-662.
Kaimal, J. C., Businger, J. A., 1963: A continuous wave sonic anemometerthermometer, J. Appl. Meteor., 2, 156-164.
Kaimal, J. C., Gaynor, J. E., 1983: The boulder atmospheric observatory, J.
Appl. Meteor., 22, 863-880.

•

•

•

•

•

•

Kanda, M., Inagaki, A., Letzel, M. O., Raasch, S., Watanabe, T., 2004: LES
study of the energy imbalance problem with eddy covariance fluxes,
Bound. Layer Meteor., 110, 381–404.
神田学, 渡辺力, マルコス・オリバー・レッツェル, ジークフリード・ラッシュ,
2002: LES によるインバランス問題に対する検討 (第 1 報) 大気境界層スケールの対流構造の影響, 水文・水資源学会誌, 15, 243-252.
Lee, X., 1998: On micrometeorological observations of surface-air
exchange over tall vegetation, Agric. For. Meteorol., 91, 39–49.
Mitsuta, Y., 1966: Sonic anemometer-thermometer for general use, J. Met.
Soc. Japan, 44, 12-24.
Mitsuta, Y., Tamagawa, I., Sahashi, K., Wang, J., 1995: Estimation of annual
evaporation from the Linze desert during HEIFE, J. Met. Soc. Japan, 73,
967-974.
Mizoguchi, Y., Miyata, A., Ohtani, Y., Hirata, R., Yuta, S., 2009: A review of
tower flux observation sites in Asia, J. For Res., 14, 1-9.

•

•

•

•

•

•

Monin, A. S. and Obukhov, A. M., 1954: Basic laws of turbulent mixing inthe surface layer of the atmosphere, Tr. Akad. Nauk SSSR Geiphiz. Inst., 24,
163-187.
Monji, N., 1973: Budgets of turbulent energy and temperature variance inthe transition zone from forced to free convection, J. Met. Soc. Japan, 51,
133-145.
Monji, N., 1975: Characteristics of the horizontal wind fluctuatioins in thesurface layer under strong convective conditions, J. Met. Soc. Japan, 53,
99-102.
Obukhov, A. M., 1971: Turbulence in an atmosphere with a non-uniform
temperature, Bound. Layer Meteor., 2, 7-29.
Ohtaki, E., 1985: On the similarity in atmospheric fluctuations of carbon
dioxide, water vapor and temperature over vegetated fields, Bound. Layer
Meteor., 32, 25-37.
Ohtaki, E., Matsui, T., 1982: Infrared device for simultaneous
measurement of fluctuations of atmospheric carbon dioxide and water
vapor, Bound. Layer Meteor., 24, 109-119.

•

•

•

•

•

•

Panofsky, H. A., Tennekes, H., Lenschow, D. H., Wyngaard, J. C., 1977: The
characteristics of turbulent velocity components in the surface layer underconvective conditions, Bound. Layer Meteor., 11, 355--361.
Paulson, C. A., 1967: Profiles of wind speed, temperature and humidity
over the sea, Ph.D. dissertation, Univ. of Washington.
Paulson, C. A., 1970: The mathmatical representation of wind speed andtemperature profiles in the unstable atmospheric surface layer, J. Appl.
Meteor., 9, 857—861.
斎藤琢, 玉川一郎, 村岡裕由, 2011: 渦相関法を用いた炭素収支評価における CO2 貯留変化量の影響, 日本森林学会誌, 93, 105-112.
斎藤琢, 玉川一郎, 村岡裕由, 小泉博, 2010: 森林における光合成・呼吸活動に伴う貯熱量, 農業気象, 66, 289-298.
Saitoh, T. M., Tamagawa, I., Muraoka, H., Koizumi, H., 2011: Energybalance closure over a cool temperate forest in steeply sloping topographyduring snowfall and snow-free periods, Journal of Agr. Met.

•

•

•

•

•

•

Saitoh, T.M., Tamagawa, I., Muraoka, H., Lee, N. -Y., Yashiro, Y., Koizumi, H.,
2010: Carbon dioxide exchange in a cool-temperate evergreen coniferous
forest over complex topography in Japan during two years with
contrasting climates, Journal of Plant Research, 123, 473-483.
Tamagawa, I., 1996: Turbulent characteristics and bulk transfer coefficientsover the desert in the HEIFE area, Bound. Layer Meteor., 77, 1-20.
Tanaka, K., Ishikawa, H., Hayashi, T., Tamagawa, I., Ma, Y., 2001: Surface
energy budget at Amdo on the Tibetan Plateau using GAME/Tibet IOP98
data, J. Met. Soc. Japan, 79, 505-517.
Tillman, J. E., 1972: The indirect determination of stability, heat andmomentum fluxes in the atmospheric boundary layer from simple scalarvariables during dry unstable conditions, J. Appl. Meteor., 11, 783-792.
Karman, Th. v., 1931: Mechanical Similitude and Turbulence, NACA Tech.
Mem., 611, 1-19.
Webb, E. K., 1970: Profile relationships: the log-linear range, and extention
to strong stability, Quart. J. Roy. Meteor. Soc., 96, 67-90.

•

•

•

•

•

•

•

Wyngaard, J. C., Businger, J. A., Kaimal, J. C., Larsen, S. E., 1982: Comments
on 'a revaluation of the Kansas mast influence on measurements of stressand cup anemometer overspeeding', Bound. Layer Meteor., 22, 245-250.
Wyngaard, J. C., Cote, O. R., 1971: The budgets of turbulent kinetic energy
and temperature variance in the atmospheric surface layer, J. Atmos. Sci.,
28, 190-201.
山口浩司, 1999: クローズドバス法とオープンバス法を用いた森林上でのCO2・水蒸気フラックスの評価, 大阪府立大学農学生命科学研究科修士
論文.
Yamamoto, G., 1959: Theory of turbulent transfer in non-neutral
conditions, J. Meteor. Soc. Japan., 37, 60-70.
Zilitinkevich, S. S., Mammarella, I., Baklanov, A. A, Joffre, S. M, 2008: The
effect of stratification on the aerodynamic roughness length anddisplacement height, Bound. Layer Meteor., 129, 179-190.
Arya, S. P., 2001: Introduction to micrometeorogy 2nd edition, Academic
Press, pp.420.
浅井冨雄, 新田尚, 松野太郎, 2001: 基礎気象学, 朝倉書店, pp.202.

•

•

•

•

•

Holton, R. J., 2004: An introduction to Dynamic Meteorology 4th ed.,Elsevier Academic Press, pp.535.
Kaimal, J. C., Finnigan, J. J., 1994: Atmospheric boundary layer flows: theirstructure and measurement, Oxford University Press, pp.289.
近藤裕昭, 2001: 人間環境の気象学, 朝倉書店, pp.156.
Lee, X., Massman, W., Law, B., 2004: Handbook of Micrometeorology: A
Guide to Surface Flux Measurement and Analysis, Kluwer Academic,
pp.268.
Lumley, J. L., Panofsky, H.A., 1964: The structure of atmosphericturbulence, Interscience Publishers John Wiley ＆ Sons, pp.239.

•

文字信貴, 2003: 植物と微気象, 大阪公立大学共同出版会, pp.148.

•

小倉義光, 1999: 一般気象学, 東京大学出版会, pp.308.

•

巽友正, 1982: 流体力学, 培風館, pp.464.

•

横山長之, 1997: 煙ー大気中における振る舞と姿, 白亜書房, pp.149.

•

•

•

•

•

•

•

•

•

Air Data and 3D Displacement Sensors for a Demanding World, <URL: http://www.spaceagecontrol.com/Main/Home>
C50 Site in Takayama, <URL: http://sateco-archive.green.gifu-u.ac.jp/C50/>
CASES-99, <URL: http://www.cora.nwra.com/cases/CASES-99.html>
Davis, <URL: http://www.davisnet.com/marine/products/marine_product.
asp?pnum=00276>
Eddy Covariance Analyzers - LI-COR Environmental, <URL: http://www.licor.com/env/products/eddy_covariance/>
FE2011 (ISHIDA) <URL: http://www.st.hirosaki-u.ac.jp/~ishida/flux/2001/i
ndex.html>
FIFE Project Page <URL: http://daac.ornl.gov/FIFE/fife.shtml>
FLUX 野郎の会改め、FLUX 野郎＆お嬢 の会, <URL: http://tama.green.gif
u-u.ac.jp/~tama/flux/>
Home of HEIFE DATA, <URL: http://ssrs.dpri.kyoto-u.ac.jp/~heife/>

•

•

•

•

•

•

•

•

I. Tamagawa's Page, <URL: http://tama.cive.gifu-u.ac.jp/~tama/COE_study_memo/Tower_Photo1/>
伊藤芳樹，林泰一，玉川一郎: 超音波風速TR-61A型プローブ特性試験,<URL: http://atmos.cc.okayama-u.ac.jp/flux/TR61A.PDF>
Open Path Gas Analyser for Eddy Covariance <URL: http://www.campbells
ci.co.uk/index.cfm?id=2117>
ソニック, <URL: http://www.u-sonic.co.jp/product/meteorology.html#b>
タワーフラックス観測マニュアル, <URL: http://www.ffpri.affrc.go.jp/labs/f
lux/manual/manual_index_j.html>
Vaisala, <URL: http://www.vaisala.co.jp/jp/meteorology/products/weathe
rradars/Pages/default.aspx>
Welcome to the BOREAS Project, <URL: http://daac.ornl.gov/BOREAS/bhs/BOREAS_Home.html>
Wikipedia, <URL: http://ja.wikipedia.org/wiki/%E3%83%A1%E3%82%A4%
E3%83%B3%E3%83%9A%E3%83%BC%E3%82%B8>