杜紅春,朱挺兵,龔進(jìn)玲,胡飛飛,陳 康,楊德國

(中國水產(chǎn)科學(xué)研究院長江水產(chǎn)研究所，農(nóng)業(yè)農(nóng)村部淡水生物多樣性保護(hù)重點(diǎn)實(shí)驗(yàn)室，武漢430223)

瀾滄江發(fā)源于青海省唐古拉山東北部，流經(jīng)西藏、云南兩省區(qū)，出中國國境后被稱為湄公河，經(jīng)老撾、緬甸、泰國、柬埔寨，于越南胡志明市注入南海，是東南亞最大的國際河流，干流全長4 880 km，總落差5 500 m，流域面積81.0萬km。瀾滄江是典型的南向流向河流。瀾滄江西藏段主要位于西藏昌都市區(qū)域內(nèi)，長317 km，年平均流量為480～ 671 m/s，流域面積為23 167 km，屬于高原溫帶氣候。瀾滄江西藏段上游海拔較高，水生生物種類及數(shù)量較少，生物多樣性相對較低。

浮游動(dòng)物是水生態(tài)系統(tǒng)的重要組成部分，不僅是生態(tài)系統(tǒng)中的初級消費(fèi)者，還是其它高營養(yǎng)級水生動(dòng)物的食物，受上行效應(yīng)和下行效應(yīng)兩種機(jī)制的影響，且由于其對環(huán)境變化敏感度較高，因此在水質(zhì)監(jiān)測中常被用作指示物種。早期關(guān)于瀾滄江水系浮游動(dòng)物的研究主要集中在青海省河段、下游云南河段和水電開發(fā)密集區(qū)域，這些研究揭示了瀾滄江浮游動(dòng)物種類較少的總體特點(diǎn)，種類上主要以輪蟲為主，枝角類和橈足類的數(shù)量較少。但目前關(guān)于瀾滄江西藏段浮游動(dòng)物的研究還未見報(bào)道。本研究于2017-2020年對瀾滄江西藏段的浮游動(dòng)物進(jìn)行了調(diào)查，分析了浮游動(dòng)物的時(shí)空變化趨勢及其與環(huán)境因子的關(guān)系，以期為瀾滄江西藏段的生態(tài)環(huán)境保護(hù)和治理提供科學(xué)依據(jù)。

1 調(diào)查方法

1.1 采樣時(shí)間與斷面

分別于2017年4月(春季)、9月(秋季)、2018年4月(春季)、10月(秋季)、2019年4月(春季)、9月(秋季)和2020年11月(秋季)對瀾滄江西藏段浮游動(dòng)物進(jìn)行了調(diào)查，其中2020年春季因新冠肺炎疫情影響，未開展調(diào)查。根據(jù)瀾滄江西藏段河流特征，參考《水域生態(tài)系統(tǒng)觀測規(guī)范》共設(shè)置8個(gè)調(diào)查斷面，其中干流(曲孜卡、如美、卡若、扎曲)和支流(麥曲、金河、昂曲、熱曲)各4個(gè)。2019年4月金河斷面因上游電站調(diào)度導(dǎo)致河道干涸，未采集到樣本。采樣斷面分布如圖1所示。

圖1 瀾滄江西藏段浮游動(dòng)物采集斷面示意圖Fig.1 Sampling station of zooplankton in the Tiber reach of the Lancang River

1.2 樣品采集

采用25號浮游生物網(wǎng)在各采樣斷面近岸區(qū)水面下呈“∞”形滑動(dòng)采集浮游動(dòng)物定性樣品。原生動(dòng)物和輪蟲樣品根據(jù)水深等量采集不同水層的水樣混合，取1 L混合樣，加入魯哥試劑固定48 h后，濃縮至30 mL，用于定量分析。枝角類和橈足類定量樣品用5 L采水器取10 L水樣，用25號浮游生物網(wǎng)過濾后，收集于100 mL瓶中，加入4%甲醛溶液固定。浮游動(dòng)物的種類鑒定以及計(jì)數(shù)在顯微鏡下進(jìn)行。參照《淡水浮游生物研究方法》進(jìn)行密度和生物量的計(jì)算。

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

1.3.1 生物多樣性和優(yōu)勢度

采用Shannon-Wiener多樣性指數(shù)(′)、Pielou均勻度指數(shù)(′)和Margalef豐富度指數(shù)()來描述浮游動(dòng)物多樣性，采用Pinkas相對重要性指數(shù)()描述群落優(yōu)勢種的優(yōu)勢度，用Jacard群落相似性系數(shù)()來分析不同斷面浮游動(dòng)物群落的相似程度。

1.3.2 排序和分析

為探討瀾滄江西藏段浮游動(dòng)物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系，首先利用Canoco 4.5軟件包對浮游動(dòng)物豐度進(jìn)行去趨勢對應(yīng)分析(Detrended correspondence analysis,DCA)，根據(jù)DCA結(jié)果，最大梯度小于3，選用冗余分析(Redundancy analysis,RDA)。將所有監(jiān)測的水環(huán)境參數(shù)和浮游植物豐度數(shù)據(jù)引入分析中。浮游動(dòng)物豐度和除pH外的水環(huán)境因子數(shù)據(jù)以及浮游植物豐度進(jìn)行l(wèi)g(+1)轉(zhuǎn)化，再進(jìn)行蒙特卡洛檢驗(yàn)，將單個(gè)因子貢獻(xiàn)率值小于0.05作為顯著影響浮游動(dòng)物群落結(jié)構(gòu)分布的環(huán)境因子。

1.3.3 聚類分析

用SPSS 24.0軟件對浮游動(dòng)物群落進(jìn)行聚類分析，聚類分析采用分層聚類，分析前將浮游動(dòng)物豐度進(jìn)行l(wèi)g(+1)轉(zhuǎn)化，以降低極值帶來的影響。

2 結(jié)果

2.1 瀾滄江西藏段浮游動(dòng)物種類組成

共鑒定出浮游動(dòng)物4大類26屬28種(除無節(jié)幼體外)，其中原生動(dòng)物11屬11種，占總種類數(shù)的39.3%；輪蟲10屬10種，占總種類數(shù)的35.7%；枝角類3屬5種，占總種類數(shù)的17.9%；橈足類2屬2種，占總種類數(shù)的7.1%(圖2)。

圖2 瀾滄江西藏段浮游動(dòng)物種類組成Fig.2 Composition of zooplankton species in the Tibet reach of the Lancang River

從圖3可知，瀾滄江西藏段浮游動(dòng)物的種類數(shù)季節(jié)變化明顯。整體而言，春季的浮游動(dòng)物種類數(shù)大于秋季，但差異不顯著(=3586，=0.064)。春季浮游動(dòng)物種類數(shù)19種，占總種類數(shù)的65.5%；秋季13種，占總種類數(shù)的44.8%。2017年僅調(diào)查到枝角類和橈足類，其中春季有1種枝角類和2種橈足類，秋季枝角類和橈足類各1種。2018年春季僅調(diào)查到原生動(dòng)物(5種)和輪蟲(7種)，秋季沒有調(diào)查到浮游動(dòng)物。2019年春季調(diào)查到浮游動(dòng)物4門7種，其中原生動(dòng)物4種，輪蟲、枝角類和橈足類各1種；秋季僅原生動(dòng)物(4種)和橈足類(1種)。2020年秋季有浮游動(dòng)物4門8種，包括原生動(dòng)物3種、輪蟲和枝角類各2種、橈足類1種。

圖3 瀾滄江西藏段浮游動(dòng)物種類數(shù)的時(shí)空分布Fig.3 Temporal and spatial distribution of zooplankton species in the Tibet reach of the Lancang River

在8個(gè)采樣斷面中，扎曲采集到的浮游動(dòng)物種類數(shù)最多，有10種；其次為如美和熱曲，均8種；卡若和昂曲斷面的種類數(shù)最少，僅2種。干流共發(fā)現(xiàn)浮游動(dòng)物19種，其中原生動(dòng)物和輪蟲均7種，枝角類2種，橈足類3種；4條支流共發(fā)現(xiàn)浮游動(dòng)物15種，其中原生動(dòng)物7中，輪蟲4種，枝角類3種，橈足類1種。干流浮游動(dòng)物的種類數(shù)多于支流。從空間上看，干流和支流浮游動(dòng)物的種類數(shù)差異不顯著(=0409，=0.801)。

2.2 瀾滄江西藏段浮游動(dòng)物優(yōu)勢種

釆用IRI對瀾滄江西藏段定量樣品中的19種浮游動(dòng)物進(jìn)行優(yōu)勢種分析(表1)。其中優(yōu)勢種0種；重要種有2種，占10.5%，分別是鐘蟲(sp.)(144.9)和美麗猛水蚤(sp.)(212.7)；常見種6種，占31.6%；一般種9種，占47.4%；少見種僅2種，占10.5%(表1)。原生動(dòng)物大多都是常見種，輪蟲大多為一般種。

表1 瀾滄江西藏段浮游動(dòng)物優(yōu)勢種Tab.1 List of dominant species of zooplankton in the Tibet reach of Lancang River

2.3 瀾滄江西藏段浮游動(dòng)物密度和生物量

瀾滄江西藏段水體流動(dòng)太快，部分季節(jié)點(diǎn)位定量樣品中未采集到浮游動(dòng)物。浮游動(dòng)物密度的變化范圍為0～ 60 ind./L，均值為(5.95±2.20) ind./L。

續(xù)表1

原生動(dòng)物、輪蟲、枝角類、橈足類的密度均值分別為(4.49±1.94) ind./L、(1.31±1.17) ind./L、(0.04±0.03) ind./L和(0.11±0.07) ind./L。浮游動(dòng)物生物量的變化范圍為0～ 0.058 0 mg/L，均值為(0.004 2±0.002 0) mg/L，原生動(dòng)物、輪蟲、枝角類、橈足類的生物量均值分別為(0.000 2±0.000 1) mg/L、(0.001 6±0.001 4) mg/L、(0.000 7±0.000 6) mg/L和(0.001 7±0.001 1) mg/L。原生動(dòng)物在浮游動(dòng)物密度上占據(jù)優(yōu)勢，但是由于其體積較小，所以在生物量中的占比最小。

除去2020年秋季外，2017-2019年春秋季浮游動(dòng)物的密度均值分別為(4.92±2.77) ind./L和(3.92±2.01) ind./L，生物量均值分別為(0.003 4±0.002 5) mg/L和(0.003 4±0.002 0) mg/L。整體而言，春季浮游動(dòng)物的密度大于秋季，除2017年外，2018-2019年春季浮游動(dòng)物的生物量均大于秋季。不同季節(jié)浮游動(dòng)物的密度和生物量沒有顯著性差異(=0083，>005；=0000，>0.05)(圖4)。

圖4 瀾滄江西藏段浮游動(dòng)物密度和生物量的時(shí)空分布Fig.4 Temporal and spatial distribution of zooplankton density and biomass in the Tibet reach of the Lancang River

干流浮游動(dòng)物的豐度均值為(7.86±1.51) ind./L，生物量均值為(0.005 2±0.003 2) mg/L，所有斷面中，如美斷面的豐度均值最大(13.57±7.60) ind./L，曲孜卡斷面的生物量均值最大(0.010 5±0.008 9) mg/L，卡若斷面的豐度和生物量均值最小[(0.14±0.12) ind./L、(0.000 2±0.000 1) mg/L]。4條支流中，熱曲的豐度均值最大(9.14±7.35) ind./L，金河的豐度均值最小(0.33±0.18) ind./L；麥曲的生物量均值最高(0.006 7±0.005 4) mg/L，昂曲的生物量均值最低(0.000 7±0.000 6) mg/L。整體而言，干流的豐度和生物量均要大于支流但差異均不顯著(=0573，>005；=0239，>0.05)(圖4)。

2.4 瀾滄江西藏段浮游動(dòng)物生物多樣性指數(shù)

瀾滄江西藏段浮游動(dòng)物Shannon-Wiener多樣性指數(shù)的變化范圍為0～1.040，均值為(0.065±0.031)；Pielou均勻度指數(shù)的變化范圍為0～0.918，均值為(0.066±0.031)；Margalef豐富度指數(shù)的變化范圍為0～1.443，均值為(0.093±0.042)。浮游動(dòng)物的生物多樣性指數(shù)均較低，說明瀾滄江西藏段浮游動(dòng)物的種類較少，各物種數(shù)量較接近，群落結(jié)構(gòu)簡單。

各采樣斷面中，干流卡若斷面采集的定量樣品中僅1種浮游動(dòng)物，多樣性指數(shù)為0，除此之外，其余斷面的Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)的變化范圍分別為0.269～ 0.868、0.245～ 1.000和0.222～ 1.443。干流Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)的均值分別為(0.500±0.158)、(0.440±0.167)和(0.559±0.232)，相比各支流的多樣性指數(shù)較小。干流與支流的Shannon-Wiener多樣性指數(shù)(=0375，=0.825)、Pielou均勻度指數(shù)(=0456，=0.767)和Margalef豐富度指數(shù)(=0381，=0.821)差異均不顯著(圖5)。

各采樣時(shí)間中，2018年秋季采集的定量樣品中沒有發(fā)現(xiàn)浮游動(dòng)物，多樣性指數(shù)無法計(jì)算。春季的Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)均值分別為(1.260±0.356)、(0.814±0.095)和(1.648±0.621)，秋季的Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)均值分別為(0.589±0.197)、(0.703±0.207)和(0.344±0.134)。整體而言，春季的多樣性指數(shù)高于秋季，春季與秋季Shannon-Wiener多樣性指數(shù)(=4711，=0.035)、Pielou均勻度指數(shù)(=4869，=0.032)和Margalef豐富度指數(shù)(=5329，=0.026)差異均顯著(圖5)。

圖5 瀾滄江西藏段浮游動(dòng)物生物多樣性指數(shù)Fig.5 Biodiversity index of zooplankton in the Tibet reach of the Lancang River

2.5 群落結(jié)構(gòu)

對瀾滄江西藏段各斷面進(jìn)行群落相似性分析，表2結(jié)果顯示，相似性系數(shù)介于0～0.167之間，各斷面間浮游動(dòng)物群落結(jié)構(gòu)表現(xiàn)為極不相似。浮游動(dòng)物聚類分析結(jié)果如圖6所示，瀾滄江西藏段浮游動(dòng)物可分為3個(gè)群落。群落Ⅰ為曲孜卡，曲孜卡斷面僅有2種浮游動(dòng)物，且其它斷面均無該2種浮游動(dòng)物；群落Ⅱ?yàn)槿缑篮蜔崆?，浮游?dòng)物種類相似，以鐘蟲最具優(yōu)勢；群落Ⅲ為卡若、金河、昂曲、扎曲和麥曲。

圖6 瀾滄江西藏段浮游動(dòng)物群落結(jié)構(gòu)聚類分析結(jié)果Flg.6 Results from the cluster analysis on zooplankton community structure in the Tibet reach of the Lancang River

表2 瀾滄江西藏段各斷面浮游動(dòng)物群落組成相似性分析Tab.2 The similarity among the sampling sites in the Tibet reach of the Lancang River

2.6 游動(dòng)物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系分析

將瀾滄江西藏段浮游動(dòng)物密度進(jìn)行去趨勢對應(yīng)分析(DCA)，結(jié)果顯示最大排序梯度小于3，選用冗余分析。將瀾滄江西藏段7次采樣調(diào)查共計(jì)55個(gè)樣點(diǎn)的浮游動(dòng)物密度與環(huán)境因子(見表3)以及浮游植物密度(Phyto)進(jìn)行冗余分析，分析結(jié)果見圖7。在RDA排序結(jié)果中，第1、2排序軸的特征值分別為0.190和0.064，兩軸物種-環(huán)境相關(guān)性系數(shù)分別為0.793和0.734(表4)。經(jīng)蒙特卡洛檢驗(yàn)分析顯示，TSS(=0.014)、WT(=0.024)和Phyto(=0.032)是影響浮游動(dòng)物分布的最主要環(huán)境因子，其中TSS、WT和Phyto均與軸1呈負(fù)相關(guān)性，相關(guān)系數(shù)分別為-0.278 1、-0.300 5和-0.239 0；Phyto與軸2呈最大負(fù)相關(guān)(=-0.290 0)，WT與軸2呈正相關(guān)性(=0.192 3)。RDA結(jié)果顯示，大多數(shù)浮游動(dòng)物的分布與WT、TSS和浮游植物密度呈正相關(guān)。

表3 瀾滄江西藏段水環(huán)理化因子(平均值±標(biāo)準(zhǔn)誤)Tab.3 Enviromental factors in the Tibet reach of the Lancang River

表4 瀾滄江西藏段浮游動(dòng)物和環(huán)境因子間RDA分析的統(tǒng)計(jì)信息Tab.4 Statistics for the axes of RDA of zooplankton and environmental factors in the Tibet reach of the Lancang River

圖7 瀾滄江西藏段浮游動(dòng)物和環(huán)境因子RDA排序圖Fig.7 Redundancy analysis (RDA) ordination diagrams of the zooplankton and environmental factors in the Tibet reach of the Lancang River浮游動(dòng)物代碼見表1

3 討論

3.1 瀾滄江西藏段浮游動(dòng)物的群落結(jié)構(gòu)特征

本研究中共鑒定出浮游動(dòng)物29種，以原生動(dòng)物和輪蟲為主，枝角類和橈足類種類數(shù)相對較少，這與瀾滄江其它江段的調(diào)查結(jié)果基本一致：青海江段共發(fā)現(xiàn)浮游動(dòng)物24種，以輪蟲為主，原生動(dòng)物次之，枝角類和橈足類較少；下游江段共發(fā)現(xiàn)浮游動(dòng)物37種，以輪蟲為主；小灣庫區(qū)江段共發(fā)現(xiàn)浮游動(dòng)物22種，以輪蟲為主。整體而言，瀾滄江浮游動(dòng)物的群落結(jié)構(gòu)特征表現(xiàn)為以小型化的輪蟲和原生動(dòng)物為主，這與其它高原河流，如拉薩河、西藏尼洋河、雅魯藏布江的組成結(jié)構(gòu)相似。

瀾滄江西藏段各種浮游動(dòng)物的豐度和分布范圍較小，沒有較為明顯的優(yōu)勢種，重要種為美麗猛水蚤和鐘蟲。瀾滄江西藏段浮游動(dòng)物的豐度均值為(5.95±2.20) ind./L，原生動(dòng)物的豐度占比較高，輪蟲次之；生物量均值為(0.004 2±0.002 0) mg/L，橈足類的生物量最大，輪蟲稍次之。在浮游動(dòng)物類群中，原生動(dòng)物個(gè)體小，而枝角類和橈足類體型較大，導(dǎo)致浮游動(dòng)物密度和生物量在大小上存在不一致性。本次調(diào)查中，春季浮游動(dòng)物的密度和生物量基本均大于秋季，造成季節(jié)差異的原因可能是環(huán)境因素(降雨等)引起的。瀾滄江西藏段每年10月至次年4月，降雨稀少，為旱季，5-9月受印度洋暖濕氣流的影響，降雨增加，為雨季，雨量的增加會(huì)加快水體流動(dòng)，同時(shí)加快水體的上下攪動(dòng)，使水體中懸浮物增加，制約浮游動(dòng)物的生長。從流域上看，瀾滄江西藏段浮游動(dòng)物的密度和生物量與青海江段、下游江段的密度和生物量相差不大，數(shù)值均較小，可能是因?yàn)檎{(diào)查的河段都是活水狀態(tài)，流速較快，不適宜浮游動(dòng)物的生長。但是與2020年小灣庫區(qū)的密度差異較大，可能是因?yàn)樾硯靺^(qū)水流緩慢，幾近靜水狀態(tài)，更適宜浮游動(dòng)物的生長(表5)。

表5 瀾滄江各江段浮游動(dòng)物群落結(jié)構(gòu)Tab.5 Zooplankton community structure in different reaches of the Lancang River

3.2 瀾滄江西藏段浮游動(dòng)物密度與環(huán)境因子的關(guān)系

本研究表明，對瀾滄江西藏段浮游動(dòng)物分布影響顯著的環(huán)境因子是Phyto、TSS和WT。浮游植物作為浮游動(dòng)物的主要食物之一，能夠通過上行調(diào)控影響浮游動(dòng)物的密度，因此浮游植物密度與浮游動(dòng)物的生長繁殖密切相關(guān)。本文顯示浮游動(dòng)物密度與浮游植物密度呈顯著相關(guān)性，這與張軍燕等對拉薩河的的研究結(jié)果相同。有研究表明，懸浮物固體含量高的水體不利于浮游動(dòng)物的生長、繁殖和攝食，浮游動(dòng)物的密度會(huì)隨著懸浮物濃度的增加而降低，懸浮物濃度達(dá)到60.172 mg/L時(shí)，能使半數(shù)的浮游動(dòng)物死亡，同時(shí)水體中懸浮物含量較高，也會(huì)使水體中的透明度下降，進(jìn)而影響浮游動(dòng)物生長。瀾滄江西藏段水體的懸浮物均值為56 mg/L，濃度較高，可能是造成水體中浮游動(dòng)物密度較低的重要原因。同時(shí)有研究表明懸浮物固體的存在也會(huì)改變浮游動(dòng)物的群落組成，當(dāng)水體懸浮物含量增加時(shí)，枝角類數(shù)量不占優(yōu)勢，以輪蟲為主，瀾滄江西藏段浮游動(dòng)物群落種類組成就符合這一特征，枝角類和橈足類等大型浮游動(dòng)物的數(shù)量較少，小型動(dòng)物原生動(dòng)物和輪蟲的數(shù)量較多。有研究表明，水溫對浮游動(dòng)物的生長、繁殖和密度等的影響極為顯著，是影響浮游動(dòng)物群落變化的重要環(huán)境因子，水溫既可以通過作用于浮游動(dòng)物的生命周期直接影響浮游動(dòng)物的生長繁殖，也可以通過影響浮游植物的生長間接影響浮游動(dòng)物的生長繁殖，瀾滄江西藏段流域內(nèi)海拔為2 300～ 6 300 m，隨著海拔的升高，水體的溫度逐漸降低，調(diào)查期間，瀾滄江西藏段的平均水溫為(11.5±0.5) ℃，水溫較低，對浮游動(dòng)物的分布有較大影響，使浮游動(dòng)物的群落結(jié)構(gòu)呈現(xiàn)明顯的季節(jié)變化，如種類數(shù)和密度等呈現(xiàn)明顯的春高秋低趨勢。