把邊江流域大型底棲動(dòng)物群落與水質(zhì)的關(guān)系

和克儉，張 琦，丁 佼

(1.云南大學(xué) 地球科學(xué)學(xué)院，云南 昆明 650500；2.云南省水利廳，云南 昆明 650021)

0 引言

大型底棲動(dòng)物是河流生態(tài)系統(tǒng)中的主要?jiǎng)游锶后w，是維持河流生態(tài)系統(tǒng)生物群落結(jié)構(gòu)和功能穩(wěn)定的重要因素[1，2]。底棲動(dòng)物在河流生態(tài)系統(tǒng)中廣泛分布，其移動(dòng)緩慢、易于采集，同時(shí)底棲動(dòng)物對(duì)水質(zhì)敏感，具有較好的水質(zhì)環(huán)境指示作用，是河流生態(tài)和水質(zhì)研究的熱點(diǎn)[1-5]。底棲動(dòng)物群落結(jié)構(gòu)和分布狀況與其生境水質(zhì)息息相關(guān)，水體總氮(TN)、總磷(TP)、無(wú)機(jī)氮(DIN)都會(huì)對(duì)底棲動(dòng)物的群落結(jié)構(gòu)與分布狀況產(chǎn)生影響[6-10]。因此，底棲動(dòng)物是目前應(yīng)用廣泛的水生生物監(jiān)測(cè)指示生物類群，厘清底棲動(dòng)物群落與水質(zhì)的關(guān)系有助于科學(xué)開(kāi)展河流生態(tài)系統(tǒng)健康評(píng)價(jià)。

目前國(guó)內(nèi)大型底棲動(dòng)物群落與水質(zhì)的相關(guān)研究主要集中在中國(guó)東部地區(qū)，如渤海灣地區(qū)[9]、長(zhǎng)江中下游地區(qū)[11]、濟(jì)南市[12]、上海市[13]、遼河流域[7]等，而對(duì)西南地區(qū)河流目前主要集中在瀾滄江等河流。紅河是中國(guó)西南地區(qū)重要的國(guó)際河流，對(duì)其上游地區(qū)底棲動(dòng)物調(diào)查數(shù)據(jù)、底棲動(dòng)物群落與水質(zhì)的關(guān)系研究極為缺乏。因此，以紅河流域上游——把邊江流域?yàn)檠芯繀^(qū)域，調(diào)查流域大型底棲動(dòng)物群落及其生境水質(zhì)狀況，分析識(shí)別影響把邊江流域大型底棲動(dòng)物群落的主要水質(zhì)因子。以補(bǔ)充完善相關(guān)研究，為科學(xué)評(píng)價(jià)目標(biāo)流域水生態(tài)系統(tǒng)，加強(qiáng)水生態(tài)系統(tǒng)的保護(hù)與管理，促進(jìn)水生態(tài)系統(tǒng)和多樣性的保護(hù)與建設(shè)提供數(shù)據(jù)支持。

1 材料與方法

1.1 研究區(qū)域概況

把邊江是紅河流域最大支流——李仙江的上游，發(fā)源于大理州南澗縣，自西北向西南流經(jīng)景東縣、鎮(zhèn)沅縣，過(guò)墨江縣和寧洱縣后與阿墨江匯合成李仙江。把邊江干流河長(zhǎng)約382 km，集水區(qū)面積約8 948 km2[14](圖1)，流域地處橫斷山脈南段，地勢(shì)西北高東南低，地形以山地為主[14]。流域氣候以亞熱帶季風(fēng)氣候?yàn)橹?，干濕分明，雨季?～10月，旱季為11～4月，年均降水量約為1 225.4 mm，年平均氣溫為19.0℃[14]。主要植被為亞熱帶常綠闊葉林，主要土壤類型為紅壤、黃壤和黃棕壤[15]。

1.2 大型底棲動(dòng)物群落及其生境水質(zhì)狀況調(diào)查

于2019年8月～9月開(kāi)展把邊江流域大型底棲動(dòng)物群落及其生境水質(zhì)狀況調(diào)查。

流域內(nèi)共布設(shè)31個(gè)采樣點(diǎn)(圖1)。采樣點(diǎn)布設(shè)盡可能均勻、全面地覆蓋流域內(nèi)的干流和支流，同時(shí)考慮樣點(diǎn)的交通可達(dá)性、安全性和生境典型性。

圖1 把邊江流域區(qū)位和采樣點(diǎn)圖

在穩(wěn)定流量狀況(48小時(shí)內(nèi)降雨小于10 mm)下采集生境水質(zhì)樣品[16]。水質(zhì)樣品采集參照《水和廢水監(jiān)測(cè)分析方法(第四版)》標(biāo)準(zhǔn)(國(guó)家環(huán)境保護(hù)總局2002)，在采樣點(diǎn)的河流中央水深0.5m處使用1L有機(jī)玻璃采水器采集水樣。實(shí)行對(duì)水溫(TEM)、電導(dǎo)率(EC)、酸堿度(pH)和溶解氧(DO)的現(xiàn)場(chǎng)測(cè)定，其余指標(biāo)采集水樣后在實(shí)驗(yàn)室利用進(jìn)行測(cè)定。水質(zhì)因子符號(hào)、單位，測(cè)量?jī)x器/測(cè)定方法詳見(jiàn)表1。

表1 水質(zhì)指標(biāo)說(shuō)明Tab.1 Description of water quality indicators

采用“復(fù)合生境采樣法”在采集生境水質(zhì)樣品同時(shí)進(jìn)行大型底棲動(dòng)物采集[17]。在每個(gè)樣點(diǎn)典型河段內(nèi)，使用D形網(wǎng)(40目紗，0.3 m寬)采集河流底質(zhì)約2 kg，采集長(zhǎng)度為3～5 m。采集的河流底質(zhì)在現(xiàn)場(chǎng)使用60目網(wǎng)篩進(jìn)行篩洗，將附著在石塊上的底棲動(dòng)物洗刷入網(wǎng)。篩洗后的底質(zhì)放入塑封袋中，置于保溫箱中低溫保存，帶回駐地進(jìn)行挑選。底質(zhì)中的底棲動(dòng)物當(dāng)日挑取并保存在75%酒精溶液中，用于進(jìn)行室內(nèi)鑒定。物種類別鑒定至科級(jí)，鑒定時(shí)對(duì)物種進(jìn)行計(jì)數(shù)，計(jì)算物種密度數(shù)據(jù)，物種密度定義為單位重量河流底質(zhì)所含底棲動(dòng)物個(gè)體數(shù)量(ind/kg)。

根據(jù)攝食對(duì)象和攝食方式分為5個(gè)攝食功能類群(Functional Feeding Groups，F(xiàn)FG)[18-23]：

(1)捕食者(Predators，PR)，直接吞食或刺食其他底棲動(dòng)物。

(2)撕食者(Shredders，SH)，食用粗顆粒有機(jī)質(zhì)。

(3)直接收集者(Gather collectors，GC)，直接收集者食物以水底的細(xì)顆粒有機(jī)物為主。

(4)濾食收集者(Filter collectors，F(xiàn)C)，主要以水流中細(xì)顆粒有機(jī)物為食。

(5)刮食者(Scrapers，SC)，以固著在底質(zhì)上的底棲動(dòng)物為食。

1.3 數(shù)據(jù)分析與處理

1.3.1 指數(shù)計(jì)算方法

在表1中，Shannon-Wiener多樣性指數(shù)的計(jì)算公式應(yīng)為：多樣性指數(shù)選用了Shannon-Weaver指數(shù)(H′)[24]和Pielou均勻度指數(shù)(J)[25]，具體計(jì)算公式如下：

(1)

J=H′/H′max

(2)

式中：S為物種數(shù)；Pi=ni/N，ni為第i種個(gè)體數(shù)，N為樣本總個(gè)體數(shù)。H′越大，物種多樣性越高。H′為實(shí)測(cè)Shannon-Weaver指數(shù)，H′max是群落中最大Shannon-Weaver指數(shù)值，反映了各物種數(shù)目分配的均勻程度。

1.3.2 數(shù)據(jù)統(tǒng)計(jì)分析

(1)相關(guān)分析：用于探究水質(zhì)因子與大型底棲動(dòng)物指數(shù)之間的相關(guān)關(guān)系(SPSS)。

(2)群落結(jié)構(gòu)分析：首先對(duì)物種—物種密度矩陣進(jìn)行四次方根轉(zhuǎn)換，以平衡稀有種在大型底棲動(dòng)物群落中的影響[26]；基于物種矩陣，使用PRIMER 6軟件計(jì)算相似性矩陣，利用等級(jí)聚類功能(CLUSTER)對(duì)標(biāo)準(zhǔn)化處理后的物種矩陣進(jìn)行聚類分析，對(duì)聚類結(jié)果進(jìn)行物種組成差異檢驗(yàn)(ANOSIM)，并進(jìn)行優(yōu)勢(shì)種分析(SIMPER)。

(3)方差分析：根據(jù)群落結(jié)構(gòu)分組結(jié)果，利用SPSS中單因素方差分析(ANOVA)功能分析不同群落生境水質(zhì)的差異，并進(jìn)行Duncan多重比較(post hoc test)分析，顯著性水平為P<0.05，判斷造成大型底棲動(dòng)物種群結(jié)構(gòu)差異的主要水質(zhì)因子。

(4)多元統(tǒng)計(jì)分析：使用CANOCO 5軟件進(jìn)行多元統(tǒng)計(jì)分析。首先使用去趨勢(shì)對(duì)應(yīng)分析(Detrended correspondence analysis，DCA)對(duì)四次方根轉(zhuǎn)換后的物種—密度矩陣進(jìn)行排序，DCA分析中前兩個(gè)排序軸的梯度值分別為3.63和3.55，因此使用典范對(duì)應(yīng)分析(Canonical correspondence analysis，CCA)分析大型底棲動(dòng)物群落、物種與水質(zhì)的關(guān)系[12]，使用冗余分析(Redundancy analysis，RDA)分析大型底棲動(dòng)物攝食功能類群與水質(zhì)因子的關(guān)系[27]。

2 結(jié)果與分析

2.1 大型底棲動(dòng)物種類組成

研究區(qū)31個(gè)樣點(diǎn)共調(diào)查到大型底棲動(dòng)物共計(jì)25個(gè)分類單元(鑒定至科級(jí)水平)。包括3門(節(jié)肢動(dòng)物門、軟體動(dòng)物門、環(huán)節(jié)動(dòng)物門)6綱(昆蟲(chóng)綱、軟甲綱、瓣鰓綱、腹足綱、蛭綱、寡毛綱)14目25科。其中水生昆蟲(chóng)18個(gè)分類單元，占總分類單元數(shù)的72%，是構(gòu)成把邊江流域底棲動(dòng)物群落的主要類群(表2)。

表2 把邊江流域大型底棲動(dòng)物群落組成Tab.2 Macrobenthos community composition in Babian River Basin 個(gè)

根據(jù)攝食對(duì)象和攝食方式，研究區(qū)底棲動(dòng)物可以分為5個(gè)攝食功能類群(Functional Feeding Groups，F(xiàn)FG)[18-23]。其中捕食者(PR)7個(gè)分類單元，占總分類單元的28%；刮食者(SC)6個(gè)分類單元，占總分類單元的24%；直接收集者(GC)6個(gè)分類單元，占總分類單元的24%；濾食收集者(FC)5個(gè)分類單元，占總分類單元的20%；撕食者(SH)1個(gè)分類單元，占總分類單元的4%。把邊江流域大型底棲動(dòng)物攝食功能類群以捕食者(PR)、刮食者(SC)、直接收集者(GC)和濾食收集者(FC)為主。

2.2 大型底棲動(dòng)物群落特征

2.2.1 聚類分析和群落差異分析

聚類分析結(jié)果表明，在10%的相似度水平上，可將把邊江流域大型底棲動(dòng)物劃分成4個(gè)群落類型(圖2)。在此聚類結(jié)果基礎(chǔ)上進(jìn)行相似性分析(ANOSIM)，發(fā)現(xiàn)組A(y0、y1)與組B(y8、y13、y18)物種組成差異不顯著(R=1，P>0.01)。據(jù)此對(duì)聚類結(jié)果進(jìn)行適當(dāng)調(diào)整(將組A與組B合并為組I)，對(duì)調(diào)整后的相似性分析(ANOSIM)結(jié)果顯示，組I、組II和組III之間物種組成差異顯著(R=0.544，P<0.01)(表3)。

圖2 聚類分析結(jié)果

表3 把邊江流域大型底棲動(dòng)物群落分組及組間相似性分析Tab.3 Grouping and ANOSIM of macrobenthosin Babian River Basin

2.2.2 優(yōu)勢(shì)種分析

基于調(diào)整后的物種聚類分組和物種-物種密度矩陣，進(jìn)行相似性百分比分析(SIMPER)，得到每組內(nèi)累計(jì)貢獻(xiàn)率>90%的物種，即優(yōu)勢(shì)種。結(jié)合現(xiàn)場(chǎng)調(diào)查結(jié)果和樣點(diǎn)空間分布(圖3)分析發(fā)現(xiàn)：組Ⅰ物種主要分布在上游和中游生境破壞較嚴(yán)重的干流河段，優(yōu)勢(shì)種為長(zhǎng)臂蝦科(Palaemonidae)、四節(jié)蜉科(Baetidae)、紋石蛾科(Hydopsychidae)和石蛭科(Herpobdellidae)；組Ⅱ物種主要分布在上游和下游城鎮(zhèn)附近的干流河段。優(yōu)勢(shì)種為紋石蛾科(Hydopsychidae)、泥甲科(Dryopidae)和顫蚓科(Tubificidae)；組Ⅲ物種主要分布在生境較好的支流河段，優(yōu)勢(shì)種為扁蜉科(Heptageniidae)、紋石蛾科(Hydopsychidae)、石蠶科(Phryganea Larva)和四節(jié)蜉科(Baetidae)。

圖3 把邊江流域大型底棲動(dòng)物群落分布

2.3 多樣性特征和密度特征

由表4可知，各采樣點(diǎn)之間把邊江流域大型底棲動(dòng)物的Shannon-Weaver指數(shù)(H′)和Pielou均勻度指數(shù)(J)差異較小，H′平均值普遍不大于1.5。根據(jù)群落分類，H′和J在流域內(nèi)規(guī)律均為：組Ⅰ<組Ⅱ<組Ⅲ，與分類單元數(shù)分布規(guī)律一致。

表4 把邊江流域大型底棲動(dòng)物及攝食功能類群Shannon-Weaver指數(shù)(H′)、Pielou均勻度指數(shù)(J)和密度平均值Tab.4 The average of Shannon-weaver Index(H′)，Pielou evenness Index(J) and density(ind/kg)of macrobenthos and functional feeding groups in the Babian River Basin

把邊江流域大型底棲動(dòng)物的密度平均值：組Ⅱ<組Ⅰ<組Ⅲ；各功能群密度均值，捕食者(PR)密度：組Ⅲ<組Ⅱ<組Ⅰ；直接收集者(GC)密度：組Ⅱ<組Ⅲ<組Ⅰ；刮食者(SC)密度：組Ⅰ<組Ⅱ<組Ⅲ；濾食收集者(FC)密度：組Ⅱ<組Ⅰ<組Ⅲ；撕食者(SH)密度：組Ⅰ=組Ⅱ<組Ⅲ。

2.4 水質(zhì)基本特征

把邊江流域氮磷含量較高，其中氨氮(NH3-N)平均值為1.493 mg/L，接近IV類水體1.5 mg/L的閾值?？偭?TP)平均值0.564 mg/L，超過(guò)V類水體0.4 mg/L的閾值，水體污染以氮磷污染為主；高錳酸鹽指數(shù)(CODMn)平均值為5.636 mg/L，為III類水體(6 mg/L)，有機(jī)污染程度一般；溶解氧較高，水質(zhì)偏堿性(表5)。其中氨氮(NH3-N)、硝氮(NO3-N)、葉綠素a(Chl-a)和總磷(TP)變異系數(shù)較大(>50%)，空間差異明顯，其較高值主要出現(xiàn)在人為活動(dòng)較為頻繁的流域北部以及干流部分，同時(shí)這些區(qū)域溶解氧值較低，水質(zhì)較差(圖4)。其它水質(zhì)因子變異系數(shù)較小(< 50%)，空間差異不明顯。

表5 水質(zhì)指標(biāo)基本統(tǒng)計(jì)特征Tab.5 Basic statistical characteristics of water quality indicators

2.5 大型底棲動(dòng)物群落結(jié)構(gòu)與水質(zhì)因子的關(guān)系

2.5.1 大型底棲動(dòng)物群落和物種與水質(zhì)因子的關(guān)系

(1)不同群落生境水質(zhì)差異

單因素方差分析結(jié)果(表6)表明，不同群落生境水體的NH3-N、NO3-N、TN和TP存在顯著差異(P<0.05)，其他水質(zhì)因子在不同的群落生境間無(wú)顯著差異。其中組I的NH3-N顯著高于組Ⅲ；組I的NO3-N顯著高于組II和組III；組I和組II的TN和TP均顯著高于組III。

表6 大型底棲動(dòng)物群落結(jié)構(gòu)與水質(zhì)因子單因素方差分析

(2)CCA排序分析

CCA排序前兩軸共解釋41.55%的物種與水質(zhì)關(guān)系，第一軸特征值為0.356，解釋率為21.98%；第二軸特征值為0.317，解釋率為19.57%。第一軸與pH、DO、NO3-N、EC、TN、TP、Chl-a、NH3-N、CODMn相關(guān)性較強(qiáng)，第一軸為污染軸。

從CCA分析的樣點(diǎn)—環(huán)境排序圖(圖5)可看出，群落Ⅰ(星號(hào))主要集中在第一和第四象限，主要分布在NH3-N、NO3-N、TN、TP較高值區(qū)域，群落Ⅱ(三角形)主要集中在第四象限，主要分布在NH3-N、TN、TP較高值區(qū)域；群落Ⅲ(圓形)主要集中在第二和第三象限，主要分布在DO較高，NH3-N、TN、TP較低值區(qū)域。

圖5 CCA樣點(diǎn)—環(huán)境排序圖

從CCA分析的物種—環(huán)境排序圖(圖6)可看出，蟌科(Coenagrionidac)、搖蚊科(Chironnmodae)、觹螺科(Hydrobiidae)和石蛭科(Herpobdellidae)等一般耐污種以及顫蚓科(Tubificidae)等耐污種多分布在污染軸(第一軸)的正方向。扁蜉科(Heptageniidae)、等蜉科(Isonychiidae)、蓋蝽科(Aplelocheiridae)等生活在清潔水域的物種多分布在污染軸負(fù)方向。沿污染軸，隨水體污染程度的增加，EPT指示種的生物量逐漸減少。

a b c d

e f g h

i j

圖6 CCA物種—環(huán)境排序圖

2.5.2 大型底棲動(dòng)物攝食功能類群與水質(zhì)因子的關(guān)系

把邊江流域大型底棲動(dòng)物攝食功能類群與水質(zhì)因子RDA分析結(jié)果(圖7)顯示，第一和第二排序軸共解釋73.77%的功能類群與水質(zhì)因子的關(guān)系，第一軸特征值為0.237，解釋率為51.67%；第二軸特征值為0.101，解釋率為22.10%。第一軸與TN、NH3-N、NO3-N、TP、EC、Chl-a正相關(guān)，與DO和pH負(fù)相關(guān)。NO3-N、TN、NH3-N(P值分別為0.002、0.002、0.02)單獨(dú)解釋率分別為15.8%、12.1%、8.9%，NO3-N、TN、NH3-N是影響攝食功能類群群落結(jié)構(gòu)的主要水質(zhì)因子。

圖7 RDA物種—環(huán)境排序圖

撕食者(SH)與CODMn負(fù)相關(guān)，直接收集者(GC)和捕食者(PR)與CODMn正相關(guān)，說(shuō)明這3種攝食功能群群落分布與水體有機(jī)污染相關(guān)。刮食者(SC)和濾食收集者(FC)與污染軸呈負(fù)相關(guān)，說(shuō)明這兩種攝食功能群群落偏好分布在清潔水域。

2.5.3 物種、攝食功能類群多樣性與水質(zhì)因子的關(guān)系

大型底棲動(dòng)物物種多樣性、攝食功能類群多樣性與生境水質(zhì)因子的Spearman秩相關(guān)分析結(jié)果(表7)表明：物種多樣性H′與NH3-N、NO3-N、TN、TP呈顯著負(fù)相關(guān)；攝食功能類群多樣性與NO3-N、TN呈顯著負(fù)相關(guān)。而物種均勻度、攝食功能類群均勻度與所有水質(zhì)指標(biāo)均沒(méi)有顯著的相關(guān)關(guān)系。

表7 大型底棲動(dòng)物及其攝食功能類群多樣性、均勻度與水質(zhì)因子相關(guān)性分析Tab.7 Correlation analysis of diversity，evenness and water quality factors of macrobenthos and their functional feeding groups

3 討論與結(jié)論

3.1 把邊江流域大型底棲動(dòng)物群落與水質(zhì)因子的關(guān)系

把邊江流域大型底棲動(dòng)物群落分異明顯，人類活動(dòng)和自然地理環(huán)境的差異可能是造成研究區(qū)大型底棲動(dòng)物群落分異明顯的原因之一[28]。群落組Ⅰ中優(yōu)勢(shì)種為紋石蛾科(Hydopsychidae)、長(zhǎng)臂蝦科(Palaemonidae)、四節(jié)蜉科(Baetidae)和石蛭科(Herpobdellidae)，均為一般耐污種[29]，組內(nèi)貢獻(xiàn)率最高的紋石蛾科相較于其他3組平均密度較低(表3)。這可能與組Ⅰ分布區(qū)域生境破壞比較嚴(yán)重，流域內(nèi)有建橋、修路等工程建設(shè)，或處于城市附近干流，人口聚集，人類活動(dòng)干擾明顯有關(guān)。而組Ⅱ優(yōu)勢(shì)種紋石蛾科(Hydopsychidae)和泥甲科(Dryopidae)為一般耐污種，顫蚓科(Tubificidae)常屬于污染嚴(yán)重地區(qū)的優(yōu)勢(shì)種[30]。這可能與組Ⅱ主要分布在農(nóng)村城鎮(zhèn)較為集中的干流河段，生活污染較為嚴(yán)重有關(guān)。組Ⅲ中優(yōu)勢(shì)種占比最高的扁蜉科(Heptageniidae)為清潔物種，組Ⅲ主要分布在生境完整，水質(zhì)較好的支流。同時(shí)，EPT在三組中的分類單元數(shù)為組Ⅰ<組Ⅱ<組Ⅲ，從側(cè)面體現(xiàn)不同群落生境水質(zhì)的差異。

方差分析結(jié)果表明把邊江流域不同大型底棲動(dòng)物群落生境水體的NH3-N、NO3-N、TN和TP存在顯著差異(P<0.05)。以往研究表明氮磷是影響大型底棲動(dòng)物群落的重要因子[31]。在CCA排序圖中，隨著NH3-N、NO3-N、TN和TP增加，底棲EPT物種逐漸減少。EPT指示種(蜉蝣目、襀翅目、毛翅目)在一定程度上體現(xiàn)了水體的清潔程度[33]。因此，NH3-N、NO3-N、TN和TP是引起把邊江流域大型底棲動(dòng)物群落分異的主要水質(zhì)因子[32]。此外，研究區(qū)EPT指示種對(duì)水質(zhì)污染比較敏感，說(shuō)明EPT指示種可以作為該區(qū)域水生態(tài)健康評(píng)價(jià)的指標(biāo)。

攝食功能類群與水質(zhì)關(guān)系的分析結(jié)果表明：撕食者(SH)密度與CODMn負(fù)相關(guān)，直接收集者(GC)密度和捕食者(PR)密度與CODMn正相關(guān)，說(shuō)明這3種攝食功能類群群落分布與水體有機(jī)污染相關(guān)。刮食者(SC)密度和濾食收集者(FC)密度與污染軸(第一軸)呈負(fù)相關(guān)，說(shuō)明這兩種攝食功能類群偏好分布在清潔水域。結(jié)合研究區(qū)底棲動(dòng)物功能類群的密度(表6)，平均密度值較高的GC可以作為評(píng)價(jià)有機(jī)污染的底棲因子，而平均密度值較高的SC和FC可以作為評(píng)價(jià)氮磷、溶解氧等水質(zhì)狀況的底棲因子。

3.2 把邊江流域大型底棲動(dòng)物多樣性與水質(zhì)因子的關(guān)系

由Spearman秩相關(guān)分析(表6)可知，NH3-N、NO3-N、TN和TP與大型底棲動(dòng)物物種多樣性指數(shù)(H’)呈顯著負(fù)相關(guān)。從 NH3-N、NO3-N、TN和TP較高值在流域主要分布在干流及靠近城鎮(zhèn)河段(圖4)。生活垃圾污染、生活污水排放和農(nóng)業(yè)面源污染是氮磷排放的主要來(lái)源，由此引起的營(yíng)養(yǎng)鹽濃度增加造成了其大型底棲動(dòng)物多樣性降低[10]。因此，建議目標(biāo)流域進(jìn)一步加強(qiáng)水污染防治和管理，強(qiáng)化城鎮(zhèn)生活污染治理，推進(jìn)農(nóng)業(yè)農(nóng)村污染防治。重點(diǎn)針對(duì)NH3-N、NO3-N、TN和TP在流域內(nèi)的排放，加強(qiáng)水污染排放的監(jiān)督管理，落實(shí)水環(huán)境治理行動(dòng)計(jì)劃，保護(hù)和改善流域水質(zhì)，推動(dòng)流域大型底棲動(dòng)物多樣性的保護(hù)和恢復(fù)。

3.3 結(jié)論

(1)把邊江流域大型底棲動(dòng)物群落分異明顯，不同群落生境水體的NH3-N、NO3-N、TN和TP存在顯著的差異。且NH3-N、NO3-N、TN和TP與大型底棲動(dòng)物物種多樣性呈顯著負(fù)相關(guān)。NH3-N、NO3-N、TN和TP是引起把邊江流域大型底棲動(dòng)物群落分異和多樣性降低的主要水質(zhì)因子。

(2)研究區(qū)EPT指示種(蜉蝣目、襀翅目、毛翅目)對(duì)水質(zhì)污染比較敏感，EPT指示種可以作為該區(qū)域水生態(tài)健康評(píng)價(jià)的指標(biāo)。

(3)撕食者(SH)密度與CODMn負(fù)相關(guān)，直接收集者(GC)密度和捕食者(PR)密度與CODMn正相關(guān)；而刮食者(SC)密度和濾食收集者(FC)密度與氮磷污染軸呈負(fù)相關(guān)。GC可以用于評(píng)價(jià)目標(biāo)流域有機(jī)污染狀況，而SC和FC可以用于評(píng)價(jià)目標(biāo)流域氮磷、溶解氧等水質(zhì)狀況。

致謝：感謝云南大學(xué)黃曉霞副教授和劉虹、北京師范大學(xué)劉琦在野外調(diào)研和室內(nèi)實(shí)驗(yàn)中給予的幫助。