地球環境研究，Vol.19（2017）「那須おろし」発生時のスコラー数の特徴
青

木

翔

太＊

渡

来

靖＊＊

キーワード：那須おろし、スコラー数、局地風

１ ．はじめにの橋梁で脱線事故が那須おろしによって引き起こされている。近年、鉄道や高速道路等の交通インフラでは那須

東日本の太平洋側各地では、主として西高東低の総観

おろしに対する対策が進み重大な事故は報告されていな場の下でおろし風と呼ばれる局地風が発生する。おろし

いが、那須おろしが原因と考えられる樹木の倒壊などの風は気象条件や周辺地形の条件によって猛烈な風をもた中小の被害は近年でもたびたび報告されている。例えばらすことがあり、おろし風によって過去に多くの重大な那須塩原市において、2008年２月23日～24日午前中にか事故が引き起こされてきた。例えば松井（2012）では、

けて発生した突風によるイチゴハウス、花ハウス等のビ

琵琶湖西岸に吹走する比良おろしが原因となった事故とニールの剥離、立木等の倒壊による家、納屋等の被害がして、1997年６月28日の比良駅構内での貨物転覆事故を報告された（栃木県農業共済組合連合会 2008）。このよ挙げている。栃木県北部の那須野ヶ原においても、冬季

うに、那須おろしは地域住民の生活に大きな影響を及ぼに「那須おろし」と呼ばれる北ないし北西の強い局地風しており、那須おろしに伴う強風発生のより正確な予測が発達する（小園 1983）。那須野ヶ原地域は大佐飛山地、のために、その発生条件について調査することは重要で那須連山の南東側に位置し、箒川、那珂川に挟まれた扇ある。しかしながら、那須おろしの発生条件については、状地である（図１）。那須地方の中でも特に風の強い地域

これまでに調査は行われていない。

となっており、1950年１月には那須野ヶ原を流れる箒川

青木（2016MS）は、2004/05年～2013/14年の10年間の寒候期（11月～４月）を対象に、農研機構畜産草地研究
所那須拠点（以下、NILGS：図１●）で観測された地上
風向風速データを基にして那須おろし（風向が西北西～
北北西、風速が５ms 以上）の積算吹走時間の経月変化－1を求め、那須おろしの出現頻度の季節変化を明らかにした。図２は、青木（2016MS）と同様の手法で2014/15年
までの11年間について調査した結果である。これによると、那須おろし吹走時間は３月が最も長く、11年間で積算600時間を超えており、この特徴は青木（2016MS）と一致している。一方で、関東平野及びその周辺地域においておろし風が発生しやすいといわれる西高東低の気圧配置の出現頻度（図３：吉野 ・ 福岡1967を参考に2004/05
年～2014/15年の期間を対象に作成）は１月にピーク（41.9％）があり、３月は15.8％と比較的少ない。つまり、
那須おろしの吹走しやすい条件として、西高東低の気圧
図１

研 究対象領域。等高線は200ｍ間隔。●は NILGS、■は黒磯 AMeDAS、◆は那須高原 AMeDAS を意味する。▲はそれぞれの山頂を示す

＊		 立正大学大学院地球環境科学研究科大学院生＊＊
立正大学地球環境科学部

配置とは別の因子が存在する可能性がある。
おろし風のような山越え気流の振る舞いを議論する際に、スコラー数と呼ばれるパラメータがしばしば用いら103

「那須おろし」発生時のスコラー数の特徴（青木・渡来）

那須地方における北西寄りの強い風の発生条件が単に気

吹走時間

1000
900
800
700
600
500
400
300
200
100
0

流の山越え条件によって決まるのであれば、強風の出現
頻度はスコラー数によって説明できるはずである。そこ
で本稿では、那須おろしと呼ばれる強い北西風が吹走する際のスコラー数の特徴について調査し、那須おろしの発生条件や水平分布、経月変化がスコラー数にどのように対応しているかを考察することを目的とする。
11

12

1

2

3

4

月

図２

 004/05年～2014/15年の11年間の那須おろしの月
2
積算吹走時間の経月変化（青木 2016MS を参考に作成）

那須連山が存在しており、それぞれの山の間に鞍部が存50.0%在している。また、北西から南東方向へ、主に箒川、蛇

出現率

40.0%

尾川、那珂川の３つの川が流れている。解析対象期間は

30.0%

2004/05年～2014/15年の11年間の寒候期（11月～４月）

20.0%とする。対象領域内におけるおろし風吹走の局地性を検

10.0%

図３

研究対象地域は栃木県北部の那須地方（図１）である。
この地域には南西から北東方向に高原山、大佐飛山地、

60.0%

0.0%

２ ．研究手法

討するために、NILGS（図１●）の地上気象観測データ
11

12

1

月

2

3

に 加 え て、 黒 磯（図 １ ■）、 那 須 高 原（図 １ ◆） の4

 004/05年～2014/15年における西高東低型気圧配2
置出現率の経月変化（気象庁 ・ 日々の天気図を元に作成）

AMeDAS データも用いた。解析した気象要素は風向風速であり、時間間隔は１時間である。本稿における那須
おろしの定義は、青木（2016MS）や青木 ・ 渡来（2016）と同様に、各地点における風向が西北西～北北西、風速れる。一般場の成層が一様であるとすると、スコラー数
l（m ）は以下の式⑴で求められる。

が５ms 以上のときとした。
－1

スコラー数の算出には、風上の日本海側に位置する気－1

象庁秋田地方気象台の高層気象データ（925hPa、900hPa、⑴

850hPa、800hPa、700hPa 各指定気圧面の気温、風向風
速、ジオポテンシャル高度）を使用した。このデータの̅ は平均温位（K）
ここで、ɡ は重力加速度（ms ）、θ、U

時間分解能は、毎日日本時間09時と21時の12時間間隔で－2

は風速（ms ）、z は高度（ｍ）、N はブラント－ヴァイ

ある。本稿が用いるスコラー数は700hPa 面以下（地上～

サラ振動数（s ）であり、各変数は山の風上側平野部の3000ｍ高度）の平均的な値とすることとし、以下の式⑵値で見積もられる。スコラー数は気流が障害物（山）をのように求めた。

－1－1

乗り越え易いか否かの指標であり、これを用いることによっておろし風が発達するか否かを判別することができ⑵

る。スコラー数が山岳地形の波数よりも小さいときには、鉛直断面上の流線の波形は山頂の両側で対称となり、上

ここで、添え字は算出に用いた気圧面の番号を意味し、方にいくほど波が減衰する。一方でスコラー数が山岳地

添え字が１の時は925hPa 面、添え字が５のときは700hPa

形の波数よりも大きいときには、山頂の両側で流線の波面を表す。Z はジオポテンシャル高度である。U には各形が非対称となり、波面が上空ほど上流側に傾く。後者

気圧面での風の北西風成分の大きさを用い、温位θは各の場合、風下側斜面に沿って流線が密となり、おろし風

気圧面の気温から求めた。なお、那須おろし吹走域の風が発生するとされる（小倉 1997）。

上に位置する高原山－大佐飛山地－那須連山はほぼ南西

式⑴から分かるように、スコラー数は風速や安定度の－北東方向に連なっており、主稜線に垂直な風が卓越す関数であるため、西高東低の気圧配置となっても常におる事例のみを対象とするため、秋田の850hPa 面での風向

ろし風の発生しやすい状態であるとは限らない。また、が西～北であった事例のみを抽出して解析した。

104

地球環境研究，Vol.19（2017）

３ ．結
3.1

51.3 ％ を 示 し て お り、 経 月 変 化 の 特 徴 は 前 報（青 木果2016MS）の結果と同様である。一方、那須高原と黒磯で

那須おろしの出現頻度の経月変化

は１月に最大値を示しており、その出現率は那須高原が図２に示した如く、前報（青木 2016MS）では箒川谷

57.0％、黒磯が40.6％である。また、黒磯における那須お口の南東に位置する NILGS において那須おろしの月積算

ろしの出現率は，全期間を通して、他の２地点より小さ

吹走時間の経月変化を調べ、那須おろしが３月に最も吹く推移している。このように、那須おろし吹走域内におきやすいことを明らかにした。この経月変化の特徴が那いて那須おろしの出現率やその経月変化には地域差が存

須おろし吹走域の他の地点でも同様に見られるのかを調在することが明らかとなった。

査するために、NILGS に加えて黒磯、那須高原の計３地点における那須おろし出現率の経月変化を調査した（図3.2

那須おろし吹走時に出現するスコラー数のヒストグ
ラムの地域差４）。図１から分かるように、黒磯■は NILGS ●から北東に約10km 離れており、大佐飛山地の南東に位置する。

図５は地点別に求めた解析対象期間内の北西風吹走時

那須高原◆は NILGS ●から北北東に約25km、黒磯■かにおけるスコラー数のヒストグラムである。階級幅0.3×

ら北に約16km 離れており、那須連山の南東側斜面に位10 m で北西風が測定された時刻の直近のスコラー数の置する。図４を見ると、NILGS は３月に出現率が最大のヒストグラムを作成した。直近６時間（前後12時間）以－3－1

内に１度でも北西風が測定されれば前後12時間内の北西
風出現回数とは無関係に度数１として処理した。最低１

60.0%

那須おろし出現率50.0%(a)NILGS

回、最大12回測定されている北西風のうちの最大風速が40.0%

５ ms を上回っている場合の度数は影を付して示してあ30.0%る。北西風が観測された全ての事例数（灰色＋白色）に20.0%

おいて、３地点のスコラー数のヒストグラムに差がある

10.0%

かどうかχ 検定を実施したところ、χ 値は0.824となり、

0.0%

－1

2

有意水準１％で３地点のヒストグラムが異なるとは言え
11

12

1

月

2

3

4

ない。また、どの地点においてもスコラー数が0.6×10－3

那須おろし出現率－3～ 0.9×10 m の範囲となる事例が約500事例（全体の約

60.0%
50.0%

2

－1

30％）で最も多かった。那須おろし発生時に相当する風(b)黒磯

速５ms 以上の北西風出現数でも３地点全てにおいてス
－1

40.0%

コラー数が0.6×10 ～0.9×10 m の範囲となる事例が最30.0%も多かった。しかしながら、北西風が観測された全事例

20.0%

数に対する風速５ ms 以上の北西風事例数の割合は、

10.0%

NILGS で43％、那須高原で42％とほぼ同じであったが、

0.0%

－3

－3

－1

－1

11

12

1

月

2

3

4

黒磯は32％とほかの地点に比べて約10％少なかった。これは図４で見られた那須おろし出現率の経月変化の特徴と一致している。また、風速５ms 以上の北西風におい
－1

60.0%

那須おろし出現率50.0%(c)那須高原てもスコラー数のヒストグラムに３地点間の差があるか否か調べるために、それぞれの２地点のヒストグラム間

40.0%においてχ 検定を行った。黒磯－那須高原間はχ 値が
2

30.0%

25.984であり、有意水準１％でこれらのヒストグラムの20.0%間に差あるとは言えない。一方 NILGS－黒磯間はχ 値が
2

10.0%
0.0%

図４

2

91.018、NILGS－那須高原間はχ 値が175.728であり、有
2

11

12

1

月

2

3

4

 NILGS、⒝ 黒磯、⒞ 那須高原における2004/05⒜
年～2014/15年の那須おろし（風向が西北西～北北－1
西、風速が ５ ms 以上）の出現率の経月変化

意水準１％でこれらのヒストグラムの間には差があると言える。よって、これら３地点のうち、NILGS は他の２地点と比較して、那須おろし吹走時のスコラー数の特徴が異なることが推測される。
105

「那須おろし」発生時のスコラー数の特徴（青木・渡来）

山越え気流に関しては、理想化された条件下での数値

大佐飛山地、那須連山をそれぞれ孤立峰として、それら実験等に基づいて、風速や大気の成層状態、地形等の条が南西－北東方向に連座する地形であると見做せる。さ件で流れの状態が複数のレジームに分類されることが詳らに、那須おろしの主風向は北西であり山塊の連座する細に調査されている。３次元のベル型孤立峰における山方角と直交するので、北西風に対しては単一孤立峰と比越え気流の振る舞いについては、Smith（1980）の線形モ

較的同等な地形条件であると考えられる。そこで本稿でデルに基づく研究や、それを非線形モデルに拡張したは、那須おろしをもたらす山の地形が上述のような三連

Smolarkiewicz and Rotunno（1989）の研究などで明らかの孤立峰と見做せるものと仮定し、Smolarkiewicz andにされている。しかしながら、理想化実験に基づくこれ

Rotunno（1989）の研究成果に基づき山越え気流の振る

らの結果を現実場にそのまま適用することは困難な場合

舞いについて検討する。

が多い。本稿においても、図１から分かるように、那須第１章で述べたように、気流が山を越えておろし風がおろし吹走域の風上にそびえる山々を単一の孤立峰と見発達するかどうかの閾値は山岳地形の波数である。Smith

做すことは困難である。しかし単純化すれば、高原山、（1980）や Smolarkiewicz and Rotunno（1989）では、座標（x,y）における標高 zh が以下の式⑶で表されるベル型
孤立峰を設定した。

600(a)NILGS

500

⑶

度数

400
300
200

ここで、（x0,y0）は山頂の座標、h は山の高さ、r は山頂

100からの水平距離である。a は山の水平スケールであり、

0

0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4

2.7 3.0 3.3 3.6 3.9 4.2 4.5 4.8

スコラー数(×10
5ms 以上5m/s以上m)

が1795ｍ、大佐飛山地が1908ｍ、那須連山が1917ｍと多

(b)黒磯

500

度数

からわかるように、a は標高 h/（2~_2）の等高線の半径に相当する。本稿への適用を考えると、まず標高は高原山

5ms 未満
5m/s未満600

少異なるが、簡単のために h はこれらを平均して h=1873ｍとする。この時、a は山頂から標高が1873/（2~_2）
≈662

400

ｍとなる地点までの水平距離であるから、地形図から南300

東側斜面の尾根部複数地点を選び距離を測定して平均し200

たところ、約8.0km であった。Smith（1980）によれば、

100

孤立峰によって励起される山岳波の波数は a の逆数に等0
0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4

2.7 3.0 3.3 3.6 3.9 4.2 4.5 4.8

スコラー数(×106005m/s以上5ms 以上)m

5m/s未満5ms 未満(c)那須高原500

度数式⑶において r=a を代入すると zh=h/（2~_2）となることしいから、おろし風が発達するかどうかの閾値に相当す
る波数 k t は、⑷となる。つまり、スコラー数が0.13×10 m 以上のとき、
－3

－1

400

おろし風が発達する。図５をみると、NILGS、黒磯、那300

須高原のどの地点においても、スコラー数0.3×10 m

200

未満で５ms 以上の強風が吹く事例はほとんど無く、そ100もそもこのような低スコラー数で北西風となる事例がほ

0

－3

－1

0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4

2.7 3.0 3.3 3.6 3.9 4.2 4.5 4.8

スコラー数(×10

図５

106

－1

5m/s以上5ms
以上

m )5m/s未満5ms 未満 NILGS、⒝ 黒磯、⒞ 那須高原における北西風⒜
吹走時のスコラー数のヒストグラム。灰色は ５－1
ms 以上の強風（那須おろし）吹走時、白色は、－1
５ ms 未満の弱風吹走時を表す

とんど発生していなかった。
また、山越え気流は、山頂を越えるような鉛直的な流れが卓越する場合と、山を水平方向に迂回するような流れが卓越する場合に大別され、このような流れの性質の違いは、フルード数 Fr

地球環境研究，Vol.19（2017）

500⑸

Rotunno（1989）によれば、三次元ベル型孤立峰の周辺

度数に よ っ て 表 さ れ る（荒 川 1975）。Smolarkiewicz and0.13～1.1

300

1.1～ (×10m

3

4の流れにおいては Fr ≈0.5を境にして、フルード数が大き200

い時には山を鉛直に越える流れが卓越する一方、フルー

100

ド数が小さい時には山を迂回する流れが卓越して山の風

0下側に渦の対が形成されるという、２つのレジームが存在する。また、フルード数が小さい時には山を迂回する

図６

両側面の山麓に、地上風の強風域が出現する。本稿における地形条件に当てはめると、この強風域は隣り合う山

1

月

2

80
70

る。式⑴、⑸より、２つのレジームの境界となる Fr=0.5

60
度数

る。

12

90

おいては鞍部の風下に強風域が出現することが予想されをスコラー数で表現すると、その値 lt は式⑹で求められ

11)

 田上空の850hPa 面高度におけるスコラー数帯
秋
別の西～北風出現度数の経月変化。白色は那須おろしが発生しやすいと考えられるスコラー数の範－3－3－1
囲（0.13×10 ～1.1×10 m ）の西～北風出現度
数を表す

流れと風下渦による流れが収束し、主風向に対する山の塊間の鞍部に位置することから、低フルード数の流れに0.0～0.13

400

50
40
30

⑹

20

青山（2006）は式⑹を用いて、松本盆地や伊那盆地で発

0

10
11

12

生する南風の性質を推定した。本稿においては、式⑹に
h=1873ｍを代入し、lt は以下の式⑺のように計算される。

NILGS

図７

2

月

黒磯

3

4

那須高原

 ILGS（黒）、黒磯（灰）、那須高原（白）の各地
N
点で前後 ６ 時間以内に那須おろしが吹走した時の、
秋田上空850hPa 面高度における西～北以外の風の出現度数の経月変化⑺
つまり、スコラー数が1.1×10 m 以上の時に気流は山
－3

1

－1

を迂回しやすく、1.1×10 m 以下のときに気流は山頂－3－1を越えやすい。すなわち本稿の研究対象領域においては、

ラー数帯別に示したものである。上記で議論した、おろ

スコラー数が0.13×10 ≤ l ≤1.1×10 の範囲にある時、し風が出現しやすいと考えられるスコラー数の範囲（0.13

おろし風が発生しやすいと考えられる。図５によると、

×10 ～1.1×10 m ）は白色で示している。この白色の５ms 以上の北西風事例に対して全体の88～91％が、おスコラー数帯での850hPa 面高度における西～北風出現度－3

ろし風が発生しやすい条件におおよそ相当する0.3×10数のピークは１月にあり、１月全出現度数の約75％に達≤ l ≤1.2×10 の範囲に含まれており、那須おろし発生時している。また11月～２月に限れば、その経月変化は図には上記のスコラー数条件が概ね満たされていることを

４に示される那須高原や黒磯における那須おろし出現率示している。また、NILGS と那須高原での５ms 以上のの経月変化とよく対応している。一方３月、４月では、

北西風の出現度数を比較すると、スコラー数0.3×10 ～白色のスコラー数帯での850hPa 面高度における西～北風

0.6×10 m では那須高原（227事例）の方が NILGS（178

出現度数が２月以前に比べて大きく減少しているのに対事例）よりも多く、0.9×10 ～1.2×10 m では逆にして、那須高原や黒磯での那須おろし出現度数は２月とNILGS（165事例）の方が那須高原（129事例）よりも多大きく変わらない。図７は、秋田の850hPa 面高度での風

くなっており、那須おろし発生時の NILGS におけるスコ

向が西～北以外であったにもかかわらず、前後６時間以

ラー数は他の地点に比べてやや大きいという特徴が見ら

内に那須おろしが観測された事例の出現度数の経月変化れる。

を観測地点別に示したものである。これを見ると、３月、

－3

－3

－1

－3

－1

－3

－3

－1

－3

－3

－1

－3

－3

－1

図６は、秋田上空850hPa 面高度（約1500ｍ）における４月各月の出現度数が70～80に達しており、11～２月各

西～北風の出現度数の経月変化を0.0×10 ～0.13×10 、月の出現度数の約２倍となっている。すなわち、３月、

0.13×10 ～1.1×10 および1.1×10 以上の３つのスコ

４月には図６の範疇外で発生した那須おろしが比較的多

－3

－3

－3

－3

－3

107

「那須おろし」発生時のスコラー数の特徴（青木・渡来）

下側では強風域が平野部に長く伸びている。本稿で地上

200

観測風を解析した３地点を図９に重ねると、NILGS は鞍

度数

150

部の風下、黒磯は山頂部の風下、那須高原は山頂部の風下側斜面に位置している。図９において NILGS と那須高100

原は強風域に位置しているのに対し、黒磯は相対的な弱50風域に位置しており、このような那須おろしの地域差が出現率の地域差をもたらしたと考えられる。

0
北

東

南

西

風向

図８

 ILGS において前後 ６ 時間以内に那須おろしが吹N走した時の、秋田上空850hPa 面高度における西～北以外の各風向間の風のヒストグラム

地上風速の分布に地域差が現れる原因に関連して、青木 ・ 渡来（2016）は、強風域の先端がハイドロリックジャ
ンプの発生と対応しており、山頂部の風下と鞍部の風下
でのハイドロリックジャンプ発生位置の違いが強風分布を決めているとしている。Saito（1992）は、鞍部が存在く存在することを示す。また図８は、図７に収録されたする山脈を越える浅水波の解として、ハイドロリックジャ

那須おろし出現時直近の秋田上空850hPa 面の風向のヒスンプを伴う流れを理論的に求めた。その結果、山脈の風トグラムである。図８は、図６の範疇外で発生した那須

下側に停滞性のジャンプを形成するような条件下におい

おろしは風上側風向が南西～西の場合が多いことを示しては、狭い鞍部の存在がジャンプの形成位置を山越え気ている。これらの事から、３月、４月には気圧の谷の通流に垂直な方向に一様な２次元山岳の場合より風下側に過によって850hPa 面高度の風向が短時間で大きく変化す

移動させ、ジャンプ手前の流速は強められることを示しるような場合が多いため，直近の秋田上空850hPa 面高度ており、これは青木 ・ 渡来（2016）のシミュレーションの風向が西～北以外となったことが示唆される。

結果と整合的である。Saito（1992）の解を現実場に適用すれば、特に鞍部の風下で強化される那須おろしの発生

４ ．考

条件をより明確に出来る可能性が高くなったが、この点

察については稿を改めたい。また、図５によると、鞍部の高原山－大佐飛山地－那須連山のおおよその地形条件

風下に位置する NILGS は、他の地点に比べるとやや大きから式⑷、⑺により求めた、那須おろしの出現し易いス

いスコラー数の時に那須おろしが発生する頻度が高いこ

コラー数の範囲は0.13×10 ≤ l ≤1.1×10 であると見積－3－3もられた。那須おろし吹走域に位置する NILGS、黒磯、那須高原の３地点における北西風の時のスコラー数のヒストグラム（図５）によれば、スコラー数が上記の条件を満たすとき常に強風となるわけではないが、風速が５ms を超える事例の約９割は上記条件を満たしている。
－1

この結果は、那須地方で観測される北西寄りの強風が、山岳による風下側斜面下降流の強化の影響を受けている、
すなわちおろし風であることを意味している。
図４や図５によると、本稿で定義された那須おろしの出現率およびその経月変化は、那須おろし吹走域内で一様ではなく、黒磯においては NILGS や那須高原に比べて10％前後低い。また、この傾向は季節（月）やスコラー数の大小には依存していないように思われる。青木 ・ 渡来（2016）の示した典型的な那須おろし事例の数値シミュ
レーション結果（図９）によると、那須おろし吹走時の山岳風下側の地上風の空間分布は非一様であり、高原山
や大佐飛山地、那須連山の山頂部の風下側を中心に強風域が山麓付近で途切れるのに対し、各山塊間の鞍部の風108

図９ 域気象モデル WRF により再現計算された那須領おろし吹走時（2015/03/10 0400JST）の地上風速分布。図中のシンボルは図 １ と同様である。等値
線は標高（100ｍ間隔）を示す（青木 ・ 渡来 2016に彩色）。左上枠内は気象庁 MSM（メソ数値予報モデル）の2015/03/10 0300JST、0600JST の秋
田上空850hPa 面高度における風向風速の ０ 時間
予報値を表す地球環境研究，Vol.19（2017）
350

変化は地上風の特徴を示しているのに対して、那須高原

300

は上空風の特徴を示す傾向があり、明瞭な日変化を示さ

積算吹走時間250

ないと推測される。日射量は季節変化するので、混合層200の成長は１、２月より３、４月の方がより活発である。

150

この混合層の季節変化が NILGS における那須おろしの季

100

節変化に影響を及ぼしているかもしれない。しかし、こ500の論理では黒磯で１月にピークとなることを十分説明で0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
時間(JST)
NILGS

黒磯

那須高原

図10 那
 須おろし積算吹走時間の日変化の比較。実線は
NILGS、点線は黒磯、破線は那須高原を表す

きず、他の要素が関わっている可能性もある。3.2節で指摘したように、図６作成のためのデータ抽出方法に起因して３、４月を過小とする要素が含まれていた可能性もあり、さらなるデータの精査が必要である。経月変化の非一様性をもたらしうる熱対流混合やその他の要因を探とが明らかになった。上述の如く、スコラー数が大きいるために、数値モデルによるシミュレーションや混合層

ほど山越え気流は迂回流の性質が強くなるので、NILGS

上端の風の観測などのさらなる研究が必要である。

で強風が発生する場合には、山頂を迂回しギャップとなっている鞍部を通過する流れにより強化されたものが相当数含まれていると考えられる。山頂を迂回した気流は鞍部を選択的に流れ、それが鞍部風下における強風帯形成の一因となっている可能性がある。

５ ．結

論

本稿では、寒候期（11～４月）に吹く風速５ms 以上－1の北西寄りの地上風として定義した那須おろしについて、

青木（2016MS）は、鞍部風下の NILGS における観測

その吹走時におけるスコラー数の特徴を調べた。高原山

結果から、那須おろしが３月に最も多く発生するという－大佐飛山地－那須連山のおおよその地形条件から、山経月変化を示した（図２）。しかし、本稿による追加調査

岳の影響によるおろし風が形成されやすい条件は、スコによって、３月にピークとなる特徴は那須おろし吹走域

ラー数 l が 0.13×10 ≤ l ≤1.1×10 の範囲であることが全体での共通した特徴ではなく、山頂部の風下に当たる

導かれた。那須おろし吹走時のスコラー数の約90％が上

黒磯や那須高原では１月がピークとなることが明らかと記の条件を満たしており、那須おろしは山岳による風下なった。図６に示した秋田上空850hPa 面高度（約1500

側斜面下降流の強化の影響を受けたおろし風であることｍ）におけるスコラー数帯別西～北風の出現度数の経月が示された。那須おろしの出現頻度には、鞍部風下の変化によると、スコラー数から予測されるおろし風発生

NILGS で多く山頂部風下の黒磯で少ないという地域差が頻度は１月にピークとなることを表しており、NILGS の見られた。これは青木 ・ 渡来（2016）による那須おろし

結果はむしろ特異である。このような、場所による経月

吹走事例のシミュレーション結果とよく対応しており、変化の相違をもたらした原因の候補の一つとして、日射個々の事例における地上風の非一様性が出現頻度に影響

による熱対流混合が挙げられる。図10は、那須おろしのしていると思われる。また、鞍部風下の NILGS で観測さ積算吹走時間の地点別の日変化である。平野部に位置すれる那須おろしの際には、他の２地点における那須おろ

る NILGS や黒磯では、14時頃をピークとした一つ山の分しの際に比べてスコラー数がやや大きい値となる傾向が

布をしており、那須おろしはその多くが日中に観測され示された。山越え気流は、スコラー数が大きくなるほどていることがわかる。これは、日中は混合層の発達によ

山頂を水平に迂回する流れがより卓越するため、NILGS

り上空の運動量が地上に運ばれやすくなるためである。

で観測された那須おろしの一部では、山頂を迂回し鞍部

一方、那須高原ではこのように明瞭な日変化は見られなを吹き抜けた流れが強風の形成に寄与していた可能性がかった。また、わずかにみられる日変化においても、ピーある。青木（2016MS）で示された、NILGS における那クが他の地点に比べて遅い時間に出現している。これは、

須おろしの出現度数が３月に最大となる特徴は那須おろ

典型的な上空風の特徴である（例えば、Kusaka et al.

し吹走域全体の特徴ではなく、山頂部風下の黒磯や那須

2011）。那須高原は那須連山の風下側斜面に位置するた高原ではむしろ１月に最大となることがわかった。スコ

め、他の２地点に比べて標高が高い（NILGS：335ｍ、黒ラー数から予測される那須おろしの出現度数は１月に最

磯：343ｍ、那須高原：749ｍ）。そのため他の２地点の日大を示しており、NILGS における結果はむしろ特異であ

－3

－3

109

「那須おろし」発生時のスコラー数の特徴（青木・渡来）る。このような経月変化の地域差の原因を理解するため

風系分布．東北地理，35，20-25．

には、さらなる研究が必要である。
謝

小園 修，1983：那須野ヶ原北西部における「那須おろし」のKusaka, H., Y. Miya, and R. Ikeda, 2011: Effects of solar radiation amount and synoptic-scale wind on the local wind

辞本研究を進めるにあたり、11年間にわたる長期間の貴重な気象観測データを提供して頂きました国立研究開発
法人農業 ・ 食品産業技術総合研究機構（農研機構）の中
村美和研究員をはじめとする農研機構畜産草地研究所那須拠点の方々に心より御礼申し上げます。匿名の査読者からは大変詳細かつ有益なコメントを頂きました。記して御礼申し上げます。

“Karakkaze”over the Kanto Plain in Japan. J. Meteor.
Soc. Japan, 89, 327-340.
松井 一幸，2012：比良おろし．内藤 正明（編），琵琶湖ハンドブック，263pp．
小倉 義光，1997：メソ気象の基礎理論．東京大学出版会，
228pp．
Saito, K., 1992 : Shallow water flow having a lee hydraulicjump over a mountain range in a channel of variablewidth. J. Meteor. Soc. Japan, 70, 775-782.
Smith, R. B., 1980 : Linear theory of stratified hydrostatic

参考文献

flow past an isolated mountain. Tellus, 32, 348-364.

青木 翔太，2016MS：領域気象モデルを用いた那須おろしの解析．立正大学平成27年度卒業論文，108pp．
青木 翔太 ・ 渡来 靖，2016：WRF を用いた那須おろしの立体構造の解析．2016年度春季大会講演予稿集，日本気象学会，
109，95．

Smolarkiewicz, P. K., and Rotunno, R., 1989：Low Froudenumber flow past three-dimensional obstacles. Part I:
Baroclinically generated lee vortices. J. Atmos. Sci., 46,1154-1164.
栃木県農業共済組合連合会，2008：NOSAI とちぎホームペー

青山 高義，2006：松本盆地の南風にみられる地形の影響について．地理学評論，79，405-422．

ジ，観察員情報，平成19年4月～．http://www.nosai-tochigi.
or.jp/kans/19/kansatu19.htm （閲覧日：2017/1/22）

荒川 正一，1975：おろし風を中心とした山越え気流論．気象研究ノート，125，51-84．

吉野正敏 ・ 福岡義隆，1967：半旬別の気圧配置ごよみ．天気，
14，250-255．

Characteristics of the Scorer parameter in the generation of the“Nasu-oroshi”
AOKI Shota ＊ and WATARAI Yasushi ＊＊＊Graduate Student, Graduate School of Geo-environmental Science, Rissho University＊＊
Faculty of Geo-environmental Science, Rissho University

Key words: Nasu-oroshi, Scorer parameter, Local wind110