曹金玲,許其功,2*,席北斗,2,高如泰,姜 磊,丁京濤,毛敬英,3,孫雪松,4 (.中國環(huán)境科學(xué)研究院,水環(huán)境系統(tǒng)工程研究室,北京 0002；2.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室,北京 0002；3.西南交通大學(xué)環(huán)境科學(xué)與工程學(xué)院,四川 成都 6003；4.重慶交通大學(xué)河海學(xué)院,重慶 400074)

第二階梯湖泊富營養(yǎng)化自然地理因素及效應(yīng)

曹金玲1,許其功1,2*,席北斗1,2,高如泰1,姜 磊1,丁京濤1,毛敬英1,3,孫雪松1,4(1.中國環(huán)境科學(xué)研究院,水環(huán)境系統(tǒng)工程研究室,北京 100012；2.中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室,北京 100012；3.西南交通大學(xué)環(huán)境科學(xué)與工程學(xué)院,四川 成都 610031；4.重慶交通大學(xué)河海學(xué)院,重慶 400074)

為了明晰影響我國地形地貌第二級階梯湖泊富營養(yǎng)化的主要自然地理驅(qū)動因素,為湖泊富營養(yǎng)化的有效治理提供參考,采用曲線回歸方法分析了位于我國第二級階梯的云南和新疆湖泊中Chl-a以及營養(yǎng)狀態(tài)綜合指數(shù)(TLI)與自然地理特征的相關(guān)性.結(jié)果表明,平均水深、海拔以及年均溫與TLI指數(shù)具有較好的相關(guān)性,擬合模型分別符合三次/對數(shù)模型、二次模型和三次模型.Chl-a與年平均日照時數(shù)具有較好的相關(guān)性,擬合模型符合二次/三次模型.分析富營養(yǎng)化效應(yīng),藻類生長對透明度的影響以及營養(yǎng)物質(zhì)對浮游藻類的生態(tài)效應(yīng)(Chl-a/TP),結(jié)果發(fā)現(xiàn),新疆湖泊中Chl-a與透明度呈倒數(shù)關(guān)系,云南湖泊中Chl-a與透明度呈冪函數(shù)關(guān)系;Chl-a與總磷的線性方程的斜率分別為177.595(云南)和222.758(新疆).影響我國第二級階梯湖泊富營養(yǎng)化的主要自然地理驅(qū)動因素是湖泊的平均水深、海拔、年均溫度及年平均日照時數(shù),透明度與Chl-a含量的相關(guān)性以及營養(yǎng)物質(zhì)TP對浮游藻類的生態(tài)效應(yīng)存在區(qū)域差異.

湖泊富營養(yǎng)化；驅(qū)動因素；富營養(yǎng)化效應(yīng)；綜合營養(yǎng)狀態(tài)指數(shù)；曲線估計；擬合模型

隨著社會經(jīng)濟的發(fā)展,人類過度開發(fā)利用湖泊資源并不斷向湖泊中排放污水,造成氮、磷等營養(yǎng)物質(zhì)過剩、江湖關(guān)系阻斷[1]、流域清水產(chǎn)流機制破壞[2],使湖泊自凈能力下降或者喪失,導(dǎo)致水體中藻類以及其他水生生物異常繁殖,水體透明度和溶解氧下降[3],加速湖泊由貧營養(yǎng)到中營養(yǎng)到富營養(yǎng)、沼澤化,最后走向消亡的過程[4-5].為了減緩日趨嚴(yán)重的富營養(yǎng)化現(xiàn)象,制定湖泊富營養(yǎng)化控制標(biāo)準(zhǔn)已成為治理湖泊富營養(yǎng)化的當(dāng)務(wù)之急,而湖泊營養(yǎng)物基準(zhǔn)是制定湖泊富營養(yǎng)化控制標(biāo)準(zhǔn)的基礎(chǔ)和依據(jù).目前,我國尚無明確的湖泊營養(yǎng)物基準(zhǔn)和富營養(yǎng)化控制標(biāo)準(zhǔn)[6],湖泊富營養(yǎng)化的管理與控制仍然參照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》或者其他發(fā)達(dá)國家的相關(guān)基準(zhǔn)值.但我國湖泊數(shù)量眾多、分布廣泛、類型多樣、成因復(fù)雜[3,7],不同區(qū)域湖泊的富營養(yǎng)化成因、類型、演變過程以及物理、化學(xué)、生物學(xué)特性等方面存在顯著差異,致使我國湖泊營養(yǎng)物水平和富營養(yǎng)化效應(yīng)具有明顯的區(qū)域差異性.因此,對處于不同區(qū)域、營養(yǎng)狀態(tài)各異的湖泊,采用同一水平的營養(yǎng)物基準(zhǔn)和富營養(yǎng)化控制標(biāo)準(zhǔn),缺乏針對性、科學(xué)性和指導(dǎo)性.因此,從分區(qū)控制和科學(xué)管理出發(fā),把我國的湖泊按其生態(tài)系統(tǒng)特征進(jìn)行科學(xué)分區(qū),針對不同分區(qū)湖泊系統(tǒng)特點、生態(tài)特征和營養(yǎng)物生態(tài)效應(yīng),結(jié)合湖泊參照狀態(tài),分別制定相應(yīng)的營養(yǎng)物基準(zhǔn)和富營養(yǎng)化控制標(biāo)準(zhǔn),是目前我國湖泊富營養(yǎng)化控制的主要任務(wù)之一.對我國湖泊進(jìn)行科學(xué)分區(qū)的基礎(chǔ)是明晰處于不同區(qū)域的湖泊富營養(yǎng)化效應(yīng)及驅(qū)動因素的差異性,然而對于全國大尺度湖泊富營養(yǎng)化區(qū)域差異性的研究資料較為罕見.本研究針對位于我國地形地貌第二級階梯的湖泊,分析影響這一區(qū)域湖泊營養(yǎng)狀態(tài)的主要自然地理因素,為探討全國大尺度湖泊富營養(yǎng)化驅(qū)動因素的區(qū)域差異性以及營養(yǎng)物生態(tài)分區(qū)提供基礎(chǔ)和依據(jù),為實現(xiàn)湖泊富營養(yǎng)化的“分區(qū)、分類、分級、分期”管理提供技術(shù)支持.

1 材料與方法

1.1 湖泊自然地理要素數(shù)據(jù)的收集

湖泊的水深(Depth)及所在區(qū)域的海拔(Altitude)、年均溫(Tem)、年平均降水量(AAP)、年平均日照時數(shù)(ASH)、年平均無霜期(AFFP)以及年平均蒸發(fā)量(AAE)7項自然地理要素及湖盆形態(tài)的數(shù)據(jù)來源于《中國湖泊志》[7]以及其他相關(guān)的文獻(xiàn).

1.2 水質(zhì)監(jiān)測

CODMn、TP、TN的測定采用《水和廢水監(jiān)測分析方法》[8]的方法進(jìn)行測定.透明度(SD)采用塞氏盤測定.Chl-a的測定采用丙酮萃取分光光度計法.

1.3 統(tǒng)計分析

利用 SPSS17.0中的曲線估計功能,分別對營養(yǎng)狀態(tài)綜合指數(shù)(TLI)與平均水深、海拔、年均溫、年平均降水量、年平均日照時數(shù)、年平均無霜期以及年平均蒸發(fā)量進(jìn)行回歸分析,探討影響該階梯內(nèi)湖泊富營養(yǎng)化的主要自然地理驅(qū)動因素.

2 結(jié)果與討論

2.1 第二階梯主要湖泊自然地理特征及營養(yǎng)狀態(tài)評價

位于中國地形地貌第二級階梯的省份主要包括:新疆維吾爾自治區(qū)、內(nèi)蒙古自治區(qū)、寧夏回族自治區(qū)、甘肅省、山西省、四川省、陜西省、重慶市、云南省、貴州省.湖泊分布比較集中的省份是新疆維吾爾自治區(qū)、內(nèi)蒙古自治區(qū)和云南省.利用綜合營養(yǎng)狀態(tài)指數(shù)法評價以上3個省份典型湖泊的營養(yǎng)狀態(tài),評價項目選取反映水體營養(yǎng)程度的5項指標(biāo):葉綠素a(Chl-a)、總磷(TP)、總氮(TN)、透明度(SD)、高錳酸鹽指數(shù)(CODMn).綜合營養(yǎng)狀態(tài)指數(shù)計算公式為:

式中:TLI(Σ)為綜合營養(yǎng)狀態(tài)指數(shù);Wj為第j種參數(shù)的營養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重;TLI(j)為第j種參數(shù)的營養(yǎng)狀態(tài)指數(shù).

以Chl-a作為基準(zhǔn)參數(shù),則第j種參數(shù)的歸一化的相關(guān)權(quán)重計算公式為:

式中:rij為第j種參數(shù)與基準(zhǔn)參數(shù)Chl-a的相關(guān)系數(shù);m為評價參數(shù)的個數(shù).營養(yǎng)狀態(tài)指數(shù)計算式為:式中:C(j)為第 j種污染物的濃度;其他污染指標(biāo)的參數(shù)值如表1所示.

采用 0～100的一系列連續(xù)數(shù)字對湖泊營養(yǎng)狀態(tài)進(jìn)行分級: TLI(Σ)＜30為貧營養(yǎng);30≤TLI(Σ)≤50為中營養(yǎng); TLI(Σ) ＞50為富營養(yǎng).評價結(jié)果如表2所示.其中內(nèi)蒙古的烏梁素海,云南的淺水湖泊滇池、星云湖、異龍湖、長橋海、大屯海和個舊湖處于富營養(yǎng)狀態(tài),其他湖泊均處于貧營養(yǎng)和中營養(yǎng)狀態(tài).

表1 營養(yǎng)狀態(tài)指數(shù)計算公式中各水質(zhì)指標(biāo)的參數(shù)值Table 1 Parameter values of water qualities in TLI equation

表2 中國第二級階梯典型湖泊自然地理特征及營養(yǎng)狀態(tài)評價Table 2 Physical geography characteristics and trophic level in the second ladder of lakes in China

2.2 TLI指數(shù)與湖盆形態(tài)(水深)的相關(guān)性分析

對位于我國地形地貌第二級階梯的19個湖泊的TLI指數(shù)與平均水深進(jìn)行曲線回歸分析,探討湖盆形態(tài)對湖泊營養(yǎng)狀態(tài)的影響.由表3和圖1可知,擬合度較高的回歸模型如下:

二次模型:

三次模型:

統(tǒng)計量對比分析:比較以上5個模型的修正R2值,由大到小的順序依次為三次模型＞二次模型＞對數(shù)模型＞冪模型＞倒數(shù)模型,由此可以判斷擬合最好的是三次模型.方差分析的F值概率均小于0.000,因此比較F值,F值由大到小的順序依次為對數(shù)模型＞冪模型＞倒數(shù)模型＞二次模型＞三次模型.通過以上判斷得出最佳模型為三次模型和對數(shù)模型.

圖1 TLI指數(shù)與平均水深的5種模型擬合Fig.1 Five fitted models of TLI and depth

2.3 TLI指數(shù)與自然地理特征的相關(guān)性分析

對位于中國地形地貌第二級階梯的19個湖泊的 TLI指數(shù)與自然地理特征進(jìn)行曲線回歸分析,結(jié)果(表4和圖2)可知,擬合度較高的回歸模型如下:

二次模型:

圖2 TLI指數(shù)與海拔的兩種模型擬合Fig.2 Two fitted models of TLI and altitude

表4 TLI指數(shù)與海拔的模型匯總和參數(shù)估計值Table 4 Model summary and coefficients of TLI and altitude

三次模型:

統(tǒng)計量對比分析:二次和三次模型的R2相等,但二次模型的F值較三次模型的大,而P值較三次模型的小,因此,二次模型的模擬效果較三次好.同理,復(fù)合、增長、指數(shù)和 Logistic模型的修正R2值及 F值相等,方差分析的 F值概率均小于0.05,擬合曲線重合為一條曲線,另外,4個模型的F值較二次模型的大,而P值較二次模型的小,因此,復(fù)合、增長、指數(shù)和Logistic模型的模擬效果最佳.

表5 TLI指數(shù)與年均溫的模型匯總和參數(shù)估計值Table 5 Model summary and coefficients of TLI and temperature

圖3 TLI指數(shù)與年均溫度的模型Fig.3 Fitted model of TLI and temperature

藻類的生長與繁殖需要適宜的溫度,因此,溫度同樣是影響湖泊營養(yǎng)狀態(tài)的重要氣候因素,通過分析TLI指數(shù)與年均溫的相關(guān)性發(fā)現(xiàn),TLI指

數(shù)與年均溫度的模擬曲線符合三次模型(表5和圖3):

而在分析 TLI指數(shù)與其他自然地理因素相關(guān)性時發(fā)現(xiàn), TLI指數(shù)與年平均降水量、年平均日照時數(shù)、年平均無霜期及年平均蒸發(fā)量無任何相關(guān)性.

2.4 Chl-a與自然地理特征的相關(guān)性分析

表6 Chl-a與年平均日照時數(shù)的模型匯總和參數(shù)估計值Table 6 Model summary and coefficients of Chl-a and ASH

圖4 Chl-a與年平均日照時數(shù)的模型Fig.4 Fitted model of Chl-a and ASH圖中二次模型與三次模型重合

湖泊水體中Chl-a的濃度可以反映藻類的生長量,從而表征湖泊的營養(yǎng)狀態(tài).對 19個湖泊的Chl-a與平均水深、海拔、年均溫、年平均降水量、年平均日照時數(shù)、年平均無霜期及年平均蒸發(fā)量進(jìn)行曲線回歸分析,發(fā)現(xiàn) Chl-a僅與平均日照時數(shù)擬合度較好,符合二次/三次模型(表 6和圖4):

加拿大學(xué)者Rawson[16]提出影響湖泊富營養(yǎng)化的主要因素有人類活動、地質(zhì)地形和緯度等因素,其中地質(zhì)地形包括湖泊的湖盆形態(tài),水深和面積等;湖泊所處的緯度決定湖泊所在區(qū)域的氣候,包括降水、溫度、日照等.由于我國地域遼闊,湖泊所處的地理位置、地形地貌特點不盡相同,導(dǎo)致富營養(yǎng)化的主要驅(qū)動因素及富營養(yǎng)化效應(yīng)存在顯著差異.通過分析位于我國地形地貌第二級階梯的19個湖泊的湖盆形態(tài)、氣候特征與湖泊中 Chl-a濃度及營養(yǎng)狀態(tài)綜合指數(shù)(TLI)之間的關(guān)系,發(fā)現(xiàn) Chl-a濃度與年平均日照時數(shù)具有較好的相關(guān)性,TLI指數(shù)與水深、海拔及年均溫的擬合模型分別符合三次/對數(shù)、二次和三次模型.說明位于我國地形地貌第二級階梯的湖泊富營養(yǎng)化主要驅(qū)動因素為水深、海拔、年均溫及年平均日照時數(shù).

2.5 中國地形地貌第二級階梯湖泊富營養(yǎng)化效應(yīng)特征

湖泊富營養(yǎng)化發(fā)生過程中會產(chǎn)生一系列的效應(yīng),包括透明度下降、溶解氧下降等,但不同的湖泊富營養(yǎng)化效應(yīng)的類型及程度存在差異,因此湖泊營養(yǎng)物基準(zhǔn)和富營養(yǎng)化控制標(biāo)準(zhǔn)的制定應(yīng)以富營養(yǎng)化效應(yīng)的區(qū)域差異性為基本依據(jù).

藻類的生長對湖泊透明度的影響是富營養(yǎng)化效應(yīng)的主要指標(biāo),通過分析2007年云南和新疆湖泊中 Chl-a與透明度之間的關(guān)系,發(fā)現(xiàn)新疆湖泊中 Chl-a與透明度呈倒數(shù)關(guān)系,而云南湖泊中Chl-a與透明度呈冪函數(shù)的關(guān)系,即透明度隨Chl-a濃度的降低而升高,模擬方程如圖5所示.

以 Chl-a/TP代表營養(yǎng)物對浮游藻類的生態(tài)效應(yīng),其生態(tài)學(xué)含義為高營養(yǎng)水平將對應(yīng)高營養(yǎng)利用效率[17].云南湖泊中Chl-a與TP線性模型的斜率為177.595而新疆的為222.758,如圖6所示,這說明,同樣位于我國第二級階梯,新疆和云南湖泊中浮游藻類利用營養(yǎng)物質(zhì)的效率存在差異.

圖5 中國地形地貌第二級階梯湖泊Chl-a與透明度之間的關(guān)系(2007年數(shù)據(jù))

圖6 中國地形地貌第二級階梯湖泊浮游藻類利用營養(yǎng)物質(zhì)的效率比較(2007年數(shù)據(jù))

3 結(jié)論

3.1 位于我國地形地貌第二級階梯的湖泊富營養(yǎng)化主要自然地理驅(qū)動因素為水深、海拔、年均溫及年平均日照時數(shù).

3.2 我國地形地貌第二級階梯湖泊,透明度與Chl-a有較強的相關(guān)性,但存在區(qū)域差異性.新疆的湖泊,Chl-a與透明度呈倒數(shù)關(guān)系,透明度隨Chl-a濃度的升高緩慢降低;云南的湖泊Chl-a與透明度呈冪函數(shù)關(guān)系,透明度隨 Chl-a濃度的升高,降低的較迅速.新疆湖泊藻類對營養(yǎng)物質(zhì)總磷的利用效率(即Chl-a/TP)高于云南湖泊.

[1] 劉永定.湖泊富營養(yǎng)化成因知多少？[EB/OL]. http://www. cenews.com.cn/xwzx/gd/qt/200910/t20091029_624023.html.2009-10-29.

[2] 金相燦,胡小貞.湖泊流域清水產(chǎn)流機制修復(fù)方法及其修復(fù)策略[J]. 中國環(huán)境科學(xué), 2010,30(3):374-379.

[3] 金相燦,劉樹坤,章宗涉,等.中國湖泊環(huán)境(第一冊) [M]. 北京:海洋出版社, 1995.

[4] 日本水污染研究會編.湖泊環(huán)境調(diào)查指南 [M]. 謝其明,王恕蓉譯.北京:中國環(huán)境科學(xué)出版社, 1989:22.

[5] 范成新,王春霞.長江中下游湖泊環(huán)境地球化學(xué)與富營養(yǎng)化[M]. 北京:科學(xué)出版社, 2007:9.

[6] 霍守亮,陳 奇,席北斗,等.湖泊營養(yǎng)物基準(zhǔn)的制定方法研究進(jìn)展 [J]. 生態(tài)環(huán)境學(xué)報, 2009,18(2):743-748.

[7] 王蘇民,竇鴻身.中國湖泊志 [M]. 北京:科學(xué)出版社, 1998.

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

[9] 王明翠,劉雪芹,張建輝.湖泊富營養(yǎng)化評價方法及分級標(biāo)準(zhǔn) [J].中國環(huán)境監(jiān)測, 2002,18(5):47-49.

[10] 金相燦,屠清瑛.湖泊富營養(yǎng)化調(diào)查規(guī)范 [M]. 北京:中國環(huán)境科學(xué)出版社, 1990.

[11] 許秋瑾,鄭丙輝,朱延忠,等.三峽水庫支流營養(yǎng)狀態(tài)評價方法[J]. 中國環(huán)境科學(xué), 2010,30(4):453-457.

[12] 金相燦,劉鴻亮,屠清瑛,等.中國湖泊富營養(yǎng)化 [M]. 北京: 中國環(huán)境科學(xué)出版社, 1990.

[13] 饒 娟,傅超英.體現(xiàn)新疆地域特色的小城鎮(zhèn)綠地系統(tǒng)規(guī)劃——以新疆達(dá)坂城鎮(zhèn)區(qū)為例 [J]. 河北農(nóng)業(yè)科學(xué), 2009,13(11):81-83.

[14] 金相燦,劉樹坤,章宗涉,等.中國湖泊環(huán)境(第二冊) [M]. 北京:海洋出版社, 1995.

[15] 石昱禎.四川省瀘沽湖地區(qū)地下水水化學(xué)特征研究 [J]. 山西建筑, 2010,36(2):146-148.

[16] Rawson D S. Some physical and chemical factors in the metabolism of lakes [J]. American Association for the Advancement of Science Public, 1939,10:9-26.

[17] 王海軍.長江中下游中小型湖泊預(yù)測湖沼學(xué)研究 [D]. 武漢:中國科學(xué)研究院水生生物研究所, 2008.

Physical geography drivers of lake eutrophication and the characteristics of eutrophication effects in the secondary topography ladder in China.

CAO Jin-ling1, XU Qi-gong1,2*, XI Bei-dou1,2, GAO Ru-tai1, JIANG Lei1, DING Jing-tao1, MAO Jing-ying1,3, SUN Xue-song1,4(1.Laboratory of Water Environmental System Engineering, Chinese Research Academy of Environmental Sciences, Beijing 100012, China；2.State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China；3.School of Environmental Science and Engineering, Southwest Jiaotong University, Chengdu 610031, China；4.School of River and Ocean Engineering, Chongqing Jiaotong University, Chongqing 400074, China). China Environmental Science, 2011,31 (11)：1849～1855

In order to investigate the key physical geography drivers of lake eutrophication in the secondary topography ladder in China, curve estimation in SPSS was used to analyze the relationships between tropic level index (TLI)/Chl-a and physical geography characteristics of lakes. relationships between TLI and depth, altitude, temperature were fit to cubic/logarithm, quadratic and cubic model, respectively. The relationship between Chl-a and annual sunshine hours (ASH) was fit to quadratic/cubic model. It indicated that the key physical geography drivers of lake eutrophication in the second geography ladder were depth, altitude, temperature and ASH. Lake eutrophication effects, the relationship between Chl-a and transparency and the slopes of linear equation for Chl-a and TP were also analyzed. Transparency was the reciprocal of Chl-a in the lakes in Xinjiang Province, and in the lakes in Yunnan Province transparency and Chl-a were in power function. The slopes of linear equation for Chl-a and TP were 177.595 and 222.758 in lakes in Yunnan and Xinjiang, respectively. This results indicated that the efficiencies of nutrients used by algae were different in the second geography ladder in China.

lake eutrophication；drivers；eutrophication effects；trophic level index；curve estimation；model fitting

X144

A

1000-6923(2011)11-1849-07

2011-03-01

國家“973”項目(2008CB418006);國家水體污染控制與治理科技重大專項(2009ZX07106-001)

* 責(zé)任作者, 研究員, xuqigong@tom.com

曹金玲(1981-),女,黑龍江齊齊哈爾人,中國環(huán)境科學(xué)研究院水環(huán)境系統(tǒng)工程研究室博士后,主要從事湖泊富營養(yǎng)化區(qū)域差異性研究.發(fā)表論文10篇.