Rossby 波動が関係する成層圏平均東西風の変動例

Presentation on theme: "Rossby 波動が関係する成層圏平均東西風の変動例"— Presentation transcript:

1 Rossby 波動が関係する成層圏平均東西風の変動例
1978 冬季 1979 東西平均の温度の時間変化（突然昇温と呼ばれるー＞８章） 東風 西風 東風 西風 １９７９年２月２６日の東西平均風 この現象に中高緯度の惑星波動（planetary波）が重要、β効果から引き起こされる線形波動の振る舞いや性質について述べておく。 北半球冬 １月の平均東西風

2 ５—１：準地衡風方程式について - 線形的なRossby波を取り出すための準備 -
いま中・高緯度を見るため sin 　 をある緯度のまわりに展開（赤道 beta 平面近似と同じやり方） (1) これを以下のように表す　 (2) 第２項の大きさは1000kmのスケールでは　10-11x106m=10-5 1/s 惑星スケールの擾乱について、中・高緯度においては地衡風近似が観測的に成り立つから第０近似の式として　　 (3) (4) これは定常の式である —＞　時間発展的にはどうなるか？　ということで、次のorderに進む 大規模な運動についての鉛直移流の項については w が小さいとしてその項を落とす。すると第１次近似として次の式になる。ここで時間微分、移流項およびβ - 項の u, v については地衡風とし、f 項には次のorderの量を考慮する　　　　　 (5) (6)

3 １の添え字 は第１次近似（地衡風からのずれ、この項があると地衡風が少しずつ運動する）の量。運動が水平的なのでそれを表すものとして渦度の鉛直成分の方程式を導く。渦度の鉛直成分（地衡風成分のみ）は
(7) (6)のｘ-微分から(5)のｙ-部分を引く。地衡風の水平発散はゼロ（(3), (4) から）であることを考慮すれば、 (8) ここで渦度に地衡風近似を使えば（ (3) , (4) ） (9) 第０近似では水平速度は地衡風で非発散であった。１次のオーダーの u1 , v1 に対応した連続の式は (10) である。ここでw にはゼロの添え字をつけた。それを考慮すると(8)の準（完全には地衡風の定常状態ではなく時間変化するのでこの様な名前がついたのだろう）地衡風の渦度方程式は (11) ここで密度（圧力）は H のスケール・ハイトで変化している。

4 次に準地衡風方程式での熱力学の方程式は以下のよう
(12) ここで温度の水平移流の速度の所に地衡風を用いた。 (11)と(12)を一つにまとめると準地衡風近似におけるPotential Vorticity 方程式が得られる。 (13) (14) この様に、準地衡風近似の方程式は１つの変数のみの時間発展の式で表される。また(3)と(4)から流線関数を導入する。 (1５) β はコリオリ項の南北微分からでたことを思い出すと(14)はさらに (16) とすると（準地衡風近似でのポテンシャル・渦度） (17) Rossbyモードの保存的な時間発展の式になる この方程式を、この章では惑星波動の鉛直・南北の伝播性を、７章の一部では中層大気の傾圧不安定や順圧不安定の問題に適用してみよう。

5 ５−２：惑星波（Planetary wave）の鉛直伝播
成層圏大規模擾乱の冬と夏の違い 中・高緯度中層大気の大規模な渦の様子：図は１０hPaの水平断面図（等圧面高度）を示す。地球規模の波動的擾乱（惑星波、Planetary wave 又はロスビー波、Rossby wave ）をみることが出来る。下に対流圏高低気圧の図をのせている。それに比べ、この図では水平スケールが大きいことにきずかれるであろう。なぜ？　これが惑星波の鉛直伝播の問題。 高 図：１９７９年１月２６日、10hPaのHeight図 低 右に夏の場を示しておこう、ほぼ丸いパターン 図：７月平均の10hPa の温度分布　ー＞ 500hPa, 1964 Nov 20の高度場

6 惑星波動（Planetary wave）鉛直伝播の１つの見方
　冬の成層圏では東西方向に一様ではない渦がみられ、夏には東西方向にほぼ一様な風のみである。成層圏平均東西風の緯度−高度断面図では、冬と夏では東西風の向きが異なる。すなわち夏では東風、冬では西風である。冬と夏の、東西方向に一様でない大規模擾乱の振舞いの違いについては、CharneyとDrazinによってはじめて示された(JGR, 1961)。 西風 東風 　線形の議論：基本流として”一定”の東西風のみが吹いていると仮定する。もし東西風が高さの関数であれば、特別の場合を除き解析的に解くのは難しい（高さの１次関数のときは合流型超幾何関数で表される）。一般には東西風は高さ及び南北方向にも変化している。そのときには球面の効果を考慮した準地衡風の方程式が必要であろう。この問題についてはMatusno(1970) を参照　−＞　あとで結果を示す ということで、一定の東西風 u0 が吹いているときの線形準地衡風の方程式は以下のよう、また、ここでN2は一定と仮定、 (18) <-- 東西、南北に波の仮定をおこない、東西に位相速度 c で動くとし、また鉛直方向には密度の効果を考慮して流線関数を以下の様に仮定する。 (19) この様に仮定すると鉛直の構造を決める方程式は以下の様になる。 (20)

7 である。m2 が正ならば波として鉛直に伝播が可能である。m2 が負ならば鉛直方向に指数関数の形になり伝播は不可能。
ここで (21) である。m2 が正ならば波として鉛直に伝播が可能である。m2 が負ならば鉛直方向に指数関数の形になり伝播は不可能。 (21)をみてきづくことは k , l が大きいと負になることである。水平のスケールが小さい（あまり小さくなると準地衡風の近似が使えなくなる、１０００ｋｍくらいの水平スケールまではいいよう）と鉛直方向に伝播しない。対流圏中・高緯度の大規模擾乱の方がスケールが小さい。これは数千ｋｍの擾乱が自励的に対流圏に生成することと矛盾しないだろう。伝播でエネルギーが逃げたら不安定になりにくいであろうから。 定常惑星波の生成として思いつくのは山岳による強制であろう。このとき位相速度はゼロである（海陸の熱の違いによる時も同様）。このとき(21)は (22) 周期 この式から東風のときm2 は負となり鉛直に伝播不可。これが夏の成層圏において定常惑星波がほとんどない理由である。また上式から西風があまり強くても伝わらない。この臨界速度をRossby critical velocityと呼ぶ。式で表すと (23) 波数 陰の部分が鉛直伝播可 ここで β = 1.6 x s-1 m-1 , Lx = 2π / k = km , Ly = 2 π / l =10000km , f = s-1, N = 2x10-2 , H = 7 km と選ぶと Uc = 28 ms-1 、南北巾１００００ｋｍの南北にstanding modeならば Uc = 58ms-1と大きい。

8 これまでは β-平面の一定東西風での議論であった。球面上の準地衡風近似の方程式を使った鉛直、南北伝播の議論をみておこう（cf
これまでは β-平面の一定東西風での議論であった。球面上の準地衡風近似の方程式を使った鉛直、南北伝播の議論をみておこう（cf. Matsuno, 1970）。 球面上で東西一様な東西風（南北と鉛直の関数とする）があるときの線形の準地衡風方程式は（Adrews et al., 1987, Middle Atmosphere Dynamicsから）： 計算に使われた東西風緯度-高度図 有効的なβ 定常な（c=0）惑星波動として（上式で時間微分の項をおとす） のような形を仮定する（東西方向には波形、鉛直南北の波の構造を決めたい） N2が一定のときは以下の式になる、 n : 屈折率と呼ばれる =0は特異点

9 観測で見積もられ、計算に使われた波に関しての下部境界(500hPa)条件、モデル下端で擾乱を強制する
有効βの分布図 振幅 位相 波数１ 緯度 波数0に対する屈折率の二乗 波数３ 定常な惑星波動の臨界層 観測で見積もられ、計算に使われた波に関しての下部境界(500hPa)条件、モデル下端で擾乱を強制する

10 水平分布の計算例（Matsuno,1970, J. Atmos. Sci
水平分布の計算例（Matsuno,1970, J. Atmos. Sci.）：500hPaでの高度場偏差（これを境界条件として与える、100m間隔）ー＞　30kmでの高度分布（波数＝１－３を足してある）、アリューシャン高気圧が見える。右下は観測による高度場。200m間隔で1967年の1月の状況、 境界条件としての500hPa高度場偏差、 左を強制することで、図のような線形波のresponse 60Nでの波数１成分の振幅と位相 観測による、1967年1月

11 波数１と２の鉛直、南北の伝播性、１の方が上層まで伝播している。
左の振幅・位相に対応した波数１惑星波動のエネルギーフラックス、上方および赤道方向に波が伝わっている。 位相 波の振幅 対応：異なる表現 準地衡風近似で、定常惑星波動の場合は、 近年、東風運動量的な流れである Eliassen-Palm フラックスによる表現が多い、惑星波動の場合、

12 波数１、上が波の振幅(decameters)、下が波の位相。左７月（南半球）、右が１月（北半球）
平均東西風、左が南半球で右が北半球、上図が１月で下図が７月 西風 振幅 西風 東風 位相 30hPa以上は Nimbus 5（ Jan. 1973からDec. 1974）, 6（ Jul. 1975からJun. 1978） 衛星データ

13 近年の衛星観測例：Garcia et al., JAS, 2005
赤道中間圏（西風）まで penetration SABERで得られた2002 June-July（夏）の平均温度、右の惑星波動の解析とは季節が異なる 14K 120km 中間圏界面 0.01mb 成層圏界面 1mb 14km 1/25-2/24, 2002：定常惑星波動（波数＝１）に伴う温度の振幅と位相

14 ５−３：球面上の自由振動について 　前節は惑星波動の強制問題（対流圏で生成された波が鉛直に伝播する問題）をあつかった。この節では等温静止大気の自由振動のはなし 地球上の自由振動なので、地表面で鉛直流ゼロが自然な選択であろう　　　　 等価深さ（鉛直構造）が決まる（どんな構造になるか？） 高度座標での境界条件は のように書かれるであろう。 Geopotentialで書き換えると 全微分　　　　は log-p 座標で線形的に（平均の風がないとする） であろうから、 がlog-pでの下部境界条件となる。等温静止大気では なので、上の式は 一方、熱力学の式は log-p 座標では だったので、２つの式から　　　　を消すと、

15 等温大気のとき　　　　　　　　　であるので、上の式を変形して、
となる。上式の時間微分を落とすと、　　　についての境界条件としては、 これが Geopotential に関する自由振動にたいする下部境界条件である この境界条件から自由振動の鉛直方向の解は変数分離形 として、上の境界条件の式に対応するものとして大気の中も のような形をもつであろう。これはあきらかに外部波である（鉛直には伝播しない形）。鉛直方向の方程式は１章の赤道波と同様に だから、計算すると ここで、　　　　　　　を用いて、

16 書き直すと、 これが自由振動の等価深さである。浅水波としたときの速度　　　　　　　　　は音速に等しい。等価深さとしての h は おおよそ１０kmとなる。スケールハイトと比較して、　（比熱比）だけ、大きくなる。 その時の変数分離した形で地球の回転を考慮した球面上の式を書くと， 　　　は自由振動の時の等価深さである． この式をいつものように（球だから経度をつかう） のように分離すると，南北の構造を決めるLaplaceの潮汐方程式がでてくる。　　　　　　として，　　　　　　　　　　 これが長周期波動に関する球面上の式である。　　　　　 ２章で示した式との違いは、東西波数sの部分と１項の　　　が異なる。

17 s=1（東西波数１）の、等価深さ h（横軸） と振動数との関係図。 図の(a), (b)はh>0のときで、(c),(d)はh<0を示す。左図は東に伝わる波で右図は西に伝わる波である。東に伝わる波はKelvin波とか慣性重力波と書いてある。西に伝わる波は惑星波とかRossby-重力波とか慣性重力波と書いてある。 　図(b)の○印のドット入りはあとの例にある自由振動の５日wave（惑星波のところ）、h=10kmに対応している。 　南北に高次のRossby波に対応しては、波の周期がゆっくりになっている。Rossby波の分散式で 等価深さ 振動数 周期 東進慣性重力波 5日 惑星波 南北波数が大きくなれば、ｃ は小さくなり、周期はゆっくりなるであろう。 ⊗は１日潮汐の第１モードを示す（１３章）

18 ５−４：現実大気中でのRossby波の自由振動
前節でみたように、自由振動は外部波の１つで、なめらかな固体境界を満たす振動である．境界条件を考慮すると、 　　　　　　　　　　　　の形なので、鉛直方向には指数関数的に大きくなる。 赤道対称 Hirooka(1992, Tenki)：s=1, 周期=5日の成層圏 1hPa での自由Rossby波、西に伝播しているようす。 この５日波が南極大陸で励起されているという話しがある、Cheong and Kimura(1997, JAS)：500mbで、軸が南北に傾いている（自由振動の解では南北に位相軸は傾かない）。また、Miyoshi and Hirooka(1999, JAS）では、湿潤過程のheatingの重要性を指摘している。 Maddenにより解析された周期５日の自由Rossby波、s=1で赤道対称な基本モード（海面気圧高度場）の水平構造、右は鉛直分布。破線は理論的な鉛直構造、それよりゆっくり振幅が増大している。

19 衛星観測結果：5day波の全球鉛直構造 equinox seasonの5day波 平均東西風、boldは高度偏差最大
Lawrence and Randel, 1996, JGR 平均東西風、boldは高度偏差最大 1977 Oct. 温度 位相 西向き 約81kmでの温度スペクトルの緯度分布 高度の振幅と位相

20 〜100km 近年の衛星観測結果から：5day wave SABER data, 2002年の春での5 day 赤道対称Rossby normal modeの温度分布 Garcia, 2005, JAS

21 水蒸気変動に見られる準５日波 （normal modeらしい）
Sonnermann et al., JGR 2008 水蒸気変動に見られる準５日波 （normal modeらしい） マイクロ波をもちいた観測とLIMAモデルとの比較 74km 2003年May/Juneでの観測された水蒸気変動、ALOMAR(69N, 16E), Norway, マイクロ波 5/1 2003 7/1 モデルの水蒸気変動の時間緯度断面図

22 極中間圏雲(Polar Mesospheric Cloud)変動に見られる準５日波
Merkel et al., JASTP 2009 CIPS, 中間圏の氷エアロのミーspacecraft上の、雲imageと粒子サイズ測器 Cloud Imaging and Particle Size Instrument on the Aeronomy of Ice in the Mesosphere spacecraft 2007年の夏、それぞれの緯度における、solstice期のzonal dayly頻度変化、5日程度の変動 PMCアルベードの波数ー振動数スペクトル解析、正の振動数が、西向き伝播をしめす、波数 s=1で５日あたりに強いシグナルが見える 経度180度での、雲アルベードの変化の時間変動

23 データ解析による、10day wave（自由振動）の構造、1981, Apr.の1hPa高度偏差、50m間隔
１０日波について データ解析による、10day wave（自由振動）の構造、1981, Apr.の1hPa高度偏差、50m間隔 東西波数１で南北に基本的な３つのモードの理論的な構造分布図、10日波＝図の(1,2)モードは反対称モードの最も簡単な南北構造をもつ。 振幅の時間的な変化の様子

24 赤道角運動量（Equatorial Atmospheric Angular Momentum, EAAM）に現れる１０日波
Feldstein, JAS, 2006, 風の寄与 質量トルク EAAMの位相（NCEP/NCAR 冬期dataから） EAAMのtendencyに線形回帰した地表面気圧ー＞赤道に関して反対称に近い構造 EAAM vectorの約10-day の変動にs=1 反対称Rossby normal modeが対応しているらしい ＜ー地表面気圧の変化によるマスバランスの違いで変動をもたらしているよう。 グリニッジにおける赤道角運動量の位相 ー＞１０日程度の振動

25 観測で求められた、16day wave自由振動、赤道に関して対称モードである。Hirooka and Hirota, 1985, J
観測で求められた、16day wave自由振動、赤道に関して対称モードである。Hirooka and Hirota, 1985, J. Atmos. Sci.　冬半球の方にシグナルが見える（非対称的）。下は４月の振幅変動（対称的） 線形計算に用いられた冬の基本風 90km 赤道対称的 計算で求められた（Salby, 1981, J. Atmos. Sci.）１６日波の緯度-高度図、冬の極に大きな振幅をもっていて、構造がいびつになっている。 1980年

26 補足：対流圏の赤道Rossby波 Dima, Wallace, and Kraucunas, 2005, Dima and Wallace, 2007, JAS： 東西に周期的な、定常のforcing 年平均の降水量（影、mm/day）と150hPaの高度（線）および水平風ー＞赤道Rossby波的構造がみえる n=1赤道Rossby波的 Kelvin波的 定常応答の解、Matsuno, 1966, JMSJ 10S-5N平均の、eddy成分の高度（色）と東西、鉛直速度（矢、10m/sと1cm/s）、ERA-40 data