大気再解析データを用いた 帯状平均場の解析

Slides:



Advertisements
Similar presentations
生物統計学・第 5 回 比べる準備をする 標準偏差、標準誤差、標準化 2013 年 11 月 7 日 生命環境科学域 応用生命科学 類 尾形 善之.
Advertisements

CMIP5 気候モデルにおける三 陸沿岸の SST の再現と将来予測 児玉安正・ Ibnu Fathrio ・佐々木実紀 (弘前大学大学院・理工学研究科)
多次元データ 解析・可視化ソフトウェア GAVE A Grid Data Analyzer and Viewer, GAVE 竹本 和彰 北海道大学理学部地球科学科 地球流体力学研究室 4年 2004 年 2 月 3 日.
過去 100 年に観測された 夏季日本の気候変動 気象研究所 遠藤洋和 第 10 回ヤマセ研究会.
CMIP5 気候モデルにおける ヤマセの将来変化: 海面水温変化パターンとの関係 気象研究所 気候研究部 遠藤洋和 第 11 回ヤマセ研究会 1.
Ruby on Rail の紹介 石渡正樹 Ruby on Rails とは? スクリプト言語 Ruby で書かれた web アプリケー ションフレームワーク 作者 –Devid Heinemeier Hansson という人だそうです ( 詳 しいことは知りません.
ITPASS Informational Training program with a spirit of self-help オプション課題の概要 高橋芳幸.
Doxygen ~ IGModel を一例にした, 数値モデルのドキュメンテーションにおける Doxygen の利用
大気科学入門 - 金星大気東西風の超回転について -
北大における Super-SINET 接続と利用: 2004 年度報告
JRA-55再解析データの 領域ダウンスケーリングの取り組み
表計算ソフトで動作するNEMUROの開発
島田照久(1) 沢田雅洋(2) 余偉明(2) 川村宏(1)
数値気象モデルCReSSの計算結果と 観測結果の比較および検討
北日本における4月と8月気温の強い相関関係とその時間変動(2)
数値モデルの出力データをどのように取り扱っているか?
電脳Rubyプロジェクト概要 堀之内 武(京大生存圏研究所),
実証分析の手順 経済データ解析 2011年度.
相関係数 植物生態学研究室木村 一也.
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
力学的ダウンスケールによる2003年東北冷夏の アンサンブル予報実験
成層圏突然昇温の 再現実験に向けて 佐伯 拓郎 神戸大学 理学部 地球惑星科学科 4 回生 地球および惑星大気科学研究室.
貧困と出産の関係.
北海道大学大学院理学研究科地球惑星科学専攻 地球流体力学研究室 M1 山田 由貴子
永井晴康、都築克紀(原研)、植田洋匡(京大防災研)
惑星大気大循環モデル DCPAM を用いた 地球大気に関する数値実験
davis / gtool4 プロジェクト その背景と野望
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
*大気の鉛直構造 *太陽放射の季節・緯度変化 *放射エネルギー収支・輸送 *地球の平均的大気循環
CMIP5マルチ気候モデルにおける ヤマセに関連する大規模大気循環の 再現性と将来変化
CMIP5マルチ気候モデルにおける ヤマセに関連する大規模大気循環の 再現性と将来変化(その2)
赤道QBOの影響の統計的有意性 ― 大標本法に基づいた評価 ―
脳活動に関するデータ データの種類 データの特徴 脳波・脳磁図・fMRI画像 脳活動とパフォーマンスの関係はきわめて冗長。
近年の北極振動の増幅 Recent Arctic Oscillation amplification
2005/06年冬季における 北半球環状モード変動の 予測可能性について
高次元データの解析 -平均ベクトルに関する検定統計量の 漸近分布に対する共分散構造の影響-
海氷が南極周辺の大気循環に与える影響 地球環境気候学研究室  緒方 香都 指導教員:立花 義裕教授.
ヤマセによる冷夏をターゲットにした アンサンブルダウンスケール予報実験
2.温暖化・大気組成変化相互作用モデル開発 温暖化 - 雲・エアロゾル・放射フィードバック精密評価
バルク法について ~deepconv を用いて地球の積雲対流の数値計算をするにあたって~
YT2003 論文紹介 荻原弘尭.
2010 年北極振動の 冬から夏への極性反転と 猛暑の連関 ―北極振動と猛暑と今年の夏―
講義ノート(ppt)は上記web siteで取得可 #但し、前日に準備すると思われるのであまり早々と印刷しない方が身の為
熱帯海上における降水特性による 降水・循環の将来変化パターンの マルチモデル間の違い 廣田渚郎、高薮縁 (東大AORI) 2011/6/9.
冬季北大西洋振動が 翌冬の日本の気候に与える影響
講義ノート(ppt)は上記web siteで取得可 #但し、前日に準備すると思われるのであまり早々と印刷しない方が身の為
南北両半球間を横断する 水蒸気輸送と降水量との関連性
2009年秋の北極海ラジオゾンデ観測によって観測された 大気の順圧不安定とメソ渦列
共生第二課題における 陸域生態系炭素循環モデルの研究計画 名古屋大学大学院 環境学研究科地球環境科学専攻 市井 和仁
気候シナリオモデルを用いた将来のヤマセ発生可能性について
気候モデルのダウンスケーリングデータにおける ヤマセの再現性と将来変化
菅野洋光 (農研機構東北農業研究センター) 渡部雅浩 (東京大学大気海洋研究所)
Johnson et al., 1999 (Journal of Climate)
CMIP3/CMIP5気候モデルにおける ヤマセに関連する大規模大気循環の再現性 ~モデル解像度による違い~
冬季、東シナ海・日本南方海域における 温帯低気圧の発生に関する気候学的研究
CMIP5気候モデルにおける ヤマセの将来変化
気候モデルのダウンスケーリングデータにおけるヤマセの再現性と将来変化2
CMIP3 マルチモデルにおける熱帯海洋上の非断熱加熱の鉛直構造 廣田渚郎1、高薮縁12 (1東大気候システム、2RIGC/JAMSTEC)
冬期極域成層圏対流圏循環の変動に おける赤道域QBOの影響の 統計的有意性
竜巻状渦を伴う準定常的なスーパーセルの再現に成功
潮流によって形成される海底境界層の不安定とその混合効果
MIROC5による将来のヤマセの再現性について(2)
2006 年 11 月 24 日 構造形成学特論Ⅱ (核形成ゼミ) 小高正嗣
地球環境気候学研究室 513M230 松本直也 指導教員 立花義裕
北極振動の増幅と転調は 何故20世紀末に生じたか? Why was Arctic Oscillation amplified and Modulated at the end of the 20th century? 地球環境気候学研究室 鈴木 はるか 513M228 立花 義裕, 山崎 孝治,
将来気候における季節進行の変化予測 (偏西風の変化の観点から)
地球環境気候学研究室 谷口 佳於里 指導教員:立花義裕 教授
夏季日本における前線帯の変動と その天候への影響
卒論中間発表 2001/12/21 赤道の波動力学の基礎 北海道大学理学部 地球科学科 4年 山田 由貴子.
K2地球システム統合モデル 成層圏拡張の進捗について
Presentation transcript:

大気再解析データを用いた 帯状平均場の解析 塚原大輔 北海道大学 大学院理学研究科 地球惑星科学専攻 地球惑星流体科学講座 地球流体力学研究室 このようなタイトルで発表させていただきます. 地球流体力学研究室 修士二年 塚原大輔です. どうぞよろしくお願いいたします. <<2004 年度 DM2 セミナー>> 2005 年 2 月 23 日

大目次 (1)データの蒐集と整備, 教科書的図表の作成 (2) Ruby による気象解析 (3)ハドレー循環の年々変動

データの蒐集と整備, 教科書的図表の作成

再解析とは 数値予報/データ同化手法を用いて過去の観測データから,大気循環場と境界条件(地表面フラックス, 放射フラックスなど)のデータセットを作成すること 1990 年代初めから欧米の気象機関が先行して作成開始 NCEP/NCAR(NRA=NCEP/NCAR Re-Analysis, NRA1) ECMWF(ERA=ECMWF Re-Analysis, ERA-15/40) 気象庁/電中研(Japanese Re-Analysis 25 years, JRA-25) 長期間の均質なデータを提供することが目標 現在, 品質の検査とデータの利用(科学研究)が再帰的に行われている ERA-40がもっとも先行, NRA1→NRA2ヘ, JRAは初期段階 NRA1 : 1991 年開始(Kalnay, et al. 1996) ERA40 : 2000 開始(http://www.ecmwf.int/research/era/Project/Plan/) JRA : 2001 年開始(http://www.jreap.org/publications/plan.html) 近年, 欧米や日本の気象機関にて再解析データと呼ばれるデータセットが作成, 配布されはじめました. 再解析とは過去の観測データを元に, 気象機関で現業的に用いられている数値予報やデータ同化手法を用いて 品質の一様な大気大循環場と境界条件のデータセットを作成することです.

現実に(もっとも)近い数値モデル出力結果 再解析データは「現実」データか? 同化に用いる数値モデルの「癖」に大きく依存 物理的/統計的につじつまが合うよう調節されてはいる 再解析データによって取り込まれている観測データが異なる よって再解析データは観測データというよりも 現実に(もっとも)近い数値モデル出力結果

「現実」を見極めるには? 再解析の製法を調べる 再解析の材料となる観測データを押さえる 異なる再解析データ, 観測データの絵を描いて相互比較 数値モデル(スキーム等) 再解析の材料となる観測データを押さえる 異なる再解析データ, 観測データの絵を描いて相互比較

以上を踏まえて… 惑星大気の標準データの蒐集 教科書的な標準図表を作成 上記プロダクトをWeb 上で公開 再解析データの相互比較のために 異なる気象機関の再解析データを統一フォーマットに 教科書的な標準図表を作成 再解析データの相互比較のために 上記プロダクトをWeb 上で公開 http://www.gfd-dennou.org/arch/dcchart/ 想定する利用者 数値モデルを用いて大気大循環の研究を行う人 基本的には林研の人々(特に森川さん?) 詳しくは上記 URL をお見せして説明します

Ruby によるデータ整備, 可視化・解析

気象業界の可視化・解析事情 標準的なツール/ライブラリ 現状に対する不満 コンパイル型言語はうっとおしい Fortran/DCL, GrADS, GMT, etc… 解析 (Fortran|C)/数値計算ライブラリ, GrADS 現状に対する不満 コンパイル型言語はうっとおしい GrADS, GMT は可視化・解析専門 データハンドリング, 文字処理, CGI など一つの言語で実現したい なんか無機質な感じがする 理想 強力な汎用スクリプト言語 製作したライブラリ/スクリプトの整理, 管理がたやすい 多くの人が「開発者」になれる敷居の低い言語

オブジェクト志向スクリプト言語 Ruby Ruby 長所 型なしスクリプト言語 言語仕様レベルでオブジェクト指向サポート 強力な言語機能 開発効率, 保守性良好. 強力な言語機能 文字処理, イテレータ, ソケットプログラミングetc 可視化/数値解析関連の外部ライブラリも登場 FORTRAN や C のライブラリを呼び出し可能 コーディングしてて楽しい! 設計思想が「プログラム作成を楽しくする言語」 Ruby 短所 言語自体が発展途上中 拡張ライブラリの数はまだ Perl や Java, Python にかなわない それ以外は思いつかない…(誰か教えてください)

気象解析に利用するライブラリ NArray Ruby/NumRu 多次元データ解析/数値計算用モジュール群, ユーティリティ群 多次元 数値 配列クラス Ruby でも大量の数値を扱う計算が、簡単かつ高速にできる Ruby/NumRu 多次元データ解析/数値計算用モジュール群, ユーティリティ群 京都大学生存圏研 堀之内武 氏を中心とするコミュニティで開発中 電脳 Ruby プロジェクト (http://dennou-k.gfd-dennou.org/arch/ruby/index-j.htm) (自己主張)塚原も参加 一部資源紹介 Ruby/DCL(Dennou Club Library): 日本の気象業界で使われているグラフィック機能中心のライブラリの Ruby ラッパー Ruby/NetCDF: NetCDF IO の Ruby ラッパー GPhys: 格子状に離散化された物理量を扱うクラス

可視化/解析プログラム netCDF 形式のデータ 算出した物理量 電脳 Ruby プロジェクト公開のライブラリを利用して可視化・解析プログラムを作成 netCDF 形式のデータ クイックビューワ, 属性書き換え, 2 項演算, 次元平均 算出した物理量 質量流線関数, EP フラックス(京大堀之内氏と共同), 南北エネルギー輸送量, etc…

ハドレー循環の年々変動

目次 過去の研究 目的 データ ハドレー循環強度 El Nino との関連 熱帯降水活動との関連 東西風ジェットとの関連 中緯度擾乱活動度との関連 まとめ

目的 再解析データを用いて長期間の帯状平均場の様子を描画する (上記の一環として)ハドレー循環強度の変動と関連する物理量の変動の関連を調べる 熱帯降水活動 ハドレー循環を駆動する積雲対流活動の指標 中緯度東西風ジェット 循環強度が強ければ角運動量保存的にジェットは強くなるはず 中緯度の擾乱の活動度 擾乱による角運動量輸送 => 擾乱によってハドレー循環が駆動

過去の研究 特にハドレー循環に着目して, 再解析データが明らかにする大気大循環の描像を得ることが目標 (今回は NRA1) ハドレー循環にまつわる最近の研究 Oort and Yienger(1996) ハドレー循環強度とENSO , 諸量の変動の関連を相関係数を用いて示した(ゾンデデータ). ハドレー循環強度以外の時系列は示していない Wang(2002) ローカルハドレーの循環強度が強いと中緯度ジェットも強くなると議論(NRA1). 帯状平均東西風ジェットについては触れていない. Dima and Wallace(2003) ハドレー循環を赤道対象・非対称成分に分けた. ただし地表面気圧より高い気圧面上の格子点の値も用いている( NRA1) きちんと現実の大気に対応する格子点のみによってみてみるとどうであろうか?

データ NCEP/NCAR Reanalysis 1(NRA1) 期間:1978年 12 月 – 2003 年12月 http://www.cdc.noaa.gov/cdc/reanalysis/reanalysis.shtml より取得 期間:1978年 12 月 – 2003 年12月 FGGE 実施以降を選択(信頼度の高い均質なデータを選択) 地表面気圧より高い気圧面上の値は使わない 現実大気と対応する格子点値のみ使用

ハドレー循環強度の平均年変動からの偏差の時系列 北半球は正, 南半球は負 それぞれ絶対値が大きいほど循環は強い ハドレー循環強度の時系列 ハドレー循環強度の平均年変動からの偏差の時系列 El Nino の発生期間を左図に描く(この際手書きでもいい) 黒線: 北半球ハドレー循環強度 赤線: 南半球ハドレー循環強度 黒線: 北半球ハドレー循環強度 赤線: 南半球ハドレー循環強度 (右)ハドレー循環強度の時系列 (右)ハドレー循環強度の平均年変動からの偏差の時系列

El Nino 発生時期(気象庁のサイトより転載) ハドレー循環強度の強い時期は概ね El Nino と一致する Oort and Yienger(1996) とは 1982/83 年を除いて整合的 NRA1 では 82/83 年の El Nino 時にも強度が強くなる El Nino 1982(MAM) - 1983(JJA) 1986(SON) - 1988(DJF) 1991(MAM) - 1992(JJA) 1993(MAM) - 1993(JJA) 1997(MAM) - 1998(JJA) 2002(MAM) - 2003(DJF) El Nino 発生時期(気象庁のサイトより転載) http://www.data.kishou.go.jp/climate/elnino/mikata/elnino_table.html

熱帯降水活動との関連 北(南)半球の循環強度が強い時 南(北)半球ITCZ 付近の降水量が増加 低緯度全体の平均降水量はあまり変動していない => 熱帯大気に対する潜熱加熱はさほど年々変動していない 南北の加熱差による強制によってハドレー循環は駆動される? 帯状平均降水量の平均年変動からの偏差 帯状平均降水量の緯度時間変動 黒線: N5-N15 平均値, 赤線: S5-S15 平均値 緑線: N30-S30 平均値

中緯度東西風ジェットとの関連 時系列を見る限り, 循環強度との対応はよく見えない Oort and Yienger(1996) の算出した循環強度とジェットの相関係数は北半球で 0.39, 南半球で –0.34. 200 hPa 面東西風の平均年変動からの偏差 黒線: N20-N40 平均値 赤線: S20-S40 平均値

EP フラックス鉛直成分の平均年変動からの偏差の時系列 フェレル循環強度の平均年変動からの偏差の時系列 擾乱の活動度との関連 時系列を見る限り, 循環強度との対応はよく見えない Oort and Yienger(1997)はフェレル循環強度が弱まると主張 しかし佐藤(1998)は,理想的な境界条件下ではハドレー循環強度とフェレル循環は比例すると結論 EP フラックス鉛直成分の平均年変動からの偏差の時系列 フェレル循環強度の平均年変動からの偏差の時系列 黒線: N30-N60 平均値 赤線: S30-S60 平均値 黒線: 北半球フェレル循環強度 赤線: 南半球フェレル循環強度

まとめ ハドレー循環が強いときの諸量との関連は? 今後の課題 Ruby を用いた可視化・解析 El Nino の時期は循環強度は強い 降水は時系列で比べても関係がわかる 降水の集中化による強制加熱 東西風ジェット, 擾乱の活動度とは時系列で見ても対応付けは見られない 今後の課題 相関解析や統計的有意性の検定を行い, 変動の有意性を吟味 複数のデータソースを解析, 比較 NRA1 に加えERA40 や ERBE などのデータと比較. データ自体のチェックを行う 循環強度の強い年の子午面構造を解析 Ruby を用いた可視化・解析 (感想として)非常に有益. まとめに際して注意点.

参考文献 CDC Data Management Group, The NCEP/NCAR Reanalysis Project at the NOAA-CIRES Climate Diagnostics Center. http://www.cdc.noaa.gov/cdc/reanalysis/reanalysis.shtml. Dima, I. M., Wallace, J. M., 2003: NOTES AND CORRESPONDENCE On the Seasonality of the Hadley Cell. J. Atmos. Sci.,60, 1522-1527. European Centre for Medium-Range Weather Forecasts, 2001:ERA-40 Project Plan. http://www.ecmwf.int/research/era/Project/Plan/Project_plan_TOC.html 気象庁, 2005: エルニーニョ/ラニーニャ現象の情報. http://www.data.kishou.go.jp/climate/elnino/mikata/elnino table.html 気象庁・電力中央研究所, 2001: JRA-25 長期再解析プロジェクト. http://www.jreap.org/index.html 増田耕一, 1989: 大気データの4 次元同化による気候系の解明. 日本の科学と技術, 30, 54-59. Oort, A.H., Yienger, J. J., 1996: Observed Interannual Variability in the HadleyCirculations and Its Connection to ENSO. J. Climate, 9, 2751-2767. 佐藤正樹, 1998: 地球大気の子午面循環と角運動量輸送. 日本惑星科学会誌「遊星 人, 7, 91-101. 塚原大輔, 石渡正樹, 小高正嗣, 高橋芳幸, 山田由貴子, 森川靖大, 林祥介, 地球流体電脳倶楽部DCCHART プロジェクト, 2003: 大気大循環. http://www.gfddennou.org/arch/dcchart/atmos global/, 地球流体電脳倶楽部. Wang, C., 2002: Atlantic Climate Variability and Its Associated Atmospheric Ciculation Cells. J. Climate, 15, 1516-1536.

以下メモ

6. 25 年平均値帯状平均場 北半球: 6.0e10 , 南半球:-7.0e10 kg/s Dima and Wallace(2003): NRA1 , 地表面以下の格子は抜いてない 北半球: 6.0e10 +few , 南半球:-7.4e10 kg/s + few ゾンデの見積もり(Oort and Yienger, 1996) は 北半球: 9.2e10 , 南半球:-7.4e10 kg/s

(脱線)なぜ Ruby は楽しいのか?

熱帯域ラジオゾンデステーション分布の一例 (Oort and Yienger, 1996 より転載)