地球での接地境界層の乱流観測から分かったこと
玉川一郎(岐阜大学流域圏科学研究センター）

超音波風速計• 音波を往復させ，その音速の差から媒体である大気等の流速を計測する装置• 超音波パルスを使った実用システムは，光田＆海上電機（現 株式会社ソニック）の開発である。
• 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次元で，乱流変動まで計測可能。航空機に使われるが，地上用に販売していたこともある。Rosemountearospace社 （ピトー管と圧力センサーの間に水滴受けをつけるなどの工夫が必要）

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

大気境界層乱流の話に戻る• 超音波風速温度計の開発以降，接地境界層乱流の直接観測が盛んに行われた。
• 乱流の統計力学としての中心的な理論である
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

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/p21902m3p8106311/)• 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_Obukhov_1954.pdf ) 最初の実験的証拠
渦粘性係数 t  u * (kz) の対数則の拡張を考えた

Monin-Obukhovの相似則
gkz  w' v 'v z


3Lu*

M-Oの安定度

u *3L
gk   w' v ' v 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 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 z/L は• MOの安定度(z/L)は，水平一様定常のもとでの，乱流運動エネルギー方程式の浮力消滅項／シアー生成項
gkz  w' v 'v z
Lu *3kz Uz m  u * zL

植物と微気象（文字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 virtualtemperature)
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.

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

ATIKaijo

長期安定は難しいFlux measurement
Sonic Anemometer Thermometer
Sonic temperature includes cross wind effect andwater vapor effect.

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

In wind tunnel

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 より

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

Ｗ不安定
安定T

Tamagawa 1996

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

不安定
安定U

Tamagawa 1996

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

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

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 とする。
地面の特性量になる。
• 温度に対しても同様にする。地面の温度を放射温度計で計測したもの等に置き換えて処理 （圧力による抵抗，ソースの位置の違い，表面値の定義の相違等でかなり違う値になる）• キャノピーモデルを使う現在でも，キャノピー直上でこの値が使われている。

粗度高ほか上、人間環境の気象学 近藤裕昭著より
数値計算に使用したパラメタ
右、基礎気象学 浅井冨雄、新田尚、松野太郎著
より

ちょっと火星を意識して多分，顕熱フラックスから推定できること

日中の混合層上端の交換率や放射冷却など，さまざまな仮定が必要であるが，温位プロファイルと下端からの熱供給で，混合層高度が推定できるかも。境界層の変化に関する情報
高度

昇温

気圧と放射冷却の推定が必要か
温位

地表面熱収支

粗度高等は，H LE のモデル化に影響Observed data (daytime)
uTv

cw

q

TKC での観測値の例

Observed data (nighttime)
uTv

cw

q

TKC での観測値の例渦相関法による乱流輸送量推定
Eddy covariance:

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

From GAME-Tibet 1998

Covariance w and T

Covariance w and mixing ratio ofwater vapor

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

Unpub. Data by me精度良く計測するには• Ｗ の方向 （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.

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

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

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

結果，顕熱と潜熱のフラックスが測れるようになったので，熱収支観測が盛んに行われるようになった。

Field Projects• 多分， FIFE (The First ISLSCP FieldExperiment ) あたりが最初(1987-1989)。
http://daac.ornl.gov/FIFE/fife.shtml• 他に，TOFA-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でも，多数の地点で乱流観測が行われた。

HEIFESTART OF MY RESEARCH
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の乱流観測の様子

GAME/Tibet
1998 Intensive Observation

安多I belonged to Nagoya univ.

Tanaka et.al. 2001

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

記録• 1990年代以前は，統計値のみを保存。ある
いは 30分前後の期間の乱流データを保存。
• それ以降は，連続的にデータを保存するのが普通。
• サンプリングは，10Hz 程度• 90年代のFlux-PAMは，10分毎に平均値と共分散を記録。

Flux-PAM

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

データ処理の流れ• データを保存 （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 laboratoryin 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
H2Oand CO2 analyzer
LI-7500(LI-COR)

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!

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

Light source & detecter

筺体の熱の影響の補正

Observed data (daytime)
uTv

cw

q

TKC での観測値の例

Observed data (nighttime)uTv

cwq

TKC での観測値の例

Eddy Covariance• Horizontally homogeneous condition
 EC is direct measurement of flux at themeasuring point• Perfectly Steady condition
 Storage between measuring point andsurface can be neglected.

REAL FLUX from surface

Eddy Covariance
AVERAGE:We must choose it to obtainstatistically stable covariance!
Representativeness of the area
We want assume that convective plumes
and other circulations pass through theobservation point

CO2 profile measuring system

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

Radiation observation

Ground observation

So many sensors for profiles

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

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

火星•

6hPa, 180 ± 30K, 自転周期 24時間37分，公転周期 687
日, 重力 0.38 ＊9.8 m/s2 大気成分CO2 （一般気象学 小倉1984より）• 乱流の統計理論には，一応，密度や組成は出てこない。
• 非圧縮仮定が崩れたら未知の世界 （航空力学にある？）• 地上のモニターとしては，地表面状態（粗度）。放射観測と地温観測があれば，熱収支。CO2の潜熱輸送？
• 顕熱から境界層高度？
• 温室効果気体なので，放射冷却も多そう。
• ダストを巻き上げるシアー応力のモニター （地球では，風送ダストのプロジェクトがあった。三上正男ら，風速の3乗でダストが飛ぶ）• そもそも風速，風向のモニター。大気擾乱。日変化，年変化。