徐雅倩， 馬 駿， 楊正健

(湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院/ 河湖生態(tài)修復(fù)與藻類(lèi)利用湖北省重點(diǎn)實(shí)驗(yàn)室， 湖北 武漢 430068)

水華暴發(fā)是一個(gè)復(fù)雜的生態(tài)動(dòng)力學(xué)過(guò)程[1-3],國(guó)內(nèi)外學(xué)者從水動(dòng)力[4-6]、水溫[7-9]、光照[10-13]、營(yíng)養(yǎng)鹽[14-18]等方面開(kāi)展了一系列野外觀(guān)測(cè)和現(xiàn)場(chǎng)圍隔實(shí)驗(yàn),取得了諸多進(jìn)展。然而,對(duì)水庫(kù)或水質(zhì)管理部門(mén)來(lái)說(shuō),僅在水華發(fā)生后跟蹤監(jiān)測(cè)遠(yuǎn)遠(yuǎn)不夠,更需要能在水華暴發(fā)前較準(zhǔn)確地預(yù)測(cè)出水華情勢(shì),提前采取干預(yù)措施。水華的數(shù)值模擬是解決這一問(wèn)題的有效手段:一方面可以預(yù)測(cè)未來(lái)一段時(shí)間內(nèi)的水華情勢(shì);另一方面也可以用于水華暴發(fā)機(jī)理研究。

國(guó)內(nèi)外學(xué)者對(duì)水華的數(shù)值模擬研究經(jīng)歷了從黑箱模型、統(tǒng)計(jì)模型到過(guò)程模型研究的過(guò)程。目前水華模型以過(guò)程模型為主,模擬精度主要取決于關(guān)鍵生態(tài)動(dòng)力學(xué)參數(shù)的取值。藻類(lèi)生態(tài)動(dòng)力學(xué)參數(shù)取值不僅與藻種有關(guān),還與環(huán)境水體物化特性有關(guān)[19-20],即使同一藻種,在不同生長(zhǎng)階段參數(shù)取值也不一樣[21]?；谂R界層理論[22],垂向梯度差異(包括溫度、光照和濃度等方面)對(duì)于藻類(lèi)長(zhǎng)勢(shì)有著決定性影響。由美國(guó)陸軍工程兵團(tuán)和波特蘭州立大學(xué)開(kāi)發(fā)的CE-QUAL-W2模型是應(yīng)用于河流、河口、湖泊和水庫(kù)的立面(縱向/垂向)二維水動(dòng)力水質(zhì)模型。該模型假定水體橫向均勻,在相對(duì)較狹窄的、縱向和垂向上存在溫度或濃度梯度的水體都有廣泛應(yīng)用[21]。CE-QUAL-W2模型在水流、水溫和營(yíng)養(yǎng)鹽方面模擬應(yīng)用較多,對(duì)水華直接模擬應(yīng)用較少(對(duì)水華模擬多為通過(guò)水動(dòng)力、營(yíng)養(yǎng)鹽間接聯(lián)系,或?qū)⒃撃Ｐ椭性孱?lèi)生長(zhǎng)動(dòng)力學(xué)模塊與其他模型耦合),而對(duì)于實(shí)驗(yàn)室尺度的藻類(lèi)生長(zhǎng)模擬幾乎沒(méi)有;現(xiàn)有宏觀(guān)尺度水華模擬應(yīng)用中對(duì)于藻類(lèi)生態(tài)動(dòng)力參數(shù)極少進(jìn)行系統(tǒng)細(xì)致的研究[23-26]。由此,目前水華模擬的難點(diǎn)在于如何合理有效地選擇藻類(lèi)的生態(tài)動(dòng)力學(xué)參數(shù),以提高模擬效果。

銅綠微囊藻(Microcystisaeruginosa)水華是目前湖庫(kù)中主要的有害藍(lán)藻水華[27-28]。因此筆者選擇銅綠微囊藻作為藻生長(zhǎng)數(shù)值模擬的研究對(duì)象,旨在證明CE-QUAL-W2模型能較好地模擬室內(nèi)養(yǎng)藻過(guò)程,并研究關(guān)鍵參數(shù)對(duì)模擬效果的影響規(guī)律以及實(shí)驗(yàn)室條件下模型關(guān)鍵參數(shù)的合理優(yōu)化取值,為后續(xù)研究藻類(lèi)數(shù)值模擬提供借鑒。

1 研究方法

1.1 藻類(lèi)生長(zhǎng)動(dòng)力學(xué)模型

CE-QUAL-W2中的生態(tài)動(dòng)力學(xué)模型[21]將藻作為水體物質(zhì)的一種進(jìn)行計(jì)算,其對(duì)流擴(kuò)散方程為

(1)

式(1)中Φ為橫向平均組分濃度,g·m-3;U為x垂向流速,m·s-1;W為z垂向流速,m·s-1;B為水面寬,m;Dx為溫度和組分縱向擴(kuò)散系數(shù),m2·s-1;Dz為溫度和組分垂向擴(kuò)散系數(shù),m2·s-1;qΦ為單位體積內(nèi)物質(zhì)橫向流入或流出的量,g·m-3·s-1;SΦ為橫向平均的源匯項(xiàng),g·m-3·s-1。

該方程中的源匯項(xiàng)如圖1中藻的生態(tài)學(xué)過(guò)程所示,增量只有藻類(lèi)生長(zhǎng),損量包括藻類(lèi)的(暗)呼吸、排泄、死亡、沉降以及被捕食。

圖1 藻的生態(tài)學(xué)過(guò)程Fig.1 Schematic diagram of the algae ecological process

源匯項(xiàng)Sa的計(jì)算公式為

(2)

Kag=γarγafλminKag,max。

(3)

式(2)～(3)中,Z為網(wǎng)格厚度,m;Kag為藻類(lèi)生長(zhǎng)速率,d-1;Kar為藻類(lèi)(暗)呼吸速率,d-1;Kae為藻類(lèi)排泄速率,d-1;Kam為藻類(lèi)死亡率,d-1;ωa為藻類(lèi)沉降速率,m·d-1;Zμ為浮游動(dòng)物凈生長(zhǎng)速率,d-1;σalg、σpom、σzoo分別為浮游動(dòng)物捕食藻類(lèi)、活性顆粒態(tài)有機(jī)質(zhì)、浮游動(dòng)物的喜好因子;Φa為藻類(lèi)質(zhì)量濃度,g·m-3；Φlpom為活性顆粒態(tài)有機(jī)質(zhì)質(zhì)量濃度,g·m-3;Φz(mì)oo為浮游動(dòng)物質(zhì)量濃度,g·m-3;γar和γaf分別為藻的生長(zhǎng)曲線(xiàn)上升段和下降段的溫度系數(shù);λmin為藻類(lèi)生長(zhǎng)限制因子(光照、營(yíng)養(yǎng)鹽)的系數(shù);Kag,max為最大生長(zhǎng)速率,d-1。

在數(shù)值模擬中,藻類(lèi)最大生長(zhǎng)率代表其生長(zhǎng)“能力”,而實(shí)際生長(zhǎng)率則受光照、溫度和營(yíng)養(yǎng)鹽等條件影響。實(shí)驗(yàn)所用藻類(lèi)培養(yǎng)液以BG-11培養(yǎng)基為基礎(chǔ),N、P營(yíng)養(yǎng)鹽濃度根據(jù)實(shí)驗(yàn)需求進(jìn)行相應(yīng)調(diào)整,未涉及硅營(yíng)養(yǎng)鹽。共選取11種動(dòng)力學(xué)參數(shù)為研究對(duì)象,探究其對(duì)銅綠微囊藻生長(zhǎng)的影響(表1)。

表1銅綠微囊藻的關(guān)鍵動(dòng)力學(xué)參數(shù)及其取值依據(jù)

Table1KeykineticparametersofMicrocystisaeruginosaandbasisoftheirvalues

1.2 PIKAIA遺傳算法的自動(dòng)率定方法

由于采用Fortran 90語(yǔ)言編寫(xiě)的PIKAIA遺傳算法以適應(yīng)度函數(shù)最大化為目標(biāo)[30],故選擇常用于CE-QUAL-W2模型評(píng)價(jià)的均方根的倒數(shù)作為適應(yīng)度函數(shù)。適應(yīng)度函數(shù)計(jì)算公式[31]為

(4)

式(4)中,Oi,j為實(shí)測(cè)值;Pi,j為模擬值;m為實(shí)測(cè)與模擬匹配的數(shù)目;wi為權(quán)重因子;q為待率定參數(shù)的個(gè)數(shù)。選擇葉綠素a(Chl-a)作為研究指標(biāo)。

2 模型構(gòu)建及相關(guān)條件的處理

2.1 銅綠微囊藻室內(nèi)生長(zhǎng)實(shí)驗(yàn)

2.1.1實(shí)驗(yàn)裝置

銅綠微囊藻室內(nèi)生長(zhǎng)實(shí)驗(yàn)裝置如圖2所示,每組工況設(shè)置2個(gè)平行樣,各水框之間水體不連通,視為互不影響。單個(gè)水框的尺寸約為70 cm×40 cm×17 cm,水深10 cm,水框上方布設(shè)一排日光燈管。

2.1.2實(shí)驗(yàn)設(shè)計(jì)

室內(nèi)生長(zhǎng)實(shí)驗(yàn)控制的單一變量為水溫,水溫由可調(diào)式絕緣電熱棒控制。根據(jù)銅綠微囊藻AT1、AT2、AT3、AT4這4個(gè)特征溫度的5個(gè)區(qū)間設(shè)定5個(gè)具有代表性的水溫工況,依次為(10±1)、(14±1)、(24±1)、(28±1)、(32±1) ℃,編號(hào)為1#～5#。實(shí)驗(yàn)藻種購(gòu)自中國(guó)科學(xué)院淡水藻種庫(kù),具體培養(yǎng)條件如表2所示。光照測(cè)定使用TES-1339照度計(jì),室溫和濕度測(cè)定采用溫濕度計(jì),pH值、營(yíng)養(yǎng)鹽和Chl-a濃度測(cè)定方法參照文獻(xiàn)[32]。

2.1.3銅綠微囊藻室內(nèi)生長(zhǎng)結(jié)果

實(shí)驗(yàn)周期約為14 d,覆蓋銅綠微囊藻從接種到死亡的整個(gè)過(guò)程。取樣時(shí)間為每天08:00,各裝置的Chl-a濃度及水溫實(shí)測(cè)值如圖3所示。

靜水實(shí)驗(yàn)的水溫由空調(diào)和電熱棒共同控制,低溫水平較為穩(wěn)定,高溫水平則較難保持恒定,在控制水平上下稍有浮動(dòng)。當(dāng)水溫為10 ℃左右時(shí),銅綠微囊藻極難生長(zhǎng),其Chl-a質(zhì)量濃度不超過(guò)2 mg·m-3;當(dāng)水溫在14 ℃左右時(shí),銅綠微囊藻直到實(shí)驗(yàn)周期最后3 d才顯現(xiàn)出微弱的增長(zhǎng)趨勢(shì);水溫為24、28和32 ℃左右時(shí),Chl-a濃度變化與藻類(lèi)生長(zhǎng)理論曲線(xiàn)基本吻合,經(jīng)歷遲緩期、對(duì)數(shù)生長(zhǎng)期、穩(wěn)定期和衰亡期4個(gè)階段;其中,在28 ℃左右條件下Chl-a濃度峰值最高。這說(shuō)明5組工況中,10和14 ℃左右并不在銅綠微囊藻的適宜溫度范圍內(nèi),而28 ℃左右為銅綠微囊藻生長(zhǎng)最適溫度。因此模型使用(28±1) ℃實(shí)驗(yàn)進(jìn)行率定,另外4組進(jìn)行驗(yàn)證。

圖2 實(shí)驗(yàn)裝置示意Fig.2 Schematic diagram of the experimental apparatus

表2室內(nèi)實(shí)驗(yàn)各條件情況及其實(shí)現(xiàn)方法

Table2Conditionsofthelaboratoryexperimentandmethodsadopted

實(shí)驗(yàn)條件 控制手段或方法 水平供試藻液BG-11培養(yǎng)后離心并除去細(xì)胞表面營(yíng)養(yǎng)鹽銅綠微囊藻915(Microcystisaeruginosa)供試水體自來(lái)水代替水庫(kù)原水、蒸餾水補(bǔ)水初始放置2d,每日補(bǔ)水光照日光燈管[光]照度2000～2300lx初始pH值氫氧化鈉溶液8.6初始ρ(TN)硝酸鉀溶液4.3mg·L-1初始ρ(TP)磷酸二氫鉀溶液0.3mg·L-1室溫空調(diào)14～17℃相對(duì)濕度空調(diào)60%～70%水溫可調(diào)式絕緣電熱棒(10±1)、(14±1)、(24±1)、(28±1)、(32±1)℃

圖3 室內(nèi)實(shí)驗(yàn)水溫及葉綠素a濃度隨時(shí)間的變化Fig.3 Measured data of water temperature and Chlorophyll a in laboratory experiment

2.2 銅綠微囊藻生長(zhǎng)數(shù)學(xué)模型

2.2.1網(wǎng)格劃分

由圖4可見(jiàn),模型網(wǎng)格長(zhǎng)度Δx=10 cm,網(wǎng)格厚度Δz=5 cm,加上虛擬層,共計(jì)(2+2)×(7+2)=36個(gè)網(wǎng)格(圖中未顯示虛擬層)。

圖4 水框模型網(wǎng)格示意圖Fig.4 Grid diagram of the tank model

2.2.2模型初邊值簡(jiǎn)化

模型的邊界條件和初始條件簡(jiǎn)化設(shè)置見(jiàn)表3。

表3模型初邊值設(shè)置

Table3Initialboundaryvaluesofthemodel

初邊界條件具體取值設(shè)定依據(jù)出/入流量0實(shí)際情況水深10cm實(shí)際情況初始ρ(TN)4.3mg·L-1實(shí)際情況初始ρ(TP)0.3mg·L-1實(shí)際情況相對(duì)濕度65%實(shí)際情況取中值太陽(yáng)輻射22W·m-2實(shí)際情況換算后取中值云量0實(shí)際情況風(fēng)速、風(fēng)向0實(shí)際情況氣溫、露點(diǎn)溫度見(jiàn)表4根據(jù)水溫工況和相對(duì)濕度換算

光照以實(shí)驗(yàn)條件為準(zhǔn),[光]照度與輻射之間的換算關(guān)系[33]為

E=1.04L/100。

(5)

式(5)中,L為[光]照度,lx;E為太陽(yáng)總輻射,W·m-2。實(shí)驗(yàn)中L為2 000～2 300 lx,換算成E為20～24 W·m-2,模型中取中值,即22 W·m-2。

模型中的水溫由太陽(yáng)輻射、氣溫和露點(diǎn)溫度共同決定,實(shí)驗(yàn)中光照和相對(duì)濕度一定,水溫通過(guò)氣溫調(diào)節(jié),露點(diǎn)溫度根據(jù)式(6)換算[34]。5個(gè)水溫工況對(duì)應(yīng)的氣溫和露點(diǎn)條件見(jiàn)表4。

(6)

式(6)中,T為氣溫,℃;Tdp為露點(diǎn)溫度,℃;HR為相對(duì)濕度,%;γ(T,HR)為中間函數(shù);b、c為常數(shù),分別取值17.271、237.7 ℃。

表4模型關(guān)鍵溫度值設(shè)置

Table 4 Critical temperature of the model ℃

3 參數(shù)率定與模型驗(yàn)證

3.1 動(dòng)力學(xué)參數(shù)敏感性分析

根據(jù)CE-QUAL-W2模型中生態(tài)動(dòng)力學(xué)模型的各項(xiàng)參數(shù)模擬藻類(lèi)生長(zhǎng)過(guò)程需要率定的動(dòng)力學(xué)參數(shù)有11個(gè),各參數(shù)取值范圍見(jiàn)表5。

所有11個(gè)參數(shù)的敏感性分析遵循控制單一變量原則。每個(gè)參數(shù)在其建議范圍內(nèi)取5～6個(gè)水平進(jìn)行試算,結(jié)果見(jiàn)圖5。根據(jù)各參數(shù)對(duì)Chl-a的影響規(guī)律,分析發(fā)現(xiàn)AG、AS、 ALGP、ALGN、ASAT和ACHLA這6個(gè)關(guān)鍵參數(shù)有較大影響,為藻類(lèi)生長(zhǎng)動(dòng)力學(xué)模型的敏感參數(shù);而AR、AE、AM、AHSP和AHSN這5個(gè)參數(shù)影響較小,為不敏感參數(shù)。

6個(gè)敏感參數(shù)對(duì)Chl-a濃度具體影響為:AG越大,Chl-a峰值越大(正相關(guān)),峰現(xiàn)時(shí)間越早(負(fù)相關(guān));AS越大,Chl-a峰值越小(負(fù)相關(guān)),峰現(xiàn)時(shí)間越晚(正相關(guān));ALGP越大,Chl-a峰值越小(負(fù)相關(guān)),峰現(xiàn)時(shí)間越早(負(fù)相關(guān));ALGN規(guī)律與ALGP相同:在10～50 W·m-2范圍內(nèi),ASAT越大,Chl-a峰值越大(正相關(guān)),峰現(xiàn)時(shí)間越早(負(fù)相關(guān)),在50～170 W·m-2范圍內(nèi),ASAT規(guī)律相反;ACHLA越大,Chl-a峰值越小(負(fù)相關(guān)),峰現(xiàn)時(shí)間不變(無(wú)關(guān))。

表5藻類(lèi)生長(zhǎng)動(dòng)力模型參數(shù)取值范圍[21]

Table5Valuerangesofthekineticparametersofthealgaegrowth

參數(shù)單位取值范圍參數(shù)單位取值范圍AGd-10.5～8.0ACHLAmg·μg-10.01～0.40ARd-10.01～0.25AHSPg·m-30.0005～0.0800AEd-10.01～0.25AHSNg·m-30～0.4AMd-10.001～0.200ALGP0.0008～0.0300ASm·d-1-0.1～10.0ALGN0.006～0.160ASATW·m-210～170

AG為最大生長(zhǎng)速率; AR為(暗)呼吸速率; AE為排泄速率; AM為死亡率; AS為沉降速率; ASAT為最適光強(qiáng)(光合速率最大時(shí)的光強(qiáng)); ACHLA為單位藻干重葉綠素a含量; AHSP為吸收磷的半飽和常數(shù); AHSN為吸收氮的半飽和常數(shù); ALGP為藻細(xì)胞磷含量; ALGN為藻細(xì)胞氮含量。

AG—最大生長(zhǎng)速率； AS—沉降速率； ALGP—藻細(xì)胞磷含量； ALGN—藻細(xì)胞氮含量； ASAT—最適光強(qiáng)； ACHLA—單位藻干重葉綠素a含量;AHSP—吸收磷的半飽和常數(shù); AHSN—吸收氮的半飽和常數(shù); AR—(暗)呼吸速率; AE—排泄速率; AM—死亡率。

3.2 敏感參數(shù)率定

模型以28 ℃水溫條件下的實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ),經(jīng)過(guò)反復(fù)自動(dòng)率定得到局部最優(yōu)解后再進(jìn)行手動(dòng)微調(diào),最終率定結(jié)果見(jiàn)圖6,Chl-a濃度峰出現(xiàn)時(shí)間與實(shí)驗(yàn)保持一致(第7天)，峰值相對(duì)誤差為0.3%,模型關(guān)鍵參數(shù)取值:AG為7.526 8 d-1,AS為 0.002 2 m·d-1,ASAT為102.774 4 W·m-2,ALGP為0.000 5,ALGN為0.041 3,ACHLA為0.125 8 mg·μg-1。

3.3 模型驗(yàn)證

基于率定所得參數(shù),通過(guò)改變模型氣象條件從而改變水溫條件,得到另外4組水溫條件下的Chl-a模擬結(jié)果,與實(shí)測(cè)數(shù)據(jù)對(duì)比情況見(jiàn)圖7。模型關(guān)于Chl-a量級(jí)的模擬準(zhǔn)確度高,10和14 ℃左右水溫對(duì)銅綠微囊藻的生長(zhǎng)都極為不利,因而Chl-a質(zhì)量濃度極低,<10 μg·L-1。在14 ℃左右水溫條件下,實(shí)驗(yàn)最后2 d銅綠微囊藻才呈現(xiàn)出微弱的增長(zhǎng)趨勢(shì),模型模擬結(jié)果與之相符。在24和32 ℃左右水溫條件下,模擬的峰現(xiàn)時(shí)間與實(shí)測(cè)一致(第7天),峰值的相對(duì)誤差分別為18%和23%,說(shuō)明模型能準(zhǔn)確模擬出峰值和峰現(xiàn)時(shí)間,基本能擬合銅綠微囊藻的緩慢生長(zhǎng)期和衰亡過(guò)程。

圖6 28 ℃水溫條件下模型率定結(jié)果Fig.6 Calibration of the model at 28 ℃(WT)

°實(shí)測(cè)值； —模擬值。當(dāng)水溫為10 ℃時(shí)，實(shí)測(cè)和模擬的ρ(Chl-a)幾乎為0。

4 結(jié)果與討論

4.1 銅綠微囊藻生長(zhǎng)動(dòng)力學(xué)參數(shù)敏感性分析及其取值優(yōu)化

CE-QUAL-W2模型對(duì)于實(shí)驗(yàn)室不同水溫條件下銅綠微囊藻生長(zhǎng)的模擬應(yīng)用效果可觀(guān),敏感性生長(zhǎng)動(dòng)力參數(shù)包括AG、AS、ASAT、ALGN、ALGP和ACHLA。已有的水華模擬研究中銅綠微囊藻的AG取值一般<2.0 d-1[29,35-36],但由圖5可知，當(dāng)AG<2.0 d-1時(shí),銅綠微囊藻的增長(zhǎng)趨勢(shì)很弱,增長(zhǎng)速度極低,除非藻本身的活性低,否則AG范圍建議為2.0～8.0 d-1。AS≥0.5 m·d-1時(shí),銅綠微囊藻幾乎無(wú)增長(zhǎng)可能;AS為-0.1和0 m·d-1時(shí),Chl-a濃度曲線(xiàn)重合,此時(shí)的生態(tài)學(xué)意義是銅綠微囊藻的偽空泡調(diào)解其細(xì)胞密度從而垂直向上遷移或懸浮于水中,但室內(nèi)實(shí)驗(yàn)過(guò)程中銅綠微囊藻沒(méi)有出現(xiàn)過(guò)上浮,因此在銅綠微囊藻室內(nèi)條件下AS取值范圍建議為0～0.5 m·d-1。由圖5可知,ASAT對(duì)Chl-a的影響大致分為2個(gè)區(qū)間,分別為10～50和50～170 W·m-2,2個(gè)區(qū)間的響應(yīng)規(guī)律相反,具體變化的閾值還需進(jìn)一步試算。ALGN、ALGP、ACHLA這3個(gè)參數(shù)均可參考表5中的取值范圍。

以水溫為28 ℃的模型為基礎(chǔ),經(jīng)過(guò)多次率定,以峰值和峰現(xiàn)時(shí)間為優(yōu)先擬合準(zhǔn)則,得到較優(yōu)的參數(shù)組合如下:AG為7.526 8 d-1,AS為0.002 2 m·d-1,ASAT為102.774 4 W·m-2,ALGP為0.000 5,ALGN為0.041 3,ACHLA為0.125 8 mg·μg-1。該組參數(shù)在溫度約為10、14、24和32 ℃的模型中得到良好驗(yàn)證。

4.2 銅綠微囊藻數(shù)值模擬的合理性分析

模擬結(jié)果與實(shí)測(cè)值的差異可能有以下2個(gè)原因:首先是水溫條件,在模型中每種工況的水溫恒定不變,而室內(nèi)實(shí)驗(yàn)中水溫會(huì)在一定范圍內(nèi)波動(dòng);其次是選取的對(duì)比指標(biāo)Chl-a,模型所得Chl-a濃度是由藻的生物量通過(guò)ACHLA的換算后得到,而實(shí)驗(yàn)測(cè)定的Chl-a濃度并不完全與藻的生物量對(duì)應(yīng),還包括部分未被分解的死亡藻細(xì)胞中的Chl-a。因此,為提高模擬精度,建議后續(xù)研究在測(cè)定Chl-a濃度的同時(shí),還需通過(guò)顯微鏡測(cè)量藻的生物量。

該銅綠微囊藻室內(nèi)生長(zhǎng)實(shí)驗(yàn)與野外最大的不同在于,室內(nèi)實(shí)驗(yàn)中的微囊藻是以單個(gè)細(xì)胞存在,不會(huì)在水面聚集成團(tuán),雖然理論上可以上下浮動(dòng),但由于培養(yǎng)環(huán)境相對(duì)穩(wěn)定,并未出現(xiàn)上浮現(xiàn)象;而野外銅綠微囊藻大量聚集成團(tuán),并且分泌胞外多糖和調(diào)解偽空泡大小適應(yīng)自然環(huán)境,具有顯著的垂向遷移規(guī)律。李小龍等[37]研究表明,銅綠微囊藻光飽和點(diǎn)在500 μmol·m-2·s-1(約260 W·m-2)附近,而光強(qiáng)達(dá)到900 μmol·m-2·s-1(約468 W·m-2)時(shí)仍無(wú)光抑制現(xiàn)象發(fā)生,因此ASAT的最優(yōu)取值可能會(huì)大于建議的取值范圍上限170 W·m-2。而AS可能取負(fù)值,其他具有主動(dòng)遷移能力的藻類(lèi)(如擬多甲藻等)AS也可能為負(fù)值,但是無(wú)主動(dòng)遷移能力的藻類(lèi)(如隱藻、硅藻等)AS取值只能為正。后續(xù)研究建議在可控的實(shí)驗(yàn)條件下,有目的地探索各類(lèi)生長(zhǎng)動(dòng)力參數(shù)與藻類(lèi)生長(zhǎng)變化之間的敏感性關(guān)系及其響應(yīng)規(guī)律,可節(jié)省大量人力、物力和時(shí)間,同時(shí)為野外水華預(yù)測(cè)提供相關(guān)經(jīng)驗(yàn)。

[1] REYNOLDS C S.The Ecology of Phytoplankton[M].New York,USA:Cambridge University Press,2006:179-236.

[2] J?RGENSEN S E,BENDORICCHIO G.Fundamentals of Ecological Modelling[M].3rd ed.Oxford，UK:Elsevier Science Ltd.,2001:236-259.

[3] 秦伯強(qiáng),楊桂軍,馬健榮,等.太湖藍(lán)藻水華“暴發(fā)”的動(dòng)態(tài)特征及其機(jī)制[J].科學(xué)通報(bào),2016,61 (7):759-770.[QIN Bo-qiang,YANG Gui-jun,MA Jian-rong,etal.Dynamics of Variability and Mechanism of Harmful Cyanobacteria Bloom in Lake Taihu,China[J].Chinese Science Bulletin,2016,61(7):759-770.]

[4] DUR G,HWANG J S,SOUISSI S,etal.An Overview of the Influence of Hydrodynamics on the Spatial and Temporal Patterns of Calanoid Copepod Communities Around Taiwan[J].Journal of Plankton Research,2007,29(Suppl. 1):97-116.

[7] DUPUIS A P,HANN B J.Warm Spring and Summer Water Temperatures in Small Eutrophic Lakes of the Canadian Prairies:Potential Implications for Phytoplankton and Zooplankton[J].Journal of Plankton Research,2009,31(5):489-502.

[8] RALSTON D K,KEAFER B A,BROSNAHAN M L,etal.Temperature Dependence of an Estuarine Harmful Algal Bloom:Resolving Interannual Variability in Bloom Dynamics Using a Degree Day Approach[J].Limnology and Oceanography,2014,59(4):1112-1126.

[9] 袁琳娜,楊常亮,李曉銘,等.高原深水湖泊水溫日成層對(duì)溶解氧、酸堿度、總磷濃度和藻類(lèi)密度的影響:以云南陽(yáng)宗海為例[J].湖泊科學(xué),2014,26(1):161-168.[YUAN Lin-na,YANG Chang-liang,LI Xiao-ming,etal.Effect of Daily Thermal Stratification on Dissolved Oxygen,pH,Total Phosphorus Concentration,Phytoplankton and Algae Density of a Deep Plateau Lake:A Case Study of Lake Yangzonghai,Yunnan Province[J].Journal of Lake Sciences,2014,26(1):161-168.]

[10] 孔繁翔,高光.大型淺水富營(yíng)養(yǎng)化湖泊中藍(lán)藻水華形成機(jī)理的思考[J].生態(tài)學(xué)報(bào),2005,25(3):589-595.[KONG Fan-xiang,GAO Guang.Hypothesis on Cyanobacteria Bloom-Forming Mechanism in Large Shallow Eutrophic Lakes[J].Acta Eologica Snica,2005,25(3):589-595.]

[11] 徐寧,陳菊芳.廣東大亞灣藻類(lèi)水華的動(dòng)力學(xué)分析:Ⅰ.藻類(lèi)水華的生消過(guò)程及其與環(huán)境因子的關(guān)系[J].海洋環(huán)境科學(xué),2001,20(2):1-6.[XU Ning,CHEN Ju-fang.Dynamic Analysises on Several Algal Bloom Events in Daya Bay of Guangdong:Ⅰ.Process of Algal Bloom and Its Relationship With Environmental Factors[J].Marine Environmental Science,2001,20(2):1-6.]

[12] CARNEIRO R L,SANTOS M E V D,PACHECO A B F,etal.Effects of Light Intensity and Light Quality on Growth and Circadian Rhythm of Saxitoxins Production inCylindrospermopsisraciborskii(Cyanobacteria)[J].Journal of Plankton Research,2009,31(5):481-488.

[13] EDWARDS K F,THOMAS M K,KLAUSMEIER C A,etal.Phytoplankton Growth and the Interaction of Light and Temperature:A Synthesis at the Species and Community Level[J].Limnology and Oceanography,2016,61(4):1232-1244.

[15] JAKOBSEN H H,MARKAGER S.Carbon-to-Chlorophyll Ratio for Phytoplankton in Temperate Coastal Waters:Seasonal Patterns and Relationship to Nutrients[J].Limnology and Oceanography,2016,61(5):1853-1868.

[16] HARDING L W,GALLEGOS C L,PERRY E S,etal.Long-Term Trends of Nutrients and Phytoplankton in Chesapeake Bay[J].Estuaries and Coasts,2016,39(3):664-681.

[17] 李哲,方芳,郭勁松,等.三峽小江回水區(qū)段 2007 年春季水華與營(yíng)養(yǎng)鹽特征[J].湖泊科學(xué),2009,21(1):36-44.[LI Zhe,FANG Fang,GUO Jin-song,etal.Spring Algal Bloom and Nutrients Characteristics in Xiaojiang River Backwater Area,Three Gorge Reservoir,2007[J].Journal of Lake Sciences,2009,21(1):36-44.]

[18] 鄭丙輝,曹承進(jìn),秦延文,等.三峽水庫(kù)主要入庫(kù)河流氮營(yíng)養(yǎng)鹽特征及其來(lái)源分析[J].環(huán)境科學(xué),2008,29(1):1-6.[ZHENG Bing-hui,CAO Cheng-jin,QIN Yan-wen,etal.Analysis of Nitrogen Distribution Characters and Their Sources of the Major Input Rivers of Three Gorges Reservoir[J].Environmental Science,2008,29(1):1-6.]

[19] 孔繁翔,馬榮華,高俊峰,等.太湖藍(lán)藻水華的預(yù)防、預(yù)測(cè)和預(yù)警的理論與實(shí)踐[J].湖泊科學(xué),2009,21(3):314-328.[KONG Fan-xiang,MA Rong-hua,GAO Jun-feng,etal.The Theory and Practice of Prevention,Forecast and Warning on Cyanobacteria Bloom in Lake Taihu[J].Journal of Lake Sciences,2009,21(3):314-328.]

[20] WHITEHEAD P G,BUSSI G,BOWES M J,etal.Dynamic Modelling of Multiple Phytoplankton Groups in Rivers With an Application to the Thames River System in the UK[J].Environmental Modelling & Software,2015,74:75-91.

[21] THOMAS M C,EDWARD M B.CE-QUAL-W2:A Two-Dimensional,Laterally Averaged,Hydrodynamic and Water Quality Model,Version 2.0[Z].Portland:Portland State University,2016.

[22] SVERDRUP H U.On Conditions for the Vernal Blooming of Phytoplankton[J].Ices Journal of Marine Science,1953,18(3):287-295.

[23] 謝濤,紀(jì)道斌,尹衛(wèi)平,等.三峽水庫(kù)不同下泄流量香溪河水動(dòng)力特性與水華的響應(yīng)[J].中國(guó)農(nóng)村水利水電,2013(11):1-6.[XIE Tao,JI Dao-bin,YI Wei-ping,etal.The Hydroynamic of Xiangxi Bay Response to Water Bloom on Different Discharges of the Three Gorges Reservoir[J].China Rural Water and Hydropower,2013(11):1-6.]

[24] 姚燁.基于抑制近壩支流水華的三峽水庫(kù)非汛期優(yōu)化調(diào)度研究[D].天津:天津大學(xué),2014.[YAO Ye.Research on Three Gorges Reservoir Operation for Non-Flood Season Based on Controlling Algal Blooms in Near-Dam Tributaries[D].Taijin:Tianjin University,2014.]

[25] 龍良紅,紀(jì)道斌,劉德富,等.基于CE-QUAL-W2模型的三峽水庫(kù)神農(nóng)溪庫(kù)灣水流水溫特性分析[J].長(zhǎng)江科學(xué)院院報(bào),2016,33(5):28-35.[LONG Liang-hong,JI Dao-bin,LIU De-fu,etal.Characteristics of Flow Velocity and Water Temperature in Shennong Bay Using CE-QUAL-W2 Model[J].Journal of Yangtze River Scientific Research Institute,2016,33(5):28-35.]

[26] 李光浩,楊霞,陳和春,等.三峽水庫(kù)香溪河庫(kù)灣水華水動(dòng)力調(diào)控初步研究[J].人民長(zhǎng)江,2017,48(10):18-23,29.[LI Guang-hao,YANG Xia,CHEN He-chun,etal.Preliminary Study on Hydrodynamic Regulation of Algal Bloom in Xiangxi River Bay[J].Yangtze River,2017,48(10):18-23,29.]

[27] 王蕾,石利利,蔡道基.藍(lán)藻水華優(yōu)勢(shì)藻高效防控銅制劑的篩選[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2011,27(1):60-64.[WANG Lei,SHI Li-li,CAI Dao-ji.Screening of Copper Algaecide to Suppress Dominant Algae of Cyanobacteria Algal Bloom[J].Journal of Ecology and Rural Environment,2011,27(1):60-64.]

[28] 方麗娟,劉德富,楊正健,等.三峽水庫(kù)香溪河庫(kù)灣夏季浮游植物演替規(guī)律及其原因[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2013,29(2):234-240.[FANG Li-juan,LIU De-fu,YANG Zheng-jian,etal.Succession of Phytoplankton in Xiangxi Bay of Three-Gorge Reservoir in Summer and Its Causes[J].Journal of Ecology and Rural Environment,2013,29(2):234-240.]

[29] 崔玉潔.三峽水庫(kù)香溪河藻類(lèi)生長(zhǎng)敏感生態(tài)動(dòng)力學(xué)過(guò)程及其模擬[D].武漢:武漢大學(xué),2017.[CUI Yu-jie.The Sensitive Ecological Dynamic Processes and Their Simulations of Algal Growth of Xiangxi Bay in the Three Gorges Reservoir[D].Wuhan:Wuhan University,2017.]

[30] CHARBONNEAU P.Genetic Algorithms in Astronomy and Astrophysics[J].The Astrophysical Journal Supplement Series,1995,101:309-334.

[31] PELLETIER G J,CHAPRA S C,TAO H.QUAL2Kw:A Framework for Modeling Water Quality in Streams and Rivers Using a Genetic Algorithm for Calibration[J].Environmental Modelling & Software,2006,21(3):419-425.

[32] 國(guó)家環(huán)境保護(hù)總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì).水和廢水監(jiān)測(cè)分析方法[M].4版.北京：中國(guó)環(huán)境科學(xué)出版社，2002：88-284.[Editorial Board of “Water and Wastewater Monitoring and Analysis Methods” of State Enviromental Protection Administration.Water and Wastewater Monitoring and Analysis Methods[M].4th ed.Beijing:China Environmental Science Press,2002:88-284.]

[33] 毛軍需.太陽(yáng)輻射度與光照度換算淺析[J].河南農(nóng)業(yè)科學(xué),1995(1):11-12.[MAO Jun-xu.Analysis of Solar Radiation Degree and Illuminance Conversion[J].Journal of Henan Agricultural Sciences,1995(1):11-12.]

[34] BOLTON D.The Computation of Equivalent Potential Temperature[J].Monthly Weather Review,1980,108:1046-1053.

[35] 馬駿.三峽水庫(kù)香溪河庫(kù)灣水動(dòng)力循環(huán)模式及其環(huán)境效應(yīng)研究[D].南京:河海大學(xué),2015.[MA Jun.Study on Hydrodynamic Circulation Patterns and Their Environment Effects in Xiangxi Bay of the Three Gorges Reservoir[D].Nanjing:Hohai University,2015.]

[36] 嚴(yán)萌.不同環(huán)流模式對(duì)香溪河庫(kù)灣水華生消影響規(guī)律研究[D].宜昌:三峽大學(xué),2016.[YAN Meng.Reserch on the Influence of Different Hydrodynamic Circulation Patterns Impacting on Formation and Disappearance of Algal Bloom in Xiangxi Bay[D].Yichang:China Three Gorges University,2016.]

[37] 李小龍,耿亞紅,李夜光,等.從光合作用特性看銅綠微囊藻(Microcystisaeruginosa)的競(jìng)爭(zhēng)優(yōu)勢(shì)[J].武漢植物學(xué)研究,2006,24(3):225-230.[LI Xiao-long,GENG Ya-hong,LI Ye-guang,etal.The Advantages in Competition Based on the Photosynthetic Characteristics ofMicrocystisaeruginosa[J].Joumal of Wuhan Botanical Research,2006,24(3):225-230.]