1
國立臺灣大學生命科學院漁業科學研究所 碩士論文
Institute of Fisheries Science, College of life Science National Taiwan University
Master Thesis
石門水庫集水區湳仔溝溪之臺灣馬口魚族群生態學研究
The population ecology of lake candidus dace (Candidia barbatus) in Nanzaigou stream of Shihmen Reservoir in Taiwan
研究生:楊芷毓 撰 Chih-Yu Yang
指導教授:李英周 博士 朱達仁 博士 Advisor: Ying-Chou Lee, Ph.D
Ta-Jen Chu, Ph.D
中華民國 100 年 1 月
January, 2011
ii
謝辭
本研究得以完成,首先要感謝指導教授 李英周老師在我修讀碩士期間在課 業、研究或是論文寫作上細心的教導與指正,及感謝共同指導教授 朱達仁老師 在湳仔溝溪資料上的全力協助與研究論文上的悉心指導,在此致上最崇高的敬意 與謝意。此外也感謝交通大學土木工程系 郭一羽老師與臺灣大學生命科學系 莊 鈴川老師,百忙之中給予論文的指正與建議。
實驗過程中,很幸運的有許多默默幫助我的人,在此一一致謝。感謝桃園縣 百吉國小志工隊及中華大學朱達仁老師實驗室的同仁們在湳仔溝溪採樣上的大 力幫助,因為你們我的論文才可以順利的進行。也感謝實驗室惟哲學長和君翰學 長當我在論文與實驗上遇到瓶頸時,熱心的幫助並給予建議,還有秉宏、承容、
榮琦、小不點學長們跟長毛在碩士過程中帶給我許多歡樂和度過適應期,阿達、
正勳、小馬、小號、孟祺、宗祐和小暐的陪伴,還有我的同學球球和小 P,因為 我們的互相打氣、一同奮鬥、歡笑、忙碌、鬱悶和互吐苦水才可以順利的完成學 業。同時也感謝大學和高中時期好友們的關心和鼓勵。
最後,僅將此文獻給摯愛的爺爺、奶奶和爸爸與乾媽一家人,感謝你們一路 上的支持和無悔的付出,我才能無後顧之憂的好好完成學業,謝謝你們。
iii
摘要
石門水庫上游集水區湳仔溝溪多年來受到水患影響,於2004年3月至2006年 11月進行整治工程,本研究針對該溪之優勢魚種臺灣馬口魚,期望透過長時間調 查湳仔溝溪之臺灣馬口魚族群探討其密度、體長、體重、肥滿度受到溪流整治之 影響及年齡與成長和生殖週期之研究,了解湳仔溝溪的臺灣馬口魚在受到溪流整 治的人為因子及環境變動下族群的結構與動態是否受到衝擊,並探討臺灣馬口魚 族群在湳仔溝溪的年齡組成和生殖週期等特性,以期更充分理解臺灣馬口魚在湳 仔溝溪之族群生態。2003年11月至2010年11月,依上游百吉四號橋、湳仔溝二號 橋、龍行橋、鳳儀僑、新福圳至下游懷德橋各測點由下游往上游以電魚法進行魚 類採樣,以one-way ANOVA比較(1)密度、(2)體長、(3)體重、(4)肥 滿度之平均值有無差異,試圖了解整治工程對臺灣馬口魚之族群有無影響;2009 年9月至2010年11月,每月以蝦籠法與電魚法在相同之六測點採樣並進行年齡形 質及與生殖週期之研究。結果顯示,溪流整治對臺灣馬口魚的影響中,密度平均 值以龍行橋施工期間最大(7.00尾/15分鐘/100公尺),但龍行橋在施工期間馬口 魚的體長及體重偏小,推測原因為龍行橋有加設「避魚柵」及「生態島」對優勢 魚種之馬口魚幼魚提供了庇護之故;在第一期施工期間受到艾莉風災之影響,上 游各測點之馬口魚密度下降,溪流整治在施工期發生颱風對施工測點馬口魚密度 影響極大,而鳳儀橋在施工期間馬口魚體長體重平均值為最大,鳳儀橋其施工工 法可能對馬口魚之成魚有助益,因此,湳仔溝溪進行溪流整治的過程中,對於湳 仔溝溪臺灣馬口魚的密度、體長、體重及肥滿度未造成極大的傷害,但還是需要 後續更長期的環境調查資料探討在施工過後湳仔溝溪台灣馬口魚族群的變動情 況;以鱗片形質來進行臺灣馬口魚年齡與成長之研究中,共捕獲192隻,體長範 圍在2.9-17.0cm,判讀輪數最高為三輪,利用邊緣成長率分析輪紋一年形成一輪,
輪紋形成期約在7-9月,與生殖季同時,推測可能為生殖輪。年齡組成分別為0歲
iv
魚23.4%,1歲魚24.7%,2歲魚41.8%及3歲魚10.1%。臺灣馬口魚族群在湳仔溝溪 面臨艾利風災與溪流整治使族群的密度降低、體重與肥滿度偏小、以及物種間食 物的競爭壓力大,使臺灣馬口魚傾向選擇發育快、體型小及壽命較短的生存策略,
導致年齡結構以提早成熟的低齡魚為主。未來對於臺灣馬口魚保育和復育湳仔溝 溪臺灣馬口魚族群,重點應該放在湳仔溝溪的環境改善,保持良好的水質與適合 棲息的棲地環境,將有助於族群恢復。
關鍵字:臺灣馬口魚、族群生態、湳仔溝溪、鱗片、年齡與成長
v
Abstract
Fish will change with and increases along with the environmental condition, the community fades and the age structure. Research in fish’s age structure, sex ratio, density, the sphere of action as well as the environment factor influence are called the ecology population. Nanzaigou stream was located at the tributaries of Shihmen Reservoir, it was often over-flooded and caused serious damage along the stream side, so it was restored by the natural ecological engineering from March 2004 to November 2006. The study intends to explore the differences of the dominant species Zacco barbata population before, during and after restoration at each site in Nanzaigou Stream, as well as the age structure and reproductive cycle of this species. Sampling was conducted along the stream from Baiji NO.4 Bridge, Nanzihgou NO.2 Bridge, Longxing Bridge, Fengyi Bridge, Sinfu Canal and Huaide Bridge by electrofishing from November 2003 to November 2010. And the mean values of dominant species densities, body length, body weight and condition factor were compared by one-way ANOVA. And used scales to assessment the age and growth of the dominant species were capture of the same sites by Shrimp cage and electrofishing from September 2009 to November 2010. The results were as follows.
The density mean value in Long Xingqiao is biggest (7.00/15min), but the fork length and body weight were small, extrapolated that maybe the Evades the fish grid and Ecology island has provided reason of the asylum to the juvenile fish. The first phase of construction with the typhoon, the densities of fish each site declined, so has the typhoon to be enormous to the construction measuring point to the fish density influence. The mean value of body length and weight were biggest at Fengyi Bridge, extrapolated that maybe the construction labor law to the adult fish is of
vi
help.Comprehensive result, when Nanzaigou stream carries on the restoration there has not caused the enormous damage on the mean value of density, body length, body weight and condition factor for Zacco barbata, But it will need the following longer-term environment inquiry material discussion to the population of Zacco barbata after construction in Nanzaigou stream. A total of 192 individuals were captured to use scales to assessment the age and growth of the dominant species. The range of fork length was 2.9-17.0 cm, the highest rings were counted to 3. The rate of marginal increment was highest from July to September, and which indicated that one ring was formed yearly. The age composition respectively was 0 year old of fish 23.4%, 1 year old of fish 24.7%, 2 year old of fish 41.8% and 3 year old of fish 10.1%.
The population in Nanzaigou stream faced with typhoon caused disaster and ecological engineering causes population density to reduce, the body weight and the condition factor is small, as well as the species food competition pressure high, causes it to choose r-selection the survival strategy, the age structure to shift to an earlier time the mature low age fish primarily.In the future regarding the species’ care and the cicada chrysalis, which should place the environment improvement with emphasis, maintains the good water quality and suitable to perch roosts the environment, will be helpful in the population restores.
Keywords: Candidia barbatus, population ecology, Nanzaigou Stream, scales, Age and growth.
vii
目錄
口試委員審定書...i
謝辭...ii
中文摘要………iii
英文摘要………....v
第一章 前言………1
1.1.河川生態系之簡介………1
1.2.魚類族群生態學簡述………2
1.3.生態工法沿革與意涵………....4
1.4.實驗魚種簡介………5
1.5 研究緣起……….6
1.6.研究目的………8
第二章 材料與方法………9
2.1.研究地點………9
2.2.採樣方法………9
2.3.樣本魚處理……….10
2.4.鱗片處理………..10
2.5.肥滿度指數之估計………..10
viii
2.6.生殖腺指數之估計………..11
2.7.性比………..11
2.8.資料處理………..11
第三章 結果………..13
3.1.溪流整治對臺灣馬口魚之影響………..13
3.1.1.族群密度………..13
3.1.2.體長………..14
3.1.3.體重………..15
3.1.4.肥滿度………..15
3.1.5.環境因子………..15
3.2.年齡與成長………..16
3.2.1.研究樣本簡介………..16
3.2.2.尾叉長與體重之關係………...16
3.2.3.年齡查定鱗片之選定………..17
3.2.4.輪徑長與尾叉長之關係………..17
3.2.5.輪紋形成時期………..17
3.2.6.輪紋形成時期………..18
3.2.7.生殖腺指數的變化………..18
ix
3.2.8.肥滿度指數的變化………..19
3.2.9.性比………..19
3.2.10.年齡組成………19
第四章 討論………..20
4.1.溪流整治對臺灣馬口魚族群之影響………..20
4.2.年齡與成長………..23
4.3.未來展望………..27
第五章 結論………..28
參考文獻………..29
x
表目錄
表 1.國內外生態工法之沿革與定義……….39
表 2.湳仔溝溪各測點施工時間與施工項目……….40
表 3.河川棲地型態分類系統表……….41
表 4 湳仔溝溪臺灣馬口魚施工期間之密度變化………..42
表 5.各測點臺灣馬口魚之平均密度……….43
表 6.各測點臺灣馬口魚之平均體長……….44
表 7.各測點臺灣馬口魚之平均體重……….45
表 8.各測點臺灣馬口魚之平均肥滿度……….46
表 9.湳仔溝溪各測點水文因子平均值……….47
表 10.湳仔溝溪各測點底石組成平均值………...48
表 11.湳仔溝二號橋環境因子主成分分析結果的特性向量與特徵值…………...49
表 12.龍行橋環境因子主成分分析結果的特性向量與特徵值………...50
表 13.鳳儀橋環境因子主成分分析結果的特性向量與特徵值………...51
表 14.新福圳環境因子主成分分析結果的特性向量與特徵值………...52
表 15.未施工測點環境因子主成分分析結果的特性向量與特徵值………...53
表 16.鱗片形質判讀二次之判讀率………...54
表 17.各輪數之平均輪徑長………...55
xi
圖目錄
圖 1.湳仔溝溪地理位置圖……….56
圖 2.湳仔溝溪六測點地理位置圖……….57
圖 3.湳仔溝溪施工期間之魚種組成……….58
圖 4.湳仔溝溪各測點臺灣馬口魚之平均密度……….59
圖 5.臺灣馬口魚 2003-2010 平均數量變化圖……….60
圖 6.臺灣馬口魚 2003-2010 平均重量變化圖……….61
圖 7.湳仔溝溪各測點臺灣馬口魚之體長頻度……….62
圖 8.湳仔溝溪各測點臺灣馬口魚之平均體長……….63
圖 9.湳仔溝溪各測點臺灣馬口魚之體重頻度……….64
圖 10.湳仔溝溪各測點臺灣馬口魚之平均體重………...65
圖 11.湳仔溝溪各測點臺灣馬口魚之肥滿度頻度………...66
圖 12.湳仔溝溪各測點臺灣馬口魚之平均肥滿度………...67
圖 13.臺灣馬口魚之尾叉長與體重頻度圖………...68
圖 14.雄魚與雌魚之尾叉長頻度圖………...69
圖 15.臺灣馬口魚尾叉長與體重之關係………...70
圖 16.鱗片選取之部位之測量軸平均值與 c.v 值………...71
圖 17.臺灣馬口魚輪徑長與尾叉長之關係………...72
圖 18.臺灣馬口魚尾叉長與輪數頻度圖………...73
圖 19.臺灣馬口魚輪徑與輪數頻度圖………...74
圖 20.臺灣馬口魚輪徑與輪紋之關係………...75
圖 21.臺灣馬口魚鱗片季別邊緣成長率頻度圖………...76
圖 22.臺灣馬口魚季別生殖腺指數(GSI)變化圖……….77
圖 23.臺灣馬口魚尾叉長與生殖腺成熟狀態之關係圖………...78
xii
圖 24. 雌、雄魚 50%性成熟尾叉長圖………..79
圖 25.雌、雄魚季別平均肥滿度(CF)變化圖………...80
圖 26. 雌、雄魚體長頻度圖……….81
圖 27.湳仔溝溪臺灣馬口魚年齡分布圖………...82
圖 28.雌、雄魚年齡分布圖………...83
xiii
照片目錄
照片 1.鱗片選取之部位圖………84
照片 2.臺灣馬口魚鱗片之結構圖………85
照片 3-a.湳仔溝溪整治項目避魚柵………86
3-b.湳仔溝溪整治項目生態島………86
1
第一章、前言 1.1 河川生態系之簡介
河川環境(river environment)廣義上可分為(1)河川的外在環境:泛指河 道外的陸域環境;(2)河川的內在環境:以河川生物棲息範圍為主的水、陸域 環境(汪 1995)。而河川的環境會隨上、中、下游不同的地質、地形、坡度、
溫度、生物種類、汙染程度、人為侵入等因子變化,而使流水生態系統(lotic ecosystem)也隨之變化。生活在河川水域中的生物,是河川環境中有生命的組 成因子,並與周遭的非生物因子理化環境,共構了河川生態系(Hynes 1970)。
河川主要有三種地形區分:(1)侵蝕區:河床坡度大、河道不穩定、河岸陡峭 而狹窄、沖蝕極為嚴重,多分布在河川的上游源頭。(2)儲藏和運送區:河道 最穩定、輸送沉積物、養分和水、沖蝕與淤積平衡,通常是在侵蝕區的下游。(3)
沉積區:面積比侵蝕區與輸送區大,坡度較為和緩,淤積問題嚴重,通常位於河 口(Faber et al. 1989;陳等 1998)。因河川流經了此三種區域,形成了許多不 同的地形地貌,如峽谷、瀑布、水潭、沙洲、潟湖及濕地等(汪 1992;高 1993)。
而河川也可依據等級和長度來做分類,例如河川序級(stream order)就是「河川 在支流等及位置正式的測量法」,意即第一級河川沒有支流,第二級河川有第一 級河川作支流,第三級河川有第一級和第二級河川作支流,依此類推(Leopold et al. 1964)。Vannote et al.(1980)提出河川續動概念(river continuum concept;
RCC),其將自然河川生態的動態平衡原理(Principle of dynamic equilibrium)
與河川序級的尺度大小融和以描述河水中動植物縱向的分布模式(Vannote et al.
1980;Lorenz et al. 1997;Aarts and Nienhuis 2003),模式主要是指不受干擾的 自然河川中,來自枯枝落葉的營養物質被河流從上游往下游帶經過連續不斷變化 的物理及化學因子的作用,形成了不同的棲息環境,同時也使棲息於其中的生物 分布呈現連續變化的分布,主要以生產/呼吸作用(P/R,P:gross primary productivity,R:community respiration)的比率分成三區:(1)低序級溪段屬 上游生產區:無脊椎的掠食者和採集者支配整個動物群,生物多樣性受到低溫、
少光和少養分的限制,因此此區的 P/R 率小於 1。(2)中序級溪段屬中游過度 區:此區光線較多,地質植物茂盛,生物多樣性高,河川生物多為自營性,如浮 游、附生性藻類,消費者組成以濾食者與刮食者(grazers)最多,所以此區的
2
P/R 率大於 1。(3)高序級溪段屬下游沉積平原區:來自上游輸入的有機質及營 養鹽充分,過濾型植物占優勢,加上此區河道寬廣、流速較慢,增加了河流棲地 的多樣性,P/R 率小於 1 (Vannote et al. 1980;Allan 1995;Mitsch and Jørgensen 1995)。臺灣位於亞熱帶地區,地形山高平原小,造成河川源短流急,落差可達 數千公尺,年平均降雨量可達2500 公厘,雨量豐沛但時空分配卻不均勻,雨、
旱季降雨量相差極大,也使的臺灣的河川豐、枯流量差異懸殊,但其淡水魚類的 多樣性極其豐富(郭 2003)。臺灣的淡水魚類可分為:(1)初級性(primary):
指終生生活在淡水水域中,多分布於湖泊或河川的中上游,如鯉科、鰍科等;(2)
次級性(secondary):大部分時間生活於淡水,亦可短期生活於海水的魚類,如 慈鯛;(3)周緣級(peripheral):為溯河(鮭科),降海(鰻、鰕虎)洄游之 魚類及偶爾可游入河川的廣鹽性魚類,如鯔等。若再將河川魚類由上游至下游區 分,其分布可分成上游為鱒魚區及鏟頷魚區;中游為石 區及中下游的鯉魚區;
下游或河口的大肚魚、吳郭魚或烏魚區(邵 1993)。由於地形上的特徵再加上 長期為了安全以及經濟上之需求,河川已成為自然與人為力量相互爭奪之地,破 壞了河川生態系的平衡。臺灣河川的問題有(1)水土保持;(2)氣候變遷;(3)
汙染控制;(4)都市開發;(5)生物棲地消失;(6)河川渠道化;(7)強勢 外來物種;(8)生態基流量不足;(9)水資源枯竭(10)災害控制等問題(郭 2003)。而臺灣河川的治理,近年來由於生態環境觀念的興起,已由早期的以防 洪工程為主要目的,演變成兼顧水土保持、環境保護、生態復育的生態工程。由 上可知,河川的生物多樣性及資源對人類的生存,扮演著舉足輕重的角色,若能 加強河川生物多樣性之生態保育及對棲息於河川中之魚類等水生生物的進行保 育、復育等工作並進行長期河川生態研究是必須的。
1.2.魚類族群生態學簡述
魚類生態學的研究首先從物種(species)開始,但是如果單以個體為單位,
魚類會隨環境條件和群體的增衰、年齡結構發生變動而變化,使出生率和死亡率 等許多生命現象無法獲得正確的說明,所以研究魚類通常不以個體為單位,而是 以種內個體棲息在同一生態環境裡所形成的組合群,即族群(population)為對 象。而研究魚類的年齡結構、性比、密度、活動範圍以及環境因子的影響稱為魚
3
類族群生態學(Population Ecology)(殷 1998)。以下敘述族群生態學的主要 研究內容:
(1)年齡:年齡評估(Age assessment)被認為是研究魚類生物學和生態學特性 的基礎,也是分析和評價魚類種群數量變動趨勢的基本依據之一(Ambrose 1983;
赤嶺、麥谷 1997;殷 1998)。年齡查定的方法主要有飼育法、標誌放流法、體 長頻度分析法與年齡形質法等;其中飼育法與標誌放流法易造成生物生長狀態與 自然情形不同;體長頻度分析法對於高齡魚之分離識別較困難,有時僅能獲得相 對年齡;而年齡形質法是最繁瑣,但卻是最能有效查定年齡的方法(黃和郭 1990)。
年齡形質法主要是利用鱗片(scale)、耳石(otolith)、棘(spine)、脊椎骨(vertebra)
與鰓蓋骨(opercular)等骨質組織上留下的輪紋標誌來判定(Ambrose 1983;赤 嶺、麥谷 1997;殷 1998)。最早作出用鰻魚的鱗片判定鰻魚的年齡的是荷蘭的 Leeuwenhoek(1684),而瑞典的 Hans Hederström(1759)利用白斑狗魚(Esox lucius)的脊椎骨最先作出可靠的魚類年齡鑑定(殷 1998)。以鱗片作為年齡鑑 定研究的歷史相當悠久,其優點是不需要犧牲魚隻、對魚體造成的傷害為最小;
淡水魚則是從1930 年代才開始運用(Summerfelt and Hall 1987)。魚類鱗片(scale)
是覆蓋體表位於真皮(dermis)的硬組織,其功能之一是保護體表。鱗片源自中 胚葉真皮結締組織中的骨質層,和表皮細胞無關。鱗片前緣的骨芽細胞分泌骨基 質促使鱗片往前生長,骨芽細胞則扁平化後消失。鱗片上方的骨芽細胞則分泌形 成隆起線,鱗片下方的骨芽細胞形成纖維芽細胞最後變成纖維層板。因此硬骨魚 類之鱗片是由外層堅硬的骨質層(bony layer)和內層的纖維層板(fibrillary layer)
兩層構成,依形狀又可區分為表面光滑的圓鱗(cycloid scale;不具突起棘)和 表面有小棘(cteni)突起的櫛鱗(ctenoid scale)。臺灣馬口魚的鱗片屬於不具 突起棘的圓鱗(Cycloid scale),屬於圓鱗。鱗片的後端露出於表皮,前端則潛 入前一片鱗片之下,形成瓦狀重疊覆蓋體表,其生長線(Growth line, Circuli)呈 同心圓排列。年輪(Annulus, Annual ring)由於生長線的不連續、分枝、切斷、
密集等而形成線狀(line)或帶狀(zone)(黃和郭 1990)。在過去二十年間,
鱗片的研究已由直接計數年輪的方法發展到由鱗片之微化學結構的層面來探討 年齡軌跡,而此法對於海、淡兩棲及底棲的魚種較為有效(Gauldie et al. 1991;
Coutant et al. 1993)。而國內在年齡成長的研究上,主要針對海水魚作探討,
4
對於淡水魚之年齡成長有彭(1986)針對臺灣鏟頷魚的鱗片作初步研究,及林等
(2003 ) 以 鱗 片 和 鱗 序 推 論 台 灣 陸 封 櫻 花 鉤 吻 鮭 ( Oncorhynchus masou formosanus)魚類年齡和成長等。
(2)性比(Sex ratio):魚類性比會隨著成長階段、棲息場所、季節變化而有 所差異,性比之多樣性為魚類對其生活環境多樣性之適應屬性(李 1990),而 繁殖是魚類生活史的重要環節,透過繁殖將資源傳遞給後代,透過了解魚類的性 腺發育程度、性週期等,有助了解其繁殖策略。
(3)密度:是指族群在某一時間、某一區域內個體的絕對數量,透過密度估計 可以了解族群的大小和利用率(殷 1998)。
(3)體長、體重關係:魚類的體長和體重有一定的相關關係。一樣長度的魚,
體重越重,表示魚體越豐滿,營養狀況和環境條件越佳;肥滿度也是魚類體長、
體重關係的另一種表達方式,常用於衡量魚體的豐滿程度、營養狀況和環境條件 的指標(殷 1998)。
1.3.生態工法沿革與意涵
生態工法之觀念源自於德國及瑞士,德國的Seifert 在 1938 年首先提出近自 然河溪整治的概念,表示能夠在完成治理河流的傳統基礎上可以達到接近自然、
廉價並保持景觀美的的一種治理方案(林等 2004)。1951 年,德國 Kruedener 提出「生物工程(Bioengineering)」直接用於溪流治理,並有最接近台灣現況 應用的定義。1962 年,Odum 等人提出將自律行為(self-organizing activities)之 生態工程概念應用於汙水處理工程之中,並且首度提到「生態工程」(Ecological engineering)(Odum 1971)一詞。1985 年,Holzmann 則認為生態治理是一種 可以創造出各式各樣的水深和流路之溪流斷面,注重工程治理和自然景觀的和諧 性,把河岸植被規畫成具有多種小生態環境之治理概念。1989 年,Mitsh 則提出
「生態工程(Ecological engineering)」是「運用生態系之自我設計(self-design)
能力為基礎,強調透過人為環境與自然環境間之互動,達到互利共生(Symbiosis)
之目的。」此外,近自然工法必須考慮到整體河道的安定,並且在容許範圍之內 設法創造出有深潭、急流或淺灘等各式各樣之棲地環境(福留脩文 2001)。國
5
內,1997 年,行政院農業委員會水土保持局提出生態工法的定義為「採用天然 資源為主要材料,以融合週邊地形自然景觀,減少造成生態環境之衝擊為理念設 計,構築可供動植物棲息之空間,創造兼具防災及復育功能之工法。」2002 年,
行政院公共工程委員會更進一步提出生態工程的定義為「基於對生態系統之深切 認知與落實生物多樣性保育及永續發展,而採取以生態為基礎,安全為導向的工 程方法,以減少對生態環境造成傷害(林等 2001)。」,國內外生態工法之沿 革與定義(如表 1)。臺灣河川治理規劃早期以滿足人類需求為主要考量,依照 人類需求及活動範圍進行河川治理工程,以確保河川周邊的產業活動或居家安全。
但是,臺灣因其地形、地質及水文條件特殊,加上地狹人稠的土地利用形態,這 種臺灣河溪的自然特性與人為干擾使溪流處於極不穩定的狀態。近年來,由於經 濟發展,人民所得提高後,防洪排水安全不再為人民追求的唯一目標,因此河溪 整治的目標,除防災、減災外,尚須考量河溪水域環境之保全。常用的於溪流整 治的生態工法包括滯洪池、堤防及護岸、固床工、丁壩工、邊坡整治、擋土工及 水域棲地營造等工法。其目標包含:(1)減低河溪洪峰流量、(2)降低河溪流 速、(3)保護溪岸防止沖刷、(4)維持河溪沖淤平衡、(5)提供水生物洄游 路徑、(6)改善水域水質、(7)提供多樣性生物棲地(巨廷工程顧問公司 2005)。
溪流整治對於生存在溪流生態系中之魚類、蝦類、藻類等生物影響重大,尤以魚 類較常被使用來作為評估指標,魚類具有不同汙染耐受性的種類,不同游泳能力 之種類,以及不同的遷移能力,故可利用此三項特性來反應特定水域即時、綜合 的影響(朱等 2004)。
1.4.實驗魚種簡介
臺灣的淡水魚有一百五十餘種以上(曾 1986),早期以分類與分布為主,
而溪流魚類的生態研究包括族群結構、棲地選擇、食性及生殖特性,到了 1980 年代才逐漸起步,主要對象以櫻花鉤吻鮭(Oncorhynchus masou formosanus)(曹 1988;莊 1988;林等 1989;林等 1990;詹等 1991;林等 1993),台灣纓口 鰍(Crossostoma lacustre)(張 1989)與臺灣石 (Acrossocheilus paradoxus)
(林及許 1990;林及張 1991;孔 2006)作過較詳盡的生殖及棲地利用之研究,
再加上探討平鰭鰍科(Homalopteridae)(曾 1981)之分類與分布,與探討魚類 不同生命週期的不同階段棲地使用與物理棲地因子之量化關係(孫及呂 2010),
6
較少有對魚類生態有更深入之研究。而本島中低海拔溪流的魚類生態,則有於北 部淡水河支流針對臺灣石 (Acrossocheilus paradoxus)(許 1991)、臺灣鏟頷 魚(Varicorhinus barbatulus)(張 1994)及臺灣馬口魚(Candidia barbatus)(顏 1993)進行生殖生物學與族群生態學之研究。為了進一步了解本島溪流魚類的整 體特性,有必要針對棲息於本島各溪流的魚類進行族群生態之研究。臺灣馬口魚
(Candidia barbatus)為臺灣特有種魚類,在系統分類學上是屬於脊椎動物亞門
(Vertabrata)、條鰭魚綱(Actinopterygii)、新鰭亞綱(Neopterygii)、鯉科
(Cyprinidae)、鬚鱲屬(Candidia)(臺灣魚類資料庫)。最早由 T. Regan(1986)
在日月潭發現。臺灣馬口魚分布於臺灣西部各河川中上游及其支流,性喜冷水域,
為雜食性魚類(詹 1989;曾 1990)。形態特徵為體延長而側扁,腹鰭基部後方 具不完全的腹稜。頭中大,吻略短,口端位,口裂中大,向後達眼前緣下方;口 裂末端有一對細小的紅色觸鬚。成熟雄魚的上頜較為前突,且略呈鉤狀。體被細 小圓鱗;側線完整,由鰓蓋的上緣向腹側明顯下彎呈弧形;側線鱗數54-57。各 鰭均無硬棘,背鰭軟條 3(不分枝軟條)+ 7(分枝軟條);臀鰭 3(不分枝軟 條)+ 9(分枝軟條)。體背側呈灰褐色,腹部灰白,體側自前鰓蓋的後緣至尾 鰭基部有一極為明顯的藍黑色縱帶;成熟雄魚的吻部有極為明顯而尖銳的追星,
其頭部下側、腹部、胸鰭及腹鰭均呈鮮紅色;雌魚的體色較淡;各鰭微黃而白色
(臺灣魚類資料庫)。目前關於台灣馬口魚之研究包括:一般性描述(詹 1989;
曾 1990),系統分類與遺傳結構(王 1995)、生殖行為(莊等 1995)、哈盆 溪臺灣馬口魚之生態學初步研究(顏 1993)及繁養殖及環境生物學研究(熊 1999)
外,尚無其他研究報告發表。
1.5. 研究緣起
湳仔溝溪位於桃園縣大溪鎮復興里與新峰里,是石門水庫上游集水區之ㄧ支 流(圖1)。湳仔溝溪由於河道狹窄,易遇雨成災,並在當地居民數次陳情下,
以野溪生態工法治理湳仔溝溪。於2004年3月至2004 年12 月進行第一期中下游 段整治工程,工程內容主要以溪河拓寬,並包含橋樑工程、護岸工程等整治工程,
其中沿河道路亦做了部分調整改道。又於2005年2月至2006年12 月進行第二期中 上游段整治工程(表2)。並自整治前2003年11月至2010年11月整治過後共7年之 時間,依湳仔溝溪上游百吉四號橋、湳仔溝二號橋、龍行橋、鳳儀僑、新福圳至
7
下游懷德橋,共六個測點(圖2),進行湳仔溝溪溪流生態監測與情勢評估,共 調查了魚類、水生昆蟲、附著性藻類、底棲動物之種類、數量、分布、種組成等 生物性資料。湳仔溝溪的魚類組成非常多樣化,於2003至2008年調查至今魚種總 計6目9科28種,包括鯉目鯉科的台灣石 ﹙Acrossocheilus paradoxus﹚、台灣鏟 頜魚(Onychostoma barbatulum)、紅鰭鮊(Culter erythropterus)、粗首鱲﹙Zacco pachycephalus﹚、平頜鱲﹙Zacco platypus)、魚密﹙Hemibarbus labeo﹚、白鱎
(Hemiculter leucisculus)、鯉(Cyprinus carpio carpio)、台灣馬口魚﹙Candidia barbatus﹚、羅漢魚(Pseudorasbora parva)、高體鰟鲏(Rhodeus ocellatus ocellatus)、 革條副鱊(Tanakia himantegus)、日本鯽(Carassius cuvieri)、鯽魚(Carassius auratus auratus)、朱文錦(Carassius auratus);鰍科的中華花鰍﹙Cobitis sinensis﹚、
泥鰍(Misgurnus anguillicaudatus);鯉齒目花魚將科的大肚魚(Gambusia affinis);
合鰓魚目合鰓魚科的黃鱔(Monopterus albus);鱸形目鳢科的斑鱧(Channa maculata);慈鲷科的吉利吳郭魚(Tilapia zillii)、雜交吳郭魚(Orechromis sp.);
及 鰕 虎 魚 亞 目 鰕 虎 科 的 褐 吻 鰕 虎 ﹙ Rhinogobius brunneus ﹚ 、 明 潭 吻 鰕 虎
(Rhinogobius candidianus)、極樂吻鰕虎(Rhinogobius giurinus)、短吻褐斑吻 鰕虎(Rhinogobius rubromaculatus)等。除大肚魚、鯉魚與吳郭魚為外來種外,
其餘多為台灣地區普遍分布的特有種或是原生種(吳等 2008)。
過去有關湳仔溝溪的調查包括:生態評估(朱等 2004;朱 2005;吳等 2008)、 應用生物指標法探討湳仔溝溪以生態工法整治影響之研究(朱等 2004)及石門 水庫上游集水區生態評估之研究:以湳仔溝溪為例(朱等 2008)。過去的研究 主要以概括性及應用生物指數的方式探討湳仔溝溪受到溪流整治之影響,較缺少 針對單一物種受到溪流整治影響之研究,魚類較易於觀察其溪流環境受到變化的 指標生物之一。從2003年11月至2008年10月,臺灣馬口魚在湳仔溝溪之魚種分布,
施工前10%、第一期施工17%、第二期施工17%及施工後13%(圖3),顯示在 這四個時期臺灣馬口魚為屬於臺灣特有種之魚種中分布比率最高之魚種,顯示臺 灣馬口魚為湳仔溝溪之優勢魚種。
過去國內對於馬口魚之相關研究有一般性描述(詹 1989;曾 1990),系統 分類與遺傳結構(王 1995)、生殖行為(莊等 1995)、哈盆溪臺灣馬口魚之生 態學初步研究(顏 1993)及繁養殖及環境生物學研究(熊 1999)外,尚無其他
8
研究報告發表。因此,本研究期望透過長時間調查湳仔溝溪之臺灣馬口魚族群探 討其密度、體長、體重、肥滿度受到溪流整治之影響及年齡與成長和生殖週期之 研究,了解湳仔溝溪的臺灣馬口魚在受到溪流整治的人為因子及環境變動下族群 的結構與動態是否受到衝擊並探討臺灣馬口魚族群在湳仔溝溪的年齡組成、性比 和生殖週期等特性,以其更充分理解臺灣馬口魚在湳仔溝溪之族群生態。
1.6. 研究目的
1、針對湳仔溝溪之臺灣馬口魚,比較正在進行施工的測點(B)湳仔溝二號橋
( 2005.2-2006.2 ) 、 ( C ) 龍 行 橋 ( 2005.2-2006.2 ) 、 ( D ) 鳳 儀 橋
(2004.3-2004.12)、(E)新福圳(2004.3-2004.12)及 A-F 測點中自 2003.11 至 2010.11 未進行施工的時期統稱為(Z)未施工測點之密度、體長、體重 與肥滿度是否有差異,以及水文因子與棲地間的關係,進而了解湳仔溝溪的 臺灣馬口魚在受到溪流整治的人為因子及環境變動下族群的結構與動態是否 受到衝擊與影響。
2、利用鱗片做為年齡形質,探討湳仔溝溪臺灣馬口魚年齡與成長之關係,了解 其年齡結構。
3、利用性比、生殖腺指數(gonadosomatic index,GSI)等的變化,瞭解臺灣馬 口魚在湳仔溝溪之生殖週期變化。
4、利用上述各項因子希望更充分理解臺灣馬口魚在湳仔溝溪之族群生態的變化 及分布。
9
第二章、材料與方法 2.1.研究地點
湳仔溝屬於石門水庫集水區,位於桃園縣大溪鎮內,發源於新峰里北側標高 470m 之山區,向南流經約 2km 之山谷,於腦窟寮南方轉向西南,至百吉後轉向 西北,復於湳仔溝再轉向西南,於龍珠灣流入石門水庫,水系面積約7.0km2主流 長約 10km,河床平均坡降約 1/100(朱等 2008)。本研究以湳仔溝溪流域六測 點,由上游依次往下游為(A)百吉四號橋、(B)湳仔溝二號橋、(C)龍行橋、
(D)鳳儀橋、(E)新福圳、(F)懷德橋等共六站,集水區範圍(圖 1)及調 查測點如(圖 2)所示。
2.2.採樣方法
從2003年11月至2010年11月,依湳仔溝溪上游百吉四號橋、湳仔溝二號橋、
龍行橋、鳳儀僑、新福圳至下游懷德橋,共六個測點,其中各測點以電魚法(電 壓8伏特)由下游至上游以Z字型方向進行15分鐘,100公尺之魚類採樣,使單位 努力捕獲量CPUE(Capture Per Unit of Effort)維持相同,捕獲之魚隻記錄其尾叉 長(fork length)準確至(1mm)及體重準確至(0.1g),最後將魚放回溪流內,
採樣頻度每月採樣壹次。鱗片及生殖腺樣本採集則從2009年9月至2010年11月,
於百吉四號橋(A)、湳仔溝二號橋(B)、龍行橋(C)、鳳儀橋(D)、新福圳(E)及 懷德橋(F)六測點,每月以蝦籠法及電魚法輔助進行魚類之採樣,每月所捕獲魚 隻依其體長及數量分層隨機採樣(6-10隻),以10%福馬林保存,帶回實驗室作 進一步研究。
水質檢測項目以穿越線法記錄採樣測站內底石種類與比例,底質分析係根據 汪(1990)所建議採用之底質分類系統來記錄調查(表3)和水質儀器測量水深、
流速、水溫、酸鹼值(pH值)、導電度、懸浮固體(SS)、溶氧(DO)、生化 需氧量(BOD)、化學需氧量(COD)、氨氮、濁度和總磷共12項。
10
2.3.樣本魚處理
將保存於福馬林中之樣本魚取出,先將魚編號並量測其體全長(total length)、
尾叉長(fork length)及體重後,隨即解剖魚體,依生殖腺辨別性別並量測生殖 腺重及採集鱗片形質。
2.4. 鱗片處理
根據後述結果所判定,鱗片標本採自背鰭後方側線之上之部位(照片1 之 C 部位),每條魚採取8 片左右,從中選取大小相近的 3 片鱗片,放置於含有去離 子水的培養皿中清洗,去除表面汙垢,再放置於離心管中,利用去離子水及超音 坡清洗機(Model 1210R-DTH, Branson)震盪 5 分鐘,去除鱗片上的黏膜、表皮 組織及其他汙物,再以滴管吸出離心管中的去離子,重複此兩步驟三到五次,視 鱗片大小及汙穢程度調整沖洗次數,再將清理乾淨的鱗片放置於兩片載玻片之間,
用標籤紙將兩端黏住,自然陰乾(廖 2001)。將製作完成的鱗片製於顯微鏡(Nikon eclipse E400)下,其物鏡的放大倍率設定4 倍,由 Canon 550D 型相機攝取鱗 片影像以供觀察判讀。
2.5.肥滿度指數之估計
肥滿度為反映魚類個體重量增長程度、生長好壞、肥瘦、營養狀況與對餌 料利用程度的一種指標,也可當作比較不同時期成長差異的依據(Bagenal and Tesch 1978),肥滿度指數(Condition factor;CF)的計算式如下(Htun-Han 1978):
CF= BW FL⁄ 10 ………..(1)
BW:體重 FL:尾叉長
(BW:g;FL:cm)
11
2.6.生殖腺指數之估計
生 殖 腺 重 量 的 變 化 表 示 生 殖 腺 發 育 狀 況 而 可 作 為 生 殖 腺 發 育 之 指 標
(Htun-Han 1978;Delahunty and De Vlaming 1980)。本研究以生殖線指數
(gonadosomatic index;GSI)來判定生殖期,生殖腺指數計算式如下:
GSI GW FL⁄ 10 ………...(2)
GW:生殖腺重 FL:尾叉長
(GW:g;FL:cm)
2.7.性比
魚類性比會隨著成長階段、棲息場所、季節變化與漁具選擇性而有所差異,
性比之多樣性為魚類對其生活環境多樣性之適應屬性(李 1990)。性比的計算 式如下(Lee et al. 1986):
Sex ratio F F⁄ M ……….(3)
(F:雌魚數;M:雄魚數)
2.8.資料處理
1. 以 One-way ANOVA 比 較 正 在 進 行 施 工 的 測 點 ( B ) 湳 仔 溝 二 號 橋
( 2005.2-2006.2 ) 、 ( C ) 龍 行 橋 ( 2005.2-2006.2 ) 、 ( D ) 鳳 儀 橋
(2004.3-2004.12)、(E)新福圳(2004.3-2004.12)及 A-F 測點中自 2003.11 至 2010.11 未進行施工的時期統稱為(Z)未施工測點,在(1)密度,(2)
體長,(3)體重及(4)肥滿度之平均值間有無差異,如有差異,再以Duncan 進行兩兩比較,其數學模式如下:
Χij=μ+τi+εij ………(4)
μ=族群均值,τi =月別效應,εij=試驗誤差。
i=1…I(測點)。j=1…J(月別)。
12
2.利用主成分分析法(Principle Component Analysis)來分析湳仔溝溪調查到的 水文因子與底石組成在環境變異上的貢獻度。主成分分析法是一種多變量統 計方法,可將許多相關的變數轉換成新的獨立變數,這些新的變數為原變數 的線性組合。分析方法為找出原有變數的線性組合,算出共變異數矩陣 S 及 相關矩陣分析R,求出 S 與 R 的特徵值(eigenvalue)、特徵向量(eigenvector)、
貢獻度(proportion)和累積貢獻度(cumulative proportion)。各主成分得點 之計算公式如下(Johnson and Wichern 1988):
Zi=∑ E S ………(5)
Zi=第 i 個主成分得點,Eij=第 i 個主成分得點之第 j 個因子的特徵向量,Sij=第 i 個主成分得點之第 j 個因子的標準化值,其中
∑ E =1
再將特徵值依大小順序排列,設分別為λ1,λ2……λp,求出對應的單位特徵向量 a1,a2……ap,其中ai‵ai=1 且 ai‵aj=1,則
y1= a1X= a11X1+ a12X2+……+ a1pXp 為第一主成分 y2= a2X= a21X1+ a22X2+……+ a2pXp 為第二主成分
………
Yp= apX= ap1X1+ ap2X2+……+ appXp 為第 p 主成分
接著再簡化多變量資料的變數個數並求出解釋的變異比例,若只取最大的 q 個主 成分替代原有的p 個變數,則這 q 個主成分解釋的變異比例為:
R2=λ λ λ ⋯ λ
λ λ λ ⋯ λ
λ λ λ ⋯ λ
………(6)
當所能夠解釋的變異比例值越高,顯示這q 個主成分所能解釋各因子與目的變數 之間的可信度越高。
本研究所有測量與計算之數據,皆輸入Excel,統計軟體使用 spss 12.0。
13
第三章、結果 3.1.溪流整治對臺灣馬口魚之影響
3.1.1. 族群密度
本研究將2003 年 11 月至 2008 年 10 月湳仔溝溪魚種組成依施工期距分成四 個時期,顯示施工前臺灣馬口魚佔10%、第一期施工 17%、第二期施工 17%及 施工後 13%,在這四個時期皆為湳仔溝溪魚類屬於臺灣特有種中分布比率最高 的魚種,表示臺灣馬口魚為湳仔溝溪之優勢魚種(圖3)。而在這四個時期各測 點以圓點的大小表示以電魚法調查之臺灣馬口魚平均密度(隻/15 分鐘/100 公尺)
變化:在未施工時期密度以最下游測點懷德橋最高(5.5 隻/15 分鐘/100 公尺),
其次為龍行橋(5.0 隻/15 分鐘/100 公尺)、百吉四號橋(3.0 隻/15 分鐘/100 公尺)、湳仔溝二號橋(2.8 隻/15 分鐘/100 公尺)、鳳儀橋(1.5 隻/15 分鐘/100 公尺),最低為新福圳(0.8 隻/15 分鐘/100 公尺);第一期施工期間,由於受 到艾利風災(2004 年 8 月)之影響,各測點密度都呈現大幅下降的趨勢,施工 之中下游測點鳳儀橋密度為(1.2 隻/15 分鐘/100 公尺)及新福圳甚至完全沒有 調查到馬口魚(0.0 隻/15 分鐘/100 公尺),而最下游測點懷德橋的密度卻達到 最高(12.1 隻/15 分鐘/100 公尺),其他測點密度依次為龍行橋(2.5 隻/15 分 鐘/100 公尺)、湳仔溝二號橋(2.2 隻/15 分鐘/100 公尺)及百吉四號橋(1.5 隻/15 分鐘);第二期施工期間,密度最高為未施工測點百吉四號橋(16.4 隻/15 分鐘/100 公尺),進行施工測點湳仔溝二號橋(2.4 隻/15 分鐘/100 公尺)及龍 行橋(7.0 隻/15 分鐘/100 公尺),進行施工之測點馬口魚密度受到施工影響不 大,而施工後之鳳儀橋(1.2 隻/15 分鐘/100 公尺)及新福圳(4.1 隻/15 分鐘/100 公尺),最下游懷德橋則為(7.3 隻/15 分鐘/100 公尺),顯示艾莉風災過後各 測點密度都有恢復的趨勢;施工過後,密度最高還是未施工之百吉四號橋(18.1 隻/15 分鐘/100 公尺),其次為龍行橋(11.5 隻/15 分鐘/100 公尺)、最下游之 懷德橋(3.1 隻/15 分鐘/100 公尺)、新福圳(1.5 隻/15 分鐘/100 公尺)、鳳儀 橋(0.7 隻/15 分鐘/100 公尺),最低為湳仔溝二號橋(0.5 隻/15 分鐘/100 公尺)。
從以上結果得知,密度在進行施工時期各測點變化皆不相同(表4)。
14
在溪流整治影響之研究中,自2003 年 11 月至 2010 年 11 月共調查到 1179 尾臺灣馬口魚,其中平均密度最大為(C)龍行橋7.00±3.83 隻/15 分鐘/100 公尺,
其次依序為(Z)未施工測點 4.79±6.03 隻/15 分鐘/100 公尺、(B)湳仔溝二號橋 2.43±1.90 隻/15 分鐘/100 公尺、(D)鳳儀橋 1.20±1.93 隻/15 分鐘/100 公尺,最 低為(E)新福圳0.00±0.00 隻/15 分鐘/100 公尺,經檢定過後平均族群密度達到顯 著差異(p=0.02<0.05),再進行兩兩比較後發現(C)、(Z)及(B)測點之平 均密度沒有差異,(Z)、(B)、(D)及(E)測點之平均密度沒有差異,但(C)
與(D)、(E)測點之平均密度有差異(表 5、圖 4)。
從2003年至2010年臺灣馬口魚的數量變化中,在第一期施工整治數量每月約 在20-40隻,艾利風災過後數月數量下降,至第二期施工整治期間數量有呈現緩 慢恢復的情況,甚至有些月份數量到達60隻,但在整治過後,數量持續下降,至 2010年每月數量只約在20隻,可見臺灣馬口魚的族群數量經過溪流整治後有下降 的趨勢(圖5)。重量上的變化則是,在施工前臺灣馬口魚的平均重量在15-20g,
第一期施工期間可能受到艾利風災與溪流整治的影響開始下降,之後臺灣馬口魚 之重量隨著整治時間逐漸下降,至2010年臺灣馬口魚的平均重量只約在5g左右
(圖6)。以上結果顯示,臺灣馬口魚族群在受到艾利風災、溪流整治與魚群之 間競爭壓力下發生變化。
3.1.2.體長
臺灣馬口魚體長範圍介於 1.8cm-14cm 之間,大部分集中在 3-4cm 之間
(n=497)。其中(B)湳仔溝二號橋體長範圍 3-13cm(n=17),以 7、9cm 所占比 例最高(41%、35%);(C)龍行橋體長範圍 3-12cm(n=49);(D)鳳儀橋體長 範圍集中在 3-6cm 及 12-14cm 之間(n=12);(E)新福圳未調查到馬口魚而(Z) 未施工測點體長範圍3-14cm(n=1101)(圖 7)。平均體長依序為(D)7.00±4.57cm
>(B)6.99±2.25cm>(Z)6.07±3.15cm>(C)5.58±2.63cm>(E)0.00±0.00cm,經檢定 過後平均體長達到顯著差異(p =0.00<0.05),再進行兩兩比較後發現(D)、(Z)、
(B)及(C)測點之平均體長沒有差異,但與(E)測點之平均體長有差異(表 6、
圖 8)。
15
3.1.3 體重
體重範圍介於 0.1-44g 之間,大部分集中在 0-10g(n=949)。除了(E)未調 查到馬口魚外,(B)、(C)、(D)及(Z)測點之體重範圍皆集中在0-10g 之間(圖 9)。
平 均 體 重 依 序 為 (D)11.81±13.30g > (B)5.97±7.26g > (Z)6.36±8.82g > (C)4.58±5.11g > (E)0.00±0.00g , 經 檢 定 過 後 平 均 體 長 達 到 顯 著 差 異 (p
=0.02<0.05),再進行兩兩比較後發現(D)、(Z)測點之平均體重沒有差異;
(Z)、(B)及(C)測點之平均體重沒有差異;(C)及(E)測點之平均體重 沒有差異,但(D)、(B)及(E)測點之平均體重有差異(表 7、圖 10)。
3.1.4 肥滿度
肥滿度範圍介於0.1~183.1 之間,大部分集中在 10-30 之間(n=909)。其中 (B)湳仔溝二號肥滿度範圍10-40(n=17),主要集中在 10-20(59%、71%);
(C)龍行橋肥滿度範圍10-70(n=49);(D)鳳儀橋肥滿度範圍集中在20-30 及 80-90 之間(n=12);(E)新福圳未調查到馬口魚而(Z)未施工測點肥滿度範圍 10 大於 100(n=1101)(圖 11)。平均肥滿度依序為(D)37.75±31.75>(C)28.94±19.21
>(Z)24.64±24.91>(B)14.56±6.44>(E)0.00±0.00,經檢定過後平均體長達到顯 著差異(p =0.00<0.05),再進行兩兩比較後發現(D)、(C)及(Z)測點之平均 肥滿度沒有差異;(C)、(Z)及(B)測點之平均肥滿度沒有差異;(B)及(C)
測點之平均肥滿度沒有差異,但(D)及(E)測點之平均肥滿度有差異(表 8、
圖12)。
3.1.5.環境因子
將(B)、(C)、(D)、(E)及(Z)測點調查到的水深、流速、水溫、
pH 值、溶氧、導電度、濁度、氨氮、COD、BOD、懸浮粒子、總磷與底石組成 共13 項環境因子進行平均值之計算(表 9、表 10),並分別在各測點進行主成 分分析。在(B)測點發現必須到第四主成分達到累積總變異的 94.1%,四個主 成分的特徵值都大於1,主成分 1 的變異因子是圓石,主成分 2 的變異因子是懸 浮固體,主成分 3 的變異因子是總磷,主成分 4 的變異因子是導電度(表11)。
(C)測點發現必須到第五主成分達到累積到總變異的 94.7%,五個主成分的特
16
徵值都大於 1,主成分 1 的變異因子是 COD,主成分 2 的變異因子是流速,主 成分3 的變異因子是水深,主成分 4 的變異因子是導電度,主成分 5 的變異因子 是氨氮(表12)。(D)測點發現必須到第五主成分達到累積到總變異的 94.1%,
五個主成分的特徵值都大於1,主成分 1 的變異因子是小漂石,主成分 2 的變異 因子是溶氧,主成分3 的變異因子是導電度,主成分 4 的變異因子是酸鹼值,主 成分 5 的變異因子是 BOD(表 13)。(E)測點發現必須到第五主成分達到累 積到總變異的92.9%,三個主成分的特徵值都大於 1,主成分 1 的變異因子是溶 氧,主成分2 的變異因子是大漂石,主成分 3 的變異因子是濁度,主成分 4 的變 異因子是圓石,主成分5 的變異因子是水溫(表 14)。(Z)測點發現必須要到 第七主成分才會累積到總變異的71.5%,七個主成分的特徵值都大於 1,主成分 1 的變異因子是小漂石,主成分 2 的變異因子是水溫,主成分 3 的變異因子是沉 積砂土,主成分4 的變異因子是導電度,主成分 5 的變異因子是 BOD,主成分 6 的變異因子是水深,主成分7 的變異因子是酸鹼值(表 15)。
3.2.年齡與成長 3.2.1. 研究樣本簡介
為了瞭解臺灣馬口魚在湳仔溝溪的年齡與成長之關係,本研究於 2009 年 9 月至2010 年 11 月利用蝦籠法及電魚法採集臺灣馬口魚,共捕獲 192 隻,尾叉長 範圍2.9-17cm,平均尾叉長 7.8cm;體重範圍 0.3-53.0g,平均體重 7.1g(圖 13)。
其中,可辨別性別臺灣馬口魚107 隻,雄魚 78 隻,平均尾叉長 9.1cm;雌魚 29 隻,平均尾叉長10.1cm(圖 14)。
3.2.2.尾叉長與體重之關係
尾叉長(FL)與體重(W)之間呈 W=aLb 的關係,經自然對數轉換後呈 lnW=lna+blnL 之直線關係(Ricker,1973),由於本研究可分辨出雌雄魚之尾數 並不多,故將雌雄魚合併計算,經迴歸係數顯著性檢定之結果,顯示迴歸直線成 立(圖15)。尾叉長與體重之關係式如下:
W 6.9 10 FL . ………(5)
(W:g;FL:cm)
17
3.2.3. 年齡查定鱗片之選定
使用鱗片作為年齡查定之形質時,為了減少鱗片形態之差異,必須決定固定 之採鱗部位。鱗片因著生的位置不同,其形狀大小均有差異(劉和曾,1972)。
本研究選取5 尾尾叉長相近之魚(平均體長 15.2±1.3),並選取 6 個部位(照片 1),每個部位採取 3 片鱗片來分析各部位鱗片之鱗徑平均值及變異系數
(coefficient of variation;CV)的變動情形,變異係數求法如下:
CV S x⁄ 100%...(6)
S:標準偏差 x:平均值
結果顯示,以 C 部位的平均鱗徑最大,變異係數最小,表示此區域鱗片的變 動情形最小,故選背鰭後方測線以上之 C 部位,即背鰭後方,側線以上之部位,
作為採鱗之部位(圖 16)。測量軸則是採以鱗片之焦點(focus)向右與邊緣之 垂直線為測量軸(照片2),量取鱗徑(R)及各輪徑(Rn)之長度,此 3 片鱗 片之輪徑及各輪徑之平均值,即為該尾魚之輪徑及各輪徑之代表值。
3.2.4. 輪徑長與尾叉長之關係
瞭解年齡形質與體長間的關係其目的在於藉由逆算(back-calculation)回推 早期各年齡的體長(Casselman 1990),且有助於了解其成長史各階段之成長狀 況(Bartlett et al. 1984)。馬口魚之輪徑長(R)與尾叉長(FL)的關係呈簡單 迴歸直線L=a+bR 之傾向(圖 17)。
其關係式為:FL=1.9215+0.0094R ………(7)
(FL:cm ;R:μm)
3.2.5.輪紋之規則性
年齡形質上輪紋數及輪紋徑若具有(1)體長愈長,輪紋數愈多,(2)同一 個體之各鱗片的輪紋數一定,且各輪紋的相對位置大略一致,(3)等輪紋數之 異個體間,各輪紋的相對位置大略一致,(4)相鄰輪紋間,各輪紋徑組成彼此 均有所分離且各平均輪紋徑均有差異、等規則性,則該形質為有效之年齡形質
18
(Tan et al. 1965; Lai and Liu 1974)。本研究結果顯示:(1)馬口魚具有隨著 輪徑越大,尾叉長也增加的現象(圖17)。(2)隨著尾叉長增加,輪紋數也會 增大,且輪數群峰形分離(圖18)。(3)隨著輪徑增加,輪紋數也會增大,且 輪數群峰形分離(圖19)。(4)2、3 輪群其輪紋與輪徑之迴歸直線彼此分離,
顯示同輪數異個體間,各輪紋形成之相對位置大致相似(圖20-a、b)。綜合上 述結果可以判定鱗片為一有效之年齡。
3.2.6.輪紋形成時期
輪紋形成期是根據鱗片的邊緣成長率(marginal increment;MI)月別平均值 之變化來推定,邊緣成長率之計算式如下(Joung and Chen 1992):
MI R r ⁄ r r ………(8)
R:輪徑
rn:最外一輪之輪紋徑 rn-1:次外一輪之輪紋徑
(R、rn、rn-1:μm)
由於本研究樣本數較少,將雌、雄魚樣本合併處理,並採取季別方式探討輪 紋形成期,結果顯示,1-3 輪魚之邊緣成長率平均值 1-3 月為 0.27,4-6 月為 0.26,
並在7-9 月達到最高值 0.30,10-12 月隨即下降至 0.24,一年只有一個高峰,因 此推測鱗片之輪紋一年形成一輪,形成期大約在7 至 9 月間(圖 21)。
3.2.7.生殖腺指數的變化
生 殖 腺 重 量 的 變 化 表 示 生 殖 腺 發 育 狀 況 而 可 作 為 生 殖 腺 發 育 之 指 標
(Htun-Han 1978;Delahunty and De Vlaming 1980),本研究以生殖線指數
(gonadosomatic index;GSI)來判定生殖期。其中可辨別出雌雄之樣本數為 107 隻,根據生殖腺指數之季別變化得知,雄魚(n=77)生殖腺指數 1-3 月為最小值 0.067,4-6 月為最大值 0.113,其值在四季中變化不大,雄魚無明顯的生殖盛期;
雌魚(n=29)生殖腺指數在 7-9 月達到最大值 0.346,10-12 月則下降至 0.112,顯
19
示雌魚的生殖季為7-9 月(圖 22)。將雌、雄魚的生殖腺成熟情形分為未成熟、
開始成熟及完全成熟三階段,觀察到雄魚的最小成熟尾叉長皆約為 7.3cm,雌魚 為 7.5cm,(圖 23),而雌、雄魚在尾叉長皆約 6-7cm 達到 50%性成熟體長(圖 24)。
3.2.8.肥滿度指數的變化
雄魚肥滿度季別變化不大,在1-3 月為最小值 9.884,之後肥滿度值持續上升 至10-12 月達到最大值 11.211;雌魚肥滿度的季別變化可知在 4-6 月為最小值 5.190,
隨後在7-9 月略為上升至 8.541,在 10-12 月達到最大值 11.062。顯示肥滿度變化 是在生殖腺指數達到最高值前一季為肥滿度指數為最小值,之後持續呈現上升的 趨勢(圖25)。
3.2.9.性比
從 2009 年 9 月至 2010 年 11 月共 192 隻,其中可辨別性別之馬口魚為 107 隻(雄魚:78;雌魚:29),性比為 0.271。經卡方分布檢定(Chi-suuare test)與 雌、雄性比1:1 有顯著性差異(X2=21.736> X2(0.01,1)=6.64)。再將尾叉長以 1cm 作組距,將雌雄魚尾叉長用 t-test 比較,雄魚平均尾叉長8.5±2.3cm,雌魚平均 尾叉長 9.3±2.39cm,兩性平均尾叉長間無顯著性差異(t=1.73,df=104,p>0.05),
再以 K-S 法比較雌雄魚的尾叉長分布後,發現兩性尾叉長分布也無顯著性差異
(p=0.55>0.05)(圖 26)。
3.2.10.年齡組成
利用鱗片形質判讀馬口魚的年齡,有 4 個年齡群(0-3 歲),年齡組成分別 為0 歲魚 23.4% ,1 歲魚 24.7% ,2 歲魚 41.8%及 3 歲魚 10.1% (圖 27)。再以 年齡別將馬口魚分成 3 個年齡組,探討各年齡組的馬口魚數量,發現 0-2 歲魚數 量雄魚多於雌魚,3 歲雌雄魚數量則相同(圖 28)。
20
第四章、討論 4.1.溪流整治對臺灣馬口魚族群之影響
臺灣馬口魚在施工前、中、後在六個測點密度的變化中,在未施工時期,其 六個測點的環境經利用生化需氧量(BOD5)、溶氧量(DO)、氨氮(NH3-N)及 懸浮固體(SS)等四項參數組成並根據其數值來對汙染程度分成:Non-impaired、
Slightly impaired、Moderately impaired、Severely impaired 四類的河川汙染指標法
(River Pollution Index;RPI),評估顯示是第一級 Non-impaired(未受汙染,稍 受汙染),顯示在施工前湳仔溝溪的環境是良好的(朱等,2008),而此時馬口 魚的密度是以最下游測點(F)懷德橋最高(5.5 隻/15 分鐘),其次為(C)龍行 橋(5.0 隻/15 分鐘)、(A)百吉四號橋(3.0 隻/15 分鐘)、(B)湳仔溝二號橋(2.8 隻/15 分鐘)、(D)鳳儀橋(1.5 隻/15 分鐘),最低為(E)新福圳(0.8 隻/15 分 鐘)的分布狀態。在第一期施工進行期間,由於受到艾利風災之影響,暴漲的溪 水會沖刷原來著生於溪床底石上的藻類及水棲昆蟲,引起魚類短期的食物缺乏,
溪床中的底石組成也可能改變,溪岸邊的濱溪植物也被沖光,使溪流的有機質數 量改變(林等 2002),可能為上游各測點馬口魚密度皆下降之原因,但是最下游 之 F 測點馬口魚密度反而上升,推測為颱風帶來的雨量將上游的魚往下游帶,或 是石門水庫區中的水位上升,讓庫區中存在的馬口魚回游至 F 測點,形成此點密 度上升之故,與應用生物整合指標法(Index of Biotic Integrity,IBI)(Karr 1981;
1991)探討湳仔溝溪以生態工法整治影響之研究指出其 F 測點在 2004 年 9 月時 有相當高的IBI 等級(Non-impaired)(朱等 2007)相同,顯示在施工期間發生 颱風使溪流及庫區水位的變動增大,也對馬口魚的密度分布影響極大。
第二期施工期間,在颱風過後環境漸漸恢復,各測點馬口魚的密度皆比受颱 風影響期間上升,而密度較高的是以未施工測點的 A 及 F 兩區。施工後,密度 以最上游之未施工測點 A 密度最高,其次依序為第二期施工測點 C、最下游之 未施工測點 F、第一期施工測點 E、D,最低為第二期施工測點。由以上結果可 知,湳仔溝溪馬口魚在施工過後其密度的分布與施工前不同,除了 C 點在施工 後未施工測點的馬口魚密度皆比施工測點高,顯示施工整治對馬口魚的族群密度 是有影響的,與探討嘉義縣赤蘭溪針對當地魚類調查發現在施工前後魚類分布的 變化有改變(莊等 2004)及探討後番子坑溪在經過生態工法整治過後之魚類變
21
化情形中粗首鱲、臺灣鏟頷魚、臺灣石魚賓以及臺灣馬口魚在工程進行中出現頻 度出現大幅下降之情形相似(朱等 2007)。
本研究為了瞭解溪流整治對馬口魚族群的影響,再將進行整治期間的四個測 點對其他未進行整治的測點(Z)針對其密度、體長、體重及肥滿度平均值四個 參數探討整治是否對馬口魚產生影響,由結果知 C 點馬口魚的密度平均值為最 大(7.00±3.83 尾/15 分鐘),但其平均體長(5.58±2.63cm)及平均體重(4.58±5.11g)
卻偏小,平均肥滿度高(28.94±19.21),推測其原因可能為 C 測點在施工中水域 型態為淺流、岸邊緩流,底石組成種類以砂土、礫石為主(吳等 2008),在整治 中有設「避魚柵」及「生態島」的整治項目(照片3-a、b),「避魚柵」其功能為 河道在低流量時提供深潭供水生動物棲息,但在中高流量或洪水時,必須提供水 生物避難之空間;「生態島」由卵塊石所組成,其空隙因有泥砂淤積,已有草本 植物生長,再以大顆粒塊石舖設防止受水流沖刷,可供鳥類、兩棲類棲息,並提 供生物安全隱避的的空間(經濟部水利署北區水資源局),而生態島上加上岸邊 種植豐富的植生提供了豐富的營養來源,而馬口魚的幼魚喜愛成群地聚集在溪流 兩岸的緩流處覓食(臺灣魚類資料庫),讓C 測點的密度及肥滿度平均值高,但 是體長及體重平均值偏小之故。
在 C 測點水質資料利用主成分分析後可發現影響 C 測點環境變異的水文因 子主要來自COD、流速、水深、導電度與氨氮,COD 大部分應用在測定河水表 面水中有機污染物的含量,可藉此得知水質的好壞,水溫的變異也會影響水中的 耗氧量,氨在水中有未解離氨(NH3)及解離氨(NH4+
),一般未解離氨是比較 有毒(Burkhalter and Kaya 1977;Armston et al. 1978),此形式為脂溶性而很容易 擴散到細胞膜(Fromm and Gillette 1958),此外,氨經由硝化細菌的作用,產生 中間產物-亞硝酸,亞硝酸會引起魚類的棕血症,將血紅素氧化成為變性血紅素,
降低血紅素的攜氧能力(熊 1999),而馬口魚是喜愛棲息於水質清澈、流速較緩 且淤積度較淺的水域中,所以C 測點中 COD、流速、水深與氨氮的變異情形可 能會影響馬口魚對棲地的選擇性。
D 測點馬口魚平均密度不高(1.20±1.93 尾/15 分鐘),但其體長(7.00±4.57cm)、
體重(11.81±13.30g)及肥滿度(37.75±31.75)平均值皆為最大,其原因可能為 D 測點施工中水域型態為淺流、岸邊緩流,底石組成種類以砂土、礫石為主,而 在施工過程中底質組成中砂土數量比例提高(吳等 2008),可能影響 D 測點馬口
22
魚環境棲地的變異因子可能為溶氧,臺灣馬口魚的致死氧量為 0.77mgO2/L(熊 1999),比臺灣石 的1-1.2mgO2/L(彭 1993)及青魚(Mylopharyngodon piceus)
的0.8-1.1 mgO2/L(廖等 1991)為低,所以水中溶氧的變化可能會造成馬口魚密 度降低之原因,而此區整治工法的「淺灘深潭營造」其淺灘可作為魚群的餌料場,
潭區則提供一個休息庇護的之地(郭 2003),適合臺灣馬口魚成魚生存之區域。
E測點在施工整治期間無法記錄到馬口魚,原因為 E 測點在施工整治期間 遭遇到艾利風災之侵襲再加上 E 測點中有1.5 公尺高防砂壩的存在,經由主成分 分析發現濁度、溶氧可能為為 E 測點影響馬口魚的環境變異因子,濁度的變化可 能為艾利風災帶來的豐沛雨量使溪中的大量沙土被帶起往下游沖,也讓水中的溶 氧量降低,使喜愛低淤積度且水質清澈棲地類型的馬口魚受到影響,而防砂壩對 水生生物的影響可能造成(1)洄游生物路逕阻隔;(2)魚類族群縮小與區隔化;(3)
魚類棲息地單調化;(4)水生生物族群基因庫縮小與區隔化,使得 E 測點其環境 可能不適於馬口魚之生存(張與林 1994)。
此外,在先前湳仔溝溪的研究「應用生物整合指標法探討湳仔溝溪以生態工 法整治影響之研究」中指出,優勢度指數(C) 數值越大顯示有明顯的優勢種出 現,湳仔溝溪之優勢度指數(C)在施工後有明顯的變動,且呈上揚趨勢,施工區 測點C、D、E 有較高數值形成優勢種出現;數值越大顯示有較豐富的種類出現,
湳仔溝溪的種歧異度指數(H’)呈現著下降的趨勢,尤以施工測點 C、D、E 等測 站較為明顯;而 IBI 值在施工區的 C、D、E 測點也有明顯的變化,呈現著下降 的趨勢:C 測點下降至Slightly impaired 及 Moderately impaired 等級之間;D、E 測點皆下降至Moderately impaired 等級(朱等 2007)。本篇研究發現,C 點進行 施工時,臺灣馬口魚佔全部魚種比例最高(25%),推測此點之優勢種可能為臺灣 馬口魚;D 測點及 E 測點在施工中魚種分布以鯉魚比例最高(35%,52%),E 測 點施工中並無調查到臺灣馬口魚之存在;臺灣馬口魚的平均密度變化也與IBI 的 變化趨勢類似,C 測點的密度皆比 D、E 兩測點高。
在湳仔溝溪實施的溪流施工整治中,施工項目「生態島」、「避魚柵」、「生態 滯洪池」及「斜坡式跌水工」等工法,其作用都是為了提高了河道水流阻力,以 淨化湳仔溝溪的水質,朱等(2010)針對湳仔溝溪應用微濕地生態工程技術去除 水質污染研究指出,「生態島」的除污型水生植物發揮效果使氨氮等營養鹽有下 降的趨勢,「滯洪池」、「斜坡式連續性跌水工」、「淺灘深潭營造」等工法發揮了
23
降低懸浮固體的功效,顯示在湳仔溝溪淨化水質上有正面的影響。但臺灣的溪流 因地形、氣候及人為的干擾下,生態工程施作的效果都只維持二、三年,湳仔溝 溪整治過後幾年陸續在艾利等颱風影響下,生態島等工程也因棲地受到干擾,無 法發揮其效果,湳仔溝溪的臺灣馬口魚在整治過後密度等因子也較難恢復至施工 前之情況,綜合結果,湳仔溝溪的臺灣馬口魚在河溪整治的過程中,馬口魚會因 棲地環境以及施工工法的不同在密度、體長、體重及肥滿度上出現不同的反應,
但溪流整治對於臺灣馬口魚族群所帶來的影響是有相當的衝擊存在,湳仔溝溪位 於石門水庫上游,再加上有防砂壩的阻隔,都會使臺灣馬口魚在洄游上受到干擾,
是否有造成棲地零碎化的效應,日後對於臺灣馬口魚族群的保育及復育更需進一 步的觀察及討論。
4.2.年齡與成長
一般常用來作為年齡查定的年齡形質有鱗片、耳石、脊椎骨、棘與尾椎骨等 硬組織(能勢等 1988)。關於淡水魚的年齡鑑定,除少數無鱗魚或鱗片上的年 輪特徵不明顯之種類外,通常都以鱗片作為年齡鑑定的材料(李等 2010),近 年來有研究報告指出以鱗片作為年齡查定的形質並不適用,主要原因是因為易產 生低估魚類年齡的誤差,這些現象不論是在海水魚或淡水魚中皆有例子,而低估 的情形主要以海水魚高於淡水魚(Barnes and Power 1984;Beamish and McFarlane 1983)。原因為就鱗片形質而言,某些魚種在先天上有有判讀上的困難,如年輪 的不明顯、斷裂、重合、偽輪或再生鱗等現象,但這並不是造成鱗片較不被接受 作為年齡查定形質的主因,1983年Beamish and McFarlane 的研究報告表示,魚 體的成長有一定的極限,當魚體成長到一定的限度後,鱗片年輪也伴隨著無法繼 續隨年齡再增長下去,以致於高齡魚之年紀不易判讀,容易造成低估,不過由於 本樣本於之鱗片形質為年齡組成以淡水魚的低齡魚為主(0-3歲),相較於高齡 魚的形質,較不易產生誤判,且鱗片形質的輪徑具有等規則性,是為良好且有效 的年齡形質。
由於輪紋判讀法在判讀過程中容易造成人為判定之誤差(Lai and Gunderson, 1987),為了瞭解誤差,本研究取標本魚鱗片樣本,每一鱗片皆判讀兩次,結果 顯示兩次判讀間輪數相同者約88.5%(表16),且其餘均相差在一輪以內。在李 氏現象(Rosa Lee’s phenomenon)方面,理論上隨著魚體的生長,身體各部位亦