Silviase/jgov_v1.4_full_ann huggingface dataset:
した昆虫は大学に持ち帰り、目レベルまで同定し個体数を記録した。

## （7） 共分散構造モデルによる鳥類種数・個体数、および昆虫個体数変動の要因に関する検討

本研究の重点調査区のように、多分野にまたがる非常に多くの項目のデータがある場合、一般的な線形モデルのように、一つまたは複数の説明変数により、一つの目的変数の振る舞いを説明する、という方法ではモデル化に限界がある。あるモデルでは目的変数となった変数が、他の変数の振る舞いを説明するのに必要な説明変数になるようなことは往々にしてある。むろん、多くの変数がより錯綜した関係におかれていることも十分にあり得る。

そこで本研究では、共分散構造分析により変数間の関係をモデル化することを試みた。分析にはSPSS社のAmos 16.0を使用した。具体的な手順は以下の通りである。

## 1 ）使用したデータ (表(2)-2)

### a. NDVI

噴火前（1994-1998）、噴火後2年以内、噴火後3～5年、噴火後6～9年の最大NDVI値（前述）。各調査区の中心から、25m以内の平均値。

### b. 固定調査区の植生と樹種構成の把握

サブテーマ1で得た植物種ごとの被度と、階層別の植被率および高さ、樹種別のDBHを利用した。被度は0、+、1～5の7段階で記録し、各被度階級における被度の中央値、すなわち0、0.005、0.055、0.175、0.375、0.625、0.875に換算した上で、DCAによる序列化を行った。第1軸に、ほとんど回復していない植生から噴火の影響がほとんど認められない植生までの種組成の変化の軸が抽出された。これと、第2軸、第3軸のサンプルスコアを、以後の解析に利用した。

毎木調査で記録されたDBHからは、これを直径とする円の面積を調査対象樹木の胸高断面積と見なして算出した。同一調査区内で、同一樹種の胸高断面積を合計し、その種の胸高断面積合計値として解析に用いた。また、常緑樹種の胸高断面積合計が全体に占める割合を、植生遷移の進行状況の指標と見なして解析に用いた。

これ以外に、植生高、上層植被率（高木層植被率と亜高木層植被率の合計）、下層植被率（低木層植被率と草本層植被率の合計）、植物種数（植物社会学的群落調査で記録された植物の種数）を解析に用いた。

### c. 土壤呼吸と枯死木量の把握

枯死木は2000年噴火の影響を受けて大量に発生し、それらが他の生物の貴重な餌資源や生育地となっている可能性が高い。さらに、大量に存在している枯死木は防災上の観点からもその動態が注目されている。サブテーマ1の調査において、まだ直立している枯死木を「立枯木」、既に倒れてしまったものを「倒木」として、重点調査区でその量を計測した結果を使用した。

土壌呼吸は、土壌から放出されるCO₂量として評価され、土壌微生物等の有機物分解に伴う呼吸量と植物の根の呼吸量の総和である。したがって、土壌有機物分解の指標にもなりうるし、また植物の量やその活性の指標にもなりうる重要な炭素の流れである。土壌呼吸量、その温度依存性、温度0度の元で推定される土壌呼吸量、および土壌呼吸に関連する平均地温の各項目について、重点調査区で調査された結果を利用した。

### d. 固定調査区における土壌断面調査

重点調査区で記録された、O層厚さ、最上位A層厚さ、2000年噴火による火山灰の厚さの3項目
を解析に用いた。

### e. 土壤動物調查

重点調査区で得られた土壌動物のデータを利用した。ツルグレン装置を用いて得られたデータ、ハンドソーティングによるミミズ個体数、ツルグレン装置を用いてとりだした土壌動物総重量を解析に用いた。ツルグレン装置を用いて採取された土壌動物の種類が多かったため、DCAを行い、出現傾向より分類群を五グループに分けた。個体数が少なかった1つを除く4つのグループそれぞれに属する分類群の個体数の合計を解析に用いた。

グループ1：カニムシ、コムシ、ササラダニ、トゲダニ、甲虫幼虫、トビムシ、ケダニ、ミミズ綱

グループ2：ヤスデ、ハエ幼虫、アリ、エダヒゲムシ、カマアシムシ、ウズムシ

グループ3：ダニspp.、クモ、ムカデ、甲虫、ワムシ

グループ4：コムカデ、線虫、カメムシ、巻き貝

### f. 分解速度の調査

重点調査区に10種類（植物5種×メッシュサイズ2通り）のリターバッグをセットし、1年後に回収した結果から算出された分解定数を利用した。10種類の分解定数の間の相関は非常に高かったため、主成分分析を適用して第1主成分得点を求め、これを分解速度の指標とした。

### g. 鳥類調査

前述の定点センサスの結果を用いた。得られた結果から、越冬期、繁殖期それぞれにおける、各地点での種数と個体数を求めた。さらに越冬期、繁殖期それぞれのデータをDCAにより序列化して、調査区のスコアを求めた。ほとんど回復していない鳥類相から噴火の影響がほとんど認められない鳥類相までの種組成の変化の軸と判断された第1軸のスコアのみを、以降の分析で使用した。

### h. 昆虫調查

鳥類の食物になり得る動物の生息量の把握という観点から、重点調査区において、地上設置式フライトインターセプトトラップ（FIT）、吊り下げ式FIT、ベイトトラップ、スイーピング、ビーティングの5つの方法による調査で捕獲された節足動物のデータを利用した。地点×分類群のデータにまとめ、DCAにより序列化した結果から、出現傾向に基づき、以下の3グループに分けた。

グループA：ヨコエビ、バッタ

グループB：アリ、ハエ、カメムシ

グループC：甲虫、ハチ、クモ

それぞれのグループに属する分類群の個体数の合計を以降の分析に用いた。

### i. 火山ガス濃度の調査

固定調査区にデータロガー付きの二酸化硫黄ガスセンサーを設置した。結果は、二酸化硫黄ガス濃度が、0.1ppm以上、0.2ppm以上、0.3ppm以上、0.4ppm以上、0.5ppm以上になった回数を1ヶ月あたりに換算して算出し、各階級の中央値（0.5ppm以上の場合は0.6ppm）と1ヶ月あたり回数の積により表現した。

## 2 ）相関分析と主成分分析
表(2)-2 分析に用いた変数一覧

調査対象	設定した変数(上段:変数名、下段:説明)
植生	LowCover	VegHeigh	HiCover	BAtotal	VegRich	BA%EVG	ps_DCA1	ps_DCA2	ps_DCA3
下層植被率	植生高	上層植被率	胸高断面積合計	植物種数	常緑樹胸高断面積比率	DCA第1軸スコア	DCA第2軸スコア	DCA第3軸スコア
土壤	O_layer	A-layer	VolcAsh
O層厚さ	A層厚さ	2000年噴火火山灰厚さ
土壤呼吸	Dojokoky	Toboku	Tatigare	Ondoizon	Odeg_kok	STavg
土壤呼吸量	倒木量	立ち枯れ木の量	土壤呼吸の温度依存性	0度時の土壤呼吸量	平均地温
腐植の分解定数	bunkaiPC	分解定数の第1主成分得点
土壤動物	DojoG5	mimizu	DojoG2	DojoWeig	DojoG4	DojoG3
コムカデ、線虫、カメムシ、巻き貝	ハンドソーティング時のミミズ個体数	カムシ、ササラダニ、トゲダニ、甲虫幼虫、トビムシ、ケダニ、ミミズ網	土壤動物総重量	ダニ-pp、クモ、ムカデ、甲虫、ワムシ	ヤスデ、ハエ幼虫、アリ、エダビゲムシ、カマアシムシ、ウズムシ
昆虫	yok_batt	arhaekam	ko_ha_ku	甲虫・ハチ・クモ
ヨコエビ・バッツ	アリ・ハエ・カメムシ	Abund_BR	SR_BR	bdbrDCA1	bdwiDCA1
鳥	Abund_WN	SR_WN	繁殖期個体数	繁殖期種数	繁殖期DCA1	越冬期DCA1	火山ガス	SO2
越冬期個体数	越冬期種数	繁殖期個体数	繁殖期種数	繁殖期DCA1	越冬期DCA1	火山ガス	SO2
GIS	DEM2001	Max_NDVI_2	Max_NDVI_3	Max_NDVI_4	DEM1983	Max_NDVI_1	火山ガス	SO2
2001年標高	噴火直後NDVI	噴火数年後NDVI	直近NDVI	1983年標高	噴火前NDVI	亜硫酸ガス濃度推定值

以上の、合計41変数の全ての組み合わせについて、ピアソンの積率相関係数を算出した。さらに、相関行列に基づく主成分分析を実行し、類似性の高い変数のグループを探索した。

## 3 ）共分散構造分析

本研究では、多くの変数が調査された一方で、調査地点数は限られている。

このような状態でも共分散構造分析を行うことはできるが、適合度の推定に問題を生じ、モデルの妥当性の検証が難しくなってしまう。

そこで本研究では、2)の相関分析と主成分分析により関係が特に強い変数を抽出した上で、冗長な変数を除き、残りの変数の一部を対象としてモデル構築を行うことを、変数の構成を変えて繰り返し主要な変数間の関係性を検討する、という作業を行って、変数間の関係の検討を試みた。

## （8） 三宅島におけるスダジイ堅果の地域間変異

三宅島内の3地点間でスダジイ堅果の大きさと形状を比較した。調査地点は、坪田林道の中間の林道上約100メートル（以下坪田）、薬師堂と御祭神社社殿に続く参道上約200メートル（以下薬師堂）、アカコッコ館駐車場から同館に至る道路上約150メートル（以下アカコッコ館）である。

## （9） 三宅島のスダジイ堅果における昆虫による食害率の地域間変異

三宅島のスダジイの堅果の大きさと形状が場所により異なる場合、大きさと、昆虫による食害率の間に関係があることが予想される。この仮説に基づき、（8）と同じ3地点においてスダジイ Castanopsis sieboldiiの堅果を採取し、昆虫による食害率を調べた。

## （１０）三宅島におけるオオバエゴノキの種子の昆虫による食害率

伊豆諸島三宅島には葉や花の大きなオオバエゴノキStyrax japonica var. kotoensisが自生しているが、その種子を食害する昆虫種についてはあまり調べられていない。そこで本研究では三宅島のオオバエゴノキ（以下、エゴノキ）の落下種子を採取し、昆虫による食害率を調べた。

## （１１）大路池周辺に植栽されたツバキ園芸種の花形態の観察：野生ヤブツバキ集団との交雑の可能性

ヤブツバキ Camellia japonicaは、日本・中国などの東アジアに分布するツバキ科の常緑広葉樹
で、伊豆諸島にも広く分布する。大きく色鮮やかな花弁をもつことから、本種をもととして数千種に及ぶ園芸品種が観賞用に作り出されてきた。また本種は火山ガスに対する耐性もあり、三宅島島内では2000年火山噴火による植生破壊を強く受けた高標高地の疎林から、低標高地の成熟した常緑樹林内にいたるまで広く生育している(Abe and Hasegawa 2007) $ ^{5} $ 。いつぽう、島の南部に位置する大路池の周辺では、10年以上前から園芸品種のツバキが多数植栽されており、自生するヤブツバキとの交雑が懸念されている。

2013年2月に、大路池北西部の林道に植栽されているツバキ園芸品種の花の形態を観察、写真撮影した。なお、園芸品種のなかには葯や雌しべの退化した八重咲きのものがあるが、これらは花粉の受け渡しができないため野生集団との交雑のおそれはないと考えられる。それに対して、葯から花粉を露出させる株や雌しべを残している株は、野生集団の遺伝子汚染を引き起こす可能性が否定できない。このため、調査では個体ごとに、葯から花粉を露出させる株や雌しべを残しているかどうかに焦点をあてて観察した。

## （１２）キイチゴ属２種の花の咲く向きと訪花昆虫

三宅島にはカジイチゴRubus trifidusとハチジョウイチゴR. ribisoideusが多く見られる。カジイチゴの花は上を向く一方、ハチジョウイチゴの花は下を向く（七島花の会 2011） $ ^{6)} $ 。2012年3月から4月にこのキイチゴ属2種についてフェノロジーおよび訪花昆虫を調査した。

## （１３）植物-訪花昆虫間の送粉ネットワークの解明

植生遷移段階が異なる島内6地点に、長さ200-300mのライントランセクトを設置した。調査は、2012年4月・6月・7月・8月・10月の計5回行った。各回の調査にあたっては、各トランセクトを日中約6時間歩いて往復し、訪花している昆虫(体長5mm以上)を見つけしだい捕虫網で捕獲した。またトランセクトにおける開花植物種を記録した。

## （１４）三宅島における大型鳥類の種子散布者としての機能評価

宇都宮大学農学部の杉田昭栄教授の研究室と共同で、ハシブトガラス飼育個体（3個体）を用いて、種子排泄実験を行った。缶詰のマンゴーなどの果肉に、ハゼノキ・カラスザンショウなどの種子を挟み込んだ擬似果実を作成し、採食させた。採食終了から種子が排出されるまでの時間を全ての種子について記録した。

## （１５）島内で分布を拡大しているメダケが鳥類の生息に与える影響の検討

ラインセンサスデータを利用し、メダケが優占していた区間とそうでない区間の間で鳥類種数が異なるかどうかを、一般線型モデルにより検討した。100mあたり鳥類種数を目的変数、樹木植被率対数値と、メダケ優占群落が区間の2割以上あれば1，そうでなければ0となるダミー変数を説明変数とするモデルを構築した。

## 4 ．結果及び考察

## （1） 現存植生図

空中写真（1999年4月、2001年2月、2006年3月、2010年7月）に基づきそれぞれの時期の現存植
図(2)-5 固定調査区(重点調査区およびそれ以外の調査区)の近傍25m内におけるNDVI値の経時変化。衛星画像から算出したNDVI値の経時変化に基づく。調査区は、①噴火や火山ガスによる被害を受けなかった（青）、②噴火の被害は軽微だがその後火山ガスによる被害があり、現在は回復過程にある（緑）、③噴火直後に植生が大きく損なわれたが、その後短期間で回復過程に入った（赤）、④噴火により植生が損なわれた上さらに火山ガスによる被害も生じたが現在は回復過程にある（橙）、⑤噴火直後に植生が大きく損なわれ、比較的最近まで回復が見られなかった（紫）、⑥噴火による被害が著しく、その後の回復も見られない（黒）、の6つに分けられた。重点調査区は、①②③が各2、④が3、⑤⑥が各1となっており、全てのパターンをカバーしている。

生図を作成した上で、現地踏査を通じて細部の修正を行った。原図は大判になるため、本報告書においては報告書の判型に縮小したものを巻末に添付する（資料1）。原図については、要請があり次第、提供可能である。

## (2) NDVI

2005年4月、2009年4月、2010年4月撮影のIKONOS画像、2011年12月撮影のQuickbird画像、2013年4月撮影のGeoEYE(IKONOS後継機)画像からそれぞれ計算されたNDVIを地図上に表示した結果を、順に巻末の資料2として示す。こちらについても、原図については要請に応じて提供可能である。

衛星画像については、2011年12月の画像は雲やガスがほとんどなく、単一の画像により当該時期の全島の状況を把握できると考えられた。それ以外の画像については、雲やガスの影響があり、またセンサーの仕様に関する情報が販売元から十分に提供されなかったこともあって、画像間の直接の比較は難しいと考えられた。そこで、別途入手した中解像度の画像4時点分×4セットから算出されたNDVI値と、2011年12月の値を比較した。それぞれの固定調査区（11の重点調査区を含

表(2)-3 検討した統計量の一覧

種別	指標	单位	英名	和名
Amplitude	Sa	nm	Roughness Average	平均粗さ
Sq	nm	Root Mean Square	二乗平均平方根粗さ
Sy	nm	Peak-Peak Height	最大高さ
Sz	nm	Ten Point Height	十点平均高さ
Ssk	Surface Skewness	表面歪度
Sku	Surface Kurtosis	表面尖度
Hybrid	Sdq	%	Root Mean Square Slope	二乗平均平方根勾配
Sdr	Surface Area Ratio	表面面積比
Functional	Sbi	Surface Bearing Index	表面ペアリングインデックス
Sci	Core Fluid Retention Index	コア溶液保持インデックス
Svi	Valley Fluid Retention Index	ハレー溶液保持インデックス
Spk	nm	Reduced Summit Height	縮小されたサミット高さ
Sk	nm	Core Roughness Depth	コア粗さ深さ
Svk	nm	Reduced Valley Depth	縮小されたパレー深さ
Sdc0_5	nm	I-h% Height Intervals of Bearing Curve	ベアリング曲線のI-h%高さ間隔 0_5
Sdc5_10	nm	I-h% Height Intervals of Bearing Curve	ベアリング曲線のI-h%高さ間隔 5_10
Sdc10_50	nm	I-h% Height Intervals of Bearing Curve	ベアリング曲線のI-h%高さ間隔 10_50
Spatial	Sdc50_95	nm	I-h% Height Intervals of Bearing Curve	ベアリング曲線のI-h%高さ間隔 50_95
Sds	1/μm²	Density of Summits	サミット密度
Std	度	Texture Direction	テクスチュアの方向
Stdi	Texture Direction Index	テクスチュア方向のインデックス
Srw	nm	Radial Wavelength	放射状の波長
Srwi	Radial Wave Index	放射状波インデックス
Shw	nm	Mean Half Wavelength	平均半波長
Sfd	Fractal Dimension	フラクタル次元
Scl20	nm	Correlation Length Parameters	相関長パラメータ 20%
Scl37	nm	Correlation Length Parameters	相関長パラメータ 37%
Str20	Texture Aspect Ratio Parameters	テクスチュア・アスペクト比パラメータ 20%
Str37	Texture Aspect Ratio Parameters	テクスチュア・アスペクト比パラメータ 37%

図(2)-6 Core Fluid Retention Index と噴火被害段階の関係（Eが最も甚大）
左が繁殖期、右が越冬期の結果を示す。種の総計には利用場所がその他のものや確認できなかった記録を含む。

む23調査区、鳥類の定点センサスで用いた調査区と同じ）の25m以内におけるNDVI値の経時変化を、以上のデータから算出した。その結果、調査対象とした調査区は6つの類型に分けられた。11の重点調査区は全類型をカバーしていることが確認された（図(2)-5）。

さらに、表(2)-3に挙げた統計量metricsについて、その性質を検討した。量的指標（NDVIなど）は、植生の単位面積当たりの量的変化を効果的に捉えることができ、噴火など自然的あるいは人為的インパクトによる植生への影響と回復の状況を類型化できるため、広く用いられている。ここで取り上げた表面形状指標（表面粗さ）を用いた場合、植生被害は、表面形状の振幅の乱れとして評価することができる。例えば、Fluid Retention特性の機能パラメータとして評価することができ、被害を受けて不揃いになった植生では値が大きくなる傾向が見られた(図(2)-6)。量的指標に比較して、植生変化への感度特性が高い傾向が認められ、両者を組み合わせることで、植被の多少だけでなく植被の質的な面も含めた評価指標として利用できるようになると考えられた。

## （３）鳥類の定点センサス

固定調査区では、繁殖期には18種、越冬期には20種（ともに不明種を除く）が記録された。噴火の被害が大きな場所ほど種が少なくなるという傾向は、噴火直後から一貫して同様であった。

個々の個体の利用空間については、分布が限られている種は常緑樹を好んで利用する傾向が認められた（表(2)-4）。繁殖期の方が越冬期よりもこの傾向がより明瞭であった。利用していた場所の高さについても、種ごとの特徴は認められたが、分布が限られている種とそうでない種の間の違いは認められなかった。

種組成を規定する要因を明らかにするため、サブグループ1により調査された固定調査区の植生に関する情報、または共通作業によって得られたNDVIを説明変数として、CCAを行った。その結果、どちらの場合でもよく似た序列が得られた（図(2)-7）。上層植被率で代表される高木・亜高木による植被の発達具合と、植物種数に代表される植生全体としての発達の程度が、相互にある程度独立して鳥類相を規定しているという結果になった。NDVIを用いた場合、噴火直後のNDVIは上層植被率とほぼ同じ意味合いを持ち、2012年のNDVIは植物種数と同様の意味合いを持つという結果になった。直近のNDVIは植被の全体としての発達を指標するのに対し、噴火直後のNDVIは噴火を生き延びた高木（特に常緑樹）の量を反映しており、これは2012年の三宅島においてはその時点で生育していた高木（特に常緑樹）の植被量の指標にもなっていたことが示唆された。

## 表(2)-4 定点センサスで記録された種が初認時に利用していた場所。
図(2)-7 定点センサスの結果得られた繁殖期の種組成データに対する CCA の結果。第1軸（横軸）と第2軸（縦軸）で張られる平面を示す。左はNDVIをセカンドマトリクスとした場合。右は植生調査結果を利用した場合。

基本的には、種組成の主要な変化は第1軸に反映されている。左が種の豊富な定点、右が種の少ない定点。ホオジロ、ハシブトガラス、モズは種の少ない＝樹林が回復していない＝所でも(or所を好んで)出現するが、コマドリ、ミソサザイ、カラスバト、ヤマガラ、アカコッコ、コグラは樹林が回復しているところでのみ出現した。

## （４）鳥類のラインセンサス

ラインセンサスは、同じ手法で2002年以降調査を継続してきた（加藤・樋口2003 $ ^{7} $ 、2011 $ ^{8} $ ）。ここでは、2002、2006、2009年の繁殖期における調査結果、および越冬期である2002年2月、2004年12月、2008年2月、2011年2月の調査結果を、本研究の実施期間に行った2012年繁殖期および2013年2月の調査結果とあわせて解析した。

噴火後の三宅島において、鳥類の種数および個体数は、樹木植被率と強い相関を示すとされてきた（加藤・樋口 2003 $ ^{7} $ ）、2006 $ ^{9} $ ）、2011 $ ^{8} $ ）。2012年繁殖期および2013年2月においても、その傾向は変わらなかった。繁殖期では、2004～2007年に回帰直線の切片と傾きがともに増大し、その後減少したとされる（加藤・樋口 2011） $ ^{8} $ 。本調査期間でもその傾向は変わっていない（図(2)-8）。紙面の都合で図は省くが、越冬期でも同様であった。

図(2)-8 ラインセンサスで記録された樹木植被率と鳥類種数の関係
2012年繁殖期および2013年2月の調査結果を、TWINSPANにより整理したところ、種組成変化のパターンは、繁殖期、越冬期ともに、典型的な入れ子構造を呈した(図(2)-9、越冬期の結果は省略)。特に限られたところにしか分布していない種、すなわちヤマガラ、コマドリ、ミソサザイ、カラスバト、コゲラは、噴火による植生被害の影響を強く受けている種と判断できる。

樹木植被率および全階層植被率(高木層、亜高木層、低木層、草本層の植被率の合計値)、またはNDVIをセカンドマトリクスとし、過去の調査結果も含めてCCAを実行した結果が図(2)-10である(越冬期の結果は省略)。セカンドマトリクスの異なる2通りの分析の間で、結果

図(2)-9 ラインセンサスで記録された鳥類種組成の場所による違い。分類はTWINSPANに基づく、

はきわめてよく似ていた。繁殖期、越冬期いずれの場合も、第1軸は植被、特に樹木植被の発達具合とよく対応し、第2軸は調査年に対応した。

第2軸についての種のスコアから、繁殖期では、最近になるほど、ウグイス、ヒヨドリ、メジロ、ホオジロが多くなる傾向があると判断された。年を経るに従って森林内部性の種が減少し、草地性のホオジロ、ササ原や藪を好むウグイス、あるいはジェネラリストとしての性質が強いヒヨドリやメジロが優占度を増していることを示唆する。これは、島の植生の多くを占めるタイプ

図(2)-10 繁殖期のラインセンサス結果に対し、CCAを適用した結果。左は異なる時期のNDVIと調査年を、セカンドマトリクスとした場合、右は現場で記録した植被率と調査年をセカンドマトリクスとした場合。

year：調査年、gy2012：2013年4月NDVI、just_aft：噴火直後のNDVI、before：噴火前のNDVI、T：樹木植被率、ALLsum：全階層植被率
が、常緑広葉樹林から草地および落葉広葉樹林(オオバヤシャブシ群落)へと移行したことに起因すると思われる。噴火直後に直ちに鳥類が減少あるいは鳥類相が変化せず、徐々に変化が進んでいるのは、いわゆる混み合い効果が作用したことと、噴火直後からしばらくの間は枯死木や衰退木から発生した昆虫が森林性の鳥類の食物として機能していた殻であろうと考えられる(加藤・樋口 2011) $ ^{8} $ 。

ラインセンサスのデータに基づき、100m 調査区間あたりの出現個体数を予測するモデルの構築を試みた。鳥類のデータとして 2002、2005、2009、2012 年の繁殖期の調査結果を用い、過去のものも含め衛星画像により算出された NDVI と、共通作業により作成された現存植生図を ArcGIS 上に読み込み各区間から 50m 以内にあるスギ・ヒノキ植林の面積割合を算出した結果を説明変数として、重回帰モデルを構築したところ、以下の結果が得られた。

2002 年 5 月

個体数 = -1.5 + 15.9 NDVI(2000-2002 年最大値) - 0.013 (1999 年 50m 範囲内スギ・ヒノキ植林比率)

$$ \mathbf{R}^{2}=0.66,\ \mathbf{N}=37 $$

2005年5月

$$ 個体数 =-3.5+18.3NDVI(2000-2002 年最大值 ) $$

$$ \mathbf{R}^{2}=0.73,\mathbf{N}=35 $$

$$ 2009\ 年 \ 5\ 月 $$

$$ 個体数 =-4.2+18.3NDVI(2003-2005 年最大值 ) $$

$$ \mathbf{R}^{2}=0.61,\mathbf{N}=37 $$

$$ 2012\ 年 \ 6\ 月 $$

個体数 = 1.1 + 13.0 NDVI(2013 年 4 月) + 0.11 （2010 年 50m 範囲内スギ・ヒノキ植林比率）

$$ \mathbf{R}^{2}=0.70,\ \mathbf{N}=37 $$

衛星画像をこのような予測に使うことの妥当性を検討するため、現場にて記録された植被率を説明変数に加えたモデルを構築したところ、以下のようになった。

2002 年 5 月

$$ 個体数 =2.9+0.12\ 上層植被率 \% $$

$$ \mathbf{R}^{2}=0.73,\mathbf{N}=37 $$

$$ 2005\ 年 \ 5\ 月 $$

$$ 個体数 =-0.42+0.057\ 上層植被率 \%+0.08\ 下層植被率 \% $$

$$ \mathbf{R}^{2}=0.80,\mathbf{N}=35 $$

$$ 2009\ 年 \ 5\ 月 $$

个体数 = 1.8 + 5.1 NDVI(2003-2005 年最大值) + 0.09 上层植被率% + 0.03 下层植被率%

$$ \mathbf{R}^{2}=0.90,\mathbf{N}=37 $$

2012年6月

個体数 = -0.27 + 3.0 NDVI(2013 年 4 月) - 0.027 （2010 年 50m 範囲内スギ・ヒノキ植林比率）

$$ +0.14\ 上層植被率 \%+0.023\ 下層植被率 \%\\\ R^{2}=0.95,N=37 $$

ステップワイズ法による変数選択の際に、2002年および2005年についてはNDVIや植林比率は説明変数として選択されなかったが、モデルの当てはまりは現地調査によるデータのみを説明変数とするモデルのほうが優れていた。他の2時期のモデルも、現地調査のデータ、特に上層植被率の影響が強く現れていた。とはいえ、決定係数0.6～0.7のモデルが得られており、実用性は高いと判断された。

衛星画像と現存植生図より、上述の重回帰モデルを用いて繁殖期の鳥類の個体密度を島の全域について予測した結果が図(2)-11である。既に報告したように、三宅島では鳥類の種組成の変化
は入れ子形式であること、個体密度と種密度はほぼ同じ要因によって同様に（係数の値は異なるが）説明されていることから、個体数が多いところがすなわち種組成が豊富なところである。年を経るに従って山頂近くまで何らかの植被は回復しているのだが、鳥類の空白域はむしろ拡大していることがわかる。また、個体数が多いところでも、推定される個体密度は少なくなっている。

図(2)-11 衛星画像から算出したNDVIを用い、重回帰式から予測された島内の鳥類種密度。左上2002年, 右上2006年, 左下2009年、右下2012年。繁殖期についての予測。

NDVI から種数を予測するモデルについて同様の検討を行い、以下のモデルが得られた。

2002 年 5 月

$$ 種数 = 1.40 + 11.0 NDVI(2000-2002 年最大值) $$

$$ \mathbf{R}^{2}=0.44,\mathbf{N}=37 $$

2006年5月

種數 = -1.20 + 16.0 NDVI(2000-2002 年最大值)

$$ \mathbf{R}^{2}=0.66,\mathbf{N}=34 $$

2009年5月

種數 = 0.21 + 11.7 NDVI(2000-2002 年最大值)

$$ \mathbf{R}^{2}=0.62,\mathbf{N}=28 $$

2012年6月

$$ \text{種数} =0.883 + 6.437(2000 - 2002 \text{ 年最大值}) + 4.27 \text{ NDVI}(2013 \text{ 年 4 月}) \quad R^{2}=0.59, \text{ N}=38 $$

今度は植林地の項目は選択されず、決定係数はやや小さな値となった。これらのモデルに基づいて、それぞれの年における島全体の状況を予測した結果を図(2)-12に示した。

CCA の結果から、年を説明変数に入れてモデルを構築することができるようにも思える。調査年を増やした方が年変動をよりよく反映できると考え、2005 年と 2008 年のデータも加えてモデルの

図(2)-12 衛星画像から算出したNDVIを用い、重回帰式から予測された島内の鳥類個体密度。左から2002, 2005, 2009年繁殖期についての予測。
構築を試み、以下のモデルが得られた。

種數 = 274 - 0.137 調查年西暦 + 13.4 NDVI(2000-2002 年最大值) R²=0.60, N=192

調査年は説明変数に取り入れられたので、現在の傾向が変わらない限り、島の広範囲で鳥類群集の衰退が続くことが予想される。

なお、2012年の種組成に基づき、定点センサスの結果も含めて島内の鳥類相の場所による違いを図示したのが図(2)-13である。これと比較することで、いずれのシミュレーションにおいても、鳥類の分布パターンの大筋は適切に反映されていることがわかる。

一連のモデルにおいて、NDVI(2000-2002年最大値)が変数に選択されている。定点センサスの項でも述べたが、噴火直後のNDVIは噴火を生き延びた高木（特に常緑樹）の量を反映し、2012年の三宅島においてはその時点で生育していた高木（特に常緑樹）の植被量の指標にもなっていたためと考えられた。

図(2)-13 2012年6月の島内における鳥類相の場所による違い。ラインセンサスと定点センサスの結果に基づく。●：鳥類が記録されなかった、●：イイジマムシクイ、ウグイス、ホオジロが記録された。●：上の3種に加えヒヨドリやメジロが記録された。●：上の5種に加えて、アカコッコが出現。●：コマドリ、カラスバト、ミソサザイなど樹林性の希少種まで出現。

図(2)-14 大路池の鳥類相の経年変化。左が繁殖期、右が越冬期の結果を示す。縦軸は個体数。

図(2)-15 御祭神社の鳥類相の経年変化。左が繁殖期、右が越冬期の結果を示す。縦軸は個体数。
## （５）噴火直後からの長期的動向の検討

繁殖期の大路池では、2004年以降、小刻みな増減はあっても個体数は概ね安定しているように見える(図(2)-14)。但し、これはヒヨドリ、メジロ、ウグイス、ホオジロが増えているから。この4種を抜いてしまうと、2008年以降漸減傾向にあるようにも見える。越冬期ではこの漸減傾向がより明瞭に認められる（同）。

御祭神社および曽利川橋では、2013年の繁殖期、2013年および2014年の越冬期に、例年よりも多数の個体が記録された。これが一時的なものであるのかどうかが問題であるが、それより前を見るとやはり2007-2008年以降漸減傾向にあるように見られる（図(2)-15、曽利川橋での結果は省略）。

## （６）鳥類の食物となり得る昆虫の調査

## 1 ）地上設置式フライトインターセプトトラップ (FIT)

ハエ目（双翅目）が多く採集されたので、ハエ目とその他の分類群を分ける。まずハエ目の結果は図(2)-16のようである。季節別にみると、4月が最も多く、6月、9月の順に少なかった。被害程度については、4月は被害の大きな調査区において個体数が多く、6月は被害が中程度の調査区において個体数が多かった。最大の個体数は4月のIZ2であった。

次に、ハエ目以外の分類群については図(2)-17に示したとおりである。ハチ目、甲虫目、カメムシ目が多く採集された。局所的にアリ類が多く採集された地点があった。4月、6月に被害の中程度の調査区（N4-M1）において多くの分類群が採集された。ハチ目は4月、6月にはヒメバチ類が多く見られた一方、9月はツチバチ（キンケハラナガツチバチ）が多くみられた。カメムシ目はヨコバイが多く採集された。

図(2)-16 地上設置式FITによるハエ目昆虫の捕獲結果

図(2)-17 地上設置式FITによるハエ目以外の捕獲結果

図(2)-18 吊り下げ式FITによるハエ目昆虫の捕獲

図(2)-19 吊り下げ式FITによるハエ目以外の捕獲結果
## 2 ）吊り下げ式FIT

ハエ目（双翅目）が多く採集されたので、ハエ目とその他の分類群を分ける。まずハエ目の結果を図(2)-18に示す。4月に被害の比較的大きな調査区において多くの個体数が採集された。一方6月はどの調査区においてもあまり多く採集されなかった。

ハエ目以外の結果を図(2)-19に示す。4月はIG9、N1においてハチ目が多く採集された。6月はN4においてアリ類が多く採集された。ハチ目はヒメバチ科が多くを占めた。

## 3 ）ベイトトラップ

ベイトトラップの結果を図(2)-20に示す。4月、6月共に局所的にある分類群の節足動物が多く採集された。4月はIZ1で甲虫目（ゴミムシ科）、AK1でアリ類、IG9およびN1で甲虫目（ゴミムシ科）、IG7でハエ目（ガガンボ幼虫）が多く採集された。6月はIZ1でヨコエビ目、N1で甲虫目（ゴミムシ科）が多く採集された。4月に被害の中程度の調査区において個体数が多く、また比較的多くの分類群が採集された。

## 4 ）スイーピング

ハエ目（双翅目）が多く採集されたので、ハエ目とその他の分類群を分ける。まずハエ目の結果を図(2)-21に示す。7月よりも6月に多く採集された。6月はIZ2、IG9で特に個体数が多かった。被害の大きな調査区（M1-OY8）ではほとんど採集されなかった。次に、ハエ目以外の分類群の結果を図(2)-22に示す。6月よりも7月に多くの個体数、分類群が採集された。7月は被害が中程度（IZ2-IG8）の調査区において多くの個体数、分類群が採集された。AK1で6、7月に多く採集されたカメムシ目はヨコ

図(2)-20 ベイトトラップによる捕獲結果

図(2)-21 スイーピングによるハエ目昆虫捕獲結果

図(2)-22 スイーピングによるハエ目以外捕獲結果

図(2)-23 ビーティングによる捕獲結果
バイ科の成虫及び幼虫である。被害が中程度以上の調査区（IG9-IG7）においてアリが多く採集された。これはハチジョウススキから多く採集されたためである。このことはカメムシ目についてもいえた。

## 5 ）ビーティング

採集された節足動物の結果を図(2)-23に示す。4月、6月ともカメムシ目およびアリ類が多く採集された。カメムシ目はヨコバイ科が多くを占めたが、IG8、IG7ではハチジョウ

回復の程度

図(2)-24 節足動物群集の噴火後の回復パターン

ススキを吸汁するナガカメムシ科もみられた。4月、6月とも、被害が中程度の調査区（IG9-IG8）において多くの個体数、分類群がみられた。

## 6 ）考察

調査を通じてハエ目昆虫が最も多く採集された。次いでハチ目、カメムシ目が多くみられた。2000年噴火の被害の程度についてみると、調査区被害が少なかった調査区（IZ1やAK1）で多くの節足動物が採集されることはあまりなかった。一方、被害が中程度の調査区（N4-IG8）では個体数、分類群ともに多くの節足動物が採集された。被害が大きな調査区では4月にハエ目昆虫が多く採集された。以上のことから、いくつかのことが示唆される。

まず、三宅島の鳥類のエサ資源になっているのはハエ目昆虫であることが示唆される。ハエ目の成虫は4月から6月に多く出現したことから、鳥類にとってこの時期にこれらのハエ目昆虫を摂食し、繁殖することが都合が良い。より詳細にみるとハエ目昆虫は4月に被害が大きな調査区から出現する一方、6月に被害が中程度の調査区から多く出現した。したがって鳥類はハエ目昆虫をエサ資源として利用する場合、春は被害の大きな地域を訪れ、その後夏にかけては被害が中程度の地域を訪れるのがよいだろう。この点については鳥類の調査結果と組み合わせて解析することによって明らかになると考えられる。

次に、噴火の被害が中程度の調査区に多くの節足動物がみられたことから、生物多様性は生態系の回復の過程の中途において最も高まることが示唆される（図(2)-24）。このことはすでに他の研究でも指摘されており（たとえば宮下・野田 2003 $ ^{10} $ ））、今回の調査もこのことを支持する。

さらに本研究の結果から、三宅島の節足動物相の回復は従来考えられているよりも早いことが示唆される。すなわち三宅島において節足動物相は10年程度でかなり回復し、10年から20年にピークに達するように思われる。一方、ハエ目は節足動物相の回復とともに個体数、多様性が減少していく傾向にあると思われる。このことは鳥類や他の生物との相互作用に重要な影響をもつだろう。今後は鳥類や植生などの調査結果を組み合わせて解析を行い、これらのことについて検証することを予定している。

## （7） 共分散構造モデルによる鳥類種数・個体数、および昆虫個体数変動の要因に関する検討

主成分分析を適用したところ、多くの変数が第1主成分と強い関係を示した。これは、変数がいくつかのグループにまとまるというよりは、噴火による被害が大きいところから小さいところへという変化を示す変数が大半を占め、一部、それとは異なる振る舞いを示す変数があることを意
図(2)-25 共分散構造モデルの結果-I （左）植生の属性と土壌の属性、火山ガス濃度の関係。（中）冗長な変数を省いたモデル。（右）繁殖期の鳥類個体数を組み込んだモデル。

植生の構造と土壌、火山ガスの関係についてのモデル化をまず行ったところ、図(2)-25左図のような結果が得られた。土壌の回復状況に関わる項目はいずれも強い関係を示していた。植生に関わる項目では、胸高断面積がやや弱い関係を示した。主な適合度指標は、\chi^{2}=8.778\left(p=0.553\right)，GFI = 0.828，AGFI = 0.519，RMSEA = 0.000であった。さらに火山ガスを土壌と対応づけようとしても、適応度が下がったり、不適解しか得られなかったりして適切なモデルは得られなかった。

次に動物に関する情報をモデルに組み込んだ。変数の数を抑制するため、前の段階で得られたモデルの簡略化を行い、植生回復、土壌回復ともに二つの変数のみを残した。適応度は上昇し、除外した変数は冗長性が高いものであったことが理解できる(図(2)-25中)。このモデルの主な適合度指標は以下の通り。 $ \chi^{2}=2.764(p=0.429) $ ，GFI = 0.902，AGFI = 0.511，RMSEA = 0.000。

直前のモデルに繁殖期の鳥類個体数を組み込んだ結果を図(2)-25右図に示す。植生のほか土壌からの影響も無視できないことが示された。これは、地表性の小動物を捕食する鳥が相応に存在することを示唆するものと思われる。このモデルに関する主な適合度指標は、以下の通りである。 $ \chi^{2}=4.204\quad(p=0.649) $ 、GFI = 0.879，AGFI = 0.577，RMSEA = 0.000

このモデルに、リターの分解速度を反映する分解定数と、甲虫、ハエ、クモの捕獲個体数を組

図(2)-26 共分散構造モデルの結果-2 （左）繁殖期の鳥類、分解定数、昆虫を加えたもの、（中）越冬期の鳥類を加えたモデル。（右）倒木量を含むモデル。
み込んだモデルも作ることができた(図(2)-26左図)。分解速度は主に植生の発達状況から説明できたが、土壌動物の重量と相互の誤差間に正の相関があるべきと診断された。甲虫、ハエ、クモの個体数は、植生よりも土壌からよく説明できた。このモデルに関する主な適合度指標は、以下の通りである。 $ \chi^{2}=13.127\ (p=0.663) $ 、GFI=0.775、AGFI=0.494、RMSEA=0.000。

越冬期の鳥類の個体数を組み込むモデルも作成した。こちらのモデルでは、下層植生の植被率を取り入れると結果の改善が見られた(図(2)-26中図)。適合度指標は $ \chi^{2}=12.632\left(p=0.396\right) $ 、GFI=0.795、AGFI=0.521、RMSEA=0.076で、RMSEAがやや好ましくない値であった。

枯死木や倒木は、他の調査項目と弱い相関しか示さなかった。そのため枯死木を含むモデルは構築できなかった。倒木を含むモデルとしては、図(2)-26右図に示すものが構築できた。適合度指標は以下の通りである。 $ \chi^{2}=12.463(p=0.490) $ 、GFI=0.761, AGFI=0.485, RMSEA=0.000。以上の結果から、特に次の点を指摘できる。

・火山ガスが現在の植生に強く影響している。今後の生物多様性の保全にあたっては、生物の生息の基礎となる植生の状況や推移の芳香性を的確に把握するため、ガス濃度のモニタリング（または植生指標等による推定）と、ガスにあわせた生態系回復手法の適用が必要である。

・夏鳥の個体数や種数、一部節足動物(甲虫、ハエ、クモ)の個体数は、土壌にも影響を受けている。これは、土壌と関連の深い食物の供給を介した関係と判断できる。

・冬鳥にとっては、樹林地内外の低木や草本植物の意味が大きくなる。

動物の生息状況が土壌と強く関係することを示したことは、重要な成果であると言える。

## （8） 三宅島におけるスダジイ堅果の地域間変異

坪田で、堅果長が有意に短かく、横幅が有意に大きかった。体積指標も坪田で有意に大きかった。縦横比は坪田で有意に小さかった（形態としては幅広）。三宅島という比較的小さい空間スケールでもスダジイの堅果の大きさと形状が地点間で異なっていることがわかった。

## （9） 三宅島のスダジイ堅果における昆虫による食害率の地域間変異

薬師堂での食害率は他の2地点と比べて有意に高かった。つまり大きな堅果ほど昆虫による食害を生き残りやすいという仮説は支持されなかった。坪田のスダジイ堅果が大きい理由は、昆虫による種子捕食では説明できなかった。

## （１０）三宅島におけるオオバエゴノキの種子の昆虫による食害率

三宅島のエゴノキの種子は体積にして本州のものよりも1.5倍大きいことがわかった。食害率は63.7%だった。食害を受けていない種子のほとんどは胚乳の発達がみられないシイナだった。これらの結果と考察から、ヤマガラが昆虫の食害を受けたエゴノキの種子やシイナを地面に落としている可能性が指摘できた。

## （１１）大路池周辺に植栽されたツバキ園芸種の花形態の観察：野生ヤブツバキ集団との交雑の可能性

植栽されたツバキの一部の個体では、葯や雌しべが退化した八重咲きの花が見られた。だが、花中に葯を残し花粉を露出している個体、通常の雌しべを保持している個体も多数見られた。こ
図(2)-27 訪花昆虫の目別個体数

図(2)-28 ハチジョウイタドリ訪花昆虫の目別個体数

のことから、動物による花粉媒介を介した野生集団への遺伝子浸透が懸念される。ヤブツバキ個体群の遺伝的な状況を保全するうえで、ツバキ園芸品種のうち葯をつけ花粉を放出しうる株や雌しべを保持した株については、島外への移植や伐採などの処置が必要であると考えられた。

## （１２）キイチゴ属２種の花の咲く向きと訪花昆虫

カジイチゴにはシラキチビケシキスイが多く訪花し、ハチジョウイチゴにはオドリバエが多く訪花する。野外操作実験の結果、シラキチビケシキスイの訪花は花の咲く向きに関係していた。

## （１３）植物-訪花昆虫間の送粉ネットワークの解明

計 5 回の調査において植物48種から計2440個体の訪花昆虫を捕獲した。捕獲昆虫について目レベルでの分類を行った(図(2)-27)。島内に広く見られるハチジョウイタドリに訪花した個体に着目したところ、環境によって同じ開花植物種でも訪花昆虫の構成や訪花頻度が大きく異なるという傾向が示唆された(図(2)-28)。

## （１４）三宅島における大型鳥類の種子散布者としての機能評価

ハシブトガラスを用いた種子の排出実験からは、果実採食から1時間後にも、排泄される種子は約50％に留まることが明らかになった(図(2)-29)。採食後1時間以内にほとんどの種子を排出するヒヨドリと比較して、かなり長い時間種子を体内に滞留させている。また、採食後7時間以上経過してから排出される種子もあることを確認した。ハシブトガラスの行動圏の広さを考え合わせると、植物種子を非常に広範な範囲に排

泄する可能性があることが示された。カラスバトについては、糞中の種子は認められなかった。

2012年1月・2月に島の北岸にあたる神着地区で採取されたハシブトガラスの糞には、ホルトノキの種子が多く含まれていた。同種の分布の拡大に寄与する可能性がある。

図(2)-29 ハシブトガラス3個体のカラスザンショウ種子排出パターン(400分まで)
## （１５）島内で分布を拡大しているメダケが鳥類の生息に与える影響の検討

越冬期では、メダケ群落で鳥類種数は減少したものの有意ではなかった $ (p=0.18) $ 。繁殖期では、樹木植被率を考慮してもメダケ群落で鳥類種数は有意に少なかった $ (p=0.002 $ 、係数値0.46)。

## 5 ．本研究により得られた成果

## （1） 科学的意義

三宅島内で噴火による被害程度が様々な11箇所の固定調査区で鳥類と昆虫の状況を把握した。三宅島全域で昆虫調査がこのように体系的に行われたことは、過去にない。また、固定調査区以外に広く鳥類の調査を行い、島の鳥類相の現状と空間的な違いを把握した。

噴火直後からの継続調査の結果とあわせて、島内における噴火後の鳥類相の変化の様相を明らかにした。前述の空間的変動の調査結果とあわせ、三宅島2000年噴火による鳥類群集への被害の状況を総合的かつ時空間的に把握することができた。

火山ガス、植生、土壌、鳥類、昆虫の相互関係を、共分散構造分析によりモデル化した。その結果、植生が他の生物や土壌の状況に大きな影響を与えていることがあらためて確認されたほか、以下の新たな知見が得られた。①一部の鳥類および一部の昆虫は、植生よりも土壌の状態により強く影響を受けている、②植生と、それに依存している他の生物は、現状では火山灰よりも火山ガスに専ら影響を受けている、③繁殖期の鳥類は植生のうち特に樹木に依存する程度が大きいが、越冬期の鳥類は下層植生や草本植物に依存する程度が大きい、④鳥類と植生の関係は、植物の種多様性よりも植生の構造に関わるものの重要性が遙かに高い。なお、共分散構造分析を適用してこのような解析を行った事例は、国内で公表されたものはほぼ皆無、海外でも非常に少数である。

島内の鳥類相が、噴火による被害が軽微な箇所も含めて全体的に変容しつつあることを示した。いっときにおける植生-鳥類の関係だけでなく、こうした長期的な変容も考慮に入れた生態系管理の必要性を示すことができた。

噴火前後の現存植生図4時期分を整備した。これらは、行政目的でも有用性が高いが、今後の学術的な調査研究の基礎資料としても広く活用され得るものである。

## （２）環境政策への貢献

## ＜行政が既に活用した成果＞

三宅村営の自然観察施設である「アカコッコ館」が実施する各種活動に協力した（2012年定例会（岸、吉川（ともに本事業に係る特任研究員）、2013年2月自然観察会（吉川（同））。

## ＜行政が活用することが見込まれる成果＞

三宅島には、カラスバト（準絶滅危惧（NT））、ヤマガラ（亜種オーストンヤマガラ、絶滅危惧IB類（EN））、コマドリ（亜種タネコマドリ、絶滅危惧II類（VU））、ミソサザイ（亜種モスケミソサザイ、絶滅危惧IB類（EN））、コゲラ（亜種ミヤケコゲラ）、アカコッコ（絶滅危惧IB類（EN））、イイジマムシクイ（絶滅危惧II類（VU））など、固有性の高い鳥類の種や亜種が生息している。以上に挙げた種はいずれも森林性で、ミヤケコゲラを除いて環境省のレッドリストに掲載されている。本研究の結果、これらは減少傾向にあり、速やかに保全の必要があることが示された。
鳥類の生息に関わる以下の事項は、島内での防災、緑化事業に反映されることが望まれる。

・森林性鳥類保全のための常緑高木林の再生の促進

・鳥類多様性保全のためメダケの繁茂の抑制

・植生の保全・再生にあたり、土壌（一部の鳥類の食物となる土壌中の小動物）も考慮することやブツバキ園芸種の野外への植栽を制限することの必要性についても、あわせて指摘された。

### 6. 国際共同研究等の状況

特に記載すべき事項はない。

### 7. 研究成果の発表状況

## （1） 誌上発表

## <論文（查読あり）>

特に記載すべき事項はない。

<査読付論文に準ずる成果発表>

特に記載すべき事項はない。

## くその他誌上発表（査読なし）>

1）岸茂樹、加藤和弘：Miyakensis、15，9–14（2013）

「三宅島におけるスダジイ堅果の地域間変異」

2）岸茂樹、加藤和弘：Miyakensis、15，15–24（2013）

「三宅島のスダジイ堅果における昆虫による食害率の地域間変異」

3）岸茂樹、加藤和弘：Miyakensis、15，25–30（2013）

「三宅島におけるオオバエゴノキの種子の昆虫による食害率」

4）吉川徹朗、加藤和弘：Miyakensis、受理

「三宅島・雄山上部の植生荒廃地で観察されたカラスバト」

5）吉川徹朗、樋口広芳：Miyakensis、受理

「大路池周辺に植栽されたツバキ園芸種の花形態の観察：野生ヤブツバキ集団との交雑の可能性」

## （2） 口頭発表（学会等）

1）岸茂樹・吉川徹朗・加藤和弘：個体群生態学会大会第28回大会（2012）

「三宅島におけるキイチゴ属2種の花の咲く向きが訪花昆虫に与える影響」

2）岸茂樹：日本生態学会第60回大会（2013）

「訪花昆虫群集をオスとメスにおける」

3）加藤和弘・吉川徹朗・樋口広芳：日本生態学会第60回大会（自由集会）（2013）

「三宅島2000年噴火後の森林における鳥類相の変遷」

4）高橋俊守・指村奈穂子・加藤和弘：日本生態学会第60回大会（自由集会）（2013）

「三宅島2000年噴火後のランドスケープの変化」

5）吉川徹朗：日本野鳥の会東京支部定例会(2013)

「鳥類と植物の複雑な関係を紐解く」
6）高橋俊守・指村奈緒子・加藤和弘：日本生態学会第61回大会（自由集会）（2014）

「2000年噴火に伴うランドスケープ構造の変化」，広島，2014年3月.

7）加藤和弘・岸茂樹・吉川徹朗・高橋俊守・上條隆志・廣田充・田村憲司・樋口広芳：日本生態学会第61回大会（自由集会）（2014）

「2000年噴火が生態系に与えた影響の総合的な把握に向けて」

8）岸茂樹：日本生態学会第61回大会(企画集会)(2014)

「見落とされてきた訪花昆虫のオス」

## （3） 出願特許

特に記載すべき事項はない。

## （4） シンポジウム、セミナー等の開催（主催のもの）

1）日本生態学会自由集会「2000年噴火後の三宅島における生態系研究 現状と今後に向けた課題」（平成25年3月5日、静岡県コンベンションアーツセンター、参加者約45名）

2）日本生態学会自由集会「2000年噴火後の三宅島島の生物はどうなったか」（平成26年3月14日、広島国際会議場、参加者約50名）

## （５）マスコミ等への公表・報道等

特に記載すべき事項はない。

## (6) その他

特に記載すべき事項はない。

### 8 \. 引用文献

1）高橋俊守・加藤和弘・上條隆志（2011）衛星リモートセンシングによる三宅島2000年噴火後の植生モニタリング．日本生態学会誌61：167-175.

2) ter Braak, C. J. F. (1986) Canonical Correspondence Analysis: A New Eigenvector Technique for Multivariate Direct Gradient Analysis. Ecology 67:1167–1179.

3) Ter Braak C. J. F. & Šmilauer P. (2012): Canoco reference manual and user's guide: software for ordination, version 5.0. Microcomputer Power, Ithaca, USA, 496 pp.

4）槙原寛・岡部宏秋(2006) 三宅島噴火4,5年後のカミキリムシ相--大発生したカミキリムシ.月刊むし 419:32-41.

5) Abe H, Hasegawa M 2007. Impact of volcanic activity on a plant-pollinator module in an island ecosystem: the example of the association of Camellia japonica and Zosterops japonica. Ecological Research 23(1): 141-150.

6）七島花の会（2011）神津島花図鑑．日本出版ネットワーク．

7）加藤和弘・樋口広芳（2003）三宅島2000年噴火後の島の森林における鳥類群集. Strix 21:81-98

8）加藤和弘・樋口広芳（2011）三宅島2000年噴火後の植生の退行・回復にともなう鳥類群集の変化．日本生態学会誌61：177-183

9）加藤和弘・樋口広芳（2006）三宅島2000年噴火・鳥類への影響と回復．森林科学46：16-19

10）宫下直・野田隆史．2003．群集生態学．東京大学出版会．
GeoEye-1 20130428