温暖化にともなうブナ林の適域の変化予測と影響評価
Prediction and impact assessment of the changes in suitable habitats for beech (Fagus crenata)forests under climate change scenarios
松井 哲哉

1＊

・田中 信行 ・八木橋
2

1＊

勉 ・小南 裕志 ・津山 幾太郎 ・高橋 潔
3

4

2

2

3

5

4

Tetsuya MATSUI , Nobuyuki TANAKA , Tsutomu YAGIHASHI , Yuji KOMINAMI ,2
5
Ikutaro TSUYAMA and Kiyoshi TAKAHASHI（独）森林総合研究所 北海道支所・（独）森林総合研究所 植物生態研究領域・（独）森林総合研究所 東北支所
45（独）森林総合研究所 関西支所・（独）国立環境研究所 地球環境研究センター
1
Hokkaido Research Center, Forestry and Forest Products Research Institute2
Department of Plant Ecology, Forestry and Forest Products Research Institute3
Tohoku Research Center, Forestry and Forest Products Research Institute4
Kansai Research Center, Forestry and Forest Products Research Institute5
Center for Global Environmental Studies, National Institute for Environmental Studies
1

2

3

摘

要

温暖化がブナ林の分布に与える影響を評価するために、分布予測モデルを用いて、現在の気候および RCM20 と MIROC シナリオによる生態的生育適域（適域）の将来予測を全国と各地域で行った。また、温度を 1℃ずつ、降水量を 10％ずつ変化させた場
合の適域面積の変化率を温暖化影響関数として算出した。現在ブナ林が分布する地域
における適域面積は、RCM20 と MIROC では現在に比べ 2031 ～ 2050 年には 47％と32％に、2081 ～ 2100 年には 21％と 4％にそれぞれ減少すると予測された。影響関数
によって、温度上昇だけではなく降水量の減少も適域の面積を低下させることがわかった。気温が 2℃上昇かつ降水量が 40％増加すると全国レベルの適域面積は現在の 9割程度にとどまるが、気温が 2℃上昇かつ降水量が 40％減少すると 2 割以下になる。
九州、四国、本州太平洋側のブナ林は、分布面積がもともと狭い上に、将来、夏期の高温にさらされることにより適域面積が狭くなると予測された ｡ 現在はほとんどが適域である白神山地世界遺産地域とその周辺では、将来、夏期の高温のために適域が大幅に縮小すると予測された。北海道では適域が現在の分布北限を越えて北へ広がるが、ブナの分布拡大速度が遅い上に、分断化された天然林の配置ではブナの移動は困難なので、温暖化のペースに追いつけないだろう。西日本ではカシ類やモミなどが、また本州の低標高域ではコナラ、クリ、ミズナラなどが、ブナの衰退後に置き換わる可能性がある。ブナの寿命は 200 ～ 400 年であり、温暖化によりブナ林がすぐに衰退する
可能性は低い。ブナの老齢木の枯死後に高木種の交代が順調に進行するかについて、
ブナ林の大きな変化や衰退が予測される地域を中心に監視していくことが重要である。

キーワード：影響関数、温帯林、脆弱性、適応策、分布予測モデル、
メッシュ植生データ

Key words：impact function, temperate forest, vulnerability, adaptive management,
predictive distribution models, mesh vegetation data

1．はじめに環境条件から野生植物種の分布を予測する統計モデルを作り、生育地を予測する研究が近年進められている。環境条件の適した場所に植物種がくまなく広がっている場合は、モデルから予測される分布確率は種にとっての生育地条件に対応し、高い確率は生態的に生育に適している場所（適地、適域）を示すと考えられる。現在の気候条件で植物種の分布する地域を環境条件から予測する統計モデルを開発する
ことにより、将来の気候条件のもとでの潜在的な分
布確率を空間的に予測することができる。また、気候変化による適域の面積変化だけでなく、将来も適域である地域、現在は適域であるが将来は外れてしまう地域、現在は適域ではないが将来新たに適域と
1）－ 3）
なりうる地域が特定できる。このような研究は、

受付；2008 年 12 月 22 日，受理：2009 年 5 月 20 日＊
〒 062-8516 北海道札幌市豊平区羊ヶ丘 7 番，e-mail：tematsui@ffpri.affrc.go.jp
2009 AIRIES

165

松井ほか：温暖化にともなうブナ林の適域の変化予測と影響評価

広域における長期間の変化予測が可能という利点がある。衰退や絶滅の危険性のある脆弱な植物種や群落とその地域が特定できるので、温暖化影響に対す
る適応的保全計画の策定資料となる。また温暖化影
響検出のためのモニタリング地を選定する際の参考になる。
欧米では、地球温暖化の植物分布への影響を予測・評価した研究が多くの種類について行われてい4）－ 6）る。アメリカ東部では 80 種の樹木種について
4）
適域の変化予測が行われ 、ヨーロッパでは約 1,400
5），6）種の植物の多様性の変化予測が行われている 。
日本においては温暖化が植生帯や個々の植物種（例えばブナ、ハイマツ、ササ類など）の分布に与える影響について、1990 年代以降研究が行われてい1）－ 3），7）－ 12）る。これらのうち日本の代表的天然林タイプの 1 つであるブナ林については、IS92a 温室効果13）
ガス排出シナリオに基づくCCSR/NIES シナリオによる 2090 年代の分布確率が、分類樹モデルを用
10）
いて全国スケールで算出された 。この研究によれば、ブナ林の適域（種の生存に適しており、分布頻度または優占度が高くなりうる環境条件を備えた地
域）は、九州、四国、本州太平洋側ではほぼ消滅す
ると予測された。さらに、世界自然遺産地域の白神
山地においては SRES-A2 排出シナリオに基づく
RCM20 シナリオによる 2081 ～ 2100 年の予測が追加され、同地域のブナ林適域の面積は 0.1％に減少12）すると予測された 。
その後も、SRES 排出シナリオに基づく新しい気候
変化モデルが発表されている。一方で、大気循環モ
デルに基づく気候変化シナリオを用いないで、気温や降水量などの気候因子を感度解析的に変化させた
14）
多数回シミュレーションを行う温暖化影響関数 の15）開発も、ブナ林をモデルケースとして進められた 。
ブナは北海道南部から鹿児島県まで分布し、その面積は日本の天然林総面積の約 17％にあたる
2
23,000 km である。ブナが林冠を優占するブナ林は
北海道南部、東北、本州日本海側に広く分布し、本州の太平洋側、四国、九州では山岳上部などに分布が限られる。ブナ林は、野生生物の生息地として生16）
態学的な重要性が認められているばかりでなく 、世界遺産の白神山地のように世界的な評価を受ける森林である。したがって、温暖化によるブナ林への影響を予測評価し、適応策について論じることは、生態系の保全管理計画を策定する上で重要である。
本稿では、ブナ林の分布を高い精度で予測できる
10）
ENVI モデル をベースとして、近年公表された
RCM20 および MIROC 気候変化シナリオに基づく
ブナ林の適域の予測および、温暖化影響関数の予測結果を紹介し、温暖化が全国と各地域におけるブナ林の分布に与える影響について、既存の研究成果を踏まえながら生態学的な考察を加えることを目的とする。
166

2．温暖化影響評価の手法
2.1 分布予測モデル
本研究で用いたブナ林の分布予測モデルである
ENVI モデルは、統計モデルの一種である分類樹モ
18）
デル である。分類樹モデルはデータをそれ以上分離しても無意味になるまで、可能な限り等質になる
ように 2 つに分離していくことによって作成される
19）モデルである 。このモデルは目的変数、説明変数
ともに数値変数と因子変数の両方を扱うことが可能な点や、モデルの構築が簡易で頑健な点、説明変数
間の交互作用関係を処理できる点、解釈が容易な点
20）で優れている 。よって、植生の分布と環境傾度の4），21）－ 23）
関係解析にしばしば用いられる。また、ブナ林の日本全国の分布を気候変数から予測する場合
では、一般化線形モデル（GLM）や一般化加法モデ
17）ル（GAM）よりも予測精度が高かった 。
ENVI モデルより、ブナ林の分布確率（自然植生
2
域の中でブナ林が出現する確率）が全国 1 km 解像10）－ 12），15），17），23）－ 26）度で算出可能となった。ブナ林の在・不在を分布確率から予測するために、モデルの27）
実際の分布に対する適合度を示す Kappa 係数 を、
分布確率 0 から 1 の間で 0.01 刻みで変化させながら求めたところ、分布確率が 0.47 の場合に Kappa
係数が最大 0.57 となり、当てはまりの良いモデル
28）と判定された 。そこで、分布確率 0.5 以上の地域
をブナ林の適域と呼称した。
本稿では、生態的な適域を面積として算出するた
めに、Kappa 係数を用いた適合度に基づいて分布確率 0.5 未満を適域外とした。そのために、例えば確率 0.3 でブナ林が分布する地域は適域外になるが、分布確率が閾値以下でも分布可能性がある場所は、
1）3）
分布辺縁域 もしくは Marginal habitats と呼ばれる。本稿ではブナ林の成立に適している地域が、温暖化によってどの程度変化するのかをわかりやすく示すことを目的とするので、適域についてのみ議論を進める。
分布予測モデルである ENVI モデルの目的変数
は、3 次メッシュ植生データベースから抽出した全自然植生データにおけるブナ林の在・不在データで29）ある 。これは、環境庁第 3 回自然環境保全基礎調査において 1979 ～ 1986 年に作成された 1：50,000
植生図を 3 次メッシュセル（緯度 30 秒、経度 45 秒の約 1 km メッシュセル）ごとにデジタル化したデータベースである。このデータベースは、日本全国をカバーする 345,000 以上のメッシュセルの植生タ
イプデータから構成されている。各メッシュセルの植生タイプは、メッシュセルの中央から半径 125 m
29）の円内に優占する植生タイプで代表させている 。
これから抽出したブナ林セル 23,432 点を含む自然
植生メッシュセル 156,804 点を目的変数として使用した。

地球環境 Vol.14 No.2

説明変数には、気候変数 4 つと環境変数 5 つを用30）いた。気候変数には暖かさの指数（WI） 、最寒月の
日最低気温の月平均値（TMC）、夏期降水量（5 ～ 9月）（PRS）、冬期降水量（12 ～ 3 月）
（PRW）
を用いた。
環境変数には、地質、地形、土壌、斜面方位、斜面
傾斜度を用いた。このうち地質、地形、土壌データ
は類型が都道府県間で異なる場合があるので、地質は 21 類型、土壌は 18 類型、大地形は 15 類型に類31）型を統一したものを利用した 。斜面方位と斜面傾
斜度は各 3 次メッシュ中心点から半径 125 m 圏内の 50 m デジタル標高データから計算した斜面方位10）と傾斜角度の値を使用した 。
気候変数では、モデル構築と現在のブナ林の分布
確率予測のために 3 次メッシュ気候値である気象庁
32）
観測平年値 を用いた。将来の 2 つの気候変化シナ
33）
34）
リオには RCM20 と MIROC の 2 期間（2031 ～2050 年と 2081 ～ 2100 年）を用いた。シナリオの年32）
平均気温は気象庁観測平年値 と比較すると、2031～ 2050 年の RCM20 と MIROC ではそれぞれ 2.2℃と 2.3℃の上昇、2081 ～ 2100 年の RCM20 とMIROC ではそれぞれ 2.8℃と 4.4℃の上昇である。
本研究では本州を太平洋側と日本海側に区分して統計値の算出を試みたが、この場合の両地域の区分
35）
は吉野（1981） によった。
2.2 ブナ林の温暖化影響関数の開発
全国と各県のブナ林適域面積率が多様な気候変化
シナリオで予測できるように、現時点で利用が可能
な複数の気候シナリオが予測する気温と降水量の変14）化幅（図 1）におけるブナ林の温暖化影響関数 を作15）成した 。温暖化影響関数とは、ブナ林の分布予測

図1

 平均気温変化および年平均降水量変化率の
年
県別・年別プロット．
□：SRES-B1 シナリオでの 2001 ～ 2010 年の 20 個のGCM の予測値
×：SRES-B1 シナリオでの 2081 ～ 2090 年の 20 個のGCM の予測値△：SRES-A2 シナリオでの 2001 ～ 2010 年の 18 個のGCM の予測値○：SRES-A2 シナリオでの 2081 ～ 2090 年の 18 個のGCM の予測値

165－174（2009）

モデル（詳細モデル）を用いて、平均気温を 1℃ずつ、降水量を 10％ずつ変化させ、感度解析的に多数回
シミュレーションを行い、その出力を全国、または県別に集計することで得られる数表の集合体（データベース）である。すなわち、詳細モデルの出力を近似的に得るための簡易影響評価モデルと考えるこ
15）とができる 。ブナ林の影響関数の場合は、分布確率 0.5 以上を有する 3 次メッシュセル面積合計の全面積に占める割合（ブナ林適域面積率）を県別に計算し、それを影響関数の出力変数とした。
3．ブナ林適域の予測と温暖化影響評価
3.1 ブナ林の成立を阻害する要因
ENVI モデルは説明変数の 87 通りの組み合わせ10）で成り立つ分類樹モデルである 。これらの組み合
わせのうちブナ林の分布確率が低い地域において、ブナ林の分布を主要に制限する要因とその閾値は以下のとおりであった（図 2）。分布確率の低い北海道の大半は、冬期の寒冷と乾燥（TMC ＜－12.45、
PRW ＜ 498.5）（図 2a）が、新潟県以南の日本海側沿岸域では夏期の高温（WI ＞ 95.15）（図 2b）が、本州
内陸部や東北の阿武隈や北上山地では冬期と夏期の乾燥（PRW ＜ 330.5、PRS ＜ 794.5）（図 2c）が、宮城県以南の本州太平洋側地域と四国および九州では夏期 の 高 温 と 冬 期 の 乾 燥（ W I ＞ 8 9 . 5 5 、 P R W ＜498.5）（ 図 2d）が、それぞれブナ林の成立を阻害す
る要因と考えられる。
3.2 全国のブナ林適域の将来変化
ENVI モデルに基づき分布予測を行った結果、ブナ林の適域の面積は、3 次メッシュセル 1 つの面積
22を 1 km とすると、現在の気候下では 26,220 km であった（表 1 下段）。このうち、実際にブナ林であ2る地域は 56％にあたる 14,579 km であった。
RCM20 と MIROC シナリオによる将来予測では、
2
適域全体の面積は、現在の 26,220 km から 2031 ～2050 年には RCM20 と MIROC ではそれぞれ 86％
2
2（22,637 km ）と 56％（14,744 km ）に、2081 ～ 2100
22
年にはそれぞれ 37％（9,618 km ）と 21％
（5,427 km ）に減少すると予測された（表 1、図 3）。これらの適域のうち、ブナ林が実際に分布している地域だけを2対象とすると、現在（14,579 km ）に比べて 2031 ～2050 年には RCM20 と MIROC ではそれぞれ 47％
2
2（6,821 km ）と 32％（4,637 km ）に、2081 ～ 2100 年
2
2にはそれぞれ 21％（3,117 km ）と 4％（544 km ）に減少すると予測された（表 1、図 3）。2081 ～ 2100 年の適域の分布は、現在の気候に比べ、本州太平洋側
や西日本でほとんどなくなり、本州日本海側でも大きく縮小すると予測された（図 3）。
温暖化影響関数に基づく分析によれば、気温が現状より 2℃上昇する場合、同時に降水量が 40％増加
すると、全国レベルでの適域の面積合計は現状の 9167

松井ほか：温暖化にともなうブナ林の適域の変化予測と影響評価（a）ターミナルノード 1
分布確率 0.0024
メッシュ数 56,676
PRW＜498.5
WI＜89.55
TMC＜−12.45（b）ターミナルノード 87
分布確率 0.0040
メッシュ数 27,195
PRW＞498.5
WI＞95.15（c）
ターミナルノード 2
分布確率 0.0082
メッシュ数 15,146
PRW＜330.5
WI＜89.55
TMC＞−12.45
PRS＜794.5（d）
ターミナルノード 32
分布確率 0.0009
メッシュ数 121,154
PRW＜498.5
WI＞89.55

図2

表1

 NVI モデルにおける，ブナ林の分布確率が低く面積の広い 4 つのE
ターミナルノードを規定する条件と対応する地域．

地域別・気候シナリオ別の適域（確率 0.5 以上）の変化．
現在の気候
RCM20
RCM20
MIROC
MIROC
（1952 ～ 1983 年） （2031 ～ 2050 年） （2081 ～ 2100 年） （2031 ～ 2050 年） （2081 ～ 2100 年）

北海道 ＊2＊1（3,089 /81,497 ）＊5本州日本海側（17,885/92,378）

1,695 /3,348＊4

1,484/11,716

413/4,062

1082/7,888

193/3,297

12,424/21,677

4,948/9,100

2,601/5,137

3,384/6,254

346/1,910

本州太平洋側（1,995/122,082）

456/1,181

384/1,809

103/419

171 /602

5/220

四 国
（223/16,415）

0/5

3/7

0/0

0/0

0/0

九 州
（240/32,741）

4/9

2/5

0/0

0/0

0/0

6,821/22,637

3,117/9,618

4,637/14,744

544/5,427

＊5

合
＊1
＊4

計

14,579/26,220

ブナ林の 3 次メッシュセル数．
ブナ林適域の 3 次メッシュセル数．

図3

168

＊3

＊2
＊5

対象地域の 3 次メッシュセル数．
＊3 ブナ林かつ適域である 3 次メッシュセル数．
35）
本州太平洋側と日本海側の区分は吉野 による．

ブナ林の分布（a）と，各気候条件において予測されたブナ林分布確率の分布（b～f）．
ブナ林の適域は，分布確率 0.5 以上（赤色）の地域である．

165－174（2009）（％）

地球環境 Vol.14 No.2

図4 気
 温と年間降水量を一定割合ずつ変化させた場
合の，ブナ林適域（確率 0.5 以上）の全国面積に対する割合（%）の等値線図（温暖化影響関数）．
現状は座標（0,0）である．

図5

気 温と年間降水量を一定割合ずつ変化させた
場合の，
ブナ林適域
（確率 0.5 以上）の県別面積割合（%）の等値線図（県別温暖化影響関数）の例．

（a）北海道，（b）富山県，（c）静岡県，（d）鹿児島県

割となるが、降水量が 40％減少すると現状の 2 割以下に減少する予測となった（図 4）。また気温が4℃上昇すると、降水量が 40％増加しても適域面積は現状の 3 割となる。このように、温度上昇だけでなく、降水量の減少がブナ林の分布確率を低下させることが今回の解析で判明した。
3.3 北海道への影響
ブナ林の天然分布北限は、北海道の渡島半島にある黒松内町、寿都町や蘭越町周辺にある。しかし、十勝地方の虫類村の更新世中期頃（約 78 万年前～
12 万年前）の地層からは、ナウマン象の化石ととも
36）にブナの殻斗や種子が出土している 。この時代は現在よりもやや温暖であったと考えられている。一方、現在よりも 5℃程度寒冷であったとされる約 2万年前の最終氷期の最寒冷期には、ブナ林の分布の中心は本州の新潟県から福島県以南にあったと考え
37），38）
られている。そして、過去 1 万年間の温暖化にともない、ブナ林は分布を北に拡大し、約 1,000 年
39）
前に現在の北限域に到達したとされる 。
このように、
ブナは過去の気候変動に対応して分布拡大縮小を繰り返してきた。現在のブナ林の北限は夏期降水量で26）
分布が制限されている可能性が指摘されている 。
2
現在の気候における適域の面積は 3,348 km である。2081 ～ 2100 年の MIROC シナリオでは
2
3,297 km に減少するものの、それ以外は 120％2
2
（4,062 km ）～ 350％（11,716 km ）
に増加する
（表 1）
。
適域は、ブナ林の分布北限を越えて北東に拡大するが、2081 ～ 2100 年には縮小に向かい断片化する（図 3e、f）。温暖化影響関数の結果も、気温上昇と同時に降水量が増加すれば適域の面積は増加する可能性を示した（図 5a）。温暖化にともなう適域の拡
大は、北海道のブナ林の分布を制約している冬期の
寒冷と乾燥（図 2a）が温暖化によって緩和されるた−12.25（温暖化シナリオ）−12.25（現在の気候 32））
−12.25（現在の気候 32））

（A）
RCM20
（2081 〜 2100 年）

図6

（B）
MIRCO
（2081 〜 2100 年）

北 海 道 に お け る ， 最 寒 月 最 低 気 温（ T M C ）が－12.25℃以上と未満の地域の変化予測．

めと考えられる。
現在の温暖化の水平移動速度は、花粉分析から推定された過去のブナ林の移動速度よりも速い。花粉分析の研究から推定される最終氷期以降のブナ林の
40）
，
41）
推定北進速度は北海道で 1.1 ～ 2.0 km/100 年である。RCM20 と MIROC シナリオの将来の北海道における最寒月最低気温（TMC）の北進距離は、例
え ば －1 2 . 2 5 ℃ の 等 温 線 を 例 と す る と 約 1 0 ～50 km/100 年である（図 6）。ブナの移動速度を本州38）の推定値である 23.3 km/100 年 と仮定しても、温暖化による等温線の移動速度に追いつけない可能性が高い。
ブナがスムースに移動するためには天然林が連続している必要があるが、現在の土地利用は人工林、農耕地、都市などが天然林を分断しているため、将来のブナの移動は阻害されると予想される。さらに
将来は、適域面積は増加しても断片化するため、ブナの分布移動は一層困難となるだろう。ブナ林の変化をより正確に予測するためには、ブナの移動速度
や天然林の分断化を考慮した予測モデルを開発する必要がある。
169

松井ほか：温暖化にともなうブナ林の適域の変化予測と影響評価
32）

図7

白神山地世界遺産地域とその周辺におけるブナ林の分布確率の変化予測．

（a）現在の気候下での分布確率，（b）RCM20 シナリオ（2031－2050 年），（c）RCM20
シナリオ（2081－2100 年），（d）MIROC シナリオ（2031－2050 年），（e）MIROC シナリオ
（2081－2100 年）．ブナ林の適域は，分布確率 0.5 以上（赤色）の地域である．

3.4 本州日本海側への影響
2
本州の日本海側では、ブナ林が 17,885 km と広く分布する（表 1）。この地域の適域の面積は現在の2
気候下では 21,677 km あるが、RCM20 と MIROC
2
2
両シナリオ下で 9％（1,910 km ）～ 42％（9,100 km ）に減少する（表 1）。また、実際にブナ林かつ適域で2ある地域は、現在の 12,424 km から将来は 3％22（346 km ）～ 40％（4,948 km ）に減少する。山陰、北陸地方では実際にブナ林かつ適域である地域はほとんど消失し、新潟から青森県にかけての地域でも縮小する（図 3）。これは、日本海側沿岸の低標高域でブナ林の成立阻害要因となっている夏期の高温（図 2b）
が、将来は高標高域にも到達するためと考えられる。
一方で、富山県の温暖化影響関数は、適域面積の減少率が比較的小さく、降水量の変動にもあまり左右されない結果となった（図 5b）。これは、温暖化後に適域が移動できる高山があることと、現在の降水量が多いことが原因と考えられる。
世界遺産で青森県と秋田県にまたがる白神山地では、適域が現在の気候下では世界遺産地域の 95％を占めるが、2031 ～ 2050 年には RCM20 で 30％、MIROC で 2％に減少し、2081 ～ 2100 年には RCM20
で 0.6 ％、MIROC では消滅する（図 7）。RCM20 シ
ナリオの 2081 ～ 2100 年には、白神岳（1,232 m）
周辺
12）の狭い地域が適域として残るだけとなる 。
自然遺産地域の施業管理計画図によると、約 8 割
が林齢 150 ～ 200 年生であるので、2100 年頃には12）
多くのブナが壮齢期から老齢期を迎える 。このため、将来は適域でなくなることを考慮すると、ブナ林下限域から、風倒等により発生する林冠ギャップの形成を契機として他樹種との競争によりブナが排除され、ブナの立木密度の低下が起こると考えられ
170

図8

白 神山地世界遺産地域とその周辺における現在の気候と気候変化シナリオ下での適域の垂直分
布変化予測．

標高値は図 7 に表示された地域に含まれる 2,470 個の 3
32）
次メッシュの標高値を用いた．
「現在」は気象庁 ，「R31」
は RCM20 シナリオ 2031－2050 年，「R81」は RCM20 シ
ナリオ 2081－2100 年，「M31」は MIROC シナリオ 2031
－2050 年，「M81」は MIROC シナリオ 2081－2100 年を示す．各気候シナリオごとに適域と判定された 3 次メッシュセルの数を「n ＝ 」で示す．箱ヒゲは，最小値，下四分位値，中央値，上四分位値，最大値，および外れ値（黒丸）を表示．

る。その場合、この地域の競合種としては、ミズナ
ラ、
コナラ、クリ、
アカシデ、
トチノキ、
ケヤキ、
ハリギ
12）
リ、シナノキなどの落葉広葉樹種が考えられる 。
白神岳におけるブナ優占林分の上限は標高1,070 m（WI≈44）であるが、ブナ単木の上限は標高
1,120 m（WI≈42）で、ミヤマナラ、ダケカンバ、ミネ
カエデ、ナナカマド等の低木類とチシマザサが混生12）する偽高山帯群落に侵入している 。温暖化によりブナ林の適域が高標高域へシフトすると（図 8）、積
雪期間の減少、生育期間の長期化、WI の上昇にともない、現在は標高 900 m 以上に分布している偽高山
12）帯群落へブナの侵入が活発になる可能性がある 。

地球環境 Vol.14 No.2

一方で、自然遺産地域の外側で、今回予測に用いた 4 通りの気候変化のいずれにおいても適域が存在するのは岩木山（1,625 m）
のみである
（図 7）
。過去に
使用された CCSR/NIES シナリオ
（2091 ～ 2100 年）
12）においても岩木山の一部は適域と判定された 。したがって岩木山は、この地域のブナ林の逃避地としての保全価値が比較的高いといえる。
3.5 本州太平洋側への影響
この地域のブナ林適域の面積は現在の気候下では
2
1,181 km と予測された（表 1）。RCM20 と MIROC 両
2
2
シナリオ下で 19％（220 km ）～ 153％（1,809 km ）に変化する（表 1）。また、実際にブナ林かつ適域であ22る地域は、現在の 456 km から将来は 1％（5 km ）～
2
84％（384 km ）に減少する。
RCM20 シナリオの 2031 ～ 2050 年では、適域の面2積が 2031 ～ 2051 年に現状の約 1.5 倍（1,809 km ）に増加するが、それらの地域のほとんどは新規に適域となる地域で、長野県と岐阜県にまたがる飛騨山脈
南部、富士山、北関東の一部分が該当する（図 3c）。
これは、RCM20 シナリオの 2031 ～ 2050 年において、
本州太平洋側の降水量が増加することに起因すると考えられる。静岡県の温暖化影響関数においても、
気温が 2℃上昇かつ年平均降水量が 20％増加すれば、適域の面積はわずかながら増加することと矛盾しない（図 5c）。
それ以外の気候シナリオでは、適域は現状の2
19％～ 51％（220 ～ 602 km ）に減少すると予測された（表 1、図 3）。これは、宮城県以南の本州太平洋側地域でブナ林の成立阻害要因となっている夏期の高温と冬期の乾燥（図 2d）が、ブナ林の分布する高標高域に到達するためと考えられる。本州内陸部や東北の阿武隈や北上山地においては夏期・冬期の乾燥がブナ林の分布を制限していると考えられるが（図 2c）、将来、降水量が増加しても、気温も同時に上昇するために適域は減少すると考えられる（表 1、図 3）。
将来、WI が 90 を超えて上昇し適域から外れると、ブナ優占林はほとんど成立しなくなる。このような気候におけるブナの存在形態としては、ブナ個体の分布密度に地形の影響が明確になることや、ブナが他樹種と混交して低密度で生存することがある。前者の例では、栃木県の低地に成立するブナ小林分（WI 103 ～ 105）で、コナラ林の中で北斜面にのみブ42）ナが生育していることが報告されている 。後者の例では、阿武隈山地における WI が 100 以上の温暖な地域のブナは、モミ、イヌブナ、コナラ、シデ類
43）などと混交することが報告されている 。このように適域から外れたブナは個体密度が徐々に低くな
り、WI が 105 ～ 110 以上となる地域では出現頻度
44）が著しく低くなると考えられている 。
低山の山頂にのみブナ林が孤立して存在する場合
は、高標高への逃げ場がないので、温暖化影響によ165－174（2009）

るブナ林の消失が起こりやすいだろう 。例えば茨城県の筑波山（877 m）では、ブナ林は頂上周辺にのみ分布しており、分布下限でアカガシと混交している。そのため、ブナ林内へのアカガシの侵入にともないブナの密度が低下する可能性がある。
3.6 四国、九州への影響
四国、九州地方のブナ林の分布下限標高は高く、
例えば九州北部で約 700 m より上部が、南部では約
46）
1,000 m より上部がブナの主要な分布域である 。
そのためブナ林の面積は小さく、現状では四国で
2
2
223 km 、九州では 240 km である（表 1、図 3a）。
これら現存するブナ林の多くは適域からはずれており、適域内に成立するブナ林は四国では存在せず、
2
九州で 4 km であった（表 1）。とはいえ四国や九州
でも、ブナ林が存在するエリアは存在しないエリアに比較して、分布確率は高くなっている。例えば、
2
四国と九州のブナ林約 463 km のうち、確率 0.3 以
2上の場所が 290 km ある。四国や九州のブナ林は人為影響によって減少した可能性があり、このようなブナ林の在データの不足によって ENVI モデルが地域的に予測精度の低下を引き起こした可能性がある
10）
ので 、結果の解釈には注意が必要である。
温暖化後は、RCM と MIROC いずれのシナリオにおいても四国と九州ではブナ林の適域はほぼ消滅すると予測されたので、温暖化に対して脆弱と考えられる。これは、本州太平洋側と同様に、四国や九州でブナ林の成立阻害要因となっている夏期の高温と冬期の乾燥（図 2d）が、ブナ林の分布する高標高域に到達するためと考えられる。ブナの南限が位置する鹿児島県の温暖化影響関数は、降水量変化の影
響は少なく、気温の上昇によって適域が減少するという結果を示した（図 5d）。
温暖化にともない、現在ブナ林が広がる標高域にまで WI 90 以上の高温域が上昇すると、アカガシやウラジロガシなどの常緑広葉樹のブナ林内への侵入が起きて、ブナの密度が減少する可能性がある。この地域のブナは、アカガシなどの常緑広葉樹やモミ47）類などの針葉樹としばしば混生する 。九州西南部にそびえる紫尾山（1,067 m）はブナの分布南西限に位置する山で、WI は約 89 であるが、頂部の北西斜面ではブナの密度が高く、反対側の斜面ではアカガ
48）
シ、ウラジロガシが成立する 。一方、分布南限の鹿児島県高隅山系（大蓑柄岳 1,236 m）においてもブ
ナは尾根の北西側で個体数密度が比較的高く、南東49）側では低いことが報告されている 。この半世紀で高隅山系の WI は 15 ～ 20 上昇したと見積もられており、ブナの年輪成長幅は減少傾向にあり、近年の50）
温暖化との関連が示唆されている 。
45）

4．温暖化対策としての植生管理と適応策
温暖化により、森林の構成種に消失、移動、加入
171松井ほか：温暖化にともなうブナ林の適域の変化予測と影響評価

といった変化が起こる。ブナ林保全のための適応策としては、ブナの移動経路を確保するための生態的回廊の設定が挙げられる。温暖化により適域から外れてブナが衰退する地域では、温暖な気候に適する
コナラやモミ・カシ類などがブナに置き換われば、森林としての一定の生態系サービスを享受することができるだろう。衰退が進行するブナ林では、変化をそのまま受け入れるのか、ブナ林を維持するための対策を行うかどうかの判断が各地で求められるよ
うになるだろう。どうしてもブナ林として維持しなければならない理由がなければ、植栽等の人為を加えず、森林の変化を見守るためモニタリングしていくほうが賢い選択だろう。モニタリングの場所としては、大きな変化が予測される感受性の高い地域や衰退が予測される脆弱な地域を中心に選定することが望ましい。長期間にわたるモニタリングにより実際の植生の変化を把握しながら、具体的な適応策を検討していくことが重要である。いずれにしても地域の実情に合わせた保全管理計画が必要である。

温暖化により適域から外れて衰退が進行するブナ林では、変化をそのまま受け入れるのか、ブナ林を維持するための適応策を行うかどうかの判断が必要である。どうしてもブナ林を維持しなければならな
い理由がなければ、植栽等の人為を加えず、森林の変化をモニタリングすることが有効である。
謝辞

本稿を作成するにあたり、査読者の方々から大変有益なご意見をいただいた。ここに御礼申し上げます。本研究は、環境省地球環境研究総合推進費（S-4）
「温暖化の危険な水準及び温室効果ガス安定化レベ
ル検討のための温暖化影響の総合的評価に関する研究」で行った。
引用文献
1） 津山幾太郎・松井哲哉・小川みふゆ・小南裕志・
田中信行（2008）本州東部におけるチシマザサの潜5．結論

在分布域の予測と気候変化の影響評価．GIS －理
論と応用，
16，
11-25.

気候と土地的変数を組み込んだブナ林の分布予測モデルである ENVI モデルに、近年公表された温暖化シナリオである RCM20 と MIROC のいずれかをあてはめると、将来のブナ林適域の総面積は減少し、2081 ～ 2100 年には RCM20 では現状の 37％、
MIROC では 21％になる。このうち、ブナ林が実際に分布している地域の適域の変化をみると、現在に比べて RCM20 で 21％、MIROC で 4％に減少する。
北海道ではブナ林の分布を制約している冬期の寒さと乾燥が緩和されることから、適域は北方へ移動しながら面積が 120％～ 350％増加するが、同時に断片化する。また、人為的土地利用による天然林の断片化の影響があるため、新たな適域へのブナの移動は困難だろう。本州の日本海側では夏期の高温によってブナ林の分布が制限され、適域が高標高地域
へ上昇しながら 9％～ 42％に減少する。実際にブナ林かつ適域である地域は、3％～ 40％に減少する。
世界遺産の白神山地周辺では適域面積の著しい減少が予測された。本州太平洋側では主に夏期の高温と冬期の乾燥がブナ林の分布を制限しているが、RCM20 シナリオでは一時的に適域が増加するものの、2081 ～ 2100 年には現状の 19％～ 35％に減少
し、実際にブナ林かつ適域である地域は、1％～ 23％に減少する。四国や九州でも主に夏期の高温と冬期の乾燥がブナ林の分布を制限している。ブナ林は高標高地域にのみ分布するため、現状でも適域面積は狭いが、将来はほぼ消滅すると予測された。しかし、ENVI モデルの地域的な予測精度不足に起因している可能性もあるため、結果の解釈には注意を要する。
1722） 津山幾太郎・松井哲哉・堀川真弘・小南裕志・
田中信行（2008）日本におけるチマキザサ節の潜在
分布域の予測と気候変化の影響評価．GIS －理論と応用，
16, 99-114.
3） Horikawa, M., I. Tsuyama, T. Matsui, Y. Kominami
and N. Tanaka（2009）
Assessing the potential impactsof climate change on the alpine habitat suitability ofJapanese stone pine（ Pinus pumila）. LandscapeEcology, 24, 115-128.
4） Iverson, L.R. and A.M. Prasad（1998）Predicting
abundance of 80 tree species following climatechange in the eastern United States. Ecological
Monographs, 68, 465-485.
5） Bakkenes, M., J.R.M. Alkemade, F. Ihle, R. Leemans
and J.B. Latour（2002）
Assessing effects of forecasted
climate change on the diversity and distribution ofEuropean higher plants for 2050. Global Change
Biology, 8, 390-407.
6） Thuiller, W., S. Lavorel, M.B. Araujo, M.T. Sykes and
I.C Prentice（2005）Climate change threats to plant
diversity in Europe. Proceedings of the National
Academy of Sciences of the United States of America,
102, 8245-8250.
7） Uchijima, Z., H. Seino and M. Nogami（1992）
Probable shifts of natural vegetation in Japan due to
CO 2-climate warming. In: M. Shiomi et al., eds.,
Ecological Processes in Agro-Ecosystems, Yokendo
Publishers, 189-201.
8） Tsunekawa, A., H. Ikeguchi and K. Omasa（1996）
Prediction of Japanese potential vegetation distribution

地球環境 Vol.14 No.2

in response to climatic change. In: K. Omasa et al.,eds., Climate Change and Plants in East Asia. SpringerVerlag, 57-65.
9） 小島 覚（2001）環日本海地域の植生と気候温暖化の影響．富山県 .

165－174（2009）

proach for predicting species distributions. Journal of
Vegetation Science, 11, 679-694.
23）八木橋 勉・松井哲哉・中谷友樹・垰田 宏・田中
信行（2003）ブナ林とミズナラ林の分布域の気候条
件による分類．日本生態学会誌，53, 85-94.

10）Matsui, T., T. Yagihashi, T. Nakaya, H. Taoda, S.

24）Matsui, T., T. Yagihashi, T. Nakaya, N. Tanaka and H.

Yoshinaga, H. Daimar u and N. Tanaka（ 2004）

Taoda（ 2004）Climatic controls on distribution of
Fagus crenata forests in Japan. Journal of Vegetation

Probability distributions, vulnerability and sensitivity
in Fagus crenata forests following predicted climatechanges in Japan. Journal of Vegetation Science, 15,605-614.
11）田中信行・松井哲哉・八木橋 勉・垰田 宏（2006）天Science, 15, 57-66.
25）Matsui, T., K. Shimada, T. Yagihashi, T. Nakaya, H.
Taoda and N. Tanaka（ 2005）Assessing potential
distribution of beech（Fagus crenata）
in central Japan.

然林の分布を規定する気候要因と温暖化の影響予Journal of Agricultural Meteorology, 60, 439-444.

測－とくにブナ林について．地球環境，
11, 11-20.

26）Yagihashi, T., T. Matsui, T. Nakaya, N. Tanaka and H.

神山地ブナ林の気候温暖化に伴う分布適域の変化

Taoda（2007）Climatic determinants of the northernrange limit of Fagus crenata forests in Japan. Plant

予測．日本森林学会誌，
89, 7-13.

Species Biology, 22, 217-225.

12）松井哲哉・田中信行・八木橋 勉（2007）世界遺産白13）Yokozawa, M., S. Goto, Y. Hayashi and H. Seino（2003）
Mesh climate data for evaluating climate change
impacts in Japan under gradually increasing
atmospheric CO2 concentration. Journal of Agricultural
Meteorology, 59（2）, 117-130.
14）肱岡靖明・高橋 潔・久保田 泉（2006）統合評価モデ
ルを用いた温室効果ガス安定化濃度目標下におけるイネ・小麦の潜在生産性変化の国別影響評価．
環境情報科学論文集，
20, 19-24.
15）高橋 潔・松井哲哉・肱岡靖明・田中信行・原沢英夫（2008）
温暖化政策支援モデルのための県別ブナ27）Cohen, J.
（1960）
A coefﬁcient of agreement of nominal
scales. Educational and Psychological Measurement,
20, 37-46.
28）Landis, J.R. and G.C. Koch
（1977）The measurement ofobserver agreement for categorical data. Biometrics,
33, 159-174.
29）アジア航測（1988）第 3 回自然環境保全基礎調査．
環境庁．
30）吉良龍夫（1948）温量指数による垂直的な気候帯のわかちかたについて－日本の高冷地の合理的利用のために．寒地農学，2, 47-77.

林影響関数の開発．地球環境研究論文集，
16, 111-119.

31）吉永秀一郎・鈴木裕一・松倉公憲・小林 守・新井

16）村井 宏・山谷孝一・片岡寛純・由井正敏（1991）ブ

正（1994）国土数値情報を用いた酸性雨に対する感

ナ林の自然環境と保全．ソフトサイエンス社．

受性分布図の作成．日本土壌肥料学雑誌，65, 565-

17）Matsui, T., T. Nakaya, T. Yagihashi, H. Taoda and N.

568.

Tanaka（2004）Comparing the accuracy of predictive
distribution models for Fagus crenata forests in32）気象庁（1996）気象庁観測平年値（CD-ROM）．
（財）気Japan. Japanese Journal of Forest Environment, 4633）気象庁（2004）気候統一シナリオ，第 2 版．気象庁（2）, 93-102.

象業務支援センター.
気候・海洋気象部・気候情報課．

18）Clark, L.A. and D. Pregibon（1992）Tree-basedmodels. In: J.M. Chambers and T.J. Hastie, eds.,

34）K-1 Model Developers（ 2004）K-1 Coupled Model（MIROC）Description. In: H. Hasumi and S. Emori,Statistical Models in S. Wadsworth & Brooks/Coleeds., K-1 Technical Repor t 1. Center for Climate

Advanced Books & Software, 377-419.
19）Breiman, L., J.H. Friedman, R.A. Olshen and C.J. Stone
Classification and Regression Trees. Chapman &（1984）
Hall/CRC.
20）De’
ath, G. and K. E. Fabricius（2000）Classification
and regression trees: A powerful yet simple techniquefor ecological data analysis. Ecology, 81, 3178–3192.
21）Michaelsen, J., D.S. Schimel, M.A. Friedl, F.W. Davisand R.C. Dubayah（1994）Regression tree analysis of

System Research, University of Tokyo.
35）吉野正敏（1981）
日本の気候地域区分．災害の研究，
12, 39-59.
36）矢野牧夫（1978）ナウマンゾウ包含層から産出した植物遺体．十勝団体研究会（編）, 十勝平野．地学団
体研究会，
389-393.
37）Tsukada, M.（1982）Late-Quaternary development of
the Fagus forest in the Japanese archipelago.
Japanese Journal of Ecology, 32, 113-118.

satellite and terrain data to guide vegetation samplingand surveys. Journal of Vegetation Science, 5, 673-686.

38）Tsukada, M.（1982）Late-Quaternar y shift of Fagus
distribution. Botanical Magazine Tokyo, 95, 203-217.

22）Vayssieres, M.P., R.E. Plant and B.H. Allen-Diaz（2000）

39）紀藤典夫・瀧本文生（1999）完新世におけるブナのClassiﬁcation trees: An alternative non-parametric ap-

個体群増加と移動速度－北海道南西部の例．第四173

松井ほか：温暖化にともなうブナ林の適域の変化予測と影響評価

紀研究，38, 297-311.

45）田中信行・八木橋 勉・杉田久志・藤田和幸・林40）五十嵐八枝子（1994）北上するブナ．北海道の林木
育種，37, 1-7.

遺伝別冊，
17, 109-118.

41）滝谷美香・萩原法子（1997）西南北海道横津岳における最終氷期以降の植生変遷．第四紀研究，36,

46）宮脇 昭（編著）
（1981）日本植生誌九州．至文堂 .
47）西尾孝佳・福嶋 司（1996）九州地方のブナ林群落における組成分化の機構．植生学会誌，
13, 73-86.

217-234.
42）小倉洋志・矢板美智夫（1989）栃木県南那須町におけ る ブ ナ の 分 布 ． 栃 木 県 立 博 物 館 研 究 紀 要 ，6 ,48）水永博己（2005）分布南限のブナ個体群．森林技術，
759, 20-23.
49）板谷明美・水永博己・山本進一（2004）ブナ南限の107-116.
43）原 正利（2005）東日本太平洋側におけるブナの生態
地理と保全．原 正利（編），レッドデータ植物群落である太平洋側低地ブナ林の分布と保全に関する研究．1-75.
44）原 正利（2006）東日本太平洋側におけるブナの分布とその下限を規定する要因について．植生学会誌，
23, 1-12.

松井哲哉

Tetsuya MATSUI

東京都出身。専門は植物生態学・植生学。東京農工大学在学中にブナ林に興味を持つ。環境コン
サルタント会社を経て、ニュージ
ーランド・オタゴ大学にて森林生態学を学ぶ。その後、森林総合研究所にて地球温暖化がブナ林に与える影響についての研究を開始。分布予測モデルによるブナ林の生態的分布適域予測のほか、北限域のブナ林の動態についても研究中。現在、森林総合研究所 北海道支所 森林育成
研究グループ主任研究員。博士
（農学）。

田中

信行

Nobuyuki TANAKA

専門は森林生態学、農学博士。
東京農工大学・東京大学大学院で学ぶ。学生時代に、樹木の分布と気候・土地要因の関係を研究テー
マの一つとしていた。森林総合研究所では、1993 年に温暖化影響評価に関する研究プロジェクトに参加して以来、温暖化影響研究に関わってきた。IPCC 第 3次（2001）と第 4 次（2007）の報告書 WG2 のレビューワー。現
在、温暖化影響総合予測プロジェクト（S-4、2005 ～ 2009 年度、環境省）の森林生態系部門の責任者。森林総合研究所主
任研究員。

八木橋

勉

Tsutomu YAGIHASHI

専門は森林生態学、造林学。農学博士。北海道大学大学院農学研
究科博士後期課程終了。農林水産省森林総合研究所、国際農林水産業研究センターを経て、現在は森林総合研究所東北支所に勤務。鳥散布種子の発芽など天然更新に関わる課題と、温暖化の影響を評価する際の土台ともなる植物の分布と気候要因の関係解析などが主な研究テーマ。

174

哲・垰田 宏
（2003）
森林生態系への影響と森林管理．

高隈山系におけるブナ分布の特徴．中部森林研究，
52, 67-68.
50）Mizunaga, H., S. Sako, Y. Nakao and Y. Shimono（2005）Factors af fecting the dynamics of the
population of Fagus crenata in the Takakuma
Mountains, the southern limit of its distribution area.
Journal of Forest Research, 10, 481-486.

小南

裕志

Yuji KOMINAMI

専門は森林水文学。森林総合研究所関西支所森林環境グループ主任研究官。アメダスと衛星データを用いた積雪の広域分布の推定、
粘性圧縮モデルを用いた積雪層の圧密過程の評価、京都府南部の山城試験地において森林の全炭素交換量の評価、等に関する研究を行っている。

津山

幾太郎

Ikutaro TSUYAMA

専門は森林植生学、
植物地理学、
景観生態学。筑波大学大学院生命環境科学研究科博士課程単位取得
満期退学後、森林総合研究所にて、地球温暖化の森林への影響評価に関する研究に従事。現在、日本の森林の林床優占種であるササ類を対象として、その分布を規定する気候要因の解明と、地球温暖化による影響評価の研究を行っている。博士（農学）。

高橋潔

Kiyoshi TAKAHASHI

1973 年山形県鶴岡市生まれ。
博士（工学）。専門は、環境システ
ム工学、地球環境モデリング。大学在学中より地球温暖化の影響対策に関する統合評価モデルの開発の研究プロジェクトに携わり、地理情報システムを活用した全球規
模の農業・水資源分野の温暖化影響評価手法の開発を担当してきた。最近では、影響予測の不確実性の定量化及びその伝達に焦点を当てて影響評価手法の高度化に取り組むとともに、予期される悪影響を軽減するための適応策の評価についても重要な研究事項と認識している。