宋云龍， 張金松*， 朱 佳， 高靜思， 王 麗， 陶 益， 張 麗

1．哈爾濱工業(yè)大學(xué)深圳研究生院土木與環(huán)境工程學(xué)院， 廣東 深圳 518055 2．深圳職業(yè)技術(shù)學(xué)院建筑與環(huán)境工程學(xué)院， 廣東 深圳 518055 3．深圳市水質(zhì)檢測中心， 廣東 深圳 518055 4．清華大學(xué)深圳研究生院， 廣東 深圳 518055 5．深圳市國家氣候觀象臺， 廣東 深圳 518040

西麗水庫葉綠素a時空分布特征及其影響因素

宋云龍1， 張金松1*， 朱 佳2， 高靜思1， 王 麗3， 陶 益4， 張 麗5

1．哈爾濱工業(yè)大學(xué)深圳研究生院土木與環(huán)境工程學(xué)院， 廣東 深圳 518055 2．深圳職業(yè)技術(shù)學(xué)院建筑與環(huán)境工程學(xué)院， 廣東 深圳 518055 3．深圳市水質(zhì)檢測中心， 廣東 深圳 518055 4．清華大學(xué)深圳研究生院， 廣東 深圳 518055 5．深圳市國家氣候觀象臺， 廣東 深圳 518040

以典型亞熱帶供水型水庫——西麗水庫為研究對象，對其2013年4月—2014年6月藻類ρ(Chla)及其相關(guān)環(huán)境因子進(jìn)行監(jiān)測，分析了藻類ρ(Chla)時空演變規(guī)律及其與環(huán)境因子的關(guān)系．結(jié)果顯示：西麗水庫ρ(TChla)(TChla為總?cè)~綠素a)為2．65～60．35 μgL，其中藍(lán)藻ρ(Chla)為0．77～30．58 μgL，硅藻ρ(Chla)為0．77～31．62 μgL，綠藻ρ(Chla)為0．77～12．71 μgL．汛期(6—9月)優(yōu)勢藻為藍(lán)藻，其他時期硅藻占優(yōu)勢，綠藻全年處于較低水平．藻類分布具有明顯的空間異質(zhì)性，ρ(Chla)整體呈西北庫區(qū)高東南庫區(qū)低的特點，汛期白芒河入庫區(qū)域ρ(Chla)偏高，存在一定的藻華風(fēng)險．汛期盛行西南風(fēng)，平均風(fēng)速為2．7 ms，有利于藻類自高濃度的西部庫區(qū)向東遷移，一定程度上降低了藻華風(fēng)險．水庫流場空間差異較大，水體流速為0．005～0．025 ms，水庫主導(dǎo)流向為東北—西南，對西北庫區(qū)的高濃度藻類起到稀釋作用．Pearson相關(guān)分析和典范對應(yīng)分析(CCA)顯示，溫度、有機(jī)物以及徑流量是影響藻類生長和分布的前三位因子，磷是藻類生長的限制性營養(yǎng)鹽．

西麗水庫； 葉綠素a； 時空分布； 環(huán)境因子； 相關(guān)分析； 典范對應(yīng)分析

水體富營養(yǎng)化引起的藻華是我國面臨的日益嚴(yán)重的水環(huán)境問題．目前藻華的研究集中在藻華嚴(yán)重的太湖、巢湖、滇池和三峽水庫等大型湖庫[1- 4]，對藻類濃度相對較低但敏感度更高的飲用水水源水庫關(guān)注較少．藻類生長和演替受水質(zhì)、氣象和水動力等多種因素共同影響[5- 6]，現(xiàn)有研究多關(guān)注其中某一方面，對多因素的綜合影響研究較少．對藻類集中爆發(fā)的春夏時期研究較多[7- 8]，對藻類風(fēng)險較低的秋冬季節(jié)研究不足，尤其是對藻類生物量年際演變規(guī)律的研究更少．藻類爆發(fā)的本質(zhì)是其生物量長期演變和聚集的最終體現(xiàn)，忽略整個演變過程而只關(guān)注藻類爆發(fā)時段則難以全面認(rèn)識和深入理解藻華機(jī)理．

傳統(tǒng)藻類研究技術(shù)手段主要有分光光度法和顯微鏡計數(shù)法兩大類．前者只能檢測總ρ(Chla)，對藻類群落結(jié)構(gòu)揭示能力不充分[9]．后者可將不同藻類鑒別開來，但由于不同種類的藻細(xì)胞體積相差可達(dá)幾百倍，顯微鏡計數(shù)法反映藻類生物量的能力欠佳．調(diào)制熒光儀根據(jù)藍(lán)藻、綠藻和硅藻等的葉綠素a熒光特性差異，可分別測得相應(yīng)的ρ(Chla)，使研究者在掌握藻類生物量的同時更深入的探究藻類群落演替規(guī)律[10- 12]．目前關(guān)于西麗水庫的研究相對較少，并且多關(guān)注水質(zhì)方面，藻類與環(huán)境因子關(guān)系的研究比較缺乏．因此課題組對西麗水庫ρ(Chla)及其相關(guān)環(huán)境因子進(jìn)行了15個月的連續(xù)監(jiān)測，研究了藻類爆發(fā)、消退、休眠、復(fù)蘇和再爆發(fā)全過程ρ(Chla)時空演變規(guī)律，識別了藻類高風(fēng)險區(qū)域，分析了藻類與其關(guān)鍵環(huán)境因子的相關(guān)關(guān)系，為預(yù)防藻類爆發(fā)，維護(hù)水生態(tài)安全提供理論依據(jù)和數(shù)據(jù)支撐．

1 材料與方法

1．1 采樣點的布設(shè)

1．1．1 采樣點

對西麗水庫10個采樣點進(jìn)行每月1～2次采樣監(jiān)測，采樣點布設(shè)見圖1，其中1號采樣點靠近麻磡河，2號采樣點靠近東江來水口，3號采樣點靠近西麗湖壩，4、5號采樣點為庫心活水區(qū)，6號采樣點靠近自來水廠取水口，7號采樣點為死水區(qū)，8號采樣點靠近農(nóng)業(yè)污染源，9號采樣點靠近鐵崗水庫取水口，10號采樣點靠近白芒河．

圖1 采樣點布設(shè)Fig．1 Location map of sample sites

1．1．2 采樣水層

每個點位分3～5層采樣，表層為水面下10 cm左右水層，透光層水樣在與透明度同深度處采取，溫躍層為水體1 m內(nèi)溫度變化超過1 ℃水層，活性層為ρ(DO)為2 mgL水層，底層為底泥上方約10 cm水層，其中表層、透光層和底層始終納入監(jiān)測范圍，溫躍層和活性層視WT(水溫)和ρ(DO)分布狀況而定，出現(xiàn)即監(jiān)測．

1．2 試驗方法

采用ZPY- 1型采水器分層采樣并分別保存，水樣采集后2 h內(nèi)運(yùn)回實驗室，4 ℃保存，當(dāng)天完成ρ(Chla)的測定，ρ(TOC)、ρ(CODMn)、ρ(TP)、ρ(TN)、ρ(NH4+-N)、ρ(NO3--N)、ρ(Fe)和ρ(Silicate)等2 d內(nèi)完成檢測．ρ(Chla)采用調(diào)制熒光儀(德國WALZ Phyto- PAM)測定，并定期與《水和廢水監(jiān)測分析方法》[13]中的丙酮萃取分光光度法校正；WT、水深、pH、ρ(DO)、濁度采用多參數(shù)水質(zhì)測定儀(美國YSI 6600V2型)現(xiàn)場測定；透明度采用塞氏盤現(xiàn)場測定；ρ(CODMn)采用酸性高錳酸鉀法測定；ρ(TP)采用鉬酸銨分光光度法測定；ρ(Silicate)采用硅鉬藍(lán)分光光度法測定；ρ(TN)、ρ(NH4+- N)、ρ(NO3-- N)采用流動分析儀(法國Alliance- Futura)測定，ρ(TOC)采用總有機(jī)碳分析儀(GE- Siever 5310C)測定；水庫流場采用SonTek RiverSurveyor M9流場儀監(jiān)測．

1．3 時期劃分

根據(jù)氣象條件將全年劃分為5個時期：前汛期為2013年和2014年4—6月，后汛期為2013年7—9月，高溫?zé)o雨期為2013年10—11月，冬季枯水期為2013年12月—2014年2月，溫度躍升期為2014年3月．

1．4 數(shù)據(jù)分析

藻類與環(huán)境因子間的雙變量關(guān)系采用Pearson相關(guān)分析在SPSS 20．0上進(jìn)行，藻類與環(huán)境因子間的綜合相關(guān)關(guān)系采用CCA在XLSTAT 2010上進(jìn)行．藻類平面分布圖是根據(jù)各點位藻類總ρ(Chla)在相應(yīng)氣候期內(nèi)的多月平均值，通過克里金插值法，在ArcMap 10．2上生成，水庫流暢分布圖在Matlab．R2012上繪制．

2 結(jié)果與討論

2．1 葉綠素a時空分布特征

相對于太湖、巢湖和滇池等大型天然湖泊[14- 17]，西麗水庫ρ(Chla)處于較低范圍，表現(xiàn)出城市供水水源水庫的鮮明特征．由圖2可見，西麗水庫藻類ρ(TChla)平均值為16．84 μgL，最高為60．35 μgL，出現(xiàn)在2014年4月10號采樣點透光層；最低為1．45 μgL，出現(xiàn)在2013年11月采集的2號采樣點表層.監(jiān)測期內(nèi)ρ(Chla)呈現(xiàn)V型變化趨勢：從2013年4月開始快速增加，并于2013年前汛期達(dá)到全年最高水平(21．85 μgL)，后汛期繼續(xù)維持高濃度(17．95 μgL)，高溫?zé)o雨期濃度開始下降(13．86 μgL)，冬季枯水期跌至全年最低水平(7．01 μgL)，2014年溫度躍升期藻類開始復(fù)蘇(8．85 μgL)，至2014年前汛期又劇增至高濃度水平(24．82 μgL)，并且比去年同期高13．59%．2013年前汛期藻類濃度高于后汛期，是因為前汛期流域地表污染物濃度高，強(qiáng)降雨對流域地表的沖刷作用將大量的營養(yǎng)鹽帶入水庫[18]，促進(jìn)藻類快速生長，而后汛期地表徑流帶入的營養(yǎng)鹽相對較少．2013年7月在暴雨后采樣，ρ(Chla)明顯降低．究其原因：①連續(xù)陰雨天氣不利于藻類生長②暴雨對藻類濃度產(chǎn)生稀釋作用．

注: 各組柱由左至右采樣點編號為1～10.圖2 ρ(TChla)隨時間變化Fig．2 Annual trend in total chlorophyll-a

由圖3可見，西麗水庫主要藻類是硅藻和藍(lán)藻，綠藻次之．硅藻ρ(Chla)平均值為7．81 μgL，最高為31．62 μgL，出現(xiàn)在2014年4月的10號采樣點，月均最高值為2014年4月的18．44 μgL，月均最低值為2014年1月的3．56 μgL．藍(lán)藻ρ(Chla)平均值為7．34 μgL，最高值為30．58 μgL，出現(xiàn)在2013年6月的10號采樣點，月均最高值為2014年6月的13．15 μgL，月均最低值為2014年1月的1．73 μgL．綠藻ρ(Chla)平均值為3．31 μgL，最高值為12．71 μgL，出現(xiàn)在2013年6月的10號采樣點，月均最高值為2014年5月的7．81 μgL，月均最低值為2013年12月的0．51 μgL．總體上東北庫區(qū)(1、2、3號采樣點)硅藻比例明顯偏高，這與東江來水硅藻含量高有關(guān)．白芒河庫灣(8、9、10號采樣點)藍(lán)藻優(yōu)勢明顯，這可能是因為入庫河流帶來的高濃度營養(yǎng)物質(zhì)適宜藍(lán)藻大量繁殖．南部庫區(qū)(4、5、6、7號采樣點)硅藻和藍(lán)藻優(yōu)勢度相當(dāng)．

注: 各組柱由左至右采樣點編號為1～10.圖3 藍(lán)藻、 綠藻、 硅藻ρ(Chla)隨時間變化Fig．3 Annual trend in chlorophyll-a of cyonophyta， bacillariophyta and chlorophyte

圖4反映了藻類種群結(jié)構(gòu)變化規(guī)律，硅藻和藍(lán)藻交替占優(yōu)勢，而綠藻全年處于較低水平．2013年4月藍(lán)藻比例為25．25%，之后快速上升，6月達(dá)到38．64%，超越硅藻成為優(yōu)勢藻類，9月達(dá)到峰值50.62%；之后逐步降低，2013年12月—2014年2月為30%左右，處于最低水平；2014年3月起隨著氣溫升高，藍(lán)藻開始復(fù)蘇，并在前汛期迅速增加，5月超越硅藻，重新成為優(yōu)勢藻類．硅藻與藍(lán)藻存在激烈的競爭關(guān)系，其生長趨勢與藍(lán)藻相反：2013年4—6月比例為44．31%～65．11%，優(yōu)勢明顯；6—9月比例降為26．55%～35．33%，低于藍(lán)藻；10月比例迅速上升，再次成為優(yōu)勢藻類，并且于12月末達(dá)到峰值63．62%，并將優(yōu)勢一直維持到2014年4月．綠藻比例為7．03%～28．18%，全年處于較低水平，呈現(xiàn)出前汛期高而秋冬低的演變規(guī)律；與2013年相比，2014年綠藻比例明顯降低，水庫藻類多樣性降低．綜合可見硅藻從10月至來年4月占優(yōu)勢，藍(lán)藻從5月至9月優(yōu)勢明顯，而綠藻濃度常年處于較低水平．

圖4 藻類種群結(jié)構(gòu)變化Fig．4 Community structure of phytoplankton

2．2 藻類高風(fēng)險區(qū)域識別

綜合考慮水環(huán)境質(zhì)量和水廠對藻類的去除能力，認(rèn)為ρ(TChla)<20 μgL無藻華風(fēng)險，20 μgL≤ρ(TChla)<40 μgL存在輕度藻華風(fēng)險，40 μgL≤ρ(TChla)<60 μgL存在中度藻華風(fēng)險，ρ(TChla)≥60 μgL存在高藻華風(fēng)險．圖5反映了藻類ρ(TChla)在各氣候期的平面分布情況．西麗水庫ρ(TChla)整體呈現(xiàn)西北庫區(qū)高東南庫區(qū)低的特點．2號采樣點是東江來水進(jìn)入西麗水庫的位置，水質(zhì)較優(yōu)，硅藻優(yōu)勢明顯，是全庫ρ(TChla)最低處(10．25 μgL)，無藻華風(fēng)險．1號采樣點是麻磡河入庫區(qū)，徑流帶入的藻類和營養(yǎng)鹽增加了水庫的藻華風(fēng)險，但由于靠近2號采樣點，受到大量東江來水的稀釋作用，風(fēng)險有所緩解．3、4、5號采樣點位于水庫活水區(qū)，水體交換迅速，流速大，ρ(TChla)為10．74～12．38 μgL，無藻華風(fēng)險．6號采樣點位于水庫最南端，是自來水廠取水口，遠(yuǎn)離麻磡河和白芒河，水體流速高，ρ(TChla)平均值為12．41 μgL，無藻華風(fēng)險．10號采樣點是白芒河入庫區(qū)域，徑流帶入大量的營養(yǎng)鹽和藻類，使得其附近的整個庫灣藻類濃度一直處于全庫最高水平(26．32 μgL)，汛期表層最高ρ(TChla)達(dá)到60．35 μgL，藻華風(fēng)險極高，高溫?zé)o雨期也存在一定風(fēng)險，容易爆發(fā)局部藍(lán)藻水華，是需要密切關(guān)注的高風(fēng)險區(qū)域．8號采樣點附近存在荔枝園，部分肥料和農(nóng)藥等面源污染進(jìn)入水庫后引起氮磷濃度升高，為藻類生長提供了營養(yǎng)物質(zhì)，存在較高的藻類風(fēng)險．9號采樣點是鐵崗水庫取水口，該區(qū)域庫形狹長，水體交換緩慢，其ρ(TChla)處于較高水平(14．82 μgL)，汛期存在藻華風(fēng)險．從時間上看，高溫?zé)o雨期、冬季枯水期和溫度躍升期各采樣點均無藻華風(fēng)險，汛期尤其是前汛期8、9、10號采樣點不同程度存在藻華風(fēng)險．

圖5 不同氣候期內(nèi)ρ(TChla)平面分布Fig．5 Chlorophyll- a chart in different climatic period

2．3 庫內(nèi)藻類輸移規(guī)律分析

2．3．1 風(fēng)對藻類遷移的影響

水體表面風(fēng)速和風(fēng)向?qū)υ孱惖妮斠朴绊戯@著[19- 21]．研究認(rèn)為，當(dāng)風(fēng)速小于3 ms時，藻類在水面順風(fēng)漂流；當(dāng)風(fēng)速大于3 ms時，風(fēng)浪會將表面漂浮的藻類帶入水下，發(fā)生垂直方向上的遷移[22- 23].由表1可見，2013年主導(dǎo)風(fēng)向為東北風(fēng)，頻率為19%，平均風(fēng)速2．1 ms．2014年主導(dǎo)風(fēng)向為東北偏東風(fēng)和東北風(fēng)，頻率均為16%，平均風(fēng)速分別為2．3、2．2 ms．監(jiān)測期間深圳主導(dǎo)風(fēng)風(fēng)速小于3 ms，主要促使西麗水庫藻類沿水面漂散，導(dǎo)致西麗水庫藻類呈現(xiàn)西北庫區(qū)高東南庫區(qū)濃度低的分布特點．2013年9月—2014年4月主要盛行東北偏北風(fēng)或東北風(fēng)，在風(fēng)力作用下，藻類從低濃度的東北部向高濃度的西南部遷移，風(fēng)加劇了藻類聚集，但其間西麗水庫藻類處于低濃度狀態(tài)，因此不會引發(fā)藻華風(fēng)險．2013年6月、8月以及2014年5月、6月主導(dǎo)風(fēng)為西南風(fēng)，促進(jìn)了水庫西部的高濃度藻類向東部遷移擴(kuò)散，有利于疏散高濃度藻團(tuán)，降低藻華風(fēng)險．2013年4月、5月、7月以及2014年4月主導(dǎo)風(fēng)為東北風(fēng)或東北偏東風(fēng)，則會促進(jìn)藻類向西南庫區(qū)聚集，加劇藻華風(fēng)險．10號采樣點位于西麗水庫西北角的狹長庫灣區(qū)域，不論是東北風(fēng)還是西南風(fēng)，均不利于區(qū)域內(nèi)藻類輸散，爆發(fā)局部藻華的風(fēng)險較高．此外，深圳地處東南沿海，汛期時有臺風(fēng)出現(xiàn)，最大風(fēng)速為8．1～12．7 ms，極大可風(fēng)速可達(dá)14．5～24．3 ms，風(fēng)向為西南偏南風(fēng)或西南風(fēng)．臺風(fēng)引發(fā)巨大波浪，促使表層藻類在垂直方向上劇烈運(yùn)動，加上暴雨的沖刷作用，使得表層聚集的藻類快速回到中下層水體，降低藻類水華風(fēng)險．

表1 西麗水庫風(fēng)速、風(fēng)向統(tǒng)計分析

2．3．2 流場對藻類遷移的影響

西麗水庫是典型的調(diào)水供水型水庫，接納大量的東江水并轉(zhuǎn)供給下游的鐵崗水庫，水體更新速度快于天然湖泊．圖6是水庫流場分布圖，受人工的供水和調(diào)水的影響，水庫流場空間差異大：在東江來水口、鐵崗水庫取水口以及河流入水口處，水體流速達(dá)0．020～0．025 ms，而死水區(qū)(7號采樣點)水體流速為0．005 ms左右．與風(fēng)驅(qū)動為主導(dǎo)因素的太湖流場0．005～0．010 ms的流速[24]相比，西麗水庫水體流速明顯偏高．低水力停留時間和高水體流速是西麗水庫保持較低藻類濃度的重要原因．

西麗水庫主導(dǎo)流向是東江進(jìn)水口(2號采樣點)—西麗湖壩(3號采樣點)—庫心(4號采樣點)—西部庫區(qū)(8號采樣點)．優(yōu)質(zhì)的東江來水自東北流向西南，確保了東半庫區(qū)水質(zhì)優(yōu)良，對西半庫區(qū)高濃度藻類起到稀釋作用．次主導(dǎo)流向是東江進(jìn)水口(2號采樣點)—西麗湖壩(3號采樣點)—庫心(5號采樣點)—水廠取水口(6號采樣點)，次主導(dǎo)流向流程短，不與河流交匯，外源污染少，確保水廠取水口水質(zhì)良好．白芒河入庫區(qū)域是西麗水庫藻類濃度最高的區(qū)域，該區(qū)域水體流速高，水體交換強(qiáng)烈，有利于藻類的擴(kuò)散，降低了局部藻類爆發(fā)風(fēng)險．麻堪河來水中含有較高濃度的藻類，入庫后受到主導(dǎo)流場的影響，藻類向西南方向疏散．整體上看，西麗水庫水由藻類濃度低的東部庫區(qū)流向藻類濃度高的西部庫區(qū)，對西部和西北部庫區(qū)的高濃度藻類起到稀釋作用，但不利西部庫區(qū)高濃度的藻類向中東部擴(kuò)散．

注：箭頭是表示該點位水體流速的矢量，箭頭大小表示流速大小，箭頭方向表示流速方向．圖6 西麗水庫流場分布Fig．6 Flow field distribution of Xili Reservior

2．4ρ(Chla)與環(huán)境因子的關(guān)聯(lián)性分析

2．4．1ρ(Chla)與水質(zhì)因子雙變量相關(guān)分析

表2反映了ρ(Chla)與水質(zhì)因子的相關(guān)關(guān)系．由表2可見，不同藻類與環(huán)境因子相關(guān)關(guān)系存在顯著差異．總藻、硅藻和綠藻的主要相關(guān)因子一致，按相關(guān)系數(shù)絕對值大小排序依次為水溫>pH>ρ(NO3--N)>ρ(Silicate)>ρ(TN)，其中ρ(Chla)與ρ(Silicate)負(fù)相關(guān)，與其他環(huán)境因子為正相關(guān)．藍(lán)藻ρ(Chla)與多個水質(zhì)因子存在較高相關(guān)性，依次是水溫>ρ(TOC)>pH>ρ(CODMn)>ρ(Silicate)>ρ(Fe)>透明度>ρ(TN)>濁度．

水溫對于藍(lán)藻、硅藻和綠藻的ρ(Chla)均是首要影響因子，這與Lü等[25]的研究結(jié)果相一致．西麗水庫汛期水溫為23～32 ℃，與藍(lán)藻的最適宜溫度范圍(25～33 ℃)重疊，容易爆發(fā)藍(lán)藻水華，因此汛期要加強(qiáng)對水溫和藍(lán)藻的監(jiān)測．水溫與硅藻和綠藻的相關(guān)性小于藍(lán)藻，原因可能是藍(lán)藻在適宜的水溫中大量增殖，在與硅藻、綠藻的競爭中取得優(yōu)勢，一定程度上抑制了硅藻和綠藻的生長．

pH與3種藻類均存在一定相關(guān)性．西麗水庫pH為6．51～8．76，汛期和旱期平均pH分別為7．89和7．25，呈現(xiàn)汛期高而旱期低的特點，這與藍(lán)藻的變化規(guī)律相一致．常淳等[26]在于橋水庫也發(fā)現(xiàn)類似規(guī)律．Berge等[27]認(rèn)為，當(dāng)pH在6．5～8．5范圍內(nèi)波動時，藻類的生長速率基本不變．西麗水庫pH變化可能是藻類增殖的結(jié)果而非原因，因為藻類的光合作用和呼吸作用等生命活動對pH影響顯著[28- 29]．白天光照充足，光合作用明顯大于呼吸作用，吸收CO2使水體pH升高．夜間以呼吸作用為主，生成CO2使pH降低．水質(zhì)和藻類的監(jiān)測是在白天進(jìn)行，因此檢測到高藻期間水體pH明顯升高．

西麗水庫ρ(TN)和ρ(TP)平均值分別為1．41、0．023 mgL，TNTP〔即ρ(TN)ρ(TP)〕約為61．相關(guān)分析顯示，ρ(Chla)與ρ(TN)存在一定相關(guān)性，與ρ(TP)無明顯相關(guān)性．西麗水庫ρ(TN)略低于天然湖泊，而ρ(TP)僅為其15左右[30]，營養(yǎng)鹽處于較低水平．吳雅麗等[31]認(rèn)為，太湖春季藻類生長的ρ(TP)閾值為0．059 mgL，而夏季藻類生長旺盛，對TP需求量更高．與天然湖泊較高的藻類葉綠素a含量相比，西麗水庫藻類濃度明顯處于低水平，原因可能是ρ(TP)較低，在一定程度上限制了藻類大量生長．

ρ(TOC)和ρ(CODMn)與藍(lán)藻相關(guān)性較高，可能的原因有：①高濃度的有機(jī)物為藻類增殖提供充足的碳源；②藍(lán)藻本身即為有機(jī)物，作為TOC和COD的一部分被檢測出來．透明度與藍(lán)藻負(fù)相關(guān)，濁度與藍(lán)藻正相關(guān)，說明藍(lán)藻的爆發(fā)是水體濁度升高、透明度降低的重要原因．

表2 藻類ρ(Chla)與水質(zhì)因子Pearson相關(guān)系數(shù)

注：*表示在0．05水平上顯著相關(guān)；**表示在0．01水平上顯著相關(guān)．

2．4．2ρ(Chla)與水動力和氣象因子雙變量相關(guān)分析

同一天各點位的水動力因子和氣象因子數(shù)值相同，數(shù)據(jù)特征與水質(zhì)因子不同，因此單獨分析其與藻類ρ(Chla)的相關(guān)性．由表3可見，總藻、硅藻和綠藻的主要相關(guān)因子一致，依次為HRT>濕度>日照時數(shù)>氣壓，其中與HRT、濕度為正相關(guān)，與日照時數(shù)、氣壓為負(fù)相關(guān)．與藍(lán)藻相關(guān)的因子較多，依次為氣壓>濕度>氣溫>降水量>徑流量>HRT，其中與氣壓負(fù)相關(guān)，與其他因子正相關(guān)．HRT與3種藻類均存在較高的相關(guān)性，硅藻、綠藻與HRT相關(guān)性更高，并且水體流速與3種藻類均有一定的負(fù)相關(guān)性，說明優(yōu)化水利調(diào)度，加快水體流動性，降低HRT，是抑制藻類生長的有效途徑．氣壓與3種藻類成負(fù)相關(guān)性，尤其是與藍(lán)藻相關(guān)系數(shù)達(dá)到-0．692；濕度與3種藻類正相關(guān)，與藍(lán)藻相關(guān)系數(shù)高達(dá)0．624．目前還沒有氣壓、濕度與藻類生長存在直接因果關(guān)系的報道．降水量、徑流與藍(lán)藻相關(guān)系數(shù)分別為0．441和0．422，是因為汛期降雨對流域的沖刷，以地表徑流的形式將大量的藍(lán)藻和營養(yǎng)物質(zhì)帶入水庫[18]．

表3 藻類ρ(Chla)與水動力和氣象因子Pearson相關(guān)系數(shù)

Table 3 Pearson correlation coefficients between chlorophyll-a and hydrodynamicmeteorological factors

表3 藻類ρ(Chla)與水動力和氣象因子Pearson相關(guān)系數(shù)

水動力∕氣象因子ρ(Chla)總藻藍(lán)藻硅藻綠藻HRT0．423??0．366??0．423??0．423??水體流速-0．154?-0．179?-0．112-0．189?氣溫0．0490．472??0．0410．041氣壓-0．200??-0．692??-0．200??-0．201??日照時數(shù)-0．253??0．088-0．253??-0．253??濕度0．348??0．624??0．348??0．348??降水量0．0870．441??0．0830．083蒸發(fā)量-0．125-0．029??-0．126-0．124徑流量0．0510．422??0．0510．052風(fēng)速0．0250．138?0．0240．025

注：*表示在0．05水平上顯著相關(guān)，**表示在0．01水平上顯著相關(guān)．

2．4．3 葉綠素a與環(huán)境因子典范對應(yīng)分析

利用CCA研究了水質(zhì)、水動力和氣象因子的綜合作用對藻類生長的影響，結(jié)果見圖7．F1特征值為0．079，解釋率為80．66%．水溫與F1軸關(guān)系最為密切，其次是氣溫，這說明F1軸主要反映了溫度的變化．溫度是影響西麗水庫藻類生長和分布的首要因子，藍(lán)藻和綠藻位于F1軸正方向，硅藻位于F1軸負(fù)方向，這說明溫度的升高有利于藍(lán)藻和綠藻在水體中占優(yōu)勢，而對硅藻不利．即水溫升高后，藍(lán)藻和綠藻競爭力高于硅藻．氣壓與F1軸顯著負(fù)相關(guān)，可能是由氣壓與溫度呈顯著負(fù)相關(guān)這一氣象規(guī)律引起的．F2特征值為0．019，解釋率為19．34%．ρ(TN)、ρ(TP)和ρ(Silicate)與F2軸最相關(guān)，因此F2軸代表了水體中的氮、磷和硅等營養(yǎng)鹽成分．ρ(TN)位于F2軸負(fù)向而ρ(TP)和ρ(Silicate)位于正向，這說明西麗水庫由氮鹽豐富態(tài)向磷鹽豐富態(tài)轉(zhuǎn)變更有利于藻類生長，磷是西麗水庫的限制性營養(yǎng)鹽．ρ(TOC)和ρ(CODMn)代表的有機(jī)物是藻類的第二大類影響因子，水體有機(jī)污染程度高的時期或區(qū)域中藻類生長繁殖速度快．徑流量和降水量是藻類的第三大影響因子，這可能是因為白芒河和麻磡河等河流污染相對嚴(yán)重，前汛期徑流將流域內(nèi)的大量藻類和有機(jī)物帶入水庫．ρ(Fe)與F2軸呈明顯的負(fù)相關(guān)性，因為Fe是藻類光合作用和新陳代謝的必要元素，對藻類生物量和生長速率均有重要影響[32- 33]．日照時數(shù)、蒸發(fā)量、pH、ρ(NO3--N)和ρ(NH4+-N)對藻類的影響也比較微弱，ρ(DO)、HRT和透明度等與藻類沒有明顯相關(guān)性．

環(huán)境因子：DO—ρ(DO)；TOC—ρ(TOC)；COD—ρ(CODMn)；TP—ρ(TP)；TN—ρ(TN)； 氨氮—ρ(NH4+-N)；硝酸鹽氮—ρ(NO3--N)；鐵—ρ(Fe)；硅酸鹽—ρ(Silicate)．圖7 ρ(Chla)與環(huán)境因子CCA分析Fig．7 Canonical correspondence analysis of environmental factors and phytoplankton

3 結(jié)論

a)ρ(TChla)時間變化規(guī)律為前汛期>后汛期>高溫?zé)o雨期>溫度躍升期>冬季枯水期．硅藻在冬春季節(jié)占優(yōu)勢，藍(lán)藻在汛期占優(yōu)勢，綠藻全年處于低濃度水平．ρ(TChla)整體呈西北庫區(qū)高東南庫區(qū)低的分布規(guī)律，汛期白芒河入庫區(qū)域ρ(TChla)偏高(以藍(lán)藻為主)，是藻華高風(fēng)險區(qū)．8號采樣點由于農(nóng)業(yè)面源污染的進(jìn)入，在汛期也存在一定藻華風(fēng)險，其他點位藻華風(fēng)險較低．

c) Pearson相關(guān)分析和CCA結(jié)果基本一致，溫度是藻類生長和分布的第一影響因子，TOC和COD代表的有機(jī)物以及徑流量和降水量分列第二位和第三位，磷是藻類生長的限制性營養(yǎng)鹽．藍(lán)藻的關(guān)鍵環(huán)境因子與硅藻、綠藻存在一定差異．

[1] 姚緒姣，劉德富，楊正健，等．三峽水庫香溪河庫灣冬季甲藻水華生消機(jī)理初探[J]．環(huán)境科學(xué)研究，2012，25(6)：645- 651． YAO Xujiao，LIU Defu，YANG Zhengjian，etal．Preliminary studies on the mechanism of winter dinoflagellate bloom in Xiangxi Bay of the Three Gorges Reservoir[J]．Research of Environmental Sciences，2012，25(6)：645- 651．

[2] WANG P F，WANG Y G，WANG C H，etal．Ecological characteristics and environmental factors of phytoplankton during different seasons and in different parts of Taihu Lake[J]．Fundamental and Applied Limnology，2015，187(1)：33- 42．

[3] LI M，LIU G J，WU L，etal．Eutrophication and water quality characteristics of Western Chaohu Lake in winter via planktonic alga diversity analysis[J]．Journal of University of Science and Technology of China，2015，45(2)：180- 188．

[4] 劉勇，朱元榮，弓曉峰，等．滇池近代富營養(yǎng)化加劇過程的沉積記錄[J]．環(huán)境科學(xué)研究，2012，25(11)：1236- 1242． LIU Yong，ZHU Yuanrong，GONG Xiaofeng，etal．Sedimentary records of accelerated eutrophication in Dianchi Lake over the recent decades[J]．Research of Environmental Sciences，2012，25(11)：1236- 1242．

[5] KEITH D，RICHARD J G，PAUL T，etal．Harmful algal blooms：How strong is the evidence that nutrient ratios and forms influence their occurrence[J]．Estuarine，Coastal and Shelf Science，2012,115：399- 413．

[6] HALLEGREF G M．A review of harmful algal blooms and their apparent global increase[J]．Phycologia，1993，32(2)：79- 99．

[7] 李春華，葉春，張詠，等．太湖湖濱帶藻密度與水質(zhì)、風(fēng)作用的分布特征及相關(guān)關(guān)系[J]．環(huán)境科學(xué)研究，2013，26(12)：1290- 1300． LI Chunhua，YE Chun，ZHANG Yong，etal．Temporal and spatial distribution of algal density and its relationship with water quality and wind factor in the littoral zone of Lake Taihu[J]．Research of Environmental Sciences，2013，26(12)：1290- 1300．

[8] NAZEER M，NICHOL J E．Development and application of a remote sensing- based Chlorophyll- a concentration prediction model for complex coastal waters of Hong Kong[J]．Journal of Hydrology，2016，532：80- 89．

[9] PYO J，HA S，PACHEPSKY Y A，etal．Chlorophyll- a concentration estimation using three difference bio- optical algorithms，including a correction for the low- concentration range：the case of the Yiam reservoir，Korea[J]．Remote Sensing Letters，2016，7(5)：407- 416．

[10] SJOLLEMA S B，BEUSEKOM SAM，GEEST H G，etal．Laboratory algal bioassays using pam fluorometry：effects of test conditions on the determination of herbicide and field sample toxicity[J]．Environmental Toxicology and Chemistry，2014，33(5)：1017- 1022．

[11] MAJARREIS J M，WATSON S B，SMITH R H．Nutrient status and its assessment by pulse amplitude modulated(PAM)fluorometry of phytoplankton at sites in Lakes Erie and Ontario[J]．Canadian Journal of Fisheries and Aquatic Sciences，2014，71(12)：1840- 1851．

[12] EDWARDS M S，KIM K Y．Diurnal variation in relative photosynthetic performance in giant kelpMacrocystispyrifera(Phaeophyceae，Laminariales)at different depths as estimated using PAM fluorometry[J]．Aquatic Botany，2014，92(2)：119- 128．

[13] 國家環(huán)境保護(hù)總局．水和廢水監(jiān)測分析方法[M]．4版．北京：中國環(huán)境科學(xué)出版社，2002：722- 723．

[14] WANG P F，WANG Y G，WANG C H，etal．Ecological characteristics and environmental factors of phytoplankton during different seasons and in different parts of Taihu Lake[J]．Fundamental and Applied Limnology，2015，187(1)：33- 42．

[15] 姜霞，王書航，鐘立香，等．巢湖藻類生物量季節(jié)性變化特征[J]．環(huán)境科學(xué)，2010，31(9)：2056- 2062． JIANG Xia， WANG Shuhang， ZHONG Lixiang，etal．Seasonal variation characteristics of algae biomass in Chaohu Lake[J]．Environtal Science，2010，31(9)：2056- 2062．

[16] LI M，LIU G J，WU L，etal．Eutrophication and water quality characteristics of western Chaohu Lake in winter via planktonic alga diversity analysis[J]．Journal of University of Science and Technology of China，2015，45(2)：180- 188．

[17] MA X X，WANG Y A，F(xiàn)ENG S Q，etal．Vertical migration patterns of different phytoplankton species during a summer bloom in Dianchi Lake，China[J]．Environmental Earth Sciences，2015，74(5)：3805- 3814．

[18] 宗棟良，方俊峰，王依林．降雨及徑流對西麗水庫水質(zhì)的影響[J]．中國農(nóng)村水利水電，2006，26(2)：44- 46．

[19] NASELLI-FLORES L，PADISAK J．Blowing in the wind：how many roads can a phytoplanktont walk down? a synthesis on phytoplankton biogeography and spatial processes[J]．Hydrobiologia，2016，764(1)：303- 313．

[20] VIDAL J，RIGOSI A，HOYER A，etal．Spatial distribution of phytoplankton cells in small elongated lakes subject to weak diurnal wind forcing[J]．Aquatic Sciences，2014，76(1)：83- 99．

[21] NEVEUX J，LEFEBVRE JP，Le Gendre R，etal．Phytoplankton dynamics in the southern New Caledonian lagoon during a southeast trade winds event[J]．Journal of Marine Systems，2010，82(4)：230- 244．

[22] ZHOU J，QIN B Q，CASENAVE C，etal．Effects of wind wave turbulence on the phytoplankton community composition in large，shallow Lake Taihu[J]．Environmental Science and Pollution Research，2015，22(16)：12737- 12746．

[23] 白曉華，胡維平，胡志新，等．2004年夏季太湖梅梁灣席狀漂浮水華風(fēng)力漂移入灣量計算[J]．環(huán)境科學(xué)，2005，26(6)：57- 60． BAI Xiaohua，HU Weiping，HU Zhixin，etal．Importation of wind-driven drift of mat-like algae bloom into meiliang bay of taihu lake in 2004 summer[J]．Environmental Sciences，2005，26(6)：57- 60．

[24] 姜恒志，崔雷，石峰．風(fēng)場、地形和吞吐流對太湖流場影響的研究[J]．水力發(fā)電學(xué)報，2013，32(6)：165- 171． JIANG Hengzhi，CUI Lei，SHI Feng．Study on influences of wind field，topography and inflowoutflow on flow in lake Tai[J]．Journal of Hydroelectric Engineering，2013，32(6)：165- 171．

[25] LV H，YANG J，LIU L M，etal．Temperature and nutrients are significant drivers of seasonal shift in phytoplankton community from a drinking water reservoir，subtropical China[J]．Environmental Science and Pollution Research，2014，21(9)：5917- 5928．

[26] 常淳，馮平，孫冬梅，等．基于逐步聚類分析的水庫浮游藻類生長預(yù)測[J]．中國環(huán)境科學(xué)，2015，35(9)：2805- 2812． CHANG Chun，F(xiàn)ENG Ping，SUN Dongmei，etal．Prediction of the alga growth in a reservoir based on the stepwise cluster analysis[J]．China Environmental Science，2015，35(9)：2805- 2812．

[27] BERGE T，DAUGBJERG N，ANDERSEN B B，etal．Effect of lowered pH on marine phytoplankton growth rates[J]．Marine Ecology Progress Series，2010，416(9)：79- 91．

[28] 王慶安，黃時達(dá)，孫鐵珩．多藻淺水體中pH值和溶解氧協(xié)同周期性變化初探[J]．四川環(huán)境，2001，20(2)：4- 29． WANG Qingan，HUANG Shida，SUN Tieheng．Study on the coordinate periodic change and the relativity between pH and DO in shallow water with alga[J]．Sichuan Environment，2001，20(2)：4- 29．

[29] 張澎浪，孫承軍．地表水體中藻類的生長對pH值及溶解氧含量的影響[J]．中國環(huán)境監(jiān)測，2004，20(4)：49- 51． ZHANG Penglang，SUN Chengjun．The influence of algae growing on pH and DO in surface water[J]．Environmental Monitoring in China，2004，20(4)：49- 51．

[30] DENG J M，QIN B Q，PAERL H W，etal．Effects of nutrients，temperature and their interactions on spring phytoplankton community succession in Lake Taihu，China[J]．South African Journal of Botany，2012，79(79)：62- 70．

[31] 吳雅麗，許海，楊桂軍，等．太湖春季藻類生長的磷營養(yǎng)鹽閾值研究[J]．中國環(huán)境科學(xué)，2013，33(9)：1622- 1629． WU Yali，XU Hai，YANG Guijun，etal．Developing the critical phosphorus threshold for spring algal growth in Lake Taihu，China[J]．China Environmental Science，2013，33(9)：1622- 1629．

[32] 姚波，席北斗，胡春明，等．鐵限制對浮游植物生長和群落組成的影響研究綜述[J]．生態(tài)環(huán)境學(xué)報，2010，19(2)：459- 465． YAO Bo，XI Beidou，HU Chunming，etal．A review of iron limitation on the growth and competition of phytoplankton[J]．Ecology and Environmental Sciences，2010，19(2)：459- 465．

[33] SUGIE K，KUMA K，F(xiàn)UJITA S，etal．Importance of intracellular Fe pools on growth of marine diatoms by using unialgal cultures and on the Oyashio region phytoplankton community during spring[J]．Journal of Oceanography，2011，67(2)：183- 196．

Spatial and Temporal Characteristics of Chlorophyll-a and Influencing Factors in Xili Reservoir

SONG Yunlong1， ZHANG Jinsong1*， ZHU Jia2， GAO Jingsi1， WANG Li3， TAO Yi4， ZHANG Li5

1．School of Civil and Environment Engineering， Harbin Institute of Technology Shenzhen Graduate School， Shenzhen 518055， China 2．School of Architectural and Environmental Engineering， Shenzhen Polytechnic， Shenzhen 518055， China 3．Shenzhen Water Quality Testing Center， Shenzhen 518055， China 4．Graduate School at Shenzhen， Tsinghua University， Shenzhen 518055， China 5．Shenzhen National Climate Observatory， Shenzhen 518040， China

A fifteen-month field observation was conducted from April 2013 to June 2014 in Xili Reservoir，a public drinking water supply reservoir in Shenzhen．The spatial and temporal distribution of environmental factors and chlorophyll-a concentration were studied，and their relationship was analyzed．The results indicated that the average ofρ(TChla)was 16．84 μgL，ranging from 2．65-60．35 μgL．The annual averageρ(Chla)of Cyanophyta was 7．34 μgL，ranging from 0．77-30．58 μgL．The annual averageρ(Chla)of Bacillariophyta was 7．81 μgL，ranging from 0．77-31．62 μgL．The annual averageρ(Chla)of Vhlorophyta was 3．31 μgL，ranging from 0．77-12．71 μgL．Cyanophyta was the primary species in the rainy season，while Bacillariophyta was the primary species in the other months．There was obvious spatial heterogeneity in the community of phytoplankton at Xili Reservoir．Chlorophyll-a concentration decreased gradually from the northwest to the southeast．In the rainy season，ρ(Chla)near Baimang River was very high，and there was a risk of algal blooming．Wind blew at 2．7 ms from the southwest in the rainy season．As a result，phytoplankton moved in the wind′s direction，reducing the risk of algal blooming in the western reservoir．The spatial heterogeneity of flow field at Xili Reservoir was obvious．The water flow rate was 0．005-0．025 ms，and the leading flow was from the northeast to the southwest．High concentrations of algae in the western part of Xili Reservoir were diluted by the leading flow．The results of Pearson correlation analysis and canonical correspondence analysis revealed that water temperature(WT)，ρ(TOC)and rain capacity were the most significant factors for algal blooms．Phosphorus was the restricting factor of algae growth．

Xili Reservoir； chlorophyll-a； temporal and spatial variations； environmental factors； correlation analysis； canonical correspondence analysis

2016- 10- 02

2016- 11- 29

水利部公益性行業(yè)科研專項(201301047)；深圳科技創(chuàng)新項目基礎(chǔ)研究(JCYJ20140508155916418)

宋云龍(1983-)，男，山東威海人，songyunlongwater@163．com．

*責(zé)任作者，張金松(1963-)，男，黑龍江哈爾濱人，教授，博士，主要從事飲用水安全保障研究，zhangjinsong@waterchina．com

X524

1001- 6929(2017)04- 0510- 10

A

10．13198j．issn．1001- 6929．2017．01．70

宋云龍，張金松，朱佳，等．西麗水庫葉綠素a時空分布特征及其影響因素[J]．環(huán)境科學(xué)研究，2017，30(4)：510-519.

SONG Yunlong，ZHANG Jinsong，ZHU Jia，etal．Spatial and temporal characteristics of chlorophyll-a and influencing factors in Xili Reservoir[J]．Research of Environmental Sciences，2017，30(4)：510-519.