趙奇,毛文
河南工業(yè)大學(xué),河南 鄭州 450001
浮游植物作為淡水生態(tài)系統(tǒng)中的初級生產(chǎn)者,是生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動的基礎(chǔ),其群落結(jié)構(gòu)會引起生態(tài)系統(tǒng)中食物網(wǎng)結(jié)構(gòu)的改變(張莉等,2019;文紅星等,2017;宋倫等,2016),從而影響水生生態(tài)系統(tǒng)的功能、結(jié)構(gòu)和信息傳遞(孫軍等,2016;邱陽凌等,2018;劉媛等,2017)。由于浮游植物對環(huán)境因子的響應(yīng)十分敏感,其種類組成、數(shù)量分布、群落結(jié)構(gòu)等是評價(jià)水環(huán)境質(zhì)量的重要指示作用,成為用來評價(jià)江河、湖泊、水庫等水質(zhì)變化的敏感指標(biāo),在水體的監(jiān)測和評價(jià)中起極為重要的應(yīng)用意義(鐘愛文等,2017;許秀梅等,2017)。浮游植物的群落結(jié)構(gòu)與生活的水質(zhì)狀況密切相關(guān),利用浮游植物來評價(jià)和監(jiān)測水質(zhì)的研究也在逐步開展(趙秀俠等,2017;劉璐等,2017)。近年來,對河流水質(zhì)指標(biāo)研究得比較多,但關(guān)于浮游植物群落結(jié)構(gòu)及水質(zhì)的季節(jié)變化較少。
雁鳴湖水域是河南省黃河沿岸濕地重要的組成部分,對于區(qū)域氣候調(diào)節(jié)、生物多樣性保護(hù)、涵養(yǎng)水源等有重要意義(金磊等,2017)。雁鳴湖也是重要的漁業(yè)基地和水利區(qū),發(fā)揮著調(diào)節(jié)氣候、保護(hù)生物多樣性、發(fā)展生態(tài)經(jīng)濟(jì)、維護(hù)生態(tài)平衡、水土保持等重要功能(劉明典等,2017)。在全球氣候變化和人類活動影響下雁鳴湖水環(huán)境發(fā)生了嚴(yán)重的改變,已面臨農(nóng)業(yè)面源污染、有機(jī)污染、生活垃圾污染等嚴(yán)峻的生態(tài)環(huán)境問題(宋慶洋等,2019;李晶等,2017;嚴(yán)廣寒等,2018)。有鑒于此,本文連續(xù)5年較系統(tǒng)地調(diào)查和分析了不同季節(jié)雁鳴湖浮游植物群落組成及水質(zhì)狀況特征,運(yùn)用灰關(guān)聯(lián)分析與生物多樣性指數(shù)相結(jié)合的方法,分析浮游植物的多樣性和季節(jié)變動,評價(jià)雁鳴湖水質(zhì)狀況,為水生態(tài)環(huán)境的保護(hù)和水生態(tài)基準(zhǔn)的建立、水生資源恢復(fù)提供科學(xué)依據(jù)。
雁鳴湖地處鄭州市區(qū),屬于中緯度地帶,半干旱半濕潤大陸性季風(fēng)氣候區(qū)。1957—2006年年平均氣溫14.4 ℃;年平均降水量614.2 mm,降水量年內(nèi)分布不均,7—9月降水量占全年降水量的60%—70%。研究區(qū)土壤以潮土、風(fēng)沙土為主,西北部耕作區(qū)主要為潮土;湖面北側(cè)主要為魚塘;趙口總干渠南側(cè)防護(hù)林區(qū)為沙崗地,主要為風(fēng)沙土。
2014—2018連續(xù)5年不同季節(jié)(春季2月、夏季5月、秋季8月、冬季11月)定期(月中)、定點(diǎn)(每個(gè)區(qū)15個(gè)采樣點(diǎn),圖1)采集樣品,現(xiàn)場測定水溫、透明度,pH和采用多參數(shù)水質(zhì)監(jiān)測儀(HI 9828)測定,定性、定量采集樣品測定浮游植物多樣性及水質(zhì)狀況。浮游植物(定性測量樣品)以國際標(biāo)準(zhǔn)的25號浮游生物網(wǎng)在水面表層呈“∞”字形緩慢來回拖取3—5 min撈取浮游植物樣品,濃縮生物網(wǎng)中的水到100 mL后1—2 mL魯哥試液固定,帶回實(shí)驗(yàn)室以備鏡檢,另取表層水樣1 L,搖勻后帶回實(shí)驗(yàn)室用于測定水質(zhì)狀況;浮游植物(定量測量樣品)使用有機(jī)玻璃采水器在表層采取水樣1 L,10—15 mL魯哥試劑固定,帶回實(shí)驗(yàn)室后經(jīng)24 h沉淀濃縮至30 mL,加入4%甲醛溶液保存以備鏡檢(張莉等,2019;文紅星等,2017)。
圖1 研究區(qū)圖Fig.1 The study area
1.3.1 浮游植物的鑒定
將定量樣品搖勻后,在顯微鏡(400倍)下進(jìn)行觀測,采用特定的浮游生物計(jì)數(shù)框(Palmer Counting Cell)計(jì)數(shù),選取20—40個(gè)視野,每個(gè)樣本重復(fù)計(jì)數(shù)3次,有效統(tǒng)計(jì)數(shù)值。計(jì)算結(jié)果為藻類細(xì)胞密度,即單位體積內(nèi)藻類細(xì)胞個(gè)體數(shù)表示,對于比較難判斷的藻類,則任選20個(gè)個(gè)體在高倍鏡下觀察,測出細(xì)胞數(shù)取均值,依據(jù)《中國淡水藻類志》鑒定藻類組成。
選用Margalef豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)、Simpson優(yōu)勢度指數(shù)和Pielou均勻度指數(shù)分析浮游植物多樣性:
Margalef豐富度指數(shù)(S):
Shannon-Wiener多樣性指數(shù)(H):
Simpson優(yōu)勢度指數(shù)(D):
Pielou均勻度指數(shù)(JP):
式中,N為浮游植物總數(shù);ni為第i種的浮游植物個(gè)體數(shù)量。
1.3.2 水質(zhì)樣品測定
水質(zhì)測定項(xiàng)目包括TN、TP、BOD5、CODCr;TN采用過硫酸鉀-紫外分光光度法;TP采用鉬銻抗分光光度法;BOD5采用稀釋接種法;CODCr采用重鉻酸鉀氧化法(劉璐等,2017)。
水質(zhì)監(jiān)測實(shí)測值及歸一化數(shù)值:為求得關(guān)聯(lián)度,最終對雁鳴湖的水質(zhì)進(jìn)行評價(jià),需先將水質(zhì)監(jiān)測實(shí)測值和與之比較的《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)相應(yīng)的指標(biāo)值按照方法進(jìn)行歸一化處理,雁鳴湖春、夏、秋季水質(zhì)監(jiān)測實(shí)測均值及歸一化數(shù)值結(jié)果見表1、2,關(guān)聯(lián)度及水質(zhì)評價(jià):根據(jù)公式對上述歸一化數(shù)據(jù)進(jìn)一步計(jì)算得雁鳴湖各采樣點(diǎn)水質(zhì)與《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》Ⅰ—Ⅴ類水質(zhì)等級關(guān)聯(lián)度見表3。
1.4.1 水質(zhì)標(biāo)準(zhǔn)的選擇
在灰關(guān)聯(lián)分析中采用《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》作為評價(jià)標(biāo)準(zhǔn)(劉璐等,2017)。
1.4.2 參考數(shù)列、比較數(shù)列的確定
設(shè)《標(biāo)準(zhǔn)》為比較數(shù)列,選擇相應(yīng)的雁鳴湖實(shí)測值為參考數(shù)列。分別用X0={X0(k)|k=1, 2, …,n}表示參考數(shù)列,Xi={Xi(k)|k=1, 2,…,n}(i=1, 2,…,m)表示比較數(shù)列。
1.4.3 歸一化處理
水質(zhì)標(biāo)準(zhǔn)的量級、單位都不盡相同,因此對參考和比較數(shù)列進(jìn)行歸一化處理。方法:將Ⅰ級水體水質(zhì)標(biāo)準(zhǔn)濃度設(shè)為1,m級水體水質(zhì)標(biāo)準(zhǔn)濃度設(shè)為0。經(jīng)過線性變換計(jì)算,1—m級之間的水質(zhì)標(biāo)準(zhǔn)在0—1之間。處理后的比較數(shù)列、參考數(shù)列分別用X0′和Xi′表示。
表1 雁鳴湖水質(zhì)理化指標(biāo)及其歸一值Table 1 Physical and chemical indexes and normalized values of Yanming Lake water quality
表2 《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)歸一化數(shù)值Table 2 Normalized value of Surface Water Environmental Quality Standard (GB 3838—2002)
表3 雁鳴湖各采樣點(diǎn)關(guān)聯(lián)度與水質(zhì)等級Table 3 Correlation degree and water quality grade of each sampling point in Yanming Lake
關(guān)聯(lián)系數(shù)的求解:用X0′和Xi′數(shù)列中的數(shù)據(jù),求解這兩數(shù)列間的灰關(guān)聯(lián)系數(shù)。公式如下(劉媛等,2017):
其中,(k)i是第k項(xiàng)指標(biāo)X0′和Xi′的絕對差;min(i)min(k)Δi(k)與max(i)max(k)Δi(k)分別是各項(xiàng)指標(biāo)絕對差中最小和最大值,Δi表示分辨系數(shù),取0.5。
關(guān)聯(lián)度的求解方法:
式中,γi代表關(guān)聯(lián)度;k代表第k項(xiàng)指標(biāo);n表示數(shù)列中k的總數(shù)。
1.4.5 關(guān)聯(lián)度的排序
關(guān)聯(lián)度的大小反映了實(shí)測數(shù)據(jù)與《標(biāo)準(zhǔn)》的相似程度,關(guān)聯(lián)度越大,此水域接近于該級的水質(zhì)標(biāo)準(zhǔn),判定此水域?yàn)樵摷壦|(zhì)。
由表1可知,雁鳴湖水質(zhì)季節(jié)性變化較為明顯,水溫變化范圍在5.11—8.13 ℃之間,歸一化指數(shù)在0.722—0.887之間;pH實(shí)測值變化范圍在7.15—7.62之間,歸一化指數(shù)在0.828—0.913之間;水溫、pH、BOD5和CODCr、TN和TP呈一致的變化規(guī)律,在夏季和秋季最高,春季和冬季較低,基本表現(xiàn)為夏季>秋季>春季>冬季,局部有所波動,透明度呈“V”字形變化規(guī)律,在春季最高,夏季最低,基本表現(xiàn)為春季>冬季>秋季>夏季。應(yīng)用灰關(guān)聯(lián)分析評價(jià)雁鳴湖水質(zhì)類別(表3),呈現(xiàn)不同季節(jié)水質(zhì)與《標(biāo)準(zhǔn)》中的Ⅰ類水質(zhì)的關(guān)聯(lián)度最大,分別為0.712、0.823、0.856、0.768,因此,判定不同季節(jié)水質(zhì)為Ⅰ類水質(zhì)。
在雁鳴湖采樣水體中,共檢出浮游植物7門231種(表4),其中浮游植物種類在夏季和秋季達(dá)到較高,春季和冬季較低,不同季節(jié)均以綠藻門(Chlorophyta)和硅藻門(Bacillariophyta)為主,分別占38.10%和24.24%;甲藻門(Euglenophyta)18種,所占比例為7.79%;金藻門(Euglenophyta)15種,所占比例為6.49%;藍(lán)藻門(Cyanophyta)29種,所占比例為12.55%;隱藻門(Euglenophyta)20種,所占比例為8.66%,裸藻門(Euglenophyta)5種,所占比例為2.16%。
根據(jù)以上初步的配合比(02號)中水膠比和砂率按《水工混凝土試驗(yàn)規(guī)程》(SL352-2006)規(guī)定進(jìn)行調(diào)整,水膠比增減0.05、砂率相應(yīng)增減1%,得出01和03號調(diào)整后的配合比進(jìn)行試配,如表6共3個(gè)配合比進(jìn)行試配,表7為3組配合比試配后得出各項(xiàng)混凝土性能指標(biāo)。
表4 雁鳴湖浮游植物群落結(jié)構(gòu)Table 4 Phytoplankton community structure in Yanming Lake
不同植物多樣性指數(shù)顯示(圖2),雁鳴湖浮游植物Shannon-Wiener多樣性指數(shù)(H)、Pielou均勻度指數(shù)(JP)、Margalef種類豐富度指數(shù)(S)在夏季和秋季較高,春季和冬季較低,隨季節(jié)呈先增加后降低趨勢,而浮游密度(D)在秋季達(dá)到最低。其中浮游植物密度變化范圍在0.67×106—1.98×106cells·L-1之間,夏季和冬季差異不顯著(P>0.05);多樣性指數(shù)變化范圍在0.95—1.79之間,春季和冬季差異不顯著(P>0.05),夏季和秋季差異不顯著(P>0.05);均勻度指數(shù)變化范圍在5.45—9.51之間,不同季節(jié)差異均顯著(P<0.05);豐富度指數(shù)變化范圍在0.55—0.89之間,春季和秋季差異不顯著(P>0.05)。
由表5可知,多樣性指數(shù)春、夏、秋、冬季均值分別為1.65、0.95、0.79、1.79,根據(jù)生物多樣性指數(shù)的評價(jià)標(biāo)準(zhǔn),雁鳴湖秋季和夏季屬于α-中污型,春季和冬季屬于β-中污型;均勻度指數(shù)春、夏、秋、冬季均值分別為0.66、0.89、0.73,0.55,根據(jù)生物多樣性指數(shù)的評價(jià)標(biāo)準(zhǔn),雁鳴湖春季、秋季和冬季屬于清潔-寡污型,夏季屬于清潔型;豐富度指數(shù)春、夏、秋、冬季均值分別為6.23、9.51、7.67、5.45,根據(jù)豐富度指數(shù)的評價(jià)標(biāo)準(zhǔn),雁鳴湖屬于清潔型。
圖2 雁鳴湖浮游植物多樣性的影響Fig.2 Influence of phytoplankton diversity in Yanming Lake
表5 Shannon-Wiener指數(shù),Margalef指數(shù)和Pielou指數(shù)的評價(jià)標(biāo)準(zhǔn)Table 5 Evaluation criterion about Shannon-Wiener, Margalef and Pielou index
由表6可知,雁鳴湖春季浮游植物多樣性與透明度呈顯著負(fù)相關(guān)(P<0.05),與CODCr濃度呈極顯著正相關(guān)(P<0.01),與BOD5濃度和TN濃度呈顯著的正相關(guān)關(guān)系(P<0.05);夏季浮游植物多樣性與透明度呈極顯著負(fù)相關(guān)(P<0.01),與TN濃度、TP濃度、BOD5濃度和CODCr濃度呈極顯著正相關(guān)(P<0.01),與水溫呈顯著的正相關(guān)關(guān)系(P<0.05);秋季浮游植物多樣性與透明度呈極顯著負(fù)相關(guān)(P<0.01),與TN濃度、TP濃度、BOD5濃度和CODCr濃度呈極顯著正相關(guān)(P<0.01),與水溫呈顯著的正相關(guān)關(guān)系(P<0.05);冬季浮游植物多樣性與透明度呈極顯著負(fù)相關(guān)(P<0.01),與TN濃度、BOD5濃度和CODCr濃度呈顯著正相關(guān)(P<0.05)。
由表7可知,前2個(gè)排序軸特征值分別為0.615和0.232,第一排序軸可反映不同生境下環(huán)境因子的梯度變化特征,浮游植物多樣性與環(huán)境因子2個(gè)排序軸的相關(guān)性均為1.000,前2個(gè)排序軸特征值占總特征值的95.12%,蒙特卡羅檢驗(yàn)分析環(huán)境因子對植物的影響達(dá)到顯著性(第一軸P=0.001,F(xiàn)=8.56;第二軸P=0.001,F(xiàn)=6.11),能夠很好地解釋環(huán)境因子對植物的影響。RDA排序圖(圖3)顯示雁鳴湖浮游植物多樣性與水質(zhì)狀況均呈顯著正相關(guān)關(guān)系(P<0.01),沿RDA第1排序軸,雁鳴湖浮游植物多樣性與水質(zhì)狀況的季節(jié)變化規(guī)律相一致,而雁鳴湖浮游植物多樣性影響最大的是水BOD5和CODCr濃度。
表6 雁鳴湖浮游植物多樣性與水質(zhì)指標(biāo)的關(guān)系Table 6 Relationship between phytoplankton diversity and water quality indexes in Yanming Lake
表7 RDA排序結(jié)果Table 7 Results from ordination
圖3 雁鳴湖浮游植物多樣性與環(huán)境因子的RDA排序圖Fig.3 RDA sequence diagram of phytoplankton diversity and environmental factors in Yanming Lake
本研究中,雁鳴湖共檢出浮游植物7門251種,對照區(qū)共檢出浮游植物7門166種,主要種類為綠藻門,可以初步推斷雁鳴湖藻類植物的群落組成主要為綠藻類型;浮游植物種類在夏季和秋季達(dá)到較高,春季和冬季較低??傮w而言,Shannon-Wiener多樣性指數(shù)(H)、Pielou均勻度指數(shù)(JP)、Margalef種類豐富度指數(shù)(S)隨季節(jié)呈先增加后降低趨勢,在春季和冬季溫度較低時(shí),雁鳴湖的浮游植物種類比較少,但仍然是綠藻門和硅藻門浮游植物占主導(dǎo),與前人的研究結(jié)果相一致(Eva,2001;Katsiapi et al.,2011;Webber et al.,2005)。夏季和秋季,水溫較高,在加上水質(zhì)中TN和TP濃度較高,導(dǎo)致浮游植物種類數(shù)豐富。本研究中浮游植物在夏季密度最大,秋季最低,參照評價(jià)水質(zhì)的標(biāo)準(zhǔn),浮游植物細(xì)胞密度≤5×105cells·L-1,水體為極貧營養(yǎng)水體,≤1.0×106cells·L-1為貧營養(yǎng),1.0×106—9.0×106cells·L-1為貧中營養(yǎng)。由此可見,雁鳴湖秋季水質(zhì)最差,為極貧營養(yǎng)水體,夏季和冬季達(dá)到貧中營養(yǎng)水平。浮游植物多樣性指數(shù)是判斷水庫營養(yǎng)狀況最常用的檢測指標(biāo),H>3為輕或無污染,3—2為β-中度污染(輕中污染),2—1為α-中度污染(重中污染),0—1為重度污染,也即H指數(shù)值越大,水質(zhì)越好,即藻類的種類多樣性指數(shù)越高,其群落結(jié)構(gòu)越復(fù)雜(Chen et al.,2009;Gharib et al.,2011;Ciudice et al.,2010)。根據(jù)生物多樣性閾值的分級評價(jià)標(biāo)準(zhǔn),夏季和秋季浮游植物種間個(gè)體數(shù)分布較為均勻,但種類豐富程度較低,水體受到一定污染,同時(shí),夏季和秋季浮游植物種間個(gè)體數(shù)分布欠均,種類豐富程度偏低,這與當(dāng)?shù)厣钗鬯呐湃胗嘘P(guān)。另外,夏季溫度較高,雨量較為豐沛,上游、支流為雁鳴湖輸入了較多的營養(yǎng)物質(zhì),適合于浮游植物的生長(Danaher et al.,2010;Dondajewska et al.,2019)。從多樣性的季節(jié)變化來看,春季和冬季生物多樣性升高、種類的豐富程度加大、群落中優(yōu)勢成分的集中程度降低,水體趨向健康(張莉等,2018;喻龍等,2017)。另一方面,夏季和秋季浮游動物對浮游植物較高的攝食壓力也影響了浮游植物的分布、數(shù)量秋季,雁鳴湖水溫回降,雖然浮游植物的多樣性仍然豐富,種類的豐富程度較高,但水體的健康程度急劇下降。
冗余分析(RDA)作為一種直接梯度分析方法,不僅維持了良好的景觀,而且環(huán)境各因子對浮游植物也表現(xiàn)出明顯的反饋?zhàn)饔?。冗余分析顯示,沿RDA第1排序軸,雁鳴湖浮游植物多樣性與水質(zhì)狀況的季節(jié)變化規(guī)律相一致,而雁鳴湖浮游植物多樣性影響最大的是BOD5和CODCr濃度,與Pearson相關(guān)分析的結(jié)果基本一致。水質(zhì)的研究結(jié)果顯示水溫、pH、TN、TP、BOD5和CODCr濃度呈一致的變化規(guī)律,在夏季和秋季最高,春季和冬季較低,基本表現(xiàn)為夏季>秋季>春季>冬季,局部有所波動,透明度呈“V”字形變化規(guī)律,在春季最高,夏季最低,基本表現(xiàn)為春季>冬季>秋季>夏季。從季節(jié)的變化規(guī)律來看,春季和冬季污染較小,隨著水溫的增加,各種水生生物和淡水藻類瘋狂生長,加上旅游開發(fā)后人類對雁鳴湖的污染,大部分水域出現(xiàn)水體富營養(yǎng)化,在夏季污染程度達(dá)到最大,導(dǎo)致水質(zhì)中TN、TP、BOD5和CODCr濃度也達(dá)到最大,春季和冬季淡水藻類和水生生物數(shù)目急劇降低,水質(zhì)狀況有所好轉(zhuǎn)(Akomeah et al.,2010)。淡水水體中N、P長期以來被認(rèn)為是與浮游植物生長密切相關(guān),其中P被廣泛認(rèn)為是淡水浮游植物演替的影響因素,但在雁鳴湖TP并非影響浮游植物群落分布的主要驅(qū)動力。
(1)雁鳴湖共檢出浮游植物7門231種,其中浮游植物種類在夏季和秋季達(dá)到較高,春季和冬季較低,雁鳴湖浮游植物Shannon-Wiener多樣性指數(shù)(H)、Pielou均勻度指數(shù)(JP)、Margalef種類豐富度指數(shù)(S)在夏季和秋季較高,春季和冬季較低,隨季節(jié)呈先增加后降低趨勢,而浮游密度(D)在秋季達(dá)到最低。
(2)雁鳴湖秋季和夏季屬于α-中污型,春季和冬季屬于β-中污型;從均勻度指數(shù)來看,雁鳴湖春季、秋季和冬季屬于清潔-寡污型,夏季屬于清潔型;根據(jù)豐富度指數(shù)的評價(jià)標(biāo)準(zhǔn),雁鳴湖屬于清潔型。應(yīng)用灰關(guān)聯(lián)分析評價(jià)雁鳴湖水質(zhì)類別,不同季節(jié)水質(zhì)為Ⅰ類水質(zhì)。
(3)相關(guān)性分析顯示雁鳴湖浮游植物群落多樣性與BOD5和CODCr濃度呈顯著或極顯著正相關(guān),與pH呈顯著負(fù)相關(guān),由此表明雁鳴湖更多受氮等營養(yǎng)元素干擾,并且污染物面積相對寬闊,需注意控制加強(qiáng)面源污染,合理規(guī)劃農(nóng)業(yè)土地利用,盡量減少化肥施用量等,有利于控制營養(yǎng)元素流入漢庫,對于浮游植被的恢復(fù)有重要意義。
(4)浮游植物多樣性對水溫、TN濃度和NH4+-N濃度反應(yīng)較為敏感,在某種程度上反映了其指示作用,能夠表征雁鳴湖浮游植物的敏感性,但不同季節(jié)植物多樣性與環(huán)境因子之間的內(nèi)在聯(lián)系還需要深入研究。