Transcript 陸域炭素循環結合モデルの開発 (kato_060619 6843KB)

気候-陸域炭素循環結合モデルの開発
加藤知道 （かとうともみち）
独立行政法人 海洋研究開発機構
地球環境フロンティア研究センター
生態系変動予測研究プログラム
共生2連絡会議2006年6月１9日
これまでの展開
2005年度まで
・モデル結合を完成
・20世紀の炭素循環を再現
2006年度
・21世紀の炭素循環の推定
・SEIB-DGVMの導入
AGCM
(大気モデル)
ua, va, Ta, qa,
Ps, Prec, Rad↓
E, H,τ
Rad↑
CO2a
MATSIRO
(陸面の熱・水)
64
C a rb o n flu x (P g C yr -1 )
出力(完了)
NPP
62
ua, va, Ta, qa,
Ps, Prec, Rad↓
出力
NEP
LAI
Ts, Tg
CO2a
Sim-CYCLE
(陸域炭素)
導入
60
HR
58
4
C o u p led
NEP
2
0
-2
N E P - L U C e fflu x
U n co u p led
-4
1900
19 20
19 40
1960
19 80
2000
ye a r
20世紀の全球炭素動態
1900
2000
2100
21世紀のCO2濃度
動的植生モデル
(SEIB-DGVM)
本日の発表内容
・20世紀の気候-炭素循環の推定結果
・投稿準備中の結果について、紹介します
・21世紀ランへむけた結合モデルの改良について
・Sim-CYCLEとMATSIRO結合についての問題点
・Land use changeデータの導入に向けた取り組み(仮)
・ECRP陸域モデルGの北大での打ち合わせの内容
・2005年度FRCGC評価委員会の評価（Aだった）
→ ほとんどが陸域モデルについての(厳しい)コメント
→ ６/8に北大地環研にて対応策を検討
地球の炭素収支
750
放出源
化石燃料消費 &
セメント生産
5.5
土地利用変化
1.1
2190
吸収源
大気
3.2
海洋
2.0
陸域生態系？
1.4
(1980-1989, PgC yr-1)
Schimel(1995)
目的： 20世紀の陸域炭素循環を再現してみる！！
20世紀における陸域炭素フラックスの変化
64
NPP
C a rb o n flu x (P g C yr -1 )
62
・NPP、HRともに5-7％程度
上昇
60
HR
58
・NEPわずかだが正
4
→陸域が炭素シンク
NEP
2
・LUCを考慮すると、陸域は
たびたび炭素ソース
0
→LUCの影響大
-2
N E P - L U C e fflu x
-4
1900
19 20
19 40
1960
19 80
2000
ye a r
全球積算 NPP(純一次生産), HR(従属栄養生物呼吸),
NEP(生態系純生産) and NEP-LUCefflux(土地利用変
化による放出)
SST, CO2コントロール実験
66
64
a) NPP
・NPPについて、CO2controlはほぼ一定
だった
N o rm a l
C a rb o n flu x (P g C yr -1 )
62
S S T -c trl
60
→ NPP上昇はCO2
の影響のみを受けて
いる
58
C O 2 -c trl
56
66
64
b) HR
S S T -c trl
62
N o rm a l
・HRについて、SSTcontrolは中間的な
動き
60
→HR上昇はCO2と
SST上昇の二つの影
響を受けている
58
C O 2 -c trl
56
1900
1920
1940
1960
ye a r
1980
2 00 0
全球のCO2濃度
1月
7月
地表面付近CO2濃度 (ppmv; 1999）
地上観測点におけるCO2濃度変化の検証
38 0
a) BRW (N71.3, W156.6)
36 0
34 0
M o d el
O b se rv atio n
32 0
Atmospheric CO2 concentration (ppmv)
38 0
b) CBA (N55.2, W162.7)
36 0
・季節変化が北方では
大きく、南方では小さ
い傾向は一致
34 0
32 0
38 0
c) MLO (N19.5, W155.6)
36 0
34 0
・年々変化はかなり近
い動き
32 0
38 0
36 0
d) SMO (S14.3, W170.6)
→グリッドサイズ等を考
慮すると、それなりに
精度は高い
34 0
32 0
38 0
e) BHD (S41.4, W174.9)
36 0
34 0
32 0
38 0
f) SPO (S89.2, W24.8)
36 0
34 0
32 0
19 7 5
19 80
1985
Year
1990
19 95
2 0 00
輸送モデルによる大気CO2濃度の年々変化
Table. 全球炭素収支
(1980-1989; PgC/yr)
370
観測値 (0.25*spo+0.75*mlo)
CO2 濃度 (ppmv)
360
本研究 Houghton (2003)
350
モデル推定値
340
330
320
化石燃料
+5.4
+5.4
海洋
-2.2
-1.7
陸域生態系
-0.5
-2.4
土地利用変化
+0.6
+2.0
Total
+3.3
+3.3
310
1955
1960 1965
1970 1975 1980
1985 1990
1995
2000
year
全球平均大気CO2濃度
(ppmv; 1959-1999）
・Totalは実測値によく一致
・陸域は大気CO2に中立的な振舞
・しかし、その内訳は文献値と異なる
→ 文献値が過大評価の可能性
2
0 .4
a)
CO2 anomaly
1
0 .2
0
0
-1
-0 .2
・CO2アノマリーは温度
アノマリーから1-2年
遅れてピークがくる
→ Keelingらの
観測と一致
Temp. anomaly
-2
-0 .4
3
Carbon flux (Pg C yr-1)
Temperature (oC)
CO2 concentration (ppmv)
気候・炭素動態アノマリーの年々変化
2
・温度アノマリーとHRの
動きは近い
HR
b)
1
・CO2アノマリーの変化
速度はNCBの動きと
近い
0
-1
NPP
-2
NCB
-3
1960
1965
1970
1975
1980
Year
1985
1990
1995
→陸域が短期的な
CO2変動を起こす
CO2アノマリーに対する
気候・炭素動態アノマリー変化の先行
0.8
a)
Temperature
anomaly
0.4
Correlation coefficient
R (C O 2 fu ll_ N P P )
0
NiNO3
-0.4
CO2アノマリーから18-21ヶ月前に、
相関係数のピークがくる
-0.8
0 .8
b)
HR
0 .4
0
NPP
-0 .4
NCB
-0 .8
-4 0
-30
-2 0
-1 0
0
Lagged months from CO2 anomaly
10
まとめ
・本研究で開発した結合モデルは、20世紀中の炭素循環
をうまく再現することができ、その推定精度は高いと考
えられる
・20世紀中のNPP、HRともに5-7％程度上昇し、差し引き
のNEPはわずかに正だった。
・土地利用変化による炭素の放出の影響は大きい。
・20世紀中の気候変化において、陸域生態系は大気CO2
に対して、長期にはほぼ中立的で影響を及ぼしていな
いが、短期には主たる変動要因であった。
21世紀ランへむけた結合モデルの改良について
今年度の計画
・21世紀の炭素循環の推定
・SEIB-DGVMの導入
土地利用変化・植生分布の変化が
炭素循環と気候に及ぼす影響を調べる
→ アルベド変化・蒸発散量変化
→ MATSIROとのさらなる融合が必要
→ これまで各モデルで並列的に扱ってきた
陸面のアルベド・水収支・蒸散量推定を
統一する必要がある。
現在
AGCM
CO2 transport
ua, va, Ta, qa,
Ps, Prec, Rad↓
CO2a
E, H,τ
Rad↑
NEP
ua, va, Ta, qa,
Ps, Prec, Rad↓
CO2a
LAI
MATSIRO
Water
Balance
Photosynthesis
Ts, Tg
Sim-CYCLE
Water
Balance
Wg
Photosynthesis
最終形
AGCM
CO2 transport
ua, va, Ta, qa,
Ps, Prec, Rad↓
CO2a
E, H,τ
Rad↑
NEP
ua, va, Ta, qa,
Ps, Prec, Rad↓
CO2a
LAI
MATSIRO
Water
Balance
Photosynthesis
gc
Ts, Tg
Wg
Sim-CYCLE
Photosynthesis
MATSIROとSim-CYCLEの親和性を高める
～鉛直構造について～
1
2
0m
0.05m
0.25m
Sim-CYCLE
0-0.1m
1
3
1m
4
SEIB-DGVM
0m
0.3m
0m
0-0.5m
1
0.5m
0.6m
0.1-2m
0.5-1.5m
2
地温・土壌水分
2m
1.5m
1.5-3m
重み付け内挿
2
～
MATSIRO
3
Rooting depth
PFT typeによって決定
2.7m
3m
5
・各モデルの鉛直構造は変更しない
4m
・地温だけでなく土壌水分もMATSIROからもらうよう
にする
MATSIROとSim-CYCLEの親和性を高める
～水平構造について～
MATSIRO
Sim-CYCLE
PFTs: 14 types
PFTs: 20 types
融合させる必要があること
・アルベド
・気孔コンダクタンス
・植生変化への対応
Land use changeデータの導入に向けた取り組み(仮)
・21世紀中のLUCにおける気候と炭素循環への影響を調べるために、将来
のLUCのデータが必要になった。
・しかし、IMAGE groupはいつまでたってもデータをくれない。
・NIES groupのデータをもらうことを考える。
→ 環境研モデルの利用可能性
→ Jamstec研究開発促進アワードとの関連
LUCデータを作っている場所
・IMAGE group (Bas Eickhout;
RIVM/Netherland)
・SAGE group (Rammunkutty and
Foley; Univ. Wisconsin)
・NIES group (山形さんら; NIES)
ECRP陸域モデルGの北大での打ち合わせの内容(6/8)
○趣旨説明
○各人の成果概要と研究計画
○外部評価内容の検討
○グループとしての目標設定
○その他、コメント
発表者：伊藤、佐藤、稲冨、加藤
参加者：甲山先生、山中先生、高田先生、
もう一人研究者、大学院生
評価委員会のコメント要約
総評： 研 究 プ ロ グ ラ ム
生態系変動予測
総合評価
研究成果
研 究 課 題 ・目 標 設 定
社 会 へ の 還 元 ・イ ン パ ク ト
Ａ
Ａ
A
A
・昨年度の評価結果を受け止め、全体の体制が見直され、研究陣
の整備が進み、研究担当者の士気の向上、メンバー間での連携
の強化など、効果的な組織運営になってきた。
・3次元放射伝達モデルなど陸域生態系モデルの開発に大きな進
展があった。モデルのラインアップがそろい、全球モデルにつ
なぐ方向が見えて、H17年度に比べ格段の進歩がある。
・特に陸域生態系のモデルを全球の水・炭素収支などに結びつけ
るデータベースに関する戦略を立てる必要がある。
・この分野での日本唯一のチームとして、更なる人材確保に努め
るべき。
個別の点へのコメント１（厳しいコメント）
1.モデルの検証が不十分
・DGVMは、現在気候の植生分布を再現しただけでは
検証にはならない。
・各モデルの平均気候特性、気候変動特性、不確実性の幅に
関して、しっかり検証を行なうべき。
2.モデル間・研究者間の連携が不十分
・陸域は空間的非一様性が高いので、地域性に重点を置く
とともに、それと全球モデルとの関係を明らかすべきだ。
・まだなんでもありの印象が強く、目標設定をもっと絞った
ほうがよいのではないか。
・さまざまな試みがなされ成果を上げているが 、全球モデル
に何が必要なのかの検討もするべき。
・モデルを応用性の高いものにしていくために、実データの
予測をモデルの改良にフィードバックする必要がある
・既存サイト等を使い、フィールド研究者との連携を確実に行う。
・プログラム間での更なる連携が期待される。
個別の点へのコメント２（励ましのコメント）
3.今後の発展の方向性
・グローバルな気象変化の、人間の生活基盤への影響を評価するこ
とがFRCGC全体の目的だとすると、生態系変動の影響評価、と
くに陸域生態系のグループの研究は、もっとも重要な部門と
考えるべきで、より強化すべきである。
・特に日本ではこのような研究グループがまったくないため、人材
の確保と育成に十分留意いただきたい。
4.アウトリーチ活動のススメ
・国際的活動への積極参加が目立つ(C4MIP?)ことは評価できる。
・論文業績もあがりつつあるが、研究成果のすべてが必ずしも出版
されていない。まだ努力の余地がある。
・地域的な現実の問題に取り組むことによって社会へのインパクト
が高まることが期待される。
・外部へのインパクトを高める意味でも、ECRPの存在意義について
の議論を関連学会やその他の場で、もたれることが望ましい。
グループ構成
生態系変動予測研究プログラム（ECRP）
陸域生態系モデルグループ（GL：和田英太郎）
・物質循環サブグループ（SL：伊藤昭彦）
ポスドク研究員：加藤知道、稲冨素子
→ Sim-CYCLEの高度化
・植生動態サブグループ（SL：甲山隆司）
ポスドク研究員：佐藤永
→ SEIB-DGVMの開発
個々人の仕事をなるべく増やすことなく、
コメントに対応する答えを探す！
１．モデルの検証（データベースの利用）
・全球スケール水＆炭素収支
・環境変動への植生動態応答
問題点：
・観測とのスケールギャップ
・長期データの不足
・観測データの信頼性
GPPDI data (Zheng et al. 2003)
2. モデル間・研究者間の連携
Sim-CYCLEのスケール間の関わり
植物動態
＋
窒素循環機能
＋
それなりの対応関係はすでにある！
プログラム内・間連携は十分か？
プログラム内：
・海洋Ｇとはほとんど研究上の交流無し
・空間Ｇとリモセンデータのモデル利用に関して共
同研究が進められている
例：高精度のPARデータをモデルで利用
プログラム間：
・プログラム横断的に地球システムモデルを構築
・大気組成Ｇ研究者へバイオマス分布データを提供
・大気組成Ｇ研究者とN2O収支について議論
既存観測サイト・データとの比較
○岐阜高山サイト：土壌からのGHG放出（稲富）
○中国青海草原サイト：土壌からのGHG放出（加藤）、
PAR（小林）
○環境省Ｓ１：フラックス観測データ
○共生３：衛星観測データ
S-1
18 sites
Asiaflux
9 sites
CREST etc. 11 sites
Japan site total 38 sites
+
Koflux, Chinaflux (~20 sites)
東アジアの
CO2フラックスタワー
データ同化例：衛星データ利用
葉面積指数（LAI）
・モデル：炭素重からLAIへ変換（変換係
数：SLA, 比葉面積＝植生ごとの固定値）
Sim-CYCLE
・衛星（MODIS）：アルゴリズム
・ナジング：衛星LAIで置き換え
・同化：衛星に合うようSLAを較正
SLA
LAI
比較
LAI
MODIS
Data by P.J.Baruah (IIS, U.Tokyo)
モデルの観測データへのリアルタイム利用：
火災延焼と炭素放出の推定
野外火災の予報システム（JSTSORST）
モデル推定における不確実性の評価
空間分布・経年変動
3. 今後の発展の方向性
物質循環モデル
・陸域生態系の炭素-窒素-水循環を統合的に扱う新モデル
・アジア陸域における温室効果ガス収支マップの作成
・国内外の観測サイトにおけるデータの積極的利用
統合モデル
・共生課題における陸域部分の完成（炭素循環+植生動態）
・20-21世紀における土地利用変化の気候的影響評価
衛星データの利用
・高精度PAR、土地被覆データを用いた陸域炭素収支評価
・MODIS火災マップを用いたCO2放出の評価 (JST-SORST)
4. アウトリーチ
提案：生態学会でのシンポジウム/自由集会
「地球環境問題と陸域生態系モデル：可能性と課題」
温暖化をはじめとする地球環境問題は陸域生態系に深刻な影響を及ぼ
し、逆にその変化は気候変動へのフィードバック効果として作用する。
このような相互作用を解析し予測する上で、モデル手法は不可欠な
ツールとなりつつある。特に最近数年間で顕著な研究の展開がみられ
たが、それが生態学にもたらす新たな可能性については十分に伝えら
れていない。一方、複雑多様な生態系を予測するモデル開発には課題
も多く残されている。ここでは、国内の陸域生態系モデル研究の最新
成果を報告し、今後の展開について議論を行いたい。