山形県におけるクマタカの生息適地推定モデルの構築

Abstract

本研究では，希少猛禽類クマタカの保全・保護対策のために, クマタカの生息適地推定モデルを構築し山形県内における生息適地分布図を作成した．本モデルでは対象地域をグリッドで分割し, 個々のグリッドに対して周辺の地形・植生状況を表現する各種指標値を与えて生息地としての適否を推定した. 従来の研究では, 約5km×5kmのグリッド内の, 平均標高・平均傾斜・最低標高等の地形指標が主に用いられていた. しかし本研究では, クマタカの生息環境をより的確に表現するため, クマタカの生息に重要とされる谷地形の分布状況を表現する指標等, 全部で13種類の地形・植生指標を対象とした. また, 詳細な生息適地分布図を作成するためグリッドサイズを50mとした. なお, 各指標が集計対象とする周辺範囲の広さによって指標の環境説明力が変化するため，既往の研究を参考に, 集計対象範囲を各グリッドから半径0.5km, 1km, 1.5km, 3kmの円内の4通りに設定した. 各指標値の算出には国土地理院の数値地図50mメッシュ(標高)および環境省の現存植生図(1/50, 000)を用いた. 最適モデルの選択においては，ロジスティック回帰分析を用いて各地形・植生指標を組み合わせたモデル式を作成し, 各々の適合度や推定精度等を比較した. その結果, 周辺半径1kmの谷地形の分布状況を表す指標と周辺半径3kmの樹林地の面積割合を表す指標を組み合わせたモデルが最適モデルとして選択された(Overall accuracy 91.8％). さらに, 最適モデルにもとづいて作成した生息適地分布図を5kmグリッド化したものとクマタカ確認情報分布図を比較すると, 概ね高い精度(Overall accuracy 89.0%)で一致し, 推定された生息適地の大部分においてクマタカの確認情報が存在していることが判明した. 以上より, 本研究で作成したモデルは高い精度で生息適地を推定できており, クマタカの生息適地推定には谷地形の分布と樹林面積に着目することが重要であるといえる.

The mountain hawk-eagle (Spizaetus nipalensis) is an endangered and an umbrella species ranking high in the food chain of the forest ecosystem. The objective of this study was to estimate the potential habitat of the mountain hawk-eagle in Yamagata prefecture, northern Japan. Habitat suitability of each grid was estimated using values of indices corresponding to terrain and vegetation features of the surroundings. Previous studies have computed terrain indices such as minimum and mean altitude and mean slope gradient within each grid-square of approximately 5km in size. To estimate potential habitat more accurately, we used 13 terrain and vegetation indices, including the valley distribution index that is important for the mountain hawk-eagle, and set a grid resolution of 50m to describe potential habitat map at fine scale. Since these indices would be sensitive to the analyzing window size placed around each grid, we examined this effect by changing radii as 0.5, 1, 2 and 3km. Values of these indices were computed in each grid with a 50m digital elevation model (50m DEM) and a rasterized 1/50, 000 digital vegetation map. By comparing Akaike’s Information Criterion (AIC) values and accuracies among logistic regression models with different combinations from a set of 13 indices, an optimal model with the lowest AIC value and high accuracy was determined. The model used the distribution of valleys within 1km radius and the area of forest within 3km radius, which indicates the importance of valleys and forest for the mountain hawk-eagle’s habitat. The potential habitat areas estimated by this model substantially corresponded (89.0%) to the areas where the mountain hawk-eagle has been observed, with the overall accuracy of 91.8%. We concluded that this model could accurately estimate the potential habitat of the mountain hawk-eagle.