まず、2003年から2014年にかけて、京都大学芦生研究林のモンドリ谷に設置した16haの固定プロットにおいて、304本のミズナラの被害状況を確認しました。前年の被害木からのカシノナガキクイムシの分散カーネルが対数正規分布に従うと仮定し、平均と標準偏差を変化させることでピーク位置とばらつき具合の異なる100種類の仮想分散カーネルを準備しました。前年被害木からの距離に応じてカシノナガキクイムシの移動分散確率を計算し、どのような分散カーネルを仮定した場合に被害発生確率を予測するモデルのあてはまりがよくなるのかを探索しました。
調査地では2004年に初めてナラ枯れ被害が確認されたのですが、以後2011年までカシノナガキクイムシの穿孔による枯死がみられました。各年各個体の被害の有無を応答変数とした一般化線形モデルを構築したところ、説明変数の中では胸高直径が最も逸脱度を説明した割合が高く、次に高かったのがカシノナガキクイムシの移動分散確率でした。前年の被害木から300mにピークがあり、比較的ばらつき具合の小さい分散カーネルを仮定した場合が最もモデルのあてはまりがよくなり、直近には飛ばないカシノナガキクイムシの移動分散様式が示唆されました。また、最適モデルには、胸高直径と移動分散確率以外に周辺10mのミズナラの密度も説明変数として採択され、局所的に集中分布する太いミズナラがカシノナガキクイムシに攻撃されやすいことが明らかとなりました (Yamasaki et al., 2016)。
