11章:赤道域中層大気における平均東西風の長周期変動 ー準2年振動と半年振動についてー

Slides:



Advertisements
Similar presentations
ケルビン波( Kelvin Wave ) 境界が無い場合、回転を感じる重力波(慣性重力波)の 最大振動周期は慣性周期であった。つまり、慣性周期 以上の時間スケールを持つ重力波は、 “ 境界が無い場合 ” には存在しない。 “ 境界があれば ” 、どうなるのであろう か? 境界があれば、慣性周期以上の時間スケールを持つ重力.
Advertisements

海洋流体力学 2014 海洋流体力学とは、海洋に関する流体力学。本講義では、 海洋のみならず、大気も含めた地球流体力学について学ぶ。 Fluid Dynamics( 流体力学 ) Geophysical Fluid Dynamics (地球流体力学) 目標 海洋・大気大循環のイメージを描けるようにする。
YohkohからSolar-Bに向けての粒子加速
大気科学入門 - 金星大気東西風の超回転について -
電磁気学C Electromagnetics C 7/27講義分 点電荷による電磁波の放射 山田 博仁.
自己重力多体系の 1次元シミュレーション 物理学科4年 宇宙物理学研究室  丸山典宏.
13章:中層大気中の大気潮汐 大気潮汐のはなしを最後にしておこう。
数値気象モデルCReSSの計算結果と 観測結果の比較および検討
熱帯太平洋における季節内スケールの 赤道波動特性の解析 5AOOM007 萩原 右理 指導  轡田 邦夫 教授.
Rossby 波動が関係する成層圏平均東西風の変動例
第10章:成層圏突然昇温など --惑星波動による平均東西風変化の例-- 10−1: EP-フラックスによる解析
高橋研究室 惑星大気 中層大気 化学モデリング 気象力学
第6章 波の平均場への作用 そして運動のエネルギーと、、  質も形も進度も位置も時間も みな因縁が…
冨川喜弘 (国立極地研究所・トロント大学)
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
Rossby 波動が関係する成層圏平均東西風の変動例
第10章:成層圏突然昇温 --惑星波動による平均東西風の変化について-- 10−1:Charney -Drazin の定理 
成層圏突然昇温の 再現実験に向けて 佐伯 拓郎 神戸大学 理学部 地球惑星科学科 4 回生 地球および惑星大気科学研究室.
渦位(Potential Vorticity)と角運動量の保存
ジェット気流が延々と吹き続けている理由を理解する
13章:中層大気中の大気潮汐 大気潮汐のはなしを最後にしておこう。
北海道大学 理学部 地球科学科 惑星物理学研究室 4年 馬場 健聡
8章:赤道域中層大気における平均東西風の長周期変動 ー準2年振動と半年振動についてー
風成海洋大循環 (準地衡流渦位方程式+エクマン層の力学)
第6章:熱帯域の不安定について(個別的?)かき混ぜか 6ー1:成層圏内の傾圧不安定(?)で起こっている例
第9章:中層大気の不安定擾乱について -不安な雲のうかび出て ふたたび明るく晴れるのは-
惑星大気大循環モデル DCPAM を用いた 地球大気に関する数値実験
電磁気学C Electromagnetics C 7/13講義分 電磁波の電気双極子放射 山田 博仁.
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
*大気の鉛直構造 *太陽放射の季節・緯度変化 *放射エネルギー収支・輸送 *地球の平均的大気循環
第5章 波の生成と シアー中の重力波について 北上山地はほのかな幾層の青い縞をつくる-
CMIP5マルチ気候モデルにおける ヤマセに関連する大規模大気循環の 再現性と将来変化(その2)
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
2.伝送線路の基礎 2.1 分布定数線路 2.1.1 伝送線路と分布定数線路 集中定数回路:fが低い場合に適用
赤道QBOの影響の統計的有意性 ― 大標本法に基づいた評価 ―
近年の北極振動の増幅 Recent Arctic Oscillation amplification
2005/06年冬季における 北半球環状モード変動の 予測可能性について
三重大学・大学院生物資源学研究科 共生環境学専攻 地球環境気候学研究室 教授 立花義裕
ロスビー波( Rossby wave) 渦度 (vorticity) 順圧非発散流(絶対渦度の保存) ポテンシャル渦度(渦位)
Fermi Bubble と銀河中心の巨大構造
YT2003 論文紹介 荻原弘尭.
川崎浩司:沿岸域工学,コロナ社 第2章(pp.12-22)
潮流によって形成される海底境界層の不安定とその混合効果
電磁気学C Electromagnetics C 7/17講義分 点電荷による電磁波の放射 山田 博仁.
講義ノート(ppt)は上記web siteで取得可 #但し、前日に準備すると思われるのであまり早々と印刷しない方が身の為
冬季北大西洋振動が 翌冬の日本の気候に与える影響
講義ノート(ppt)は上記web siteで取得可 #但し、前日に準備すると思われるのであまり早々と印刷しない方が身の為
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
渦位(Potential Vorticity)と角運動量の保存
電磁気学Ⅱ Electromagnetics Ⅱ 8/11講義分 点電荷による電磁波の放射 山田 博仁.
Iida, O., N. Kasagi and Y. Nagano
速度ポテンシャルと 流線関数を ベクトルで理解する方法
Johnson et al., 1999 (Journal of Climate)
CMIP3 マルチモデルにおける熱帯海洋上の非断熱加熱の鉛直構造 廣田渚郎1、高薮縁12 (1東大気候システム、2RIGC/JAMSTEC)
冬期極域成層圏対流圏循環の変動に おける赤道域QBOの影響の 統計的有意性
竜巻状渦を伴う準定常的なスーパーセルの再現に成功
川崎浩司:沿岸域工学,コロナ社 第4章(pp.58-68)
地球フロンティア・モデル統合化領域 渡辺真吾
潮流によって形成される海底境界層の不安定とその混合効果
落下水膜の振動特性に関する実験的研究 3m 理工学研究科   中村 亮.
北海道大学 理学部 地球科学科 惑星物理学研究室 B4 近藤 奨
2006 年 11 月 24 日 構造形成学特論Ⅱ (核形成ゼミ) 小高正嗣
北極振動の増幅と転調は 何故20世紀末に生じたか? Why was Arctic Oscillation amplified and Modulated at the end of the 20th century? 地球環境気候学研究室 鈴木 はるか 513M228 立花 義裕, 山崎 孝治,
将来気候における季節進行の変化予測 (偏西風の変化の観点から)
卒論中間発表 2001/12/21 赤道の波動力学の基礎 北海道大学理学部 地球科学科 4年 山田 由貴子.
冨川喜弘 (国立極地研究所・トロント大学)
K2地球システム統合モデル 成層圏拡張の進捗について
磁気リコネクションによる Alfven波の発生
従来研究 本研究 結果 南極大型大気レーダーPANSYで観測された大気重力波の数値モデル再現実験による力学特性の解明
第3章:赤道波(重力波も含めて) -擾乱表現の1つの方法- <ー赤道成層圏は波で満ちあふれている?
Presentation transcript:

11章:赤道域中層大気における平均東西風の長周期変動 ー準2年振動と半年振動についてー 11章:赤道域中層大気における平均東西風の長周期変動    ー準2年振動と半年振動についてー 西風 赤道域半年振動 10hPa 東風 北半球冬 100hPa 1月の平均東西風 赤道域下部成層圏準2年振動(10m/s間隔)

11−1:Eliassen-Palm の定理の一般化 波による平均東西流加速の一般論を、Andrews and McIntyre (1976, J. Atmos. Sci. ) から引用しておく 擾乱の式として(ブシネスク流体近似、β平面、静力学平衡)、 東西平均流の変化の式は以下のように書かれる。 ここで、 右辺:unspecified forcing terms この方程式には、重力波および、Rossby波動が含まれる。 ー>この式を変形すること: 南北熱フラックス収束が子午面循環をつくり、それにコリオリ力がかかり運動量を変化させる ー>波が担っている運動量という考え方 平均東西風は緯度、および高さの関数であり、対応した温度場も緯度、高度依存性をもつ。

Eliassen-Palm flux: の項を東西風の変化の式にくりこむと、以下の式になる。 Eliassen-Palm flux をpseudo-運動量フラックスと呼ぶこともある <ーそれの収束が東西風の変化に対応する の項(Eliassen-Palm flux divergence)の変形から、平均東西流の加速として以下の式が導かれている。ここで、平均東西流の式で *のついた項が、小さい近似である(波が定常で、散逸やcritical levelがないときはゼロになる) ここで、 は波動に伴う流体粒子の南北変位をあらわし、     である。 波に対しての散逸や外力(1項や2項)、transienceの時(3項)、またcritical level (2、3項)のところで東西風が変化することを示している。 Richardson数が大きく、赤道β平面で、位相速度 cをもつ赤道にtrapする波動についての近似では:

Eliassen-Palm定理の破れの最も簡単な例 のようになり平均流が時間とともに変化していく。このとき西風を生成可能。  波は定常ではあるが、散逸されつつある場合;定常で散逸されつつあるのだから、例えば対流圏で常に強制されている。散逸として、同じ係数 a のRayleigh摩擦と Newtonian冷却を考え、散逸は小さいとする。  ブシネスク近似で、静力学平衡をみたす2次元重力波の鉛直波数 m について を解いた例:1つの東に伝わる波のみの、平均東西風の時間発展の様子をみたもの。Plumb, 1977, J. Atmos. Sci. 左図はフラックスの時間変化。 いま、波は平均流に対して東に動いているとしている。このとき上方に伝播する波の解はWKB近似的に以下のように表される(Lindzen, Dynamics in Atmospheric Physics, 1990) 時間 鉛直座標や時間は無次元化されている 鉛直EPフラックスは今の場合、 基本流に対して西向きの波の運動量フラックスは負である。 となり高さの関数。また、 このとき  東向きの波と西向きの波が両方あるとどうなる? 条件によっては西風と東風で振動しそう

波が散逸しつつあるとき、東西平均流が変化することを述べた。その典型的な例が赤道域の下部成層圏に存在する準2年振動と考えられている。 11—2:準2年振動(QBO)のありよう  波が散逸しつつあるとき、東西平均流が変化することを述べた。その典型的な例が赤道域の下部成層圏に存在する準2年振動と考えられている。 QBOに関する観測結果をいくつか述べる(cf. Andrews et al. ,1987) 西風と東風の繰り返し、上から伝播してくる(40kmくらいからか) 準2年振動は年振動と関係があるらしい。QBOの西風が下降するとき、季節的振動である半年周期振動の西風(equinoxのとき)と同期している時もあるよう 71-72 59-60 準2年振動の周期は22ヶ月から34ヶ月と一定ではない。平均の周期は28ヶ月くらいPlumb(1984)より 東風 西風持続at44hPa 西風 Pascoe et al., JGR, 2005では、太陽 Minで西風がより持続;東風がおりる(20-44hPa へ)時間が太陽 Maxで2ヶ月短い という統計的結果となっている。 半年振動(約48kmの高さ)と準2年振動、Wallace(1973)より、5m/sごと、

赤道からすこしはずれると、 が正のとき(西風が高さとともに大きい時)、 赤道からすこしはずれると、   が正のとき(西風が高さとともに大きい時)、 北半球で  が負だから赤道の方が温度が高い。このとき、熱力学の式から(T’>0として)    赤道にtrapされた現象である のようであろうから    は下降流となり、西風shearのときは鉛直移流により、はやくQBOは下降する。 QBO 西風と東風の下方伝播の違い: 下方伝播の速さは約1km/月で西風の下方伝播の方が幾分速い。これは、子午面循環の違いで説明されるであろう。 西風 shear warm 下降流 地衡風近似と静力学平衡からくる温度風の関係と、熱力学の式におけるNewton冷却と断熱鉛直運動のバランスの式: から 図:Plumb and Bell, 1982, QJRMSの2Dモデルより

11−3 準2年振動の力学的説明  赤道下部成層圏の準2年振動を波と平均流の相互作用の考え方でモデル化してみる。まずは、赤道上のみを取り扱う。東西方向に一様な風(平均流)を支配する運動方程式は右の式: QBOを生成しているといわれる波動について:準2年振動の西風(上層)が下りてくる時期で、周期15日程度の擾乱がある。これは東向きの波で西風運動量をもっており、散逸するとき西風を生成する。Wallace-Kousky wave(1968, J. Atmos. Sci. )とも呼ばれ、対流圏で生成された赤道ケルビン波といわれている。  西風 東風 対応 上方伝播ケルビン波の位相関係 赤道下部成層圏のケルビン波の時間−高度断面図(上が東西風で下が温度)。1963年の夏、場所はカントン島(南緯3度)

赤道下部成層圏にケルビン波はあって、波数1で振幅が最大で10msー1くらいはあるらしい。 波の生成は対流と大規模波動がcoupleして出来たものらしいが、明確ではない。わかり易い考えとして波動と第2種不安定(台風のメカニズム)を結びつけたWave-CISKを使ったHayashi(1970)があるが、この理論も潜熱放出パラメーターに強く依存する  赤道下部成層圏にケルビン波はあって、波数1で振幅が最大で10msー1くらいはあるらしい。 時間 東西 4000km 21km高度、1958, Apr. 15-30、ほぼ赤道上、影は南風成分のところに Wallace and Kousky, 1968, JASから 西向きの波について:図はYanai and Maruyama(1966, J. M. S. J.) により発見されたRossby-重力波の伝播の様子を示したもの。東西波数4くらいで、位相速度は25msー1程度、振幅は2〜3msー1の振幅をもっている。観測されているRossby-重力波の振幅はそれほど大きくない、この波は散逸するとき東風を生成する。 Holton and Lindzen(1972)はこの2つの波を使って準2年振動をモデルで再現したが、RG波の振幅を大きく与えている。 大循環モデルで表現されたRossby-重力波、Hayashi and Golder, 1994, J. Met. Soc. Japan. 波の振幅はv=0.5m/s程度である。

Holton-Lindzen(1972)のモデル:Kelvin波とRossby-重力波を使い、ニュートン冷却で波を減衰 方程式は、南北には積分された式で、 となる.HLでは、上層の半年振動 を与え,28km以上で,           はKelvin波とRossby-重力波の運動量フラックス   は Kelvin波: QBOを再現するためには赤道上のRossby-gravity waveの南北風振幅は下部境界で6msー1 程度与えている。 HLの1次元モデルで再現されないものとして西風の下方伝播が東風より速いことがある。前に述べたように鉛直と南北の2次元子午面循環を考慮すれは説明可能であろう(cf. Plumb and Bell, 1982 ) Rossby-重力波については分散式から Rossby-gravity重力波の場合,鉛直運動量フラックスは、     ではなくて の南北平均である 高度 位相速度30msー1のケルビン波及びRossby-gravity波を使ったのはQBOの南北スケールと波の南北スケールが1500km程度と同じくらいということである。 例えば、Kelvin波として 年 から c=30m/s として le は1000km程度になる。 10m/s間隔、shadeが西風

Plumb(1984)による、位相の下方伝播と振動 2つの東西に伝播する内部重力波を用いた振動の仕方: 波の波長は40000km(波数1の赤道ケルビン波に相当)、位相速度は30msー1(東向き、及び西向き)と仮定。ここで、約6msー1の東西風の振幅を仮定する。水平スケールが大きいと、この程度の振幅が必要である。それに対応して下部境界での運動量フラックスはHLと同じく、 という値を用いる。 計算結果が図に示してある。周期約1000日程度の準2年振動的な構造になっている。

3次元のmechanistic model で再現した例: 大振幅のKelvin波とRossby-重力波を下部境界で与えるとQBOは再現される。西風の方が早く下降している。 T=1800 days でのKelvin波の東西風。振幅が観測に比べて大きい(15m/sくらい)、2.5m/s間隔。 Takahashi and Boville, 1992, JAS T=1500 days でのRossby-重力波の南北風。振幅が観測に比べて非常に大きいこと、 2.5m/s間隔

GCMの中のQBO: Takahashi(1999, GRL) 現実的なQBOが再現されている。モデルでは、様々な重力波でQBOが生成されている モデルは、水平分解能が60kmであり、約200km以上の重力波が直接表現される。このモデルの範囲内で(対流のパラメータで重力波の生成が異なるであろう)、 西風シアー:東向き赤道波動は、25-50%の寄与をもち、内部重力波は50-75%程度の寄与をになっている。 東風シアー:西向き赤道波動は10%程度、中緯度からのRossby波動は10-25%程度、主に内部重力波が寄与 最近のモデルQBO, Kawatani et al., 2010, JAS、色はEP-fluxの発散、周期は15ヶ月

QBOの南北スケールに関して: QBOの南北スケールは、1500km程度である。赤道波動のみでは説明できていたが、重力波が主要因とすると、別の考え方も必要であろう。 ERA-40 dataからのQBO振幅分布、Pascoe et al., 2005, JGR 準2年振動の振幅(実線)と位相(破線)の緯度−高度断面図、Wallace(1973)より

Haynes(1998, Q. J. R. M. S.)による説明: の式を思い出そう(2章) p, NやNewton冷却の高さ依存を落とすと 数値実験による確認:南北に広いforcing(上図)にも関わらず、生成される東西風は赤道域のみとなっている(下図)。 ここで、準2年の変動に比べて、Newtonian damping係数は大きいのでおもな応答は 応答の南北スケールを L とし、時間のスケールをT、forcingの鉛直スケールをDとして左辺1項のz微分項に反映するとする。1項と2項で0に近いものが応答しやすいだろう。スケール的には のようであろう。f=βL とすれば、上式は のように南北スケールが決まる。   =2x3.14/2/3x107=10-7   =10-6   =0.1、Dを10kmとすると   2x10-2 x104/2x10-11=1013     3000km        --> L=1500km程度で、観測の値に近い値となる。

11—4:QBOに関係した幾つかの話題 QBO-likeな、流れの交代する実験: Plumb and McEwan (1978, JAS)、流体力学的(相似性)に興味深い。 t=150mで左の方への流れ、t=170mで右の方の流れが見える。膜の振幅等が下図とは異なるが、 実験装置:下でStanding波を作る。 左が実験で得られた振動、右が理論の結果

中間圏のQBO: 中間圏QBOが見つかっている、Burrage et al. (1996, J. G. R.) 中間圏QBO 成層圏QBO 重力波をパラメータ化したモデルで再現された中間圏QBO(y-z 2Dmodel)、 Mayr et al. (1997, J. G. R.)、振幅は大きくない

QBOと惑星波動との関係:Yamashita et al QBOと惑星波動との関係:Yamashita et al., 2011, JGR  QBOの西風位相のとき、東風位相のときと比べ、惑星波動の活動として、高緯度成層圏でフラックス偏差としては下向きで、発散的(西風加速)である。 QBOは中緯度成層圏に影響を及ぼしているよう(Holton and Tan, 1980, 図はYamashita et al., 2011, JGRから ) 下部成層圏QBOが西風のとき、冬の極夜Jetの西風が統計的に強くなっている、JRA-25 dataから 矢はEP flux偏差(QBO西風ー東風)、線は発散の偏差、色のあるところは有意性をしめす。

QBOの対流圏への影響、Crooks and Gray, J. Climate, 2005から 上図のQBOに対応して、統計的に有意なシグナルが見られる。北半球対流圏にQBOと関係するanomalyが見え、赤道20kmの高度でのanomalyがつながっているようにみえる。 下部成層圏QBOの西風位相ー東風位相の偏差、Pascoe et al., JGR, 2005から 18

補足:対流圏との関係 Maruyama and Tsuneoka ( 1988 )は ENSO と QBO の関係を調べている。ENSO のときケルビン波の活動度が強まり西風の下降が早まっているようだと述べている(1987年のENSOの時,東風の持続が短かった) 熱帯域の深い対流(OLRと対応)と下部成層圏の東風shear(低温、上昇流)とが関係あるという話しと矛盾しない、という論文もある (Collimore et al., 1998, GRL ) 冬季惑星波動の振る舞いの違いにより、QBOの東風位相のとき、北アジア域の下層大気がwarming anomalyとなっている (Chen and Li, 2007, JGR)

QBOは物質変動にも存在する(下はオゾンQBOの例をしめす)、Randel and Wu, 1996, JASから Lu et al., 2008, JGR 影 - positive 中緯度のオゾンQBO(赤道QBOと逆位相) 全オゾン 高度緯度パターン

木星の準4年振動 赤道域の標準温度からの偏差の時間変化、□(実線)が赤道、ダイアモンド記号が14S, △が14Nである。 高度緯度の2次元モデルによる準4年振動の再現、この図は、Kelvin波とRossby重力波のforcingを与えている。 7.8μmのbrightness温度の時間変化、おおよそ20hPa高度

11−5:赤道域成層圏の半年振動 1:成層圏界面付近の半年振動 Dunkerton, JAS,1978: 基礎方程式はこれまでと同様に とする.  は,半年振動の東風成分のみ生成するようにしてある. 10m/s間隔 西風を加速するKelvin波について、観測でみつかっているような位相速度c=50m/s,東西波数は1を選ぶ, 東風加速について:非線型の子午面移流      ,中緯度からの惑星波動の効果,重力波が考えられている.どの程度の割合かはまだ分かっていないよう。 西風加速について,重力波が大事であるといわれている.   NCAR GCMの半年振動:西風はおもにKelvin波と書いてある,西風が弱いよう−>たぶん重力波が足りない <-対流のパラメータのせいであろう。 Sassi, F., R. R. Garcia and B. A. Boville, 1993: The stratopause semiannual oscillation in the NCAR community climate model. J. Atmos. Sci., 50, 3608-3624.

2:中間圏界面付近の半年振動 成層圏界面の半年振動とは位相が逆転している。成層圏の半年振動の風をかんじて、逆方向の重力波が80kmまで伝わっていきそこで、波が壊れて逆位相に半年振動が生成されているらしい GFDL- GCMの中の半年振動。この場合は西風がよく再現されている。<ー 対流のパラメータが異なる。対流調節が用いられており、調節が瞬間的におこり、そのため多くの重力波が生成されているようである。 Hamilton and Mahlman, 1988, J. Atmos. Sci.

Solsticeの西風は20m/s程度、equinoxの東風は -30m/s程度の振幅をもっている。 レーダー等で評価された中間圏半年振動: Garcia et al., 1997, JGR Solsticeの西風は20m/s程度、equinoxの東風は -30m/s程度の振幅をもっている。 80km 東風 クリスマス島(2N)のレーダーで評価されたMSAO 東風 CCSR/NIES/FRCGC GCM(T213L256)での半年振動、Watanabe et al. , 2008, JGR MSAOの東風はある程度は再現、ただ、このモデルでのtop境界近くである HRDI衛星データからのMSAO

Antonita et al., 2008, JGRでは、中間圏の半年振動は重力波がmainly causedであると言っている インド、Trivandrm(8.5N, 77E)にある流星レーダーをもちいて評価している 月平均東西風 50m/s/3month〜17m/s/month 東西風の加速、実線は月平均のSAO風変化から、dotted はestimated された加速、2004年6月から2007年5月、1W:1年目の西風 2-3時間周期の短周期重力波にともなうu’w’の季節変化 ー> 右の加速の評価

低分解能気候モデルの結果:Richter and Garcia, 2006, GRL 水平2度の分解能、重力波の効果をパラメータ化して入れてあるモデル ad. EPD G モデルの中間圏半年振動 による Solstice(西風位相)では、重力波に加えて、子午面循環、分解されている波動によるEP-flux dovergence(特に2日波)が寄与、equinox(東風位相)は全forcingが小さい EP-flux Divergence, 黒は全成分、赤は2日波(7章)、青は1日潮汐の寄与