黃鑫 王菲 吳丹 孔凡青 王茜



摘 要:為探究白洋淀淀區大型底棲動物群落結構及其與環境因子的關系,于2020年10月在白洋淀淀區共設置10個采樣點位,對其進行大型底棲動物群落結構調查及測定pH值、電導率、總磷、總氮等9個水質理化指標。結果表明:白洋淀淀區共采集到大型底棲動物3門,18科,28屬,35種。其中,以節肢動物門為主,共11科,16屬,17種;軟體動物門次之,共5科,8屬,14種;環節動物門最少,2科,4屬,4種。基于PCA和RDA結果分析表明,COD、CODMn、TOC和pH值是影響白洋淀大型底棲動物群落結構的主要環境因子。
關鍵詞:白洋淀;大型底棲動物;群落結構;PCA分析;RDA分析
大型底棲動物主要包括扁形動物、環節動物、軟體動物、甲殼類和水生昆蟲等,是能量流動和物質循環的主要載體,在水環境生態系統中有重要作用[1-2]。 因其對環境變化較為敏感且耐受,可以綜合反映水體的污染程度,從而可以評價水質狀況[3],現已被廣泛應用于水生態系統的健康評價[4]。
主成分分析(Principal Component Analysis,PCA)最早由美國統計學家威爾遜[5]提出。是從多指標分析出發,運用統計分析原理與方法提取少數幾個彼此不相關的綜合性指標而保持其原指標所提供的大量信息的一種統計方法[6]。近年來,PCA分析在水源水質評價中已被廣泛應用[7-8]。冗余度分析(Redundancy Analysis,RDA)是一種直接梯度排序方法,能夠評價一組變量和另外一組變量數據之間的關系[9]。與其他排序方法相比較,RDA分析具有可以同時結合多個環境因子,結果更直觀、明確,且包含的信息量大等優點[10]。RDA分析可以很好地表現出環境因子對生物的影響程度,在生態研究領域具有較為重要的意義[11]。
白洋淀位于河北省雄安新區,由143個淀泊組成,總面積達366 km2,是中國華北地區最大的淺型淡水湖泊,也是中國北方典型的草型湖泊之一[12]。由于近年來生態補水、污染排放管控、入淀河流濕地建設等一系列的水生態環境保護修復措施的實施,白洋淀的水生態環境得到比較顯著的改善,水生態環境質量持續向好的趨勢基本形成。然而,水生態系統與水生生物的恢復是一個緩慢的過程。目前,白洋淀的水生生物多樣性在獲得恢復的趨勢下,水生生物狀況與有歷史記載的生物多樣性最好水平相比仍有差距。大型底棲動物較其他物種有較長的生活史,且遷移能力相對較差,因此可以作為環境指示物,其群落結構和空間布局也有助于反映水生態系統質量狀況。
1 材料與方法
1.1 采樣點布設
本研究根據白洋淀地理位置、生態環境、水文水資源狀況、湖泊面積與底質等因素,結合國家監測控制點位分布情況,于2020年10月在白洋淀淀區共設置10個采樣點位,對其進行水質理化指標及大型底棲動物采樣監測調查。采樣點分布見圖1。
1.2 樣品采集及分析檢測
使用規格為1/16的Peterson采泥器采集大型底棲動物定量樣品,過40目網篩篩洗干凈后放入白瓷盤中,肉眼挑撿出盤內所有大型底棲動物,于75%的酒精中保存。及時帶回實驗室在體視鏡及顯微鏡下進行計數和鑒定,絕大部分鑒定到種,少數種類鑒定到屬或更高的分類單元。根據個體數計算出大型底棲動物密度(ind./m2)。
大型底棲動物種類鑒定參考《中國經濟昆蟲志》[13]、《中國動物志》[14-17]、《淡水生物學》[18]、《中國經濟動物志》[19]、《中國北方搖蚊幼蟲》[20]以及《中國蜉蝣概述》[21]等。
水質理化指標的測定分為現場測定和實驗室檢測。現場使用多參數測定儀(YSI)測定pH、溶解氧(DO)和電導率(EC)等參數,測定方法參考《水和廢水監測分析方法》(第四版)。另外,采用有機玻璃采水器在水下0.5 m處取水樣,于車載冰箱中冷藏保存,并及時帶回實驗室進行前處理及分析,測定總磷(TP)、總氮(TN)、氨氮(NH3-N)、化學需氧量(COD)、高錳酸鹽指數(CODMn)和總有機碳(TOC)等參數,測定方法參考《水質 磷酸鹽和總磷的測定 連續流動-鉬酸銨分光光度法》(HJ 670-2013)、《水質 總氮的測定 連續流動-鹽酸萘乙二胺分光光度法》(HJ 667-2013)、《氮的測定 連續流動-水楊酸鈉分光光度法》(HHLY-LD 007-2013)、《水質 化學需氧量的測定 連續流動-重鉻酸鉀分光光度法》(HHLY-LD 002-2013)、《水質 高錳酸鹽指數的測定 連續流動-高錳酸分光光度法》(HHLY-LD 001-2013)和《水質 總有機碳的測定 燃燒氧化-非分散紅外吸收法》(HJ 501-2009)。
1.3 分析方法
1.3.1 優勢度及優勢種
采用大型底棲動物種類的優勢度值確定,Y>0.02的種類為優勢種,計算公式如下:
Y=(ni/N)×fi
式中,Y為優勢度;ni為種i的個體數;N為大型底棲動物總個體數;fi為種i出現的頻率。
1.3.2 多樣性指數計算及統計分析
采用Shannon-Wiener生物多樣性指數(H)和Evenness均勻度指數(J)2個指標分析大型底棲動物群落特征。Shannon-Wiener生物多樣性指數(H)包含的信息豐富,受物種的種類數影響較大,而Evenness均勻度指數(J)更多側重區分各種類密度的均勻性[22]。
Shannon-Wiener生物多樣性指數(H)的計算公式為:
式中,H為多樣性指數;n為大型底棲動物總個體數;S為大型底棲動物種類數;ni為種i的個體數。
Evenness均勻度指數(J)的計算公式:
J=H/log2S
式中,J為均勻度指數;H為香農指數;S為物種數。
1.3.3 主成分分析和冗余度分析
對白洋淀淀區10個采樣點的9項水質理化指標進行PCA分析,分析在SPSS 23中完成,以確定主要環境梯度。對大型底棲動物群落各類群的密度與主要環境因子進行RDA分析,RDA分析在Canoco4.5中完成。
2 結果
2.1 水質理化指標
白洋淀水體的理化指標檢測結果見表1。
2.2 大型底棲動物群落結構
2.2.1 種類組成及空間分布
在白洋淀大型底棲動物的監測調查中,共檢測到大型底棲動物3門,18科,28屬,35種。其中,以節肢動物門為主,共11科,16屬,17種,占總門類的48.57%;軟體動物門次之,共5科,8屬,14種,占總門類的40.00%;環節動物門最少,2科,4屬,4種,占總門類的11.43%,具體見圖2。其中,在監測調查的10個點位中,安新橋的物種數最少,僅為2種,且均為軟體動物;留通的物種數最多,為16種,以節肢動物占優勢。白洋淀大型底棲動物種類組成見圖2。
2.2.2 密度及空間分布
在白洋淀的監測調查中,10個監測點位的大型底棲動物的密度變化范圍較大,在39.60~339.90 ind./m2之間,平均密度為157.34 ind./m2,其中,安新橋密度最小,王家寨密度最大,白洋淀大型底棲動物密度變化見圖3。從物種的分布來看,中國圓田螺(Cipangopaludina chinensis)總密度最大,為202.62 ind./m2;中華顫蚓(Tubifex sinicus)、球圓田螺(Cipangopaludina amapulliformis)、尖口圓扁螺(Hippeutis cantori)、大臍圓扁螺(Hippeutis umbilicalis)、高翔蜉屬(Epeorus sp.)、沼蝦屬(Macrobrachium sp.)以及克氏原鰲蝦(Procambarus clarkii)總密度最低,均為3.33 ind./m2。
2.2.3 優勢種及多樣性指數
依據上述公式計算優勢度(Y)以確定優勢種,結果表明:白洋淀大型底棲動物的優勢種為水絲蚓屬(Limnodrilus sp.)、中國圓田螺、東北田螺(Viviparus chui)、槲豆螺(Bithynia misella)、紋沼螺(Parafossarulus striatulus)、七紋尾蟌(Cercion plagiosum)以及米蝦屬(Caridina sp.)。
基于上述公式計算Shannon-Wiener多樣性指數(H)和Evenness均勻度指數(J)。白洋淀大型底棲動物的Shannon-Wiener多樣性指數(H)介于0.41~3.31之間;Evenness均勻度指數(J)介于0.41~0.94之間。其中,郭里口、圈頭、留通和王家寨Shannon-Wiener多樣性指數(H)和Evenness均勻度指數(J)均較高。各點位多樣性指數結果見表2。
2.3 白洋淀水環境因子的PCA分析
對白洋淀9個水環境因子進行KMO檢驗和巴特利特球形度檢驗,得出KMO值為0.573,巴特利特球形度檢驗顯著性為0。一般認為KMO值大于0.7,則適合主成分分析,小于0.5以下則不合適,介于二者之間的為比較適用。巴特利特球形度檢驗的原假設是相關系數矩陣為單位陣(不相關),如果不能拒絕原假設,則不適合進行主成分分析[23]。因此,實驗數據可以進行PCA分析。
對白洋淀水環境因子進行PCA分析,基于特征根大于1的篩選原則,白洋淀水環境因子主要由3個成分組成,方差百分比依次為49.40%、30.12%和13.63%,累計貢獻率為93.15%。其中,F1和F2的累積貢獻率達到79.52%,可以表征白洋淀水環境狀況,故取前兩個主成分作為數據分析的有效成分。白洋淀PCA分析變量解釋見表3。白洋淀各環境因子在各主成分上的負載值及得分系數見表4。
根據各環境因子載荷量絕對值的大小以及各環境因子在成分矩陣中所占的得分系數,選擇對各個成分貢獻率大的指標作為該成分的因子[24]。對F1貢獻較大的環境因子包括COD、CODMn和TOC;對F2貢獻較大的環境因子為pH,由此確定影響白洋淀水體狀況的主要環境因子為COD、CODMn、TOC和pH。
2.4 白洋淀大型底棲動物群落結構與環境因子的關系
對白洋淀水環境因子與大型底棲動物的密度進行RDA分析。在RDA分析中,寡毛類和軟體動物密度與環境因子呈正相關關系;蜉蝣幼蟲與COD、pH和CODMn呈負相關關系;蛭綱、搖蚊幼蟲、蜻蜓目幼蟲和軟甲動物與環境因子呈負相關關系。白洋淀RDA分析結果見圖4。
3 討論
3.1 大型底棲動物群落結構特征
湖濱帶底質及人類活動對白洋淀物理生境的壓力較顯著,單一的污泥性質的底質影響著大型底棲動物群落結構,使其呈現單一化趨勢,優勢種也多以一般耐污種或耐污種為主。武志鑫等[25]在白洋淀淀區的兩次調查發現,大型底棲動物物種組成變化較小,主要為一般耐污的椎實螺科和耐污的田螺科、搖蚊幼蟲以及顫蚓科幼蟲為主;孫梟瓊等[26]對淀區30個點位進行調查發現,能在較高污染和缺氧環境下生存的搖蚊幼蟲在所有點位都有不同程度的分布,且以裸須搖蚊幼蟲為優勢物種。而本研究調查發現,白洋淀大型底棲動物主要以腹足綱動物和節肢動物為主,且主要以耐污種類田螺科為主,這與武志鑫等[25]、孫梟瓊等[26]和陳澤豪等[27]調查結果相似。
在監測調查的點位中,腹足綱動物和節肢動物出現頻率較高,幾乎所有點位都有分布,這主要與白洋淀生境特征有關。有研究認為水生植物是影響大型底棲動物分布的重要原因[28-29],這些水生植物不僅可以為節肢動物和腹足綱動物提供棲息環境,還可以富集浮游生物、有機碎屑等餌料供其生長發育。本次調查中,在留通及郭里口發現了對水質變化較為敏感的蜉蝣科稚蟲。從監測點位的地理位置來看,留通及郭里口水面開闊,靜水積蓄,較少受到人為干擾,是輕度干擾區,因此水質狀況相對較好,加上2018、2019、2020年度白洋淀持續大量補水,白洋淀水質狀況有所好轉[3,30],水生態系統結構變得復雜,也使得耐污種占比降低,敏感物種出現。
白洋淀大型底棲動物密度存在顯著的空間差異性。本次調查中,安新橋和圈頭物種多樣性低,密度也較低。安新橋和圈頭的透明度不足1 m,從而影響光的垂直分布,導致水體光照減弱,水生植物光合作用減弱,水體初級生產力隨之降低,大型底棲動物食物來源減少,進而導致密度降低。另外,通常情況下,大型底棲動物群落密度隨水深增加而呈遞減趨勢[31]。由于水深增加,水底光照強度減弱,水體初級生產力也隨之降低,進而導致食物來源較少[32]。但也有研究表明,大型底棲動物密度隨水深增加而呈遞增趨勢。研究發現,隨水深的增加,被淹沒的水生植物分解得越徹底,越有利于大型底棲動物的攝食[33]。依據王豆豆[34]在2019年秋季對白洋淀淀區的水深調查結果,白洋淀水深變化整體不超過1.5 m。由于白洋淀淀區水深變化范圍較小,缺乏相應的數據,大型底棲動物密度與水深的關系還有待于進一步研究。
3.2 多樣性指數
大型底棲動物活動能力較差,且對外界環境變化較為敏感,因此,自然因素和人類活動等均可影響大型底棲動物群落結構及物種多樣性。物種多樣性指數可以反映一個群落或生境中生物數量大小以及群落的分配狀況,可以直接或間接了解水生態系統的穩定程度及生境差異等[35-36]。留通和王家寨受人為干擾較小,且有機污染程度較低,多樣性指數較高,大型底棲動物種類數及多樣性也均較高;端村、采蒲臺、泥李莊、郭里口和光淀張莊受旅游業、運輸船只等人為影響較大,并在淀區可發現農業耕種用地,存在使用化肥等情況,可能造成大型底棲動物種類和密度較低;而安新橋和泥李莊不僅受人為影響較大,也受自然因素影響,其位于府河濕地入淀口附近,在暴雨季節仍存在府河水流入白洋淀的現象,進而導致其大型底棲動物種類較少,多樣性較低。另外,本研究中選擇的兩個多樣性指數的側重點不同,Shannon-Wiener生物多樣性指數(H)受物種的種類數影響較大,而Evenness均勻度指數(J)更側重區分各種類密度的均勻性。因此,在分析大型底棲動物群落結構時,可選用多個指標從物種的豐富度、均勻度、多樣性等多個方面進行評價,以此確保評價結果的全面性。
3.3 水環境因子對大型底棲動物的影響
水環境因子對大型底棲動物的影響相對復雜。不同類型的大型底棲動物的影響因素不一定相同,相同類型的大型底棲動物在不同區域也可能受到不同因素的影響。如郝韻等[37]認為,在湟水河流域高濃度氮、磷鹽對敏感物種有一定的限制作用;張遠等[38]認為,電導率是影響寡毛類分布的主要影響因素;池仕運等[39]認為,水溫、鹽度、總溶解固體和電導率是影響寡毛類分布的主要影響因素。本研究結果表明:影響白洋淀水體狀況的主要環境因子為COD、CODMn、TOC和pH。其中,COD、CODMn和TOC是評價水體有機污染程度的重要指標,數值越高,表明水體受有機污染越嚴重,進而影響大型底棲動物的分布、多樣性和密度等。pH值可通過影響大型底棲動物的繁殖能力,從而導致底棲無脊椎動物多樣性的顯著降低[40]。有研究表明,當pH值低于4或高于9時,大型底棲動物密度顯著降低[41]。
本研究中對污染耐受程度較高的寡毛類幼蟲和軟體動物密度隨有機污染物濃度增加而呈現升高趨勢,而其他類群則相反,這與本研究的前期分析結果相一致。在有機污染程度相對較高的泥李莊,其軟體動物密度顯著高于其他點位,同樣受有機污染的端村,其寡毛類密度也高于其他點位。環境因素影響大型底棲動物群落結構的變化并不是由單一因素決定,而是由多種因素(包括物理因素、生物因素、人為因素等)協同作用。因此,在探究影響大型底棲動物群落結構的關鍵影響因子時,除了需要考慮微觀影響因素(如水生植物、化學因子、生物因子等)外,還需要從宏觀的角度(如海拔、水文水資源狀況、土地利用類型等)綜合研究環境因素對大型底棲動物的影響。
參考文獻:
[1] 馮建社.白洋淀浮游植物與水質評價[J].江蘇環境科技,1999(2):27-29.
[2] 金磊,李林鈺,周楊,等.白洋淀三大典型水域浮游植物群落及水質評價[J].河北大學學報(自然科學版),2017,37(3):329-336.
[3] 易雨君,林楚翹,唐彩紅. 1960s以來白洋淀水文、環境、生態演變趨勢[J].湖泊科學,2020,32(5):1333-1347+1226.
[4] 朱浩,曹坤,陳曉龍,等.生態懸床技術對白洋淀水環境修復效果的研究[J].漁業現代化,2020,47(6):42-48.
[5] PEARSON K. On lines and planes of closest fit to systems of points in space[J].PhilosMag,1901,6(2):559-572.
[6] 王思遠,張增祥,趙曉麗,等.遙感與GIS技術支持下的湖北省生態環境綜合分析[J].地球科學進展,2002,17(3):426-431.
[7] 李萍.主成分分析在洱海水質評價中的應用[J].環境科學導刊,2023,42(3):87-90.
[8] 趙恰,趙明杰,馬學軍,等.秦皇島海港區主要河道水質分析及綜合評價[J].山東化工,2023,52(9):251-254.
[9] BRAAK T F J C.Canonical correspondence analysis:a new eigenvector technique for multivariate direct gradient analysis[J].Ecology,1986,67(5):1167-1179.
[10] 劉存歧,孔祥玲,張治榮,等.基于RDA的白洋淀浮游植物群落動態特征分析[J].河北大學學報(自然科學版),2016,36(3):278-285.
[11] 佟麗麗.淺談CANOCO4.5中RDA分析在水生態評價的應用[J].科技風,2018(2):47.
[12] 夏慧琪,龔家國,王浩,等.白洋淀生境空間格局演變分析[J].水利水電技術(中英文),2023,54(8):134-143.
[13] 中國科學院.中國經濟昆蟲志[M].北京:科學出版社,1997.
[14] 楊潼.中國動物志:環節動物門,蛭綱[M].北京:科學出版社,1996.
[15] 李新正.中國動物志:無脊椎動物,甲殼動物亞門,十足目,長臂蝦總科[M].北京:科學出版社,2007.
[16] 梁象秋.中國動物志:無脊椎動物,甲殼動物亞門,十足目,匙指蝦科[M].北京:科學出版社,2004.
[17] 任先秋.中國動物志:無脊椎動物,甲殼動物亞門,端足目,鉤蝦亞目[M].北京:科學出版社,2006.
[18] 大連水產學院.淡水生物學:上冊[M].北京:農業出版社,1982.
[19] 劉月英.中國經濟動物志:淡水軟體動物[M].北京:科學出版社,1979.
[20] 王俊才,王新華.中國北方搖蚊幼蟲[M].北京:中國言實出版社,2011.
[21] 周長發,蘇翠榮,歸鴻.中國蜉蝣概述[M].北京:科學出版社,2015.
[22] 蔡琨,陸維青,牛志春,等.洮滆水系湖泊春季浮游植物群落結構和水質生物學評價[J].環境監測管理與技術,2018,30(3):37-41.
[23] 郭芬. 遼河流域水生態與水環境因子時空變化特征研究[D].北京:中國環境科學研究院,2009.
[24] 熊春暉,張瑞雷,季高華,等.江蘇滆湖大型底棲動物群落結構及其與環境因子的關系[J].應用生態學報,2016,27(3):927-936.
[25] 武志鑫,王璽,劉瑀璇,等. 白洋淀流域大型底棲無脊椎動物群落結構分析[J].河北大學學報(自然科學版),2022,42(4):424-432.
[26] 孫梟瓊,李宏昌,廖晨延,等.白洋淀底棲無脊椎動物調查及水環境分析[J].河北漁業,2019(8):40-44+50.
[27] 陳澤豪,楊文,王穎,等.白洋淀大型底棲動物群落結構及其與環境因子關系[J].生態學雜志,2021,40(7):2175-2185.
[28] 張思強.貴州鎮寧縣洞穴大型底棲動物群落結構與環境因子的關系[D].貴陽:貴州師范大學,2022.
[29] CYR H,DOWNING J A.The abundance of phytophilous invertebrates on different species of submerged macrophytes[J].Freshwater Biology,1988,20:365-374.
[30] 崔志清.實施白洋淀生態補水 確保雄安新區水安全[J].中國水利,2020(15):62-63+39.
[31] 熊偉唯.瑪瑙河大型底棲動物時空分布特征及水環境評價研究[D].宜昌:三峽大學,2022.
[32] HU Y Y,XUE R Z,JIAN H D,et al.Spatial and temporal relation rule acquisition of eutrophication in Da'ning River based on rough set theory[J].Ecological Indicators,2016,66:180-189.
[33] CORTI R,DATRY T.Terrestrial and aquatic invertebrates in the riverbed of an intermittent river:parallels and contrasts in community organisation[J].Freshwater Biology,2016,61(8):1308-1320.
[34] 王豆豆.白洋淀沉水植物分布特征、凈水效果及重建策略[D].北京:北京林業大學,2019.
[35] ZHANG N N,Zang S Y. Characteristics of phytoplankton distribution for assessment of water quality in the Zhalong Wetland,China[J].International Journal of Environmental Science and Technology,2015,12(11):3657-3664.
[36] WANG H W.,XING Q C,ZHAO C L,et al. Investigation of phytoplankton and evaluation on water quality of Taizi River in Liaoyang area[J].Advanced Materials Research,2014,955-959: 403-406.
[37] 郝韻,尚光霞,丁森,等.湟水河流域秋季大型底棲動物群落結構及關鍵影響因子識別[J].水生態學雜志,2023,44(2):64-72.
[38] 張遠,徐成斌,馬溪平,等.遼河流域河流底棲動物完整性評價指標與標準[J].環境科學學報,2007(6):919-927.
[39] 池仕運,王瑞,汪鄂洲,等.三峽庫區生態牧場底棲動物群落特征及影響因素分析[J].環境保護科學,2022,48(3):20-28.
[40] 李永剛,胡慶杰,曲疆奇,等.北京密云水庫底棲動物群落結構及其時空變化[J].水生態學雜志,2018,39(5):31-38.
[41] BEREZIN A N. Influence of ambient pH on freshwater invertebrates under experimental conditions[J].Russian Journal of Ecology,2001,32(5):343-351.
The community structure of macrobenthos and its response to environmental factors in Baiyangdian Lake
HUANG Xin1,2,WANG Fei3,WU Dan3,KONG Fanqing2,WANG Qian1
(1.Key Laboratory of Aquatic-Ecology and Aquaculture of Tianjin,College of Fisheries,Tianjin Agricultural University,Tianjin 300384, China;2.Center of Eco-environmental Monitoring and Scientific Research,Administration of Ecology and Environment of Haihe River Basin and Beihai Sea Area,Ministry of Ecology and Environment of Peoples Republic of China,Tianjin 300170,China;3.Ecological Environment Monitoring Center of Hebei Province,Shijiazhuang 050037,China )
Abstract:To explore the community structure of macrobenthos in the Baiyangdian Lake and its relationship with environmental factors,10 sampling points were set up in Baiyangdian Lake in October 2020 to investigate the community structure of macrobenthos and determine 9 water quality physicochemical indicators such as pH value,conductivity,total phosphorus,and total nitrogen. The results showed that a total of 3 phyla,18 families,28 genera,and 35 species of macrobenthos were collected from the Baiyangdian Lake. Among them,arthropoda was the main phylum,consisting of 11 families,16 genera,and 17 species;The phylum mollusca took the second place,consisting of 5 families,8 genera and 14 species;The phylum annelida was the least,consisting of 2 families,4 genera and 4 species. Based on the analysis of PCA and RDA results,it was shown that COD,CODMn,TOC and pH value were the main environmental factors affecting the community structure of macrobenthos in Baiyangdian Lake.
Key words:Baiyangdian Lake;macrobenthos;community structure;PCA analysis;RDA analysis
(收稿日期:2023-09-13)
基金項目:國家自然科學基金(編號:31672264)。
作者簡介:黃鑫(1997-),女,碩士在讀,研究方向:水生態環境保護。E-mail: huangxin456280@163.com。
通訊作者:王茜(1971-),女,博士,教授,研究方向:水生動物及昆蟲分子系統學。E-mail: Wqgt1999@163.com。