楊靜，史子莫，華杰，占陽，銀旭紅，劉婧，李云娟，胡茂林

（1.南昌大學(xué)生命科學(xué)學(xué)院，江西 南昌 330031；2.江西省水產(chǎn)技術(shù)推廣站，江西 南昌 330046；3.安義縣農(nóng)業(yè)農(nóng)村局，江西 南昌 330500）

魚類是水生生態(tài)系統(tǒng)食物鏈的頂級(jí)或近頂級(jí)群落，易鑒定、個(gè)體較大、壽命長、對(duì)水質(zhì)反應(yīng)敏感[1]，是河流生態(tài)系統(tǒng)的重要生物類群[2]。然而，在生物因素和非生物因素影響下，其資源量和補(bǔ)充量不斷下降，魚類物種組成和分布也發(fā)生改變[3]。魚類群落的空間分布特征是魚類與局域生境在長期進(jìn)化過程中形成的一種適應(yīng)性特征，它反映了環(huán)境對(duì)物種空間資源利用的影響[4]。研究魚類的群落結(jié)構(gòu)特征及其與環(huán)境因子的關(guān)系是保護(hù)魚類多樣性的前提[5]。

修河為鄱陽湖水系五大河流之一，位于江西省的西北部，其水位洪枯變化明顯。潦河為修河一級(jí)支流，河系發(fā)達(dá)，支流眾多。龍安河系潦河一級(jí)支流，自北向南注入潦河，流域面積305 km2，主河道長57.3 km[6]。1985 年修河記載有魚類75 種[7]，1995年記錄魚類34 種[8]，2008—2009 年調(diào)查到魚類41種[9]，主要原因是水庫阻隔（柘林水庫）、不合理的捕撈方式（電魚）、挖沙作業(yè)及水質(zhì)污染等造成魚類資源減少，同時(shí)小型魚類如類成為了該流域的優(yōu)勢(shì)種群，而大型魚類（如四大家魚）資源量下降[9]。這些資料主要分析了魚類的生態(tài)類型及區(qū)系組成，而對(duì)魚類群落結(jié)構(gòu)特征及其與環(huán)境之間的關(guān)系調(diào)查甚少。本研究基于2020 年7—8 月潦河永安段（圖1）環(huán)境和漁獲物的調(diào)查數(shù)據(jù)，主要闡明：（1）魚類物種組成與空間分布特征；（2）魚類群落結(jié)構(gòu)的空間差異；（3）魚類群落與環(huán)境因子之間的關(guān)系。

圖1 潦河永安段采樣點(diǎn)分布圖Fig.1 Distribution of sampling sites in the Yong′an section of the Liaohe River

1 材料與方法

1.1 研究區(qū)域概況

龍安河、北潦河北支及潦河是永安段的主要河流，其中龍安河為潦河的一級(jí)支流，北潦河北支為潦河的二級(jí)支流，共同匯入潦河，對(duì)安義、永修兩縣的經(jīng)濟(jì)發(fā)展和城市供水起重要作用。該流域?qū)賮啛釒Ъ撅L(fēng)氣候區(qū)，4—6 月為雨季，暴雨集中，平均降水量1 778.0 mm。

1.2 樣點(diǎn)設(shè)置與標(biāo)本采集

本次調(diào)查共設(shè)11 個(gè)采樣點(diǎn)（圖1），2020 年7—8 月調(diào)查了每個(gè)采樣點(diǎn)的漁獲物和環(huán)境因子。其中，S1～S6 采樣點(diǎn)位于龍安河上游，為山區(qū)性河流，河谷深切，河道彎曲，河面寬一般小于50 m,河床多礫石、粗沙；S7 位于龍安河下游，其水槽逐漸開闊，水流趨于平緩，河面寬一般80～150 m,河床多細(xì)沙、淤泥。S8 處于龍安河與潦河匯合處，S9 和S10 處于潦河中下游河段，S11 位于北潦河北支下游，這四個(gè)采樣點(diǎn)河道較寬淺順直，水流平緩，河床多粗沙、卵石[6]。水深小于1.5 m 的采樣點(diǎn)（S1-S6）用電捕法（捕撈時(shí)長15 min）采集，水深1.5 m 以上的采樣點(diǎn)（S7-S11）用絲網(wǎng)（長70 m、高1.5 m、網(wǎng)目1 cm，捕撈時(shí)長1 h）捕撈；同時(shí)，使用便攜式多參數(shù)水質(zhì)分析儀（AP-2000）現(xiàn)場(chǎng)監(jiān)測(cè)各樣點(diǎn)的經(jīng)緯度、海拔、溶解氧、水溫等環(huán)境因子，河寬用測(cè)距軟件測(cè)量?，F(xiàn)場(chǎng)鑒定魚類物種，記錄每種魚的尾數(shù)和質(zhì)量[10-12]。

1.3 數(shù)據(jù)處理

魚類的生態(tài)類型劃分參考文獻(xiàn)[13-16]，魚類的瀕危狀況和中國特有種依據(jù)FishBase（Search Fish-Base）劃分。采用相對(duì)重要性指數(shù)（IRI）確定魚類優(yōu)勢(shì)種，將IRI＞1 000 劃分為優(yōu)勢(shì)種[17]；運(yùn)用Shannon-Wiener 指數(shù)（H）、Margalef 指數(shù)（D）和Pielou 指數(shù)（J）分析魚類群落多樣性[18]。

在Bray-Curtis 相似性的基礎(chǔ)上，運(yùn)用聚類和無度量多維排序（NMDS）分析魚類群落空間分布特征[19,20]，并用相似性分析（ANOSIM）檢驗(yàn)魚類群落空間差異的顯著性（P＜0.05），利用相似性百分比（SIMPER）分析引起群落差異的主要貢獻(xiàn)種和分歧種[21]，多度/生物量比較曲線（ABC 曲線）評(píng)價(jià)群落受干擾的程度[22,23]，典范對(duì)應(yīng)分析（CCA）找出影響魚類群落空間分布的主要環(huán)境因子[17,24]。

使用Primer5.0 軟件進(jìn)行群落多樣性、NMDS、ANOSIM、SIMPER 及ABC 曲線分析與制圖，Canoco5.0 軟件進(jìn)行CCA 分析，Arcgis10.2 及Photoshop 進(jìn)行采樣點(diǎn)繪制。

2 結(jié)果與分析

2.1 物種組成

本次調(diào)查共鑒定魚類4 目9 科31 屬34 種，其中鯉形目2 科23 屬26 種、占總物種數(shù)的76.47%，鲇形目3 科4 屬4 種、占11.76%，鱸形目3 科3 屬3 種、占8.82%，鯡形目僅刀鱭1 種；鯉形目又以鯉科魚類為主，共21 屬24 種，占總物種數(shù)的70.59%（附表1）。本次調(diào)查發(fā)現(xiàn)中國特有魚類8 種、瀕危物種1 種（刀鱭Coilia nasus）、易危2 種（南方擬Pseudohemiculter dispar、月鱧Channa asiatica），其中月鱧為江西省保護(hù)魚類（附表1）。

聚類分析表明，在19.46%的水平上將11 個(gè)采樣點(diǎn)劃分為2 組（圖2）。A 組為采樣點(diǎn)S1～S6，山溪型河流，河道曲折狹窄，共發(fā)現(xiàn)魚類3 目7 科14種；B 組為S7～S11，河道寬闊平直，共采集魚類4 目6 科27 種。

圖2 潦河永安段魚類群落聚類分析圖Fig.2 Cluster analysis of fish communities in the Yong′an section of the Liaohe River

2.2 生態(tài)類型

夏季魚類生態(tài)類型組成主要以雜食、中下層、定居性以及產(chǎn)漂流性和沉性卵的魚類為主（附表1）。A、B 兩組魚類生態(tài)類型組成存在差異：（1）兩組均以雜食性魚類為主，但A 組無植食性的魚類；（2）A 組以底棲魚類為主，B 組以中下層魚類為主；（3）兩組都以定居性魚類為主，但A 組無洄游性魚類；（4）A 組產(chǎn)沉性卵的魚類最多，而B 組產(chǎn)漂流性卵的魚類最多（圖3）。

圖3 A、B 兩組魚類生態(tài)類型組成Fig.3 The composition of fish ecological types in groups A and B

2.3 群落結(jié)構(gòu)特征

2.3.1 優(yōu)勢(shì)種

本次調(diào)查共采集魚類399 尾、重10.91 kg，鯽（Carassius auratus）出現(xiàn)頻率、個(gè)體數(shù)和生物量都最高，為優(yōu)勢(shì)種（表1）。A 組（S1～S6）采到漁獲物151尾、重2.18 kg，魚類平均體質(zhì)量14.41 g，優(yōu)勢(shì)種為寬鰭鱲（Zacco platypus）、側(cè)條光唇魚（Acrossocheilus parallens）、鯽和馬口魚（Opsariichthys bidens）；B 組（S7～S11）采集漁獲物248 尾、重8.73 kg，魚類平均體質(zhì)量為35.20 g，鯽和草魚（Ctenopharyngodon idella）為優(yōu)勢(shì)種（表1）。

表1 A、B 兩組魚類優(yōu)勢(shì)種組成Tab.1 The composition of dominant fish species in group A and group B

2.3.2 群落多樣性

夏季魚類群落的Shannon-Wiener 指數(shù)（H）、Margalef 指數(shù)（D）和Pielou 指數(shù)（J）分別為2.48、3.55和0.7。其中，A 組的Shannon-Wiener 指數(shù)和Margalef 指數(shù)均低于B 組，而Pielou 指數(shù)A 組大于B組（圖4）。經(jīng)SPSS 顯著性檢驗(yàn)，A、B 兩組的Shannon-Wiener 指數(shù)和Pielou 指數(shù)差異顯著（P＜0.05），而Margalef 指數(shù)差異不顯著（P＞0.05）。

圖4 A、B 兩組魚類群落物種多樣性指數(shù)Fig.4 Species diversity indices of fish communities in group A and group B

2.3.3 群落穩(wěn)定性

雖然B 組的生物量曲線一直位于豐度曲線上方，W 統(tǒng)計(jì)值為正值（圖5-c）；但A 組的W 統(tǒng)計(jì)值為負(fù)值，豐度曲線位于生物量曲線的上方（圖5-b）,導(dǎo)致夏季魚類群落的W 統(tǒng)計(jì)值較?。▓D5-a），表明潦河永安段夏季魚類群落受到了一定的干擾。

圖5 潦河永安段（a）、A 組（b）、B 組（c）魚類群落ABC 曲線Fig.5 ABC curve of fish communities in Yong′an section（a）,Group A（b）and Group B（c）of Liaohe River

2.3.4 群落差異性

NMDS 排序結(jié)果顯示，stress＜0.1，夏季魚類群落能較好地分成兩組（圖6）。經(jīng)ANOSIM檢驗(yàn)，兩組之間存在極顯著的差異，群落完全分離、不重疊（R=0.64，P＜0.05）。

圖6 潦河永安段魚類群落空間無度量多維排序圖Fig.6 The nonmetric multidimensional scaling of fish communities in Yong′an section of the Liaohe River

SIMPER 分析顯示：A 組的平均相似性是37.67%，主要貢獻(xiàn)種為寬鰭鱲、側(cè)條光唇魚、鯽和馬口魚，累計(jì)貢獻(xiàn)率為95.67%；B 組的平均相似性是25.22%，主要貢獻(xiàn)種有鯽、草魚、黃顙魚（Pelteobagrus fulvidraco）、達(dá)式鲌（Chanodichthys dabryi）和銀（Xenocypris macrolepis），累計(jì)貢獻(xiàn)率達(dá)90.23%。

A、B 兩組的平均相異性達(dá)88.79%，主要分歧種為寬鰭、側(cè)條光唇魚、鯽、馬口魚、草魚、黃顙魚和達(dá)氏鲌，累計(jì)貢獻(xiàn)率為59.97%（表2）。

表2 A、B 兩組群落內(nèi)相似性及組間相異性貢獻(xiàn)率（＞5%）Tab.2 Contribution rate of similarity and dissimilarity of communities in groups A and B（＞5%）

2.3.5 群落與環(huán)境的關(guān)系

夏季潦河永安段各采樣點(diǎn)環(huán)境因子中，海拔、河寬、氧化還原電位和電導(dǎo)率結(jié)果變化較大，而其他環(huán)境因子變化小。海拔變化范圍為20～99 m，A、B組分別為（71±25.4）m 和（23.6±4.6）m；河寬范圍2.54～470 m，A、B 組分別為（4.39±1.5）m 和（289±166.9）m；氧化還原電位為94.1～201.1 mv，A、B 組分別為（118.1±20）mv 和（177.3±54.4）mv；電導(dǎo)率范圍為35～92 us/cm，A、B 組分別為（58.5±16.9）us/cm 和（78±15.4）us/cm（表3）。

表3 潦河永安段各采樣點(diǎn)環(huán)境因子Tab.3 Environmental factors in each sampling site in Yong′an section of the Liaohe River

夏季魚類群落與環(huán)境因子的CCA 分析顯示，前兩個(gè)排序軸的特征值分別為0.6706 和0.5291，前四個(gè)排序軸對(duì)所有變量的解釋度可達(dá)64.70%（表4）。在本次調(diào)查中，取前兩個(gè)軸做CCA 排序，發(fā)現(xiàn)海拔（Alt）和河寬（Width）對(duì)物種具有顯著性的影響（P＜0.05），兩者對(duì)群落差異的解釋度為37.3%，可認(rèn)為是驅(qū)動(dòng)群落差異的主要因子。其中海拔與第一排序軸呈最大正相關(guān)關(guān)系，河寬與第一排序軸呈最大負(fù)相關(guān)關(guān)系。除河寬、酸堿度和氧化還原電位與第二排序軸呈負(fù)相關(guān)關(guān)系外，其他環(huán)境因子與第二排序軸呈正相關(guān)關(guān)系（圖7）。

表4 夏季潦河永安段排序軸特征值及其與環(huán)境變量的相關(guān)系數(shù)Tab.4 The eigenvalue of axis and their correlation coefficient with environmental variables in Yong′an Section of Liaohe River in summer

圖7 潦河永安段夏季魚類群落和環(huán)境因子的關(guān)系Fig.7 The relationship between fish communities and environmental factors in the Yong′ an section of the Liaohe River in summer

3 討論

本實(shí)驗(yàn)中，物種的空間分布顯示：A 組的魚類物種數(shù)、個(gè)體數(shù)、生物量都低于B 組，主要原因是受河流大小和空間位置的影響[25,26]。A 組主要以山區(qū)為主，河道窄，水流急，河流中營養(yǎng)物質(zhì)少；B 組以寬淺且水流平緩的河、湖為主，水流較緩慢，有利于水中浮游生物的繁殖，為魚類提供充足的餌料。從魚類的生態(tài)類型看，該流域以中下層、雜食性的魚類為主，而草食性的魚類最少，說明了該水域中腐殖質(zhì)和有機(jī)的碎屑較多，而植物性餌料缺乏。魚類的生態(tài)類型也能從側(cè)面反映出水域生態(tài)環(huán)境[27,28]，即在環(huán)境質(zhì)量下降的情況下，食性廣的魚類更易適應(yīng)環(huán)境[29]。本研究中，因雜食性物種占比最高，污染耐受性高，如鯽、南方擬等，反映出該流域的水質(zhì)量有所下降。在本次調(diào)查中，瀕危物種刀鱭僅捕撈到4尾，江西省級(jí)保護(hù)魚類月鱧僅2 尾，僅在1 個(gè)采樣點(diǎn)捕撈到，同時(shí)鮈亞科魚類也捕撈較少。魚類優(yōu)勢(shì)種群在空間組成上明顯不同，鯽為該流域的主要優(yōu)勢(shì)種群，主要是由于其食性雜，對(duì)食物的選擇廣泛，能適應(yīng)各種水體。A 組的優(yōu)勢(shì)種主要為山溪性魚類，適宜生活在流速快，砂石底質(zhì)的河流中；而B 組的優(yōu)勢(shì)種為草魚，它主要生活在水流緩慢、水體開闊的河段；其次A 組以電捕為主，對(duì)魚類沒有選擇性，B 組以絲網(wǎng)捕撈為主，對(duì)魚的大小具有選擇性，因此優(yōu)勢(shì)種的變化受捕撈方式、河寬和魚類棲息環(huán)境有關(guān)。

群落中物種豐富度高，物種分布均勻，就容易形成一個(gè)較復(fù)雜的群落，對(duì)環(huán)境的干擾有較大的緩沖能力。多樣性指數(shù)測(cè)量群落結(jié)構(gòu)和功能的復(fù)雜性結(jié)果顯示[30]，夏季魚類群落的Shannon-Wiener 指數(shù)（H）、Margalef 指數(shù)（D）和Pielou 指數(shù)（J）分別為2.48、3.55 和0.7，表明該流域魚類物種分布較均勻，種類較豐富，同時(shí)1≤H≤3，反映出水體受到中等干擾[31]。該結(jié)果與ABC 曲線結(jié)果一致。ABC 曲線結(jié)果顯示，以A 組為代表的（采樣點(diǎn)S1～S6）受到了中等程度的干擾，可能是由于氣候原因，如干旱和洪水強(qiáng)烈的影響魚類棲息地[32]。A 組（S1～S6）為山溪性溪流，水流交換速度快。采樣時(shí)間在暴雨集中后的一段時(shí)間，洪水對(duì)A 組的魚類群落結(jié)構(gòu)影響較大，所以呈現(xiàn)出中等干擾的狀態(tài)；B 組（S7～S11），水面寬廣，水流交換速度慢，水體自凈能力強(qiáng)，群落受到的干擾較小。其次，優(yōu)勢(shì)種個(gè)體大小也會(huì)影響群落穩(wěn)定性[33]。A 組中的優(yōu)勢(shì)種以小型魚類為主，其優(yōu)勢(shì)種（寬鰭鱲、側(cè)條光唇魚、馬口魚和鯽）的重量百分比和數(shù)量百分比分別為81.46%和69.27%，豐度比生物量更具優(yōu)勢(shì)地位,生物量曲線位于豐度曲線下方。A 組多位于支流，捕魚方式以電捕為主，但其位置偏遠(yuǎn)、分散，管理難度大，也可能導(dǎo)致該河段魚類恢復(fù)緩慢。因此魚類群落穩(wěn)定性受水文、優(yōu)勢(shì)種和人為管理的影響。但W 統(tǒng)計(jì)值發(fā)現(xiàn)，都趨近于0，反映出該流域魚類群落結(jié)構(gòu)受到一定程度的干擾。

Cluster 和NMDS 的分析結(jié)果表明，夏季魚類群落的排序結(jié)果較好。通常，將Cluster 和NMDS分析的結(jié)果相互對(duì)照，可以得出更具生態(tài)意義的結(jié)論[34]。在該流域中，11 個(gè)采樣點(diǎn)被分為2 組，每組的典型種為各組的優(yōu)勢(shì)種或重要種。所以，優(yōu)勢(shì)種數(shù)量的相對(duì)變化也是造成群落空間差異的重要因素[35]。從CCA 圖可知，該流域魚類群落結(jié)構(gòu)的變化與環(huán)境中的海拔和河寬密切相關(guān)。在較寬的水域中，水較深，能容納更多的魚類，更適合大型魚類生長，且離岸較遠(yuǎn)，受干擾小，因此魚類豐富度更高，種類更多[36,37]。如草魚（Ctenopharyngodon idella）、大鰭鱊（Acheilognathus macropterus）及黃顙魚（Tachysurus fulvidracog）與河寬呈正相關(guān)，主要分布在河道寬的水域。其中泥鰍（Misgurnus anguillicaudatus）、鰻尾（Liobagrus anguillicauda）和盎堂擬鲿（Pseudobagrus ondon）與海拔呈正相關(guān)，而與河寬呈負(fù)相關(guān)，主要分布在海拔高、水流較急的區(qū)域，如S1 和S2 山溪型河段。潦河永安段夏季魚類群落空間的分布特征受很多因素影響，其他因素及具體影響機(jī)制和群落結(jié)構(gòu)的時(shí)間變化還需進(jìn)一步研究。

總結(jié)

潦河永安段主要以鯉形目鯉科魚類為主，鯡形目種類最少，這與修河流域魚類組成較一致，但不合理捕撈導(dǎo)致魚類物種數(shù)和資源量不斷下降；生態(tài)類型以雜食性魚類居多，這與水質(zhì)下降息息相關(guān)；歷史優(yōu)勢(shì)種群由之前的類、黃顙魚類、鯉鯽類轉(zhuǎn)變?yōu)轳R口魚、寬鰭鱲、鯽和鱊類等，魚類個(gè)體小型化[9]。因此，應(yīng)加強(qiáng)漁政執(zhí)法力度，制定相關(guān)的漁業(yè)資源法律法規(guī)，加大宣傳和普及魚類資源的重要性；其次，污水應(yīng)處理后再排放；最后，進(jìn)行魚類人工增殖放流活動(dòng)，增加魚類種群結(jié)構(gòu)中低齡和幼齡魚類的數(shù)量，擴(kuò)大種群規(guī)模，從根本上解決天然魚類資源不足的問題。