植生における水循環

降水が植生を介して大気や河川に移動する順番で解説します

（1）遮断蒸発

（2）浸透・流出・土壌浸食

（3）水ポテンシャルと根系における吸水

（4）通導組織による水の移動

（5）気孔コンダクタンスと蒸散

植生における水循環

主なネタ本

はじめに： 地球システム・生態システムにおける水循環の重要性

（1） 水蒸気と雲は、主要な温室効果物質である。

（2） 地表面から大気への乱流によるエネルギー輸送（潜熱熱＋顕熱）

において、その約77％を潜熱による輸送が占める。

（3） 植生分布は降水量の多寡により強く規定されるなど、水循環は生 物地球化学過程も大きく制御する。

（4） 現在、人は普通にアクセスできる淡水のうちの半分を利用している

地球放射の大 気中の減衰

遮断蒸発（Canopy Interception）

樹冠には、典型的には、針葉樹林では降水の15％程度 を蓄え、広葉樹林では降水の5～10%を蓄えることが できる。

遮断蒸発の割合は、針葉樹・広葉樹の別の他に、葉 量・樹形・樹皮テクスチャ・樹冠構造・雨の降り方な どに応じて異なる。

例えば、樹冠構造が入り組んで、乱流の発生しやすい 植生では、一般に遮断蒸発量が高くなるが、このよう な条件における遮断蒸発速度は、正味放射量や飽差と 強く相関せず、また昼夜の差も小さい。

樹冠に遮断される降雪の割合は、降雨のそれの倍にも 達することがしばしばである。ツンドラのように雪面 が常に直射に晒される場所、または大陸性の亜寒帯林 のように降水量と風速の低い場所では、それぞれ約 30%と約50%の冬期降水が昇華によって失われる。

実測例：様々なユーカリ属における葉面積と樹冠 に貯蔵可能な水量との関係

図：加藤知道監訳「生態系生態学第2版」森北出版

遮断蒸発の扱い例：

LPJ-DGVMモデル(Gertenら2004による水文過程改良版)の場合

E_I = min(α*E_q, I_sc)

E_I

遮断による一日の蒸発量

係数

E_q Priestley-Taylor

式により算出される平衡蒸発散量

遮断蒸発貯留量（

）

I_sc = P × LAI × i × f_v

P

一日の降水量 （雨と雪は区別しない

）

葉面積指数

i

植生タイプごとのパラメーター（右の表を参照）

f_v

植被率（フェノロジーの効果込み

）

Priestley-Taylor式から算出される水制約の無い時の蒸発量 樹冠に貯留できる最大水量。これから溢れた降水のみが地表に到達する

針葉と広葉とで3倍もの 差があるが、これは Biomeごとの雨の降り 方の違いの効果も暗に 表現しているため

Gerten et al. (2004) J. Hydrol. 286

「森林伐採に伴う土壌流出は、人類の歴史において繰り返し文明を崩壊させてきた」

ジャレド・ダイアモンド“文明崩壊 滅亡と存続の命運を分けるもの”より

地表面の状態は、地中への水の入りやすさに影響する

浸透できなかった降水の一部は地表流出し、しばしば土壌浸食を起こす

地表面から地中までの水の流れと土壌浸食（Infiltration）

米国南部：不適切な耕作が土壌浸食をもたらした例

出典：片岡夏実訳「土の文明史」築地書館

植被は、地表面における保水力維 持と、土壌への有機物の投入を通 じて、一般に降水の土壌への浸透 を促進する。

土壌生物による土壌の団粒構造や パイプ構造の形成を促進すること で、土壌の透水性と保水力を高く する。

森林は、遮断蒸発と蒸散のために、河川流出量を減少させる。

しかし、森林は、降水時に地表流となって一時に流出する水量を抑え、土壌に保水した水 を中間流や地下水流として徐々に河川に流出させるために、河川流量を平滑化させる。

図の出典：藤森隆郎「森林生態学」全国林業改良普及協会

時間

図の出典：塚本良則「森林水文学」文永堂出版

日本中部と同緯度のアメリカ南部アパラチャ山 系にあるCoweeta試験地で行われた流域実験の 結果。7年間の平均値。5月初旬から10月中旬 までが蒸散期間。降雨は年間を通じてほぼ一定。

実際のハイドログラフ

蒸散が最大の時期と流出量 が増加する時期との間に明 らかなギャップがある

森林の下層植生バイオマス・リター 量と土壌への浸透能力との関係

間伐の遅れたヒノキ林では、土壌の浸透能力が大幅に低下 し、土壌流出量が高くなる。

土地利用の放棄が土壌浸食を生じさせるという、

過去に例を見ないタイプの土壌流出が現在日本で進行している

様々な森林における土壌浸食量

地表面が下層植生や枯死物で厚く覆われて いるほど、土壌への浸透能力は高まる

ヒノキ林ではリター層 が著しく発達しにくい

余談：放棄人工林における土壌流出

放棄された人工林がもたらす問題

保水力低下と 土壌流出

台風や豪雨による 崩壊リスク

放棄された人工林

東海豪雨による森林崩 壊（沢抜け）の例

間伐している人工林

↑ ←写真と図作成：蔵治光一郎教授

写真：台風21・22号土木学会調査団

日本の森林は40％以上が人工林であるが、間伐が間伐対象林齢の森 林の半数程度しか実施されていない。放棄人工林は、いずれ広葉樹 から成る自然林に遷移するが、その過程で公益上の諸問題が生じる

CREST課題 「荒廃人工林の管理による流量増加と河川環境の 改善を図る革新的な技術の開発」の成果

様々な野外操作実験により、林床植生やリターの量と、降水の土壌浸透能

力・土壌流出量との定量的関係が明らかになっている。

強度の間伐を行う事が放置人工林対策になる

強度間伐後の下層植生の回復 間伐後の林床照度変化

強度の間伐により林床が明るくなり、下層植生が回復する。これに伴って、林床の最大浸

透速度が増大し（保水力の増加）、土壌流出速度が低下する（土壌流出の阻止）。そのよ

うな森林は徐々に温帯広葉樹林へ戻り、やがてメンテナンスフリーとなると期待される

土壌流出の広域シミュレーション

Revised USLE* (RUSLE) に よる土壌移動速度

= R × L × S × K × C × P

R

：表層流出量の関数

：斜面長の関数

S

：斜面傾斜の関数

K

：土壌の流れやすさの関数

：植被による関数

P

：管理の関数

*土壌流出量の代表的な経験モデルUSLE（the Universal Soil Loss Equation）の改訂版

出典：佐藤ら (2018) 日本土壌肥料学雑誌89 USLEを組み込んだVISITモデルによる、現在の土壌流出量の分布

出典：Ito (2007) GRL34 気候変動

＋

土地利用変化

気候変動のみ

（MIRIC@A1Bシ ナリオ）

土地利用変化のみ

21世紀中の変化のシミュレーション

RUSLEを使用した全球シミュレーション例：

表土流出に伴う炭素流出量＠現在気候

アジアモンスーン域のように総水量の増大が 予測される地域にて、顕著な増大が予測され ている

温帯域においては、現在も生じている植生回 復トレンドが持続し、それが減少トレンドを 生じさせている？

表の出典: 佐藤・安成(2014) 渡邊ほか編「臨床環境学」2章、名古屋大学出版会

値の出典: 久馬(2005) 「土とはなんだろうか？」、 京都大学学術出版会

ある推定によれば、世界の陸地に おいて1 haあたり土壌生成速度は 平均570

これは厚さにすると、概ね0.057

mm/年

これを遙かに上回る速度で、土壌 の流出が生じている。その主な理 由は、人の土地利用と考えられる。

現在、人の土地利用が土壌流出 量を大幅に増大させている

*

土壌ー植物ー大気間の水の流れ

水は、水ポテンシャル(ψt)が、高い方から低い方へ向かって移動 する

Ψt = Ψp + Ψm + Ψo + Ψg

Ψt : 合計水ポテンシャル。地面直上に置いた純水をΨt=0と定義する Ψp : 水力（Hydrostatic pressure）ポテンシャル

Ψm : 基質（Matrix）ポテンシャル Ψo : 浸透（Osmetic）ポテンシャル Ψg : 重力（Gravity）ポテンシャル

土壌ー植物ー大気間の水 ポテンシャル勾配

（図中の値は、おおよその範囲）

図の出典：森林立地学会編「森のバランス」

土壌から葉への吸水速度＝土壌水ポテンシャル－葉の水ポテンシャルー重力ポテンシャル 通水抵抗

基質ポテンシャルと重力ポテンシャルで決まる。

-6kPaより弱い張力で保持されている土壌水は、重力ポテンシャルで速や かに排水される。-150kPaの毛管連絡切断点より強い張力で保持されてい る土壌水（難有効水）は、植物は根の周辺のみから吸水できる

細胞質の溶質のモル濃度で決まる浸透ポテンシャル（マイナス値）と、細胞 壁が水を押し出そうとする圧ポテンシャル（プラス値）の和で決定される

重力ポテンシャル = 樹高(m) × 重力加速度(約9.8 N/Kg)

＝ 樹高(m) ×9.8 kPa

図の出典：

森林立地学会編「森のバランス」

葉の含水率と水ポテンシャル

水ポテンシャルが0のときの細胞 内の水分量を1とした時の相対値 圧ポテンシャルが0となる相対含水率よりもさら

に含水率が低下すると原形質分離が起きる

水ポテンシャル（実線）は、

圧ポテンシャルの分だけ、浸 透ポテンシャル（破線）より も高い（吸水力が低い）

写真の出典：Wikipedia”原形質分離”

図の作成：大槻恭一

吸引圧

の常用対数

吸引圧

33ｃm H₂O 1000ｃm H₂O 15900cm H₂O

植物の生長や水の移動において 重要な土壌状態を示す値

土壌含水率の低下とともに

基質ポテンシャルは高まる

水分定数は土壌の種類に強く依存する

図の作成：大槻恭一

土壌の種類ごとの圃場容水量

と萎れ点（

）

土壌含水率 小 大

圃場容水量：水頭が-0.033MPa （約-0.33mH₂O，pF=2.5） 永久萎れ点：-1.5MPa （約-150mH₂O，pF=4.2）

表の作成：大槻恭一

水ポテンシャルは様々な単位で表される

吸引圧(cm)は、符号を入れ替えてhPaと読み替えても、大体合う。

ので、ザックリと、－10^(g)hPa

表の作成：大槻恭一

空気侵入ポテンシャ ル（ψ_mの最大値。

Ψ_mがこの値の時に、

土壌空隙は全て水 で満たされている）

Ψ_m=-33cmとなる体積 容水率(圃場容水量)

Ψ_m=-1500cmとなる体積容水 率 (成長阻害水分点) 水力コン

ダクタンス

体積容水率

とマトリックスポテ

ンシャル

との関係

ところで、先の表ではψ_mが-33kPaのシルト質土の含水率θ_-33は0.33m³m^-3であり上の 結果と一致しない。表中の値は、土性別の土壌について測定された平均値であり、パラ メーターの平均値を用いて推定した含水率とは必ずしも一致しない。

作成：大槻恭一

先の表のシルト質壌土のパラメーターを用いて、マトリックスポテンシャルψ_mが-33kPa と-10kPaにおける体積含水率θを求めよ。なお、飽和含水率θ_sは0.5m³m^-3とする。

例題

手順1）

θとψ_mの関係式を変形 手順3）手順1の式に、各パラメーター値を代入

手順2）シルト質壌土のψ_eの単位を変換

−0.21 m

×

土壌に入った水は、圃場容水量 (Field capacity; 重力と基質ポテンシャルが釣り合った状態) に達するまで、または母岩や凍結層に達するまで、下方に移動する。

その際の移動速度は、 水力コンダクタンス（Hydraulic conductance）に規定される。

K_sat ：飽和時の水力コンダクタンス Z ：深さ（深いほど大きな正の値）

F ：正の値を持つ定数

粒子の粗い土壌ほど、圃場容水量は低く、また水力コンダクタンスは高い。また、水力コン ダクタンスは、土壌深度と共に急速に低下する（下式のようにモデル化される）。

Js = ΔΨt ×K / l

l : 経路長

K : 水力conductance

土壌中の水の移動速度（Js）は、水ポテン シャルの勾配と、水力コンダクタンスで決まる

また、水力コンダクタンスは、土壌の乾燥と共に急速に小さくなる。これはmatrix potentialが強くなるのと、土壌中の水の流路が大幅に伸びるため。

土壌には通水抵抗があるため水ポテンシャル勾配の解消には時間がかかる

シベリアのカラマツにおける、吸水速度と蒸散速 度の時系列変化。明け方は、幹に蓄えた水で蒸散 を行い、失われた水は夜間にrechargeされる。

木本の幹には、多くの場合、5～10日分の蒸散量 に相当する水が含まれている。

植物体には通水抵抗があるため水ポテンシャル勾配の解消には時間がかかる

根の中の通水経路には、細胞壁や細胞間 隙を通るアポプラスト経路と細胞内を通 るシンプラスト経路があるが、根内には 菌類等の侵入を抑えるカスパリー線があ るため、アポプラスト経路であっても1 度は細胞内を通過する

図：加藤知道監訳「生態系生態学第2版」森北出版

図の出典：Lee et al. (2005) PNAS 102

Hydraulic Redistribution

根系を介した土壌水の再配分のこと。土壌深部 から土壌浅部への、土壌水の持ち上げのみ現す 際にはHydraulic Liftとも呼称される。

これらは気孔が閉じているとき、根系が2形性

（例えば、直根と側根の組み合わせなど）、土 壌粒子が細かくて水が浸透しにくい、などの条 件下で発生しやすい。

Hydraulic liftで再配分される土壌水量は、蒸散

土壌浅部は深部よりも栄養塩量も土壌生物量も 多いため、Hydraulic liftは生態系全体の生産性 を増加させることに寄与するようである。

直根

側根 降雨があったタイミング

図の出典：Lee et al. (2005) PNAS 102(49)

アマゾン盆地における影響評価 蒸発散速度 地表気温

全球における影響結果

蒸発散速度 地表気温

年平均緯度帯ごと 季節変化

Hydraulic Redistributionは、全球水循環に大きな影響を与えている 可能性が指摘されている

大気循環モデルCAM2と陸面過程モデルCLMとの結合モデルを用いた影響評価シミュレーション

しかし、このメカニズムを取り込んだ動的全

球植生モデルは、多分存在しない

根の鉛直分布

生活型ごと Biomeごと

隣接した草地と灌木地帯間に おける、水ポテンシャルの鉛 直分布。草本よりも木は、より 深くに根を伸ばすため。土壌 の深いところまで乾燥させる

草本よりも木本、木本よりも灌木で、深い根を持つ傾向あり。ただ し、同じ生活型やBiomeであっても、膨大な多様性がある

広域植生モデルにおいては、草本と木本間のみに鉛直方向の水

競争を仮定することが一般的。

Root Depth

に関するメタ解析の結果

図の出典：Fan et al. (2017) PNAS 114 (40)

A.

年降水量

B.

土壌テクスチャ

土壌の厚さ

D.

植物の生活型

植物の分類群

地下水位

影響あり

顕著な影響あり

2200以上の観察データを用い たメタ解析の結果。メタ解析 とは、複数の研究結果を統合 した、より総合的な見地から 分析のこと。

図の出典：Fan et al. (2017) PNAS 114 (40) 種群に関わらず、Root depthは、年降水量ではなく、地下 水位によって制御されていることが分かった

最も顕著な影響を与えた地下水位について、更なる解析を行 うために、種群（genera）ごとに、root depthと年降水量と の相関、 root depthと地下水位深との相関を、それぞれプ ロットした。

根系の深さは、環境条件に応じた表現型可塑性の 程度がとても大きい。

幾つかのDGVMでは、土壌乾燥度に応じて、地上 部と地下部のNPP配分比率を変えてやることで、

この表現型可塑性を表現している。

図の出典：Fan et al. (2017) PNAS 114 (40)

先のメタ解析とConsistent な観測事例。同じ気候区に 生育している同じ種であっ ても、地下水面の深さに応 じて根の深さが異なる

シベリアニレ

Plains cottonwood (ポプラの一種)

地下水位 の位置

図の出典：Fan et al. (2017) PNAS 114 (40)

地下水位が、ある程度浅い（5ｍ程度まで）場合に、

根深と地下水位の関係は、最も顕著に生じる

Fisher et al. (2018) GCB 24

←パッチ間で異なる水資 源プールを持つ

（ED2, LPJ-Guessなど）

水資源を巡る競争の植生モデルにおける表現

←半乾燥帯などでは、木本は広い根系を持ち、広 い面積から土壌水を集めることがあるが、このよ うな機構を扱う広域植生モデルは無い

←全てのパッチが同一の 水資源プールを共有

（CLM, SEIBなど）

実際の状況は、これらの中間。

パッチやギャップの大きさ、

根系の水平分布に応じて、適 切なサイズが決まるはず

他方、殆どの植生モデルが土壌水の鉛直方向の競

争は考慮する（でも、きちんと検証してない）

幹における水の流れ

樹木は根から水分と養分を吸収し、幹を上昇して樹冠部に達し、葉で太陽 光線を受けて光合成を行う。この光合成産物が木部と樹皮部の間の形成 層に送られ、ここで細胞分裂が生じ、樹木は成長する。

初め未着色である年輪はある程度古くなる と着色される。未着色の部分を辺材といい、

着色された部分を心材という。

辺材の細胞は生きていて、水分を通したり 養分を貯蔵したりしているが、心材の細胞 は既に死んでおり、幹の機械強度を保つ 以外の機能は喪失している。心材には、耐 朽性を保つための物質が蓄積されている。

きれいな正相関を持つ。

傾きは、種や環境によって異なる。

北米に分布する３つの針葉樹 種において観察された、辺材 断面積と葉面積との間の相関

心材と辺材

幹内の通導組織の水中は、水の凝集力と壁面への付着力に よって維持される。

よって、水コンダクタンス（通導性）とキャビテーション

（通導組織の水中が切れてしまうこと）のリスクとの間に は、正のトレードオフが生じる。

寒冷地の木本は、凍結によるキャビテーションのリスクに 晒されるため、細い導管を多数持つ種が多く、また導管の refillが出来る種もある。

←

キャビテーションが生じる植物体の水ポテンシャル

（横軸）と、生育地の土壌における最低水ポテンシャル

（縦軸）との関係。正の相関がある。

通導組織の導水性とキャビテーション

気孔と蒸散

気孔は広葉の場合は葉の裏側に分布する。葉面の気孔以 外の部位は、クチクラという不透水層で覆われている。

したがって、顕熱の交換は葉の両面で生じるが、樹幹が 乾燥している時の潜熱の交換は、ほぼ葉の裏面でのみ生 じる。

図：熊谷朝臣

気孔の数は、1 cm²あたり約一万個。気孔面積の葉面積 に対する比率は0.3～3%を占めるにすぎない。水の蒸散 は、この面積比に比例するわけでは無く、一種のオアシ ス効果（水面の点在効果）により、葉面全体が濡れてい る時の10%程度の効率で、蒸発が生じる。

植物が気孔を開く理由はCO₂の取り込み。気孔からの水 蒸気の放出は、一般にコストと考えられている* 。従っ て、その開口度は、植物の水ポテンシャルに応じて制御 される。

葉肉細胞間隙。その葉内におけ る体積割合は、種や葉によって 異なり20～80％

*葉を冷やす事で、葉の温度を光合成に適した範囲に収める機能を持つとの指摘もあるようです

気孔コンダクタンス

気孔のガスの通りやすさを気孔コンダクタンスと呼ぶ。

気孔コンダクタンスの逆数が気孔抵抗である

G_c

と

がゴッチャに扱われている気がするので 注意。水蒸気は

よりも分子量が小さく、より 拡散しやすいので、以下の関係が成り立つ。

G_h = 1.56×G_c

光合成速度

と気孔コンダクタンス

の間には 次の関係が成り立つ。

A = G_c × (C_a - C_i)

また、蒸散速度

と気孔コンダクタンス

の間 には次の関係が成り立つ。

E = G_h × (H_i - H_s)

記号 説明 ^{よく使われる単位}

G_c 気孔コンダクタンス（CO₂） mol CO₂ m^-2 s^-1 G_h 気孔コンダクタンス(水蒸気) mol H₂Om^-2 s^-1

C_s 葉表面のCO₂濃度 μmol mol^-1 C_i 細胞間隙のCO₂濃度 μmol mol^-1 H_s 葉表面における水蒸気濃度 mmol mol^-1 H_i 細胞間隙の水蒸気濃度 mmol mol^-1

A 光合成速度 μmol CO₂ m^-2 s^-1

E 蒸散速度 mmol CO₂ m^-2 s^-1

葉内は湿度100%と仮定できるので、H_iは葉温から算出できる。

CO₂濃度は、

ppmvがそのまま使える

大気

200

1分子のCO

を取り込むために、例えば 200分子

もの水分子が、葉から流出する

1

*以下の環境条件下の概算値 気温30℃、相対湿度60%、

大気中CO₂濃度400PPM、

気孔内CO₂分圧が大気の70%、

水蒸気とCO₂の拡散係数比は1.56とした

葉内

実習：

上の200分子の根拠を示せ。なお、30℃

の飽和水蒸気圧は約40hPa、大気圧は、

海抜0m付近で約1000hPaである

気孔コンダクタンスについての経験的知見

葉面からの蒸散速度は、飽差と風速、そして気孔コンダクタンスにより制御されている。

気孔コンダクタンスには、次の経験的知見が得られている。

(1) 暗黒下ではほぼゼロ。光強度と共に上昇し、やがて飽和する。

(2) 空気が乾燥すると下がる。

(3) CO

濃度が上がると下がる。

(4) 土壌水分が不足すると下がる。

(5) 葉の水ポテンシャルが下がると下がる。

(6) 夜がくると（明るくても）下がる種もいる。

このように気孔コンダクタンスの反応は複雑で、制御の生理メカニズムも

完全には分かっていない事もあり、半経験的モデルにより扱われる

現在広く利用される気孔コンダクタンスモデルの構造

) / 1

( ) CO

( _{2 3}

2

min k k

g

g 大気中 濃度 飽差

光合成速度







 



濃度 大気中

相対湿度 光合成速度

2 1

min CO

 

 g k g

Jarvis

型

BWB

型

Leuning

型

g 気孔コンダクタンス

g_max 最大の気孔コンダクタンス g_min 最低の気孔コンダクタンス k_n 経験的な定数

Γ CO₂補償点 大気の乾燥度を表す変数として、相対湿度と飽

差の2種類が混在している。 Jarvis型モデルは、光合成速度を必要としない（観測した気象 要素とCO₂濃度のみの関数）ため、観測系の研究者に人気

BWB型モデルは、光合成速度の関数とすることで、

より実際の生理機構を反映させた。これはシンプル ながらも実測値と良く合うために、多用されている。

Leuning型モデルは、相対湿度の代わりに飽差を用いる。これは、実際の気

孔コンダクタンスが相対湿度では無くて飽差に依存しているため。ついでに、

CO₂補償点を考慮に入れ、低CO₂濃度下の挙動を正確に表現した。

g = g_max × f₁(光強度) × f₂(飽差) × f₃(大気中CO₂濃度) × f₄(気温) × ‥

Big leafモデル

Ball型モデルの検証例

出典：Ball et al. (1987) 様々な環境条件下における光合成速度(A、

横軸)と気孔コンダクタンス(g、縦軸)との関係

h_s : 相対湿度 c_s : CO₂濃度 Ball型関数の出力（横軸）と

気孔コンダクタンス（縦軸）との関係

大気乾燥度の変化予測*は、飽差と相対湿度とで大きく異なる

20世紀末

21世紀末

20世紀末

21世紀末

相対湿度 飽差

東経20度 南北縦断ライン

21 世 紀 末 20 世 紀 末

東経20度 南北縦断ライン

21 世 紀 末 20 世 紀 末

* IPCC-AR4A1B CO2排出シナリオの元におけるMIROC-AGCM出力 † 最後の10年間の平均値

このような、「飽差が増大しても、相対湿度は殆ど変化しな

い」という予測は、多くの

から得られている

42.4 hPa

@30℃ RH 70%

@35℃ 相対湿度 70%

56.2 hPa

大気中のCO₂分子数が変わらなけれ ば、VPDは上昇、RHは下降、しかし 変化の割合は前者の方が大

VP=

29.7 VPD=

12.7

VP=

39.3 飽差=

16.9 飽和水蒸気圧は気温と共に上昇

飽差と相対湿度の挙動が異なる理由の一つ

飽和水蒸気圧

飽和水蒸気圧

気温

29.7 VPD=

26.5

@35℃ RH 52%

実際には、高いVPDは高い蒸発散効 率をもたらすので、RHは比較的安 定する（しかしVPDの拡大分は全て 埋まらない）。

短縮形一覧（これら は業界で一般的です）

VP:水蒸気圧

（Vapor Pressure） VPD：飽差 (Vapor Pressure Deficit) RH ：相対湿度 (Relative Humidity)

相対湿度を用いた気孔コンダクタン スモデルは、現在においても多くの 気候モデルや地球システムモデルで 利用されている

Land Surface Model GCM Nation

CLM4 NCAR-CCSM4.0 USA

CLM3 NCAR-CCSM3.0 USA

LSM NCAR-PCM1 USA

Friend and Kiang GISS-ER USA

LM2 GFDL-CM2.0 USA

LM2 GFDL-CM2.1 USA

MOSES-2 UKMO-HadGEM1 UK

MOSES-2 UKMO-HadGEM3 UK

CTEM-CLASS2.7 CCCMA-CGCM3.1 Canada

BCC-AVIM BCC-CM1 China

CLM3 FGOALS-s2 China

ISBA BCCR-BCM2.0 Norway

ISBA CNRM-CM3 France

ORCHIDEE IPSL-CM4 France

JSBACH/BETHY ECHAM5/MPI-OM Germany Roeckner et al. MIUB-ECHO-G Germany/K

orea Based on SiB MRI-CGCM2.3.2 Japan MATSIRO MIROC3.2(hires) Japan MATSIRO MIROC3.2(medres) Japan Gordon et al. CSIRO-MK3.0 Australia Gordon et al. CSIRO-MK3.5 Australia

Sim-CYCLE(VISIT) Japan

SEIB-DGVM Japan

ED USA

NOAH-GEM USA

IBIS USA

IBIS-2 USA

SiB2.5 (SiB3?) USA

VISA USA

作成： 高橋厚裕

地球システムモデル・気候モデル (GCM・ESM)陸面過程モデル (LSM)

飽差

(Leuning type model)

相対湿度

(BWB type model)

Stomata conductance is a function of

SEIB-DGVMを用いた感度分析＠アフリカ大陸

今世紀末にかけて、全体に植物生産力が強まるという予測。

気孔コンダクタンスモデル間の差は、地図上では目立たないが、、

MIROC-GCM出力

@IAR4-A1Bシナリオ

*

飽差の関数

（

モデル）

相対湿度の関数

（

モデル）

20世紀末の気候条件において、

大陸全体のGPPが大体一致す るようにパラメーターを推定

State at the end of

Sato et al. (2015) JGR-Bioge. 120(5)

21世紀中GPP変化量の、モデル間の差 分(Leuningモデル – BWBモデル)

いずれの気孔コンダクタンスモデルの元 でも上昇したが、LeuningモデルはBWB モデルよりも上昇幅が大きかった

LeuningモデルはBWBモデルに 比べ、特にサヘル域において生 産性の増加量が目立った

年GPP の今世紀中の変化（大陸平均）

なぜ、このような差が生じたのだろう、、、

(KgC m-2yr-1) (KgC m-2yr-1)

Leuningモデル BWBモデル

大陸全体の平均値

世紀中の変化

飽差

モデル 相対湿度

モデル

年

年蒸散量

展葉期間

日

日

いずれの気孔コンダクタンスモデルの元でも、変動気候 条件のもとにおいて、GPPは増加し蒸散量は低下した

植物生産力が主に降水量によって制限されているアフリカ大陸において、飽差モデルの

「より賢い」水利用が、高い植物生産性を示したと解釈できる

成長期間の、あるスナップショットだけを見ると、温暖化条件でより気孔を開放するRHモデルの方が、光合成速 度も高かったりする。しかしこれは水資源を、より多く消費することで、土壌含水率を低下させて、展葉期間を 短くしてしまう。

Big leafモデル

光合成速度と気孔コンダクタンスとの連立

SEIB-DGVMの場合。

細胞間隙CO₂濃度が平衡に達するまで、反復計算。

光合成速度は気孔コンダクタンスの関数で、同時に気孔コンダクタンスは光合成速度の関数

（Jarvis型モデルを除いて）。そこで、繰り返し計算を行うことで、両者を収束させる必要がある。

CLM4.5の場合 蒸散量に応じて葉温を制御している ため、細胞間隙CO₂濃度と葉温の両 者が平衡に達するまで、反復計算。

葉温（T_i）を収束させるループ

細胞間隙CO2濃度（c_i） を収束させるループ

光合成に関する各種パラ メーターは葉温（C_i）の関数 で与えられている

出典：Bonan et al. (2014) Geosci. Model Dev. 7

これら正負の効果 の合計は通常負

土壌の乾燥による光合成速度のダウンレギュレーション

土壌が乾燥すると、植物の水ポテンシャルが低下し、

これにより気孔コンダクタンスの低下、続いてCO₂ 不足による光合成速度の減速が生じる。

しかし、Jarvis型を除いて、気孔コンダクタンスモ デル土壌含水率を制御因子として考慮しないため、

土壌含水率による光合成速度のdown regulationが 別途必要になる。

図：加藤知道監訳「生態系生態学第2版」森北出版

有効土壌含水率

（萎れ点で0、圃場容水量で1）

c water

SEIB-DGVMでは以下の関数を光飽和光合成速度に乗してる

𝑐_{𝑤𝑎𝑡𝑒𝑟} =

2 × 𝑠𝑡𝑎𝑡_{𝑤𝑎𝑡𝑒𝑟} − 𝑠𝑡𝑎𝑡_{𝑤𝑎𝑡𝑒𝑟}²

stat_water

図の出典：彦坂幸毅著「植物の光合成・物質 生産の測定とモデリング」共立出版

葉肉抵抗

気孔全開状態の気孔抵抗とほぼ同等とのこと。明示的に扱 われている広域植生モデルの例を知らない。光合成のパラ メーターに吸収させている？

境界層抵抗

コンダクタンスで

とか値が大いので、通常無視さ れる

気孔抵抗

CO₂

取り込みの最大の障害。光合成速度をパラメーターに持 つ気孔コンダクタンスモデルを用いる場合には、光合成モデ ルと連立させる必要がある。

CO

取り込みに関わる3つの抵抗

実は気孔抵抗だけ ではありません

気孔コンダクタンスモデルは、実は今も発展途上である

) / 1

( ) CO

( _{2 3}

2

min k k

g

g 大気中 濃度 飽差

光合成速度







 



Leuning型

経験的なパラメーター

最適理論型モデルは、「気孔コンダクタンスは、失われる水蒸気あたりの光合成量を 最大にするように制御されている」という仮定の下で、解を得ている。

最適理論型















 





 ^飽差 

濃度 大気中

濃度 気孔内部の

光合成速度

1.6 1 CO

CO

2 2

g

炭素量に換算した

単位蒸散量あたりのコスト

最適理論モデルのパラメーターは１つ。しかも、その値の上限・下限も理論的に得られ る。また、このモデルは高

においても的確に挙動することが確認されている

この他にも葉の水ポテン シャルをパラメーターとして 取り入れた気孔コンダクタン スモデルも存在するらしい

グローバルな水循環

図の出典：加藤知道監訳「生態系生態学」

土壌中に存在し植物がアクセス可能な水は、地球 上の水のわずかに0.01%。これが陸上生物を支え ている。

陸は地球表面の約30%を覆っているが、そこから の蒸発散は、グローバルな蒸発の15%を占めてい る。→陸面の蒸発効率は海面の約半分。

大気に常時存在する水量は、蒸発と蒸散によって 年間に循環する水の総量の2.6%。よって大気中の 水の回転率は、約10日。←数時間から数週間オー ダーで生じる蒸発散と大気循環が、その時空間変 動を支配する。

土壌水分の平均回転率は約2ヶ月（ただし地域差が 大きい）。←Seasonalな時間スケールで生じる降 水と蒸発散とが、その時空間変動を支配する