國立台灣大學生物資源暨農學院昆蟲學研究所 碩士論文
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
誌 謝
天
中文摘要
本研究以標識再捕法分析台灣地區青斑蝶族群的組成與周年動向,並就陽明 山大屯山區的青斑蝶再捕獲資料進行族群數量估算。2000-2008年間,大屯山區青 斑蝶族群數量最低為418隻,最高為165322隻,年度間數量變化差異極大。就長時 間的連續變化而言,此區族群的數量變動可能以3-4年為一週期。就異地再捕獲的 長距離移動與各地區的出現時間紀錄,青斑蝶有季節性、固定方向的移動行為,
但短距離的移動則似無方向規律。冬季時無群聚越冬現象,雌蝶亦無生理滯育,
於多處地區可觀察到成蝶活動與幼蟲的世代紀錄,其越冬方式可能為多策略並行。
關鍵詞:青斑蝶、族群變動、移動行為、標識再捕法、族群估算、陽明山國家公 園
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
目 錄
口試委員會審定書...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
vi
(二) 觀音樣區青斑蝶周年分布...15
二、族群估算...16
三、再捕獲移動紀錄...36
(一) 島外移動記錄...36
(二) 島內移動記錄...36
伍、討論...42
一、台灣地區青斑蝶族群組成與季節動向...42
(一) 陽明山地區的青斑蝶族群...42
(二) 南澳地區的青斑蝶族群...43
二、台灣地區青斑蝶的越冬狀態...44
三、族群估算與影響因素...45
四、台灣地區青斑蝶移動與遷移的可能性...46
(一) 標識再捕紀錄的方向性...47
(二) 遷移的生理條件...48
五、青斑蝶標識研究與推廣教育的關連...49
陸、結論...51
柒、參考文獻...52
捌、附錄...59
圖 次
圖一、陽明山國家公園大屯山研究樣區...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
viii
圖二十七、南向再捕獲紀錄 1...37
圖二十八、南向再捕獲紀錄 3...37
圖二十九、南向再捕獲紀錄 4 與 5...38
圖三十、南向再捕獲紀錄 6 與 7...38
圖三十一、島內再捕獲記錄 1...38
圖三十二、島內再捕獲記錄 2...39
圖三十三、島內再捕獲記錄 3...39
圖三十四、16 方位角示意圖...39
表 次
表一、南澳鄉觀音產業道路標定寄主植物之 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
1
壹、緒言
青斑蝶 (Parantica sita niphonica Moore) 屬於蛺蝶科斑蝶亞科,為台灣地區常 見種類的蝴蝶,廣泛分布於全島各地。因其食草為蘿藦科 (Asclepiadaceae) 的有毒 植物,故蟲體內具有毒性作為防禦之用。所以成蝶天敵較少,亦較不畏懼人為干 擾,是適合的賞蝶對象。
在台灣北部陽明山國家公園境內,每年春夏季節均有大量的青斑蝶出現,尤 以大屯山區的大發生景象最為壯觀。大發生期間蝴蝶數量在短時間內極劇的增 加,而期間過後卻又急速減少,2000 年曾首次証實此區的蝴蝶有跨海長距離飛行 至日本的行為 (Lee and Yang,2000),其後陸續的再捕獲紀錄也證實台日之間青斑 蝶的族群交流。
而以往對於此種蝴蝶的研究多屬於描述其生物學的範疇,族群的數量變化與 青斑蝶的飛行移動資料較少,雖有學者曾針對族群數量研究,但缺乏長期且持續 的觀察,對欲了解台灣地區青斑蝶族群的變動資料仍嫌不足。本研究以標識再捕 法為主,對台灣地區的青斑蝶族群進行調查,希望對青斑蝶的族群變動建立基本 資料,並對青斑蝶在台灣島嶼內及島嶼間的季節動向與飛行路徑能有初步了解。
貳、往昔研究
一、青斑蝶概述
青斑蝶 (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 年,由九州鹿兒島縣的鹿兒島昆蟲同好會所發起;隔年五月,日 本關西地區的三重縣再捕獲由九州南方的種子島所標記的青斑蝶,此為首例青斑
3
蝶北上飛行的再捕獲紀錄;而首例南下飛行再捕獲紀錄則為同年十一月,於九州 南方的奄美大島再捕獲由鹿兒島縣內標記的青斑蝶個體。
日本青斑蝶標放研究中主要的研究團體,為大阪地區的大阪青斑蝶調查會,
以及九州鹿兒島縣青斑蝶大移動研究會分別成立於 1983 與 1984 年;大阪青斑蝶 調查會並於 1990 年更名為青斑蝶調查會。1994 年起,利用網際網路架設網站於留 言版系統,加速青斑蝶標識資訊的流通與擴展參與標識研究的人員階層。當時最 接近台灣地區的再捕獲紀錄,為 1995 年沖繩縣與那國島所再捕獲,由關西大阪市 境內所邊放的青斑蝶,與台灣本島的宜蘭縣最近距離僅約 110km。
日本地區青斑蝶的移動現象可分為三個時期:
(一)、春季北上
春季到夏初 (三~七月) 西南季風盛行,由越冬地區新羽化的青斑蝶世代,
可能利用西南氣流由九州南方各島開始往日本北方移動,但因日本大部分地區 中,春季時非青斑蝶的主要活動季節,因此標識再捕的紀錄較少。
(二)、夏季山區間的移動
七月~八月下旬是日本各地青斑蝶主要活動時間,因為幼蟲主要的食草與 成蝶的蜜源植物生長於山區,因此成蝶出現地點也以山區或高地為主;一般認 為,春季北上的個體向北遷移後也會往山區移動。夏季在各山區標識再捕的結 果,鮮少有長距離移動的情形,成蝶僅在數公里之內的山系間移動。
(三)、秋季南下
夏末秋初,日本各地青斑蝶均有明顯移動現象。小地區範圍內,青斑蝶由 山區向平地移動;廣域地區範圍,本島的青斑蝶族群順著當時季風風向,往南 方朝九州方向飛行。南向飛行的地點不僅止於九州本島,而是沿著東亞島弧向 南方前進,沿薩南諸島、奄美諸島、琉球群島,最遠可到達八重山群島中的与 那國島,離台灣本島極近。因此,亦推論台灣可能為青斑蝶長距離長距離的停
留點。
遷移行為發生的可能原因,包括短日照、低溫、季節氣流方向、寄主植物老 化枯死,以及蜜源植物減少等環境因子,日方研究者並就台灣與日本的地理位置,
推測台灣地區的青斑蝶族群亦有遷移行為 (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),亦提到秋冬季節青斑蝶群聚於南部六龜地區。
但在南部諸多紫斑蝶類與青斑蝶類越冬的谷地,青斑蝶均無出現在越冬集團
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),於台灣北海岸觀察到斑蝶成群移動現象,
但其中亦並無青斑蝶紀錄。
参、材料與方法
一、研究地點
(一)、陽明山國家公園大屯山區
調查地點為大屯主峰登山車道沿線。
據 Wei (1995),大屯主峰登山車道屬於陡升車道地形,海拔高度由 872 公尺上升至 1092 公尺,車道寬度約 4 公尺。前半段為紅楠與台灣矢竹草原植 物社會,後段為台灣芒草與台灣矢竹草原植物社會。本段路面植被遮蔽度少,
車道旁大量孳生斑蝶類主要的菊科蜜源植物-島田氏澤蘭。車道兩旁每年均 有固定進行除草,除草影響蜜源植物生長,而近年除草的時間已延後至斑蝶 類大發生過後進行。
圖一、陽明山國家公園大屯山研究樣區。
Fig. 1. Study area in Datung Mt., Yangmingshan National Park.
7
自 2000 年~2008 年,視各年青斑蝶發生狀況而略有差異,於四月下旬
~五月上旬,青斑蝶開始出現起,每週兩次,每次間隔二~三日,於大屯山 登山車道沿線進行青斑蝶成蝶標識。標識進行持續到青斑蝶族群數量明顯減 少為止,約持續進行至七月中旬。其它月份則每月進行一次調查。
(二)、南澳觀音產業道路
調查區域位於省道台九線蘇花公路 144.5K 處,據楊與李(1995),為觀音 海岸野生動物重要棲息環境範圍內。此區地形南北狹長,東向陡坡急降入海,
溼度高,全年多雨。此區林相主要為台灣赤楊、大葉楠、茄冬、樟科與榕屬 的植物組成;在裸露或陡坡面則為五節芒的優勢族群。
於調查全島青斑蝶分布狀況時,發現本區有大量的青斑蝶寄主植物台灣 牛嬭菜 (Marsdenia formosana),且多位於道路兩旁的林班地內,少遭砍伐,
觀察容易,有穩定的成蟲與幼蟲數量,故於此以 GPS 標定 16 株台灣牛嬭菜,,
如表一所列。並於 2005 年十一月至 2007 年九月止,每周一次,調查並記錄 標定寄主植物上的青斑蝶幼生期,包括卵、幼蟲、蛹的數量與寄生狀況。
成蝶的調查時間則由 2005 年十一月至 2007 年九月止。除了此區,亦同 時將東澳粉鳥林、烏石鼻自然保留區、南澳金岳村等地的成蝶觀察紀錄一並 列入。
(三) 其他青斑蝶活動環境
於 2000-2008 年間,每月一次,於島內青斑蝶主要的分布環境,調查青斑 蝶成蝶數量並進行標記。
表一、南澳鄉觀音產業道路標定寄主植物之 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)、所有個體被捕獲的機率相同,與它們在棲地中的分布位置無關。
9
4、取樣須在分離的時間段進行,且標識進行的時間與研究進行時間相較 之下須佔極少比例。
標識方式以油性黑色細字簽字筆,於青斑蝶後翅中室書寫標記代碼後原地 釋放。標識所使用的標記代碼由四部份組成,如圖二所示:
1、編號:第一隻青斑蝶標識為 001,第二隻標識為 002,之後依序類推。
2、標識地點:以標識地的英文縮寫表示,如陽明山國家公園 (Yang-Ming-Shan National Park)縮寫表示為 YMS。
3、標識單位:以標識者所屬單位的英文縮寫表示,如台灣大學 (National Taiwan University) 的縮寫表示為 NTU。
4、標識日期:日期以四位代碼表示,如 4 月 13 日則記為 0413。但在陽 明山地區,因大發生期間青斑蝶數量過多,為節省時間,首次標識日,
日期代碼改註記為「A」,第二次標識日註記代碼為「B」之後依序類 推,以圖一為例,D表該年度第四次的標識日。
圖二、青斑蝶標識代號說明 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
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)。
封閉族群的估算假設前提為族群在研究進行當中不會有任何變動,任何個 體均無法離開研究樣區或由樣區外進入;而開放族群則會因研究期間的出生、
11
死亡、遷入、遷出等因素而產生變動。
而常用於開放式族群的估算有三種方法:
1、Fisher and Ford trellis method 2、Bailey`s triple catch
3、Jolly-Seber stochastic method
其中 Jolly-Seber stochastic method 為最常使用的方法,基本公式為:
i i i
i
N M m
=
r
Ni 為 i 日的族群估值。
Mi 為 i 日標記總數量的估值。
ri 為 i 日再捕獲標記個體的總數量。
ni 為 i 日捕獲的總數量。
而 Mi 的估算公式為: i i i i
i
M a Z r
=
R
+ ai 為 i 日所標記的個體數。Ri 為 i 日所標記,而於其後標示日被再捕獲的總數。
Zi 為在 i 日前標記,但不在 i 日被捕獲,而在之後被再捕獲的總數。
而 i 日中,標記個體在族群裡的比例為 i i
i
r α
=n
所以原來的估算公式可表示為: i ii
N M
=
α
Mi 帶入後估算公式變成為: i i i i i i
i i i
M n a Z
N n
α R r
= = +
其中,var
( )
i i(
i i)
i i i 1 1 1 ii i i i
M r a
N N N n
M R a r
⎧ − + ⎡ ⎤ −
α
⎫⎪ ⎪
= − ⎨ ⎢ − ⎥+ ⎬
⎪ ⎣ ⎦ ⎪
⎩ ⎭
楊 (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
13
圖四、氣象測候站位置示意圖。
Fig. 4. Locations of meteorological stations.
肆、結果
一、青斑蝶個體資料
(一) 青斑蝶族群組成
表三為 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%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 次調查,卵與幼蟲各齡期紀錄列於
附錄三,其中 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 就陽明山大屯山區 進行的標識再捕法所得到的再捕獲資料,進行青斑蝶族群的絕對豐度估算結果。
各年內與各年度間估算數量差異頗大。圖六至圖十四為各年標記數量與族群估值 的比較圖,各年度內,標記數量與該次族群估值並無明顯相關。
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
表七、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.0Jun. 19
1050 1042 8 37520.0 17437.7Jun. 22
569 550 19 26648.2 8920.4Jun. 26
744 712 32 36216.9 12948.1Jun. 29
466 446 20 16053.7 6031.3Jul. 3
48 41 7 2367.4 1425.1Jul. 11
452 431 21 21656.0 10568.0Jul. 14
83 77 6 9136.9 7576.9Jul. 18
129 119 10 13944.9 15102.3Jul. 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.
19
表八、2001 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 8. Population estimation of P. sita niphonica in Yangmingshan National Park, 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 7. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2001.
表九、2002 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 9. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.3Jun. 6
86 83 3 - -Jun. 10
40 37 3 1520.0 1440.5Jun. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 8. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2002.
21
表十、2003 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 10. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.2Jun. 5
577 566 11 3702.2 2499.6Jun. 10
237 224 13 32906.5 35686.3Jun. 19
267 263 4 29525.8 25770.2Jun. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 9. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2003.
表十一、2004 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 11. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.5Jun. 2
265 257 8 3305.4 1528.5Jun. 7
520 502 18 10188.1 3550.6Jun. 12
186 164 22 2959.1 1040.4Jun. 14
555 522 33 13723.6 7264.1Jun. 21
75 70 5 3575.0 2723.3Jun. 25
112 101 11 2168.7 2578.4Jun. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 10. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2004.
23
表十二、2005 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 12. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.1May 27
366 363 3 70485.5 47753.2May 30
6 4 2 - -Jun. 1
645 635 10 21123.8 8585.8Jun. 3
265 253 12 13497.6 5131.5Jun. 6
534 507 27 15374.5 4696.5Jun. 10
479 460 19 17101.1 5383.6Jun. 13
678 656 22 37501.3 11512.8Jun. 15
425 406 19 27669.7 9216.8Jun. 20
477 452 25 47909.9 18605.1Jun. 24
62 60 2 39122.0 35292.3Jun. 27
678 643 35 30380.4 10900.8Jun. 30
462 437 25 22381.9 10535.0Jul. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 11. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2005.
表十三、2006 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 13. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.5Jun. 4
25 24 1 2275.0 2500.0Jun. 5
8 7 1 - -Jun. 6
108 106 2 45900.0 46783.8Jun. 12
688 676 12 9989.8 4982.8Jun. 16
268 256 12 27425.3 16703.2Jun. 19
460 442 18 56937.8 63171.9Jun. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 12. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2006.
25
表十四、2007 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 14. Population estimation of P. sita niphonica in Yangmingshan National Park, 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.0May 24
85 80 5 6885.0 7975.4May 28
116 113 3 1754.5 1441.1Jun. 1
211 201 10 635.1 462.2Jun. 11
208 204 4 8458.7 5799.8Jun. 14
420 404 16 420.0 0.0Jun. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 13. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2007.
表十五、2008 年陽明山大屯山區青斑蝶標識再捕與族群估算結果。
Table 15. Population estimation of P. sita niphonica in Yangmingshan National Park, 2008.
Date Captured Marked Recaptured Population size Eatimated S.E.
May 22
136 136 - - -May 27
259 258 1 3296.4 4335.4May 30
342 332 10 11696.4 5652.4Jun. 2
165 157 8 14736.6 9421.6Jun. 6
626 614 12 25264.7 12490.4Jun. 9
525 516 9 38834.1 19890.5Jun. 13
332 326 6 38413.6 21686.6Jun. 16
423 413 10 105242.4 72930.0Jun. 19
309 304 5 43251.2 27118.5Jun. 23
244 236 8 36234.0 25868.7Jun. 26
127 122 5 8183.9 5468.8Jun. 30
124 115 9 4877.3 3488.4Jul. 4
122 113 9 4717.3 5489.8Jul. 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 年陽明山大屯山區青斑蝶族群估值與標記數量的比較。
Fig. 14. Compare estimated population numbers with marked numbers of P. sita niphonica in Yangmingshan National Park, 2008.
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
) 為資料比對的來源。而比較的結果顯示,三項氣候因子對 大屯山區青斑蝶的族群數量的變動並無明顯相關性,顯示大屯山區的氣候變化可 能不是影響青斑蝶族群變動的主因。。圖十六、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
29
圖十七、2001 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 17. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2001. a: temperature; b: rainfall; c: humidity.
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
圖十八、2002 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 18. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2002. a: temperature; b: rainfall; c: humidity.
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
31
圖十九、2003 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 19. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2003. a: temperature; b: rainfall; c: humidity.
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
圖二十、2004 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 20. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2004. a: temperature; b: rainfall; c: humidity.
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
33
圖二十一、2005 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 21. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2005. a: temperature; b: rainfall; c: humidity.
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
圖二十二、2006 年六月份,陽明山大屯山區青斑蝶族群估值與氣候資料的比較。
a. 溫度;b. 雨量;c. 溼度。
Fig. 22. Compare estimated population numbers with climate factors during June in Yangmingshan National Park, 2006. a: temperature; b: rainfall; c: humidity.
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