Developing a new dynamic global vegetation model (DGVM) for global changing predictions

23
Developing a new dynamic global vegetation model (DGVM) for global changing predictions Hisashi SATO (FRSGC) Akihiko ITOH (FRSGC) Takashi KOHYAMA (Hokkaido Univ.)

description

Developing a new dynamic global vegetation model (DGVM) for global changing predictions. Hisashi SATO (FRSGC) Akihiko ITOH (FRSGC) Takashi KOHYAMA (Hokkaido Univ.). Toward developing the integrated land-surface-model. 原図䌊藀昭圊. Land surface physical process model. Land surface carbon - PowerPoint PPT Presentation

Transcript of Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Page 1: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Hisashi SATO (FRSGC)

Akihiko ITOH (FRSGC)

Takashi KOHYAMA (Hokkaido Univ.)

Page 2: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Toward developing the integrated land-surface-model

Vegetation dynamics model

Land surface physical process model

原図䌊藀昭圊

定着、死亡、攪乱、競争これらの結果ずしおの怍生倉動

Land surface carbon cycle model

「怍生動態モデル」ず「陞域炭玠埪環モデル」ずを分離しお蚭蚈するこずが困難であったため、この䞡者を結合したモデルを䜜成した。なお、怍生動態郚分の構造は、 LPJ-DGVM ず BIOME3 を参考にした。

Page 3: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

General description of the DGVM

Major advances from the previous DGVMs

(1) Individual Based Model               (except for herbaceous PFTs)

(2) Explicitly simulate spatial structures of vegetations

† Plant Functional Types

埓来のず比范した利点

†

空間構造を明瀺的に衚した個䜓ベヌスモデルであるため、パラメヌタヌ掚定に個䜓矀生態孊のデヌタ朚本密床、サむズ分垃、暹霢構成分垃をそのたた甚いるこずが出来る。たた、個䜓間競争の衚珟が劥圓であり、したがっお気候倉動に䌎った怍生倉動の速床を、より的確に予枬できるこずが期埅される。

Page 4: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

ケッペンの怍生地図ケッペンの怍生地図

幎平均気枩℃ 

幎降氎量

mm



Page 5: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Plant Characteristics (1)

Plant Functional Types (PFTs)

1. Tropical broad-leaved evergreen2. Tropical broad-leaved raingreen3. Temperate needle-leaved evergreen4. Temperate broad-leaved evergreen5. Temperate broad-leaved summergreen

䞋蚘 PFT1  8 が朚本、 PFT9 ず 10 が草本である。蚈算区画においおは、朚本 PFT ず草本 PFT ずが同時に分垃するず仮定した。朚本に関しおは同䞀区画内に耇数の PFT が混圚可胜であるが、草本 PFTs は同䞀区画内にはタむプしか分垃できないずした。

怍生景芳の求め方

䟋えば、雚緑暹 (PFT#2) ず C4 草本 (PFT#10) ずが優占しおいる区画においおは、朚本の個䜓密床が高ければ「雚緑暹林」ずし、朚本の個䜓密床が䜎く草本バむオマスが高ければ「サバンナ」、朚本の個䜓密床も草本バむオマスも䜎ければ「砂挠地垯」ずする。

6. Boreal needle-leaved evergreen 7. Boreal needle-leaved summergreen 8. Boreal broad-leaved summergreen 9. Temperate herbaceous (C3)10. Tropical herbaceous (C4)

各区画の怍生景芳は、優占する朚本 PFT の個䜓密床、及び優占する草本 PFT の単䜍面積圓バむオマスによっお決定する。

Page 6: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Crown : biomass, diameter, depth Stem : biomass, height, sapwood & heartwood DBH Root : biomass Stock-resource per individual tree

Woody PFTs

Herbaceous PFTs

Foliage : biomass in a unit areaRoot : biomass in a unit areaStock-resource in a unit area

草本は、葉ず根のみから構成され、それぞれは単䜍面積圓たりのバむオマスのみで扱われる。なお、草本は倚幎性であるず仮定し、根は森林火灜などの攪乱が生じない限り、葉が枯れおも残存するずした。

Plant Characteristics (2)

Individual base 。各朚本は、 Crown 、 Stem 、 Root の噚官から構成され、各噚官は次の倀を有する。 ( 右図も参照 )

Page 7: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

定着 (annualy computation)Woody PFTs

草本は定着過皋を明瀺的に扱わず、垞に存圚するず仮定した。たた、 C3 経路ずC4 経路の双方が䜿えるず仮定し、そのどちらかを䜿うかは前幎の気候条件に応じお決定させた。すなわち前幎の coldest-month temperature が 15.5 床 LPJ の基準以䞊ならば C4 経路、それ未満の堎合には C3 経路に切り替えられる。

Scenario 1: 指定した皮類の PFT のみが定着するパラメヌタヌ掚定甚モヌド。

Scenario 2: 珟圚どのような PFT が分垃しおいるのかに関わらず、その環境条件で定着可胜な党おの朚本 PFT が等確率で定着する infinite seed dispersal mode   。

Scenario 3: スピンアップでは Scenario 2 を適甚。その埌のシミュレヌションでは、各朚本 PFT のバむオマスに比䟋しお定着比率が決たる no seed dispersal mode   。 Herbaceous PFTs

各朚本 PFT には定着可胜な気枩範囲 maximum coldest-month temperature ず sum of growth degree-day が䞎えられおおり、その範囲に各グリッドの最近 20 幎間の平均気候が収たる堎合、その朚本 PFT が定着可胜であるずした。たた、前幎の降雚量が100mm を䞋回る堎合は、いかなる朚本 PFT も定着できないずした。

シミュレヌションは以䞋の぀のシナリオに基づいお行う。

Page 8: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

光環境の求め方 (daily & monthly computation)

Estimate light intensity on the top of the crown by using canopy location within the forest stand (SORTIE like)

To avoid ‘edge effect’, this scanning is performed among replicated forest stands, which surround the examining area.

Estimate light distribution within canopy using leaf area concentration and light attenuation index

Estimated light intensity

朚本の葉は均䞀に分垃しおいるこずを仮定し、地衚面に降り 泚ぐ光 量を求め、これを草本 PFT が利甚できるずした。

Grass layer

Woody PFTs

Herbaceous PFTs

たた高さ 3m 以䞋の朚本の光環境に぀いおも、同様の求め方をした。

Page 9: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Leaf phenology (daily computation)

各 PFT には萜葉性か垞緑性かの属性が䞎えられおおり、萜葉性の PFTs では以䞋の芏則によっお、展葉期ず䌑眠期が移り倉わる。その刀断には気枩ず土壌氎分に ぀いおの最近 1週間の running mean を甚いる。移り倉わり基準、すなわち「基準枩床」 ず「基準氎 ポテンシャル」は、 LPJ-DGVM より埗た。

●展葉期→䌑眠期

条件 気枩 < 基準枩床  .OR. 氎利甚 ポテンシャル < 基準氎 ポテンシャルむベント それたで展葉させおいた葉が党おリッタヌに組み蟌たれる。

●䌑眠期→展葉期

条件 気枩 > 基準枩床  .AND. 氎利甚 ポテンシャル > 基準氎 ポテンシャルむベント 貯蔵資源が利甚可胜ずなり、展葉が始たる

●䞊蚘アルゎリズムに぀いおの既知の問題 未解決  湿最で枩暖な環境䞋では、萜葉性 PFT が萜葉しなくなる。実際にはそんな事ない。 

Page 10: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

展葉  (daily computation)

Woody PFTs の最倧展葉量

(1) 暹冠サむズ、 (2)蟺材による氎 茞送量 * 、 (3) 根からの絊氎量 * 、による制玄のうち、最も匷い制玄が最倧展葉量を芏定する

Herbaceous PFTs の最倧展葉量

根からの絊氎量 * をちょうど満たす葉量

)GDD

GDD,0.1min(最倧展葉たでに必芁な

その日たでの   最倧展葉量展葉量

展葉量は、 GDD  Growing Degree Day, 日平均気枩が 5℃ を超えた日に぀いお、その超えた気枩を冬至の日より積算した倀の蓄積に埓い、以䞋の匏に埓っお増加する。

* 䞀週間の running mean

* 䞀週間の running mean

Page 11: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Woody PFTs の成長 (monthly computation)

1) 暹冠の枯れ䞊げ

盎前ヶ月の平均 PAR が補償点を䞋回る暹冠レむダヌを枯れ䞊げさせる。なお、䞀床枯れ䞊げた暹冠レむダヌは、たずえ光環境が改善されたずしおも、再び葉を぀けるこずは出来ないずした。

3) 貯蔵資源の補完4) 繁殖

2) 根の補完

この時点で残る利甚可胜資源の䞀埋が繁殖に甚いられる。これはリッタヌに組み入れられる。

5) 幹ず暹幹の成長

最倧展葉時に必芁ずなる根のバむオマスを、盎前ヶ月の平均氎利甚 ポテンシャルを甚いお求める。珟圚の根バむオマスが、これに満たなければ、䞍足分を補完する。

展葉量 貯蔵資源 量ずなるたで、貯蔵資源を補完

残りの利甚可胜資源は党お蟺材の成長にあおる。たた、蟺材の成長に䌎い暹高ず暹冠幅も成長させる。これらの算出には以䞋の関係匏を甚いる。

䜆し、朚本間で暹冠は重ならないずしおいるため、暹高、暹冠断面積共に、呚囲の朚本の配眮によっお最倧倀が芏定される堎合がある。

最倧暹高胞高盎埄定数

 最倧暹高暹高 exp1

1.6胞高盎埄 定数暹冠断面積

2.0胞高盎埄暹高 定数幹のバむオマス

Page 12: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

1) 暹冠の枯れ䞊げ

盎前ヶ月の平均 PAR が補償点を䞋回る暹冠レむダヌを枯れ䞊げさせる。なお、䞀床枯れ䞊げた暹冠レむダヌは、たずえ光環境が改善されたずしおも、再び葉を぀けるこずは出来ないずした。

3) 貯蔵資源の補完4) 繁殖

2) 根の補完

この時点で残る利甚可胜資源の党おが繁殖に甚いられる。これはリッタヌに組み入れられる。

5) 幹ず暹幹の成長

最倧展葉時に必芁ずなる根のバむオマスを、盎前䞀週間の平均氎利甚 ポテンシャルを甚いお求める。珟圚の根バむオマスが、これに満たなければ、䞍足分を補完する。

展葉量 貯蔵資源 量ずなるたで、貯蔵資源を補完

残りの利甚可胜資源は党お蟺材の成長にあおる。たた、蟺材の成長に䌎い暹高ず暹冠幅も成長させる。これらの算出には以䞋の関係匏を甚いる。

䜆し、朚本間で暹冠は重ならないずしおいるため、暹高、暹冠断面積共に、呚囲の朚本の配眮によっお最倧倀が芏定される堎合がある。

最倧暹高胞高盎埄定数

 最倧暹高暹高 exp1

1.6胞高盎埄 定数暹冠断面積

2.0胞高盎埄暹高 定数幹のバむオマス

Herbaceous PFTs の成長 (daily computation)

Page 13: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

æ°Ž 収支 (daily computation)

倧郚分は Sim-CYCLE ず同䞀巊図参照であり、この蟺りは、 MATSIRO ず結合埌には、詳现なモデルず倉わるはずである。

Sim-CYCLE 同様に、土壌は局ず仮定する。䜆し、 Sim-CYCLE では root uptake が䞋局からのみ生じるず仮定しおいるが、このDGVM では、 PFT 間の氎を 巡る競争を衚珟するために、 PFT毎に根バむオマスの垂盎分垃を定矩し、それに応じお利甚できる土壌レむダヌの比率を倉えた。

すなわち、ある PFT の根バむオマスの 80が䞊局にあるず定矩した堎合、このPFT の氎利甚 ポテンシャル (0.0  1.0) は、次の通りである。0.8× 䞊局含氎 率  + 0.2× 䞋局含氎 率

この蚈算の際、含氎 率含 氎量 / その土壌の最倧含 氎量には、 最近䞀週間のrunning mean を甚いた。

原図䌊藀昭圊

Page 14: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

攪乱 (annualy computation)

怍生の攪乱芁因ずしおは火灜のみを考えた。火灜モデルには、 Thonicke et al. (2001) が開発し、 Sitch et al. (2003) が LPJぞの組み蟌みに際しお簡略化したものを甚いた。

このモデルで火灜は、 fuel load バむオマスリッタヌが 200 g C/m2 以䞊蓄積しおいる堎合にのみ、 fuel load ずリッタヌ含氎 率によっお求められる確率によっお発生する。

火灜の生じた区画では、焌死した朚本の党バむオマス、生き残った朚本の葉矀バむオマス、草本の党バむオマス、そしお党リッタヌが CO2 ずしお攟出されおしたうず仮定した。

䜆し、䞀般的には火灜跡には倚くの炭化朚片が残り、これらは生分解されにくいため炭玠ストックずしお長く地䞭に留たる。たた、それらの炭化朚片はアルベドを䜎䞋させる効果を持぀ずも考えられるが、これらの効果に぀いおは、このモデルには含たれおいない。

火灜時における朚本の生存率には、 PFT毎に特有の倀が䞎えられおいる。

䜆し、 Sitch et al. (2003) は、これらの倀の出所を匕甚しおいないので、どの皋床信頌のおける倀であるかは䞍明である。

Page 15: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

死亡 (annualy computation)

死亡率 = Max [1.0, mort_greff + mort_heat + mort_limit]

mort_greff Background mortality from current year’s growth efficiency

mort_heat

mort_limit If 20yrs running mean value of bioclimatic valiables of air-temperature fall outside a PFT limits for survival

最倧展葉時の葉面積量幎間のバむオマス増加

const2

const1

1

300

PFT,0max

,1min day

毎の基準枩床その日の気枩

( 成長効率 )

(熱ストレス )

( 生存可胜な気枩範囲 )

Otherwise 0.0

0.1

Page 16: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

光合成

呌吞

Turnover

Sim-CYCLE のルヌチンを、そのたた䜿甚

維持呌吞ず成長呌吞を考える。基本的には Sim-CYCLE ず同じだが、維持呌吞に芁する資源量に利甚可胜な資源が満たない時には、たず貯蔵資源が甚いられ、次に葉ず根のバむオマスが甚いられ、それでも足りない堎合には死亡するずした。

LPJ のアルゎリズムずパラメヌタヌを、そのたた䜿甚

土壌呌吞 Sim-CYCLE のルヌチンを、そのたた䜿甚する予定 実はただ組み蟌んでいない

(daily computation)

(daily computation)

(daily computation)

(daily computation)

Page 17: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

幎間の動態デモンストレヌション熊本県氎 俣垂

前半 100 幎間この堎所で定着可胜な党おの PFT が等確率で定着埌半 100 幎間 PFT毎のバむオマス比に応じお定着枩垯性萜葉広葉暹ず枩垯性垞緑広葉暹の混合林 → 枩垯性萜葉広葉暹林

※諞

パラ

メヌ

タヌ

は未調敎

Page 18: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

幎間の動態デモンストレヌション Pasoh, Malaysia

前半 100 幎間この堎所で定着可胜な党おの PFT が等確率で定着埌半 100 幎間 PFT毎のバむオマス比に応じお定着

熱垯性垞緑広葉暹林

※諞

パラ

メヌ

タヌ

は未調敎

Page 19: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

モデル党䜓の入出力

気枩地枩 ( 地衚 , 地䞋 0.1m ず2m)

降雚量湿床雲 量颚速CO2濃床

怍生 status 優占皮、陞䞊バむオマス、土壌炭玠含 量

二酞化炭玠収支 光合成、呌吞、分解、火灜

æ°Ž 収支  runoff 、蒞散、蒞発

攟射収支 アルベド

党お Sim-CYCLE ず同じDaily time-step Daily time-step

Yealy time-step

Page 20: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

炭玠埪環の抂略

点線内陞域生態系点線倖倧気

Page 21: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Simulation will be conducted on the T42 global grid (128×64), each of which includes 10 replication forest stands.

Thus, assuming 1/3 of the earth surface is terrene, about 27000 independent forest stand will be independently simulated.

To date, this would be the most complex ecosystem model that have ever made.

Simulation procedure

小サむズの林分を耇数シミュレヌトさせる䞻な理由ずしおは、攪乱の問題があげられる。䟋えば寒垯林で頻発する森林火灜は、䞀床生じるず、シミュレヌトしおいる林分の倧きさが 30×30だろうが 1haだろうが、その殆ど党おが壊滅しおしたう。このように機䌚的に倧きく倉動する単䞀の林分をもっお、グリッドの代衚倀ずさせるこずは適圓ではない。

Page 22: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

Estimate parameters and algorithm of a tree growthso that tree-form are reasonably simulated for each PFT

Estimate dynamics parameters (Establishment, Mortality, Disturbance):so that density and age distribution of tree are reasonably simulated when only one PFT composes the forest

Estimate metabolic parameters (Photosynthesis, Respiration, Allocation):So that GPP, biomass, LAI, and distribution of DBH are reasonably simulated. This will be conducted on forest that was composed of only one PFT.

By repeating above (2) and (3), convergence parameters

Conduct test run on global gridthen examine that distribution of vegetation and GPP at equilibrium are reasonably simulated.

Procedure for parameter estimation and tuning

(1)

(2)

(3)

(4)

(5)

Page 23: Developing a new dynamic global vegetation model (DGVM) for global changing predictions

点蚈算版のコヌドに぀いおは、昚幎末にいちおう曞き䞊がった。珟圚はコヌドチェックず仕様曞論文マテメ゜の叩き台を䜜成しおいる。その埌、ベクタラむズを行う。

地点蚈算版モデルを甚いお、パラメヌタヌ掚定ずチュヌニングを行う (on SX5) 。

党球蚈算ができるようにコヌドを盎し、党球グリッドにおけるオフラむン実隓 (on the Earth-simulator) 、および解析。

以䞊を、なるべく新幎床明け早々たでには終えお、その埌、早めに論文執筆ぞ移りたい。

 今埌の䜜業手順

月䞋旬

ヶ月

ヶ月