440 likes | 667 Views
大気大循環モデルで表現される赤道域降水活動の構造解析. 北海道大学 大学院理学研究科 地球惑星科学専攻 山田 由貴子 DM2semi 2006/02/09. 0.もくじ. はじめに 研究背景 (観測と数値実験) 研究目的とこれまでの取り組み 時空間スペクトル解析 赤道波理論 , 時空間スペクトル 解析結果 鉛直構造解析 wave-CISK 理論 アニメーション解析 スペクトルフィルターコンポジット解析 まとめ. 研究目的とこれまでの取り組み. 1. 赤道域降水活動に対する認識. Nakazawa (1988) 赤道域の降水階層構造
E N D
大気大循環モデルで表現される赤道域降水活動の構造解析大気大循環モデルで表現される赤道域降水活動の構造解析 北海道大学 大学院理学研究科 地球惑星科学専攻 山田 由貴子 DM2semi 2006/02/09
0.もくじ • はじめに • 研究背景 (観測と数値実験) • 研究目的とこれまでの取り組み • 時空間スペクトル解析 • 赤道波理論, 時空間スペクトル • 解析結果 • 鉛直構造解析 • wave-CISK 理論 • アニメーション解析 • スペクトルフィルターコンポジット解析 • まとめ
1.赤道域降水活動に対する認識 • Nakazawa (1988) • 赤道域の降水階層構造 • クラウドクラスター, スーパークラスター, MJO
2. APE; Aqua-planet experiment project AGUforAPE CSIRO 水惑星比較実験 T42程度の解像度の結果 2005/04/20-22, Reading GSFC LASG K1JAPAN MGO NCAR UKMO
3. これまでの取り組み • 水惑星条件におけるモデルパラメタ依存性調査 • 潤沢な?計算機資源 (スパコン) を使って, とにかくたくさん?の実験をする • パラメタ空間マップを作る • まずはどんな降水分布パターンを取りうるのか? • そのパターンはどのくらいのパラメタ空間で存在するのか • その上で, 理論的考察をして整理する (できたらいいな) • キーワード: wave-CISK, 水蒸気移流 • 水平解像度依存性実験 • モデル: AFES (AGCM for Earth Simulator; CCSR NIES AGCM ver 5.4.02 の並列化効率を高めた), 積雲スキーム使用せず • 解像度:T39L48 (水平格子 330km), T159L48 (水平格子 80 km)
4.水平解像度依存性 (降水赤道上 x-t 図 ) T39L48_non T159L48_non
6.赤道波 (Matsuno,1966) • Matsuno (1966) • 赤道ベータ面浅水方程式系 44.7 m/s 1300 km
6.赤道波 (Matsuno,1966) • Matsuno (1966); 分散曲線と解の構造 慣性重力波 混合ロスビー重力波 ケルビン波 ロスビー波
7.時空間スペクトル(Wheeler and Kiladis 1999) • Wheeler and Kiladis (1999) • NOAA: Twice-daily, 2.5grid, 1979/01--1996/08, 15N-15S • 「バックグラウンドスペクトルとの比」を取ることで, ケルビン波シグナルと MJO シグナルを分離 • 等価深度 • 鉛直第一モード, 鉛直第二モード, … 30 m/s 22.3 m/s 15 m/s 11 m/s 9 m/s 2500km
8.時空間スペクトル (降水、比湿) 比湿 (下層: σ=0.85) 降水
9.時空間スペクトル (温度、東西風) 比湿 (上層: σ=0.25) 比湿 (上層: σ=0.25)
10.時空間スペクトル(Wheeler and Kiladis 1999) • Wheeler and Kiladis (1999) • NOAA: Twice-daily, 2.5grid, 1979/01--1996/08, 15N-15S • 「バックグラウンドスペクトルとの比」を取ることで, ケルビン波シグナルと MJO シグナルを分離 • 等価深度 • 鉛直第一モード, 鉛直第二モード, … Lz=19 km 14 km 10 km 6.8 km 5.6 km
11.時空間スペクトル(Wheeler and Kiladis 1999) Lz=19 km c =30 m/s 10 km 15.8 m/s 14 km 22.3 m/s 6.8 km 10.9 m/s 5.6 km 8.9m/s
12.moist-Kelvin • Reduced gravity • 浅水系で考えると… • 重力波の速度 • 加熱を受けると遅くなる
13.wave-CISK モデル 鉛直モード 加熱分布 • Yoshizaki(1991a), 保坂 (199?) • 赤道ベータ面、ブジネスク系、線型 変数の座標位置 加熱の パラメタリゼーション 境界条件
19.偏東風移流 -7 m/s
20.理論から得られる予想 • wave-couple していれば • 位相速度は第1モードと第2モードの間くらい • 鉛直構造の位相は進行方向に向かって逆に傾く(温度、風速、ジオポテンシャル) • wave-couple していなければ • 位相速度は第1モードより遅い? • 鉛直構造の位相の傾きは見られない • 背景風を考慮すると? • 位相速度は、遅くなる?速くなる?
? ? ? ? • 解析上の問題 • 参照点の取り方の任意性が大きい • 抽出したい構造がぼやける、イメージがつかみにくい 25.コンポジット解析 • 解析上の重大な問題 • とにかく面倒くさいっ!!! • たくさんの実験ケース …やってらんないっ!!!
26.スペクトルフィルターコンポジット解析26.スペクトルフィルターコンポジット解析 • Wheeler and Kiladis (2000) • Data: NOAA OLR, NCEP Reanalysis • スペクトル空間のある領域だけ取り出して実空間に戻す • 上記データに対して, OLR の極大値を参照点としてコンポジットする • 我々の実験結果に対して • 批判もあるだろうが、イメージ作りのためにやってみる価値あるかも
27.スペクトルフィルターコンポジット解析27.スペクトルフィルターコンポジット解析 1.0 -4 -1 3 gravity filter advect filter kelvin filter -28 -15 15
38.まとめ • まとめ • 。。。。 • 課題 • 解析上のこまごま • wave-CISK の説明の精度, 簡単モデルから複雑モデルへ • 積雲パラメタリゼーション問題 • wave-CISK ではない可能性の検討
39.移流を考慮した分散曲線 -7m/s とした場合