樹木の蒸散と森林水循環：Penman-Monteith式と森林の水収支

森林の水循環は、樹木の生理学（蒸散）・気象学（大気水蒸気）・水文学（流出・涵養）・気候学（雲形成・降水循環）が交差する重要分野です。本稿では蒸散の物理学、Penman-Monteith式の数理、グローバル水循環における森林の役割、樹液流計測と渦相関法、樹種別の蒸散特性、気候変動とAIモデリング、林業実務への含意までを最新の数値と出典で10K規模に詳述し、最後にFAQ 10項目で実務者の疑問に答えます。

クイックサマリ：蒸散と水循環

項目	内容・典型値
蒸散（transpiration）	樹木が葉から水蒸気を放出する生理過程、年400〜1,200 mm
蒸発散（evapotranspiration、ET）	蒸散＋遮断蒸発＋林床蒸発の総和、PM式で算定
主要規定要因	VPD（飽差）、気孔コンダクタンス g_s、葉面温度、風速、日射 R_n
定量モデル	Penman-Monteith式（FAO-56）、Priestley-Taylor式、Shuttleworth-Wallace 二層モデル
森林の蒸散割合	降水量の40〜60%（日本の森林、FFPRI 2020）
遮断蒸発（canopy interception）	10〜25%（針葉樹）、5〜15%（広葉樹）
1個体の蒸散量	20〜400 L/日、年100〜200 m³（大径木）
林分蒸発散量	2〜6 mm/日、年500〜1,200 mm（日本）
水利用効率（WUE）	2〜8 g-C/kg-H₂O、樹種で2〜4倍差
地球規模 ET	年65,000 km³（陸域）、太陽エネルギーの約25%消費

蒸散の物理学：気孔・水ポテンシャル・SPAC

蒸散は樹木が能動的に維持する生理現象です。葉内細胞間隙は相対湿度ほぼ100%で、外気との水蒸気圧較差（VPD）が駆動力。土壌（−0.03 MPa）→根→幹→葉（−1.5〜−2.5 MPa）→大気（−50〜−200 MPa）の連続体（SPAC：Soil-Plant-Atmosphere Continuum）をポテンシャル勾配で水を吊り上げる「凝集張力説」が現代の標準理論です（Tyree & Zimmermann 2002）。

1. 気孔の制御：気孔密度は樹種で50〜500個/mm²、気孔コンダクタンス g_s は0.05〜0.4 mol/m²/s。光、CO₂濃度、ABA（アブシジン酸）、葉水ポテンシャルにより秒〜分スケールで動的制御される。

2. 水ポテンシャル勾配：水流 J = (Ψ₁−Ψ₂)/R（van den Honert 1948）に従い水が流れる。樹高80 mのセコイアでは梢頂で約−2.0 MPaの追加負荷で樹高上限要因（Koch et al. 2004 Nature）。

3. 葉温調節：蒸散1 mm（=1 L/m²）で気化熱2.45 MJ/m²を奪い葉温を3〜10°C低下。都市林業のヒートアイランド削減効果（前A10記事）の物理的根拠。

4. 光合成とのトレードオフ：WUE = A/E（光合成速度／蒸散速度）は2〜8 g-C/kg-H₂OがC3植物のレンジ。CO₂濃度上昇でWUEは10〜25%向上（Ainsworth & Long 2005）。

5. キャビテーション：木部に過剰張力が掛かると道管にエアブロックが生じ通水阻害（前D12記事）。被子植物 P50 = −1.5〜−3 MPa、針葉樹 −3〜−6 MPa（Choat et al. 2012 Nature）。

Penman-Monteith式：物理ベースの蒸発散モデル

世界標準はPenman-Monteith式（FAO-56準拠）。1948年Penmanの蒸発式に1965年Monteithが気孔・空気力学抵抗を組み込み、1998年FAO-56として国際標準化。

λET = (Δ(R_n − G) + ρ_a c_p (e_s − e_a)/r_a) / (Δ + γ(1 + r_s/r_a))
（λ：気化潜熱 2.45 MJ/kg、Δ：飽和水蒸気圧曲線の傾き、γ：乾湿計定数 ≈0.066 kPa/K）

パラメータ	意味	典型値（森林）
R_n（純放射）	下向き−上向き放射収支	50〜200 W/m²（夏日中）
G（地中熱フラックス）	地表面下方向熱流	R_n の5〜10%
(e_s − e_a)（VPD）	飽和－実際水蒸気圧差	0.5〜3 kPa
r_a（空気力学抵抗）	境界層・粗度長	10〜100 s/m
r_s（表面抵抗）	気孔抵抗の樹冠平均	50〜500 s/m

FAO-56は基準作物（草地、r_s=70 s/m）のET₀標準で、森林適用には作物係数 K_c（針葉樹1.0〜1.2、広葉樹1.0〜1.3）またはr_s直接設定が必要。Shuttleworth-Wallace二層モデル（1985）は疎林・若齢林の精度を向上させます。

森林水収支の構成要素：5フローの定量配分

日本の典型的森林（年降水量1,500〜2,500 mm）の水収支は以下のように配分されます。林分構造・樹種・気候帯で変動するため、施業設計の出発点となります。

水フロー	割合（年降水比）	典型量（年）	主要規定要因
降水（input、P）	100%	1,500〜2,500 mm	気候、地形、海洋距離
樹冠遮断蒸発（I）	10〜25%（針葉樹） 5〜15%（広葉樹）	150〜500 mm	樹冠貯留容量 S（1〜3 mm）、降雨強度
樹幹流（stemflow）	1〜10%	15〜200 mm	樹皮粗度、樹冠形態
林内雨（throughfall）	70〜90%	1,200〜2,000 mm	遮断ロスの残差
蒸散（T、樹木→大気）	40〜60%	500〜1,200 mm	VPD、g_s、LAI、根系深度
林床蒸発（E_soil）	5〜15%	75〜300 mm	林床日射、土壌湿度
表流水流出（Q_s）	5〜10%	75〜200 mm	傾斜、土壌透水性
地下水涵養（Q_g）	15〜30%	225〜600 mm	下層浸透能、岩盤亀裂
樹体内貯留（ΔS_tree）	≦1%	15〜25 mm相当	木部含水量、季節変動

森林の蒸発散合計（ET = T + I + E_soil）は降水の50〜75%を占め、流出量を大きく規定します。北海道・東北の冷温帯林ではET比率がやや低く（45〜60%）、九州・四国の暖温帯林では高い（60〜75%）傾向があります（FFPRI 各地試験地長期データ）。

グローバル水循環における森林の役割

地球規模では陸域からの蒸発散は年65,000〜70,000 km³、太陽放射の約25%を消費し（IPCC AR6 WG1, 2021）、森林がその40〜50%を担います。

1. 降水リサイクル比：陸域降水のうち再蒸発由来の割合。アマゾン流域25〜50%、コンゴ30〜45%、東アジアモンスーン域20〜35%、ユーラシア内陸40〜60%（van der Ent et al. 2010）。

2. 大気の川（atmospheric river）：アマゾンの「飛ぶ川」は南米南部降水の30%を供給（Marengo et al. 2018）。

3. 森林破壊の遠隔影響：熱帯林の30%以上を伐採すると下流圏で降水が10〜20%減少するシミュレーション（Lawrence & Vandecar 2015 Nat. Clim. Change）。日本も東アジアモンスーン水循環に組み込まれ、流域単位の水資源確保で森林保全が重要。

4. 同位体水文学：水の δ¹⁸O・δD 比から起源・経路を追跡し「いつの雨水を樹木が吸っているか」を定量可能（Sprenger & Allen 2020）。

樹液流計測：個体スケールの蒸散測定

個体・林分スケールで蒸散を実測する標準法は樹液流（sap flow）計測。Granier法（HD）を中心に多様な手法が確立しています。

計測法	原理	精度	対象
Granier法（HD）	加熱・非加熱プローブ温度差	±10〜20%	DBH 5 cm以上
Heat Pulse法（HPV）	パルス加熱の伝播速度	±5〜15%	低流速・逆流可
HFD（Heat Field Def.）	放射状熱場の歪み	±10%	大径木、深さ別
Stem Heat Balance	幹熱量バランス	±5〜10%	小径木・若齢樹

Granier法は1985年発表で世界最普及。流速 V = α·((ΔT_max − ΔT)/ΔT)^β（α=119, β=1.231）で算出、辺材深度2〜3 cm、長期連続観測可能。FFPRIは富士吉田・札幌・宝川等で30年以上のデータを蓄積し、スギ・ヒノキ・ブナのアロメトリ式を公開しています。

渦相関法：林分スケールのフラックス観測

森林ETを直接実測する標準手法が渦相関法（Eddy Covariance、EC）。FLUXNET・AmeriFlux・AsiaFluxの500地点以上で連続観測が稼働中。

1. 観測原理：観測タワー（樹高1.5〜2倍）に3次元超音波風速計（10〜20 Hz）と赤外水蒸気分析計。鉛直風速 w’ と水蒸気濃度 q’ の共分散 LE = ρλ·⟨w’q’⟩ から潜熱フラックスを直接算出。

2. 計測範囲：風上200〜2,000 mの林分平均。

3. AsiaFlux：富士吉田スギ林・高山落葉広葉樹林・苫小牧カラマツ林等で20年以上の連続データ。

4. データ処理：エネルギー保存（R_n = LE + H + G）検証、ギャップ充填、Webb-Pearman-Leuning補正等。FLUXNET Tier 1データは品質管理済み。

5. 衛星検証：MODIS MOD16・GLEAM・PML-V2 等のグローバルET推定アルゴリズムは FLUXNET を検証基準とする。誤差±15〜30%。

樹種別の蒸散特性：保守的 vs 寛容的戦略

樹種により蒸散戦略は大きく異なります。乾燥応答による分類はisohydric（等水ポテンシャル型、保守的）とanisohydric（変動水ポテンシャル型、寛容的）の二極で整理されます（Tardieu & Simonneau 1998）。

樹種	個体蒸散量 （L/日）	林分ET （mm/日）	気孔戦略	P50（MPa）
スギ（Cryptomeria japonica）	30〜120	2.5〜4.0	isohydric	−3.5〜−4.5
ヒノキ（Chamaecyparis obtusa）	25〜90	2.0〜3.5	isohydric	−4.0〜−5.0
カラマツ（Larix kaempferi）	40〜180	2.5〜4.5	中間	−3.0〜−4.0
アカマツ（Pinus densiflora）	35〜130	2.0〜3.8	isohydric	−4.5〜−5.5
ブナ（Fagus crenata）	50〜250	3.0〜5.0	中間	−2.5〜−3.5
ミズナラ（Quercus crispula）	60〜220	3.5〜5.5	anisohydric	−2.0〜−3.0
コナラ（Quercus serrata）	60〜200	3.5〜5.5	anisohydric	−2.0〜−3.0
ケヤキ（Zelkova serrata）	80〜300	3.5〜5.5	anisohydric	−1.8〜−2.8
ユーカリ（Eucalyptus grandis）	100〜400	4.0〜7.0	anisohydric	−1.5〜−2.5

1. Isohydric戦略：スギ・ヒノキ・マツ類は土壌乾燥時に早期気孔閉鎖で Ψ を一定保持。短期乾燥に強いが長期乾燥で炭素枯渇リスク。

2. Anisohydric戦略：ミズナラ・コナラ・ユーカリ等は気孔を開け続け生産性を維持するが、Ψ_leaf 低下でキャビテーションリスク増大。

3. 樹齢・サイズ効果：個体蒸散はDBHの0.7〜1.5乗にスケール（Meinzer et al. 2005）。DBH 30 cmスギで100〜150 L/日、DBH 60 cmブナで200〜300 L/日。

4. ユーカリの水資源圧迫：豪州・南米・地中海で個体蒸散300〜500 L/日に達し地域水資源の圧迫が問題化。本格導入には水文影響評価が不可欠。

5. 混交林効果：単一樹種より混交林の方が季節・年変動を平均化し蒸発散が安定（Forrester 2014）。

樹齢・林齢構造と水収支の動態

同一林分でも林齢で水収支は大きく変動します。日本のスギ・ヒノキ人工林を例にした典型的動態：

林齢	LAI	個体蒸散	林分ET	遮断率	流出比
幼齢林（〜20年）	1〜3	10〜30 L/日	1.5〜3.0 mm/日	5〜10%	30〜50%
壮齢林（20〜50年）	4〜7	100〜250 L/日	3.0〜5.0 mm/日	15〜25%	15〜25%
成熟林（50〜80年）	5〜8	150〜300 L/日	3.5〜5.5 mm/日	20〜30%	10〜20%
老齢林（80年〜）	4〜7	150〜250 L/日	3.0〜4.5 mm/日	15〜25%	15〜25%

老齢林で蒸散がやや減少するのは「水力学的制約仮説（hydraulic limitation hypothesis）」（Ryan & Yoder 1997）で説明される。樹高が増すと水の輸送距離・重力負荷が増え、最大気孔コンダクタンスが低下するためです。

気候変動と森林蒸散：相反する効果のせめぎ合い

気候変動は森林蒸散に複数の相反する影響を与え、地域・樹種で結果が異なります。

1. VPD上昇：気温＋1°CでVPDが約7%上昇（Clausius-Clapeyron関係）。グローバルVPDは1979〜2013年に0.05〜0.15 kPa/decadeで上昇し、森林衰退・乾燥誘発死亡の主要因（Yuan et al. 2019）。

2. CO₂施肥効果：CO₂濃度（産業革命前280 ppm→2026年430 ppm）で気孔開度5〜25%低下、蒸散抑制（Ainsworth & Long 2005, FACE実験メタ解析）。

3. 降雨パターン変化：日本は極端豪雨頻度が10年で5〜10%増加、無降水日数も増加傾向（気象庁2024）。樹木水ストレス頻度上昇。

4. 干ばつ誘発樹木枯死：Allen et al. 2010 以降、米国西部・地中海・豪州・アマゾン等で「forest die-off」報告増加。

5. ハイドロリック故障：世界樹種の70%が現在の最大渇水時に安全マージン1 MPa以下（Choat et al. 2012 Nature）。

6. 高緯度の蒸散増加：北欧・カナダ・シベリア等で融雪期前倒し・成長期延長で年蒸散量増加。亜熱帯・地中海性気候では減少（Jung et al. 2010 Nature）。

AIモデリングと衛星観測の最前線

近年は機械学習・衛星観測・地上観測の統合により、森林ETの予測精度が大きく向上しています。

1. 機械学習ベースのET推定：Random Forest、XGBoost、LSTMがFLUXNET2015＋衛星データ（MODIS、Sentinel-2、Landsat 9）を結合してET推定。FLUXCOM v1（Jung et al. 2019 Nature）でGRACE水収支との整合が向上。

2. 物理＋ML ハイブリッド：PM式出力をMLで補正する「PM-ML」（Pan et al. 2020）。

3. 衛星ETプロダクト：MODIS MOD16（500 m）、GLEAM v3（25 km）、PML-V2、ECOSTRESS（70 m、ISS）等。日本ではGCOM-C「しきさい」。

4. デジタルツイン森林：3D LiDARスキャン＋気象＋生理パラメータをリアルタイム結合。林野庁・FFPRIで国有林試験運用、森林経営計画への組み込みが進行。

森林施業と水資源：5つの介入軸

森林管理は流域水資源を5つの軸で操作可能です。各軸のトレードオフを定量的に把握することが、水源涵養と木材生産の両立に不可欠です。

施業	蒸散変化	流出変化	洪水ピーク	水質
皆伐＋再造林	−40〜−80%（伐採直後）	+30〜+200 mm/年（10年間）	+10〜+30%	濁度悪化、栄養塩流出
強度間伐（30〜50%）	−15〜−30%	+50〜+150 mm/年（5年）	+5〜+15%	軽微
弱度間伐（10〜20%）	−5〜−10%	+10〜+50 mm/年	+0〜+5%	影響小
長伐期化（80年〜）	安定	安定	低減	清浄
樹種転換（針→広）	−10〜+10%（樹種依存）	±100 mm/年	変動	季節変動増
水源涵養林化	ET最適化	ベースフロー増	洪水ピーク低減	清浄、安定

日本では「水源林造成事業」（森林整備センター、1965年〜）が、急傾斜地・水源地の森林を国費で整備する制度。年予算約100億円、累計対象面積40万ha。皆伐・強度間伐は短期的に流出量を増やし、洪水・土砂災害リスクを高めるため、急傾斜地での適用には慎重な設計が必要です（前D03記事「砂防ダム」、D08記事「森林の役割と気候緩和」参照）。

都市林業と蒸散冷却効果

都市環境では蒸散による潜熱輸送が「蒸散冷却」を発揮し、前A10記事と連携する重要メカニズムです。

1. 物理：蒸散1 mmで気化潜熱2.45 MJ/m²、地上1 m層で約2°C/h の冷却ポテンシャル。

2. 都市樹冠：成熟ケヤキ・クスノキ単独樹で日中-2〜-5°Cの局所気温低下。

3. ストリートツリー蒸散量：都市では森林の60〜80%、年50〜100 m³/個体（DBH 30 cmケヤキ）。

4. 緑地計画：i-Tree等で蒸散冷却を貨幣換算。国交省「みどりの基本計画」で導入実績増加中。

水源涵養林の経済評価

日本学術会議2001年答申は森林の多面的機能の経済価値を年74兆円と評価し、うち水源涵養関連は約30兆円を占めます。

機能	評価額（年・全国）	評価根拠
洪水緩和	6.5兆円	ダム代替コスト
水資源貯留（地下水涵養）	8.7兆円	ダム貯水量代替
水質浄化	14.6兆円	下水処理代替
土砂流出防止	28.2兆円	土砂ダム費用

流量平準化（基底流量増・洪水ピーク減）・水質浄化（窒素・リン除去率30〜70%）・土砂安定（表面流出抑制90%以上）が主要評価軸。近年はウォータークレジットが国内外で試行されています。

関連分野・前後記事との位置付け

• 蒸散の生理学（D10）：個体スケール蒸散の詳細は前D10記事「樹木の蒸散」を参照。
• キャビテーション（D12）：水ストレス時の通水阻害は別途D12で詳述。
• SPACモデル（D13）：本記事の理論的基盤、土壌-植物-大気連続体の数式モデル。
• 都市林業（A10）：蒸散冷却の都市環境応用。
• 森林炭素（A02・A08）：水収支とカーボンストックの統合評価。
• 砂防ダム（D03）：森林水収支と土砂災害防止。
• 気候緩和と森林（D08）：地球規模水循環と気候の相互作用。
• カラマツ（D14）・ブナ（D15）等樹種記事：樹種別蒸散特性の詳細。

よくある質問（FAQ）

Q1. 1本の樹木は1日にどれだけ蒸散しますか？

A. 樹種・サイズで異なります。DBH 30 cmスギで30〜120 L/日、DBH 60 cmブナで150〜250 L/日、ユーカリ大径木で300〜500 L/日。年間では大型樹木で100〜200 m³に達します。

Q2. なぜPenman-Monteith式が世界標準なのですか？

A. (1) 観測気象パラメータから直接計算可能、(2) 気孔抵抗 r_s で植生効果を組込可、(3) FAO-56として1998年国際標準化、(4) FLUXNET2015の500地点以上で誤差±10〜15%の検証実績、という理由から標準採用されています。

Q3. 森林伐採は水資源を増やしますか減らしますか？

A. 短期（5〜10年）は蒸散減で流出増（年降水量の5〜15%増）。しかし長期では洪水ピーク増・渇水期流量減・土砂流出・水質悪化のトレードオフで「使える水資源」は減少。急傾斜地・水源地の伐採には慎重な設計が必要です。

Q4. 都市の街路樹も蒸散冷却の効果がありますか？

A. はい。成熟街路樹（DBH 30 cmケヤキ等）単独樹で日中-2〜-5°Cの局所気温低下。樹冠下のヒートインデックス改善はさらに顕著で、ヒートアイランド緩和の主要メカニズム。前A10記事参照。

Q5. 気候変動で森林蒸散はどう変わりますか？

A. (1) 気温上昇＋VPD増（10年で+0.05〜0.15 kPa）で促進、(2) CO₂上昇で気孔閉鎖・抑制（−5〜−25%）、(3) 干ばつで水ストレス・場合により樹木枯死、(4) 高緯度では成長期延長で増加。日本ではモデル予測で正味+5〜+15%（IPCC AR6 RCP4.5）。

Q6. 樹液流計測で個体蒸散量はどこまで正確に測れますか？

A. Granier法で±10〜20%、HFD法で±10%、Stem Heat Balance法で±5〜10%。個体スケール実測の最高精度群です。樹液流から個体蒸散への換算には辺材面積・樹冠投影面積のアロメトリ式が必要。

Q7. 日本の森林の年間蒸散量は具体的にどれくらいですか？

A. FFPRI・AsiaFlux長期観測の集計では年蒸発散量500〜1,200 mm/年。森林面積2,500万haを掛けると年約1,750億〜3,000億トンで、国内年間水使用量（約810億トン）の2〜4倍に相当。

Q8. 雨水のうち地下水になる割合はどう決まりますか？

A. 森林タイプ・地質・地形・降雨強度の4要因。日本の森林で年降水の15〜30%。風化花崗岩・第三紀層で高く、堅固火山岩で低い。広葉樹は針葉樹より涵養率が3〜5%高い場合あり。林床土壌構造の保全が重要。

Q9. グローバル水循環における森林の独自性は？

A. 森林は陸域蒸発散の40〜50%を担い、(1) 深根系（5〜30 m）で乾季の蒸発を持続、(2) 高LAI（4〜8）で年蒸発散量が農地・草地より20〜50%多い、(3) 樹冠遮断で土壌侵食を防止、(4) 「飛ぶ川」で内陸降水を支える、という点で他生態系と一線を画します。

Q10. 自分の流域の水収支データを調べるには？

A. (1) 気象庁アメダス、(2) 国交省水文水質DB、(3) FFPRI試験地、(4) AsiaFlux/FLUXNET、(5) JAXA GCOM-C「しきさい」。簡易にはPM式（FAO-56）にアメダスデータを代入してExcel計算可能です。

まとめ：森林水循環の理解と実務への展開

森林水循環は、樹木の蒸散という生理現象を起点に、流域の水資源・気候・生態系を結ぶ中核プロセスです。本稿で示した通り、Penman-Monteith式という物理ベースのモデル、Granier法・渦相関法という実測技術、衛星観測・AIモデリングという最新ツールが整い、定量的な施業判断が可能になっています。気候変動下では蒸散・降水・流出のバランスが急速に変化しており、林業実務者・水資源管理者にとって最新データに基づく適応計画が不可欠です。前後の関連記事（D10〜D15、A02・A08・A10、D03・D08）と併せて、流域単位での統合的森林管理の参考としてください。

関連記事

• 植林技術カテゴリトップ
• 樹木の蒸散メカニズム（D10）
• キャビテーションと水ストレス（D12）
• SPACモデル（D13）
• 都市林業ヒートアイランド削減（A10）
• 世界森林炭素吸収量再計算（A08）
• 砂防ダムと森林水文学（D03）