メンバー

気候モデルにおける重要要素に関する基礎的検討と非専門家への翻訳，さらにその予測精度の評価に向けた検討会
2015年6月3日
吉田英生
京都大学 工学研究科 航空宇宙工学専攻
sakura@hideoyoshida.com【主査・幹事】
吉田英生
京都大学工学研究科航空宇宙工学専攻教授
【副主査】
草野完也
【委員】
江守正多

国立環境研究所地球環境研究センター
気候変動リスク評価研究室長
増田耕一
海洋研究開発機構地球環境変動領域
地球温暖化プログラム 主任研究員
大西 領
海洋研究開発機構横浜研究所地球シミュレータセン
ターマルチスケールモデリング研究グループ研究員
桑原不二朗 静岡大学大学院機械工学専攻教授
富村寿夫
熊本大学自然科学研究科産業創造工学専攻教授
村田 章
東京農工大学工学府 機械システム工学専攻教授
吉田篤正
大阪府立大学工学研究科機械工学分野 教授

申請書より抜粋 （１）






もともと地球全体で１℃程度という平均温度の増減を過去に遡って確認すること自体が極めて困難な課題であり、IPCCの主論拠は気候モデル（General Circulation Model、GCM）による将来予測に基づいている。
しかし、熱流体の数値シミュレーションの実際をよく知っている本会会員の視点からは、モデル化に不可避の任意性を直感的にも理解できよう。
機械系に軸足を置く本会の大部分の会員の研究対象と比較すると、そのスケールと複雑さの相違は著しいが、現象の各要素や素過程はまさに本会会員の専門とするところでもある。
そこで、気象学者にその検討を全面的に任せるのではなく、本会会員の視点から中立的な第三者としてメスを入れることは極めて重要であると考えられる。

名古屋大学太陽地球環境研究所総合解析部門教授

申請書より抜粋 （２）

気候モデルは、太陽からのふく射をエネルギー源として、地球の海洋・大陸・大気中におけるさまざまな熱物質輸送過程をモデル化して総合した極めて複雑な数値シミュレー
ションコードである。
総じて、数値シミュレーション結果は一般市民にも強いインパクトで迫るが、その精度・不確かさに関して厳正なチェックの裏付けとともに透明に提示されることは皆無といっても過言ではない。
本研究は、そのような数値シミュレーションの代表例の一つである気候モデルを対象として、それを言わば「まっ裸」にして一般社会への橋渡しする最初の試みとしてもその意義を強調したい。

パリティ 2012年7月号

気候モデルにおける検討要素マトリックス【たとえば最近話題になっているhiatus 】
20世紀後半以降、地球全体の地表気温（以下、全球平均地表気温）は上昇の傾向を示しており、2001年以降の10年間の平均気温は、1961～1990年の平均に比べて約0.5℃高くなっています。しかし、21世紀に入ってからの気温上昇率は10年あたり0.03℃とほぼ横ばいの状態を示しており、こうした温暖化の停滞状態はハイエイタスと呼ばれています。ハイエイタスの要因には諸説ありますが、その原因は解明されていません。 http://www.aori.u-tokyo.ac.jp/research/news/2014/20140901.html

大気大循環モデル AGCM(1995)の解読
大気大循環モデル (Atmospheric General Circulation Model)
東京大学気候システム研究センター (CCSR) と国立環境研究所（NIES) の共同で1995年に開発されたモデルは、海洋とのカップリン
グや二酸化炭素循環は考慮せず、これらを境界条件として所与のものとするものであるが、気候モデルの骨格を理解する入口としては好適である（とはいっても門外漢には非常に複雑であるが）ので、そのマニュアルを読み解くところから着手している。
第1章 モデルの概要
第2章 力学過程
表面気圧で正規化した圧力で鉛直方向に２０層に分割して数値計算

第3章 物理過程○積雲対流過程 (3.2節)○大規模凝結過程 (3.3節)○放射過程 (3.4節)○鉛直拡散過程 (3.5節)○地表フラックス (3.6節)○地表面・地中過程 (3.7節)○重力波抵抗 (3.9節)

全球3次元大気大循環モデルの概要
方程式系
予報変数
水平離散化
鉛直離散化
放射雲過程
積雲対流
鉛直拡散
地表flux地表面過程
重力波抵抗

全球プリミティブ方程式系
風速、温度、地表気圧、比湿、雲水量
スペクトル変換法
s系 (Arakawa and Suarez, 1983)
2ストリームDOM/adding法(Nakajima and Tanaka, 1986)
雲水予報 (Le Treut and Li, 1991)
簡略化Arakawa-Schubertスキーム
Mellor and Yamada (1974) level2Louis (1979) バルク式
Miller et al. (1992) の対流効果
多層熱伝導バケツモデル
McFarlane (1987)

方程式系
地表気圧

東西風

u
 Fx D  Fx Pt
南北風

気温

物理過程

v
 Fy D  Fy PtT Q D  Q Ptp S M D  M Pt

○積雲対流過程 (3.2節)
 凝結熱が大気の成層を不安定化させる場合数値モデルの空間分解能より小さな積雲が発生する
 積雲群全体の効果をパラメーター化して表現する必要があるq S D  S Pt

○大規模凝結過程 (3.3節)
 層雲系の降水など凝結を起こしても成層が不安定化しない場合
 通常は特別なパラメータ化をせず、過飽和を飽和に変え余分な水蒸気を降水として潜熱を大気に与える

比湿地中温度Tgt Q g D  Q g P

力学過程
大気運動の移流など (添字 D の項)
物理過程
雲・放射などの各プロセス (添字 P の項)

積雲対流過程 （１） 積雲は空気に対して軽い水蒸気を含むので浮力で上昇しつつ周囲の空気を側面から巻き込む（エントレインメント）
 逆に浮力を失うと空気を周囲の放出する（デトレインメント）
 また、エントレインメントの周辺には反作用で下降流が生じる
 したがって、解像度の粗い数値モデルでは、相対湿度の格子点代
表 (平均) 値が100％以下でも降水が生じる→このため、積雲群の統計的知見に基づいて格子点値から降水量を決める（積雲対流のパラメータ化におけるクロージャー問題）

Arakawa andSchubert, 1974

○放射過程 (3.4節)○鉛直拡散過程 (3.5節)○地表フラックス (3.6節)○地表面・地中過程 (3.7節)○重力波抵抗 (3.9節)

積雲対流過程 （２） 鉛直運動にともなって発生する熱 (対流有効位置エネルギー) の大部分は積雲スケールの運動エネルギーに変換される。
 これをすべて格子点スケールの有効位置エネルギーに変換すると、モデルは過剰な運動を表現→凝結熱の鉛直配分が大規模な運動エネルギーに関わる有効位置エネルギー生成効率をモデル化する必要

要するに
 どれだけの水蒸気を凝結させるか（前スライドの問題）
 その熱をどのように鉛直に配分するか（上記の問題）
 水蒸気をどのように、雲・降水に再分配するか（上記の問題）

放射過程 （１）

放射過程 （２）dI d   ka  ks I   ka  B T   ks  P ,    I   4  dz4

熱放射を評価するためには、放射伝達方程式をプランク分布の波長帯域にわたっても積分する必要があるが、不連続的な大気成分の吸収帯を有限の波長域に分割して各波長域で平均した値で近似する
大気圏外から測定した
熱帯晴天域での外向きの長波放射スペクトル
晴天時の太陽放射スペクトル100m

放射過程 （３）

10m

5m

放射過程と鉛直拡散過程： ポイント
放射過程でのポイント （検討課題）
 気体ふく射のモデルと雲やエアロゾルによる散乱の観点から検討
 二酸化炭素の増加によるふく射への直接的な影響は、吸収係数に現れる
 二酸化炭素の増加によるふく射への関節的な影響は、水蒸気量の変化を通して雲生成に影響を与え、散乱係数に変化が現れる
鉛直拡散過程でのポイント （検討課題）
 鉛直方向へのサブグリッドスケール乱流拡散を拡散係数を用いて与える レイノルズ応力、乱流熱流束の基礎式を代数方程式に簡略化して用いる
 密度変化による大気の不安定成層による寄与をRichardson数の関数として拡散係数を算出