楊朝霞，李暢游，史小紅，孫標(biāo)，杜丹丹，全棟

內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院，內(nèi)蒙古 呼和浩特 010018

近年來，隨著經(jīng)濟(jì)的發(fā)展、城鎮(zhèn)化進(jìn)程的加劇及湖泊水資源不合理開發(fā)利用導(dǎo)致湖泊水體生態(tài)環(huán)境惡化，水體富營養(yǎng)化程度日趨嚴(yán)重，湖泊水體富營養(yǎng)化使得水體中藻類以及其它水生生物異常繁殖，破壞了湖泊水生生態(tài)系統(tǒng)平衡，造成水生生物資源浪費(fèi)，進(jìn)而影響湖泊水體功能（王蘇民等，1998）。湖泊水體的富營養(yǎng)化一直是全球最為突出的水環(huán)境問題之一，國內(nèi)外很多專家學(xué)者致力于湖泊富營養(yǎng)化研究（姜忠峰等，2011；田偉東等，2016；Coutinho et al.，2019；吳鋒等，2012；張志勇等，2017；張志紅等，2015），主要基于富營養(yǎng)化形成機(jī)理（孔繁翔等，2005），采用支持向量機(jī)（張成成等，2013）、依靠介質(zhì)和參數(shù)的分析（Lermontov et al.，2009；Purandara et al.，2012）與從經(jīng)濟(jì)發(fā)展因素對(duì)湖泊水質(zhì)影響角度出發(fā)（張民等，2015）進(jìn)行湖泊富營養(yǎng)化的評(píng)價(jià)研究。呼倫湖是我國北方高緯度寒旱區(qū)大型草原型湖泊，對(duì)我國北方經(jīng)濟(jì)的發(fā)展及生態(tài)環(huán)境的改善起著重要的作用，由于水文氣候條件的變化及人類活動(dòng)的雙重影響，呼倫湖近年來出現(xiàn)一定程度的水質(zhì)惡化與富營養(yǎng)化問題，直接威脅著湖泊生態(tài)健康發(fā)展。趙慧穎等（2007）對(duì)呼倫湖1959—2006 年的鹽度分析評(píng)價(jià)指出呼倫湖水體已處于富營養(yǎng)化狀態(tài)；梁麗娥等（2016）分析呼倫湖2006—2015 年富營養(yǎng)化趨勢(shì)指出呼倫湖水體富營養(yǎng)化程度呈現(xiàn)中度—重度—中度—輕度的變化趨勢(shì)。

本文以呼倫湖為研究區(qū)，采用綜合營養(yǎng)狀態(tài)指數(shù)法對(duì)2015—2017 年冰封期及非冰封期呼倫湖水體營養(yǎng)狀態(tài)進(jìn)行評(píng)價(jià)，并以各監(jiān)測(cè)點(diǎn)綜合營養(yǎng)指數(shù)均值作為標(biāo)準(zhǔn)，將綜合營養(yǎng)指數(shù)進(jìn)行離散化建立呼倫湖富營養(yǎng)化影響因子的非線性概率Probit 模型，分析湖泊水體富營養(yǎng)化主要水環(huán)境影響因子，并結(jié)合研究區(qū)水環(huán)境特征、社會(huì)狀況及人類活動(dòng)等影響因素進(jìn)行水體富營養(yǎng)化成因討論，本研究旨在明確呼倫湖冰封期及非冰封期水體的營養(yǎng)狀態(tài)，并對(duì)富營養(yǎng)化成因及其主要水環(huán)境影響因子進(jìn)行分析，以便對(duì)湖泊水體富營養(yǎng)化治理提供一定的理論參考。

1 材料與方法

1.1 研究區(qū)概況

呼倫湖（116°58′—117°48′E，48°33′—49°20′N）地處內(nèi)蒙古自治區(qū)滿洲里市、新巴爾虎左旗及新巴爾虎右旗之間（如圖1），是寒旱區(qū)極為罕見的具有生物多樣性和生態(tài)多功能的草原型湖泊濕地生態(tài)系統(tǒng)。湖區(qū)呈以東北至西南方向?yàn)檩S的不規(guī)則斜長方形，長度為93 km，最大寬度41 km，最大水深達(dá)18 m，蓄水量達(dá)120×108m3。氣候類型屬中溫帶大陸性氣候，流域多年平均氣溫0.16 ℃，多年平均降水量264.3 mm，多年平均蒸發(fā)量1 411 mm；湖泊受西北冷空氣與蒙古高壓控制，冬季寒冷而漫長，冰封期長達(dá)6 個(gè)月；湖泊水量補(bǔ)給除降水外，主要來自克魯倫河與烏爾遜河的外源補(bǔ)給（李暢游等，2016）。

圖1 研究區(qū)及采樣點(diǎn)分布 Fig. 1 Location of Hulun Lake and sample points

1.2 數(shù)據(jù)來源及檢測(cè)方法

根據(jù)對(duì)呼倫湖多年連續(xù)監(jiān)測(cè)，選取最具代表性的13 個(gè)監(jiān)測(cè)點(diǎn)，于2015 年、2016 年與2017 年1月冰封期及7 月非冰封期對(duì)湖泊水體進(jìn)行采樣檢測(cè)分析，測(cè)定指標(biāo)包括水深（H）、水溫（T）、pH值、電導(dǎo)率（EC）、溶解氧（DO）、氧化還原電位（OPR）、溶解性固體總量（TDS）、鹽度（S）、總氮（TN）、總磷（TP）、化學(xué)需氧量（COD）、葉綠素（Chl.a）及透明度（SD）。其中總磷（TP）采用鉬銻抗分光光度法測(cè)定，總氮（TN）采用堿性過硫酸鉀消解-紫外分光光度法測(cè)定，化學(xué)需氧量（COD）采用硫酸亞鐵銨滴定法測(cè)定，葉綠素（Chl.a）采用丙酮萃取分光光度計(jì)法測(cè)定，透明度（SD）采用塞氏盤測(cè)定，其余常規(guī)水質(zhì)指標(biāo)由瑞士梅特勒多參數(shù)儀測(cè)定（國家環(huán)境保護(hù)總局，2002）。每一種指標(biāo)重復(fù)測(cè)定3 次，取平均值作為測(cè)定值。

1.3 研究方法

1.3.1 水體富營養(yǎng)化評(píng)價(jià)

呼倫湖水體富營養(yǎng)化程度采用綜合營養(yǎng)狀態(tài)指數(shù)法進(jìn)行評(píng)價(jià)（張壘等，2015），綜合營養(yǎng)狀態(tài)指數(shù)法以葉綠素作為基準(zhǔn)參數(shù)，選取水質(zhì)參數(shù)的相關(guān)權(quán)重，通過計(jì)算各水質(zhì)參數(shù)的營養(yǎng)狀態(tài)指數(shù)并最終加權(quán)平均求得水體的綜合營養(yǎng)狀態(tài)指數(shù)。根據(jù)綜合營養(yǎng)狀態(tài)指數(shù)計(jì)算值，對(duì)營養(yǎng)狀態(tài)評(píng)價(jià)采用0—100 的一系列連續(xù)數(shù)字進(jìn)行分級(jí)，其中TLI(Σ)＜30為貧營養(yǎng)，30≤TLI(Σ)≤50 為中營養(yǎng)，50＜TLI(Σ)≤60為輕度富營養(yǎng)，60＜TLI(Σ)≤70 為中度富營養(yǎng)，TLI(Σ)＞70 為重度富營養(yǎng)，同一營養(yǎng)狀態(tài)下，其綜合營養(yǎng)指數(shù)值越大，富營養(yǎng)程度越重。其中綜合營養(yǎng)指數(shù)計(jì)算公式如式（1）所示：

式中，TLI(Σ)為綜合營養(yǎng)狀態(tài)指數(shù)；Wj為第j種參數(shù)的營養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重；TLIj為第j 種參數(shù)的營養(yǎng)狀態(tài)指數(shù)，m 為評(píng)價(jià)參數(shù)的個(gè)數(shù)。其中，各營養(yǎng)狀態(tài)指數(shù)計(jì)算如公式（2）—（6）所示：

式中，Chl.a 單位為mg·m-3，SD 單位為m，其他指標(biāo)單位均為mg·L-1。Chl.a 作為基準(zhǔn)參數(shù)，則第j 種參數(shù)歸一化的相關(guān)權(quán)重計(jì)算公式如式（7）所示：

式中，rij為第j 種參數(shù)與Chl.a 的相關(guān)系數(shù)，本文選取的水質(zhì)參數(shù)為Chl.a、TP、TN、COD 和SD，相關(guān)系數(shù)如表1 所示。

1.3.2 富營養(yǎng)化影響因子模型

呼倫湖富營養(yǎng)狀態(tài)的水環(huán)境影響因子采用非線性概率Probit 模型進(jìn)行研究，將湖泊富營養(yǎng)化程度作為因變量并進(jìn)行離散化，其余各水環(huán)境影響因子作為自變量建立模型，通過對(duì)模型進(jìn)行邊際效應(yīng)分析，計(jì)算各水環(huán)境影響因子變化時(shí)相應(yīng)因變量的響應(yīng)概率，進(jìn)而分析各影響因子對(duì)湖泊營養(yǎng)狀態(tài)的影響程度（金光，2013；許羅丹等，2014）。Probit模型是二元離散選擇模型，其因變量只有兩個(gè)虛擬數(shù)值0 和1，代表因變量的兩種狀態(tài)，模型估計(jì)采用最大似然法估計(jì)，其模型分布為式（8）：

表1 中國湖泊（水庫）部分參數(shù)與Chl.a 的相關(guān)關(guān)系 Table 1 Table 1 The correlation between Chl.a and the some parameters of lake (reservoir) water in China

式中，Pi為因變量取值為1 時(shí)的概率， Xi*為自變量，β 為自變量的回歸系數(shù)。

2 結(jié)果與分析

2.1 呼倫湖冰封期與非冰封期營養(yǎng)狀態(tài)特征

圖2 2015—2017 年呼倫湖冰封期營養(yǎng)狀態(tài)分級(jí)圖 Fig. 2 Eutrophication distribution during icebound period in 2015-2017

圖3 2015—2017 年呼倫湖非冰封期營養(yǎng)狀態(tài)分級(jí)圖 Fig. 3 Eutrophication distribution during non-icebound period in 2015-2017

對(duì)2015 年、2016 年與2017 年1 月冰封期及7月非冰封期湖泊水體進(jìn)行富營養(yǎng)化狀態(tài)評(píng)價(jià)，湖泊水體營養(yǎng)狀態(tài)分級(jí)如圖2 和圖3 所示。由圖2 可知，呼倫湖冰封期水體富營養(yǎng)化程度自2015—2017 年總體呈現(xiàn)明顯下降趨勢(shì)，其中2015 年1 月全湖82.8%區(qū)域呈現(xiàn)輕度富營養(yǎng)狀態(tài)，僅17.2%區(qū)域呈現(xiàn)中度富營養(yǎng)狀態(tài)且主要集中于湖區(qū)西南部；2016 年1 月全湖84.2%區(qū)域呈現(xiàn)輕度富營養(yǎng)狀態(tài)，僅15.8%區(qū)域呈現(xiàn)中度富營養(yǎng)狀態(tài)且主要集中于湖區(qū)西南部；2017 年1 月全湖89.4%區(qū)域呈現(xiàn)輕度富營養(yǎng)狀態(tài)，僅10.6%區(qū)域呈現(xiàn)中度富營養(yǎng)狀態(tài)且主要集中于湖區(qū)東北部。由圖3 可知，呼倫湖非冰封期水體富營養(yǎng)化程度自2015—2017 年總體呈現(xiàn)先加重后下降趨勢(shì)。其中2015 年呼倫湖全湖呈現(xiàn)輕度富營養(yǎng)狀態(tài)；2016 年全湖98.7%區(qū)域呈輕度富營養(yǎng)狀態(tài)，且在湖區(qū)東北部1.3%區(qū)域富營養(yǎng)化程度加重呈現(xiàn)中度富營養(yǎng)狀態(tài)；2017 年全湖98.8%區(qū)域呈現(xiàn)輕度富營養(yǎng)狀態(tài)，且在E8 點(diǎn)附近1.2%區(qū)域內(nèi)水體無富營養(yǎng)化，呈現(xiàn)中營養(yǎng)狀態(tài)。對(duì)比呼倫湖冰封期與非冰封期富營養(yǎng)狀態(tài)可知，湖泊水體富營養(yǎng)化程度冰封期較非冰封期整體嚴(yán)重，主要原因?yàn)楹捶鈨龊?，冰凍濃縮作用使?fàn)I養(yǎng)鹽向下排至水體造成水體中氮、磷等營養(yǎng)鹽濃度偏高，造成水體富營養(yǎng)化程度加重；此外，非冰封期呼倫湖由于克魯倫河、烏爾遜河和地下水對(duì)湖泊的充沛補(bǔ)給，湖泊蓄水量增加使?fàn)I養(yǎng)物質(zhì)濃度降低，水體富營養(yǎng)化程度因此明顯下降。

呼倫湖水體富營養(yǎng)化主要集中在湖區(qū)西南部（G2，I2，H3）以及湖區(qū)東北部，主要由于湖區(qū)西南部有一條主要入湖河流克魯倫河，如圖4 所示克魯倫河多受畜牧業(yè)所產(chǎn)生的氮、磷等有機(jī)物污染，且流經(jīng)人類聚居區(qū)，同時(shí)沿岸又有熟皮廠和硝礦，向水體排放的氨氮、農(nóng)藥、氟化物、砷、汞、酚等污染物成為呼倫湖水體富營養(yǎng)化的主要外源輸入；湖區(qū)東北部緊靠扎賚諾爾區(qū)，同時(shí)湖區(qū)東北端有小河口漁場(chǎng)，也是呼倫湖重要的旅游景點(diǎn)，人類活動(dòng)頻繁對(duì)湖泊水質(zhì)造成一定影響，成為湖泊水體富營養(yǎng)化的重要影響因素（付堯，2016；丁燕青等，2011；薛洋，2015）。

圖4 克魯倫河2016 年6 月總氮總磷濃度圖 Fig. 4 Kluron River June 2016 total nitrogen total phosphorus concentration map

2.2 呼倫湖富營養(yǎng)化影響因子模型

呼倫湖水體富營養(yǎng)化是多重水環(huán)境因子共同作用的結(jié)果，其中氮、磷等有機(jī)物和營養(yǎng)鹽是湖泊水體富營養(yǎng)化的直接影響因子，而湖泊水體理化性質(zhì)、水環(huán)境因子與湖泊環(huán)境特征也是造成水體富營養(yǎng)化的重要影響因子?，F(xiàn)采用綜合營養(yǎng)指數(shù)對(duì)2015—2017 年呼倫湖水體進(jìn)行富營養(yǎng)化評(píng)價(jià)并將其作為因變量，并以水深、水溫、pH 值、電導(dǎo)率、溶解氧、氧化還原電位、溶解性固體總量、鹽度、總氮、總磷、化學(xué)需氧量等11 各水環(huán)境指標(biāo)為自變量，研究呼倫湖水體富營養(yǎng)化與水環(huán)境因子之間的響應(yīng)關(guān)系?，F(xiàn)計(jì)算2015—2017 年呼倫湖各監(jiān)測(cè)點(diǎn)綜合營養(yǎng)指數(shù)均值為57.24，并將此值作為標(biāo)準(zhǔn)，將2015—2017 年各監(jiān)測(cè)點(diǎn)富營養(yǎng)化程度離散化，其中監(jiān)測(cè)點(diǎn)綜合營養(yǎng)指數(shù)大于均值表示水體富營養(yǎng)化程度較重，取值為“1”，監(jiān)測(cè)點(diǎn)綜合營養(yǎng)指數(shù)小于均值表示水體富營養(yǎng)化程度較輕，取值為“0”。其余各水環(huán)境因子均為連續(xù)變量，取2015—2017 年各監(jiān)測(cè)點(diǎn)實(shí)測(cè)值，通過Probit模型建立呼倫湖營養(yǎng)狀態(tài)影響因子模型，結(jié)果如式（9）所示：

Xi*β =-1.838+0.031 4T+0.019 9TDS+0.825S- 0.916pH+0.024OPR+0.495H+0.198DO-0.012 4EC+ 0.052 8COD+0.369TN+5.485TP

式中，Pi表示湖泊水體呈現(xiàn)富營養(yǎng)化程度較重狀態(tài)時(shí)的概率。

對(duì)營養(yǎng)狀態(tài)影響因子模型進(jìn)行檢驗(yàn)，擬合優(yōu)度為0.458，擬合結(jié)果較好，LR-statistic 檢驗(yàn)值為32.84，表明模型總體顯著，建模結(jié)果比較理想。對(duì)模型進(jìn)行邊際效應(yīng)分析，求解自變量變化時(shí)相應(yīng)因變量變化的響應(yīng)概率，結(jié)果如圖5 所示。由圖4 可知，營養(yǎng)狀態(tài)影響因子模型中各影響因子的響應(yīng)概率 大 小 排 序 為 TP＞pH＞S＞H＞TN＞DO＞COD＞T＞ ORP＞TDS＞EC。其中總磷（TP）響應(yīng)概率為65.07%，對(duì)湖泊水體富營養(yǎng)化程度影響最大，說明磷是呼倫湖水體富營養(yǎng)化的直接限制因子；pH 值與鹽度的響應(yīng)概率分別為10.87%和9.79%，對(duì)湖泊水體富營養(yǎng)化程度影響次之；水深、總氮和溶解氧響應(yīng)概率分別為5.87%、4.38%和2.35%，對(duì)湖泊水體富營養(yǎng)化程度也有一定影響，而其余水環(huán)境因子則對(duì)湖泊水體富營養(yǎng)化程度影響相對(duì)較小。

圖5 Probit 模型影響因子響應(yīng)概率 Fig. 5 Influence factor response probability of Probit model

3 討論

通過對(duì)呼倫湖2015—2017 年冰封期與非冰封期水體富營養(yǎng)狀態(tài)評(píng)價(jià)，并通過呼倫湖營養(yǎng)狀態(tài)影響因子Probit 模型分析湖泊水體富營養(yǎng)化的主要水環(huán)境影響因子，并結(jié)合呼倫湖的湖泊環(huán)境特征及實(shí)際狀況可知引起湖泊水體富營養(yǎng)化的主要影響因素為以下幾個(gè)方面。

3.1 冰封期與非冰封期對(duì)水體富營養(yǎng)化的影響

通過對(duì)呼倫湖2015—2017 年冰封期及非冰封期水體營養(yǎng)狀態(tài)對(duì)比可知，從冰封期到非冰封期水體各富營鹽增大，因?yàn)殡S著冰封期湖面冰體形成，營養(yǎng)鹽從冰體不斷向水體中遷移，造成水質(zhì)惡化，富營養(yǎng)化程度不斷加重（高寧，2018）。此外，非冰封期呼倫湖由于有降雨及克魯倫河、烏爾遜河和地下水對(duì)湖泊的充沛補(bǔ)給，湖泊蓄水量增加使?fàn)I養(yǎng)物質(zhì)濃度降低，水體富營養(yǎng)化程度因此明顯下降。

3.2 水環(huán)境因子對(duì)富營養(yǎng)化影響

通過對(duì)呼倫湖水體營養(yǎng)狀態(tài)影響因子建模分析可知，總磷（TP）、pH 值、鹽度（S）、水深（H）、總氮（TN）和溶解氧（DO）等6 個(gè)水質(zhì)指標(biāo)對(duì)呼倫湖水體營養(yǎng)狀態(tài)影響程度最大，對(duì)各水環(huán)境因子與綜合營養(yǎng)指數(shù)建立相關(guān)關(guān)系如圖5 所示，由圖6可知總磷（TP）、pH 值、鹽度（S）、水深（H）、總氮（TN）和溶解氧（DO）與綜合營養(yǎng)指數(shù)間的相關(guān)系數(shù)分別為0.746 5、0.773 9、0.841 1、0.775 6、0.707 5 和0.682 1，各水質(zhì)指標(biāo)與綜合營養(yǎng)指數(shù)之間存在較強(qiáng)的線性相關(guān)關(guān)系。根據(jù)各水質(zhì)指標(biāo)與綜合營養(yǎng)指數(shù)間變化趨勢(shì)關(guān)系可知，綜合營養(yǎng)指數(shù)隨總磷、pH 值、鹽度與總氮的值增加而呈現(xiàn)增大趨勢(shì)，即總磷、pH 值、鹽度與總氮含量的增加，水體富營養(yǎng)化程度加重，主要由于總磷和總氮是水體富營養(yǎng)化的直接影響因子，且磷對(duì)呼倫湖水體富營養(yǎng)化起著限制性作用，pH 值與鹽度是水體理化性質(zhì)的直接反映，其中pH 值決定水體酸堿程度，鹽度反映湖泊水體離子狀況，水體理化指標(biāo)對(duì)水生動(dòng)植物、微生物生長、繁殖均產(chǎn)生重要影響，也對(duì)水體水質(zhì)狀況及營養(yǎng)程度有著重要影響（丁燕青，2011）。而綜合營養(yǎng)指數(shù)隨水深和溶解氧的值增加而呈現(xiàn)減小趨勢(shì)，即水深增大和溶解氧含量增加，湖泊水體富營養(yǎng)化程度減輕，主要由于溶解氧反映水體富氧程度，溶解氧含量的增加，促進(jìn)了水生植物的光合作用、降低厭氧型水生微生物繁殖，有利于水體水質(zhì)的改善；水深則反映湖泊水位變動(dòng)與蓄水量變化情況，水動(dòng)力條件的變化和蓄水量的變化對(duì)湖泊污染物的釋放、遷移、稀釋與轉(zhuǎn)化產(chǎn)生著重要影響，進(jìn)而影響湖泊水體的營養(yǎng)程度（陳義中，2018）。

3.3 人類活動(dòng)對(duì)富營養(yǎng)化影響

圖6 水環(huán)境因子與綜合營養(yǎng)指數(shù)關(guān)系 Fig. 6 Relationship between water environment factors and comprehensive nutrition index

呼倫湖流域內(nèi)生活著大約400 戶左右的牧民，生活中大量使用的洗滌劑等大多含有磷酸鹽，由此而產(chǎn)生的生活污水未經(jīng)水質(zhì)處理，直接排入呼倫湖引起湖泊水體含磷量的增加，造成水體富營養(yǎng)化。另外，呼倫湖湖區(qū)東北部的扎賚諾爾區(qū)也是呼倫湖的主要污染源，捕魚業(yè)是扎賚諾爾區(qū)經(jīng)濟(jì)發(fā)展的第一產(chǎn)業(yè)，漁業(yè)所占比重高達(dá)70%以上，當(dāng)?shù)鼐用袷芙?jīng)濟(jì)利益的驅(qū)使，通常采用高投入、高產(chǎn)出的方式來實(shí)現(xiàn)高收入，忽視對(duì)養(yǎng)殖水域的水環(huán)境生態(tài)平衡的保護(hù)，以致養(yǎng)殖水域環(huán)境快速受到污染并進(jìn)步惡化，導(dǎo)致呼倫湖水體出現(xiàn)營養(yǎng)物質(zhì)的過量富集而使水體呈現(xiàn)富營養(yǎng)化現(xiàn)象（薛洋，2015）。其次為農(nóng)業(yè)和采礦業(yè)也是扎賚諾爾區(qū)經(jīng)濟(jì)發(fā)展的重要產(chǎn)業(yè)，扎賚諾爾區(qū)有耕地面積117 hm2，研究發(fā)現(xiàn)磷在農(nóng)業(yè)污染中的貢獻(xiàn)率為75%，農(nóng)田中的磷通過降水和排水造成的地表徑流和土壤侵蝕從土體進(jìn)入水體（王森等，2013）；扎煤公司煤炭年產(chǎn)量為1567 萬噸，煤炭在開采過程中向環(huán)境排放大量廢棄物通過干、濕沉降與地表徑流作用，將可溶性氮、磷污染物帶入地表水體；或通過降水淋洗滲入地下水，地下水與湖水間通過補(bǔ)給過程，引起湖泊水體氮、磷含量升高，導(dǎo)致水體富營養(yǎng)化（李云鵬等，2017）。

3.4 不同類型放牧草場(chǎng)過度放牧對(duì)富營養(yǎng)化影響

呼倫湖地處呼倫貝爾大草原，流域內(nèi)畜牧業(yè)發(fā)達(dá)，牧民的過度放牧引起草場(chǎng)面積的減少、草原的退化與生態(tài)平衡的破壞。通過對(duì)比不同放牧類型草場(chǎng)中的磷元素流失情形發(fā)現(xiàn)，當(dāng)降雨強(qiáng)度為75 mm·h-1時(shí)，傳統(tǒng)放牧草場(chǎng)的徑流總磷流失負(fù)荷為5.69 mg·m-1草庫倫草場(chǎng)的徑流總磷流失負(fù)荷為2.73 mg·m-1，不放牧草場(chǎng)的徑流總磷流失負(fù)荷為0.45 mg·m-1，其中傳統(tǒng)放牧草場(chǎng)和不放牧草場(chǎng)徑流中磷元素以顆粒態(tài)流失為主，而草庫倫草場(chǎng)以溶解態(tài)磷流失為主，并且降雨強(qiáng)度越大，其溶解態(tài)磷元素流失所占的比例越大，說明放牧的草場(chǎng)磷元素水土流失較大，進(jìn)而造成地表水體富營養(yǎng)化加重（趙偉等，2011）。

4 結(jié)論

通過綜合營養(yǎng)指數(shù)法對(duì)呼倫湖2015—2017 年冰封期和非冰封期水體富營養(yǎng)狀態(tài)評(píng)價(jià)分析可知：

（1）呼倫湖2015—2017 年冰封期內(nèi)水體富營養(yǎng)化程度呈現(xiàn)逐年增加變化特征，而非冰封期內(nèi)水體富營養(yǎng)化程度總體呈現(xiàn)2016 年富營養(yǎng)化程度最大2017 年最小的特征，且冰封期內(nèi)水體富營養(yǎng)狀態(tài)程度較非冰封期嚴(yán)重。

（2）通過建立呼倫湖富營養(yǎng)化影響因子Probit模型分析可知：呼倫湖水體富營養(yǎng)化的主要水環(huán)境影響因子是總磷、pH 值、鹽度、水深、總氮和溶解氧，且總磷因子的響應(yīng)概率為65.7%，對(duì)水體富營養(yǎng)化程度影響最大，說明磷指標(biāo)是呼倫湖水體富營養(yǎng)化的重要限制因子。

（3）通過水質(zhì)指標(biāo)與綜合營養(yǎng)指數(shù)間相關(guān)性分析可知：總磷、pH 值、鹽度、水深、總氮和溶解氧與綜合營養(yǎng)指數(shù)間的相關(guān)系數(shù)分別為0.746 5、0.773 9、0.841 1、0.775 6、0.707 5 和0.682 1，相關(guān)關(guān)系顯著，表明可以通過控制湖泊水環(huán)境因子平衡來改善水體營養(yǎng)狀態(tài)。