台灣地區青斑蝶族群變動與移動行為之探討

(1)

國立台灣大學生物資源暨農學院昆蟲學研究所碩士論文

Graduate Institute of Entomology College of Bioresources and Agriculture

National Taiwan University Master Thesis

台灣地區青斑蝶族群變動與移動行為之探討 Study on the population dynamics and moving behavior

of the chestnut tiger butterfly (Parantica sita niphonica Moore) in Taiwan

李信德 Hsin-Te Lee

指導教授：楊平世博士

Advisor: Ping-Shih Yang, Ph.D.

中華民國 98 年 6 月

June, 2009

(2)

口試委員會審定書

(3)

ii

誌謝

天

(4)

中文摘要

本研究以標識再捕法分析台灣地區青斑蝶族群的組成與周年動向，並就陽明山大屯山區的青斑蝶再捕獲資料進行族群數量估算。2000-2008年間，大屯山區青斑蝶族群數量最低為418隻，最高為165322隻，年度間數量變化差異極大。就長時間的連續變化而言，此區族群的數量變動可能以3-4年為一週期。就異地再捕獲的長距離移動與各地區的出現時間紀錄，青斑蝶有季節性、固定方向的移動行為，

但短距離的移動則似無方向規律。冬季時無群聚越冬現象，雌蝶亦無生理滯育，

於多處地區可觀察到成蝶活動與幼蟲的世代紀錄，其越冬方式可能為多策略並行。

關鍵詞：青斑蝶、族群變動、移動行為、標識再捕法、族群估算、陽明山國家公園

(5)

iv

英文摘要

Study on the population dynamics and moving behavior of the chestnut tiger (Parantica

sita niphonica Moore) in Taiwan

by Hsin-Te Lee Abstract

Chestnut tiger butterfly（Parantica sita niphonica Moore）population dynamic was studied by using mark-release-recapture method of Yangmingshan national park, Datun mountain. The results from 2000 to 2008 indicated that the estimated chestnut tiger butterfly population size ranged from 418 to 165322. Thus the population size of the butterfly varied diffentially year by year. However, there seems to be a tendency for the population size of chestnut tiger butterfly to undergo a 3-4 year cycle. According to the records of recapture in different sites, distances and areas, chestnut tiger butterfly may have seasonal and directional movement. The short distance movements, on the other hand, had no specific orientation. The chestnut tiger butterfly may employ various strategy to overwinter. The chestnut tiger butterfly did not gather in the winter, and all stages of larvae and adult could be found in many areas.

Key words: Parantica sita niphonica Moore、population dynamics、moving behavior、

mark-release-recapture method、population estimation、Yangmingshan national park

(6)

口試委員會審定書...i

誌謝...ii

中文摘要...iii

英文摘要...iv

目錄...v

表次...vi

圖次...vii

壹、緒言...1

貳、往昔研究...2

一、青斑蝶概述...2

二、日本地區的青斑蝶研究...2

(一) 春季北上...3

(二) 夏季山區的移動...3

(三) 秋季南下...3

三、台灣地區青斑蝶的周年活動與分布...4

参、材料與方法...6

一、研究地點...6

(一) 陽明山國家公園大屯山區...6

(二) 南澳觀音產業道路...7

(三) 其他青斑蝶活環境...7

二、研究方法...8

(一) 標識再捕法...8

(二) 青斑蝶族群絕對豐度估算...10

(三) 島內與島嶼間青斑蝶移動追蹤...12

(四) 氣象資料收集...12

肆、結果...14

一、青斑蝶個體資料...14

(一) 青斑蝶族群組成...14

(7)

vi

(二) 觀音樣區青斑蝶周年分布...15

二、族群估算...16

三、再捕獲移動紀錄...36

(一) 島外移動記錄...36

(二) 島內移動記錄...36

伍、討論...42

一、台灣地區青斑蝶族群組成與季節動向...42

(一) 陽明山地區的青斑蝶族群...42

(二) 南澳地區的青斑蝶族群...43

二、台灣地區青斑蝶的越冬狀態...44

三、族群估算與影響因素...45

四、台灣地區青斑蝶移動與遷移的可能性...46

(一) 標識再捕紀錄的方向性...47

(二) 遷移的生理條件...48

五、青斑蝶標識研究與推廣教育的關連...49

陸、結論...51

柒、參考文獻...52

捌、附錄...59

(8)

圖次

圖一、陽明山國家公園大屯山研究樣區...6

圖二、青斑蝶標識代號說明...9

圖三、青斑蝶翅膀磨損等級區分...10

圖四、氣象測候站位置示意圖...13

圖五、2005 年十一月至 2007 年九月觀音樣區青斑蝶成蝶數量與卵數...17

圖六、2000 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...18

圖七、2001 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...19

圖八、2002 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...20

圖九、2003 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...21

圖十、2004 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...22

圖十一、2005 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...23

圖十二、2006 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...24

圖十三、2007 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...25

圖十四、2008 年陽明山大屯山區青斑蝶族群估值與標記數量的比較...26

圖十五、2000 年-2008 年間，陽明山大屯山區青斑蝶族群最高估值...27

圖十六、2000 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較...28

圖十七、2001 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較...29

圖十八、2002 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較...30

圖十九、2003 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較...31

圖二十、2004 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較...32

圖二十一、2005 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較....33

圖二十二、2006 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較....34

圖二十三、2007 年六月，陽明山大屯山區青斑蝶族群估值與氣候資料的比較....35

圖二十四、北向再捕獲紀錄 1 與 2...36

圖二十五、北向再捕獲紀錄 3...36

圖二十六、北向再捕獲紀錄 4...37

(9)

viii

圖二十七、南向再捕獲紀錄 1...37

圖二十八、南向再捕獲紀錄 3...37

圖二十九、南向再捕獲紀錄 4 與 5...38

圖三十、南向再捕獲紀錄 6 與 7...38

圖三十一、島內再捕獲記錄 1...38

圖三十二、島內再捕獲記錄 2...39

圖三十三、島內再捕獲記錄 3...39

圖三十四、16 方位角示意圖...39

(10)

表次

表一、南澳鄉觀音產業道路標定寄主植物之 GPS 位置...8

表二、氣象測候站之 GPS 位置...12

表三、2000-2008 年，台灣地區青斑蝶標識數量與雌雄組成比例...14

表四、2000-2008 年，陽明山地區青斑蝶標識數量與雌雄組成...15

表五、2000-2008 年，陽明山以外地區青斑蝶標識數量與雌雄組成...15

表六、2000-2008 年，陽明山大屯山區青斑蝶標識再捕回收率...16

表七、2000 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...18

表八、2001 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...19

表九、2002 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...20

表十、2003 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...21

表十一、2004 年陽明山大屯山區青斑蝶標識再捕與族群估算結果... 22

表十二、2005 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...23

表十三、2006 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...24

表十四、2007 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...25

表十五、2008 年陽明山大屯山區青斑蝶標識再捕與族群估算結果...26

表十六、2000-2008 年青斑蝶異地再捕獲整理...40

表十七、2005 年六月 13 日至七月 2 日間，大屯山與二格山的連續風向變化...41

(11)

1 壹、緒言

青斑蝶 (Parantica sita niphonica Moore) 屬於蛺蝶科斑蝶亞科，為台灣地區常 見種類的蝴蝶，廣泛分布於全島各地。因其食草為蘿藦科 (Asclepiadaceae) 的有毒植物，故蟲體內具有毒性作為防禦之用。所以成蝶天敵較少，亦較不畏懼人為干擾，是適合的賞蝶對象。

在台灣北部陽明山國家公園境內，每年春夏季節均有大量的青斑蝶出現，尤以大屯山區的大發生景象最為壯觀。大發生期間蝴蝶數量在短時間內極劇的增加，而期間過後卻又急速減少，2000 年曾首次証實此區的蝴蝶有跨海長距離飛行至日本的行為 (Lee and Yang，2000)，其後陸續的再捕獲紀錄也證實台日之間青斑蝶的族群交流。

而以往對於此種蝴蝶的研究多屬於描述其生物學的範疇，族群的數量變化與青斑蝶的飛行移動資料較少，雖有學者曾針對族群數量研究，但缺乏長期且持續的觀察，對欲了解台灣地區青斑蝶族群的變動資料仍嫌不足。本研究以標識再捕法為主，對台灣地區的青斑蝶族群進行調查，希望對青斑蝶的族群變動建立基本資料，並對青斑蝶在台灣島嶼內及島嶼間的季節動向與飛行路徑能有初步了解。

(12)

貳、往昔研究

一、青斑蝶概述

青斑蝶 (Parantica sita) 屬於鱗翅目蛺蝶科斑蝶亞科。據 Shirozu (1982)，本種 分為四個亞種，其中 Parantica sita niphonica Moore 分布在韓國、日本及台灣，為 斑蝶亞科種類中分布最北者。

據楊 (1989)、Hsu (1999)、Lee (1990)、Yang (1990)，青斑蝶 (P. sita niphonica) 成蟲前翅長約 50-60mm，前翅略成鈍三角形，底色為黑色；後翅近似圓形，後翅為暗紅褐色。卵成白色，砲彈型，上有淺刻孔及縱稜。幼蟲分為五齡，一齡幼蟲頭黑色而體成青灰色，成長的幼蟲體表底色為黑色，其上佈有白色斑點，各體節左右各有兩個大黃斑。中胸及第 8 腹節背部左右有一對黑色，細長肉質突起。蛹為翠綠色懸蛹，第三腹節背側有數個黑圓點排成一環。中胸背側、翅基與頭部均有銀色斑紋。

Lee (1990)，記錄青斑蝶生活史：25℃下，卵期約 3 日，幼蟲期約 13 日，蛹期約 8 日。

二、日本地區的青斑蝶研究

據福田 (2000)，Miyatake et al. (2003): 長嶺邦雄曾於 1962 年提出沖繩地區包 含青斑蝶在內的六種蝶類，數量在夏季有突然減少的現象。而後，田中洋根據 1963 與 1964 年在九州南部鹿兒島市內週年的紀錄，認為青斑蝶有移動的可能性。1980 年，長嶺邦雄於沖繩本島觀察到族群數量為數百～上千的青斑蝶群，於兩週內數量急遽減少到僅有個位數的紀錄，亦因此提出青斑蝶有移動行為的假設。

日本地區最初的青斑標放起於 1975 年，中村耕人於九州宮崎縣山區標放四百餘隻青斑蝶，嘗試追蹤其垂直移動路徑。而正式的青斑蝶標識再捕研究則始於五年之後的 1980 年，由九州鹿兒島縣的鹿兒島昆蟲同好會所發起；隔年五月，日本關西地區的三重縣再捕獲由九州南方的種子島所標記的青斑蝶，此為首例青斑

(13)

3

蝶北上飛行的再捕獲紀錄；而首例南下飛行再捕獲紀錄則為同年十一月，於九州南方的奄美大島再捕獲由鹿兒島縣內標記的青斑蝶個體。

日本青斑蝶標放研究中主要的研究團體，為大阪地區的大阪青斑蝶調查會，

以及九州鹿兒島縣青斑蝶大移動研究會分別成立於 1983 與 1984 年；大阪青斑蝶調查會並於 1990 年更名為青斑蝶調查會。1994 年起，利用網際網路架設網站於留言版系統，加速青斑蝶標識資訊的流通與擴展參與標識研究的人員階層。當時最接近台灣地區的再捕獲紀錄，為 1995 年沖繩縣與那國島所再捕獲，由關西大阪市境內所邊放的青斑蝶，與台灣本島的宜蘭縣最近距離僅約 110km。

日本地區青斑蝶的移動現象可分為三個時期：

(一)、春季北上

春季到夏初 (三～七月) 西南季風盛行，由越冬地區新羽化的青斑蝶世代，

可能利用西南氣流由九州南方各島開始往日本北方移動，但因日本大部分地區中，春季時非青斑蝶的主要活動季節，因此標識再捕的紀錄較少。

(二)、夏季山區間的移動

七月～八月下旬是日本各地青斑蝶主要活動時間，因為幼蟲主要的食草與成蝶的蜜源植物生長於山區，因此成蝶出現地點也以山區或高地為主；一般認為，春季北上的個體向北遷移後也會往山區移動。夏季在各山區標識再捕的結果，鮮少有長距離移動的情形，成蝶僅在數公里之內的山系間移動。

(三)、秋季南下

夏末秋初，日本各地青斑蝶均有明顯移動現象。小地區範圍內，青斑蝶由山區向平地移動；廣域地區範圍，本島的青斑蝶族群順著當時季風風向，往南方朝九州方向飛行。南向飛行的地點不僅止於九州本島，而是沿著東亞島弧向南方前進，沿薩南諸島、奄美諸島、琉球群島，最遠可到達八重山群島中的与那國島，離台灣本島極近。因此，亦推論台灣可能為青斑蝶長距離長距離的停

(14)

留點。

遷移行為發生的可能原因，包括短日照、低溫、季節氣流方向、寄主植物老化枯死，以及蜜源植物減少等環境因子，日方研究者並就台灣與日本的地理位置，

推測台灣地區的青斑蝶族群亦有遷移行為 (Wei, 1995)。

三、台灣地區青斑蝶的周年活動與分布。

青斑蝶在臺灣的分布記錄散見於諸多報告書籍，早期詳細分佈資料為

Yamanaka (1973) 「台灣產蝶類的分布」，據其記載，青斑蝶分布於全台各地，主要出現時間為春季與夏季。

在青斑蝶的垂直分佈方面，由平地至海拔 3000 公尺以上的環境仍有活動記錄，但在 200 公尺以下較常見。雖在台灣各地皆有分布，但北部的數量較中南部為多，春夏季節多處可見青斑蝶活動，其中尤以陽明山國家公園境內的青斑蝶大發生現象最受矚目。青斑蝶大發生時間在五月～六月間，實際族群數量於 1989～

1990 年的估算，大屯山區的青斑蝶族群量可由數萬至十餘萬隻不等，每年的發生量均有差異。

大屯山與面天山的森林中，有青斑蝶幼蟲的寄主植物，可能是主要的幼蟲繁殖地；也被認為是陽明山青斑蝶大發生的族群來源。而 1989～1990 年，統計青斑蝶的翅膀狀況，大發生期間，翅膀為全新或是輕微磨損者佔了族群的大部分，顯示此時青斑蝶多為新生個體，且並無因為長時間或長距離飛行而造成翅膀磨損的情形。因此，亦得到陽明山區大發生時的青斑蝶族群，多數應為當地區所繁衍，

而非由較遠的外地遷移進入的相同結論。(林等，1986；陳，1981；楊，1989；魏與楊，1991；Wei，1995；Wei and Yang，1990；Yang，1990)

秋冬季節，十一月至隔年三月是青斑蝶於南部的主要發生時期(Lee，1990)；

Hamano (1987)，亦提到秋冬季節青斑蝶群聚於南部六龜地區。

但在南部諸多紫斑蝶類與青斑蝶類越冬的谷地，青斑蝶均無出現在越冬集團

(15)

5

內。Ishii and Matsuka (1990)，調查台灣西南部山區的斑蝶類成蝶越冬族群中，紀錄到 10 種斑蝶亞科種類，其中無青斑蝶紀錄。Lee and Wang (1997)，調查台灣地區冬季的蝴蝶越冬谷地，所記錄到的九種斑蝶亞科種類中也無觀察到青斑蝶。Chao and Chen (2005、2007)、Chao et al. (2007)，於台東大武地區越冬斑蝶族群的調查 裡，亦無發現青斑蝶成蝶有越冬現象。但 Hsu et al. (1998)、Hsu and Yang (2000)，

於台東紅葉村台東蘇鐵自然保留區，在 1988 便曾發現此處有斑蝶的越冬族群，且其中有青斑蝶混雜，但並未描述當時青斑蝶個體的活動狀況。

陳 (2000)、Chen (2005a)，於 1999-2000 年間冬季時，於陽明山國家公園觀察到低溫情況下幼蟲世代仍能繁衍，因此認為青斑蝶並非如同其他越冬斑蝶類是僅有成蝶越冬狀態，而可能為多樣化的輻射適應。

2000 年兩例大屯山標記個體於日本再捕獲的記錄，証實除了台灣島內的移動，青斑蝶亦有跨海長距離移動飛行的能力，此點與日方研究的推論相同

(Fukuda，2000；Lee and Yang，2000；Uchida and Lee，2004)。2001 年至 2008 年，

台日間陸續的再捕獲記錄亦証明兩地青斑蝶族群有密切的交流，且飛行的方向與 時間似乎具有關聯 (Fujii，2007；Lee and Yang，2005；Ohshima et al.，2002) 但標識再捕紀錄之外，青斑蝶直接的動向觀察資料極少，除了 Chen (2002、

2003、2004、2005b、2006、2007、2008b) 於玉山國家公園觀察到青斑蝶秋冬季時，

由塔塔加北部水里方向翻越塔塔加鞍部，而往南邊的楠梓仙溪飛行外，其它地區均少有青斑蝶的動向觀察。吳 (2001)，於台灣北海岸觀察到斑蝶成群移動現象，

但其中亦並無青斑蝶紀錄。

(16)

参、材料與方法

一、研究地點

(一)、陽明山國家公園大屯山區

調查地點為大屯主峰登山車道沿線。

據 Wei (1995)，大屯主峰登山車道屬於陡升車道地形，海拔高度由 872 公尺上升至 1092 公尺，車道寬度約 4 公尺。前半段為紅楠與台灣矢竹草原植物社會，後段為台灣芒草與台灣矢竹草原植物社會。本段路面植被遮蔽度少，

車道旁大量孳生斑蝶類主要的菊科蜜源植物－島田氏澤蘭。車道兩旁每年均有固定進行除草，除草影響蜜源植物生長，而近年除草的時間已延後至斑蝶類大發生過後進行。

圖一、陽明山國家公園大屯山研究樣區。

Fig. 1. Study area in Datung Mt., Yangmingshan National Park.

(17)

7

自 2000 年～2008 年，視各年青斑蝶發生狀況而略有差異，於四月下旬

～五月上旬，青斑蝶開始出現起，每週兩次，每次間隔二～三日，於大屯山登山車道沿線進行青斑蝶成蝶標識。標識進行持續到青斑蝶族群數量明顯減少為止，約持續進行至七月中旬。其它月份則每月進行一次調查。

(二)、南澳觀音產業道路

調查區域位於省道台九線蘇花公路 144.5K 處，據楊與李(1995)，為觀音海岸野生動物重要棲息環境範圍內。此區地形南北狹長，東向陡坡急降入海，

溼度高，全年多雨。此區林相主要為台灣赤楊、大葉楠、茄冬、樟科與榕屬的植物組成；在裸露或陡坡面則為五節芒的優勢族群。

於調查全島青斑蝶分布狀況時，發現本區有大量的青斑蝶寄主植物台灣 牛嬭菜 (Marsdenia formosana)，且多位於道路兩旁的林班地內，少遭砍伐，

觀察容易，有穩定的成蟲與幼蟲數量，故於此以 GPS 標定 16 株台灣牛嬭菜,，

如表一所列。並於 2005 年十一月至 2007 年九月止，每周一次，調查並記錄標定寄主植物上的青斑蝶幼生期，包括卵、幼蟲、蛹的數量與寄生狀況。

成蝶的調查時間則由 2005 年十一月至 2007 年九月止。除了此區，亦同時將東澳粉鳥林、烏石鼻自然保留區、南澳金岳村等地的成蝶觀察紀錄一並列入。

(三) 其他青斑蝶活動環境

於 2000-2008 年間，每月一次，於島內青斑蝶主要的分布環境，調查青斑蝶成蝶數量並進行標記。

(18)

表一、南澳鄉觀音產業道路標定寄主植物之 GPS 位置

Table 1. GPS coordinates of host plants (Marsdenia formosana) in Nan-ao Township.

No. of sample Elevation (m)

Coordinate (TW67)

1 471 329075 2698464

2 503 329035 2698276

3 503 329035 2698276

4 565 328391 2698364

5 568 328407 2698382

6 568 328407 2698382

7 451 329153 2698336

8 451 329153 2698336

9 451 329153 2698336

10 471 329075 2698464

11 541 328650 2698325

12 409 329217 2698349

13 572 328413 2698388

14 549 328391 2698346

15 558 328478 2698369

16 471 329075 2698464

二、研究方法 (一)、標識再捕法

據 Southwood (2000)，使用標識再捕法估算族群有四項假設：

1、標識個體不因標識影響其行為或壽命，且標識記號不會脫落。

2、標識個體可完全混合於族群中。

3、標識個體與族群中其他未標識個體的被捕獲機率相同。換言之，標識進行為隨機取樣，與標識對象的狀態、年齡、性別等無關。而「相同再捕機率」亦有兩項要點：

(1)、不同年齡與性別的所有個體能依照其所佔比例而被選取標識。

(2)、所有個體被捕獲的機率相同，與它們在棲地中的分布位置無關。

(19)

9

4、取樣須在分離的時間段進行，且標識進行的時間與研究進行時間相較之下須佔極少比例。

標識方式以油性黑色細字簽字筆，於青斑蝶後翅中室書寫標記代碼後原地釋放。標識所使用的標記代碼由四部份組成，如圖二所示：

1、編號：第一隻青斑蝶標識為 001，第二隻標識為 002，之後依序類推。

2、標識地點：以標識地的英文縮寫表示，如陽明山國家公園 (Yang-Ming-Shan National Park)縮寫表示為 YMS。

3、標識單位：以標識者所屬單位的英文縮寫表示，如台灣大學 (National Taiwan University) 的縮寫表示為 NTU。

4、標識日期：日期以四位代碼表示，如 4 月 13 日則記為 0413。但在陽明山地區，因大發生期間青斑蝶數量過多，為節省時間，首次標識日，

日期代碼改註記為「A」，第二次標識日註記代碼為「B」之後依序類推，以圖一為例，Ｄ表該年度第四次的標識日。

圖二、青斑蝶標識代號說明 a. 069 為個體編號；D 為標識日代碼；NTU 為標識者代碼。b 為標識地陽明山國家公園的縮寫。

Fig. 2. Marking codes of Parantica sita niphonica Moore.

a. 069: serial numbers; D: date; NTU: code name of researcher.

b. YMS: code name of study area.

標識再捕進行時，亦同時紀錄青斑蝶個體資料。包括：

a b

(20)

1、前翅長度：由前翅基部的白點量至翅末端。

2、翅膀磨損程度：依翅膀破損與鱗粉掉落的程度分為新、中、舊三等級 (圖三)。翅膀的磨損程度與蝴蝶的壽命有關，新羽化的年輕個體，翅膀與鱗粉構造完整；而存活較久個體，則會因為飛行或受到天敵攻擊等因素，導致鱗粉掉落或是翅膀破損。以青斑蝶的翅膀磨損情況可以反應出族群的年齡組成結構。

3、行為紀錄：青斑蝶被捕捉當時的行為，如飛行、停棲、訪花等，若為訪花則亦同時紀錄訪花蜜源種類。

圖三、青斑蝶翅膀磨損等級區分。a. 新鮮(N)；b. 稍有磨損(M)；c. 磨損嚴重(O)。

Fig. 3. Wing condition of P. sita niphonica. a. new (N)、b. median (M)、c. old (O).

再捕獲紀錄：於各研究樣區進行標識時，無論同地或異地再捕獲，皆於蝶前翅追加標記再捕獲當日的日期或日期代碼後釋放，並以拍攝照片作為紀錄，若有需要，則將再捕獲的標記青斑蝶攜回製成乾燥標本保存。所有的青斑蝶再捕獲資料，整理後進行族群估算與推測飛行路徑。

(二) 青斑蝶族群估絕對豐度估算

族群估算乃根據標識再捕法所得到的數據，估算該族群的絕對豐度。據 Southwood (2000)，族群估算區分為封閉式族群 (closed populations) 與開放式族群 (open populations)。

封閉族群的估算假設前提為族群在研究進行當中不會有任何變動，任何個體均無法離開研究樣區或由樣區外進入；而開放族群則會因研究期間的出生、

(21)

11

死亡、遷入、遷出等因素而產生變動。

而常用於開放式族群的估算有三種方法：

1、Fisher and Ford trellis method 2、Bailey`s triple catch

3、Jolly-Seber stochastic method

其中 Jolly-Seber stochastic method 為最常使用的方法，基本公式為：

_i ⁱ ⁱ

i

N M m

=

r

Ni 為 i 日的族群估值。

Mi 為 i 日標記總數量的估值。

ri 為 i 日再捕獲標記個體的總數量。

ni 為 i 日捕獲的總數量。

而 Mi 的估算公式為： _i ⁱ ⁱ _i

i

M a Z r

=

R

+ ai 為 i 日所標記的個體數。

Ri 為 i 日所標記，而於其後標示日被再捕獲的總數。

Zi 為在 i 日前標記，但不在 i 日被捕獲，而在之後被再捕獲的總數。

而 i 日中，標記個體在族群裡的比例為 _i ⁱ

i

r α

=

n

所以原來的估算公式可表示為： _i ⁱ

i

N M

=

α

Mi 帶入後估算公式變成為： _i ⁱ _i ⁱ ⁱ ⁱ

i i i

M n a Z

N n

α R r

= = +

其中，^var

( )

ⁱ ⁱ

(

ⁱ ⁱ

)

ⁱ ⁱ ⁱ ¹ ¹ ¹ ⁱ

i i i i

M r a

N N N n

M R a r

⎧ − + ⎡ ⎤ −

α

⎫

⎪ ⎪

= − ⎨ ⎢ − ⎥+ ⎬

⎪ ⎣ ⎦ ⎪

⎩ ⎭

(22)

楊 (1989)、Yang (1990)、Wei (1995)、Wei and Yang (1990)等，以此法針對 陽明山青斑蝶族群進行估算；Chao et al. (2008)亦以用此法對大武地區越冬的紫 斑蝶類群聚進行數量估算。

(三) 島內與島嶼間青斑蝶移動追蹤

調查進行期間，所標識之青斑蝶若於他地被再捕獲者，則紀錄再捕獲地點、再捕獲日期、再捕獲者，以及前述各項個體資料。

(四) 氣象資料收集

由中央氣象局申購氣象資料，測候站的選擇的方式：

1、標識地點周邊測候站，用以分析氣候因子對當地族群數量的影響。

2、標識地點與再捕獲地點直線距離之間的測候站，用以分析氣候因子對青斑蝶個體飛行的影響。

收集的風速、風向、雨量、溫度等氣象資料，所使用的測候站 GPS 座標與地理位置如表二與圖四所示。

表二、氣象測候站之 GPS 位置

Table 2. GPS coordinates of meteorological stations.

Code Toponym County Elevation (m) Coordinate (TW67)

466910 鞍部臺北市 825.8 302428 2786547

466920 臺北臺北市 5.3 301131 2770270

466930 竹子湖臺北市 607.1 304059 2784169

C0A540 四堵臺北縣 401 324508 2754235

C0A530 坪林臺北縣 300 320834 2759293

C0A570 桶後臺北縣 387 309796 2749302

C0A9B0 石牌臺北市 7 300982 2778958

C0A9C0 天母臺北市 63 303334 2779121

C0A9F0 內湖臺北市 20 307050 2774890

C0A9G0 南港臺北市 42 310004 2772287

C0A9H0 信義臺北市 51 306400 2769318

C0AC40 大屯山臺北市 1098 301825 2785515

(23)

13

圖四、氣象測候站位置示意圖。

Fig. 4. Locations of meteorological stations.

(24)

肆、結果

一、青斑蝶個體資料

(一) 青斑蝶族群組成

表三為 2000 至 2008 年間台灣地區青斑蝶標記數量的統計，研究期間共標記 25561 隻青斑蝶成蝶，其中雄蝶 23779 隻，雌蝶 1782 隻。陽明山地區標記 24359 隻青斑蝶，佔全部的 95.3％。

表四的統計結果顯示，青斑蝶族群中，雄性佔絕對多數。而陽明山青斑蝶族群中，雄性比例高達 93.74％。陽明山以外的青斑蝶族群，雄性僅佔 78.54％，

雌性則占 21.46％，陽明山以外的青斑蝶族群，雌性所佔比例略有提升。

表三、2000-2008 年，台灣地區青斑蝶標識數量與雌雄組成比例。

Table 3. Marked numbers and sex ratio of marked Parantica sita niphonica from 2000 to 2008 in Taiwan.

Year Numbers of total marked

Male

Numbers Percent

Female

Numbers Percent

2000

3994 3682 92.19% 312 7.81%

2001

3199 3046 95.22% 153 4.78%

2002

957 892 93.21% 65 6.79%

2003

1711 1582 92.46% 129 7.54%

2004

2118 1953 92.21% 165 7.79%

2005

5597 5104 91.19% 493 8.81%

2006

2387 2255 94.47% 132 5.53%

2007

1373 1254 91.33% 119 8.67%

2008

4225 4011 94.93% 214 5.07%

Total

25561 23779 93.03% 1782 6.97%

(25)

15

表四、2000-2008 年，陽明山地區青斑蝶標識數量與雌雄組成。

Table 4. Marked numbers and sex ratio of marked P. sita niphonica from 2000 to 2008 in Yangmingshan National Park.

Year Numbers of total marked

Male

Numbers Percent

Female

Numbers Percent

2000

3994 3682 92.19% 312 7.81%

2001

3166 3013 95.17% 153 4.83%

2002

957 892 93.21% 65 6.79%

2003

1519 1436 94.54% 83 5.46%

2004

1966 1830 93.08% 136 6.92%

2005

5499 5038 91.62% 461 8.38%

2006

2252 2150 95.47% 102 4.53%

2007

1275 1189 93.25% 86 6.75%

2008

3731 3605 96.62% 126 3.38%

Total

24359 22835 93.74% 1524 6.26%

表五、2000-2008 年，陽明山以外地區青斑蝶標識數量與雌雄組成。

Tab. 5. Marked amounts and sex ratio of marked P. sita niphonica from 2000 to 2008 in other places.

Year Numbers of total marked

Male

Numbers Percent

Female

Numbers Percent

2000

- - - - -

2001

33 33 100.00% 0 0.00%

2002

- - - - -

2003

192 146 76.04% 46 23.96%

2004

152 123 80.92% 29 19.08%

2005

98 66 67.35% 32 32.65%

2006

135 105 77.78% 30 22.22%

2007

98 65 66.33% 33 33.67%

2008

494 406 82.19% 88 17.81%

Total

1202 944 78.54% 258 21.46%

(二) 觀音樣區青斑蝶周年分布

圖五為卵與成蝶的調查結果。自 2005 年十一月～2007 年九月.，於南澳鄉觀音產業道路的食草標定樣區，共進行 53 次調查，卵與幼蟲各齡期紀錄列於

(26)

附錄三，其中 2006 年五月 23 日～2006 年九月 13 日，及 2007 年五月 4 日～2007 年九月 12 日間，於此區寄主植物上並無幼蟲與卵的記錄。

成蝶調查時間自 2005 年十一月～2007 年九月，共進行 42 次調查，其中 2006 年五月 2 日～2006 九月 13 日，及 2007 四月 13 日～2007 年九月 12 日無成蝶的紀錄。

二、族群估算

表六為 2000 至 2008 年間，陽明山區青斑蝶標識再捕的回收率，最低為 2001 年 2.18％，最高為 2004 年 5.20％。九個年度的標識再捕回收率均為偏低情形，沒有達到能建低估算誤差的 9％門檻。

表六、2000-2008 年，陽明山大屯山區青斑蝶標識再捕回收率。

Table 6. Recapture ratio of P. sita niphonica from 2000 to 2008 in Yangmingshan National Park.

Year Numbers of total marked

Total recapture

Numbers Percent

2000

3994 134 3.36%

2001

3165 69 2.18%

2002

957 23 2.40%

2003

1513 38 2.51%

2004

1944 101 5.20%

2005

5499 228 4.15%

2006

2252 63 2.80%

2007

1275 42 3.29%

2008

3700 96 2.59%

Total 24299 794 3.27%

表七至表十五為 2000 至 2008 年間，以 Jolly-Seber Method 就陽明山大屯山區進行的標識再捕法所得到的再捕獲資料，進行青斑蝶族群的絕對豐度估算結果。

各年內與各年度間估算數量差異頗大。圖六至圖十四為各年標記數量與族群估值的比較圖，各年度內，標記數量與該次族群估值並無明顯相關。

(27)

17 0 10 20 30 40 50 60

2005 Nov. 2005 De c. 2006 Jan. 2006 Fe b. 2006 Ma r. 2006 Apr. 2006 Ma y 2006 Jun. 2006 Jul. 2006 Aug. 2006 Se p. 2006 Oc t. 2006 Nov. 2006 De c. 2007 Jan. 2007 Fe b. 2007 Ma r. 2007 Apr. 2007 Ma y 2007 Jun. 2007 Jul. 2007 Aug. 2007 Se p.

No. of adults

0 20 40 60 80 100 120 140

No. of eggs

No. of adults No. of eggs

圖五、2005 年十一月至 2007 年九月觀音樣區青斑蝶成蝶數量與卵數。

Fig. 5. Numbers of adults and eggs from Nov. 14, 2005 to Sep. 13, 2007 in Nan-ao Township.

17

(28)

表七、2000 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Table 7. Population estimation of P. sita niphonica in Yangmingshan National Park, 2000.

Date Captured Marked Recaptured Population size Eatimated S.E.

Jun. 15

70 70 - - -

Jun. 17

418 417 1 418.0 0.0

Jun. 19

1050 1042 8 37520.0 17437.7

Jun. 22

569 550 19 26648.2 8920.4

Jun. 26

744 712 32 36216.9 12948.1

Jun. 29

466 446 20 16053.7 6031.3

Jul. 3

48 41 7 2367.4 1425.1

Jul. 11

452 431 21 21656.0 10568.0

Jul. 14

83 77 6 9136.9 7576.9

Jul. 18

129 119 10 13944.9 15102.3

Jul. 24

99 89 10 - -

-10000 0 10000 20000 30000 40000 50000 60000

15 -J un 17 -J un 19 -J un 22 -J un 26 -J un 29 -J un 3- Jul 11 -J ul 14 -J ul 18 -J ul 24 -J ul

N o. of e stim at ed popula tio n

-160 40 240 440 640 840 1040 1240

N o. of m ar ke d P. si ta niphonic a

No. of estimated population No. of marked butterflies

圖六、2000 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

Fig. 6. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2000.

(29)

19

表八、2001 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

Apr. 11 23 23 - - -

May 7 6 6 0 - -

May 14 20 20 0 - -

May 21 32 32 0 - -

May 28 81 81 0 - -

Jun. 4 172 172 0 - -

Jun. 5 230 228 2 20623.3 17604.4 Jun. 7 492 485 7 34580.6 20480.3 Jun. 11 124 123 1 43338.0 47521.1 Jun. 16 647 643 4 165322.0 101893.5 Jun. 18 902 881 21 20818.3 8861.2 Jun. 21 333 318 15 20099.9 11387.8

Jun. 26 94 81 13 - -

Jun. 29 18 17 1 - -

Jul. 9 52 49 3 - -

Jul. 13 6 5 1 - -

Jul. 18 2 1 1 - -

-50000 0 50000 100000 150000 200000 250000 300000

11 -A p r 7- M a y 14 -M a y 21 -M a y 28 -M a y 4- Jun 5- Jun 7- Jun 11 -J u n 16 -J u n 18 -J u n 21 -J u n 26 -J u n 29 -J u n 9- Jul 13 -J u l 18 -J u l

N o. of e st im at ed p op ul at io n

-175 25 225 425 625 825 1025

N o. o f m ar ke d P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖七、2001 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(30)

表九、2002 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

May 10

4 4 - - -

May 13

15 15 0 - -

May 16

56 56 0 - -

May 23

33 33 0 - -

May 27

303 303 0 - -

May 30

309 300 9 7666.1 4920.3

Jun. 6

86 83 3 - -

Jun. 10

40 37 3 1520.0 1440.5

Jun. 14

17 14 3 - -

Jun. 18

38 36 2 - -

Jul. 1

79 76 3 - -

0 2000 4000 6000 8000 10000 12000 14000

10 -M a y 13 -M a y 16 -M a y 23 -M a y 27 -M a y 30 -M a y 6- Jun 10 -J u n 14 -J u n 18 -J u n 1- Jul

N o. of e st im at ed po pu la ti on

0 50 100 150 200 250 300 350

No .o f m ark ed P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖八、2002 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(31)

21

表十、2003 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

May 13

20 20 - - -

May 20

76 76 0 - -

May 26

55 54 1 - -

Jun. 3

132 131 1 5319.6 6177.2

Jun. 5

577 566 11 3702.2 2499.6

Jun. 10

237 224 13 32906.5 35686.3

Jun. 19

267 263 4 29525.8 25770.2

Jun. 23

72 71 1 - -

Jun. 24

78 71 7 - -

Jul. 10

37 37 0 - -

-10000 0 10000 20000 30000 40000 50000 60000 70000 80000

13 -M a y 20 -M a y 26 -M a y 3- Jun 5- Jun 10 -J u n 19 -J u n 23 -J u n 24 -J u n 10 -J u l

N o. of e st im at ed p op ul at io n

-55 45 145 245 345 445 545 645

N o. o f m ar ke d P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖九、2003 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(32)

表十一、2004 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

May 10

7 7 - - -

May 19

106 106 0 - -

May 28

186 185 1 13420.6 14891.5

Jun. 2

265 257 8 3305.4 1528.5

Jun. 7

520 502 18 10188.1 3550.6

Jun. 12

186 164 22 2959.1 1040.4

Jun. 14

555 522 33 13723.6 7264.1

Jun. 21

75 70 5 3575.0 2723.3

Jun. 25

112 101 11 2168.7 2578.4

Jun. 28

33 30 3 - -

-5000 0 5000 10000 15000 20000 25000 30000

10 -M a y 19 -M a y 28 -M a y 2- Jun 7- Jun 12 -J u n 14 -J u n 21 -J u n 25 -J u n 28 -J u n

N o. o f es ti m at ed p op ul at io n

-100 0 100 200 300 400 500 600

N o. o f m ar ke d P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖十、2004 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(33)

23

表十二、2005 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

Apr. 18

5 5 - - -

May 20

45 45 0 - -

May 22

150 150 0 - -

May 23

202 201 1 19149.6 20752.1

May 27

366 363 3 70485.5 47753.2

May 30

6 4 2 - -

Jun. 1

645 635 10 21123.8 8585.8

Jun. 3

265 253 12 13497.6 5131.5

Jun. 6

534 507 27 15374.5 4696.5

Jun. 10

479 460 19 17101.1 5383.6

Jun. 13

678 656 22 37501.3 11512.8

Jun. 15

425 406 19 27669.7 9216.8

Jun. 20

477 452 25 47909.9 18605.1

Jun. 24

62 60 2 39122.0 35292.3

Jun. 27

678 643 35 30380.4 10900.8

Jun. 30

462 437 25 22381.9 10535.0

Jul. 4

43 39 4 - -

Jul. 7

192 171 21 - -

Jul. 15

13 12 1 - -

-20000 0 20000 40000 60000 80000 100000 120000 140000

18- Apr 20- M ay 22- M ay 23- M ay 27- M ay 30- M ay 1- Jun 3- Jun 6- Jun 10- Jun 13- Jun 15- Jun 20- Jun 24- Jun 27- Jun 30- Jun 4- Jul 7- Jul 15- Jul

No . of e stim ate d pop ulatio n

-100 0 100 200 300 400 500 600 700

N o. of ma rke d P. si ta ni pho ni ca

No. of estimated population No. of marked butterflies

圖十一、2005 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(34)

表十三、2006 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

Apr. 19

7 7 - - -

May 5

23 23 0 - -

May 23

21 21 0 - -

May 26

492 492 0 - -

Jun. 3

108 97 11 1332.5 697.5

Jun. 4

25 24 1 2275.0 2500.0

Jun. 5

8 7 1 - -

Jun. 6

108 106 2 45900.0 46783.8

Jun. 12

688 676 12 9989.8 4982.8

Jun. 16

268 256 12 27425.3 16703.2

Jun. 19

460 442 18 56937.8 63171.9

Jun. 23

15 15 0 - -

Jun. 27

67 63 4 - -

Jun. 30

25 23 2 - -

-20000 0 20000 40000 60000 80000 100000 120000 140000

19 -A p r 5- M a y 23 -M a y 26 -M a y 3- Jun 4- Jun 5- Jun 6- Jun 12 -J u n 16 -J u n 19 -J u n 23 -J u n 27 -J u n 30 -J u n

N o. o f es ti m at ed p op ul at io n

-110 -10 90 190 290 390 490 590 690 790

N o. of m ar ke d P . s ita ni phoni ca

No. of estimated population No. of marked butterflies

圖十二、2006 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(35)

25

表十四、2007 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

Apr. 30

9 9 - - -

May 17

17 17 0 - -

May 22

111 110 1 111.0 0.0

May 24

85 80 5 6885.0 7975.4

May 28

116 113 3 1754.5 1441.1

Jun. 1

211 201 10 635.1 462.2

Jun. 11

208 204 4 8458.7 5799.8

Jun. 14

420 404 16 420.0 0.0

Jun. 18

140 137 3 - -

-2000 0 2000 4000 6000 8000 10000 12000 14000 16000

30 -A p r 17 -M a y 22 -M a y 24 -M a y 28 -M a y 1- Jun 11 -J u n 14 -J u n 18 -J u n

N o. o f es ti m at ed p op ul at io n

-45 5 55 105 155 205 255 305 355 405 455

N o. of m ar ke d P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖十三、2007 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(36)

表十五、2008 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。

Date Captured Marked Recaptured Population size Eatimated S.E.

May 22

136 136 - - -

May 27

259 258 1 3296.4 4335.4

May 30

342 332 10 11696.4 5652.4

Jun. 2

165 157 8 14736.6 9421.6

Jun. 6

626 614 12 25264.7 12490.4

Jun. 9

525 516 9 38834.1 19890.5

Jun. 13

332 326 6 38413.6 21686.6

Jun. 16

423 413 10 105242.4 72930.0

Jun. 19

309 304 5 43251.2 27118.5

Jun. 23

244 236 8 36234.0 25868.7

Jun. 26

127 122 5 8183.9 5468.8

Jun. 30

124 115 9 4877.3 3488.4

Jul. 4

122 113 9 4717.3 5489.8

Jul. 7

62 58 4 - -

-20000 0 20000 40000 60000 80000 100000 120000 140000 160000 180000 200000

22 -M a y 27 -M a y 30 -M a y 2- Jun 6- Jun 9- Jun 13 -J u n 16 -J u n 19 -J u n 23 -J u n 26 -J u n 30 -J u n 4- Jul 7- Jul

N o. of e st im at ed p op ul at io n

0 100 200 300 400 500 600 700

N o. o f m ar ke d P . si ta ni ph on ic a

No. of estimated population No. of marked butterflies

圖十四、2008 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。

(37)

27

圖十五為 2000 年至 2008 年間，各年度最高族群估值。各年間最高估值仍舊相差甚大，但大約間隔三至四年，如 2001 年、2005 年、2008 年的青斑蝶標識中便出現了較多量的族群估值。

-50000 0 50000 100000 150000 200000 250000 300000

2000 yr. 2001 yr. 2002 yr. 2003 yr. 2004 yr. 2005 yr. 2006 yr. 2007 yr. 2008 yr.

N o. of hi ghe st e st im at ed popul at ion

圖十五、2000 年-2008 年間，陽明山大屯山區青斑蝶族群最高估值。

Fig. 15. From 2000 to 2008, the highest estimated population numbers of P. sita niphonica in Yangmingshan National Park.

圖十六至圖二十三為 2000 年至 2007 年間，就大屯山區六月間青斑蝶的族群估值，分別與鄰近測候站所紀錄之溫度、雨量、濕度等氣候因子資料進行比較。

由於大屯山測候站 (站號

C0AC40

) 為 2006 年方始架設，故以大屯山底部的鞍部測候站(站號

466910

) 為資料比對的來源。而比較的結果顯示，三項氣候因子對大屯山區青斑蝶的族群數量的變動並無明顯相關性，顯示大屯山區的氣候變化可能不是影響青斑蝶族群變動的主因。。

(38)

圖十六、2000 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

a. 溫度；b. 雨量；c. 溼度。

Fig. 16. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2000. a: temperature; b: rainfall; c: humidity.

2000yr.

0 10000 20000 30000 40000 50000 60000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p op ul at ion

15 16 17 18 19 20 21 22 23 24 25

T emp er at u re (cen ti gr ad e)

No. of estimated popilation Temperature

2000yr.

0 10000 20000 30000 40000 50000 60000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o . o f es ti m at ed p op ul at io n

0 10 20 30 40 50 60 70 80 90 100

R ai n fa ll (m m )

No. of estimated popilation Rainfall

2000yr.

0 10000 20000 30000 40000 50000 60000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o . o f e st im at ed p o lu la tio n

60 65 70 75 80 85 90 95 100

H u m id ity ( % )

No. of estimated popilation Humidity

b

c

a

(39)

29

圖十七、2001 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2001yr.

-50000 0 50000 100000 150000 200000 250000 300000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p opu la tio n

15 16 17 18 19 20 21 22 23 24 25

T em per at ur e (cen tig rad e)

No. of estimated popilation Temperature

2001yr.

-50000 0 50000 100000 150000 200000 250000 300000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st in at ed pop ul at ion

-10 0 10 20 30 40 50 60

R ai nf al l (mm)

No. of estimated popilation Rainfall

2001yr.

-50000 0 50000 100000 150000 200000 250000 300000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p op ul at io n

60 65 70 75 80 85 90 95 100

H u m id ity ( % )

No. of estimated popilation Humidity

b

c

a

(40)

圖十八、2002 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2002yr.

0 500 1000 1500 2000 2500 3000 3500

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p o pu lat io n

15 16 17 18 19 20 21 22 23 24 25

T emp er at ur e (cen tig rad e)

No. of estimated popilation Temperature

2002yr.

0 500 1000 1500 2000 2500 3000 3500

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p opu la tio n

0 10 20 30 40 50 60 70 80

Ra in fa ll (m m )

No. of estimated popilation Rainfall

2002yr.

0 500 1000 1500 2000 2500 3000 3500

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. o f e stim ate d p o pu la tio n

60 65 70 75 80 85 90 95 100

H u m id ity ( % )

No. of estimated popilation Humidity

b

c

a

(41)

31

圖十九、2003 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2003yr.

-10000 0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p opu la tio n

15 16 17 18 19 20 21 22 23 24 25

T em per at ur e (cen tig rad e)

No. of estimated popilation Temperature

2003yr.

-10000 0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p opu la tio n

-16 4 24 44 64 84 104 124

Ra in fa ll (m m )

No. of estimated popilation Rainfall

2003yr.

-10000 0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p op ul at io n

60 65 70 75 80 85 90 95 100

H um id ity ( % )

No. of estimated popilation Humidity

b

c

a

(42)

圖二十、2004 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2004yr.

-2000 0 2000 4000 6000 8000 10000 12000 14000 16000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o . o f es tim at ed p op ul at io n

15 16 17 18 19 20 21 22 23 24 25

T emp er at ur e (cen ti gr ad e)

No. of estimated popilation Temperature

2004yr.

-2000 0 2000 4000 6000 8000 10000 12000 14000 16000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o o f es tim ate d p o pu la ti on

-3.5 1.5 6.5 11.5 16.5 21.5 26.5

R ain fa ll ( m m )

No. of estimated popilation Rainfall

2004yr.

-2000 0 2000 4000 6000 8000 10000 12000 14000 16000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o . o f e stim ate d p op ula tio n

60 65 70 75 80 85 90 95 100

H u m id ity ( % )

No. of estimated popilation Humidity

b

c

a

(43)

33

圖二十一、2005 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2005yr.

0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. o f es ti m at ed p op ul at io n

15 16 17 18 19 20 21 22 23 24 25

T emp er at ur e (cen ti gr ad e)

No. of estimated popilation Temperature

2005yr.

0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p op ul at io n

0 20 40 60 80 100 120 140

Ra in fa ll (m m )

No. of estimated popilation Rainfall

2005yr.

0 10000 20000 30000 40000 50000 60000 70000 80000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p opu la tion

60 65 70 75 80 85 90 95 100

Hu m d ity ( % )

No. of estimated popilation Humidity

b

c

a

(44)

圖二十二、2006 年六月份，陽明山大屯山區青斑蝶族群估值與氣候資料的比較。

2006yr.

-20000 0 20000 40000 60000 80000 100000 120000 140000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o. of e st im at ed p op ul at io n

15 16 17 18 19 20 21 22 23 24 25

T emp er at ut re ( cen ti gr ad e)

No. of estimated popilation Temperature

2006yr.

-20000 0 20000 40000 60000 80000 100000 120000 140000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

N o . of e st im at ed p op ul at io n

-23 -3 17 37 57 77 97 117 137 157

R ai nf al l (mm)

No. of estimated popilation Rainfall

2006年yr.

-20000 0 20000 40000 60000 80000 100000 120000 140000

1- Jun

2- Jun

3- Jun

4- Jun

5- Jun

6- Jun

7- Jun

8- Jun

9- Jun

10- Jun

11- Jun

12- Jun

13- Jun

14- Jun

15- Jun

16- Jun

17- Jun

18- Jun

19- Jun

20- Jun

21- Jun

22- Jun

23- Jun

24- Jun

25- Jun

26- Jun

27- Jun

28- Jun

29- Jun

30- Jun

台灣地區青斑蝶族群變動與移動行為之探討

國立台灣大學生物資源暨農學院昆蟲學研究所 碩士論文

Graduate Institute of Entomology College of Bioresources and Agriculture

National Taiwan University Master Thesis

台灣地區青斑蝶族群變動與移動行為之探討 Study on the population dynamics and moving behavior

of the chestnut tiger butterfly (Parantica sita niphonica Moore) in Taiwan

李信德 Hsin-Te Lee

指導教授：楊平世 博士

Advisor: Ping-Shih Yang, Ph.D.

中華民國 98 年 6 月

June, 2009

ii

iv

sita niphonica Moore) in Taiwan

vi

viii

1

壹、緒言

貳、往昔研究

3

5

参、材料與方法

7

Coordinate (TW67)

9

a b

11

N M m

r

M a Z r

R

r α

n

N M

α

M n a Z

N n

α R r

( )

(

)

M r a

N N N n

M R a r

α

466910 鞍部 臺北市 825.8 302428 2786547

466920 臺北 臺北市 5.3 301131 2770270

466930 竹子湖 臺北市 607.1 304059 2784169

C0A540 四堵 臺北縣 401 324508 2754235

C0A530 坪林 臺北縣 300 320834 2759293

C0A570 桶後 臺北縣 387 309796 2749302

C0A9B0 石牌 臺北市 7 300982 2778958

C0A9C0 天母 臺北市 63 303334 2779121

C0A9F0 內湖 臺北市 20 307050 2774890

C0A9G0 南港 臺北市 42 310004 2772287

C0A9H0 信義 臺北市 51 306400 2769318

C0AC40 大屯山 臺北市 1098 301825 2785515

13

肆、結果

Year Numbers of total marked

Male

Female

2000

2001

2002

2003

2004

2005

2006

2007

2008

Total

15

Year Numbers of total marked

Male

Female

2000

2001

2002

2003

國立台灣大學生物資源暨農學院昆蟲學研究所碩士論文

指導教授：楊平世博士

466910 鞍部臺北市 825.8 302428 2786547

466920 臺北臺北市 5.3 301131 2770270

466930 竹子湖臺北市 607.1 304059 2784169

C0A540 四堵臺北縣 401 324508 2754235

C0A530 坪林臺北縣 300 320834 2759293

C0A570 桶後臺北縣 387 309796 2749302

C0A9B0 石牌臺北市 7 300982 2778958

C0A9C0 天母臺北市 63 303334 2779121

C0A9F0 內湖臺北市 20 307050 2774890

C0A9G0 南港臺北市 42 310004 2772287

C0A9H0 信義臺北市 51 306400 2769318

C0AC40 大屯山臺北市 1098 301825 2785515