王纖纖 劉樂樂 楊學(xué)芬 楊瑞斌* 劉海平

(1.華中農(nóng)業(yè)大學(xué)水產(chǎn)學(xué)院, 武漢 430070; 2.西藏自治區(qū)農(nóng)牧科學(xué)院水產(chǎn)科學(xué)研究所, 拉薩 850002; 3.西藏自治區(qū)動物疫病預(yù)防控制中心(畜牧總站), 拉薩 850000)

周叢藻類(Periphyton), 又名著生藻類、底棲藻類等[1,2], 是水生態(tài)系統(tǒng)中主要的生物類群之一, 營自養(yǎng)生活, 其通過光合作用將太陽能轉(zhuǎn)化為化學(xué)能儲存起來[3], 并產(chǎn)生和釋放氧氣, 推動水生態(tài)系統(tǒng)中的物質(zhì)循環(huán)、能量流動和信息傳遞, 是水生態(tài)系統(tǒng)中食物鏈的開端。周叢藻類擁有較高的初級生產(chǎn)力, 其種類豐富, 附著生長, 生活周期短, 對水環(huán)境變化敏感[4,5], 能吸附水體中的污染物, 凈化水質(zhì)[6],且其種類組成和數(shù)量多寡, 又會隨水環(huán)境中水質(zhì)理化因子和污染物的變動和而發(fā)生相應(yīng)的改變, 加之,周叢藻類的附著生活使其易于采集, 所以, 在急流水域周叢藻類是重要的水質(zhì)監(jiān)測指示物種, 常用于水生態(tài)系統(tǒng)的健康評價中[7—9]。

河流作為生命的搖籃, 是自然生態(tài)系統(tǒng)中一種主要的形態(tài), 隨著人類對自然干預(yù)力度的增加, 許多河流面臨天然水生態(tài)系統(tǒng)發(fā)生改變[10]、水文情勢發(fā)生變化[11]、河流連通性改變[12]、水環(huán)境污染[13]、河道萎縮[14,15]和河道斷流[16,17]等問題, 這些問題直接作用于水環(huán)境中的動植物, 使得水體中的魚、蝦、蟹、螺、浮游動植物、底棲動植物等的群落結(jié)構(gòu)、生物多樣性和種群的穩(wěn)定性受到威脅[18], 降低河流的健康狀態(tài)和生態(tài)服務(wù)功能。故對河流水生態(tài)系統(tǒng)的健康狀態(tài)進(jìn)行科學(xué)的評估, 逐漸成為了國內(nèi)外研究者關(guān)注的焦點, 但若僅依靠水質(zhì)理化因子對河流進(jìn)行健康評價, 局限性較多, 而依據(jù)生物完整性指數(shù)(Index of Biological Integrity, IBI)體系, 則可從水生態(tài)系統(tǒng)結(jié)構(gòu)與功能的角度, 分析河流的健康現(xiàn)狀及其原因[19], 在李喆等[7]對松花江、劉園園等[8]對三峽庫區(qū)、張遠(yuǎn)等[19]對遼河流域等的研究中, 均運用了生物完整性指數(shù)對水域進(jìn)行健康評價。

雅魯藏布江流域平均海拔高、水溫低、紫外輻射強(qiáng), 是西藏境內(nèi)最為主要的自然水生態(tài)系統(tǒng),能調(diào)節(jié)西藏的氣候特征, 具有眾多低海拔河流所不具有的特殊地貌及生態(tài)條件, 其生態(tài)環(huán)境十分敏感脆弱, 極易受到全球氣候變化和人類活動干擾等事件的影響, 保證其水環(huán)境擁有較高的健康狀態(tài), 是守護(hù)西藏綠水青山的重要基礎(chǔ), 研究雅魯藏布江流域水生態(tài)系統(tǒng)中有機(jī)生命體的組成和無機(jī)生活環(huán)境的狀態(tài)等, 將對整個流域的可持續(xù)發(fā)展有積極的推動作用。目前, 對雅魯藏布江流域進(jìn)行健康評價的研究報道較為少見, 故本文將以雅魯藏布江流域的周叢藻類為研究對象, 利用周叢藻類來構(gòu)建適合于雅魯藏布江流域的生物完整性指數(shù)(Periphytic index of biological integrity, P-IBI)綜合評價指標(biāo)體系[20], 以期對研究區(qū)域的水生態(tài)系統(tǒng)健康狀況進(jìn)行科學(xué)評價, 這將為西藏土著魚類食性和攝食機(jī)制等的研究提供餌料數(shù)據(jù)基礎(chǔ), 也能科學(xué)評價該區(qū)域河流水生態(tài)系統(tǒng)的健康狀況, 是保障西藏水生態(tài)安全,實現(xiàn)資源可持續(xù)發(fā)展與環(huán)境保護(hù)等作用的基礎(chǔ)。

1 材料與方法

1.1 研究區(qū)域和采樣時間

雅魯藏布江流域位于岡底斯-念青唐古拉山脈和喜馬拉雅山脈之間[21], 起源于中國西藏自治區(qū)日喀則市仲巴縣境內(nèi)的杰馬央宗曲, 依據(jù)流域地形特征與氣候類型對河道進(jìn)行分區(qū): 杰馬央宗曲至拉孜為上游, 是水流平緩的高原寬谷地帶, 河流兩岸以新月形沙丘居多, 為高原寒溫帶半干旱氣候,年平均降水量<300 mm, 流域面積2.67×104km2, 長約268 km,總落差1190 m, 占干流總落差的21.9%; 拉孜至派鎮(zhèn)為中游, 河谷寬窄交錯, 河流兩岸以河漫灘與階地居多, 多為高原溫帶半干旱氣候,年平均降水量300—600 mm, 流域面積16.5×104km2, 長約1293 km, 總落差1520 m, 占干流總落差的28.1%; 派鎮(zhèn)至巴昔卡為下游, 河道曲折迂回, 水流湍急, 河谷峽深, 高溫多雨,年平均降水量>4000 mm, 流域面積5.03×104km2, 長約496 km, 總落差2725 m, 占干流總落差的50.0%[22,23]。帕隆藏布江、尼洋河、拉薩河和年楚河是雅魯藏布江流域下游至上游沿岸的一級支流(圖1)。雅魯藏布江流域?qū)俑咴撅L(fēng)溫帶半干旱氣候, 有明顯的干濕季, 氣候垂直差異明顯, 其流域常年低溫,年平均最低氣溫為–3℃,年平均最高氣溫為16.2℃, 晝夜氣溫溫差可達(dá)14.7—16.5℃, 每年5—9月為集中降雨期, 降水量可達(dá)全年的80%[24],是西藏境內(nèi)最為主要的自然水域。

圖1 研究區(qū)域采樣點分布圖Fig.1 Distribution map of sampling points in the study area

本研究在雅魯藏布江流域干支流上設(shè)置了26個采樣點, 并先后于2013年10月、2014年4月和2014年7月, 開展了3次采樣工作。每個采樣點均用GPS定位, 采樣前需記錄每個點位的海拔、坐標(biāo)、河道底質(zhì)類型及構(gòu)成比例、河岸植被類型、水電工程等環(huán)境數(shù)據(jù)(表1), 隨后進(jìn)行周叢藻類樣品的采集。每次采樣的同一采樣點均需進(jìn)行3次平行采樣, 樣點分布如下:

干流采樣點設(shè)置: 上游2個采樣點, 中游6個采樣點(中上游、中下游各設(shè)置3個采樣點), 下游3個采樣點;支流采樣點設(shè)置:年楚河2個采樣點, 拉薩河4個采樣點、尼洋河4個采樣點、帕隆藏布江5個采樣點(表1)。

1.2 周叢藻類采集及物種鑒定

根據(jù)河道底質(zhì)類型選擇不同的采樣方法, 若底質(zhì)為黏土, 在采樣點河水中提前半個月放入20 cm×20 cm的玻璃板, 采樣時刷取玻璃板上的周叢藻類,若底質(zhì)以卵石為主, 各樣點隨機(jī)選取3—5塊河流里附有藻類的形狀規(guī)整、易于測量表面積的石塊, 用直徑為3.5 cm的塑料蓋蓋在石塊表面, 并用拇指固定石塊和蓋子, 先用硬毛刷和蒸餾水將蓋子以外的周叢藻類清除, 再用新的刷子將選取表面積的周叢藻類刷下。用蒸餾水多次沖洗玻璃板或石塊, 待沖凈, 將周叢藻類刷液全部轉(zhuǎn)入500 mL采樣瓶中, 立即用甲醛溶液固定, 帶回實驗室避光靜置48h后, 濃縮至50 mL樣品瓶中, 以待后續(xù)分析[25]。

周叢藻類種類鑒定與計數(shù)。周叢藻類的鑒定主要參照《中國淡水藻類》[26]和《中國西藏硅藻》[27]。取處理好的樣品, 于面積為20 mm×20 mm,容量為0.1 mL的計數(shù)板中, 在10×40倍光學(xué)顯微鏡下進(jìn)行觀察和計數(shù), 每片計數(shù)要大于300個, 每個樣品進(jìn)行2次平行計數(shù), 2次計數(shù)結(jié)果之差應(yīng)小于15%,否則重新計數(shù), 有效統(tǒng)計數(shù)值取平均值, 所得均值即為該片周叢藻類的數(shù)量[20]。

1.3 周叢藻類生物完整性指數(shù)(P-IBI)評價體系的構(gòu)建

P-IBI體系的建立: 采樣點數(shù)據(jù)資料收集-參考點和污染點的確定-備選指標(biāo)的篩選與確定-備選指標(biāo)賦分-P-IBI體系健康評價。

(1)采樣點數(shù)據(jù)資料收集: 對周叢藻類的每一個采集地所屬的河道特征和周邊環(huán)境數(shù)據(jù)進(jìn)行記錄和拍照留影;

(2)參考點和污染點的確定: 人類對西藏河流水環(huán)境的開發(fā)和利用, 相比于內(nèi)地的其他河流程度較低, 依據(jù)參考文獻(xiàn)[28, 29]的方法, 并結(jié)合西藏實際,采樣區(qū)域參考點需滿足如下三個條件: ①采樣點無人為干擾因素; ②采樣點無明顯人類活動干擾的跡象, 其上游無點污染源, 且上、下游5 km之內(nèi)無村莊, 兩側(cè)100 m寬5 km長的范圍內(nèi)無農(nóng)田; ③采樣點無明顯工程建設(shè), 或工程建設(shè)未對河岸兩側(cè)植被、附屬沼澤濕地、河流水質(zhì)等自然生境產(chǎn)生明顯影響。上述數(shù)據(jù)主要通過野外實地考察及實驗室數(shù)據(jù)分析獲得。在本研究的26個樣點中, 結(jié)合表1可知, 有21個樣點符合無干擾和干擾極小樣點被用作參照點, 剩余的5個采樣點為污染點(表2)。

表1 研究區(qū)域采樣點分布及信息Tab.1 Sampling point distribution and information

表2 采樣點中參考點和污染點的劃分Tab.2 The division of reference points and pollution points in sampling points

(3)備選指標(biāo)的篩選與確定: 根據(jù)已有的文獻(xiàn)[20], 選取12個應(yīng)用廣泛且入選率較高的指標(biāo)作為備選指標(biāo)(表3)。并對這12個備選指標(biāo)依次進(jìn)行:指數(shù)值分布范圍分析、判別能力分析、Pearson相關(guān)性分析, 以此篩選出合適的評價參數(shù)指標(biāo)[29]。①指數(shù)值分布范圍分析: 分析各指數(shù)值的分布范圍時, 需刪除分布范圍太窄、可預(yù)測環(huán)境變化值的范圍太小、分布散亂、標(biāo)準(zhǔn)差大、不穩(wěn)定的指標(biāo), 經(jīng)上述分析, 備選指標(biāo)全部適合無需剔除指標(biāo)。②判別能力分析: 將監(jiān)測期間的26個采樣點, 依據(jù)參考點和污染點進(jìn)行區(qū)分, 隨后以參考點和污染點為分組原則, 對12項備選指標(biāo)繪制箱線圖(Box plots), 并對12個箱線圖中參考點與污染點箱體25%—75%分位數(shù)范圍, 即箱體IQ(Interquartile ranges)的重疊部分進(jìn)行比較分析, 當(dāng)參考點與污染點的箱體沒有重疊部分, 或有部分重疊, 但參考點與污染點的中位數(shù)值都不在對方箱體之內(nèi), 則可進(jìn)入下一步計算分析(圖2), 反之, 則要剔除該項備選指標(biāo)。③Pearson相關(guān)性分析: 對通過判別能力分析的指標(biāo)進(jìn)行Pearson相關(guān)性分析, 依據(jù)相關(guān)性系數(shù)|r|>0.75, 進(jìn)行指標(biāo)選擇, 即: 當(dāng)幾個指標(biāo)的相關(guān)性>0.75時, 只需選出一個指標(biāo)代替其余指標(biāo)。

圖2 符合判別能力篩選條件的2類箱線圖[29]Fig.2 Types of box plots that meet the screening conditions of discriminant ability

表3 P-IBI 備選指標(biāo)描述及對干擾的反映Tab.3 P-IBI alternative indicator description and its response to interference

(4)備選指標(biāo)賦分:①三分制法: 對隨干擾減小的指標(biāo), 高于參考點的25%分位數(shù)值(對隨干擾增大的指標(biāo), 低于參考點的 75%分位數(shù)值), 計6分; 低于25%分位數(shù)值(高于75%分位數(shù)值)的分布范圍進(jìn)行二等分, 分別計3分和0分[30]。②四分制法: 對隨干擾減小的指標(biāo), 以所有點的95%分位數(shù)值為最佳期望值(對隨干擾增大的指標(biāo), 以所有點的5%分位數(shù)值為最佳期望值), 低于(高于)此值的分布范圍進(jìn)行三等分, 依次計6、4、2、0分[20]。③比值法: 對隨干擾增強(qiáng), 指數(shù)值減小的指標(biāo), 95%分位數(shù)值為最佳期望值, 各指標(biāo)的分值等于指標(biāo)值除以所有點的95%分位數(shù)值; 對隨干擾增大的指標(biāo), 以所有點的5%分位數(shù)值為最佳期望值進(jìn)行計算, 方法為: (最大值–指數(shù)值)/(最大值–5%分位數(shù)值), 該法規(guī)定, 經(jīng)計算后分值的分布范圍為0—1, 若大于1, 則記為1[29]。

(5)P-IBI體系健康評價: 以參考點P-IBI值分布的25%分位數(shù)作為健康評價的標(biāo)準(zhǔn), 若采樣點的PIBI值大于25%分位數(shù)值, 則表示該樣點受到的干擾很小, 水環(huán)境是健康的; 對小于25%分位數(shù)值的分布范圍進(jìn)行四等分, 分別確定為健康、亞健康、一般、較差、極差五個健康等級[29]。

1.4 數(shù)據(jù)分析

采用物種豐富度(Species Richness, SR), 總豐度(Total Abundance, TA), Shannon-Wiener多樣性指數(shù)(Shannon-Wiener Index, SH,H′)和Pielou均勻度指數(shù)(Pielou index, PI,J)對周叢藻類的群落穩(wěn)定性和分配均勻度等進(jìn)行分析。本研究根據(jù)各采樣點周叢藻類出現(xiàn)的頻率和相對豐度, 采用Mcnaughton優(yōu)勢度指數(shù)(Y)來判定雅魯藏布江流域周叢藻類優(yōu)勢種組成。具體如下: (1)物種豐富度(Species Richness, SR): 指某一采樣點或者季節(jié)出現(xiàn)的物種個數(shù)[31]。(2)總豐度(Total Abundance, TA): 指某一采樣點或者季節(jié)單位面積出現(xiàn)的周叢藻類細(xì)胞個體數(shù)(cells/cm2)[31]。(3)Shannon-Wiener指數(shù)H′[32]

式中,H′值0—1為重污, 1—3為中污(1—2為α-中污型, 2—3為β-中污型), >3為輕污或無污;pi為周叢藻類i占全部周叢藻類的比例(N%)。

(4)Pielou均勻度指數(shù)J[33]

式中,J值0—0.3為多污型, 0.3—0.4為α-中污型,0.4—0.5為β-中污型, >0.5輕污或清潔;S為周叢藻類種類總數(shù)。

(5)出現(xiàn)頻率: 指某一周叢藻類出現(xiàn)的次數(shù)占所有調(diào)查樣點數(shù)的百分比[31]。

(6)本研究采用Mcnaughton優(yōu)勢度指數(shù)(Y)來判定雅魯藏布江流域周叢藻類優(yōu)勢種組成, 將調(diào)查時間內(nèi)周叢藻類的優(yōu)勢種根據(jù)每個種類的優(yōu)勢度值(Y)來確定, 即將Y>0.02的種類判定為調(diào)查水域的優(yōu)勢種[34]。公式如下:

Y=(ni/N)×fi

式中,ni為第i種周叢藻類的個體數(shù),N為所有周叢藻類的總個體數(shù),ni/N為第i種周叢藻類的個體數(shù)占所有藻類總個體數(shù)的比值,fi為第i種周叢藻類出現(xiàn)的頻率。

本研究中的數(shù)據(jù)用Excel進(jìn)行統(tǒng)計, 同一采樣點所用周叢藻類數(shù)據(jù)均為3次采樣平均值。用R 3.6.1計算物種豐富度(Species richness, SR)、總豐度(Total abundance, TA)、Shannon-Wiener多樣性指數(shù)(H′)、Pielou均勻度指數(shù)(J)和Mcnaughton優(yōu)勢度指數(shù)(Y), 并用Spss22.0對數(shù)據(jù)進(jìn)行One-way ANOVA單因素方差分析, 差異顯著水平為P<0.05, 箱線圖的繪制和Pearson相關(guān)性的分析也使用Spss22.0。

2 結(jié)果

2.1 周叢藻類群落結(jié)構(gòu)特征分析

周叢藻類的種類組成和優(yōu)勢種分析采樣期間監(jiān)測結(jié)果顯示, 雅魯藏布江干流及其四大支流共鑒定出周叢藻類70屬(種), 分別隸屬于7門10綱23目34科。如表4所示, 包括硅藻門(Bacillariophyta)、綠藻門(Chlorophta)、藍(lán)藻門(Cyanophyta)、黃藻門(Xanthophyta)、裸藻門(Euglenophyta)、隱藻門(Cryptophyta)和金藻門(Chrysophyta), 其中硅藻門2綱4目9科25屬, 綠藻門2綱9目13科23屬, 藍(lán)藻門2綱5目7科14屬, 黃藻門1綱2目2科2屬, 裸藻門1綱1目1科4屬, 隱藻門和金藻門各1綱1目1科1屬。整個監(jiān)測期所有采樣點的周叢藻類優(yōu)勢門類所占比例(圖3), 依次為硅藻門(35.71%)、綠藻門(32.86%)、藍(lán)藻門(20.00%)、裸藻門(5.71%)、黃藻門(2.86%)、隱藻門(1.43%)和金藻門(1.43%), 表明周叢藻類的群落結(jié)構(gòu)整體上均屬于硅藻-綠藻-藍(lán)藻型水體。

圖3 雅魯藏布江流域周叢藻類的種類組成Fig.3 Species composition of periphyton in the Yarlung Zangbo River

以Mcnaughton優(yōu)勢度指數(shù)Y>0.02確定為優(yōu)勢種[34], 如表4所示, 整個調(diào)查期間, 調(diào)查水域周叢藻類的優(yōu)勢種有8屬, 且均歸屬于硅藻門, 包括脆桿藻科的3屬, 分別為針桿藻屬(0.092)、脆桿藻屬(0.020)和等片藻屬(0.132), 舟形藻科的舟形藻屬(0.072), 橋彎藻科的橋彎藻屬(0.040), 異極藻科的異極藻屬(0.031), 曲殼藻科的曲殼藻屬(0.165), 菱形藻科的菱形藻屬(0.193)。

表4 雅魯藏布江干流及其四大支流周叢藻類組成Tab.4 The periphyton algae composition of the sampling period of the main stream of the Yarlung Zangbo River and its four major tributaries

續(xù)表4

周叢藻類群落參數(shù)的空間分布對26個采樣點進(jìn)行監(jiān)測發(fā)現(xiàn), 周叢藻類物種豐富度(SR)的空間變化如下: 如圖4可知, 從干流到支流(A1—A26)周叢藻類的SR逐漸增大, 帕隆藏布江為雅魯藏布江干流下游的一級支流, 其(A12—A16)周叢藻類的SR明顯大于雅魯藏布江干流中下游(A2—A9)周叢藻類的SR, 尼洋河(A17—A20)、拉薩河(A21—A24)和年楚河(A25—A26)3條中上游支流周叢藻類的SR明顯大于雅魯藏布江干流上游(A10—A11)周叢藻類的SR。如圖5可知, 周叢藻類總豐度(TA)的空間變化如下: 帕隆藏布江(A12—A16)周叢藻類的TA明顯大于雅魯藏布江干流下游(A1—A3)周叢藻類的TA, 尼洋河(A17—A20)和拉薩河(A21—A24)周叢藻類的TA明顯大于雅魯藏布江干流中游(A4—A9)周叢藻類的TA,年楚河(A25—A26)周叢藻類的TA明顯大于雅魯藏布江干流上游(A10—A11)周叢藻類的TA。所以從流域整體來看, 雅魯藏布江流域干流周叢藻類的SR和TA均稍低于其相應(yīng)支流周叢藻類的SR和TA。

圖4 雅魯藏布江流域周叢藻類物種豐富度(SR)空間變化趨勢Fig.4 The spatial variation trend of periphyton algae species richness in the Yarlung Zangbo River Basin

圖5 雅魯藏布江流域周叢藻類總豐度(TA)空間變化趨勢Fig.5 The spatial variation trend of the total abundance of periphyton algae in the Yarlung Zangbo River Basin

監(jiān)測期間發(fā)現(xiàn), 26個采樣點周叢藻類的Shannon-Wiener多樣性指數(shù)(H′)的變幅為4.84—9.41, 均值為7.38(圖6)。H′值空間變化如下: 除雅魯藏布江流域干流中游(A4—A6、A8)的H′值有極為明顯的向下波動外, 其余水域的H′值波動范圍較平緩, 且雅魯藏布江干流下游墨脫背崩鄉(xiāng)(A1)的H′值最大,為9.41, 雅魯藏布江干流中游桑日縣古雪村(A6)的H′值最小, 為4.84。由圖7可知, Pielou 均勻度指數(shù)(J)的變幅為2.41—4.80, 均值為3.30。J值空間變化如下: 雅魯藏布江流域中游干流部分點位(A4—A6、A8)及支流尼洋河(A17—A20)的J值稍低, 但J值均>0.5, 此外, 其余水域均有較高的J值; 且雅魯藏布江干流下游的墨脫縣城上游(A3)的J值最大,為4.80, 尼洋河的巴河尼洋河交匯處(A19)的J值最小, 為2.41。故周叢藻類的Shannon-Wiener多樣性指數(shù)(H′)和Pielou 均勻度指數(shù)(J)共同表明: 整個雅魯藏布江流域的清潔狀態(tài)較好, 但干流中游水域的清潔狀態(tài)稍低于其余調(diào)查水域。

圖6 雅魯藏布江流域周叢藻類Shannon-Wiener多樣性指數(shù)(H′)空間變化趨勢Fig.6 The spatial variation trend of Shannon-Wiener diversity index of periphyton algae in the Yarlung Zangbo River Basin

圖7 雅魯藏布江流域周叢藻類Pielou 均勻度指數(shù)(J)空間變化趨勢Fig.7 Spatial variation trend of the Pielou evenness index of periphyton algae in the Yarlung Zangbo River Basin

參考點與污染點周叢藻類群落參數(shù)分析如表2所示, 26個采樣點被劃分為21個參考點和5個污染點, 將參考點和污染點周叢藻類的物種豐富度(SR)、總豐度(TA)、Shannon-Wiener多樣性指數(shù)(H′)和Pielou均勻度指數(shù)(J)進(jìn)行單因素方差性分析(表5)。結(jié)果顯示, 參考點周叢藻類的SR顯著大于污染點周叢藻類的SR(P=0.014<0.05), 參考點周叢藻類的TA小于污染點周叢藻類的TA, 但差異不顯著(P=0.326>0.05), 參考點周叢藻類的H′顯著大于污染點周叢藻類的H′(P=0.005<0.05), 參考點周叢藻類的J大于污染點周叢藻類的J, 且差異不顯著(P=0.730>0.05)。

表5 參考點與污染點周叢藻類的群落參數(shù)(平均值±標(biāo)準(zhǔn)差)Tab.5 Periphyton algae community parameters at reference point and pollution point (mean±SD)

2.2 雅魯藏布江流域健康評價

確定備選指標(biāo)判別能力分析(圖8)結(jié)果顯示, 非硅藻屬類單元數(shù)(M2)、Shannon-Wiener多樣性指數(shù)(M3)、周叢藻類物種種數(shù)(M6), 共3個參數(shù)指標(biāo)符合篩選規(guī)則, 可進(jìn)入下一步分析。

圖8 各入選 P-IBI 指標(biāo)在參考點和污染點的箱線圖Fig.8 Box plot of each selected P-IBI indicator at the reference point and pollution point

對上述步驟篩選出的3項指標(biāo)進(jìn)行Pearson相關(guān)性檢驗(SPSS 22.0; 表6), 結(jié)果顯示, 非硅藻屬類單元數(shù)(M2)、Shannon-Wiener多樣性指數(shù)(M3)和周叢藻類物種種數(shù)(M6)的相關(guān)性系數(shù)兩兩之間均為|r|<0.75, 符合備選要求。

表6 備選P-IBI指標(biāo)的Pearson相關(guān)性分析Tab.6 Pearson correlation analysis of alternative P-IBI indicators

依據(jù)以上篩選方法對12項周叢藻類的生物參數(shù)進(jìn)行選擇, 最終P-IBI指數(shù)由非硅藻屬類單元數(shù)(M2)、Shannon-Wiener多樣性指數(shù)(M3)和周叢藻類物種種數(shù)(M6)共3項指數(shù)構(gòu)成。

P-IBI 體系指標(biāo)分值計算及計分標(biāo)準(zhǔn)將篩選出的指標(biāo)分別應(yīng)用3分制法、4分制法、比值法計算各指標(biāo)分值(表7—9)。

根據(jù) 3 分制法的計分標(biāo)準(zhǔn)(表7), 將篩選所得指標(biāo)的分值進(jìn)行加和, 即獲得P-IBI值。以參考點的25%分位數(shù)為最佳期望值, 確定P-IBI指標(biāo)體系評價標(biāo)準(zhǔn)(表10), 進(jìn)而得出各采樣點周叢藻類生物完整性評分即水環(huán)境健康狀況評價, 結(jié)果如下: 在26個采樣點中, 18個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為健康, 5個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為亞健康, 2個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為一般, 1個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為較差; 三分制法對污染點的水環(huán)境健康狀態(tài)評價結(jié)果為: 較差-亞健康; 三分制法對參考點的水環(huán)境健康狀態(tài)評價結(jié)果為: 亞健康-健康; 三分制法對26個采樣點整體的水環(huán)境健康狀態(tài)評價結(jié)果為較差-健康(表11)。

表7 三個P-IBI 指標(biāo)在參照點中的分布及三分制法計分標(biāo)準(zhǔn)Tab.7 The distribution of 3 P-IBI indicators in the reference point and the scoring standard of the 3-point system

表10 三分制法P-IBI體系評價標(biāo)準(zhǔn)Tab.10 3-point method P-IBI system evaluation standard

表11 三分制法P-IBI體系評價結(jié)果Tab.11 Evaluation results of the P-IBI system based on the 3-point system

根據(jù)四分制法的計分標(biāo)準(zhǔn)(表8), 將篩選所得指標(biāo)的分值進(jìn)行加和, 即獲得P-IBI值。以參考點的25%分位數(shù)為最佳期望值, 確定P-IBI指標(biāo)體系評價標(biāo)準(zhǔn)(表12), 進(jìn)而得出各采樣點周叢藻類生物完整性評分即水環(huán)境健康狀況評價, 結(jié)果如下: 26個采樣點中, 16個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為健康, 3個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為亞健康, 4個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為一般, 3個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為較差;四分制法對污染點的水環(huán)境健康狀態(tài)評價結(jié)果為:較差-亞健康; 四分制法對參考點的水環(huán)境健康狀態(tài)評價結(jié)果為: 一般-健康; 四分制法對26個采樣點整體的水環(huán)境健康狀態(tài)評價結(jié)果為: 較差-健康(表13)。

表8 三個P-IBI 指標(biāo)在所有樣點中的分布及四分制法計分標(biāo)準(zhǔn)Tab.8 The distribution of 3 P-IBI indicators in all samples and the scoring standard of the 4-point system

表12 四分制法P-IBI體系評價標(biāo)準(zhǔn)Tab.12 4-point method P-IBI system evaluation standard

表13 四分制法P-IBI體系評價結(jié)果Tab.13 Evaluation results of the P-IBI system based on the 4-point system

根據(jù)比值法的計分標(biāo)準(zhǔn)(表9), 將篩選所得指標(biāo)的分值進(jìn)行加和, 即獲得P-IBI值。以參考點的25%分位數(shù)為最佳期望值, 確定P-IBI指標(biāo)體系評價標(biāo)準(zhǔn)(表14), 進(jìn)而得出各采樣點周叢藻類生物完整性評分即水環(huán)境健康狀況評價, 結(jié)果如下: 在26個采樣點中, 16個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為健康, 7個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為亞健康, 3個采樣點水環(huán)境的健康狀態(tài)評價結(jié)果為一般; 比值法對污染點的水環(huán)境健康狀態(tài)評價結(jié)果為一般-亞健康; 比值法對參考點的水環(huán)境健康狀態(tài)評價結(jié)果為亞健康-健康; 比值法對26個采樣點整體的水環(huán)境健康狀態(tài)評價結(jié)果為: 一般-健康(表15)。

表9 比值法計算3個 P-IBI 指標(biāo)分值的公式Tab.9 The formula for calculating the scores of 3 P-IBI indicators by ratio method

表14 比值法P-IBI體系評價標(biāo)準(zhǔn)Tab.14 Ratio method P-IBI system evaluation standard

表15 比值法P-IBI體系評價結(jié)果Tab.15 Evaluation result of P-IBI system by ratio method

采用三分制法、四分制法和比值法對雅魯藏布江流域的21個參考點和5個污染點分別進(jìn)行評分,并對3類評分方法的參考點和污染點得分, 進(jìn)行判別能力分析[29](圖9), 結(jié)果顯示, 以上3類評分方法的判別能力都很高, 均滿足箱線圖法的最佳規(guī)則,即采用任何一種方法均適用于本調(diào)查水域的健康評價。

圖9 三類評分法判別能力分析Fig.9 Analysis of discriminant ability of three types of scoring

3 討論

3.1 雅魯藏布江流域周叢藻類的組成特點

本研究水域共鑒定出周叢藻類7門10綱23目34科70屬(種), 周叢藻類的種類數(shù)稍多于學(xué)者們在尼洋河調(diào)查發(fā)現(xiàn)的6門30科49屬[31]、謝通門江段的6門67屬[35], 這是因為尼洋河和謝通門江段均為雅魯藏布江流域的部分水面, 其調(diào)查區(qū)域范圍稍小于本研究。但與內(nèi)地水域相比本調(diào)查區(qū)域周叢藻類的種數(shù)較少, 這可能與調(diào)查水域常年低溫有關(guān), 雅魯藏布江流域年平均氣溫為–3—16.2℃[24], 冬季水溫可低至2.8℃[35], 梁霞[36]指出: 大多數(shù)周叢藻類在水溫降到15℃以下, 藻墊和著生基質(zhì)間便開始發(fā)生脫落, 周叢藻類就開始大量凋亡, 只有少數(shù)周叢藻類能在低溫水域長期生長繁殖, 此外, 高原環(huán)境中超強(qiáng)的紫外輻射會損傷藻類的細(xì)胞結(jié)構(gòu), 也會使水環(huán)境中的藻類減少[37]。

調(diào)查水域周叢藻類的優(yōu)勢種包括針桿藻、脆桿藻、等片藻、舟形藻、橋彎藻、異極藻、曲殼藻和菱形藻, 優(yōu)勢種均隸屬于硅藻門周叢藻類, 這與硅藻喜低溫水環(huán)境生活有關(guān)[38]。雅魯藏布江流域周叢藻類中的硅藻門、綠藻門和藍(lán)藻門為前3個優(yōu)勢門類, 故該水域?qū)俟柙?綠藻-藍(lán)藻型結(jié)構(gòu), 是典型的高海拔河流水體類型[39], 這可能與河流主要依靠冰山融雪補給, 使得河流水質(zhì)清瘦[35], 水體中營養(yǎng)鹽含量較低, 供給周叢藻類的養(yǎng)分較少有關(guān)[40]。

3.2 周叢藻類對雅魯藏布江流域的健康評價

藻類物種多樣性的變化和分布特征, 既能直接反映水體中藻類的多樣性、均勻度和群落結(jié)構(gòu)穩(wěn)定性等的狀態(tài), 還能間接對水環(huán)境的健康狀態(tài)進(jìn)行評價[41—43], 所以進(jìn)行周叢藻類的生物多樣性指數(shù)特征性分析, 對水生態(tài)系統(tǒng)的發(fā)展意義重大。依據(jù)Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)數(shù)值分布規(guī)律, 數(shù)值越大說明周叢藻類的種類越豐富,周叢藻類分布越均勻, 水質(zhì)狀態(tài)越好[44], 反之, 數(shù)值越小說明周叢藻類的種類越稀少, 周叢藻類分布越不均勻, 水質(zhì)狀態(tài)越差[45]。在調(diào)查期間, Shannon-Wiener多樣性指數(shù)(H′)的變幅為4.84—9.41, 均值為7.38, 最小值>3(輕污或無污型水體)[32], Pielou 均勻度指數(shù)(J)的變幅為2.41—4.80, 均值為3.30, 最小值>0.5(輕污或清潔型水體)[33], 說明, 2013—2014年, 雅魯藏布江流域干支流整體上為輕污型-清潔型水體。

本研究基于周叢藻類的群落結(jié)構(gòu), 利用周叢藻類的群落參數(shù)(Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù))和P-IBI評價指標(biāo)體系對雅魯藏布江流域水生態(tài)系統(tǒng)進(jìn)行健康評價研究。結(jié)果顯示: 周叢藻類群落的Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和P-IBI評價指標(biāo)體系對雅魯藏布江流域的評價結(jié)果相一致。就周叢藻類群落的Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)而言, 污染點的Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)均較參考點小, 說明污染點的清潔狀態(tài)較參考點的清潔狀態(tài)差。在P-IBI評價指標(biāo)體系中, 污染點未有判定為健康狀態(tài)的點位, 這也表明污染點的清潔狀態(tài)較參考點的清潔狀態(tài)差。由表2可知, 劃定的污染點全部分布于雅魯藏布江流域干流中游, 故從流域整體來看, 雅魯藏布江流域干流上下游及其4大支流水生態(tài)系統(tǒng)的健康狀態(tài)優(yōu)于雅魯藏布江流域干流中游的健康狀態(tài)。

周叢藻類群落參數(shù)和P-IBI評價指標(biāo)體系的評價結(jié)果, 基本反映了雅魯藏布江流域的水環(huán)境和生境質(zhì)量特征。雅魯藏布江流域干流下游水域設(shè)置的采樣點周邊無明顯的水利水電工程建設(shè), 無采礦和挖沙等的破壞, 人為干擾較小, 水質(zhì)潔凈, 健康狀態(tài)好。雅魯藏布江流域干流中游水域設(shè)置的采樣點的周邊區(qū)域人類開發(fā)利用較大, 出現(xiàn)了諸多如渡口、采砂廠和水電站等工程建設(shè), 且河流沿岸兩側(cè)的生境出現(xiàn)了植被減少, 沙化嚴(yán)重的現(xiàn)象, 雖周叢藻類的Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)仍符合輕污-清潔型水體標(biāo)準(zhǔn), 但數(shù)值顯著減小, 說明周叢藻類的生境已受采樣點沿岸水利水電作業(yè)的影響。雅魯藏布江流域干流上游水域設(shè)置的采樣點的周邊區(qū)域也建設(shè)有水文站, 但采樣點的水質(zhì)清澈, 生境質(zhì)量較好, 故沒有將其劃分為污染點。帕隆藏布江周叢藻類Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)均在A13號和A16號采樣點較低, 與P-IBI評價指標(biāo)體系對A13號和A16號采樣點的評價結(jié)果出現(xiàn)亞健康狀態(tài)相一致, 這與A13號采樣點水流極其湍急, 不利于周叢藻類的附著生長有關(guān), A16號采樣點在比值法賦分的P-IBI評價指標(biāo)體系中呈亞健康狀態(tài), 在三分法和四分法賦分的P-IBI評價指標(biāo)體系中呈健康狀態(tài), 這可能與人工自河流引水入湖有關(guān), 且說明該采樣點的人為引水對周叢藻類的生境產(chǎn)生了影響, 但影響不大, 需隨時監(jiān)測該點的水環(huán)境健康狀態(tài)。尼洋河沿岸上設(shè)置的4個采樣點周邊存在大型施工機(jī)器和施工采石車輛, 但采樣點河岸兩側(cè)植被覆蓋率高, 自然生境保存較為完整, 拉薩河沿岸上設(shè)置的4個采樣點的周邊生境質(zhì)量較好, 河岸兩側(cè)有大量的植被, 沼澤濕地生態(tài)較好, 水質(zhì)清澈, 故尼洋河和拉薩河均不劃分為污染點, 且P-IBI評價指標(biāo)體系對尼洋河和拉薩河流域的評價結(jié)果均為健康狀態(tài), 與采樣點周叢藻類的Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)均符合輕污-清潔型水體標(biāo)準(zhǔn)相一致, 但尼洋河和拉薩河采樣點周叢藻類的Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)均有向下波動, 這可能與沿岸施工或該區(qū)域人口密集有關(guān), 建議加強(qiáng)對尼洋河和拉薩河流域水生態(tài)系統(tǒng)的監(jiān)測。年楚河沿岸設(shè)置了2個采樣點, 因A25號采樣點在城鎮(zhèn)周邊, 可能受到城鎮(zhèn)生活用水的影響, 為亞健康狀態(tài),A26號采樣點河畔植被豐富, P-IBI評價指標(biāo)體系表明其健康狀態(tài)較好, 這與該采樣點擁有較高的周叢藻類Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)相一致。

文獻(xiàn)[46—48]指出: 人為干擾會影響周叢藻類的群落結(jié)構(gòu), 從而指示水域不同等級的健康狀態(tài)。雅魯藏布江流域是西藏境內(nèi)最為主要的淡水來源,保證其水環(huán)境的健康將成為守護(hù)西藏青山綠水的基礎(chǔ), 研究雅魯藏布江流域水生態(tài)系統(tǒng)中有機(jī)生命體的組成和無機(jī)生活環(huán)境的狀態(tài)等, 將對整個流域的健康和穩(wěn)定有積極的作用。

4 結(jié)論

(1)2013—2014年, 在雅魯藏布江流域干支流共觀察到周叢藻類7門10綱23目34科70屬(種), 硅藻門(35.71%)、綠藻門(32.86%)和藍(lán)藻門(20.00%)為前3個優(yōu)勢門類, 優(yōu)勢種以硅藻門周叢藻類為主, 該水域?qū)俟柙?綠藻-藍(lán)藻型水體。(2)整個調(diào)查區(qū)域中,雅魯藏布江流域干流周叢藻類的物種豐富度(SR)和總豐度(TA)均稍低于支流; Shannon-Wiener多樣性指數(shù)的變幅為4.84—9.41, Pielou 均勻度指數(shù)的變幅為2.41—4.80, 則調(diào)查水域整體為輕污型-清潔型水體。(3)雅魯藏布江流域周叢藻類的生物多樣性特征與P-IBI評價體系共同說明: 雅魯藏布江流域干流上下游及其4大支流水生態(tài)系統(tǒng)的健康狀態(tài)優(yōu)于雅魯藏布江流域干流中游的健康狀態(tài)。