地球環境研究，Vol.18（2016）

立正大学 ･ 熊谷キャンパスにおけるドップラーライダによる水平風の観測
高 咲 良 規＊中 川 清 隆＊＊

吉 﨑 正 憲＊＊
蓜 島 徹 也＊＊＊

渡 来武 井靖＊＊
祐 興＊＊＊

キーワード：ドップラーライダ、強風日、弱風日、大気境界層

１ ．はじめに２ ．DL の観測原理

気象学において風（風向 ・ 風速）は重要な観測項目でDL の用語のうちライダ（LIDAR）は Light Detectionある。その場における直接観測として、回転軸の周りにAnd Ranging の略であり、光源としてレーザ光を使うこカップ状の羽をおいてその回転数から風の強さを評価すとを意味する。大気中にレーザ光を発射すると、光は大るロビンソン風速計や、超音波によるドップラー効果を気中にある微小な物体（主にエアロゾル）にぶつかり散利用した超音波風速計（光田 ・ 水間 1964，岡本 1966）乱されるが、散乱された光の一部はレーザ光の発射場所などがある。一方、遠方場をリモートで観測する方法とを発したところに受信波として戻ってくる（図1）。散乱

して、ドップラーレーダやドップラーライダ（今後 DLさせる物体が視線方向に動いているとき、ドップラー効と略記）などのリモートセンシングの手法がある。

果により送信した周波数とは異なる周波数となる。送信DL に関して、水谷ほか（2014）に詳細なレビューが波の周波数を f、受信波の周波数を f ’、光の速さを c、物ある。そこでは、飛行機からの観測（Bibro et al. 1984）

体の速さ v とすると、受信波の周波数は次式で表される。

や地上からの観測（Fujiwara et al. 2011，後藤 ・ 大塚
2012）を挙げている。DL による地上からの観測では大f′=

c−v
fc+v

(1)

気境界層の対流運動は一つの観測対象であり、対流パター

これから、散乱体がレーザ光発射源から遠ざかっていくンとして、風の鉛直シアが弱いときには四角形から六角とき（v > 0）には f ’ は小さく、レーザ光発射源に近づ形の不定形、鉛直シアが強いときには筋状のロール形といてくるとき（v < 0）には f ’ は大きくなる。
（₁）式かなることが示された（Fujiwara et al. 2011）。また、dust

ら２つの周波数が分かると、c の大きさおよび c >> v をdevil や海洋上で見られる竜巻（waterspout）と思われる

使うと、光を散乱させる物体の流れ（v）は（₂）式のよ100m スケールの小さな渦の報告もある（Fujiwara et al.

うに計算できる。

2011，2014）。このように、DL は100m スケールの細かい大気現象を捉えるには非常に有効な測器である。
本稿では、DL の概略と観測された晴天日における強v=c( f − f ')2f(2)通常エアロゾルはミクロン（10－６m）以下のサイズであ風日と弱風日の事例について報告する。晴天日、強風日、弱風日の定義については後述する。晴天日で強風日の場合には、DL で一般風を求める手法について解説する。
また、晴天日で弱風日の場合には低高度角でみると強風－弱風のペアが複数確認されたので、それについて議論する。

図１

 に流されるエアロゾルによって、送信されるレー風ザ光が散乱されて受信波として戻ってくる様子。
そのとき周波数（あるいは波長）は変化している。
左側の三角形は発信機および受信機を表す。

  立正大学 ・ 地球環境科学部 ・ 大学院生＊＊
 立正大学 ・ 地球環境科学部＊＊＊
 立正大学 ・ 学生＊

49

立正大学・熊谷キャンパスにおけるドップラーライダによる水平風の観測（高咲・吉﨑・渡来・中川・蓜島・武井）

(a）

(b）

図２

図３

DL システムの概要。⒜ １ 式の DL の写真。⒝ DL と制御 PC の模式図。

 熊谷キャンパスおよびその周辺の航空写真。赤枠は⒝の範囲、赤点は DL の設置場⒜所である。⒝観測域内の主要な建物（黒実線）、グラウンド（緑色）、道路（二重実線）。
一号機はほぼ北側、二号機はほぼ南側を観測するように配置してあり、太破線は二つの DL の観測範囲の境界を表す。DL からの水平距離は200m ごとに描いてある。

るから、観測される流れは大気の流れと考えられるので、装置からなる（今後、１号機、２号機と呼ぶ）。二台の測こうして風が計測される。

器は水平方向に−90 − 90度、鉛直方向に0−90度の範囲

DL1台の装置による観測手法として、視線方向のを走査できるので、１号機を（ほぼ）北向き、２号機をPOINT 走査、高度角を固定してスキャナーを方位方向に（ほぼ）南向きに配置することにより全天をカバーするこ回転させる PPI（Plane Position Indicator）走査、方位角とができる。測器は３脚で支えられ、光アンテナ装置のを固定してスキャナーを鉛直方向に回す RHI（Range部分は発信機であるとともに受信機である（図２）。発信Height Indicator）走査、上空の風を推定するには VADする光の周波数は eye-safe の1.55μm 帯であり、降水に（Velocity-Azimuth Display）がある（豊田ほか，2009）。

吸収されるため観測は快晴から曇りのときに制限される。

ある距離が離れた２台の DL を使って同じターゲットを観測範囲は30−600m、75−1500m、150−3000m の３通り

観測する手法としてデュアルドップラー法があり、鉛直が設定することができ、解像度は観測範囲を20分割した流が無視できるようなときには水平風を求めることでき大きさである。

る。

本稿で紹介する観測事例は、観測範囲30−600m、解像度30m で、観測の時間間隔として６分とし、高度角は0３ ．立正大学 DL システムの概要
立正大学の DL システムは、２式の測器とデータ処理50−70度の範囲を仰角10度刻みに観測を行った。図３に DLの観測範囲の地図を示す。DL 大学構内およびその周辺の600m の範囲が観測対象である。以下の時間はすべて

地球環境研究，Vol.18（2016）

図４

図５

 015年⒜ ５ 月31日 ９ 時と⒝ ７ 月26日 ９ 時の地上天気図（気象庁 ・ 日々の天気）。
2□は立正大学 ・ 熊谷キャンパスの位置を表す。

 象観測露場における2015年 ５ 月31日 ６ 時から17時までの⒜風速（実線）・ 風向（丸印）、⒝温度（実気線）・ 相対湿度（破線）、⒞日射量（実線）の時系列。

日本時とする。

ほぼ一日風が弱い（平均1m/s 以下）日とした。本稿では、晴天日で強風日として５月31日、弱風日として７月

4 ．観測データおよび観測事例
DL の近傍の地上観測サイトとして、約400m の距離にある立正大学 ・ 熊谷キャンパス内にある気象観測露場（気温、風向 ・ 風速、相対湿度、日射量の10分値）を利用し26日を選び、それらの事例を詳しく解析した。

５ ．解析結果
５ ．１

強風日の観測－５ 月31日た。DL の測器を設置した高さは海抜高度約100m であ

図４⒜に５月31日９時の地上天気図を示す。太平洋側

り、気象観測露場（海抜高度55m）より45m 高い。また、に前線が停滞しており、関東地方では北西風が卓越して立正大学 ・ 熊谷キャンパスから約5.2km 北東方向にあるいた（図略）。図５に、その日の気象観測露場における風

熊谷地方気象台（気温、風向 ・ 風速、相対湿度の10分値）

向 ・ 風速、気温、相対湿度、日射量の時系列を示す。こも利用した。さらに、周辺の環境場を把握するために、れを見ると、９時以降、風速は2.5～3.5ms－1で風向は北気象庁メソ解析モデル（MSM）の1000hPa 面の解析値を～北西であった。また、最高気温は31.0℃と５月として用いた。MSM の格子間隔は5km×5km であり、時間は高温となり相対湿度は20% を切るなど乾燥していた。

間隔は３時間である。

５月31日10時00分における DL による高度角０度の観晴天日は、気象観測露場における日照時間が10時間以測結果を図６に示す。通常、１台の DL で観測される視上の日とした。観測は2015年５月26日から７月27日まで線方向の風速データから、風向と風速（あるいは東西風連続的に行ったが、この期間晴天日は５月30－31日、６と南北風）の２成分を同時に決めるのはむずかしい。と

月１日、７月10－15日。７月21－22日、７月25日の計12ころが、本事例のような強風のときは、観測全域で同じ日であった。気象観測露場の地上風の観測から、強風日ような風が吹いていると仮定することにより、二つの成はほぼ一日強風（平均2.5m/s 以上）が吹く日、弱風日は分を同時に決めることができる。図６から、DL の東北51

立正大学・熊谷キャンパスにおけるドップラーライダによる水平風の観測（高咲・吉﨑・渡来・中川・蓜島・武井）

図７
図６

2015年 ５ 月31日10時00分における高度角10°で観測された DL 視線方向の風速の水平分布。右下の黒矢印と値は、気象観測露場における同時刻の風向と風速を示す。北西方向にデータがないのはエレベータ棟による遮蔽のためである。

2 015年 ５ 月31日10時の強風日における風向の決定法。の DL で観測された風を示す。まず風の時間変動を見ると、７時から８時にかけて北北西から１～3ms－1の風が吹いていた。時間とともに風向は南風に変わり、12時に東方向に視線成分が符号を変えるのがわかる。これは、は風速は3ms－1以上になった。このような風の変動は、

図７の OP 方向にあたるとする。OP 方向に垂線を引く

気象観測露場における風の変動と一致した。しかし、DLと、その方向が観測全域の風の向きに相当する。しかし、で観測された風の時間変動はどのくらいの水平の広がり北北西か南南東かの風向はまだわからない。ここで一方を持つのか－大学構内規模の現象なのか、より大きいメ向の風向を決定するために、P 点を通る円を描き、P 点ソスケールの現象なのか－に関心がある。それを調べるの両側の円の上の点を Q 点と R 点とする。そこで風向がために、まず熊谷地方気象台における風の変動を調べる

P 点と同じとすると、Q 点と R 点の風ベクトルは円を横と（図８）、７～８時にかけては北北西の風、そのあとは切ることになるが、それぞれの視線成分の符号は逆にな南風となっていて、気象観測露場における風の変動と同る。この符号の正負を見ることにより風向がわかる。こじであった。さらに広域な場を調べるために、1000hPaれから、図６の場合は北北西の風向となる。そして、北面の MSM の６時から15時までの３時間ごとの風と気温北西からの風の強さは、図６の南南東方向の視線成分との場をみた（図11）。MSM では、６時には北関東域の山なる。

岳から発する北西寄りの風があり、９時から12時にかけこうして求めた風は、10時00分には北北西の風向で約ては関東地方全域で風が弱くなるとともに気温が上昇し12ms の強さであった。一方、気象観測露場で観測された。15時になると、相模湾方面から南風が吹いてきて、た同時刻の風向はほぼ同じであったが、風速は3ms で北関東域とその西部の山岳地帯は高温となった。これかあった。DL で求めた風速とは約9ms

ら、DL で捉えた北北西→東→南東風の風の変動は、北－1－1－1の風速差があったが、この違いは高さ方向の風の強さの違いからきたも関東スケールの風の変動であったことが分かる。
次に、図10に見られた風の水平分布について調べた。のと考えられる。

図10⒜と⒝には、北北西－南南東の方向には強風と弱風
５ ．２

弱風日の観測－７ 月26日のペアが並ぶのが見られた。強風域から強風域までの幅

図４⒝に７月26日９時の地上天気図を示す。気圧分布はほぼ100m であった。このパターンは、DL の視線方向をみると本州は気圧傾度が極めて弱く、そのため関東地の方位変動によるものではなく風の強さの変動と考えら方も弱風であった。この日は関東地方では猛暑日であり、れる。その鉛直構造を見るために、８時00分から04分ま熊谷地方気象台では15時50分に38.2℃の最高気温を記録

での間の高度角別の風の分布を図12に示す。高度角10度－1した（図８）。また気象観測露場においても風速は1.5msの図12⒜には、DL の南側の水平距離150～250m 付近にぐらいと弱かった。また最高気温は15時10分で36.5℃を北北西－南南東の走向をもつ強風と弱風のペアが複数見記録した（図９）。

られた。しかし、このようなペアは高度角20度ではぼや図10に、高度角０度の７時から12時までの１時間ごと52けてゆき、高度角40度では見えなくなった。DL は周辺地球環境研究，Vol.18（2016）

図８

 谷地方気象台における2015年 ７ 月26日 ６ 時～17時の⒜風速（実線）・ 風向（丸印）と⒝温度（実線）・
熊
相対湿度（破線）の時系列。

図９

図10

気象観測露場における2015年 ７ 月26日 ６ 時～17時の時系列。ほかは図 ５ と同じ。

高度角 ０ 度で観測された2015年 ７ 月26日⒜ ７ 時00分から⒡12時00分までの １ 時間ごとの DL 視線方向の風速の水平分布。図の上の数値の並びは、例
えば、⒜の0726_07:00:00-00:40_el00は ７ 月26日 ７時00分00秒から00分40秒まで高度角00度であることを意味する。右下の黒矢印と値は、気象観測露場における同時刻の風向と風速を示す。の地表から50m ぐらい高い建物の屋上にあるので、高度

図11

1000hPa 面における2015年 ７ 月26日⒜ ６ 時、⒝ ９時、⒞12時、⒟15時における気温（℃、カラー）と風（ベクトル）の気象庁 MSM 出力。立正大学と熊谷地方気象台の位置はそれぞれ□と○で示す。

いうことになる。

角０度とは約50m の高さの水平面に相当する。DL から地上では弱風といっても一般に上層にゆくと風の強さの水平距離を200m とすると、高度角10度では地上72m、は大きくなる。当日の MSM の北関東域の風の鉛直分布

高度角20度では地上122m の高さになる。これから、強をみると、1000hPa 面では1ms－1以下の北西風、975 hPa

風－弱風のペアの鉛直スケールは100m のオーダーとな

面では約2ms－1の北西風、950hPa 面では約3ms－1の北る。そのため、この風のペアの縦横比は１のオーダーと西風であり、少なくとも地上近傍では風の鉛直シアがあっ53立正大学・熊谷キャンパスにおけるドップラーライダによる水平風の観測（高咲・吉﨑・渡来・中川・蓜島・武井）た（図略）。鉛直シアのある場において卓越する対流は鉛al.（2011）に報告されているが、彼らの場合は風の鉛直直シアの方向に沿う２次元ロール解であることを、Asai

シアが強く対流の循環のサイズも数 km と対流混合層と（1970，1972）は理論的に示した。冬季日本海上に発生すして十分発達したものだった。それに対して、今回の事る筋状雲がその例である。DL で見られた風の強弱の分例は同じ２次元ロール状の対流であるけれども、100m ス

布は対流によってもたらされたものであり、地上付近でケールの縦横比１の対流であり対流混合層としては発生鉛直シアを持つ一般風が吹くとき、対流の上昇流域では初期であったと考えられる。より詳細なそうした実態や下から弱い水平風、下降流域では上から強い水平風が移その後の発展等に関しては、数値モデルの適用を含めて流されて形成されたと考えられる。

今後の課題である。

図13に風の立体構造に関する模式図を示す。このような２次元ロール状の対流パターンはすでに Fujiwara et

６ ．まとめ2015年３月に立正大学 ・ 熊谷キャンパスの学生寮屋上に DL システムを設置し観測を開始した。本稿では、晴天日で強風日である場合と晴天日で弱風日である場合に焦点をあて、それぞれ１事例を解析した。
５月31日の強風日の場合には、観測域では同じような風が吹くという仮定を用いることにより、水平風の２成分（風向と風速）を求めることができた。
７月26日の弱風日の場合には、DL では７時から８時にかけては北北西からの約１～3ms－1の風が吹いていたが、時間とともに南風に変わってゆくのが見られた。こうした時間変動は、気象観測露場や熊谷地方気象台の観測データにも見られ、気象庁メソ気象モデル（MSM）の1000hPa 面の風の空間分布から北関東スケールで起きていたことが分かった。 また、早朝から９時ぐらいまで北図12 015年 ７ 月26日 ８ 時00分から04分にかけ
2ての高度角⒜10度、⒝20度、⒞30度、⒟40度における DL 視線方向の風速の水平分布。

北西－南南東の走向に強風と弱風のペアが複数並ぶのが観測された。強風と強風の間の幅は約100m であり、またその鉛直スケールも100m あり、風のペアの縦横幅は１のオーダーであった。こうした風のペアを作ったのは、地上から約100m の高さの大気境界層内における対流であったと考えられる。
謝辞

DL の購入 ・ 運用に関して、松井秀郎 ・ 地球環境科学部長、
川野良信 ・ 前環境システム学科主任、河野忠 ・ 環境システム学科主任、地球環境科学部事務室の岩崎秀男氏をはじめ、環境システム学科の先生方にはお世話になりました。心より感謝します。

図13 2 015年 ７ 月26日 ８ 時において観測された風の概念
図。 ２ 次元ロールは対流の循環を表し、水平に一様に並んでいる。面 A は高さ約50m の建物の DLを通る面であり、面 B は面 A から200m ほど離れている。赤点線は DL からの高度角 ０ 度の線を表す。面 B の白抜きの矢印は高度角 ０ 度で見たときの水平風を表す。
54

参考文献
Asai, T.，1970: Three-dimensional features of thermal convection in a plane Couette flow. J. Meteor. Soc. Japan, 48，
18-29.
Asai, T.，1972: Thermal instability of a shear flow turning

地球環境研究，Vol.18（2016）
the direction with height. J. Meteor. Soc. Japan, 50，525532.の開発と性能評価，高層気象台彙報，70，23－30
光田寧，水間満郎， 1964：超音波風速計とその試作．天気，

Bilbro, J. W.，G. Fichtl, D. Fitzjarrald, M. Krause and R. Lee,
1984: Airborne Doppler lidar wind field measurements.
Bull. Amer. Meteor. Soc.，65，348-359.
Fujiwara, C.，and Y. Fujiyoshi, 2014: Detection of“invisiblewaterspout”using 3D scanning Doppler lidar. SOLA，10，
127-130.

11，33－40.
水谷耕平 ・ 石井昌憲 ・ 森一正 ・ 藤吉康志 ・ 藤原忠誠，2014：
高層気象観測の発展と現状．気象研究ノート（水野量 ・ 上窪哲郎 ・ 定村努編集），229，83－122．
岡本雅典， 1966：超音波風速 ・ 温度計による測定された乱れの顕熱フラックス．天気，13，367－368．

Fujiwara, C.，K. Yamashita, M. Nakanishi and Y. Fujiyoshi,

豊田康嗣，中屋耕，橋本篤，松宮央登，田中伸和，2009 ； ドッ2011: Dust devil-like vortices in an urban area detected by

プラーライダを用いた風速観測手法についての検討，研究

a 3-D scanning Doppler lidar. J. Appl. Meteor. Climatol-

報告 N08032，電力中央研究所報告，32．

ogy，50，534-547.
岩渕真海 2012：ドップラーライダー観測データ品質管理手法

Doppler Lidar Observation of Horizontal Winds at Fine Weather Days atKumagaya Campus of Rissho University
TAKASAKI Yoshinori ＊ , YOSHIZAKI Masanori ＊＊ , WATARAI Yasushi ＊＊ ,NAKAGAWA Kiyotaka ＊＊ , HAISHIMA Tetsuya ＊＊＊ , and TAKEI Yu-ki ＊＊＊＊Graduate Student, Faculty of Geo-Environmental Sciences, Rissho University＊＊
Faculty of Geo-Environmental Sciences, Rissho University＊＊＊
Student, Rissho University

Key words:	Doppler lidar，Strong wind day，Weak wind day，Atmospheric boundary layer

55