ブナ科堅果の結実豊凶がツキノワグマの個体群動態に及ぼす影響の解明

本研究では、ツキノワグマ（以下、クマ）の個体群動態が、食物資源量の年変動によって受ける影響を解明することを目的とする。令和3年度は、令和2年度分のデータに追加で取得した情報を加え、クマの出生と死亡に関わる生活史パラメータの推定と、ブナ科堅果類の結実豊凶による出生に関する繁殖パラメータへの影響の検証を行った。2005～2019年に有害捕獲された4歳以上のメス132個体の歯の年輪幅を用いて初繁殖成功年齢、繁殖成功間隔を算出した。それぞれ、5.44±0.22歳（平均±標準誤差）、2.38±0.48年と推定された。2009～2019年に群馬県で有害捕獲されたメス88個体の子宮の胎盤痕を用いて一腹産子数を算出した。1.58±0.09頭と推定された。さらに、2003～2021年に栃木県日光足尾山地にて学術捕獲されたメス43個体の情報を用いて、mark recapture recovery modelにより自然死亡率と人為死亡率を推定した。それぞれ10.8%（95%信頼区間：6.4－17.4）、0.5%（0.2－2.1）と推定され、人為死亡率が低いという結果となった。クマと人間との軋轢が少ない本調査地の傾向が反映された可能性がある。また、同調査地における2000～2019年までの、メス10個体の0歳の子連れの有無と子の生存状況の観察記録に基づき、一般化線形モデルにより出生後半年における死亡率を推定した。23.5%（95%CI：0.8－46.5）と推定された。2004～2019年に長野県、群馬県と栃木県で取得されたブナ科4樹種の堅果の結実状況の記録から、カテゴリーベース（豊作または凶作）とエネルギーベースの豊凶指標を算出した。豊凶指標が各繁殖パラメータへ及ぼす影響を一般化線形モデル及び一般化線形混合モデルにより推定したが、いずれのパラメータについても顕著な影響は見られなかった。

本研究では、ツキノワグマ（以下、クマ）の個体群動態が、食物資源量の年変動によって受ける影響を解明することを目的とする。令和3年度は、令和2年度分のデータに追加で取得した情報を加え、クマの出生と死亡に関わる生活史パラメータの推定と、ブナ科堅果類の結実豊凶による出生に関する繁殖パラメータへの影響の検証を行った。2005～2019年に有害捕獲された4歳以上のメス132個体の歯の年輪幅を用いて初繁殖成功年齢、繁殖成功間隔を算出した。それぞれ、5.44±0.22歳（平均±標準誤差）、2.38±0.48年と推定された。2009～2019年に群馬県で有害捕獲されたメス88個体の子宮の胎盤痕を用いて一腹産子数を算出した。1.58±0.09頭と推定された。さらに、2003～2021年に栃木県日光足尾山地にて学術捕獲されたメス43個体の情報を用いて、mark recapture recovery modelにより自然死亡率と人為死亡率を推定した。それぞれ10.8%（95%信頼区間：6.4－17.4）、0.5%（0.2－2.1）と推定され、人為死亡率が低いという結果となった。クマと人間との軋轢が少ない本調査地の傾向が反映された可能性がある。また、同調査地における2000～2019年までの、メス10個体の0歳の子連れの有無と子の生存状況の観察記録に基づき、一般化線形モデルにより出生後半年における死亡率を推定した。23.5%（95%CI：0.8－46.5）と推定された。2004～2019年に長野県、群馬県と栃木県で取得されたブナ科4樹種の堅果の結実状況の記録から、カテゴリーベース（豊作または凶作）とエネルギーベースの豊凶指標を算出した。豊凶指標が各繁殖パラメータへ及ぼす影響を一般化線形モデル及び一般化線形混合モデルにより推定したが、いずれのパラメータについても顕著な影響は見られなかった。