2019年夏期ヤクシマザル調査報告
2019.9.23. 2019ヤクザル調査隊
1 はじめに
ヤクシマザル(Macaca fuscata yakui)は屋久島の亜熱帯性海岸林から標高1800mを越える山頂部まで連続的に分布することが知られている (丸橋他、1986; Hanya et al., 2004)。屋久島の西部海岸域では、1970年代から餌付けに頼らない長期にわたる継続調査が行われており、 多くの成果が得られている(Yamagiwa & Hill, 1998)。中高度・上部域に生息するヤクシマザルについては調査が遅れていたが､ 1990年代に入って、西部海岸部以外での調査が本格化し､上部域のニホンザルの密度､群れサイズ､食性､活動時間配分など、 基本的な生態学的情報が、すこしずつ明らかになってきた(Hanya et al., 2003a; Hanya, 2004a, 2004b)
「ヤクザル調査隊」は1989年以降、毎年夏にヤクシマザルの分布調査を行ってきた。このうち、1994、1995、1997年には、 自然植生の垂直分布が連続して保存されている西部域で調査を行い、以下のことが明らかになった。 (1)標高300mまでの海岸部には、1km2あたり約4群と、きわめて高密度で分布している。(2)標高300mから800mの照葉樹林帯、 標高800mから1,200mの照葉樹林・ヤクスギ林移行帯、標高1200m以上のヤクスギ帯では、群れ密度は1km2あたり 1.2群から1.6群と、ほぼ一定である。(3)群れサイズについては正確な情報は得られなかったが、標高800m以上の移行帯およびヤクスギ帯では、 それ以下の標高帯に比べ小さいようであった(Yoshihiro et al., 1999; Hanya et al., 2004)。
さらに1998年からは、西部域の大川林道終点地域で継続調査を行うことにした。具体的には、調査域内の群れの正確な分布と、 対象とした群れのサイズと性年齢構成を明らかにすることを目的として行われた。これらの資料を蓄積することによって、屋久島中高度域 におけるヤクシマザルの個体群動態のメカニズムを明らかにすることが最終的な目的である。1998年の調査では、大川林道終点地域の、 瀬切川の両岸で調査を行い、1群の構成を確定したほか、3群について、大まかな構成を確認した。群れサイズはどれも20頭を越えるものではなかった。 2000年には、HR群を全頭を個体識別し、ほかにもPE、OM、SY群の合計4群が識別できるようになった。2001年の調査では、 さらにBR群が識別されるようになった。BR群は2003年以降は調査域内で観察されなくなったが、2005年にはSY群がSS群、 YY群の二つの群れに分裂し、現在調査域内では5群を識別している。これらの群れの子持ち率(または粗出産率、 観察されたオトナ・ワカモノメスのうち、アカンボウを持っていたメスの割合)は年によって4%-43%のあいだで変動した。 2000年から、これら人口学的資料の収集と並行して、ブロック分割定点調査法による集団密度の推定を行った。 この地域は自然林と伐採地が混在しており、しかも伐採された年数が異なっている。このような地域で集団密度の推定を行うことで、 ニホンザルが植生の撹乱の程度に応じてどのように土地利用のやり方を変化させていくのかを明らかにし、さらにその調査を毎年継続して行うことで、 伐採後の植生の変化に応じて土地利用のあり方がどう変化していくかを明らかにして行くことが目的である。2000年から2003年までの結果については、Hanya et al. (2005)に発表した。
さらに2003年度からは､定点調査中にヒトリザルおよびシカについてもサルの集団と並行して資料を収集するようにした｡ シカ調査は、2008年までは定点で音声が聞こえる頻度を記録していたが、この方法で得られた結果の信頼性に疑問があったため、2008年度からは、糞塊法による調査を開始した。 定点の周辺にシカ糞調査用の調査区を設置し、シカ密度の年変動を恒久的に調査できるようになった。サルの分布と関連すると思われる、 森林の果実量調査は、自然林では1999年から、伐採地では2002年から継続して行っている。また、2010年から、定点でのヒルの調査も行った。 これらの結果の概要についても、この報告書で報告する。
2 調査地域と植生
調査地は、屋久島西部の瀬切川上流の地域である(図1)。 標高は750mから1350mにあたる。調査地域の面積は7.5km2である。
調査域は全体として照葉樹林・ヤクスギ林移行帯にあたる。イスノキ、ウラジロガシなどの照葉樹が、スギ・モミ・ツガなどの針葉樹と混交している。 林床はハイノキを中心とする常緑低木に広く覆われている。調査域は大川林道を中心として植生の撹乱があり、伐採跡はスギの幼樹、ヒサカキやハイノキなどの低木に覆われている。 伐採後の年数は大川林道の終点近くや12支線沿線では22-33年程度で、西(入口方向)へ行くほど年数が経っている。もっとも西側の5班の調査域付近ではおよそ40年程度である。 伐採後の更新方法には2種類ある。一つは伐採後自然の更新過程に任せる天然更新で、12支線沿いや大川林道の終点付近からおおむね3班のdf付近まではこの方法で行なわれている。 もう一つは伐採後スギを植林する人工更新で、それより西の、3cdおよび4班と5班の調査域付近で行なわれている。
3 調査方法
参加者は、53人(前期27人、後期30人)である(表2)。1班と2班の調査は前期に、4班と5班の調査は後期に行った。
サルの集団の調査
調査は、ブロック分割定点調査法によった。調査域を500m四方のメッシュに分け、メッシュ内の音声の聞き取りやすい場所を定点に選び、そこに一人の定点調査者を配置した。 4-8つのメッシュをあわせて一つの班を作り、1-3人の統括者を配置した。定点調査者が群れを発見した場合はトランシーバーで統括者に連絡し、統括者が群れを追跡した。 群れが林道などを横切る時に群れを数え、群れの構成を調査した。
ヒトリザル
ヒトリザルを目撃した場合、時刻、定点からの距離と方角を記録した。距離は10m単位で記録し、20mより近いときには、1m単位で記録した。 最初に目撃したときから5分経過した場合、ないしはヒトリザルが10m以上動いた場合、時刻と定点からの距離、方角を再度記録し、ヒトリザルを見失ったときは、最後に確認した時刻と、定点からの距離、方角を記録した。
シカ調査
サル調査の定点近傍に、4メートル×50メートルのシカ調査区画を設け、その中に落ちているシカの糞塊数を、それぞれの糞塊の中の糞粒数とともに記録した。 この調査は、統括者が交替で行った。また、糞はすべて採取して、ここの糞粒の長径と短径を計測した。
ヒル調査
定点調査者は、1日に2回、体についているヒルの数をチェックした。一回目は、定点に到着した直後で、2回目は、1回目のヒルチェックから3時間後に行った。 ヒルがついていることを確認したら、チェックシートに時間、天候、数、体の部位、吸血の有無を記録した。観察者の皮膚に触れる前に発見したヒルは、DNA分析のために採取してアルコールに保存した。 観察者の皮膚に触れたヒルは、生きたまま採取して飼育し、その後の生存を記録した。
果実量調査
これまでの調査で、5メートル四方の植生調査区を、伐採地に26個、自然林に10個設置してある。 その中に生えている、サルが食べる液果をつける樹木(ヒサカキ、ハイノキ、オニクロキ、サカキ)について、結実数を直接計数した。伐採地の調査は前期調査中、自然林の調査は9月3日に行った。
植生調査
上記果実量調査を行っている伐採地の植生調査区の植生を調査した。4 集団密度の推定方法
集団密度はHanya et al. (2003b)にしたがって推定した。
まず、語句を以下のように定義する。｢集団｣とは、まとまって一緒に遊動しているサルの集まりのことをいい、具体的には、｢集団とは、500m以上近接して一緒に遊動しているサルの集まり｣であると定義する。 一方、｢群れ｣とは、｢同じ『集団』内で遊動する可能性のあるサルの社会的なまとまり｣であるとする。つまり、一つの群れは、一つの集団になってまとまっていることもあり、ばらばらになって二つ以上の集団を作っていることもある。 また、｢時間帯｣とは、｢6時0分から6時59分、10時0分から10時59分までのように、0分から59分までの1時間のこと｣をいう。
資料の分析に際しては、複数頭の音声か、オトナメスまたはコドモの目視情報のみを集団の発見とみなし、単発の音声、 具体的には前後1時間の間なにも情報のない音声情報や単独のオトナ･ワカモノオスの目撃情報は分析から省いた。
まず、各時間帯の各々の定点での集団の発見数を数えた。その時間帯の間に、その定点調査者が、500m以上離れた2か所からそれぞれ複数頭の音声を聞いた場合は2集団を発見したと数えた。 それ以外の場合については1集団を発見したと数えた。調査時間が30分未満の時間帯の資料は分析から除いた。次に、6時間以上調査した定点について、その日の各時間帯の発見数の平均を計算した(図3-1〜図3-2)。 最後に、全調査期間についてこの値の平均を求めた(図3-3)。
集団の発見数をn、集団密度をDとする。このふたつには---(1)
という関係がある。 λは定点からの距離yの地点にいる集団の発見率g(y)をというモデル(Half-normal model, Buckland, 1993)で近似したときの値である。この値は統括者による集団追跡と定点調査の資料を比較することでもとめることができる。今回の報告書にはこれらのパラメータを求めるための分析は間に合わなかったため、2000-2003年度の調査で求めた値(Hanya et al., 2005)を使用した。調査域を自然林(1a1b1c1d1f2a2b3a3b3c3d4a4b)、天然更新地1(若い天然更新地; 1e2f2g2h)、天然更新地2(古い天然更新地; 2c2d2e3f)、人工更新地1(若い人工更新地; 3d3e)、人工更新地2(古い天然更新地; 4c4d4e4f5a5b5c5d)の5つに分け、それぞれ個別に発見率を推定した。
5 調査日程
8/10 午前中、一部調査員は西部林道でヤクシマザル観察実習。前期集合、安房の屋久島環境文化研修センターで講習会
8/11 前期調査(日帰り)、一部調査員は西部林道でヤクシマザル観察実習
8/12-13 前期調査(日帰り)
8/14,15 台風のため停滞
8/16 入山、キャンプ設営
8/17,18 前期調査
8/19 悪天候のため下山
8/20 前期調査(日帰り)、定点調査は行わず、OM群の追跡のみ。打ち上げ、全体ミーティング
8/21 前期解散
8/22 後期集合。安房の屋久島環境文化研修センターで講習会
8/23 西部林道でヤクシマザル観察実習。午後、一部調査員でシカ調査(日帰り)
8/24-26 後期調査(日帰り)
8/27 後期調査(日帰り)。シカ調査、植物調査を行い、サル調査は休み。
8/28,29 後期調査(日帰り)
8/30 休養日
8/31 午後のみ、後期調査(日帰り)、悪天候のため定点調査は行わず、YY群の追跡のみ。
9/1 後期調査(日帰り)、打ち上げ、全体ミーティング
9/2 解散
6 調査結果
各々の定点の集団発見数の平均を図3に、統括者も含めたそれぞれの日の調査記録を図4に、長時間追跡できた集団の遊動図を図5に、それらをもとにして識別できた群れの遊動図を図6に示した。 また、群れのカウント例の一覧と、それをもとにして分かった、重複を除いた最低の群れサイズ、性年齢構成を表1に、識別できた個体の記録を図7に挙げた。
以下に、調査域全体の集団密度の傾向、シカ・ヒトリザル・ヒルの分布、果実生産の傾向、班ごとの植生、確認した群れの遊動域と遊動パターン、サイズと構成などの情報について述べる。
なお、それぞれの日に確認された集団については、1-9-aのように番号をつけた。最初の数字1は班名を、真中の数字9は日付(8月9日は9、8月11日は11、以下同様)を意味し、 最後のアルファベットはその班のその日についての通し記号である。また、調査域内ではHR群、PE群、OM群、SS群､YY群の5つの群れが識別されている。 それぞれの日に確認された群れがこれらの群れのうちのいずれかであることが確かめられた場合、1-11-b(HR)のように表記した。 ひとつの集団が同じ日に二つ以上の班にまたがって出現した場合、1-11-a(2-11-b)のように二つ以上の名前がつくことがある。
集団密度
調査域全体の集団発見数の95%信頼区間は0.10±0.03集団/定点/時間、最大は5cの0.24集団/時間、最低は1f,3c,3dの0集団/時間であった。 (1)式に代入して集団密度を計算すると、95%信頼区間が0.77±0.22集団/km2、最大が2.5集団/km2、最低が0集団/km2だった。 集団密度の植生による違いは、すべての組み合わせで有意でなかった (Tukey-Krammer法)。
図A 集団発見数の植生による変異。y軸の単位は集団数/定点/時間。平均+標準偏差を示す。
図B 集団密度の植生による変異。y軸の単位は集団数/km2。平均+標準偏差を示す。
シカ糞塊数
トランセクトあたり(50m*4m)の平均糞塊数の95％信頼区間は0.8±0.4である。最大値は4e(人工更新)の5、最小値は0であった。調査地では2010年4月から2011年3月にかけて、シカの捕獲が行われていた。 糞塊数は2009年は3.1±1.34、2010年は1.2±0.58、2011年は1.1±0.52、2012年は2.5±1.1、2013年は4.50±1.81、2014年は1.60±0.78、2015年は1.00±0.47、2016年は0.6±0.46、2017年は1.2±2.57、2018年は0.7±1.05だった。 糞塊数の植生による違いは、すべての組み合わせで有意でなかった(Tukey-Krammer法)。
図C シカ糞塊密度の植生による変異。y軸の単位は個/50m*4m。平均+標準偏差を示す。
ヒトリザルの発見頻度
単独行動のオトナオスまたはワカモノオスを目撃したヒトリザル発見頻度の95％信頼区間は(2.53±1.3)×10-3回/10分、最大は1a(一次林)の12.7×10-3回/10分、最小は0回/10分である。 図Dに植生ごとの発見頻度の違いを示した。すべての組み合わせで有意差は見られなかった(Tukey-Krammer法)。今回は調査中に合計17回発見した。 2003年から2018年までのヒトリザルの発見回数は、平均で各年19.1回(最小5回、最大30回)だったので、今年のヒトリザルの発見数は例年より下回った。
図D ヒトリザルの発見頻度。y軸の単位は回/10分。平均+標準偏差を示す。
ヒル調査
一回目のチェックでの総発見数は38匹で、一回のチェックで発見したヒルの平均数は0.27匹であった。二回目のチェックでの総発見数は7匹で、平均0.05匹であった。 一回目のチェックでは定点までの道のりで体についたヒルを取り除くことになるため、二回目よりも総じて発見数は多くなっていた。
DNA検査用に、71匹のヒルを採取し、これとは別に、飼育用に109匹のヒルを採取した。
果実量
5m四方の調査区内で結実していた液果の果肉重量の植生による変異を図Eに示した。果実生産量は天然更新1で最大で、天然更新2、自然林、天然更新3がそれに次ぐ。 人工更新地では果実生産が見られなかった。これは、過去の年の調査と同様の傾向である。すべての組み合わせで有意差は見られなかった(Tukey-Krammer法)。
図E 果実生産量の植生による違い。平均+標準偏差を示す。単位は果肉部分の湿重で、g/25m2を示す。
1班
1班の調査域のうち、12支線終点付近の1eには伐採地があるが、それ以外は自然植生である。伐採されたのは1995年ごろで、ヒサカキ、ハイノキの幼木が多い。 昨年の調査後に、12支線の終点付近と、そこから東300mくらいまでの範囲で除伐が行われ、これらの低木の多くが伐採された。自然林ではスギ、モミ、ツガ、ヤマグルマ、ハリギリなどが多い。
1班の調査域内では、HR群とPE群が確認され、そのほかにも群れの存在が推定された。
HR群はオトナを中心にほぼ全頭が個体識別されている群れである。12支線終点･1161mピーク、大川林道終点を結ぶ付近を西の端とし、1d付近を東の端として利用していることが確認された。 2019年4月に調査したときに構成を確認しており、それに今回確認したアカンボウを加えると、HR群の構成は、オトナオス5頭、オトナメス9頭、コドモ6頭、アカンボウ2頭の合計22頭だった。
PE群は、西は2e定点東の鞍部から、東は12支線終点付近、南は2h定点付近までの非常に狭い範囲でのみ発見された。 ただし、2h定点の南から声が聞こえることがあったので、おそらく調査域よりもさらに南に広がっていると考えられる。 ここ数年とは異なり、12支線の道上と道下の両方を利用していたため、定点調査を行った5日間の前期調査中、毎日出現が確認された。 確認された個体は、オトナオス1頭、オトナメス3頭、コドモ2頭、アカンボウ1頭の合計7頭だった。
ほかに、1b付近で5日間の調査中2日1日集団が出現した。1aでは集団の発見はなかった。
2班
2班の調査域は大別して、自然林、天然更新林、スギ人工林という3種の森林から構成されている。瀬切川右岸にあたる北側にはヤクスギ自然林が残され、左岸にあたる南側には伐採後の天然更新林とスギ人工林がモザイク状に見られる他、尾根沿いと南方の大川近くには伐採されなかった自然林が残っている。 自然林にはスギをはじめとした大径木が林立するが林床は比較的藪が濃く、地形の急峻さもあいまって集団の追跡が大変な場所である。 天然更新林は伐採後の年数によるが、ヒサカキ、オニクロキ、ハイノキ等の低木が密生する場所が多い。これらの低木は近年人の背丈を超える程度に成長しており、一時期のように全く歩けない程ではないものの、依然見通しの悪い状態になっている。 スギ人工林の大部分は継続的な間伐等の管理により見通しは良いが、一方で放置された間伐木により非常に歩き難い状態になっている。 これら3種の森林を、東西に走る大川林道本線と東に口をあけるU字型の12支線が、東西方向に3回横切っていることも調査域の特徴である。 集団の観察、カウント、追跡はこれらの林道に加え 林道間を結ぶ調査道を中心に行なわれた。
2班の調査域にはOM群、HR群、PE群の3つの識別群が現れた。OM群は、大川林側本線上では2c定点付近から3d定点付近で観察され、また12支線上では大曲周辺でOM群が観察された。 6日間の調査中5日OM群は確認されたが、人に対する警戒心が強く例年と同様にカウントが困難だった。今年度OM群が遊動した東端は2c定点の西側谷で、西端は3dの西側の谷であった。 一方で北端は不明であるが、少なくとも2cと3d定点間の北側を遊動しているようだ。南端は大曲周辺だが、3f定点でOM群と思わしき音声や目視があり、大曲よりも西でさらに南へ遊動している可能性がある。 構成はオトナオス3頭、オトナメス5頭、コドモ5頭、アカンボウ0頭の計13頭であった。昨年確認されていたアカンボウ3頭のうち、今年確認された1歳は1頭のみであった。
HR群は2班領域の左岸東側の2d定点の主に北側で確認された。また、PE群は2e定点東の鞍部より西側や南西を遊動し、2f定点や2h定点周辺で出現した。PE群と考えられる群れは調査5日中すべての日で2班領域に出現した。HR群とSS群の詳細は1班参照。
最後に、2g定点周辺でPE群とは異なる群れの音声が確認された。また、2a定点周辺でサルの群れが出現した。これらの群れについては統括者の人員を割くことが出来ず、本年も識別は行われていない。今年度2b定点ではサルの情報はなかった。
3班
3班の調査域は、自然林、天然更新林、スギ人工林という3種の植生帯に大別される。瀬切川右岸には伐採を経験していない自然林が広がり、スギ、ハイノキ、ヒサカキ、ヤマグルマなどが多く分布している。 この地域は大径木が多く、林床は比較的疎であるため見通しは良い場合が多い。大川林道と瀬切川に挟まれた部分には自然植生と人工林が、大川林道より南側には天然更新林と人工林がそれぞれモザイク状に入り混じっている。 人工林では見通しは良いが、最近まで間伐が行なわれていたため、林床に放置された間伐材のために歩き難い場所が散見される。天然更新林ではスギの幼樹やハイノキ・ヒサカキが密生して藪を形成している。 本年の調査では昨年までと同様瀬切川右岸で少なくとも1つの群れが、左岸および右岸でSS群が確認された。またOM群が3班域の東側で観察された。
SS群は大川林道上の3d定点付近から4d定点付近の間で観察された。3fへ向かう林道入口近くの道下に大きなカナクギノキがあり、そこでSS群がまとまって採食することがあった。 加えて、3ac4ab定点に向かう林道入り口の西側で樹上を渡たり、カウントされる時があった。SS群の遊動域の東端は3d定点付近であり、西端は4f定点付近である。 北端は不明だが瀬切川右岸を利用しており、3b定点の南や3c定点付近でSS群の音声が聞かれ、また、時系列的にSS群と思われる群れが3b定点で目視されることもあった。南端は3eよりも南側であるが今年度は確かめられなかった。 SS群は6日の調査中すべての日でカウントがおこなわれ、構成はオトナオス3頭、オトナメス7頭、コドモ5頭、アカンボウ1頭の計16頭であった。昨年確認されていたアカンボウ3頭のうち1頭のみ確認できた。 昨年度に比べて2頭オトナメスが減少(ワカモノメスを除く)したが、今年度は識別個体のオトナメスの生存が確認できず、また、オトナメスの最大頭数をカウントした3-29-fの際に、アカンボウ持ちのオトナメスが確認されていなかったため、死亡やカウント漏れの可能性が十分にある。
SS群以外では、右岸の3a,3b付近で不明群が出現している。OM群と考えられる群れは3班領域の左岸側の東で出現した。詳細は2班参照。
4班
4班の調査域には2種類の植生が認められる。右岸では江戸時代に行われた伐採以降は人為的撹乱を受けていないヤクスギの原生林が存在する。 その一方で左岸から分水嶺までは谷筋を中心として伐採後のスギ植林が見られ、尾根筋には原生林が少し残されている。また、2011年及び2015年に行われた間伐に伴い、無数の伐採道が大川林道から設置されている。
4班域の左岸ではSS群とYY群の2群が観察された。詳細な群れ構成については3班及び5班の所見を参照。8月26日、29日に4d定点周辺及び以東の大川林道上でSS群が出現した。 8月24日に4dおよび4fで出現した群れもSS群の可能性がある。4班領域にYY群が出現したのは、ほとんどが4班と5班の境界付近だったが、8月26日には、4cの定点付近でYY群が出現した。 この日は、5班ではYY群は出現しておらず、YY群は5班領域を主要な行動圏にしているものの、4班領域までまたがった行動圏を持っていることが確認できた。また、4e定点付近の森の中で、YY群と思われる音声が聞こえることもあった。
5班
調査地は全調査域の最西端にあたり、瀬切川左岸から分水嶺に位置する4区画である。植林後39-43 年経過した胸高直径30cmを越えるスギ林の林床、林道脇にイスノキ、ウラジロガシなどの照葉樹、ハイノキ、ヒサカキなどの落葉樹が、まだらな林床低木層を構成している。
5班ではYY群が確認できた。5班の調査域内には、7日間の調査期間中、6日間出現した。5班に出現していなかった8月26日には、4班に出現していた。5班での出現場所は4・5班境界付近から5c定点付近で、一昨年頻繁に出現していた万里の長城付近には現れなかった。 確認された構成は、オトナオス4頭、オトナメス5頭、コドモ2頭、アカンボウ1頭の、合計12頭だった。コドモの観察が極めて少なかったが、オスの小集団が出現した際に、コドモが一緒にいることがあり、群れ本体から離れていたのかもしれない。
謝辞
この調査を行うにあたって、屋久島森林生態系保全センターには、調査を許可していただきました。屋久島環境文化財団には、講習会場を提供していただきました。調査の準備段階で、多くの調査隊OB、OGの方にも御支援を頂きました。 今年の調査に先立つ2019年3月から5月に実施したクラウドファンディングでは、多数の方にご寄付をいただき、その資金で新しい車を購入することができました。支援者の方々のうち、お名前の掲載に同意してくださった方を、以下に明記いたします。 厚くお礼を申し上げます。
山本寛樹、向井真那、瀧本宏昭、金井塚務、野本弘子、早石周平、来栖任、上田陽子、鈴木洋平、安齋春加、松原始、日野貴文、保坂和彦、岸本真弓、冨田茉帆、滝澤信、上原昭宏、若尾美穂、加藤正吾、大坂桃子、本郷峻、中平楓、太田史生、岡村弘樹、桑島正充、佐藤みゆき、山嘉信、齊藤加奈子、東郷えりか、金岡雅浩、湯澤菜穂子、揚妻直樹、林響子、林幸平、牛崎学、久後貴寛、森嶋佳織、越川智瑛、村松佳子、ジロー、丸山晃央、河村正二、中林雅、角田史也、酒井治、館ソニン、林みな、田中徹、市川志野、阿部真也、松原幹、井上しのぶ、人見琢也、和田浩伺、中澤紗季、石部雄紀、服部正道、三上博民、神谷一彰、安藝真理、久世濃子、伊藤毅、白井啓、本田逸郎、米澤拓央、金井賢一、衣斐彼方、衣斐くるみ、大井徹、井野輝太、古閑丸斉良、丹羽百合子、清水(稲継)理恵、山岡尚平、黒江美紗子、山田容子、砥綿夕貴子、徳山豪、戸田恵美、北川訓子、足立隼、角田まみ、古田原真理、小出智矢、小貫志帆、蓑部節子、近藤健、旅するミシン店(順不同、敬称略)
引用文献
Buckland ST, Anderson DR, Burnham KP & Laake JL (1993) Distance Sampling: Estimating Abundance of Biological Populations. Chapman and Hall, London.
Hanya G, Noma N & Agetsuma N (2003a) Altitudinal and seasonal variations in the diet of Japanese macaques in Yakushima. Primates 44: 51-59.
Hanya G, Yoshihiro S, Zamma K, Kubo R & Takahata Y (2003b) New method to census primate groups: estimating group density of Japanese macaques by point census. American Journal of Primatology 60: 43-56.
Hanya G (2004a). Diet of a Japanese macaque troop in the coniferous forest of Yakushima. International Journal of Primatology 25: 55-71.
Hanya G (2004b). Seasonal variations in the activity budget of Japanese macaques in the coniferous forest of Yakushima: effects of food and temperature. American Journal of Primatology 63: 165-177.
Hanya G, Yoshihiro S, Zamma K, Matsubara H, Ohtake M, Kubo R, Noma N, Agetsuma N & Takahata Y (2004). Environmental determinants of the altitudinal variations in relative group densities of the Japanese macaques on Yakushima. Ecological Research 19: 485-493.
Hanya G, Zamma K, Hayaishi S, Yoshihiro S, Tsuriya Y, Sugaya S, Kanaoka MM, Hayakawa S & Takahata Y (2005) Comparisons of food availability and density of Japanese macaques in primary, naturally regenerated and plantation forests. American Journal of Primatology 66: 245-262.
丸橋珠樹,山極寿一,古市剛史 1986. 屋久島の野生ニホンザル. 東海大学出版会.
Yamagiwa J and Hill DA (1998). Intraspecific variation in the social organization of Japanese macaques: past and present scope of field studies in natural habitants. Primates, 39: 257-273.
Yoshihiro S, Ohtake M, Matsubara H, Zamma K, Hanya G, Tanimura Y, Kubota H, Kubo R, Arakane T, Hirata T, Furukawa M, Sato A and Takahata Y (1999). Vertical distribution of wild Yakushima macaques (Macaca fuscata yakui) in the western area of Yakushima island, Japan: preliminary report. Primates 40: 409-415.
