大気データの4次元同化による気候系の解明: 図表

図1: 質量・エネルギー・角運動量の収支解析の説明図。ある地域について、大気全層(大気の上端から地表面まで)を含む「箱」を考える。保存量である質量、エネルギーまたは角運動量の箱の中の収支を考えると、箱の中の保存量Xの時間変化ΔX/Δtは、上端から(Ftop)、地表面から(Fsfc)および横方向から(Fside)のその量の正味の供給の和に等しい。このことを使い、4項のうち3項をデータから知ることができれば、残りの項が推定できる。 [本文へ]

fig1.png


図2: 大気大循環統計量を求める2つの方法 (文献1より)。 たとえば、観測された風の南北成分vと比湿(大気中で水蒸気の占める質量比)qに基づき、時間平均した水蒸気の南北輸送q'v'を規則的に配置された格子点で求めることを考える。ただし上線は時間平均を示す記号である。方式Aでは観測地点で時系列的統計計算をしてq'v'を求め、その結果を空間内挿する。方式Bでは、まず空間内挿で格子点でのqとvを求め、それから時系列的統計を計算する。それぞれの長所、短所は本文参照。 [本文へ] [表1の関連個所の本文へ]

fig2.png


図3: 断続型および連続型の4次元同化の手順の説明 (文献2より)。 断続型(a)では、モデルの力学・物理過程に基づく時間積分(予報)と、データの空間内挿とが別のプロセスになっており、交互に行なわれる。たとえば6時間予報の結果を仮の推定値として空間内挿を行なう。連続型(b)では、時間積分の過程の各時間ステップで、モデルの積分結果に、データを空間内挿したものを線形結合で混ぜていく。なお、図の中にあるOIとは「最適内挿法」といい、観測値および予報値の線形結合による空間内挿法で、観測および予報の誤差を考慮して結果の誤差を最小にするように線形結合の重みを決めるものである。また、NMIとは、「ノーマル モード イニシャリゼーション」といい、大気の運動を重力波モードとロスビー波モードに分解し、固有周期の短い重力波モードの成長率を小さくするように風速・温度・気圧の場を調節する手続きである。 [本文へ]

fig3.png


図4: 大気と海洋による年平均の北向きエネルギー輸送 (文献4より)。 1978年12月から1979年11月までの1年間の観測データによる。 単位は1015ワット。 大気と海洋の合計輸送量(実線)は、ニンバス7号衛星による大気上端での放射収支観測から求めたもの、大気による輸送量(破線)は、ECMWF(ヨーロッパ中期天気予報センター)が4次元同化の方法で作ったFGGE IIIbデータから計算したもの、海洋による輸送量(点線)はこの2つの差として求めたものである。 [本文へ]

fig4.png


表1: 大陸・大洋別の年間の降水と蒸発の差。
大気の水蒸気収支による方法水文学的方法
地域 面積 1012 m2 FGGE ECMWF FGGE GFDL Bryan and Oort 1984年 Baumgartner and Reighel 1975年 Korzun 編 1974年 Lvovitch 1973年
アジア47 4.44.7*11.512.214.113.2
ヨーロッパ11 1.81.0*-2.02.83.03.1
アフリカ30 1.910.00.23.44.64.2
北アメリカ26 5.95.84.25.88.26.7
南アメリカ18 7.615.36.011.012.210.4
オーストラリア9 0.51.9-3.62.42.52.0
南極14 1.41.50.62.02.32.2
大陸155 23.640.36.539.747.041.7
北太平洋80 3.411.83.516.7@-9.7
南太平洋94 -1.3-9.36.9-0.9
北大西洋50 -9.9-17.4#-2.4-17.8@-31.7
南大西洋45 -8.2-15.4-10.5-18.7
インド洋78 -8.8-11.5-4.1-19.5-7.6
北極海8 1.31.4#0.42.1
海洋355 -23.6-40.3-6.5-39.7-47.0-41.7
水の体積としての総量、単位は1012 m3/年。 FGGEと書いた2列は1978年12月から1979年11月までの1年間の観測データによるもので、それぞれECMWFとGFDLの2か所で4次元同化の方法で作られたFGGE IIIbデータから筆者が計算したもの。他の列はこれまでに論文や本に発表されているもの。大陸、大洋の地域わけは論文ごとに少しずつ違うが、FGGEの2列はBaumgartner・Reichelのものにほぼ一致させてある。


2001-04-21, 2006-11-08
増田 耕一
[本文の最初へ]