江西省中小型湖泊水體營養狀態評價及其驅動因子研究

江立文，陳楊，林暾，王亞，施俊譯，阮磊磊，曹金香，李梅

華東交通大學土木建筑學院，江西 南昌 330013

湖泊是地球上重要的水資源之一，不僅滿足人類生產生活需求，還具有調節氣候、保持生態平衡和防洪灌溉等重要功能。隨著人口增長、工農業發展和城市化進程加快，過量營養元素和有機物被排入湖中，致湖泊富營養化問題日益突出。江西省位于長江中下游南岸，地處 24°29'N ～30°04'N，113°34'E～118°28'E，面積 16.69 萬 km2［1］。全省 1 km2以上湖泊有55個，總面積為3882.7 km2，主要分布在江西省西北部，除鄱陽湖為特大型湖泊外，90%的湖泊面積小于 50 km2［2-3］。

目前，江西省湖泊的研究對象大多集中在鄱陽湖區［4-7］，關于中小型湖泊水體的營養狀態卻鮮有報道。筆者利用綜合營養狀態指數(TLI)評價江西省內25個中小型湖泊水體的營養狀態，采用主成分分析法研究了湖泊水體富營養化的驅動因子，以期為江西省湖泊水體富營養狀態的防治提供理論依據。

1 數據來源

1.1 湖泊監測概況

于2010年初開始，選取江西省境內分布于長江干流水系、修河水系、環鄱陽湖水系的25個中小型湖泊，開展了為期一年的水質監測工作，監測湖泊總面積近800 km2。江西省中小型湖泊水質監測數據如表1所示。

表1 江西省中小型湖泊監測概況Table 1 General situation of small or medium size lakes in Jiangxi Province

1.2 布點及分析

湖泊監測點位設置采用網格布點法，25個湖泊共設176個取樣點，監測頻率為每季一次。采樣方法為水深小于5 m時，水面下0.5 m處取樣;水深大于5 m時，混合取樣。采集的水樣經預處理后，置于暗處4℃恒溫保存，并于當日進行水質分析。監測指標的選取及分析方法以《水和廢水監測分析方法》(4 版)［8］為依據。

使用Excel2003和PASW Statistics 18對所得監測數據進行統計分析。

2 研究方法

2.1 水體營養狀態評價

目前，我國湖泊水體營養狀態評價方法有Carlson營養狀態指數(TSI)法、修正Carlson營養狀態指數(TSIM)法、綜合營養狀態指數(TLI)法、營養度指數法和評分指數法［9］。筆者選取TLI法對湖泊水體的營養狀態進行評價。以透明度(SD)、高錳酸鹽指數(CODMn)、總氮(TN)濃度、總磷(TP)濃度和葉綠素a(Chl-a)濃度對所監測湖泊進行營養狀態評價。計算公式如下:

式中，TLI為綜合營養狀態指數;Wj為第j種指標參數的營養狀態指數的相關權重;TLI(j)為第j種指標參數的營養狀態指數;rij為第j種指標參數與基準指標參數Chl-a濃度的相關系數，中國湖泊部分指標參數與Chl-a濃度相關系數rij及其rij2值如表2所示;m為評價指標參數的個數。

表2 中國湖泊部分指標參數與Chl-a濃度的相關系數［10］Table 2 Correlation coefficients between chlorophyll-a and index parameters of some lakes in China

由表2得出:

TLI(SD)=10(5.118-1.94ln SD)

TLI(CODMn)=10(0.109+2.66ln CODMn)

TLI(TN)=10(5.453+1.694ln TN)

TLI(TP)=10(9.436+1.624ln TP)

TLI(Chl-a)=10(2.5+1.086ln Chl-a)

式中，SD，m;Chl-a 濃度，mg/m3;CODMn，TN 濃度和TP濃度，mg/L。

湖泊營養狀態分級:TLI＜30為貧營養;30≤TLI≤50為中營養;50＜TLI≤60為輕度富營養;60＜TLI≤70為中度富營養;TLI＞70為重度富營養［10］。

2.2 驅動因子分析

主成分分析方法是通過降維，將多個變量轉換成較少的幾個綜合變量［11］，在最大限度保留原始數據信息的同時，客觀地確定各變量的權重［12］，該方法簡單有效，已廣泛應用于水環境質量評價［13-20］。

3 結果與討論

3.1 湖泊營養狀態的評價

采用綜合營養狀態指數法計算了25個湖泊的TLI，如圖1所示。

圖1 江西省中小型湖泊綜合營養狀態指數Fig.1 Comprehensive eutrophication state index of small or medium size lakes in Jiangxi Province

從圖1可以看出，25個調查湖泊的TLI均大于30，已達到中營養狀態以上。八里湖(55.82)、陳家湖(51.05)、大 湖 (51.90)、芳 湖 (50.27)、南 北 湖(57.74)、太泊湖(53.57)、瑤崗湖(52.43)和瑤湖(56.44)等8個湖泊的TLI超過了50，水體呈現富營養化狀態，其中，芳湖最低，南北湖最高;8個湖泊的TLI均未超過60，為輕度富營養，其占監測湖泊總數的32%。其余的17個湖泊均為中營養湖泊，其中3個湖的 TLI較低，分別是赤湖(34.13)、軍山湖(34.31)和珠湖(37.73);而內青嵐湖(49.98)和新妙湖(49.50)的TLI相對較大，處在富營養化臨界狀態。

3.2 富營養化驅動因子的分析

以富營養化的8個湖泊(八里湖、陳家湖、大湖、芳湖、南北湖、太泊湖、瑤崗湖和瑤湖)為考察對象，以水深(DW)，水溫(TW)，SD，pH，溶解氧(DO)濃度，CODMn，TN 濃度，硝酸鹽(NO3-N)濃度，TP 濃度，溶解性磷酸鹽(DP)濃度和Chl-a濃度等11項監測指標為分析變量，研究了湖泊水體富營養化的驅動因子，8個富營養化湖泊監測指標參數的年均值統計結果如表3所示。

運用PASW Statistics 18對變量進行方差最大旋轉，得出主成分在各變量上的載荷分布，如圖2所示;湖泊水體富營養化驅動因子的主成分分析結果，如表4所示。

表3 8個富營養化湖泊監測指標參數的年均值統計結果Table 3 Statistical results of annual average monitoring index about 8 eutrophication lakes

圖2 主成分因子在各變量上的載荷分布Fig.2 The loading distribution of principal components on variables

表4 主成分法分析結果Table 4 The results of principal component analysis(PCA)

由表4可知，四類主成分的特征值都大于1，且其累計貢獻率達到了87.194%，可以表征8個富營養化湖泊11項監測指標的大部分信息。

從圖2可以看出，利用主成分分析法，可將江西省中小型湖泊水體富營養化的驅動因子分成四類:第一類驅動因子(第一主成分)是以TN(0.95)和NO3-N(0.96)為代表的含氮營養鹽指標，以及pH(0.56)指標，在該類型中的pH反映湖泊的酸堿度，經計算，所研究湖泊的pH平均為7.95，符合富營養化湖泊水質呈弱堿性的特征［9］;第二類驅動因子(第二主成分)是以TP(0.93)和DP(0.82)為代表的含磷營養鹽指標;第三類驅動因子(第三主成分)是與湖泊自身特性有關的自然特征指標，如 DW(0.91)，SD(0.80)和TW(0.75);第四類驅動因子(第四主成分)是包括Chl-a和CODMn在內的生物和有機物特性指標。

目前，在湖泊水體富營養化的機理研究中，對適合水生動植物生長的溫度、有助于藻類等浮游生物生長的水流速度和用以提供水生生物生長的氮磷等營養鹽濃度［21］三個影響因素為大多數研究者所認同。從表4可以看出，第一主成分的貢獻率最大，為24.864%，第二主成分和第三主成分的貢獻率與第一主成分相差不大，說明水體中氮磷營養鹽濃度及其自然特征變化情況在驅動江西省中小型湖泊水體富營養化方面所起的作用相當，且較大;而第四主成分的貢獻率最小，為15.759%，由于Chl-a是表征水體中浮游植物生物量和生產力的重要指標，該項貢獻率反映了CODMn與水體富營養化間的關系，說明江西省中小型湖泊水體富營養化受有機污染物的影響不大。分析表明，江西中小型湖泊水體中氮磷營養鹽濃度、溫度等自然特征變化情況與湖泊水體富營養化的關系十分密切。

4 結論

江西省25個中小型湖泊水體的營養化狀態達到中度營養狀態，其中8個湖泊的水體已達到輕度富營養化狀態，分別是八里湖、陳家湖、大湖、芳湖、南北湖、太泊湖、瑤崗湖和瑤湖。

中小型湖泊水體富營養化狀態的驅動因子分成四類:第一類是以TN、NO3-N為代表的含氮營養鹽以及pH等指標;第二類是以TP、DP為代表的含磷營養鹽指標;第三類是以DW、TW等為代表的與湖泊自身特性有關的自然特征指標;第四類是以Chl-a和CODMn為代表的生物和有機物特性指標。結果表明，水體中氮磷營養鹽濃度、溫度和水深等自然特征指標是江西省中小型湖泊水體富營養化主要驅動因子。

［1］江西省統計局，國家統計局江西調查總隊.江西統計年鑒(2010)［M］.北京:中國統計出版社，2010.

［2］楊桂山，馬榮華，張路，等.中國湖泊現狀及面臨的重大問題與保護策略［J］.湖泊科學，2010，22(6):799-810.

［3］江西省水文局.江西水系［M］.武漢:長江出版社，2006:12.

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

［5］王毛蘭，周文斌.鄱陽湖水體無機氮時空分布特征研究［J］.人民長江，2010，41(1):88-91.

［6］向蘇林，周文斌，聶發輝.鄱陽湖沉積物間隙水中磷的含量及其分布特征［J］.湖泊科學，2011，23(6):868-872.

［7］向蘇林，周文斌.鄱陽湖沉積物中磷的賦存形態及分布特征［J］.湖泊科學，2010，22(5):649-654.

［8］國家環境保護總局.水和廢水監測分析方法［M］.4版.北京:中國環境科學出版社，2002.

［9］王明翠，劉雪芹，張建輝.湖泊富營養化評價方法及分級標準［J］.中國環境監測，2002，18(5):47-49.

［10］金相燦.中國湖泊環境:第一冊［M］.北京:海洋出版社，1995.

［11］盧紋岱.SPSS統計分析［M］.4版.北京:電子工業出版社，2010:472.

［12］劉總堂，李春海，章鋼婭.運用主成分分析法研究云南湖庫水體中重金屬分布［J］.環境科學研究，2010，23(4):459-466.

［13］史小紅，李暢游，賈克力.烏梁素海污染現狀及驅動因子分子［J］.環境科學與技術，2007，30(4):37-39.

［14］伊元榮，海米提·依米提，王濤，等.主成分分析法在城市河流水質評價中的應用［J］.干旱區研究，2008，25(4):497-501.

［15］余進祥，劉婭菲，鐘曉蘭，等.鄱陽湖水體富營養化評價方法及主導因子研究［J］.江西農業學報，2009，21(4):125-128.

［16］PERONA E U，BONILLA I，MATEO P.Spatial and temporal changes in water quality in a Spanish river［J］.Sci Total Environ，1999，24(1/2/3):75-90.

［17］KARIM B，MARHABA T F.Using principal component analysis to monitor spatial and temporal changes in water quality［J］.J Hazard Materi，2003，100(1/2/3):179-195.

［18］SOLIDOROA C，PASTRES R，COSSARINI G，et al.Seasonal and spatial variability of water quality parameters in the lagoon of Venice［J］.J Marine Systems，2004，51(1/2/3/4):7-18.

［19］BU H M，TAN X，LI S Y，etal.Temporal and spatialvariations of water quality in the Jinshui River of the South Qinling Mts，China［J］.Ecotoxicology Environmental Safety，2010，73(5):907-913.

［20］de CEBALLOS B S O，KONING A，de OLIVERA J F，et al.Dam reservoir eutrophication:a simplified technique for a fastdiagnosis of an environmental degradation［J］.Water Res，1998，32(11):3477-3483.

［21］尹真真.國內外水體富營養化機理研究歷史與進展［J］.微量元素與健康研究，2006，23(3):46-47.○