鄱陽湖區(qū)藻類豐度分布及其與環(huán)境因子的關(guān)聯(lián)分析

程新+黃林+李昆太

摘要：于2012年10月至2013年7月在江西省鄱陽湖4個采樣點進行了采樣，并對表層水體中藻類豐度與環(huán)境因子之間的關(guān)系進行研究。結(jié)果表明，藻類豐度隨季節(jié)變化差異顯著（P<0.05），夏季最高[（2.58±1.85）×106 CFU/L]，冬季最低[（4.3±3.2）×105 CFU/L]，水體理化因子也隨季節(jié)和采樣點的不同而有很大的變化。通過相關(guān)分析和GAM分析表明，鄱陽湖藻類豐度與水溫和葉綠素a呈極顯著的正相關(guān)（P<0.01），與氨氮、總磷呈極顯著負相關(guān)（P<0.01）。從總體情況看，鄱陽湖水體藻類豐度呈上升趨勢，水質(zhì)正在逐步惡化，需要加以關(guān)注。

關(guān)鍵詞：鄱陽湖；藻類豐度；廣義可加模型；環(huán)境因子

中圖分類號：Q949.1；X173 文獻標識碼：A 文章編號：0439-8114（2017）15-2865-05

DOI：10.14088/j.cnki.issn0439-8114.2017.15.017

Abstract： Surface water samples from 4 sites of Poyang Lake in Jiangxi Province were collected from October 2012 to July 2013. Environmental variables of each sampling site were measured using standard methods，and the relationship between phytoplankton abundance and environmental parameters was explored. The results showed that the phytoplankton abundance varied significantly with season（P<0.05），with the highest value being recorded in summer [（2.58±1.85）×106 CFU/L] and the lowest in winter[（4.3±3.2）×105 CFU/L]. The results of correlation analysis and GAM showed that phytoplankton abundance was positive correlation with temperature and chlorophyll a（Chl-a） significantly （P<0.01），and was negative correlation with total phosphorus and amino nitrogen significantly（P<0.01）. Overall，phytoplankton abundance in Poyang Lake was on the rise，and water quality was gradually worsening in recent years，which should be paid attention.

Key words： Poyang Lake； phytoplankton abundance； generalized additive models； environmental parameters

鄱陽湖（115°49′-116°46′N，28°24′-29°46′E）位于江西省北部，長江中下游交接處南岸，是中國最大的淡水湖泊和十大生態(tài)功能保護區(qū)之一，也是具有國際性保護意義的淡水濕地和世界自然基金會劃定的全球重要生態(tài)區(qū)之一，對維護長江下游地區(qū)的水生態(tài)安全具有重要作用[1]。

浮游藻類作為湖泊的初級生產(chǎn)者，在維持水生生態(tài)系統(tǒng)的平衡方面起著重要的作用， 而浮游藻類對水體的污染有著明顯的指示作用[2，3]。近年來隨著工業(yè)化、城鎮(zhèn)化進程的推進，大量污廢水進入鄱陽湖，致使其水質(zhì)每況愈下。目前，鄱陽湖區(qū)水體氮、磷含量已接近富營養(yǎng)化水平，且有進一步惡化的趨勢[4]，鄱陽湖水體已經(jīng)具備了藍藻爆發(fā)的營養(yǎng)條件。

鄱陽湖水質(zhì)問題已經(jīng)引起人們的廣泛關(guān)注，但大多數(shù)研究集中在對氮、磷污染或其他營養(yǎng)因子的分析或?qū)υ孱惷芏鹊恼{(diào)查[4，5]，而對于湖區(qū)不同季節(jié)的藻類豐度及其與環(huán)境因子關(guān)系的研究尚未見報道。本研究通過對鄱陽湖4個采樣點浮游藻類的豐度及水體理化因子的測定，利用廣義可加模型（Generalize addictive model，GAM）闡述鄱陽湖藻類豐度與環(huán)境理化因子的響應(yīng)關(guān)系，為鄱陽湖的環(huán)境保護和生態(tài)修復提供數(shù)據(jù)支持。

1 材料與方法

1.1 樣品采集

根據(jù)鄱陽湖水面分布情況，共設(shè)置4個采樣點，分別為蚌湖、星子、棠蔭和湖口。藻類樣品采集方法參照文獻[5]，采樣時間為一個水文周期，分別為2012年10月和2013年1月、4月和7月?，F(xiàn)場采集5 L水樣，理化測定用水樣經(jīng)0.45 μm的微孔濾膜過濾后，加H2SO4并置于冰上保存，3 h內(nèi)運回實驗室。采樣期間鄱陽湖水體溫度變化見圖1。

1.2 藻類物種鑒定、計數(shù)和水質(zhì)理化參數(shù)測定

浮游藻類種類鑒定參照文獻[6]；浮游藻類數(shù)量的測定參照文獻[7]，采用顯微鏡視野法計數(shù)。

水質(zhì)理化參數(shù)的測定包括pH、溶解氧（DO）、溫度（T）、透明度（SD）、五日生化需氧量（BOD5）、化學需氧量（COD）、總磷（TP）、總氮（TN）和氨氮（NH4+-N），測定方法參照《水和廢水監(jiān)測分析方法》[8]。多樣性指數(shù)計算參照文獻[2]。

1.3 數(shù)據(jù)處理和模型建立

數(shù)據(jù)分析采用軟件R 2.15.1完成，浮游藻類豐度與環(huán)境因子的相關(guān)性分析采用Pearson相關(guān)系數(shù)法。選用廣義可加模型（GAM）探討藻類豐度與環(huán)境因子的相互關(guān)系，GAM模型的建立由統(tǒng)計軟件R中的mgcv包（Version 3.0.0）完成。endprint

2 結(jié)果與分析

2.1 鄱陽湖水體主要理化特性

從圖2可以看出，不同采樣地點采集的水樣在理化性質(zhì)上存在一定的差異，但大體趨勢相同。總氮、總磷、氨氮濃度在冬季較高，而在夏季下降；葉綠素a和BOD5在夏季較高而冬季較低；COD數(shù)值的變化趨勢隨采樣點的不同而呈現(xiàn)出不同的變化趨勢。

2.2 鄱陽湖藻類豐度的季節(jié)變化

鏡檢結(jié)果表明，鄱陽湖水域藻類的種群主要由硅藻門、綠藻門和藍藻門等構(gòu)成。鄱陽湖表層水體藻類豐度隨季節(jié)變化差異顯著（P<0.05），冬季藻類豐度最低，平均為（4.3±3.2）×105 CFU/L；夏季最高，平均為（2.58±1.85）×106 CFU/L。在不同采樣點之間，藻類豐度也存在顯著差異（P<0.05），其中夏季星子取樣點的藻類豐度最高，達到（5.6±0.23）×106 CFU/L，與此同時，棠蔭和蚌湖取樣點的藻類豐度只有（1.1±0.24）×106 CFU/L，差距較為明顯（圖3）。

根據(jù)鄱陽湖水文局在2009年7月所進行的調(diào)查顯示，鄱陽湖區(qū)棠蔭、蚌湖、湖口3個監(jiān)測點的藻類豐度均低于1.0×105 CFU/L，星子監(jiān)測點的藻類豐度為1.08×105 CFU/L，總體藻類豐度較低。本研究的結(jié)果顯示，鄱陽湖區(qū)夏季4個監(jiān)測點的藻類豐度在4年間已發(fā)生了較大變化，其數(shù)值均已大幅度提高，湖區(qū)的藻類防治工作不容忽視。

通過對湖體中藻類的鑒定結(jié)果可以看出，藍藻（Cyanobacteria）、硅藻（Bacillariophyta）和綠藻（Chlorophyta）是鄱陽湖水體中的優(yōu)勢藻，其中藍藻與硅藻的含量較高，并分別在不同季節(jié)占有優(yōu)勢，且二者占浮游植物總豐度的比例呈相反的變化趨勢（圖4）。溫度較低的季節(jié)，優(yōu)勢藻類主要是硅藻，如根管藻屬（Rhizosolenia）、針桿藻屬（Synedra）、平板藻屬（Tabellaria）、羽紋藻屬（Pinnularia）、直鏈藻屬（Melosira）等；而在溫度較高的季節(jié)，湖體中的優(yōu)勢藻類主要為藍藻，以微囊藻屬（Microcystis）為主。

分別采用Shannon多樣性指數(shù)、Margalef指數(shù)對水體中藻類多樣性進行評價，結(jié)果如圖5所示。從圖5可以看出，鄱陽湖區(qū)全年的Shannon多樣性指數(shù)變化范圍為1.17～2.69，星子、蚌湖、棠蔭和湖口4個采樣點的平均Shannon多樣性指數(shù)分別為1.62、1.93、1.56和1.61，最高和最低值分別出現(xiàn)在4月的蚌湖采樣點和1月的棠蔭采樣點。Margalef指數(shù)變化范圍為0.38～1.59，星子、蚌湖、棠蔭和湖口4個采樣點的平均Margalef指數(shù)分別為0.96、0.78、0.57和0.73，最低和最高值分別出現(xiàn)在10月的蚌湖采樣點和4月的星子采樣點?？傮w而言，春夏季節(jié)鄱陽湖區(qū)的水體生物多樣性較好，而氣溫降低后生物多樣性較差。從多樣性指數(shù)來看，鄱陽湖區(qū)水體屬于中度污染，在部分時間段接近重度污染。

2.3 鄱陽湖藻類豐度與理化因子的相關(guān)分析及GAM模型

藻類的生長與水體環(huán)境理化因子具有密切的關(guān)系，特別是氮、磷等營養(yǎng)物質(zhì)的含量起著關(guān)鍵作用。為了研究藻類豐度與環(huán)境理化因子的關(guān)系，從而確定影響藻類豐度的關(guān)鍵因子，采用Pearson相關(guān)系數(shù)法對所得結(jié)果進行分析。結(jié)果（表1）顯示，鄱陽湖藻類豐度與水溫和葉綠素a呈極顯著的正相關(guān)（P<0.01），而與氨氮、總磷呈極顯著負相關(guān)（P<0.01）。

為了進一步分析藻類豐度與各環(huán)境因子之間的關(guān)系，本研究選用廣義可加模型（GAM）探討藻類豐度與環(huán)境因子的相互關(guān)系。以藻類豐度為因變量，所監(jiān)測環(huán)境因子作為解釋變量建立GAM模型，GAM模型中各自變量之間的關(guān)系彼此獨立，可以得出各個因子與藻類豐度之間的關(guān)系。

水體中藻類的生長受很多因素的影響。楊麗標等[9]研究結(jié)果表明，總磷是影響巢湖藻類物種分布格局的重要因素。沈會濤等[10]對白洋淀藻類群落與環(huán)境因子的分析結(jié)果表明，pH和總磷是影響該區(qū)域浮游藻類分布的重要因素。從本研究GAM模型結(jié)果（圖6）可以看出，水溫對鄱陽湖藻類豐度的影響最大，藻類豐度隨著水溫的升高而增大，并且在水溫高于20 ℃時，增長速率迅速增加，曲線開始變得陡峭；葉綠素a、總磷和氨氮與藻類豐度呈良好的線性關(guān)系。

3 討論

3.1 鄱陽湖水質(zhì)有逐步下降的趨勢

從本研究的試驗結(jié)果可以看出，鄱陽湖氮、磷營養(yǎng)鹽濃度整體水平較高，營養(yǎng)鹽都是在夏季出現(xiàn)最低值，冬季出現(xiàn)最高值，說明枯水期水體營養(yǎng)程度更高。與2009年江西省鄱陽湖水文局的調(diào)查及其他相關(guān)文獻相比，鄱陽湖區(qū)水體的富營養(yǎng)化和藻類豐度均有一定程度的提高，這說明隨著環(huán)鄱陽湖區(qū)工農(nóng)業(yè)的發(fā)展和人口密度的提升，鄱陽湖區(qū)的水體生態(tài)環(huán)境有惡化的趨勢。

根據(jù)湖泊營養(yǎng)類型評價的藻類生物學標準[3]，鄱陽湖區(qū)的藻類情況總體較好，在藻類豐度最高的夏季，水體基本處于中營養(yǎng)化水平；而根據(jù)葉綠素水平評價，湖體水質(zhì)也在中營養(yǎng)化水平以下。但值得注意的是，湖體中藻類的區(qū)系可以反映水體的營養(yǎng)狀況，中營養(yǎng)型湖泊中常以甲藻、隱藻、硅藻類占優(yōu)勢，而富營養(yǎng)型湖泊則常以綠藻、藍藻類占優(yōu)勢[11，12]。根據(jù)此次的調(diào)查及計算結(jié)果，鄱陽湖水體中藍藻和綠藻的比重之和已經(jīng)超過50%，說明鄱陽湖的水體營養(yǎng)化程度正在逐步加深。李榮昉等[13]的研究結(jié)果也表明鄱陽湖水體的質(zhì)量有逐年下降的趨勢。考慮到鄱陽湖在長江中下游地區(qū)的重要生態(tài)和經(jīng)濟地位，水體污染問題亟需得到重視。

3.2 藻類豐度與環(huán)境因子之間的關(guān)聯(lián)

水體中藻類豐度與環(huán)境因子之間的關(guān)聯(lián)分析一直是研究的熱點問題，目前經(jīng)常采用的有典范對應(yīng)分析（CCA）[14]、多元回歸分析[15]等。由于藻類豐度和環(huán)境的關(guān)系錯綜復雜，單純的線性分析有時并不足以解釋二者之間的相互關(guān)系。所以本研究選用廣義可加模型（GAM）探討藻類豐度與環(huán)境因子的相互關(guān)系。GAM為廣義線性模型的非參數(shù)化擴展，其優(yōu)點是無需預定的參數(shù)模型，能直接處理響應(yīng)變量與多個解釋變量之間的非線性關(guān)系，對變量的數(shù)據(jù)類型和統(tǒng)計分布特征適應(yīng)性更廣，是目前較為先進的研究物種與環(huán)境之間關(guān)系的方法[16]，但在微生物學研究領(lǐng)域應(yīng)用較為少見。endprint

溫度可以通過影響生物體新陳代謝相關(guān)酶活性而改變環(huán)境中微生物的區(qū)系和密度。本研究結(jié)果表明，鄱陽湖水體中藻類豐度受水溫的影響較大。隨著水溫的升高，藻類豐度總體上也呈現(xiàn)上升的趨勢，因此水溫是調(diào)節(jié)鄱陽湖水體中藻類豐度季節(jié)性變化的最主要因素。

營養(yǎng)鹽也是影響細菌豐度的主要環(huán)境因子。研究表明，細菌會以提高豐度對富營養(yǎng)化過程做出響應(yīng)。在超寡營養(yǎng)湖泊至超富營養(yǎng)湖泊中，細菌的豐度一般在105～107 CFU/mL。根據(jù)本研究測定結(jié)果，盡管從總體情況來看，鄱陽湖水體還基本處于中度營養(yǎng)化階段，但近年來鄱陽湖水體中的營養(yǎng)物質(zhì)如總磷、氨氮等，其濃度都有逐步上升的趨勢，因此水體中的藻類豐度已經(jīng)超過一般的中營養(yǎng)化湖泊，如博斯騰湖[17]等，而接近一些富營養(yǎng)化湖泊如太湖、烏梁素海[18]等。

葉綠素a也是影響鄱陽湖藻類豐度的重要環(huán)境因子。葉綠素a含量代表了浮游植物的生長密度，而浮游植物與藻類的生長繁殖有著密切關(guān)系[19]。浮游植物光合產(chǎn)物——溶解性有機物（DOM）可以為藻類生長提供所需的碳源，進而促進藻類的生長繁殖，而藻類也能通過自身代謝為浮游植物的生長提供必要的有機營養(yǎng)和生長因子，并調(diào)節(jié)浮游植物生長的微環(huán)境，因此浮游植物也可以作為水體中藻類豐度的指標之一。

參考文獻：

[1] 劉倩純，余 潮，張 杰，等.鄱陽湖水體水質(zhì)變化特征分析[J].農(nóng)業(yè)環(huán)境科學學報，2013，32（6）：1232-1237.

[2] 李守淳，柴文波，葉付粥，等.南昌市艾溪湖浮游藻類的多樣性調(diào)查與評價[J].江西師范大學學報（自然科學版），2013，37（3）： 316-318.

[3] 況琪軍，馬沛明，胡征宇，等.湖泊富營養(yǎng)化的藻類生物學評價與治理研究進展[J].安全與環(huán)境學報，2005，5（2）：87-91.

[4] 胡春華，周文斌，王毛蘭，等.鄱陽湖氮磷營養(yǎng)鹽變化特征及潛在性富營養(yǎng)化評價[J].湖泊科學，2010，22（5）：723-728.

[5] 陳曉玲，張 媛，張 琍，等.豐水期鄱陽湖水體中氮、磷含量分布特征[J].湖泊科學，2013，25（5）：643-648.

[6] 胡鴻鈞，魏印心.中國淡水藻類——系統(tǒng)，分類及生態(tài)[M].北京：科學出版社，2006.

[7] 黃祥飛.湖泊生態(tài)調(diào)查觀測與分析[M].北京：中國標準出版社，2000.

[8] 國家環(huán)境保護總局.水和廢水監(jiān)測分析方法[M].北京：中國環(huán)境科學出版社，2002.

[9] 楊麗標，韓小勇，孫 璞，等.巢湖藻類組成與環(huán)境因子典范對應(yīng)分析[J].農(nóng)業(yè)環(huán)境科學學報，2011，30（5）：952-958.

[10] 沈會濤，劉存歧.白洋淀浮游植物群落及其與環(huán)境因子的典范對應(yīng)分析[J].湖泊科學，2008，20（1）：773-779.

[11] 龐清江，李白英.東平湖水體富營養(yǎng)化評價[J].水資源保護，2003，19（5）：42-44.

[12] 汪 星，鄭丙輝，劉錄三，等.洞庭湖典型斷面藻類組成及其與環(huán)境因子典范對應(yīng)分析[J].農(nóng)業(yè)環(huán)境科學學報，2012，31（5）：995-1002.

[13] 李榮昉，張 穎.鄱陽湖水質(zhì)時空變化及其影響因素分析[J].水資源保護，2011，27（6）：9-13.

[14] 李亞蒙，趙 琦，馮廣平，等.白洋淀硅藻分布及其與水環(huán)境的關(guān)系[J].生態(tài)學報，2010，31（17）：4559-4570.

[15] 吳 珺，李 浩，曹德菊，等.巢湖東半湖藍藻水華暴發(fā)時空動態(tài)及成因[J].農(nóng)業(yè)環(huán)境科學學報，2013，32（10）：2035-2041.

[16] KAWAKITA M，MINAMI M，EGUCHI S et al. An introduction to the predictive technique AdaBoost with a comparison to generalized additive models[J].Fisheries Research，2005，76（3）：328-343.

[17] 王博雯，湯 祥，高 光，等.博斯騰湖細菌豐度時空分布及其與環(huán)境因子之間的關(guān)系[J].生態(tài)學報，2014，34（7）：1812-1821.

[18] 孫鑫鑫，劉惠榮，馮福應(yīng)，等.烏梁素海富營養(yǎng)化湖區(qū)浮游細菌多樣性及系統(tǒng)發(fā)育分析[J].生物多樣性，2009，17（5）：490-498.

[19] YUAN X C，HE L，YIN K D，et al. Bacterial distribution and nutrient limitation in relation to different water masses in the coastal and northwestern South China Sea in late summer[J].Continental Shelf Research，2011，31（11）：1214-1223.endprint