気候変動が森林動態に与える影響と将来予測：
八甲田山のオオシラビソを例として
Impacts of climate change on forest dynamics of Abies mariesii（Pinaceae）in Hakkoda Mountains
田中 孝尚 ・嶋崎 仁哉・黒川 紘子・彦坂 幸毅・中静＊＊透Takahisa TANAKA , Masaya SHIMAZAKI, Hiroko KUROKAWA,
Kouki HIKOSAKA and Tohru NAKASHIZUKA東北大学大学院生命科学研究科
Graduate School of Life Sciences, Tohoku University

摘要

青森県八甲田山のオオシラビソ林を対象に，約 30 年前および最近の空中写真を利用してここ数十年間の温暖化が個体群に与えた影響を検出し，今後の分布変化を予測した。オオシラビソの個体群は，ここ約 30 年間で，分布の下限標高付近で減少し，上限付近では逆に増加しており，分布が高標高域へ移動していることが明らかになった。また，湿原の周辺でも個体数は増加していた。さらに，オオシラビソ個体群の分布と環境条件の関係を 50 m メッシュと 1 km メッシュで解析し，その結果を用いて温度上昇に伴う分布適域の変化を予測した。1 km メッシュでは，温度条件のみで分布が説明されたが，50 m メッシュでは，地形条件や湿原との距離などによって，適域の分布が影響を受けていることが分かった。これは，解像度を上げることによって見落としていた微地形などの適応的機能から局所的な適域を検出できたことによると考えられた。温度上昇にともない湿原の分布適域も減少するが，湿原が衰退しない場
合にはオオシラビソの分布適域の減少も少ないことが分かった。個体群動態の結果と合わせると，
湿原の周囲がオオシラビソのレフュジアとなる可能性が示された。また，
空間解像度の違いは，分布適域の推定に影響し，気候変動が森林動態におよぼす影響についての具体的適応策を考えるうえで，高い解像度による解析が必要なことが明確となった。

キーワード：オオシラビソ，山岳地域，地球温暖化，標高傾度，レフュジア
Key words：Abies mariesii, mountain region, global warming, altitudinal gradient,
refugia

1．はじめに世界の年平均気温は 19 世紀後半以降で 100 年あたり約 0.7℃の割合で上昇しており，温暖化の主要
因は，人間の産業活動にともなって排出された温室1）
効果ガスによるものと考えられている 。このような地球温暖化の影響は，世界的な規模での平均気温の上昇，極端な高温の増加や低温の減少，大雨の増
加，降雨の変化，乾燥化，海氷の融解や氷河および氷床減少による海面上昇など，さまざまな面で表れている。日本においても，気温の変化のみならず，
降雨の年変動，豪雨の頻度，積雪の増減などが起こると考えられ，その影響は森林生態系に生育する植物の環境や，これにともなう土壌微生物の活動などを含めた森林土壌の環境変化として現れる可能性がある。
将来的には，世界平均気温は 2090 ～ 2099 年にかけて 1.1℃～ 6.4℃上昇すると予測されており，この気候変動によって，植生の垂直分布帯の一部が減少，あるいは完全に消滅してしまう可能性がある。
特に，高緯度，高標高地域の森林は脆弱であること
1）
から ，気候変動に対する適応策が望まれている。
こうした適応策のなかでは，気候変動による植物の逃避地（レフュジア）となりうる環境を探ることが重要となる。
気候変動は，樹木の個体群動態に直接的，または
撹乱レジームの変化を介して間接的に影響をおよぼ2）
すのみならず ，その影響を通じた生態系の変化によって，生態系サービスにも影響をおよぼしてい
1），3）。近年の気候変化によって，植物種の分布がる

受付；2013 年 12 月 25 日，受理：2014 年 2 月 17 日＊
〒 980-8578 仙台市青葉区荒巻字青葉 6-3，e-mail：tanaka@m.tohoku.ac.jp
2014 AIRIES47

田中ほか：気候変動が森林動態に与える影響と将来予測

極方向あるいは高標高域へシフトしていることは明4）確であるが ，樹木に関する詳細な測定は稀であ5）る 。ほとんどの研究は，年輪や樹齢構成を調べる
ことによって，高木が生育できなくなる限界である
5）－ 8）。高木限界付近に
高木限界に焦点を当てている
おける樹木個体群の動態は気候変動の良い指標であ6）
るが ，年輪成長から個体群動態を調査することは難しい。なぜなら，年輪は生きている木からのみ採取され，死亡した個体の情報は得られないからである。大規模かつ長期的な野外調査を行うことによって，植物の種の有意な上方への移行を示した例もあ9）るが ，種の分布域全体における個体群動態に関する研究はほとんどない。過去の航空写真や衛星画像
によって植生の変化を把握することは可能である10）－ 12），リモートセンシングデータはデータ解像が度に限界があり，多くは個体群の動態解析には用いられてこなかった。
一方，野生植物の分布や植生に対する気候変動の影 響 予 測 の ほ と ん ど は， 空 間 解 像 度 が 粗 く（1 ～13），14）50 km） ，温暖化の影響を左右しうる局所，もしくは景観スケールで空間的に不均質な要因はあまり考慮されていない。こうした解像度の粗さにより，山岳植物の地球温暖化に適応する能力が過大評
15）
価されている可能性も示唆されている 。日本では， ブ ナ（Fagus crenata）の 分 布 変 化 を 約 1 km ×
1 km メッシュセルを使用して予測している例があ16）
るが ，この解像度では依然として地形のような局所的要因の影響を組み込むことは難しく，地域（地方自治体レベルを想定）または局所スケール（単独の事業レベルを想定）での実践的な適応策の計画と管17）
理には不十分である 。
八甲田山は本州最北部に位置する山岳地帯であり，高山帯，ハイマツ帯，亜高山帯や低層～高層湿原などの多様な環境を有した地域である。日本固有種のオオシラビソ（Abies mariesii）は東北地方の亜高山帯を代表する常緑樹となっており，さらに近縁種のシラビソ（Abies veitchii）も北東北には分布していないことからも，航空写真による解析が容易となっている。本研究の結果や方法は，航空写真を入手可
能な多くの亜高山帯で応用することができるため，温暖化によって減少すると考えられているオオシラ
ビソの適域を予測することは，今後の亜高山帯の森
林や生態系の有効な保全対策を考えるうえで重要である。こうした研究の動向を踏まえて，本稿では，八甲田山の亜高山帯に分布しているオオシラビソに焦点を当て，25 m あるいは 50 m のメッシュという高い解像度で，現在までの気候変動の影響と今後の予測を試みる。空中写真を利用して，直近の約 30 年間の個体群動態と現在の分布を分析することにより，変化の予測や潜在的なレフュジアの同定を行った研究を紹介し，取りうる適応策について検討してみたい。
48

2．調査地域の概要（標高経度による八甲田山の植生）
調査地の青森県八甲田山は本州最北部の火山群であり，山岳群の総称となっている。主峰である大岳は標高 1,584 m であり，標高 3,000 m 級の中部地方の日本アルプスなどと比べて標高が低いものの，この山岳群は高緯度（大岳：北緯 40 度 39 分 31 秒）に位置しているために寒冷であり，高山帯と亜高山帯を有している。また，八甲田山は最大積雪深が 3 ～5 m ほどの豪雪地であり，特に 2013 年 2 月には八
甲田山の酸ケ湯（標高 890 m）で 5 m 66 cm の積雪深を 記 録 し， 気 象 庁 の 地 域 気 象 観 測 シ ス テ ム（AMeDAS；Automated Meteor ological Data
Acquisition System）測定開始以来，全国 1 位の積雪記録となっている。八甲田山の酸ケ湯気象観測所では最大積雪深，年平均気温，年間雨量は，それぞれ
3 ～ 5 m，4℃～ 6℃，1,300 ～ 2,300 mm である（気象庁 http://www.jma.go.jp/jma/index.html）。また，調査地である八甲田山の周辺では，1977 年～ 2003 年までに年平均気温は約 0.7 ± 0.5℃上昇した。
八甲田山の植生は，標高 300 m の山麓から標高900 m あ た り ま で ブ ナ 科 ブ ナ 属 の ブ ナ（Fagus
crenata）を主体とする落葉広葉樹が分布し，標高900 m あたりから標高 1,400 m あたりまでマツ科モミ属の常緑針葉樹であるオオシラビソが分布する。
標高 900 m 付近では両者の混交林が形成されている。また，標高 1,400 m を超えるあたりからハイマツ（Pinus pumila）帯が発達した高山帯を形成しており，このような場所では雪田や風衝地も分布し，高山植物群落が見られる。標高 1,300 m 以下には，湿原がパッチ状に見られるが，その多くは更新世（新生代第四紀前期）の大規模な噴火で堆積した火山灰
18），19）。オオシラビのくぼみに発達したものである
ソは連続的な分布域の下限より標高が低い湿原周辺
20）にも分布しており ，湿原周辺が局所的なレフュジ
アとなる可能性がある。なお，この地域ではいくつ
か の 農 場 や 植 林 地 が 存 在 し て い る も の の， 標 高
700 m 以上の土地の大半は国立公園に指定されており，自然植生への人間の影響は，特に亜高山帯では極めて小さい。
3．高解像度航空写真を用いたオオシラビソ個体群の 30 年後との比較
ここでは，過去および最近の空中写真を利用した，オオシラビソ個体群の 30 年後との比較と，そ21）の変化を説明する要因の解析 について述べる。
解 析 に は， 林 野 庁 に よ っ て 1967 年 9 月 および2003 年（一部 1996 年）に撮影された空中写真を用いている。広葉樹の紅葉時期に写真撮影が行われたため，視覚的に容易に常緑のオオシラビソを識別することができた。すべての航空写真は 0.5 m グラウン

地球環境 Vol.19 No.1

おける樹冠数の差，平均樹冠高を用い，ランダム効
果として 3 次メッシュセル（約 1 km × 1 km）内の50 m メッシュセルを入れ子にすることで，空間自己相関の影響がなるべく出ないようにした。
その結果，過去と 30 年後の樹冠数の変化を最も良く説明できる変数は，標高，斜面方位，湿原からの距離，過去の樹冠数であった（図 1）。図 1 は，縦軸に目的変数として個体数変化を表す樹冠数をとり，横軸に各説明変数をとったものであるが，オオ
シラビソの樹冠数は，過去の樹冠数の多い場所で減少する（P ＜ 0.001，図 1（d））
。また，標高 1,000 m 以
下で減少，標高 1,300 m 以上で増加し（P ＜ 0.001，
図 1（a）），斜面方位ではやや南東斜面で減少するほか（P ＜ 0.05，図 1（b）），湿原からの距離が遠いほど直線的に減少していた（P ＜ 0.01，図 1（c））。
30 年後における樹冠高の変化（成長量）は，斜面
方位，日射指数，湿潤度，湿原からの距離が最もよく説明できる変数であった（図 2）。南東斜面で成長量が小さく（P ＜ 0.01，図 2（a）），太陽放射の強い場
所で大きい傾向があった（P ＜ 0.001，図 2（b））。湿
潤度は有意な相関を示したものの，顕著なパターンは検出されず（P ＜ 0.01，図 2（c）），湿原の近くでは
ゆっくりと成長する傾向が認められた（P ＜ 0.001，
図 2（d））。

s（標高, 5.28）s（斜面方位, 3.17）

ドサンプル距離（GSD；Ground Sample Distance）で幾何補正（オルソ補正）し，SOCET SET ver. 5.0.0（BAE
Systems, London, UK）を用いて 1 m GSD デジタル
表面モデル
（DSM；Digital Surface Model）を作成した。 オ ル ソ 補 正 の た め，1 m 以 下 の 位 置 精 度 で
SOKKIA GIR1600 ディファレンシャル GPS 受信機（株式会社トプコン，日本）を用いて，調査地域内に14 の地上基準点を設けている。オオシラビソが現在分布している地域全体にわたって，50 m メッシュ解像度で 178 の区画のセルをランダムサンプリングした。各 50 m のメッシュセルを，四つの 25 m のメ
ッシュセルに分割し（全体で 712 プロット），立体視
しながら区別可能なオオシラビソの林冠数と平均林冠高を測定した。
環境要因として，標高，傾斜，斜面方位，曲率，
日射指数，および湿潤度を，国土地理院の 1：25,000
地形図から作った 10 m 標高モデル（DEM；Digital
Elevation Model）から，ArcGIS ver. 9.2 ソフトウェ
ア（ESRI, Redlands, USA）の空間解析ツールを使用
して生成した。10 m メッシュデータは，25 m メッ
シュデータに平均化した。また，各セルから最も近い湿原にユークリッド距離を追加した。個体群動態と 環 境 要 因 の 関 係 は， 一 般 化 加 法 混 合 モ デ ル（generalized additive mixed ef fect models；
22）
GAMMs） を使用した。目的変数として，2 期間に

47−55（2014）

斜面方位 (°)

s（過去の樹冠数, 1）s（湿原からの距離, 1）

標高 (m)

湿原からの距離 (m)図1

過去の樹冠数

約 30 年間におけるオオシラビソの樹冠数の変化に関する偏残差プロット．

縦軸は（a）: 標高，（b）: 斜面方位，（c）: 湿原からの距離，（d）: 過去の樹冠数に対する偏残差（他の効果を取り除いたときの21）．
特定の説明変数の効果）

49

s（日射指数, 3.91）s（斜面方位, 3.97）

田中ほか：気候変動が森林動態に与える影響と将来予測

日射指数s（湿潤度, 2.6）s（湿原からの距離, 2.27）

斜面方位(°)

湿潤度
図2

湿原からの距離 (m)約 30 年間におけるオオシラビソの樹高変化に関する偏残差プロット．

縦軸は（a）: 斜面方位，（b）: 日射指数，（c）: 湿潤度，（d）: 湿原からの距離に対する偏残差（他の効果を取り除いたときの21）．
特定の説明変数の効果）

4．気候変動後のオオシラビソ個体群の分布予測
ここでは，空中写真を用いたオオシラビソの分布
環境の局所モデルと，温暖化にともなう分布変化予23）
測 について述べる。現在までの潜在的な適域の解
析は，通常地球規模の粗い気候データを用いて 1 ～50 km 程度の解像度で行われていたために，微地形や局所的な気候の影響を考慮できず，局所モデルによる解析が必要とされていた。特に，湿原周辺がオ
オシラビソの局所的なレフュジアとなる可能性に注
目し解析を行う。
調査地と分布データは，八甲田山の広範囲にわたる 20 km × 15 km の範囲で（図 3），解像度は 50 m
（117,261 セル）と 1 km
（294 セル）メッシュの 2 段階で解析した。多くの場合，1 km メッシュでの解析が行われているが，これを基に地域または局所的なスケールでの実用的な温暖化に対する適応策を考えるには不十分であり，加えてこの解像度に地形のような地理的要因の影響を組み込むことは難しい。そこで，より実用的で詳細な解析を行うために解像度を上げて，50 m メッシュを用いて解析を行った。
オオシラビソの分布図は最近の航空写真（2003 年および 1996 年に林野庁が撮影）に基づいて作成した。
人為的に影響を受けたエリアは除外し，50 m メッシュのグリッド上の座標で表現されるラスタデータ
50に変形し，それを 1 km メッシュの解像度にアップスケールした。
環 境 要 因 と し て， 暖 か さ の 指 数（WI；warmthindex）と地形および地形環境因子に関する 6 変数
（標 高， 傾 斜， 斜 面 方 位， 曲 率， 日 射 指 数，湿潤度），湿原に関する 2 変数（湿原からの距離，セル内での湿原の有無）を用い，Decision Tree Model を両
解像度（50 m，1 km メッシュ）で作成した（図 4）。
四つの温度上昇シナリオ（＋ 1℃，＋ 2℃，＋ 3℃，＋ 4℃）と湿原に関する二つのシナリオ（気候変動に合わせて縮小する，縮小しない）を組み合わせた計
8 シナリオに対し，オオシラビソの潜在生育域を予測した。湿原の縮小シナリオは，オオシラビソの分布に関する解析と同様に Decision Tree Model をあてはめ，それぞれの温度域での湿原分布を推定した（図 5）。
1 km メッシュモデルでは，WI が唯一の説明変数だった（図 4（a））。50 m メッシュモデルでも WI が最も重要な説明変数であり，二番目に重要な説明変数は傾斜であった。斜面が穏やかなほどオオシラビ
ソの分布確率が高い傾向が示され，オオシラビソの分布は斜面上の雪が地面に沿って移動するプロセス
24）（匍行）に関係している可能性がある 。雪の匍行
は，傾きが 15 度以上の斜面で起こり，樹木に大きな圧力がかかる。50 m メッシュモデルにおける傾地球環境 Vol.19 No.1

斜の分岐は 17.5 度であり，雪が滑る圧力がオオシ
ラビソの分布に影響を与えることが予想される。
三番目に重要な説明変数として斜面方位が選択され，オオシラビソの分布確率は東斜面で低かった。
北西からの冬の季節風によって形成される雪だまり
25）の影響があると考えられる 。
50 m メッシュでの予測結果を見ると，＋ 1℃～

47−55（2014）＋ 2℃のシナリオでは，オオシラビソに適した生育
適地はより高所に移行すると同時に，生育適地が二つに分断されることが予測された。さらに，＋ 4℃の上昇では，生育適地がほぼ消失する。湿原が縮小するシナリオと縮小しないシナリオの差は大きくなかったが，どの温度域でも，オオシラビソの生育適地の減少は湿原縮小シナリオのほうが大きいことが標高 (m)

図3

オオシラビソの分布．

（a）オオシラビソの水平分布．中心部の実線がオオシラビソの分布の外周を示す．斜線部は人為的な改変を受けた場所として解析から除外した．色の濃さは標高を示す．（b）標高に応じた分布．黒い点はオオシラビソが存在したセルを示す 23）．

(a)

(b)
WI ＜ 50.7

Yes

No

WI ＜ 46.5

Yes

0.323

(n=226)(n=89175)

1.000(n=40)A

WI ＜ 50.5

傾斜角度 ＜ 17.5°

WI ＜ 44.7
0.014

0.835

(n=141)C

0.689

WI＜ 44.6

湿原の有無 ：無 0.632
0.693

0.736

0.315

(n=13795) (n=1605)
AE(n=45)
B

図4

No

0.204

0.529 斜面方位 :NE E SE

0.437

0.346

0.138

WI ＜ 43.3
0.464

0.026

湿原から
の距離 ＜ 225m
0.102

0.001

(n=4814) (n=4127) 0.561
0.127 0.384
0.068 (n=43294)
C(n=5656) (n=1633) (n=1530) (n=12721)
F
I
B
G
D
H

オオシラビソの分布に関する Decision tree．

（a）1 km メッシュサイズで解析した場合．（b）50 m メッシュで解析した場合．変数名は以下のとおり．
暖かさの指数（WI），傾斜角度，斜面方位，湿原の有無，湿原からの距離 23）．例として（a）では，暖かさの指数が 50.7 未満であり，かつ暖かさの指数が 44.7 より小さい分布の可能性は 1.000 である．また，暖かさの指数が 50.7 未満であり，かつ暖かさの指数が 44.7 以上の分布の可能性は 0.689 である．また，各条件で分けた終端ノードを A ～ I のアルファベットで表した．

51

田中ほか：気候変動が森林動態に与える影響と将来予測
傾斜角度

WI

WI
傾斜角度

傾斜角度

傾斜角度

WI

図5

WI

湿原の分布に関する Decision tree．変数名は，暖かさの指数（WI），傾斜角度 23）．
各条件で分けた終端ノードを A ～ I のアルファベットで表した．

1℃

温度上昇

WI

1.5℃

2℃

2.5℃

3℃

湿原の縮小がないシナリオ

分布域外
分布可能域

湿原の縮小を考慮したシナリオ

分布適域

図6
表1

温度上昇と湿原の縮小シナリオごとのオオシラビソの分布変化 23）．

異なる色は，分布確立の違いを表す．

シナリオごとのオオシラビソの分布適域面積．
湿原が縮小しないシナリオ
現在

＋ 1℃

＋ 2℃

＋ 3℃

＋ 4℃

＋ 1℃

＋ 2℃

＋ 3℃

＋ 4℃

分布適域（ha）

645.4

311.7

126.1

1.5

0.0

284.8

76.4

1.2

0.0

分布可能域（ha）

506.1

316.5

128.2

24.1

0.0

343.5

177.9

24.4

0.0

分布不適域（ha） 1780.1

2303.32677.32905.92931.52303.32677.32905.92931.5分かった（図 6）。例えば，オオシラビソに適した生育地は湿原が縮小しない＋ 2℃のシナリオで 126.1 ha
であるのに対し，湿原が縮小する＋ 2℃のシナリオ
では 76.4 ha と差が大きくなる場合もある（表 1）。
解像度の違いは全体のパターンに大きな影響はな
いが，特に＋ 2℃のシナリオにおいて，1 km メッシュモデルはオオシラビソの適域を部分的に過小評価していた。50 m メッシュモデルで予測したいく
52

湿原が縮小するシナリオ

つかの生育適地は，1 km メッシュモデルが予測した適地の外側に位置していた。これは，おそらく1 km メッシュモデルでは使用されなかった変数の効果によるものである。つまり，粗い解像度のモデ
ルでは，微地形などに関連する種の適応的機能を見落とし，気候変動後の生存分布を正確に予測できないことを示唆している。

地球環境 Vol.19 No.1

5．考察
今回の解析では，この数十年間にオオシラビソの個体群の個体数は低標高で減少し，高標高で増加していることが明らかになった（図 1（a））。それは個体群の移動がその分布範囲内で進行中であることを示唆しており，気候変動の影響が表れている可能性がある。一方，花粉分析では，八甲田山のオオシラ
ビソ林はわずか 600 年ほど前に形成されたことが示26）
唆されている 。また，温度条件だけから推定する
27）分布の上限は，現在の分布よりもはるかに高い 。
これらの点を考えると，高所への分布拡大は，単に潜在的な分布域全体にまだオオシラビソが拡大しきっていない可能性など，気候変動以外の理由による可能性もある。しかし，高標高地の温度上昇と降雪量の減少が，これまでオオシラビソの生育に適さなかった場所への分布拡大を促進しているという報告
28）もあり ，積雪など温度以外の要因の変化が影響している可能性も否定できない。一方で，低標高域での個体群衰退は，温度上昇が主な原因である可能性が高い。暖かい環境では，通常ブナがオオシラビソ
20）に競争で勝つことがわかっている 。
湿原近くのオオシラビソの増加（図 1（c））や，温暖化による分布変化の予測結果は，湿原の周囲が潜在的なレフュジアとなり得ることを示唆していた。
オオシラビソは，より温暖環境で優占種となるブナに比べて，湿潤（そしておそらくより酸性）な土壌環
29）
境に適応した形質を持っている 。ヨーロッパブナ（Fagus sylvatica）とノルウェートウヒ（Picea abies）
30）との間にも同様の関係が観察されている 。湿原周辺部では，ブナの成長は強酸性の鉱質土壌によって阻害される可能性がある。八甲田山においても湿原周辺ではブナの成長が抑制され，その結果として，
オオシラビソに潜在的な生育地が提供された可能性がある。湿原そのものの分布条件の解析によって，温暖化による湿原の衰退が予想されるが，湿原の富栄養化のような人為的要因を除去することにより湿原の衰退を抑制することが，オオシラビソ林を保全する適応策としても有効である。
気候変動の局所スケールの影響は十分に検討され31）ていない 。なぜなら気候データは，通常粗い（例えば，50 km や 100 km の）地球規模での気候の変化を表現する全球気候モデル（GCM；Global Climate
Model）を用いており，また温度が均一に上昇していくことを前提とした理想的シナリオであることに32）よる 。本研究では，今までにほとんど用いられていない地域スケールの気候変動シナリオについて
50 m メッシュモデルを使用して，地球温暖化の樹
木の分布への影響を予測した。この解像度での解析により，1 km メッシュモデルで検出することができなかった局所的な適域を特定することができた。
ヨーロッパでも，局所スケールのモデル（25 m メッ

47−55
（2014）

シュ解像度）は，粗い分解能（16 km メッシュ）のモ
デルによって予想するのに不向きな種の生息域を予33）測できることが示唆されている 。今回，八甲田山のオオシラビソ林にとって湿原の周囲が分布の下限に近い局所的な適地として識別されたことは，気候要因に加えて他の説明変数がなければ実現しなかった。高解像度のモデルには，地形や湿原からの距離といったいくつかの非常に局所的な環境勾配の影響が重要となる可能性があるが，それらを低解像度のモデルに組み込むことはできない。それぞれの種に応じた高解像度モデルに様々な要因を取り入れることで，より多くの種の潜在的なレフュジアを見出すことができるだろう。
本研究から，オオシラビソ林の動態パターンは，全体的には解像度のレベルによる影響は少ないが，粗い解像度では見落とされる生育適地があることが明らかとなった。具体的なレフュジアの特定や適応策などの計画を考えると，より解像度の高い予測が必要である。しかし，高解像度では降水量や積雪量などの予測の難しい気候条件があり，高解像度の気候変動予測モデルを特定の地域に対して開発するな
どの手段が重要になるだろう。今回の解析方法では，常緑樹というオオシラビソの利点によって空中写真から個体を識別でき，解像度の高いデータを生成できた。こうした高解像度のデータは，生物や生態系の情報としても重要である。
謝

辞

本研究は，環境省地球環境研究総合推進費（D0904），三井物産環境基金，および文部科学省気候変動リスク情報創生プログラム（D-iii）の助成を受けて行った研究をまとめたものである。
引用文献
1）Intergovernmental Panel on Climate Change（IPCC）（2007）Climate Change 2007: Impacts, Adaptation
and Vulnerability. Contribution of Working Group II
tothe

FourthAssessment

ReportoftheIntergovernmental Panel on Climate Change. Parry,
M.L., O.F. Canziani, J.P. Palutikof, P.J. van der Linden
and C.E. Hanson, eds., 976pp, Cambridge University
Press, Cambridge, UK.
2）Overpeck, J. T., D. Rind and R. Goldberg（1990）
Climate-induced changes in forest disturbance andvegetation. Nature, 343, 51-53.
3）Gonzalez, P., R. P. Neilson, J. M. Lenihan and R. J.
Drapek（2010）Global patterns in the vulnerability of
ecosystems to vegetation shifts due to climatechange. Global Ecology and Biogeography, 19, 755
-768.

53

田中ほか：気候変動が森林動態に与える影響と将来予測

4）Parmesan, C.（2006）Ecological and evolutionar yfor beech（Fagus crenata）forests under 110 hypo-

responses to recent climate change. Annual Review

thetical climatic change scenarios in Japan. Appliedof Ecology Evolution and Systematics, 37, 637-669.

Vegetation Science, 12, 328-339.

5）Peñuelas, J., R. Ogaya, M. Boada and A. S. Jump

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

（2007）Migration, invasion and decline: changes in

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

recruitment and forest structure in a warming-linked

Probability distributions, vulnerability and sensitivity

shift of European beech forest in Catalonia（NEin Fagus crenata forests following predicted climate

Spain）. Ecography, 30, 829-837.

changes in Japan. Journal of Vegetation Science, 15,6）Wang, T., Q.-B. Zhang and K. Ma（2006）Treeline605-614.

dynamics in relation to climatic variability in the18）Muraoka, H. and S. Takakura（1988）Explanator ycentral Tianshan Mountains, northwestern China.

text of the geological map of Hakkoda Geothermal

Global Ecology and Biogeography, 15, 406-415.

Area. Geological Survey of Japan, Tsukuba, Japan.

7）Kullman, L.（2007）Tree line population monitoring19）小池一之・田村俊和・鎮西清高・宮城豊彦（2005）of Pinus sylvestris in the Swedish Scandes, 1973-

日本の地形 3，東北．東京大学出版，東京．

2005: implications for tree line theor y and climate

20）山中三男・菅原 啓・石川慎吾（1988）南八甲田山の山地帯にみられるアオモリトドマツ林の変遷．日change ecology. Journal of Ecology, 95, 41-52.
8）Batllori, E. and E. Gutiérrez（2008）Regional treeline dynamics in response to global change in the

本生態学会誌，38，147-157．
21）Shimazaki, M., T. Sasaki, K. Hikosaka and T.
Nakashizuka （2011） Environmental dependence

Pyrenees. Journal of Ecology, 96, 1275-1288.
9）Lenoir, J., J. C. Gégout, P. A. Marquet, P. de Ruffrayof population dynamics and height growth of aand H. Brisse（2008）A significant upward shiftsubalpine conifer across its vertical distribution:in plant species optimum elevation during the 20thancentury. Science, 320, 1768-1771.

photographs. Global Change Biology, 17, 3431-3438.

approachusing

high-resolution

aerial10）Kadmon, R. and R. Harari-Kremer（1999）Studying22）Wood, S. N.（2006）Low-Rank Scale-Invariant Tensor

long-term vegetation dynamics using digital process-Product Smooths for Generalized Additive Mixeding of historical aerial photographs. Remote Sensing

Models. Biometrics, 62, 1025-1036.
23）Shimazaki, M., I. Tsuyama, E. Nakazono, K. Nakao,

of Environment, 68, 164-176.
11）Witt, G. B., J. Luly and R. J. Fairfax（2006）How

M. Konoshima, N. Tanaka and T. Nakashizuka

the west was once: vegetation change in south-west（2012）Fine-resolution assessment of potential

Queensland

refugia for a dominant fir species（Abies mariesii）of

from

1930

to

1995.

Journalof

Biogeography, 33, 1585-1596.
12）Beckage, B., B. Osborne, D. G. Gavin, C. Pucko, T.

subalpine coniferous forests after climate change.
Plant Ecology, 213, 603-612.

Siccama and T. Perkins（2008）A rapid upward shift24）Leitinger, G., P. Höller, E. Tasser, J. Walde and U.

of a forest ecotone during 40 years of warming in the

Tappeiner（2008）Development and validation

Green Mountains of Vermont. Proceedings of theof a spatial snow-glide model. Ecological Modelling,

National Academy of Sciences, 105, 4197-4202.

211, 363-374.

13）Iverson, L. R. and A. M. Prasad（1998）Predicting25）四手井綱英（1956）裏日本の亞高山地帯の一部に針

abundance of 80 tree species following climate

葉樹林帯の欠除する原因についての一つの考えか

change in the eastern United States. Ecological

た．日本林学会誌，38，356-358．

Monographs, 68, 465-485.

26）Morita, Y.（2000）The vegetation histor y of the14）Horikawa, M., I. Tsuyama, T. Matsui, Y. Kominami

subalpine zone in Japan since the Last Glacial period:and N. Tanaka（2009）Assessing the potential

Were the forest zones higher than they are at

impacts of climate change on the alpine habitat

present during the Climatic Optimum period？

suitability of Japanese stone pine（Pinus pumila）.
Landscape Ecology, 24, 115-128.
15）Trivedi, M. R., P. M. Berry, M. D. Morecroft and T. P.

日本の植生．科学，61，36-49．

Dawson（2008）Spatial scale affects bioclimate model

28）Daimaru, H. and H. Taoda（2004）Effect of snow

projections of climate change impacts on mountainpressure on the distribution of subalpine Abies

plants. Global Change Biology, 14, 1089-1103.

mariesii Forests in northern Honshu Island, Japan.

16）Matsui, T., K. Takahashi, N. Tanaka, Y. Hijioka, M.

54

Japanese Journal of Historical Botany, 9, 3-20.
27）野上道男・大場秀章（1991）暖かさの指数からみたJournal of Agricultural Meteorology, 60, 253-261.

Horikawa, T. Yagihashi and H. Harasawa（2009）29）Sugita, H.（1992）Ecological geography of the range

Evaluation of habitat sustainability and vulnerabilityof the Abies mariesii forest in northeast Honshu,

地球環境 Vol.19 No.1

47−55
（2014）

Japan, with special reference to the physiographic
conditions. Ecological Research, 7, 119-132.
30）Murach, D. and B. Ulrich（1988）Destabilizationof forest ecosystems by acid deposition. GeoJournal,

黒川

17, 253-259.
31）del Barrio, G., P. A. Harrison, P. M. Berry, N. Butt,
M. E. Sanjuan, R. G. Pearson and T. Dawson（2006）
Integrating multiple modelling approaches to predictthe potential impacts of climate change on species'
distributions in contrasting regions: comparison andimplications for policy. Environmental Science &Policy, 9, 129-147.
32）Beaumont, L. J., A. J. Pitman, M. Poulsen and L.
Hughes（2007）Wher e will species go ？
Incorporating new advances in climate modelling
into projections of species distributions. Global
Change Biology, 13, 1368-1385.

紘子

Hiroko KUROKAWA

1976 年，広島県生まれ。東北大学大学院生命科学研究科助教。横浜市立大学理学部卒。京都大学大学院理学研究科博士課程修了。専門は植物生態学。マレー
シアやタイの熱帯林で人為撹乱に対する森林群集や生態系機能の応答，ニュージーランドで外来植物が生態系に与える影響，青森県八甲田で気候変動に対する生態系の物質循環の応答などを研究。主な著書に，『外来生物－生物多様性と人間社会への影響－，2 章：外来植物による生態系改変とその影響－機能形質に着目したアプローチ－』（裳華房），『生態適応科学－自然のしくみを活かし，持続可能な未来を拓く－，1 章：生態系の適応力』（日経ビジネス，
編著）
など。

彦坂

幸毅

Kouki HIKOSAKA33）Randin, C. F., R. Engler, S. Normand, M. Zappa,
N. E. Zimmermann, P. B. Pearman, P. Vittoz, W.
Thuiller and A. Guisan（2009）Climate change andplant distribution: local models predict high-elevation persistence. Global Change Biology, 15,1557-1569.

田中孝尚

Takahisa TANAKA

1970 年，埼玉県生まれ。2006 年，東北大学大学院生命科学研究科博士課程修了，博士（生命科学）。専門は植物生態学，
集団遺伝学，考古植物学。八甲田山やヒマラヤなどの環境の厳しい山岳地帯で標高傾度と植物との関係，有用植物の栽培化や起源などの植物と人間との関係を研究している。

嶋崎仁哉

Masaya SHIMAZAKI

東北大学大学院生命科学研究科生態シ
ステム生命科学専攻博士課程前期 2 年修了。気候変動に対する山岳生態系の脆弱性評価を空中写真をもとに行う。

東北大学理学部卒業，東京大学大学院
理学系研究科修了。博士（理学）。2010
年より東北大学生命科学研究科教授。植物の生態について，分子から生態系まで幅広い興味で研究を行っている。特に興味があるのは光合成系の適応と順化，環境適応の種間・集団
間変異，植物形質の適応的意義など。野外・実験・理論といった異なるアプローチの融合を目指している。

中静透

Tohru NAKASHIZUKA

1956 年，新潟県生まれ。1978 年千葉大学卒業。1983 年理学博士（大阪市立大学）。
森林総合研究所主任研究官，国際農林水産業研究センター主任研究官，京都大学
生態学研究センター教授，総合地球環境学研究所教授を経て 2006 年より東北大学生命科学研究科教授（現職）。専門は森林生態学，生物多様性科学で，熱帯林および温帯林の動態と更新，林冠生物学，森林の持続的管理と生物多様性などを研究。国際生物多様性研究計画（DIVERSITAS）
科学委員，国際森林研究機関連合（IUFRO）
理事，日本生態学
会長などを歴任。主な著書に，『モンスーンアジアの生物多様性』（岩 波 書 店， 共 著 ），『Diversity and Interaction in a
Temperate Forest Community. Ogawa Forest Reser ve ofJapan』（Springer-Verlag，共編・共著），『森のスケッチ』（東海
大学出版会）など。

55