１３章：中層大気中の大気潮汐 大気潮汐のはなしを最後にしておこう。

Presentation on theme: "１３章：中層大気中の大気潮汐 大気潮汐のはなしを最後にしておこう。"— Presentation transcript:

1 １３章：中層大気中の大気潮汐 大気潮汐のはなしを最後にしておこう。 大気中の加熱（基本的には、太陽による１日加熱）によって、大気がどのように応答するかの問題ということになる、さらにそれの非線形な振る舞いを議論する。 １３−１：対流圏中の潮汐をすこし 対流圏における潮汐の様子を眺めてみよう。ECWWF dataの解析から得られたものである。Hsu and Hoskins, 1989, Q. J. R. Met. Soc., 115, 夏 矢の長さが１日潮汐の地表圧力振動の大きさを、矢の向きが最大振幅のlocalな時間を示し、北向きが00時、東向きは06時、南向きが12時を示す。赤道の振幅が大きいが、矢の向きは東向きで、圧力最大はほぼ6時になっている。夏半球側が振幅は大きい。西に進み、１日で一周する成分がdominant。地（海）表面の影響により、１日潮汐は局所性が高い。 冬 1mb〜10m

2 対流圏１日潮汐変動の鉛直と緯度の構造 高度場振幅と位相の、緯度高度断面図、 赤道での高度場の様子、経度高度断面図、 H L -90° 90°
等値線は2m間隔 H L -90° 90° 経度 DJF 86/87年、00GMTと06GMTの状況で、等値線は5m、shadeは負の偏差、 位相の等値線は30°ごと

3 １３−２：古典的潮汐理論 太陽放射の大気加熱による潮汐の問題を、最も簡単にしたものが古典的な潮汐理論であろうから、線形の方程式はこれまでたびたび引用してきたものと同じであろう 太陽放射加熱として、時間と東西方向は　　　　　　　　　　のような形を仮定しよう。 太陽と同期する（migrating)１日潮汐はs=1, σ=1/2、半日潮汐はs=2, σ=1に対応 １章や５章と同じように、変数分離の方法を用いて、水平方向の式と鉛直方向の式に分離し（東西波数s、振動数は与える）、南北方向の方程式（Laplace’s tidal equation）として、 が得られる。 波数、潮汐の振動数は与えるとして、そのもとで、固有値 hn と対応した固有関数を求めることになる。赤道β面で、南北モードnで与えられていたように、個々の固有関数（Hough functions)の区別は指標nで構造の違いが与えられる。

4 ５章の図で、s=1が与えられていて、例えば１日周期の西に伝播する潮汐を決めたい時、１日に交差するものが個々のモードと解釈される。
古典的潮汐理論（続） ５章の図で、s=1が与えられていて、例えば１日周期の西に伝播する潮汐を決めたい時、１日に交差するものが個々のモードと解釈される。 西向きの波 東向きの波 １日 ⊗はmigrating１日潮汐の第１モードを示す h < 0 １日

5 １日の大気潮汐の第１モードは、 に対応した構造をもち、上の図から想像されるように、赤道域に集中した構造をしている。
1日潮汐モードの南北構造と固有値 図は、西に伝播する1日潮汐モードの圧力偏差の赤道対称南北構造、等価深さの大きさにより、南北の広がりが異なっている。具体的な計算は、Lindzen, 1966, MWRを参照 第１モード 鉛直波長として、おおよそ28km 負の等価深さ １日の大気潮汐の第１モードは、　　　　　　　　　　　　　　に対応した構造をもち、上の図から想像されるように、赤道域に集中した構造をしている。

6 モードの違いに相応して、鉛直方向の構造は下記の式で与えられる。
鉛直方向の構造を決める式： モードの違いに相応して、鉛直方向の構造は下記の式で与えられる。 連続の式 熱力学の式 ー＞ ここでJnは大気加熱のnモード成分を表す、オゾンや水蒸気による大気加熱を評価することになる。 整理すると、鉛直方向の構造方程式 加熱がない場合には、これまでたびたび議論したように鉛直に伝播する　　　　　　かそれともexponential 的かの方程式になっている。ちなみに、１日潮汐の第１モードの鉛直波長は約28kmに対応する。 下部の境界条件は Jが与えられて、南北にはHough modeに展開、鉛直構造は上式により解かれる、Lindzen, 1967, QJRMS参照 最終的にはこれらの重ね合わせであり、各モードの鉛直伝搬性が異なり、構造は複雑になっている。

7 古典論の１つの解： 仮定される加熱Jの鉛直および緯度分布
オゾン オゾン 水蒸気 赤道 仮定される加熱Jの鉛直および緯度分布 計算結果（点線）と観測（実線）の１日潮汐、８°Sでの北風の振幅（左）、位相（右）, Reed et al., MWR, 1967 春秋分における１日潮汐の温度振幅の緯度依存性、Lindzen, 1967, QJRMSから

8 １３−３：Migrating diurnal Tideの観測結果
衛星データからの中層大気中の１日大気潮汐について 4Nでの温度鉛直分布、52Wのlocal noonを基準にしてある。破線は衛星データのs=1 成分のみ、実線はnon-migrating成分も含めたもの、xはForbesモデル(1982)の結果、○はロケット観測から評価、 Lieberman, JAS,1991 S=1西向き１日潮汐波の温度の緯度高度図、1979, 3/4 から5/4の期間、1hPaの高度で2K程度である、1991 Lieberman, JAS

9 下部熱圏高度のmigratingな１日潮汐
衛星データの解析から導出された、１日潮汐南北風の緯度高度断面図、Hays et al., 1994, JASから T 平均東西風や散逸過程の入った、線形モデルによるmigrating１日潮汐の高度分布、Forbes, 1982, JGR より精緻化された（tidal forcingなど)モデルによる4月の１日潮汐の南北風振幅分布、Hagen et al., 1995, GRL 下図はモデルの計算結果から

10 金星大気中に存在するであろう、１日（117地球日）潮汐温度シグナルの観測結果：
Schofield and Taylor, 1983, QJRMS 金星大気中に存在するであろう、１日（117地球日）潮汐温度シグナルの観測結果： 温度の経度緯度分布：波数s=1とs=2成分が混在している構造となっている、軸が鉛直方向に傾いている70-90kmの高度で鉛直方向に伝播している。 高度 z（km） 赤道から30Nまでの平均 東西方向（x） 潮汐波動が金星大気の平均東西風生成に重要な役割を果たしているという話は後にのべる。 10

11 火星大気の１日潮汐 Wilson and Hamilton, JAS, 1996 火星大気の大循環モデルで潮汐を表現をしてみる 太陽同期のs=1, 1日潮汐の温度振幅（線）と位相（色）、Banfield et al., 2000, JGR 0.1hPa 6.1hPa 北半球冬の、s=1 西に進む１日潮汐の温度振幅 東西平均温度場 モデルの東西平均温度場

12 １３−４：大気潮汐波による平均東西風生成について
大気潮汐波も、これまでの議論から推測されるように、平均東西風を生成することが可能である。 大気加熱によって、潮汐波が生成されている状況を示している。図では、東方向の波が生成される。これまでの議論から、東に伝播する波は西風運動量をになって鉛直方向に伝播し、波動が消散されるところで、西風運動量を生成する。今、加熱のみであり、力は加われていないので、波の消散されるところで、西風運動量を生成すると、加熱されるところで、逆向きの東風運動量を生成することで、正味の運動量＝０となる。Fels and Lindzen, Geophys. Fluid Dynamics, 1974 波による加速 平均東西風 講義では、波が鉛直に伝播しながら散逸している領域での、波の伝播方向の東西風生成の話し（１：下部熱圏における赤道域の東風）と、波が生成されている領域での伝播とは反対方向の東西風生成の話し（２：金星大気における成層圏での高速東西風）をしておこう。 ー＞　大気加熱の動き(東向きとする）

13 １：地球大気下部熱圏における大気潮汐波による平均東西風生成
大気中の１日潮汐波の非線形効果を計算した結果。赤道域は東風になっている。観測の東風と対応しているよう。一方中緯度では西風が生成されている。 平均東西風の緯度ー高度断面図において、下部熱圏に注意してほしい。赤道域で東風が吹いている。 Miyahara, 1978, J. M. S. J.

14 Miyahara、JMSJ, 1978 によって解かれた問題
潮汐の線形方程式として 潮汐南北風の振幅と位相の高度分布、緯度３５度の場所 の高度 z=0で１日潮汐の第１モードのみを強制している（z=0は現実大気では94kmに対応していて、下部熱圏の潮汐が散逸する場所で解かれている） 上付き１日潮、下付き第１モード 南北風の振幅の緯度構造

15 潮汐の散逸によって生成される平均流の式が時間積分される
潮汐（ ’印）による平均流へのforcingは以下のよう 潮汐は南北に構造をもっており、それに対応して、Fuは赤道域における東風加速、および中緯度での西風加速 数値結果としての平均風　

16 ２：金星大気の東西平均温度および東西風： ー＞
左図は東西に地球一周平均した温度分布。右図の東西風の評価には旋衡風（cyclostrophic wind）バランス+静力学平衡を用いている。赤道域で、約70kmの高度に100m/sを超える東西風が吹いている。また、中緯度付近には140m/sの高速風が見える。これは41.5kmの高度で観測の風を仮定して、求めてある。 (Newman et al., 1984, JAS)。 ー＞ 温度 東西風 16

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

18 １３−５：non-migrating（太陽の動きと同期しない）潮汐のいくつか
s=1の東向き１日潮汐波の緯度高度図、1979, 3/4 kara5/4の期間、1991 Lieberman, JAS、non-migratingの eastward（日の動きとは逆向き）波 東向きs= -1のHough modeによる解析、この時期は、南北指数 3, 4がもっとも卓越していたようである

19 非同期の１日振動：SABER観測の例（高波数の例）
Garcia et al., 2005, JAS fがー、東向き 〜120km 〜87km 緯度 〜70km 時間 波数２の１日振動の時間変動 １日振動の東方向に伝播する、波数２（左）と波数３（右）の振幅と位相の緯度高度図、2002, Jun Jul.をもとに解析、s=3は大きな振幅をもつ

20 非同期の１日振動、波数５と６ Garcia et al., 2005, JAS 中間圏界面 0.01mb １日振動の東方向に伝播する、波数５（左）と波数６（右）の振幅と位相の緯度高度図、2002, Jun Jul.をもとに解析

21 対流活動で励起されるnon-migrating潮汐
Hagen and Forbes, JGR 2002 対流活動で励起されるnon-migrating潮汐 緯度高度の2次元線形モデルで、基本場は変化しないとして、東西波数、振動数を固定して解いてある。Forcingとして観測された対流圏におけるlatent heat releaseをあたえる。対流にともなう潜熱放出なので、東向き成分も含むであろう 計算の為に与えられるHeatingの鉛直構造 13 wave numberでcompositeしたもの、s=6 のW6 から E6までの範囲

22 潮汐の理論計算結果の例 Hagan and Forbes, 2002, JGR
線形の方程式でs=6までとってあるモデルで、対流活動によるforcingの成分と基本場の中の波の鉛直伝播性、および散逸に依存して、上層での波の構造が決まるであろう 4月の１日振動の対流活動に対するモデルの温度応答（上図）をしめす。W1が西に進むs=1のmigrating潮汐の応答であり、E3はs=3の東向き潮汐の応答を示す。下図は南北風応答を示す。E2, E3は観測と同程度の振幅のよう

23 熱圏高度での潮汐の理論計算例（続） E3モードの115km高度での、温度振幅の季節変化
115km高度の４月における１日振動の南北風と温度の水平構造、波数３の構造が見えている。黒が12UTでgrayが06UT時、

24 観測結果の文章として、（Banfield et al., 2000, JGR）
１３−６：潮汐波のcoupling Ls=270北半球冬の東向きの１日潮汐の温度振幅（モデルの結果）ー＞東に伝播していて、Kelvin波が第１モード、似ている所もある　Wilson and Hamilton, 1996, JAS １：火星大気の１日東向き潮汐波 観測結果の文章として、（Banfield et al., 2000, JGR） clearな空間構造が見えにくいので図には示さないが、東向き１日潮汐の振幅が２K程度はある、 西向き 東向き 西向き波数１の１日潮汐と波数２の山岳で、東向き波数１の１日潮汐がつくられるschematic図、Forbes et al., JGR, 2002から、丁寧な説明はZurek, JAS, 1976を参照

25 Palo et al., GRL 2007 ２：東に伝わる準２日波
s＝２で東向き２日波：７章で議論したs=3の西向き２日波（Rossby-gravity波）とは異なるシグナルの生成について 温度 西向き Jan22-29 西向きs=3の波（Rossby-gravity波）とs=1の西向き1日潮汐波のcouplingで、s=2で東向きの２日波が生成と考えられている、下記のような理屈であろう ー＞解析的により詳しくみているようである 上は80kmでの、s=3, 2日波とs=2, 2日波の時間緯度断面図、下は40Sでの時間高度断面図 緯度分布としての、波の対応をみていて、相互作用をしているかの確認の１つになるのだろう