香溪河庫灣春季水華縱向分布對水層結(jié)構(gòu)的響應(yīng)

劉 流 劉德富 黃鈺鈴 王從鋒

（三峽大學(xué) 水利與環(huán)境學(xué)院，湖北 宜昌 443002）

春季水華在富營養(yǎng)化水體的普遍暴發(fā)是一個(gè)世界性難題，各國學(xué)者已做了大量研究.Sverdrup［1］在1953年提出了解釋海洋春季水華暴發(fā)的經(jīng)典臨界層理論，首次提出了水溫分層誘發(fā)春季水華的概念模型，并成功應(yīng)用于解釋湖泊和水庫中水華暴發(fā)和浮游植物群落結(jié)構(gòu)演替過程［2－4］.三峽水庫支流庫灣是典型的深水河道型庫灣，對其春季水華暴發(fā)的特點(diǎn)和誘因，前人也已開展了大量研究.易仲強(qiáng)［5］發(fā)現(xiàn)香溪河庫灣春夏期間存在顯著的水溫分層，徐耀陽［6］和王嵐［7］發(fā)現(xiàn)香溪河庫灣水華的暴發(fā)具有空間異質(zhì)性，并根據(jù)縱向差異將香溪河庫灣自下游往上游劃分為干流區(qū)、湖泊區(qū)、過渡區(qū)和河流區(qū)［8］.生態(tài)分區(qū)可以為支流水華的防控提供重要參考.王玲玲［9］得出庫灣沿程浮游植物生物量與上游來流流量呈顯著負(fù)相關(guān)，提出通過調(diào)節(jié)上游來流流量來控制支流水華.

然而，他們都忽略了香溪河庫灣分層異重流的存在［10］.分層異重流作為一個(gè)重要驅(qū)動力，極大地改變了庫灣水體的理化分布格局［5，11］，其在支流庫灣的潛入深度、潛入厚度和影響范圍將隨水溫分層的季節(jié)性變化和三峽水庫的水位調(diào)度而改變，反過來導(dǎo)致支流庫灣水溫分層的縱向差異，從而直接或間接影響水華暴發(fā)情勢［11］.因此，不能對香溪河庫灣進(jìn)行簡單的縱向分區(qū).為驗(yàn)證以上理論，本研究于2010年春季水華期間對香溪河庫灣的水動力、水溫、營養(yǎng)鹽等環(huán)境因子進(jìn)行了現(xiàn)場監(jiān)測，研究分層異重流對庫灣水溫分層的影響，并試圖解釋庫灣春季水華浮游植物縱向分布差異的主導(dǎo)因素.

1 材料與方法

香溪河是三峽水庫庫首區(qū)域最大的一條支流，流域總面積3 095km2，長94km，多年平均流量為47.4 m3／s［12］.它起源于位于湖北省西北部的神農(nóng)架國家森林公園，向南流經(jīng)興山縣和秭歸縣，在香溪鎮(zhèn)注入長江.2003年6月三峽水庫初次蓄水后，距香溪河與長江干流交匯點(diǎn)26km范圍以內(nèi)即被回水淹沒，形成一個(gè)典型的深水河道型庫灣.本研究于2010年4月18日對香溪河庫灣進(jìn)行了現(xiàn)場監(jiān)測，三峽水庫干流設(shè)置一對照點(diǎn)CJXX，庫灣設(shè)點(diǎn)XX00～XX10，在未被庫灣回水淹沒的天然河道設(shè)置一對照點(diǎn)XXYT，如圖1所示.

采用馬駿［13］提出的三峽水庫支流庫灣低流速條件下的測流方法測香溪河庫灣流速.獲得的流速數(shù)據(jù)包括3個(gè)方向：東（Ve），北（Vn），垂向（Vu）.水平向的流速由Ve和Vn投影在河流主泓線上合成，本研究忽略了垂向上的流速.用水質(zhì)多參數(shù)儀（Hydrolab，美國HACH公司）測量水深、水溫和電導(dǎo)率剖面.采用薩氏透明度板測量水體透明度，IL1400A水下光量子儀測量光合作用有效輻射（400～700nm）.現(xiàn)場取350mL水樣2瓶帶回實(shí)驗(yàn)室參考《水和廢水分析方法》（第4版）［14］測定水體總氮（TN）、總磷（TP）和溶解性硅酸鹽（D－Si）及葉綠素a（Chl.a）濃度.

圖1 香溪河庫灣采樣點(diǎn)分布圖

真光層深度（Zeu）一般定義為水下光合作用有效輻射為水面的1%處對應(yīng)的深度［15］.為簡化計(jì)算或彌補(bǔ)水下光強(qiáng)數(shù)據(jù)缺失，也可通過建立Zeu與水體透明度（SD）的關(guān)系式來計(jì)算Zeu.本文采用Zeu＝1.802 1SD＋1.818 2作為計(jì)算Zeu的公式［16］.水體表面混合層深度（Zmix）是水體上層垂直混合充分的水層.本文采用水溫垂向梯度不超過0.2℃／m時(shí)對應(yīng)的深度計(jì)算.

2 結(jié)果分析

2.1 分層異重流特性

三峽水庫干流水體的化學(xué)成分與支流庫灣存在顯著差別，電導(dǎo)率作為水體中溶解質(zhì)濃度的表征，可用來示蹤倒灌異重流的運(yùn)動過程.圖2給出了4月18日香溪河庫灣的異重流沿縱剖面分布（白色箭頭表示水流流入庫灣，黑色箭頭表示水流流出庫灣；背景為電導(dǎo)率等值線，白色表示高電導(dǎo)率，黑色表示低電導(dǎo)率）.整體來看，香溪河庫灣的流速極低.三峽水庫干流水體由表層30m處倒灌進(jìn)入香溪河庫灣，在距干支流交匯點(diǎn)約2km處表層異重流開始下潛為中層異重流，至距干支流交匯點(diǎn)約2km處異重流分為兩部分，大部分水體以中層異重流形式繼續(xù)下潛，小部分仍以表層異重流形式運(yùn)行直至庫灣回水末端.異重流兩部分厚度均沿程減小.上游來流以底部異重流形式潛入庫灣，在庫灣中游約20km處出現(xiàn)分叉，一部分形成中層異重流，另一部分仍以底部異重流形式流出庫灣.

圖2 香溪河庫灣異重流縱剖面分布

2.2 水溫及光學(xué)分布特性

異重流水體一般混合比較均勻，因此異重流潛入深度的水溫結(jié)構(gòu)會發(fā)生改變.香溪河庫灣水溫垂向分布如圖3所示.

圖3 香溪河庫灣水溫垂向分布

由圖3可以看出，三峽水庫干流水體基本不存在水溫分層（最大水溫垂向梯度僅為0.07℃／m），而香溪河庫灣下游XX00水溫分層極其微弱，表層異重流水體水溫幾乎均勻分布（最大水溫垂向梯度僅為0.09℃／m）；XX01～XX02河段呈現(xiàn)相同的垂向水溫結(jié)構(gòu)，由于異重流在該河段自中層潛入，異重流的存在導(dǎo)致水體沿垂向在表層出現(xiàn)強(qiáng)制性分層（XX01處10 m以下最大水溫垂向梯度僅為0.06℃／m）；XX03自10m以下水溫分層不顯著（最大水溫垂向梯度僅為0.14℃／m），10m以上的表層水溫呈現(xiàn)較大梯度（最大水溫垂向梯度為0.35℃／m），XX04出現(xiàn)同樣的垂向水溫結(jié)構(gòu)，水溫垂向梯度轉(zhuǎn)折點(diǎn)出現(xiàn)在15m；庫灣中游XX05～XX06河段受分層倒灌異重流影響，水溫沿垂向出現(xiàn)雙溫躍層結(jié)構(gòu)；XX07～XX08河段已不受中層倒灌異重流的影響，沿垂向呈現(xiàn)顯著的水溫分層（平均水溫垂向梯度為0.27℃／m），但受表層倒灌異重流影響，表層5m以上呈現(xiàn)不同的水溫梯度（5m以下平均水溫垂向梯度僅為0.11℃／m），庫灣上游XX09～XX10河段呈現(xiàn)同樣的水溫結(jié)構(gòu).受分層異重流影響，香溪河庫灣呈現(xiàn)獨(dú)特的水溫分層結(jié)構(gòu)，缺乏一般湖庫水溫分層出現(xiàn)的典型溫躍層.

海洋和湖泊水體混合層深度一般取決于水溫分層及風(fēng)的 季節(jié)性 變化［18－19］.盧 晨晨［20］對三峽庫區(qū)河陸風(fēng)規(guī)律的研究發(fā)現(xiàn)，三峽庫區(qū)地面平均風(fēng)速小于2 m／s；且香溪河庫灣地形狹窄，風(fēng)的最大傳播距離不超過2km.因此，風(fēng)對香溪河庫灣的Zmix影響不顯著，香溪河庫灣Zmix的變化主要取決于水溫分層的發(fā)育.圖4給出了2010年4月18日香溪河庫灣的Zmix縱向分布.可以看出，由于水溫分層不顯著，三峽水庫干流和干支流交匯點(diǎn)附近的Zmix很大，幾乎到底，自XX01至庫灣回水末端Zmix呈現(xiàn)近似V字型分布，庫灣中游XX05處Zmix達(dá)到11.97m.這是由于原有的水體混合層深度與上游異重流潛入深度相接，加深了水體混合層.圖4也給出了Zeu的縱向分布，Zeu范圍為2.72～8.49m，自庫灣下游往回水末端遞減.

圖4 香溪河庫灣Zmix和Zeu縱向分布

2.3 營養(yǎng)鹽分布特性

圖5給出了2010年4月18日三峽水庫干流和香溪河庫灣的主要營養(yǎng)鹽濃度縱向分布.可以看出，干支流TN濃度值整體較高，在0.68～1.90mg／L范圍內(nèi)波動；TN的縱向分布大致呈V字型，兩端高、中上游XX07～XX09河段較低，最大值出現(xiàn)在XX10；而庫灣上游來流TN濃度卻低于XX10.TP整體值也偏高，在0.09～0.53mg／L范圍內(nèi)波動；TP大致自下游往上游呈遞增趨勢，回水末端XXYT最高.

氮磷比是衡量浮游植物生長是否受到氮或磷相對限制的一個(gè)重要指標(biāo)，TN／TP＝16時(shí)浮游植物生長不受氮或磷的相對限制［21］.由圖5可以看出，干流及庫灣回水氮磷比均小于16，在5.21～13.56范圍內(nèi)波動，沿縱向大致呈遞減態(tài)勢，在上游回水末端略有上升.D－Si濃度在三峽水庫干流較高，而在香溪河庫灣中下游XX00～XX04河段均呈現(xiàn)較高濃度（大于7.00mg／L），自XX05至庫灣回水末端整體呈U字型分布，D－Si濃度在0.16～5.10mg／L范圍內(nèi)波動，中上游XX06～XX08河段D－Si最低（均小于0.5 mg／L）；庫灣上游來流D－Si濃度（7.00mg／L）與中下游河段處于同一水平.

圖5 三峽水庫干流和香溪河庫灣主要營養(yǎng)鹽縱向分布

可以看出，三峽水庫干支流的TN、TP濃度整體均較高，滿足浮游植物生長需要，而其空間分布差異性較大.正如吉小盼［11］指出的，倒灌異重流對TN逆向補(bǔ)給效應(yīng)顯著.受昭君鎮(zhèn)生活污水排放影響，回水末端XX10河段TN濃度顯著高于中上游和上游來流.而由于香溪河庫灣流域?yàn)楦涣椎V區(qū)，上游來流TP濃度極高（0.56mg／L），TP在庫灣的縱向分布受上游來流影響顯著.氮磷比的縱向分布顯示，整個(gè)干支流均存在潛在的氮限制，而在庫灣中上游尤為顯著.上游來流與水庫干流D－Si濃度基本持平，中上游河段低D－Si濃度表明硅藻水華繁殖消耗了大量溶解性硅酸鹽.

3 討 論

圖6給出了干支流Zeu／Zmix和Chl.a的縱向分布，Zeu／Zmix是一個(gè)判別光照在垂向的分布和水體在垂向上的摻混對藻類生長的影響的綜合指標(biāo)［22］，Chl.a常被用來表征浮游植物生物量［23］.可以看出，Zeu／Zmix在水庫干流及干支流交匯點(diǎn)附近均很?。?.12），下游河段 XX01～XX04呈現(xiàn)較大值，其中XX01、XX02均大于11；庫灣中游XX05處出現(xiàn)Zeu／Zmix極小值（0.38），中上游河段XX06～XX10的Zeu／Zmix值在1.41～1.88范圍內(nèi)波動.Chl.a濃度的縱向分布存在顯著差異，三峽水庫干流Chl.a濃度極低，香溪河庫灣中下游河段XX00～XX04濃度值范圍為0.94～2.10mg／m3，而自 XX05至回水末端 Chl.a濃度沿程劇增，最大值高達(dá)161.41mg／m3，暴發(fā)了嚴(yán)重的水華.

圖6 三峽水庫干流和香溪河庫灣

關(guān)于浮游植物生物量對水溫分層的響應(yīng)，已有學(xué)者做了大量研究.在營養(yǎng)鹽、水溫等因子滿足浮游植物生長條件下，若水體出現(xiàn)較強(qiáng)溫度分層，Zeu／Zmix＞1時(shí)，浮游植物生長不受光照限制，這時(shí)候就會暴發(fā)高強(qiáng)度水華；反之，若Zeu／Zmix＜1，即在水溫分層不顯著時(shí)，浮游植物顆粒就會被摻混至真光層以下，其接受光照的機(jī)會大大減小，生長就受到光限制［24］.但若Zeu／Zmix?1，即Zmix極小，營養(yǎng)鹽受水溫分層阻隔無法進(jìn)入水體混合層，Zmix內(nèi)的營養(yǎng)鹽會很快消耗殆盡，無垂向遷移能力的浮游植物生長會受到營養(yǎng)鹽限制［25］；與之相反，Zeu／Zmix存在一個(gè)閾值，約處于0.2～0.35［24］之間（春季取0.35），小于這個(gè)閾值，由于浮游植物生物量的積累小于損失，則不可能暴發(fā)水華.因此，Zeu／Zmix＝1是水華暴發(fā)的最有利條件［25］.

香溪河庫灣的營養(yǎng)鹽條件對水華暴發(fā)無絕對限制，只在庫灣中上游存在潛在的氮、硅限制，然而中上游卻暴發(fā)了高強(qiáng)度的水華，顯然無法用王嵐［8－9］對支流庫灣的分區(qū)理論來解釋.而庫灣Zeu／Zmix的縱向分布差異可以很好的解釋這一現(xiàn)象，三峽水庫干流和干支流交匯點(diǎn)附近Zeu／Zmix＜0.35，因此Chl.a濃度極低；而XX01～XX04河段Zeu／Zmix?1，由圖7可以看出，該河段以硅藻、綠藻等無浮力或中性浮力藻類為主，這些藻類在高強(qiáng)度水溫分層（Zmix很?。l件下，極易發(fā)生沉降損失，因此Chl.a濃度也極低；庫灣中游XX05點(diǎn)Zeu／Zmix＞0.35，之后中上游XX06～XX10河段Zeu／Zmix＞1，暴發(fā)了高強(qiáng)度的硅藻水華.圖7顯示，甲藻只在香溪河庫灣干支流交匯點(diǎn)附近XX00和回水末端XX10出現(xiàn)，這與甲藻喜高摻混水流密切相關(guān)；綠藻在三峽水庫干流至庫灣中游均存在，而硅藻主要出現(xiàn)在XX04～XX10，并在XX06～XX10河段占優(yōu).這也解釋了圖5中庫灣中上游D－Si濃度偏低的現(xiàn)象.

圖7 三峽水庫干流和香溪河庫灣浮游植物群落結(jié)構(gòu)縱向分布

而由圖4可知，Zeu／Zmix值主要由Zmix決定，因此，香溪河庫灣水華情勢縱向分布的差異主要由水溫分層的差異決定，也在很大程度上受分層異重流的潛入深度和運(yùn)行距離影響.而本文限于篇幅，只研究了香溪河庫灣的一次春季水華，因此，未來有必要對香溪河庫灣的分層異重流運(yùn)動規(guī)律和浮游植物分布及演替過程進(jìn)行長系列觀測，探討不同水庫運(yùn)行期及水庫調(diào)度過程下浮游植物縱向分布及群落演替規(guī)律，以建立水庫調(diào)度與支流庫灣水華縱向分布的響應(yīng)關(guān)系，為三峽水庫的生態(tài)調(diào)度提供支撐.

4 結(jié) 論

1）春季香溪河庫灣出現(xiàn)顯著的分層異重流，水庫干流水體以表中層異重流形式倒灌進(jìn)入庫灣，而上游來流首先以底部異重流形式潛入庫灣，隨后形成中層異重流.

2）分層異重流極大地改變了香溪河庫灣的水溫結(jié)構(gòu)，水庫干流水溫不分層，而支流香溪河庫灣中受復(fù)雜的分層異重流影響出現(xiàn)獨(dú)特的水溫垂向分層結(jié)構(gòu)，缺乏典型的溫躍層，且在中游出現(xiàn)雙溫躍層形式.

3）三峽水庫干流和香溪河庫灣水華暴發(fā)均不存在營養(yǎng)鹽絕對限制，水庫干流通過倒灌異重流對庫灣的TN進(jìn)行補(bǔ)給，而庫灣上游來流是庫灣高TP濃度的重要原因，庫灣中上游受氮相對限制較顯著.

4）香溪河庫灣春季以硅藻水華為主，受水溫分層影響，香溪河庫灣Zeu／Zmix值沿縱向分布呈現(xiàn)極大差異性，從而導(dǎo)致浮游植物生物量和群落結(jié)構(gòu)的縱向分布差異，而分層異重流導(dǎo)致這種差異的根本驅(qū)動力，因此，弄清異重流運(yùn)行規(guī)律及特性對庫灣水華的研究至關(guān)重要.

［1］ Sverdrup H.On Conditions for the Vernal Blooming of Phytoplankton［J］.ICES Journal of Marine Science，1953，18（3）：287－295.

［2］ Ganf G，Oliver R.Vertical Separation of Light and A－vailable Nutrients as a Factor Causing Replacement of Green Algae by Blue－Green Algae in the Plankton of a Stratified Lake［J］.The Journal of Ecology，1982，70（3）：829－844.

［3］ Berger S A，Diehl S，Stibor H，et al.Water Temperature and Mixing Depth Affect timing and Magnitude of Events During Spring Succession of the Plankton［J］.Oecologia，2007，150（4）：643－654.

［4］ Lindenschmidt K E，Chorus I.The Effect of Water Column Mixing on Phytoplankton Succession，Diversity and Similarity［J］.Journal of Plankton Research，1998，20（10）：1927－1951.

［5］ 易仲強(qiáng)，劉德富，楊正健，等.三峽水庫香溪河庫灣水溫結(jié)構(gòu)及其對春季水華的影響［J］.水生態(tài)學(xué)雜志，2009，2（5）：6－11.

［6］ Xu Y，Cai Q，Ye L，et al.Asynchrony of Spring Phytoplankton Response to Temperature Driver Within a Spatial Heterogeneity Bay of Three Gorges Reservoir，China［J］.Limnologica－Ecology and Management of Inland Waters，2011，41（3）：174－180.

［7］ Wang L，Cai Q，Tan L，et al.Longitudinal Differences of Phytoplankton Community during a Period of Small Water Level Fluctuations in a Subtropical Reservoir Bay（Xiangxi Bay，Three Gorges Reservoir，China）［J］.International Review of Hydrobiology，2011，96（4）：381－396.

［8］ Wang L，Cai Q，Zhang M，et al.Longitudinal Patterns of Phytoplankton Distribution in a Tributary Bay under Reservoir Operation［J］. Quaternary International，2011，244（2）：280－288.

［9］ 王玲玲，戴會超，蔡慶華.香溪河生態(tài)調(diào)度方案的數(shù)值模擬［J］.華中科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，37（4）：111－114.

［10］紀(jì)道斌，劉德富，楊正健，等.三峽水庫香溪河庫灣水動力特性分析［J］.中國科學(xué)：G輯，2010，40（1）：101－112.

［11］吉小盼，劉德富，黃鈺鈴，等.三峽水庫泄水期香溪河庫灣營養(yǎng)鹽動態(tài)及干流逆向影響［J］.環(huán)境工程學(xué)報(bào)，2010，4（12）：2687－2693.

［12］黃真理，李玉梁.三峽水庫水質(zhì)預(yù)測和環(huán)境容量計(jì)算［M］.北京：中國水利水電出版社，2006.

［13］馬駿，劉德富，紀(jì)道斌，等.三峽水庫支流庫灣低流速條件下測流方法探討及應(yīng)用［J］.長江科學(xué)院院報(bào)，2011，28（6）：30－34.

［14］王心芳，魏復(fù)盛，齊文啟.水和廢水監(jiān)測分析方法［M］.4版.北京：中國環(huán)境科學(xué)出版社，2002.

［15］Kirk J T O.Light and Photosynthesis in Aquatic Eco－systems［M］.Cambridge University Press，1994.

［16］Liu L，Liu D F，Johnson D M，et al.Effects of Vertical Mixing on Phytoplankton Blooms in Xiangxi Bay of Three Gorges Reservoir：Implications for Management［J］.Water Research，2012，46（7）：2121－2130.

［17］Lawrence I，Bormans M，Oliver R，et al.Physical and Nutrient Factors Controlling Algal Succession and Biomass in Burrinjuck Reservoir［R］.Land ＆ Water Resources Australia，2000.

［18］Kling G W.Comparative Transparency，Depth of Mixing，and Stability of Stratification in Lakes of Cameroon，West Africa［J］.Limnology and Oceanography，1988，27－40.

［19］Farmer D M，Carmack E.Wind Mixing and Restratification in a Lake Near the Temperature of Maximum Density［J］.Journal of Physical Oceanography，1981，11（11）：1516－1533.

［20］魯晨晨.三峽壩區(qū)河陸風(fēng)的理論與觀測研究［D］.南京：南京信息工程大學(xué)，2009：5.

［21］Goldman J C，Mccarthy J J，Peavey D G.Growth Rate Influence on the Chemical Composition of Phytoplankton in Oceanic Waters［J］.1979，279：210－215.

［22］Jensen J，Jeppesen E，Olrik K，et al.Impact of Nutrients and Physical Factors on the Shift from Cyanobacte－rial to Chlorophyte Dominance in Shallow Danish Lakes［J］.Canadian Journal of Fisheries and Aquatic Sciences，1994，51（8）：1692－1699.

［23］Boyer J N，Kelble C R，Ortner P B，et al.Phytoplankton Bloom Status：Chlorophyll a Biomass as an Indicator of Water Quality Condition in the Southern Estuaries of Florida，USA［J］.Ecological Indicators，2009，9（6）：S56－S67.

［24］Oliver R，Hart B，Olley J，et al.The Darling River：Algal Growth and the Cycling and Sources of Nutrients［R］.Murray Darling Basin Commission，Project M386 1999，2000.

［25］Reynolds C.Succession and Vertical Distribution of Phytoplankton in Response to Thermal Stratification in a Lowland Mere，with Special Reference to Nutrient A－vailability［J］.The Journal of Ecology，1976，64（2）：529－551.

［26］Lawrence I，Bormans M，Oliver R，et al.Factors Controlling Algal Growth and Composition in Reservoirs：Report of Reservoir Management Workshops January 2000［R］.LWRRDC，Canberra，2000.