陸建芳

?

山仔水庫浮游植物群落特征及水體富營(yíng)養(yǎng)化研究

陸建芳

福州市環(huán)境監(jiān)測(cè)中心站

對(duì)山仔水庫2016年7～9月藍(lán)藻水華發(fā)生期間浮游藻類及部分理化指標(biāo)進(jìn)行了監(jiān)測(cè)，采用克里金空間插值得到葉綠素a空間分布，采用浮游植物細(xì)胞密度指示法和TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)法對(duì)庫區(qū)2016年7～9月營(yíng)養(yǎng)狀況進(jìn)行評(píng)價(jià)，探討了浮游植物群落、葉綠素a含量與水體富營(yíng)養(yǎng)化之間的關(guān)系。結(jié)果表明，2016年7～9月山仔水庫浮游植物種類以綠藻門為主，藍(lán)藻門的微囊藻和銅綠微囊藻在種群數(shù)量上占絕對(duì)優(yōu)勢(shì)，葉綠素a空間分布結(jié)果和TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)結(jié)果都顯示，2016年7月七里和皇帝洞、2016年9月皇帝洞水質(zhì)為輕度富營(yíng)養(yǎng)化，其它監(jiān)測(cè)時(shí)間水體為中營(yíng)養(yǎng)，監(jiān)測(cè)點(diǎn)位營(yíng)養(yǎng)化程度為：皇帝洞>七里>日溪>庫心>壩前。

山仔水庫 浮游植物群落 富營(yíng)養(yǎng)化 葉綠素a 空間插值

1 概述

湖泊水庫富營(yíng)養(yǎng)化問題是當(dāng)前世界范圍內(nèi)普遍關(guān)注的水環(huán)境問題之一，特別是作為飲用水源的湖泊水庫的富營(yíng)養(yǎng)化，直接威脅到人們的飲水安全[1]。富營(yíng)養(yǎng)化過程實(shí)際是水體中營(yíng)養(yǎng)物質(zhì)增加，水生生物特別是浮游藻類大量增殖, 并在水中建立優(yōu)勢(shì)的過程[2]。所以浮游藻類的種類組成和數(shù)量變化是水體富營(yíng)養(yǎng)化的重要指標(biāo)之一[3]，水體中葉綠素 a是浮游植物現(xiàn)存量的重要指標(biāo)，其含量的高低能夠反映水體的營(yíng)養(yǎng)狀況。因此，葉綠素a常作為湖泊富營(yíng)養(yǎng)化調(diào)查的一個(gè)主要參數(shù)，在水體富營(yíng)養(yǎng)狀況評(píng)價(jià)中起關(guān)鍵性作用[4]。目前，被廣泛采用的水體富營(yíng)養(yǎng)狀況評(píng)價(jià)方法，主要是從水質(zhì)、水生生物以及沉積物 3 個(gè)方面來進(jìn)行，所采用的評(píng)價(jià)方法可概括為六大類型，即特征描述法、參數(shù)法、生物指示評(píng)價(jià)法、磷收支模型法、營(yíng)養(yǎng)狀態(tài)指數(shù)法和數(shù)學(xué)分析法[5-7]，其中生物指示評(píng)價(jià)法簡(jiǎn)單、方便、可以直觀評(píng)價(jià)水體的富營(yíng)養(yǎng)化狀況[8]，而基于葉綠素a等多因子的TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)法是目前最常用的湖泊水庫富營(yíng)養(yǎng)化評(píng)價(jià)方法[9]。

本文研究對(duì)象山仔水庫位于敖江中游，水庫呈長(zhǎng)條不規(guī)則狀分布，水庫總庫容1.723億m3，調(diào)節(jié)庫容1.064億m3，為福州市第二飲用水源。近年來，由于人為活動(dòng)對(duì)庫區(qū)生態(tài)系統(tǒng)的破壞，使得該水庫水體富營(yíng)養(yǎng)化污染特征明顯，夏秋季時(shí)有藍(lán)藻水華發(fā)生，給城市供水質(zhì)量造成安全隱患[2]。本研究選擇在藍(lán)藻水華易發(fā)的夏季（7～9月）對(duì)庫區(qū)浮游藻類及部分理化指標(biāo)進(jìn)行監(jiān)測(cè)，其中，7月增加監(jiān)測(cè)點(diǎn)位，補(bǔ)充監(jiān)測(cè)葉綠素a含量，應(yīng)用地學(xué)統(tǒng)計(jì)中的空間插值方法，直觀呈現(xiàn)葉綠素a的空間分布。探討浮游植物群落、葉綠素a含量及水體富營(yíng)養(yǎng)化之間的關(guān)系，旨在為山仔庫區(qū)藻類及水體富營(yíng)養(yǎng)化研究提供基礎(chǔ)資料，以期為水華預(yù)警和防治提供參考。

2 研究方法

2.1 采樣點(diǎn)位和頻次

在山仔水庫布設(shè)5個(gè)固定采樣點(diǎn)位，分別為壩前（出口）、庫心、七里（進(jìn)口）、日溪（支流進(jìn)口）、皇帝洞。于2016年7月5日、8月1日和9月2日采樣，其中7月5日除固定5個(gè)點(diǎn)位外增加10個(gè)補(bǔ)充監(jiān)測(cè)點(diǎn)。

2.2 樣品采集和處理

樣品采集方法依據(jù)《湖泊生態(tài)調(diào)查觀測(cè)與分析》[10]進(jìn)行。

水質(zhì)監(jiān)測(cè)項(xiàng)目主要有pH、水溫（WT）、溶解氧（DO）、葉綠素a（Chl a）、透明度（SD）、總氮（TN）、氨氮（NH4-N）、硝酸鹽氮（NO3-N）、總磷（TP）、可溶性正磷酸鹽（PO43--P）、高錳酸鹽指數(shù)（CODMn），分析方法依據(jù)《水和廢水監(jiān)測(cè)分析方法（第四版）》[11]。

浮游植物定量采集0.5m水深的表層水，取1L水樣置于塑料瓶中，現(xiàn)場(chǎng)用魯哥氏液固定，樣品帶回實(shí)驗(yàn)室后靜置沉降48 h，沉淀濃縮后用浮游植物計(jì)數(shù)框計(jì)數(shù)。

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

3.1 綜合營(yíng)養(yǎng)狀態(tài)指數(shù)

采用國(guó)家環(huán)境保護(hù)部推薦方法，以chl的狀態(tài)指數(shù)TLI（chl a）為基準(zhǔn)，選取TN、TP、SD、CODMn等4個(gè)參數(shù)的狀態(tài)指數(shù)，同TLI（chla）一起進(jìn)行相關(guān)加權(quán)綜合。

綜合營(yíng)養(yǎng)狀態(tài)指數(shù)計(jì)算式：

TLI(chl a)＝10(2.5+ 1.086 lnchl a)；

TLI(TP)＝10(9.436+1.624 lnTP)；

TLI(TN)＝10(5.453+1.694 lnTN)；

TLI(SD)＝10(5.118-1.94 lnSD)；

TLI(CODMn)＝10(0.109+2.661 lnCODMn)。

結(jié)果參照湖泊營(yíng)養(yǎng)狀態(tài)分級(jí)指標(biāo)：TLI(∑)<30表示貧營(yíng)養(yǎng)，30≤TLI(∑)≤50表示中營(yíng)養(yǎng)，TLI(∑)>50表示富營(yíng)養(yǎng)，5070重度富營(yíng)養(yǎng)，在同一營(yíng)養(yǎng)狀態(tài)下，指數(shù)值越高，其營(yíng)養(yǎng)程度越重。

3.2葉綠素a含量空間插值

本文選取克里金插值法進(jìn)行山仔水庫葉綠素a的空間分布研究?？死锝穑↘ringing）插值法是空間插值最優(yōu)方法[12]，它是基于空間結(jié)構(gòu)分析后的理論變差函數(shù)模型，認(rèn)為推測(cè)值和已知值具有空間上的相關(guān)性，這種相關(guān)性隨空間屬性（距離和方向）的變化而變化。當(dāng)確定一個(gè)已知值的有效影響距離后，確定有效影響范圍，基于變差函數(shù)模型，根據(jù)有效范圍內(nèi)的已知值（實(shí)測(cè)葉綠素 a 濃度值）來推測(cè)待插值點(diǎn)的葉綠素 a 的濃度。這種方法尋求最小方差，是空間插值方法中最優(yōu)的選擇[13]。

本實(shí)驗(yàn)所有數(shù)據(jù)均采用Excel 2003錄入、圖表處理和結(jié)果統(tǒng)計(jì)分析。葉綠素a空間插值應(yīng)用ArcGIS地理統(tǒng)計(jì)分析 Geostatistical Analyst 模塊完成。

4 結(jié)果與分析

4.1 浮游植物種群特征

2016年7月～9月，庫區(qū)各采樣點(diǎn)共鑒定出浮游植物7個(gè)門41屬57種, 門類涵蓋藍(lán)藻門、硅藻門、裸藻門、隱藻門、綠藻門、甲藻門和金藻門。裸藻門和金藻門分別只在8月、9月出現(xiàn)。從種類上看（圖1），綠藻門的種類遠(yuǎn)高于其他門類，共鑒定出22種；其次為硅藻門和藍(lán)藻門，分別鑒定出9種和7種；甲藻門、裸藻門和金藻門種數(shù)較少。各月份種類組成略有不同，優(yōu)勢(shì)種以藍(lán)藻門的銅綠微囊藻、魚腥藻出現(xiàn)時(shí)間最長(zhǎng)。

圖1 山仔水庫浮游植物種類組成

從浮游植物種群數(shù)量組成上看，監(jiān)測(cè)各月，藍(lán)藻門均占絕對(duì)優(yōu)勢(shì)，其次是綠藻門和硅藻門。由圖2可知，7月，藍(lán)藻門數(shù)量占到了99.9%，其中以微囊藻和銅綠微囊藻的數(shù)量最多。8月和9月綠藻門比重增大，分別為4.36%和4.78%，優(yōu)勢(shì)種為鼓藻和美麗網(wǎng)球藻，但此時(shí)浮游植物總數(shù)僅為5.18×105個(gè)/L和5.01×105個(gè)/L?？傮w來看，夏季各月山仔水庫雖然綠藻門出現(xiàn)的種類數(shù)最多，但其在數(shù)量上并不占優(yōu)勢(shì)，所以山仔水庫夏季浮游植物組成是以藍(lán)藻為主的類型。

圖2 各門類浮游植物種群數(shù)量百分比

4.2 富營(yíng)養(yǎng)化等級(jí)評(píng)價(jià)

本文分別采用浮游植物細(xì)胞密度評(píng)價(jià)法和TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)法，評(píng)價(jià)了監(jiān)測(cè)期間5個(gè)固定點(diǎn)位的水體富營(yíng)養(yǎng)化狀況，結(jié)果見表1、表2?？梢钥闯?，兩種評(píng)價(jià)方法的結(jié)果存在差異，各點(diǎn)位應(yīng)用TLI(∑)評(píng)價(jià)法的評(píng)價(jià)等級(jí)均高于浮游植物細(xì)胞密度指示評(píng)價(jià)法，TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)法是基于葉綠素a的多因子評(píng)價(jià)指數(shù)，較單因子的浮游植物細(xì)胞密度評(píng)價(jià)法而言，能更為客觀、全面地反映水體水質(zhì)。兩種評(píng)價(jià)方法中，庫區(qū)各點(diǎn)位間綜合營(yíng)養(yǎng)狀態(tài)指數(shù)等級(jí)比較，總體趨勢(shì)均為：皇帝洞>七里>日溪>庫心>壩前。

表1 浮游植物細(xì)胞密度指示評(píng)價(jià)法

表2 綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)法

一般認(rèn)為，水庫富營(yíng)養(yǎng)化同庫區(qū)外源性污染物輸入和內(nèi)源性污染物釋放關(guān)系密切。山仔水庫每年夏季藍(lán)藻水華高發(fā)，主要由于季節(jié)變換時(shí)期，水溫回升，水中微生物分解速率加快，向水體釋放更多的營(yíng)養(yǎng)元素，促使藻類進(jìn)入生長(zhǎng)旺季。不同空間點(diǎn)位富營(yíng)養(yǎng)化程度存在差異，則可能與汛期降雨徑流影響，導(dǎo)致上游村莊的外源污染物入庫有關(guān)。七里、皇帝洞位于水域上游淺水區(qū)，氮磷等外源營(yíng)養(yǎng)物沖刷入庫后使水中藻類迅速繁殖，所以較壩前等出口位置更具富營(yíng)養(yǎng)化風(fēng)險(xiǎn)。為防控山仔水庫富營(yíng)養(yǎng)化污染，保證水源地安全，福州市開展整治畜禽養(yǎng)殖污染、控制工業(yè)企業(yè)污染、加強(qiáng)庫區(qū)周邊農(nóng)村環(huán)?；A(chǔ)設(shè)施等行動(dòng)，并于2009年起采取生態(tài)養(yǎng)魚防控藻類，每年投放鰱魚、鳙魚，在山仔水庫內(nèi)形成藻類消納氮、磷等物質(zhì)，魚吃藻的良性生態(tài)循環(huán)，優(yōu)化水庫水生系統(tǒng)?？陀^上看，山仔水庫近幾年藻類水華暴發(fā)頻次降低，富營(yíng)養(yǎng)化程度有所改善，但是根據(jù)歷年監(jiān)測(cè)數(shù)據(jù)顯示，庫區(qū)藻類數(shù)量通常從5月開始增加，7月達(dá)到頂峰，個(gè)別點(diǎn)位依然伴有藍(lán)藻水華發(fā)生，9月后隨著水溫降低，藻類繁殖力降低。2016年夏季監(jiān)測(cè)結(jié)果也符合這一規(guī)律。用兩種富營(yíng)養(yǎng)化評(píng)價(jià)法評(píng)價(jià)，結(jié)果均表明，7月水體富營(yíng)養(yǎng)化程度高于8、9月。TLI(∑)指數(shù)顯示，七里和皇帝洞點(diǎn)位在7月已達(dá)到輕度富營(yíng)養(yǎng)化水平，所以入夏時(shí)節(jié)對(duì)水庫藻類水華暴發(fā)仍需重點(diǎn)關(guān)注。

4.3 葉綠素a濃度的反演與分析

應(yīng)用克里金空間插值法對(duì)7月山仔水庫葉綠素a含量進(jìn)行插值模擬，得到當(dāng)日山仔水庫葉綠素a的空間分布，如圖3所示。可以看出，葉綠素a含量在空間分布上具有一定的分異性?？傮w上看，進(jìn)口的皇帝洞和七里葉綠素a含量高，庫心次之，出口壩前葉綠素a含量最低。對(duì)照4.2中，7月七里和皇帝洞的浮游植物細(xì)胞濃度均為107個(gè)/L數(shù)量級(jí)，而綜合營(yíng)養(yǎng)狀態(tài)均為輕度富營(yíng)養(yǎng)級(jí)別，葉綠素a含量最低的壩前浮游植物細(xì)胞濃度及綜合營(yíng)養(yǎng)狀態(tài)指數(shù)也均最低，表明空間插值法得到的葉綠素a含量空間分布較好地吻合了庫區(qū)各點(diǎn)位富營(yíng)養(yǎng)化的狀態(tài)。

圖3 山仔水庫7月葉綠素a含量空間分布

5 結(jié)論

（1）監(jiān)測(cè)期間（2016年7～9月），共鑒定出7個(gè)門41屬57種浮游植物。其中，綠藻門種數(shù)最多，硅藻門和藍(lán)藻門次之；藍(lán)藻門種群數(shù)量占絕對(duì)優(yōu)勢(shì)，優(yōu)勢(shì)種為微囊藻和銅綠微囊藻。

（2）浮游植物細(xì)胞密度指示法和TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)評(píng)價(jià)法對(duì)庫區(qū)的營(yíng)養(yǎng)狀況評(píng)價(jià)存在差異，但庫區(qū)各監(jiān)測(cè)點(diǎn)位的富營(yíng)養(yǎng)化程度總體趨勢(shì)均為：皇帝洞>七里>日溪>庫心>壩前。TLI(∑)指數(shù)顯示夏季庫區(qū)已出現(xiàn)輕度富營(yíng)養(yǎng)化情況，需加大藻類水華關(guān)注力度。

（3）各監(jiān)測(cè)點(diǎn)位葉綠素a含量總體趨勢(shì)為：七里>皇帝洞>庫心>日溪>壩前，較好地吻合了庫區(qū)各點(diǎn)位富營(yíng)養(yǎng)化狀態(tài)的分布?？臻g插值得到的葉綠素a含量空間分布比單純采樣數(shù)據(jù)更直觀、有效，便于分析和研究，但如何更精確地確定插值的系列參數(shù)，使得插值得出的預(yù)測(cè)表面與真實(shí)值更為接近，如何更好地應(yīng)用地學(xué)統(tǒng)計(jì)方法，更精確地描述葉綠素a等水質(zhì)參數(shù)的空間分布狀況，為水華預(yù)警和防治提供科學(xué)參考，需要進(jìn)一步的探究。

（4）采用克里金空間插值得到葉綠素a空間分布結(jié)果與TLI(∑)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)結(jié)果都顯示，2016年7月七里和皇帝洞、2016年9月皇帝洞水質(zhì)為輕度富營(yíng)養(yǎng)化，其它監(jiān)測(cè)時(shí)間水體為中營(yíng)養(yǎng)。

[1] 蘇玉萍,林慧,鐘厚璋,等. 富營(yíng)養(yǎng)化山仔水庫沉積物微囊藻復(fù)蘇的受控因子[J]. 生態(tài)學(xué)報(bào), 2011,31(20):6167-6173.

[2] 蘇玉萍,陳娜蓉,林婉珍,等. 福建省山仔水庫浮游植物特征與水體富營(yíng)養(yǎng)狀況分析 [J]. 亞熱帶資源與環(huán)境學(xué)報(bào),2006,1(2):48-54.

[3] 江源,彭秋志,廖劍宇,等. 浮游藻類與河流生境關(guān)系研究進(jìn)展與展望[J]. 資料科學(xué),2013，35(3):461-471.

[4] 吳潔,錢天鳴,虞左明. 西湖葉綠素 a 周年動(dòng)態(tài)變化及藻類增長(zhǎng)潛力試驗(yàn) [J]. 湖泊科學(xué)，2001,13（2）:143-148.

[5] Carlson R E. A trophic state index for lakes [J]. Limnol Oceanogr , 1977, 22 ( 2):361-369.

[6] Aizaki M. Application of modified carlsons trophic state index to Japanese lakes and its relationships to other parameters related to trophic state [J]. Res Rep Natl Inst Environ Study , 1981( 23):13-31.

[7] Thomas J R, Victor J B. A preliminary modeling analysis of water quality in Lake Okeechobee, Florida:calibration results [ J] .Water Research , 1995 , 29 ( 12):2755 -2766.

[8]王朝暉,韓博平,胡韌,等. 廣東省典型水庫浮游植物群落特征與富營(yíng)養(yǎng)化研究[J]. 生態(tài)學(xué)雜志. 2005, 24(4):402-405.

[9]周楊,許振成,蔡立梅,等. 高州水庫水體富營(yíng)養(yǎng)化評(píng)價(jià)與控制對(duì)策研究[J] .環(huán)境科學(xué)與管理, 2011, 36(2):179-183.

[10] 黃詳飛,陳偉民,蔡啟銘. 湖泊生態(tài)調(diào)查觀測(cè)與分析[M]. 北京：中國(guó)標(biāo)準(zhǔn)出版社,1999.

[11] 國(guó)家環(huán)境保護(hù)總局. 水和廢水監(jiān)測(cè)分析方法[M].4版. 北京：中國(guó)環(huán)境科學(xué)出版社,2002.

[12] 劉瑞民,王學(xué)軍,王翠紅,等. 應(yīng)用地統(tǒng)計(jì)學(xué)方法研究湖泊中葉綠素a的空間分布[J]. 農(nóng)業(yè)環(huán)境保護(hù),2001(5): 308 - 310.

[13] 胡光道,王仁鐸. 線性地質(zhì)統(tǒng)計(jì)學(xué)[M]. 北京: 地質(zhì)出版社, 1988.