褚帆， 劉憲斌，劉占廣， 劉新蕾，張文亮，劉暢

（1.天津科技大學(xué) 天津市海洋資源與化學(xué)重點實驗室，天津 300457；2.渤海海洋監(jiān)測監(jiān)視管理基地籌備處，天津 300457；3.國家海洋信息中心，天津 300171）

渤海是我國唯一半封閉型陸架邊緣海，三面環(huán)陸，水深較淺。環(huán)渤海地區(qū)是我國北方經(jīng)濟活動的中心之一，渤海也是受沿岸人為活動影響最為明顯、污染最嚴(yán)重的海域之一，入海河流逐漸成為渤海近岸海域的主要污染源（王琪等，2009；Grimvall et al，2000；劉瓊瓊 等，2013）。

富營養(yǎng)化是指在一定的水體范圍內(nèi)，水體隨著營養(yǎng)物質(zhì)的增加而導(dǎo)致營養(yǎng)狀況發(fā)生變化的過程（Meyer-Reil et al，2000）。目前，已有許多學(xué)者（秦昌波 等，2006；Lundberg et al，2005） 提出了對海水富營養(yǎng)化評價的方法，如模糊綜合評價法、主成分-因子分析法、營養(yǎng)指數(shù)評價法、集對分析方法和灰色關(guān)聯(lián)度評價法等。劉蕓等（2014） 針對渤海海水水質(zhì)年際變化情況進行了調(diào)查分析；孫元敏等（2010） 和張繼民等（2008） 分別針對南亞熱帶海島和黃河口附近海域進行了富營養(yǎng)化評價。而對天津渤海灣海域水質(zhì)富營養(yǎng)化主要驅(qū)動因子進行分析的文章卻鮮見報道。

本研究通過對天津主要河流入?？谥車逗Ｓ虻暮Ｋ|(zhì)因子進行采樣調(diào)查，運用營養(yǎng)狀態(tài)質(zhì)量指數(shù)（Nutrient Quality Index，NQI） 和富營養(yǎng)化指數(shù)法（Eutrophication Index Method，EI） 對天津近岸海域的富營養(yǎng)化現(xiàn)狀進行了評價，并用主成分分析法（Principal Component Analysis） 得出影響海水富營養(yǎng)化的主要驅(qū)動因子，為控制該區(qū)域海水富營養(yǎng)化提供科學(xué)依據(jù)。

1 材料與方法

1.1 調(diào)查時間和區(qū)域

于2013年7月對天津主要入海河流（永定新河、海河、獨流減河、子牙新河） 入海口周圍海域海水進行采樣，采樣站位的坐標(biāo)范圍為117.5°E-118.2°E、38.5°N-39.2°N，站位分布如圖1 所示。樣品采集、貯存與運輸?shù)劝凑铡逗Ｑ蟊O(jiān)測規(guī)范》（GB17378-2007） 和《海洋調(diào)查規(guī)范》 （GB/T 12763-2007） 等要求進行。

1.2 水質(zhì)因子調(diào)查

本次調(diào)查的海水水質(zhì)因子包括鹽度、葉綠素a（Chl-a）、總氮（TN）、總磷（TP）、化學(xué)需氧量（CODMn）、銨氮（NH4-N）、亞硝酸鹽氮（NO2-N）、硝酸鹽氮（NO3-N）、活性磷酸鹽（PO4-P） 和水溫等10 個。監(jiān)測方法按照《海洋監(jiān)測規(guī)范》（GB17378-2007） 進行。

1.3 評價與分析方法

1.3.1 營養(yǎng)狀態(tài)質(zhì)量指數(shù)法

營養(yǎng)狀態(tài)質(zhì)量指數(shù)法（NQI） （陳于望 等，1999；邊佳胤等，2013） 是目前我國使用較多的評價方法，計算公式如下：

式中，CCOD、CDIN、CPO4-P、Chl-a 分別為海水CODMn、無機氮、活性磷酸鹽、葉綠素a 的監(jiān)測濃度；C′COD、C′DIN、C′PO4-P、Chl-a′ 分別為3.0 mg/L、0.3 mg/L、0.03 mg/L、5 mg/m3。

圖1 調(diào)查站位示意圖

根據(jù)營養(yǎng)質(zhì)量指數(shù)值將海域營養(yǎng)水平劃分為三級：NQI ≤2 時，營養(yǎng)等級為Ⅰ級，為貧營養(yǎng)化狀態(tài)；2 ＜NQI ≤3 時，營養(yǎng)等級為Ⅱ級，為中等營養(yǎng)化狀態(tài)；NQI ＞3 時，營養(yǎng)等級為Ⅲ級，為富營養(yǎng)化狀態(tài)（郭衛(wèi)東等，1998）。

1.3.2 富營養(yǎng)化指數(shù)法

目前，國內(nèi)較常用富營養(yǎng)化指數(shù)法（EI） （鄒景忠等，1983；于春燕等，2013） 對海水水質(zhì)進行分析評價，以便于與我國其它海域進行比較，公式如下：

式中，EI 為富營養(yǎng)化指數(shù)，CODMn為化學(xué)需氧量，DIN 為無機氮（NO2-N、NO3-N 和NH4-N 的總和），DIP 為無機磷。EI ≥1 時，海水達到富營養(yǎng)化，富營養(yǎng)化程度隨EI 值增高而增高：其中1 ≤EI ≤3 為輕度富營養(yǎng)化，3 ＜EI ≤9 為中度富營養(yǎng)化，EI ＞9 為重度富營養(yǎng)化（國家海洋局，2013）。

1.3.3 富營養(yǎng)化驅(qū)動因子的主成分分析

利用統(tǒng)計分析軟件SPSS 19.0 對海水水質(zhì)進行主成分分析，對結(jié)果進行KMO 檢驗以及Bartlett球形檢驗。當(dāng)KMO 檢驗數(shù)大于0.6，Bartlett 球形檢驗顯著性概率P 小于0.05 時，該組因子間相互獨立，可以進行主成分分析。

2 結(jié)果與討論

2.1 海水水質(zhì)現(xiàn)狀

海水水質(zhì)因子監(jiān)測結(jié)果如圖2。站位間水質(zhì)變化趨勢幾乎一致，調(diào)查海域約80.95%站位的無機氮和61.90 %站位的活性磷酸鹽為劣四類海水水質(zhì)，這些站位幾乎均位于入?？谥車Ｓ?。

對不同海域海水水質(zhì)因子對比見表1，與國內(nèi)其它海域相比，本研究海域的水質(zhì)污染水平高于江蘇、廈門、廣西等地近岸海域。對天津近岸海域年度水質(zhì)變化進行比較（表2），可以看出水質(zhì)有污染加劇的趨勢（闞文靜等，2013；尹翠玲等，2014），由于本文調(diào)查海域靠近河流入?？冢廴据^其它海域嚴(yán)重。

2.2 營養(yǎng)狀態(tài)質(zhì)量指數(shù)評價

運用營養(yǎng)狀態(tài)質(zhì)量指數(shù)法對海水水質(zhì)進行評價，NQI 值為2.54～12.26，均值是5.83。由圖3 可知，NQI 值由北塘入?？诤妥友佬潞尤牒？谙蜻h岸逐漸遞減，71.43 %站位海域達到富營養(yǎng)化水平。本研究結(jié)果與曹佳蓮等（2007） 結(jié)果一致，主要是受海水的旋轉(zhuǎn)流和往復(fù)流影響，營養(yǎng)鹽從近岸入?？谔幭蜻h岸擴散。

2.3 富營養(yǎng)化指數(shù)評價法

運用富營養(yǎng)化指數(shù)評價法對海水水質(zhì)進行評價，結(jié)果發(fā)現(xiàn)調(diào)查海域95.24%的站位處于重度富營養(yǎng)化狀態(tài)。如圖4 可以看出，EI 由北塘入?？诤妥友佬潞尤牒？谙蜻h岸海域逐漸變低。指數(shù)大于50 的站位分布于永定新河、海河以及子牙新河入海口周圍海域。

由圖3、4 可以得出，營養(yǎng)狀態(tài)指數(shù)評價和富營養(yǎng)化狀態(tài)指數(shù)評價分析結(jié)果一致，均是在北塘入?？谥車Ｓ蛭廴据^為嚴(yán)重，其次為子牙新河入海口周圍海域。2013年中國海洋環(huán)境狀況公報（國家海洋局，2013） 顯示，該區(qū)調(diào)查站位海水水質(zhì)狀況為四類或劣四類。劉學(xué)海（2010） 對塘沽-北塘-漢沽近岸海域進行調(diào)查，認(rèn)為北塘入?？谂c子牙新河入?？谒|(zhì)均嚴(yán)重超標(biāo)。朱文俊等（2013）對永定新河河口周圍海域水質(zhì)監(jiān)測結(jié)果表明永定新河河口附近的近岸海域有些站點達到重度污染等級。上述與本文調(diào)查結(jié)果一致。

2.4 主成分分析結(jié)果

對葉綠素a（Chl-a）、總氮（TN）、總磷（TP）、化學(xué)需氧量（CODMn）、銨氮（NH4-N）、亞硝酸鹽氮（NO2-N）、硝酸鹽氮（NO3-N）、活性磷酸鹽（PO4-P） 和水溫等9 個因子進行主成分分析。本次調(diào)查中監(jiān)測數(shù)據(jù)的KMO 統(tǒng)計量為 0.611，Bartlett 球形檢驗相伴概率為0.000，可以對導(dǎo)致海水富營養(yǎng)化的成因進行主成分分析。提取出主成分1 和2，方差貢獻率分別為56.06%和37.58%，累計貢獻率為93.64%。

這2 個主成分的變量載荷見表3。由表中可以看出主成分1 中CODMn、銨氮、亞硝酸鹽氮、硝酸鹽氮、總氮和水溫占的載荷較大，說明第一主成分反映的信息主要是溫度、氮營養(yǎng)鹽以及CODMn，營養(yǎng)鹽促進海洋浮游植物生長，CODMn代表的是水體中有機物，有機物含量高也會直接促進海洋浮游植物繁殖，而適宜的水溫是促進其生長的重要條件；主成分2 中葉綠素a 占的載荷較大，葉綠素a 指標(biāo)間接反映了海域的浮游植物量。

圖2 水質(zhì)因子：（A） CODMn、（B） 無機氮、（C） 總氮、（D） 總磷、（E） 活性磷酸鹽和（F） 葉綠素a

表1 國內(nèi)不同海域水質(zhì)參數(shù)比較

表2 天津近岸海域海水水質(zhì)參數(shù)年際變化

圖3 營養(yǎng)狀態(tài)質(zhì)量指數(shù)平面分布

對各水質(zhì)參數(shù)與主成分相關(guān)系數(shù)進行分析得出水體富營養(yǎng)化的主、次要驅(qū)動因子（表4）?？梢钥闯隹偟c主成分1 的相關(guān)系數(shù)最高，為影響該海域海水富營養(yǎng)化的主要驅(qū)動因子?？偭着c主成分2的相關(guān)系數(shù)最高，為影響該海域海水富營養(yǎng)化的次要驅(qū)動因子。鹽度與主成分1 相關(guān)性為-0.910，呈明顯負相關(guān)性，說明鹽度與營養(yǎng)鹽具有一定相關(guān)性，鹽度低的海域即入海口處海水富營養(yǎng)化程度相對較高，可能入海河流帶來的營養(yǎng)鹽含量較高，對近岸海水富營養(yǎng)化造成影響。

圖4 富營養(yǎng)化指數(shù)平面分布

表3 主成分分析結(jié)果

周斌等（2010） 對山東半島南部近岸海域海水水質(zhì)進行評價，總無機氮是影響海水富營養(yǎng)化的主要驅(qū)動因子，與本文的富營養(yǎng)化驅(qū)動因子為總氮的結(jié)論一致。2013年上半年天津市海洋環(huán)境質(zhì)量狀況報告（天津市海洋局，2013） 指出5月份永定新河、薊運河和潮白新河的水質(zhì)均劣于第Ⅴ類地表水水質(zhì)標(biāo)準(zhǔn)，說明可能由于入海徑流對入?？诤Ｋ柠}度進行了稀釋，且其攜帶的營養(yǎng)鹽對海水富營養(yǎng)化造成了較大影響。此結(jié)果與Redfield A C（1958）研究一致，鹽度作為影響海水富營養(yǎng)化的主要的物理因子，對海域富營養(yǎng)化狀況有著重要的影響。總之，天津近岸海域的海水富營養(yǎng)化是由各種因素共同作用的結(jié)果，而入海河流攜帶的營養(yǎng)鹽為影響富營養(yǎng)化狀況的最主要因素。

表4 各水質(zhì)參數(shù)與主成分的相關(guān)性分析

3 結(jié)論

（1） 調(diào)查海域80.95%的站位無機氮以及61.90%的站位活性磷酸鹽為劣四類海水水質(zhì)，主要位于入?？谥車Ｓ?。

（2） 營養(yǎng)狀態(tài)質(zhì)量指數(shù)評價與富營養(yǎng)化指數(shù)評價的結(jié)果趨勢相近，在北塘入?？谥車Ｓ蛑笖?shù)值最高，其次為子牙新河和海河周圍海域。

（3） 主成分分析結(jié)果表明，影響調(diào)查海域富營養(yǎng)化的驅(qū)動因子為總氮和總磷；而鹽度與主成分因子呈顯著負相關(guān)性，鹽度低的區(qū)域海水富營養(yǎng)化程度高，原因是入海河流水質(zhì)超標(biāo)，引起海水富營養(yǎng)化。

Grimvall A, 2000. Time scales of nutrient losses from land to sea a European perspective.Ecological Engineering,14 (4) :363-371.

Lundberg C, Lnnroth M, Von Numers M, et al, 2005. A multivariate assessment of coastal eutrophication. Examples from the Gulf of Finland, northern Baltic Sea. Marine Pollution Bulletin, 50 (11) :1185-1196.

Meyer-Reil L A, Koster M. Eutrophication of marine waters: effects on benthic microbial communities.Marine Pollution Bulletin,2000,41(1/6) :255-263.

Redfield A C, 1958. The biological control of chemical factors in the environment.American Scientist,46:561-600.

邊佳胤，袁林，王瓊，等，2013.洋山深水港海域水質(zhì)變化趨勢分析及富營養(yǎng)化評價.海洋通報，32（1）：107-112.

曹佳蓮，劉憲斌，劉占廣，等，2007.天津港南部海區(qū)水體富營養(yǎng)化水平評價.海洋湖沼通報，（S）：121-130.

陳于望，王憲，蔡明宏，1999.湄洲灣海域營養(yǎng)狀態(tài)評價.海洋環(huán)境科學(xué)，18（3）：39-42.

郭衛(wèi)東，章小明，楊逸萍，1998.中國近岸海域潛在性富營養(yǎng)化程度的評價.臺灣海峽，17 （1） ：64-69.

國家海洋局，2013.2013 中國海洋環(huán)境質(zhì)量公告.北京：國家海洋局.

姜晟，李俊龍，李旭文，等，2012.江蘇近岸海域富營養(yǎng)化現(xiàn)狀評價與成因分析.環(huán)境監(jiān)測管理與技術(shù)，24（4）：26-29.

闞文靜，張亞楠，張秋豐，等，2013.天津近岸海域水體富營養(yǎng)化特征及健康狀況分析.海洋湖沼通報，3：147-151.

賴俊翔，陳憲云，柯珂，等，2013.廣西防城港市近岸海域營養(yǎng)鹽分布及富營養(yǎng)化研究.海洋技術(shù)，32（3）：63-69，91.

劉瓊瓊，邵曉龍，劉紅磊，等，2013.天津市主要河流水質(zhì)及污染物入海通量.天津師范大學(xué)學(xué)報（自然科學(xué)版），33（2）：56-59.

劉學(xué)海，2010.渤海近岸水域環(huán)境污染狀況分析.環(huán)境保護科學(xué)，36（1）：14-18.

劉義峰，2013.三沙灣海水中磷的分布和富營養(yǎng)化研究.環(huán)境保護科學(xué)，39（4）：43-47.

劉蕓，馮劍豐，2014.天津近岸海域2004-2012年水化學(xué)環(huán)境的變化趨勢分析.海洋通報，33（2）：215-221.

歐陽玉蓉，王翠，李青生，等，2014.廈門灣海域營養(yǎng)鹽時空分布與富營養(yǎng)化狀況分析.福建農(nóng)業(yè)學(xué)報，29（1）：88-93.

秦昌波，鄭丙輝，秦延文，等，2006.渤海灣天津段海岸帶水環(huán)境質(zhì)量灰色關(guān)聯(lián)度評價.環(huán)境科學(xué)研究，19（6）：94-99.

施曉來，2013.近15a 來三門灣海域氮、磷營養(yǎng)鹽分布與富營養(yǎng)化狀況的變化分析.海洋學(xué)研究，31（4）：63-67.

孫元敏，陳彬，馬志遠，等，2010.南亞熱帶海島周邊海域富營養(yǎng)化評價及原因分析.海洋通報，29（5）：572-576.

天津市海洋局，2013.2013年上半年天津市海洋環(huán)境質(zhì)量狀況報告.天津，天津市海洋局.

王琪，紀(jì)朝彬，2009.渤海環(huán)境綜合治理的制度安排.中國海洋大學(xué)學(xué)報（社會科學(xué)版），（2）：1-3.

尹翠玲，張秋豐，劉洋，等，2014.2012年渤海灣赤潮監(jiān)控區(qū)赤潮發(fā)生及富營養(yǎng)化狀況.海洋湖沼通報，2014，1：137-142.

于春艷，梁斌，鮑晨光，等，2013.渤海富營養(yǎng)化現(xiàn)狀及趨勢研究.海洋環(huán)境科學(xué)，32（2）：175-177.

張繼民，劉霜，張琦，等，2008.黃河口附近海域營養(yǎng)鹽特征及富營養(yǎng)化程度評價.海洋通報，27（5）：65-72.

張萬磊，張永豐，張建樂，等，2014.北戴河赤潮監(jiān)控區(qū)營養(yǎng)鹽變化及富營養(yǎng)化狀況分析.海洋湖沼通報，1：143-147.

周斌，王悠，王進河，等，2010.山東半島南部近岸海域富營養(yǎng)化狀況的多元評價研究.海洋學(xué)報，32（2） ：128-137．

朱文俊，劉憲斌，李廈，等，2013.永定新河河口周圍水域環(huán)境質(zhì)量現(xiàn)狀評價.安徽農(nóng)業(yè)科學(xué)，41（5）：2243-2245.

鄒景忠，董麗萍，秦保平，1983.渤海灣富營養(yǎng)化和赤潮問題的初步探討.海洋環(huán)境科學(xué)，2（2）：41-54.