第６章　波の平均場への作用 そして運動のエネルギーと、、　　質も形も進度も位置も時間も みな因縁が…

Presentation on theme: "第６章　波の平均場への作用 そして運動のエネルギーと、、　　質も形も進度も位置も時間も みな因縁が…"— Presentation transcript:

1 第６章　波の平均場への作用 そして運動のエネルギーと、、　　質も形も進度も位置も時間も みな因縁が… これまで波は線形で一般風（平均流）を基本の状態とした。 ー＞重力波の作用（働き）の話し。線形の波が、基本の流れを変形する（作用を及ぼす）。それが大気中で興味ある現象を引き起こす。 ６−１：Eliassen-Palmの定理（１） 上の名前の定理を述べる。それは波に伴うエネルギー・フラックスと運動量フラックスとの関係です。 前章で述べたように東西方向の線形の運動方程式は（x-z２次元）、 (1) ここで 　定常な波（一定の位相速度をもち、波の振幅は変化しない）を考え、以下の波形のような、きれいな波を仮定する　 ＜ー 現実はきれいな波は少ないと思う： ここら すると(1)式は (2) 鉛直方向の運動方程式は静力学平衡の近似 (3)

2 連続の式として (4) 熱力学の方程式は(3)を用いて (5) 以上が定常な波にたいする線形波動方程式である（散逸などはない）。 次に(2)式を以下のように変形する。 (6) 図からわかるように、u’w’や p’w’は正の相間関係 ここで u0 は高さのみの関数として偏微分を全微分に置き換えた。この式に左辺の第一項のx の偏微分の中の変数を左から掛けると (7)  線形の式を変形して波の２次の量を評価すること。 　波の２次の量のうち、これからはおもに東西に平均した量を議論する。そこで上式に、東西方向に１波長（ Lx = 2 π / k ）の平均操作を適用してみる。式で書けば そして１波長の平均操作を over bar で表す。 　波の一次の量の１波長の平均はゼロになるが積の１波長の平均は一般的にはゼロにはならない。例えばcos k x の１波長平均はゼロだが、cos 2 kx = (1+cos2kx)/2 の１波長平均はゼロではない。

3 なので、(7)式に平均操作をすることにより以下の式が導かれる。 (8)
一般に なので、(7)式に平均操作をすることにより以下の式が導かれる。 (8) Lindzen（1990）によると、これが Eliassen-Palm の第一定理と呼ばれる（波に伴うエネルギー・フラックス（左辺）と運動量フラックス（右辺のu とw の相関、uという運動量が鉛直に流れるとして）の関係を示すものを第一定理としている、１章に出ていた項）。ただし、非粘性の線形定常波で位相速度がはっきりした波についての関係式。 GFDLの大循環モデルで得られた重力波に伴う運動量フラックスの緯度−高度断面図を示す。成層圏中緯度の東風のところで運動量フラックスは正、西風のところで負になっている（以下で大事）。赤道の成層圏では正になっている。　ただし、これは全ての擾乱成分であり非定常部分も含む。 西風 東風 正 負 北半球冬 左が平均東西風で右がそのときの運動量フラックス。Miyahara et al. (1986)より。

4 ６−２：波のエネルギー方程式について（ここでの式はRossby波も含んでいる）
　波のエネルギー方程式を導く。前節と同様に定常の波を仮定する。ここでは南北方向も考慮する。普通のエネルギー方程式の導出と同様に東西方向の式にu’ を掛け、南北方向の式にv’をかけて足すと（南北シアーもあり）、 (9) 次に熱力学の方程式から（南北基本温度差も含める） (10) この式に 　　を掛けて (11) (9)と(11)を足して (12) 連続の式(13) を使って

5 (14) 静力学平衡を仮定しているので鉛直成分の運動エネルギーが、また音波を落としているので弾性エネルギーが(14)にはない。 運動エネルギー： Potential エネルギー： 波が時間的に成長するような成分を持ち、東西に平均したエネルギーの式は以下のように書かれる。 第２項の微分の中は圧力と南北方向の速度の積であるが圧力によってなされる仕事を示している（ランダウの流体力学６節参照）。第３項の微分の中は圧力によってなされる鉛直方向の仕事を示している。またエネルギーフラックスとよばれる。最後の３つの項は、ここでは基本場から（または基本場へ）のエネルギー変換を表している。

6 いくつかの解析例： 重力波に伴うPotential Energy、 の図（Tsuda et al., 2000, J. G. R. )。５月から８月の平均で、高度は20-30kmの領域での全球分布である。赤いところが重力波のPotential Energy の高いところ。Global Positioning System データから得られたもの。 5-8月のOLR図：対流の強さの指標である。大西洋のPotential Energyの大きいところは強い対流とずれている。 GCMを用いた実験での短周期重力波に伴うPEの図：大西洋に重力波にともなうPotential Energyの大きなところがある（20-30kmの高度）。これは６月の結果

7 衛星から見積もられた中層大気全体の重力波シグナル（microwave limb sounder）
Wu and Waters, 1996, GRL １月 ７月 １月 80km 48km ７月 33km 北半球 但し、重力波に伴う温度偏差（K2）、重力波のpotential エネルギーに対応したもの 緯度高度断面図

8 衛星から見積もられた重力波シグナル（CRISTAという測器）
Ern et al., 2004, JGR 97年8月、25km高度における、重力波に伴う温度偏差（K2）ー＞　慣性重力波として鉛直運動量フラックスの絶対値を見積もっている

9 ゾンデ観測： 18-25kmでの重力波の全エネルギーの時間的変化、Vincent and Alexander, JGR, 2000, 場所はCoros Islands (12S, 97E)、インド洋でのラジオゾンデ観測を解析、1月はwet season 重力波の水平波長頻度 実線：見積もられた月平均 kmの高度 月平均の東西風

10 波数 数値モデルから： 対流からの運動量フラックスとエネルギーフラックス：Eitzen and Randall, JAS, 2005 周波数 u’w’のスペクトル, 14km高度、redが正の 成層圏では上向きのエネルギーフラックス T=4hのおける対流（対流圏）と重力波（成層圏）のようす：細い線は等温位線、太い線は雲>0.1g/kg 鉛直エネルギーの流れ（左図に対応しては、実線の方）

11 ６−３：Eliassen-Palmの定理（2）
条件としては、 　（i）波が定常であること 　（ii）Forcing （例えば thermal forcing ） 　　または Damping がない 　（iii）critical level ( u0 - c = 0 ) がない を満たすときである。 　次にEliassen-Palmの第二定理を述べる。東西／鉛直２次元のエネルギーの式に東西に１波長平均の操作を施す。するとx の偏微分の項は消える。残りを書き表すと、 (15) 補足注：線形波動として物理量が 前の(8)式は 密度を掛けて のように表されているとする。ここでAとBは複素数とする。このとき積の量の平均値のみを問題にするときには、 上式をz微分すると、 z u0(z) で、(15)式を用いると、 なので、下記の条件をみたすときEliassen-Palmの第二定理が導かれる。 (16) ＊じるしはcomplex conjugateを示す。例えばランダウの電磁気学の45節参照

12 u0 ＊critical level ( u0 - c = 0 ) がないWKB近似解の場合: 　WKB近似解が(16)式を満たすことを示しておく。もう一度書き下すと（この波は基本流に対して東に進む波である）、 c 　話しの簡単化のために連続の式としてBoussinesq 近似の連続の式を使い、波は鉛直に平面波的とすれば、 z 最終的に (17) u0(z) このようでも運動量fluxは一定である。 ここで ρ00 は地表面での密度である。Eliassen-Palm の定理が導かれた。 WKB近似解の w’ の中の分母に m1/2 の factor があったが、物理的には波の運動量の保存則（Eliassen-Palm の定理）を満たすように摂動の変動が基本流の中でおこっているといっていいであろう。

13 となる。critical levelの上下で差があることに注意 ー＞ 波が吸収された分の差であることになる。
　前章において重力波の critical level の議論をした。その結果を用いて運動量フラックスのとびを計算する（上の条件（iii）の破れの場合である）。 z>0で（上）、 上 ＋ ＋ ー z<0では（下）、 となる。critical levelの上下で差があることに注意　ー＞　波が吸収された分の差であることになる。 この差により、平均東西風への加速がおこることになる。

14 の別の見方： 線型近似で、鉛直変位と鉛直流との関係は 東西に波の形を仮定すれば， となる．これに圧力偏差をかけて１波長の平均をとると， となる．これを式 　　　　　　　を使って変形すると， のように表される． ＞０の状況 図は位相速度c＞０で山を動かしている状況 斜の矢羽根は風速を示しており、　　　＞０ 　 ＞０　 ＞０のところでは　　　　　　 の 式から負になっており（u0=0）、 shadeの部分を山の所までもっていったところが 　　　　　 に対応している。 　　　　　　　　　の形から　　　　　　　　　で c＞０なので、u>0 のところは 　　　　　　　から、 図のよう に p＞０ のようになっている。 図の影の部分 は 　　 　のところで山がおしていて圧力偏差＞０の状況となっており、ζをhと見なせば山が流体に加える力（圧力の次元）と見做す事が出来る。その力が波動として上方に伝わる。その力の鉛直差で流体が加速される

15 ３次元への拡張 ２次元内部重力波について Eliassen-Palm の定理を導いた。この定理は３次元 stationary（c=0を議論してあるがｃがあっても同様）の長波（ｆも含む） についても拡張されている( Eliassen and Palm,1961) 。 論文の孫引きですが（ p-座標 ）： 基本状態としては 温度風は stationaryな波の式は 慣性重力波に伴う南北、鉛直energy flux：Kawatani et al., GRL, 2003 は安定度 このとき、第１定理は（エネルギーフラックスと運動量フラックスの関係）ー＞　 南北方向 鉛直方向 保存則（第２定理）は がなりたつ。　　の中が Eliassen-Palm fluxと呼ばれる（南北の成分および熱フラックスを含む）。定常で保存的、および臨界層のない波の場合はEliassen-Palm fluxの発散＝０となることがEliassen-Palmの定理である。

16 ６−４：平均東西風（帯状流）の変化について １章において、東西に平均した物理量と波動成分に分離して (18) 東西方向の運動方程式から以下のような式を導いた。 (19) Eliassen-Palmの定理を導いたときと同じように東西と高度のみの２次元運動だけを考えると、左辺の子午面循環やコリオリ項、および右辺の１項は落ちて、 (20) この式により、もし右辺がゼロでなければ東西平均流が変化していくことを示している。はじめ線形の波動方程式を議論していたときは０次の基本場と仮定して線形の波動擾乱を議論していたわけであるが、今や線形の波により基本場が変化していくことがわかる。 前に述べた、Eliassen-Palmの定理は、非常に特別な場合（定常な波で散逸などがない）に(20)の右辺がゼロになることを示している

17 Eliassen-Palm の定理の一般化（破綻したときはどのようになる？）: Andrews and McIntyre (1976, J
Eliassen-Palm の定理の一般化（破綻したときはどのようになる？）: Andrews and McIntyre (1976, J. Atmos. Sci. ) から 擾乱の式として（ブシネスク流体近似、β平面、静力学平衡）、 東西平均流の式は以下のように書かれる。 ここで、 右辺：unspecified forcing terms この方程式には、重力波だけでなく、Rossby波動も含まれる。 ー＞この式を変形すること： （南北熱フラックスが子午面循環をつくり、それにコリオリ力がかかり運動量を変化させるので、そこを一緒に解釈する方法） 平均東西風は緯度、および高さの関数であり、対応した温度場も緯度、高度依存性をもつ。

18 Eliassen-Palm fluxを見て想像されるように、南北熱フラックスがf（コリオリ項）をとおして運動量フラックスとからむので、その項を東西風の変化の式にくりこむと、以下の式になる。 Eliassen-Palm flux をpseudo-運動量フラックスと呼ぶこともある。ー＞それの収束が東西風の変化に対応。 の項（Eliassen-Palm flux divergence）の変形から、平均東西流の加速として近似的に以下の式が導かれている。平均東西流の式の*のついた項は、小さい近似である（定常で散逸やcritical levelがないときはゼロになる）。 ここで、　　　は擾乱に伴う南北変位をあらわし、　　　　 で定義される。 波に対しての外力（１項や２項）、transienceの時（３項）、critical level （２、３項）のところで東西風が変化することを示している。  Eliassen-Palmの定理がなりたたない状況で東西風が変化していく。 Richardson数が大きく、赤道β平面で、cをもつtrap赤道波動について：

19 中層大気における平均東西風の変動について（言葉の羅列ですが）
１月の平均東西風 西風 赤道域半年振動 中間圏弱風層 東風 赤道 突然昇温、惑星波動による（９章） 重力波も関係 赤道域下部成層圏準２年振動

20 熱圏下部の平均東西風に大気潮汐波（全球的な重力波）が寄与をしている話がある
大気中の１日潮汐波（お日さまの加熱と一緒に西に伝播）の非線形効果を計算した結果。赤道域は東風になっている。観測の東風と対応している。一方中緯度では西風が生成されている。 平均東西風の緯度ー高度断面図において、下部熱圏に注意してほしい。赤道域で東風が吹いている。 Miyahara, 1978, J. M. S. J.

21 金星大気の成層圏における高速の平均東西流に波動による運動量輸送が重要な役割を果たしている。
CCSR/NIES/FRCGC AGCMで得られた高速風の実験結果である（Ikeda et al.,2006)。金星大気の成層圏では高速の風が再現されている。 温度観測から推測された平均東西風 しかしながら、この結果では対流圏の中では東西風は再現されていない。 観測されている東西風 赤道 45ºN 東西風 0

22 金星大気の成層圏における高速の平均東西流に寄与する波動の輸送：

23 運動量バランスには大気潮汐波動が重要である、 今の場合は、潮汐が生成されている場所で高速風が生成されている。 地球大気の下部熱圏の東風は潮汐波動がつぶれている場所で作られている。（80km近傍の減速と対応している） 赤道での各項の加速