謝林環(huán)，江濤，李銳,唐常源，，陳建耀，曹英杰,劉小波,曾紅平

(1.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275;(2.中山大學(xué)環(huán)境科學(xué)與工程學(xué)院，廣東 廣州 510275)

飛來峽水庫氮磷營養(yǎng)鹽形態(tài)及分布特征

謝林環(huán)1，江濤1，李銳2,唐常源1，2，陳建耀1，曹英杰2,劉小波1,曾紅平1

(1.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275;(2.中山大學(xué)環(huán)境科學(xué)與工程學(xué)院，廣東 廣州 510275)

以飛來峽水庫為研究對象，于2015年6月(豐水期)和12月(枯水期)對庫區(qū)內(nèi)9個采樣點進行水樣采集。通過測定各形態(tài)氮、磷濃度,結(jié)合Spearman相關(guān)性檢驗結(jié)果和氮磷比，分析水庫氮、磷營養(yǎng)鹽分布特征、結(jié)構(gòu)以及來源。結(jié)果表明：① 硝酸鹽氮是氮營養(yǎng)鹽的主要存在形態(tài)，磷營養(yǎng)鹽的主要存在形態(tài)為顆粒磷和溶解態(tài)磷；② 總氮、總磷濃度總體表現(xiàn)分別為豐水期<枯水期和豐水期>枯水期，氮濃度受生活污染影響較大，磷濃度變化具備面源污染特征；③ 氮、磷元素在豐水期具有同源性，枯水期輸入源不完全一致；④ 對比3條主要入庫河流(北江干流及滃江、連江等支流)的入庫通量，枯水期北江總氮通量最大，為1 036.0 g/s，豐水期滃江總氮通量最大，為1 128.4 g/s，豐、枯水期北江總磷通量均最大，分別為85.4和70.3 g/s；⑤ DIN/TP介于4.0～47.1，水庫屬于磷營養(yǎng)限制類型。

氮和磷形態(tài)；Spearman相關(guān)性檢驗；氮磷比；飛來峽水庫

氮、磷是構(gòu)成水體初級生產(chǎn)力重要的生源要素，也是造成湖泊、水庫富營養(yǎng)化的關(guān)鍵限制因子[1-2]，其賦存形態(tài)和含量制約著水中浮游植物的生長繁殖[3-7]。亞硝酸鹽氮、硝酸鹽氮和氨氮簡稱“三氮”，其在水體中的含量及分布特征一直是研究的重點[8-12]。磷元素同樣作為水體富營養(yǎng)化的關(guān)鍵因子，其質(zhì)量濃度的變化也會影響藻類的生長[13]。蔡龍炎等[14]通過分析32篇已公開發(fā)表的學(xué)術(shù)論文的湖泊富營養(yǎng)化有效監(jiān)測數(shù)據(jù)，結(jié)果表明我國湖泊富營養(yǎng)化空間分布特征與我國磷素和氮磷比的空間變化趨勢存在密切關(guān)系。傳統(tǒng)觀點認為氮磷比以TN/TP表示[15]，但是不同研究水體得出的比值有很大差異[16-17]。近年來，文獻[18-19]提出以DIN/TP作為營養(yǎng)限制類型的指標，當(dāng)DIN/TP從1.5～3.4變化時，湖泊、水庫生態(tài)系統(tǒng)從氮營養(yǎng)限制類型轉(zhuǎn)變?yōu)榱谞I養(yǎng)限制類型。

隨著人口的增加和社會、經(jīng)濟的發(fā)展，我國水資源供需矛盾逐年加劇，水庫供水成為緩解城市供水壓力的主要途徑[20]。飛來峽水庫作為廣東省北江干流上一座具有可調(diào)控性能的大型水庫，是珠江三角洲地區(qū)重要的后備水源，承擔(dān)著保障北江下游廣州、佛山、清遠等城市供水安全的戰(zhàn)略地位，是區(qū)域水資源安全和經(jīng)濟社會可持續(xù)發(fā)展的重要保障[21]。隨著珠江三角洲地區(qū)產(chǎn)業(yè)向東西兩翼和粵北地區(qū)轉(zhuǎn)移，飛來峽水庫流域經(jīng)濟和社會得到快速發(fā)展，人口不斷增加，庫區(qū)流域污染也隨之加大，水庫水質(zhì)明顯惡化，部分庫灣出現(xiàn)季節(jié)性藍藻、甲藻水華現(xiàn)象[22-23]，水庫水質(zhì)保護顯得尤為迫切。而目前關(guān)于飛來峽水庫氮、磷營養(yǎng)鹽的研究較少，為此，本研究擬通過豐水期和枯水期對水庫水體的采樣分析，探討水庫氮、磷營養(yǎng)鹽賦存形態(tài)、分布特征、結(jié)構(gòu)及來源，為水庫的富營養(yǎng)化控制及水資源保護提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

飛來峽水庫地處北江中游，屬典型的河道型水庫，主要由北江干流及滃江、連江等支流匯入組成。滃江于英德市江頭咀匯入北江，連江于連江口匯入北江。水庫總庫容19.04億m3，庫區(qū)北江干流全長73 km，壩址位于北江干流清遠市飛來峽鎮(zhèn)，控制流域面積34 097 km2，占北江流域面積的73%。水庫上游回水與白石窯水利樞紐水電站尾水銜接，白石窯水電站是北江干流梯級開發(fā)中的第三級，地處英德市望埠鎮(zhèn)獎家州。

圖1 采樣點分布Fig.1 Location of sampling sites

飛來峽水庫流域橫跨清遠市清城區(qū)以及英德市兩個縣級區(qū)，包括飛來峽鎮(zhèn)、黎溪鎮(zhèn)、連江口鎮(zhèn)、大站鎮(zhèn)、英紅區(qū)、英城街道、望埠鎮(zhèn)、石灰鋪鎮(zhèn)、西牛鎮(zhèn)以及水邊鎮(zhèn)等十個鎮(zhèn)級行政區(qū)劃，總面積為2 450 km2，2015年總?cè)丝诩s44萬。庫區(qū)流域內(nèi)除了少量洗礦工業(yè)外，生產(chǎn)活動基本以種植、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖等農(nóng)業(yè)生產(chǎn)為主。

1.2 樣品采集

以北江干流白石窯獎家州-飛來峽水利樞紐河段為主要研究對象，共設(shè)9個采樣點(圖1)，其中：B1(白石窯獎家洲)、B2(英德大橋)、B3(連江口)、B4(大樟庫彎)、B5(黎溪)、B6(大湖庫彎)和B7(飛來峽壩前)采樣點位于北江干流；L1、W1分別為北江干流的一級支流——連江和滃江的入河口。于2015年豐水期(6月份)和枯水期(12月份)進行2次采樣，每個點位只采集表層0.5 m深的水樣。

1.3 樣品分析與數(shù)據(jù)處理

利用SPSS22.0軟件對同一水期各形態(tài)營養(yǎng)鹽質(zhì)量濃度進行Spearman相關(guān)性分析，并進行雙側(cè)顯著性檢驗。

2 結(jié)果與討論

2.1 氮、磷營養(yǎng)鹽質(zhì)量濃度基本特征

豐水期水樣DP占TP的比例為0～99.25%，枯水期所占比例為18.33%～97.62%，豐、枯水期DP占TP比例的平均值為55.43%；豐水期水樣PP占TP的比例為0.75%～100%，枯水期所占比例為2.38%～81.67%，豐、枯水期PP占TP比例的平均值為44.57%?？菟诔^50%的水樣DP占TP的比例大于50%，豐水期超過50%的水樣PP占TP的比例大于50%，即枯水期DP為TP的主要存在形態(tài)，豐水期PP為TP的主要存在形態(tài)。

2.2 豐、枯水期氮、磷營養(yǎng)鹽質(zhì)量濃度變化

從豐水期和枯水期各形態(tài)氮、磷質(zhì)量濃度分布(圖2)可以得知：除了中游B3、B5采樣點，其他采樣點TN質(zhì)量濃度表現(xiàn)為豐水期<枯水期。豐水期、枯水期TN平均質(zhì)量濃度分別為1.59和2.03 mg/L,TN質(zhì)量濃度隨著汛期來水量增加而降低，枯水期總體水質(zhì)較差。TP總體含量較低，除了L1采樣點，其余采樣點TP質(zhì)量濃度均表現(xiàn)為豐水期>枯水期，豐水期、枯水期TP平均質(zhì)量濃度分別為0.14和0.10 mg/L。

TP和TN質(zhì)量濃度在豐、枯水期的不同特征表明水庫氮、磷元素來源不完全一致。TN在枯水期有較高的質(zhì)量濃度，說明其來源輸入較穩(wěn)定，揭示了TN受生活污染的影響可能更大。TP質(zhì)量濃度變化具有面源污染的特征，受降雨徑流污染影響較大。研究河段干流下游B5采樣點豐水期TP質(zhì)量濃度最高，其主要原因是PP質(zhì)量濃度較高，PP質(zhì)量濃度與水體中懸浮泥沙含量關(guān)系密切[25],豐水期在暴雨徑流的作用下，流域地表泥沙顆粒被徑流挾帶進入庫區(qū)水體，從而使PP質(zhì)量濃度增加。

2.3 主要入庫河流氮、磷營養(yǎng)鹽通量

外源輸入是水庫氮、磷營養(yǎng)鹽的主要來源，控制入庫污染物總量，是目前防治庫區(qū)水體富營養(yǎng)化的重要手段之一，而入庫污染物通量是外來點源、非點源污染的綜合表現(xiàn)[26]。利用B1、W1和L1等3個采樣點氮、磷質(zhì)量濃度測定值，結(jié)合北江干流白石窯水電站下泄流量及滃江長湖水庫壩下、連江高道水文站對應(yīng)月份流量數(shù)據(jù)，估算飛來峽水庫主要入庫河流TN、TP及其主要營養(yǎng)鹽存在形態(tài)的物質(zhì)通量，見表1。

豐水期TN通量呈現(xiàn)滃江>北江>連江，最高值1 128.4 g/s出現(xiàn)在滃江，占主要入庫河流總?cè)霂炝康?6.8%，這主要與滃江氮營養(yǎng)鹽濃度較高有關(guān)。滃江下游長湖水庫庫區(qū)內(nèi)有大量的水產(chǎn)養(yǎng)殖和餐飲服務(wù)，造成采樣點W1的氮營養(yǎng)鹽濃度較高?？菟赥N通量呈現(xiàn)北江>滃江>連江，最高值1 036.0 g/s出現(xiàn)在北江，占主要入庫河流總?cè)霂炝康?2.2%，這主要與白石窯下泄流量大有關(guān)。

豐水期和枯水期TP通量最高值均出現(xiàn)在北江，分別為85.4和70.3 g/s，占主要入庫河流總?cè)霂炝康?2.9%和52.2%，這主要與白石窯水電站下泄流量較大有關(guān)。

圖2 豐水期、枯水期氮、磷營養(yǎng)鹽質(zhì)量濃度分布Fig.2 The nitrogen and phosphorus concentration distribution in different water seasons

表1 主要入庫河流營養(yǎng)鹽平均通量Table 1 The average flux of nutrients in main rivers (g·s-1)

2.4 氮、磷營養(yǎng)鹽間的相關(guān)性

2.5 氮、磷營養(yǎng)鹽結(jié)構(gòu)

表2 豐水期和枯水期各項因子相關(guān)關(guān)系1)Table 2 Correlationship between nutrients in different seasons

1)**相關(guān)性非常顯著(P<0.01)；*相關(guān)性顯著(P<0.05)

表3 采樣點氮磷相對豐度DIN/TPTable 3 DIN/TP ratio of different sampling sites

3 結(jié) 論

通過對飛來峽水庫2015年6月(豐水期)和12月(枯水期)水體氮、磷營養(yǎng)鹽的采樣監(jiān)測，結(jié)合Spearman相關(guān)性檢驗和氮磷比分析，得到以下結(jié)論：

2)TN濃度總體表現(xiàn)為豐水期<枯水期；TP濃度總體表現(xiàn)為豐水期>枯水期。TN濃度主要受生活污水排放的影響，磷濃度變化情況具備面源污染特征。

3)三個主要入庫河流中，豐水期滃江TN入庫通量最大，為1 128.4 g/s，占主要入庫河流總?cè)霂炝康?6.8%；枯水期北江干流TN入庫通量大，為1 036.0 g/s，占主要入庫河流總?cè)霂炝康?2.2%；豐水期和枯水期TP入庫通量最大值均為北江干流，分別為85.4和70.3 g/s，占主要入庫河流總?cè)霂炝康?2.9%和52.2%。

4)豐水期氮、磷元素具有同源性，受生活污水排放影響的同時可能也受到農(nóng)業(yè)面源污染的影響；枯水期氮、磷元素輸入源不完全一致，揭示了氮元素受點源生活污水排放影響可能更大，磷元素可能主要受農(nóng)業(yè)面源污染影響。

5)豐水期、枯水期水樣的DIN/TP在4.0 ～47.1之間，均大于3.4，即飛來峽水庫營養(yǎng)鹽結(jié)構(gòu)為磷營養(yǎng)限制類型。

[1] 吳豐昌，金相燦，張潤宇，等.論有機氮磷在湖泊水環(huán)境中的作用和重要性[J].湖泊科學(xué),2010,22(1):1-7.

WU F C, JIN X C, ZHANG R Y, et al. Effects and significance of organic nitrogen and phosphorous in the lake aquatic environment[J]. J Lake Sci, 2010,22(1):1-7.

[2] 江濤，劉祖發(fā)，陳曉宏，等.廣東省水庫富營養(yǎng)化評價[J].湖泊科學(xué),2005(4):378-382.

JIANG T, LIU Z F, CHEN X H, et al. Assessment of reservoir eutrophication in Guangdong Province[J]. J Lake Sci, 2005,17(4):378-382.

[3] LI M T, XU K Q, WATANABE M, et al. Long-term variations in dissolved silicate, nitrogen, and phosphorus flux from the Yangtze River into the East China Sea and impacts on estuarine ecosystem[J]. Estuarine, Coastal and Shelf Science,2007,71(1/2):3-12.

[4] WIEGNER T N, SEITZINGER S P, GLIBERT P M, et al. Bioavailability of dissolved organic Nitrogen and Carbon from nine rivers in the eastern United States[J]. Aquatic Microbial Ecology,2006, 43(3):277-287.

[5] BERMAN T, BRONK D A. Dissolved organic nitrogen: a dynamic participant in aquatic ecosystems. Aquatic Microbial Ecology, 2003, 31(3):279-305.

[6] 葉琳琳，吳曉東，趙冬悅，等.崇明島河網(wǎng)浮游植物和無機、有機氮的時空分布特征[J].湖泊科學(xué),2016,28(3):528-536.

YE L L, WU X D, ZHAO D Y, et al, Temporal and spatial distributions of phytoplankton and inorganic and organic nitrogen in Chongming Island[J]. J Lake Sci, 2016(3):528-536.

[7] 楊柳，章銘，劉正文.太湖春季浮游植物群落對不同形態(tài)氮的吸收[J].湖泊科學(xué),2011,23(4):605-611.

YANG L, ZHANG M, LIU Z W. Uptake of various forms of nitrogen by phytoplankton community in spring in Lake Taihu[J]. J Lake Sci, 2011,23(4)：605-611.

[8] 楊驥，姜樹春.松花江氮轉(zhuǎn)化規(guī)律與氮污染容量研究[J].地理科學(xué),1989(3):232-241.

YANG J, JIANG S C. The study on transformation principle of nitrogen and the capacity of nitrogen pollution in the Songhua River[J]. Scientia Geographica Sinica, 1989(3):232-241.

[9] 劉敏，林莉，董磊，等.長江下游干流水體中氮的空間分布特征[J].長江科學(xué)院院報,2015,32(06):65-69.

LIU M, LIN L, DONG L, et al. Spatial distribution characteristics of nitrogen in mainstream lower Yangtze River[J]. Journal of Yangtze River Scientific Research Institute, 2015,32(06):65-69.

[10] 陳曉玲，張媛，張琍，等.豐水期鄱陽湖水體中氮、磷含量分布特征[J].湖泊科學(xué),2013,25(5):643-648.

CHEN X L, ZHANG Y, ZHANG L, et al. Distribution characteristic of nitrogen and phosphorus in Lake Poyang during high water period[J]. J Lake Sci, 2013,25(5):643-648.

[11] 王書航，王雯雯，姜霞，等.蠡湖水體氮、磷時空變化及差異性分析[J].中國環(huán)境科學(xué),2014,34(5):1268-1276.

WANG S H, WANG W W, JIANG X, et al. Spatial-temporal dynamic change of nitrogen and phosphorus and difference analysis in water body of Lihu Lake[J]. China Environmental Science, 2014,34(5):1268-1276.

[12] 朱愛萍，陳建耀，江濤，等.深圳茅洲河口營養(yǎng)鹽的時空分布特征[J].熱帶地理,2015,35(3):298-305.

ZHU A P, CHEN J Y, JIANG T, et al. Spatial and seasonal variation of nutrients transportation in the Maozhou Estuary of Shenzhen, Guangdong[J]. Tropical Geography, 2015,35(3):298-305.

[13] 豐茂武，吳云海，馮仕訓(xùn)，等.不同氮磷比對藻類生長的影響[J].生態(tài)環(huán)境,2008,17(5):1759-1763.

FENG M W, WU Y H, FENG S X, et al. Effect of different N/P ratio on algal growth[J]. Ecology and Environment, 2008,17(5): 1759-1763.

[14] 蔡龍炎，李穎，鄭子航.我國湖泊系統(tǒng)氮磷時空變化及對富營養(yǎng)化影響研究[J].地球與環(huán)境,2010,38(2):235-241.

CAI L Y, LI Y, ZHENG Z H. Temporal and spatial distribution of nitrogen and phosphorus of lake systems in China and their impact on eutrophication[J]. Earth and Environment, 2010,38(2): 235-241.

[15] DOWNING J A, McCAULEY E. The nitrogen: phosphorus relationship in lakes[J]. Limnol & Oceanogr,1992, 37(5):936-945.

[16] 李哲，郭勁松，方芳，等.三峽水庫小江回水區(qū)不同TN/TP水平下氮素形態(tài)分布和循環(huán)特點[J].湖泊科學(xué),2009(4):509-517.

LI Z, GUO J S, FANG F, et al. Potential impact of TN/TP ratio on the cycling of nirogen in Xiaojiang backwater area, three Gorges Reservior[J]. J Lake Sci, 2009(4):509-517.

[17] 聶澤宇，梁新強，邢波，等.基于氮磷比解析太湖苕溪水體營養(yǎng)現(xiàn)狀及應(yīng)對策略[J].生態(tài)學(xué)報,2012,32(1):48-55.

NIE Z Y, LIANG X Q, XING B, et al. The current water trophic status in Tiaoxi River of Taihu Lake watershed and corresponding coping strategy based on N/P ratio analysis[J]. Acta Ecologica Sinica, 2012,32(1):48-55.

[18] ANN-KRISTIN B. The use of TN:TP and DIN:TP ratios as indicators for phytoplankton nutrient limitation inoligotrophic lakes affected by N deposition[J]. Aquatic Sciences,2010,72(3):277-281.

[19] PTACNIK R, ANDERSEN T, TAMMINEN T. Performance of the redfield ratio and a family of nutrien limitation indicators as thresholds for phytoplankton Nvs. P limitation[J]. Ecosystems,2010,13(8):1201-1214.

[20] 韓博平.中國水庫生態(tài)學(xué)研究的回顧與展望[J].湖泊科學(xué)，2010,22(2):151-160.

HAN B P. Reservoir ecology and limnology in China: a retrospective comment[J]. J Lake Sci, 2010,22(2):151-160.

[21] 宋濤，黃煥坤.飛來峽水庫為珠三角壓咸補淡調(diào)度有關(guān)問題的分析[J].廣東水利水電,2007(3):39-41.

SONG T, HUANG H K. The study on scheduling of Feilaixia Reservoir for estuarine fresh water providing on Pearl River Delta[J].Guangdong Water Resources and Hydropower, 2007(3):39-41.

[22] 譚超，黃本勝，邱靜，等.飛來峽水利樞紐庫區(qū)納污能力及限制排污總量研究[J].水利水電技術(shù),2015(12):81-84,89.

TAN C, HUANG B S, QIU J, et al. Study on pollutant bearing capacity and limitation of pollutant discharge gross quantity reservoir area of Feilaixia Water Control Project[J]. Water Resources and Hydropower Engineering, 2015(12)：81-84,89.

[23] 譚超，黃本勝，洪昌宏，等.飛來峽水利樞紐庫區(qū)動態(tài)納污能力計算研究[J].水文,2014,34(5):51-56.

TAN C, HUANG B S, HONG C H, et al. Analysis and calculation of dynamic water environmental capacity in Feilaixia Reservoir[J]. Journal of China Hydrology, 2014,34(5):51-56.

[24] GB3838-2002. 地表水環(huán)境質(zhì)量標準[S]. 北京：中國環(huán)境科學(xué)出版社,2002.

[25] 陳野.長江中游水體顆粒磷與懸沙定量關(guān)系及其通量變化研究[D].武漢：長江科學(xué)院,2014.

[26] 楊尚，羅瀲蔥，翟海濤，等.大沙河水庫主要入庫河流氮磷營養(yǎng)輸入分析[J].生態(tài)科學(xué),2013(2):158-164.

YANG S, LUO L C, ZHAI H T, et al. Analysis on inputs of nitrogen and phosphorous from the main inflow rivers of Dashahe Reservoir[J]. Ecological Science, 2013(2):158-164.

[27] 李健，金中武，楊文俊.顆粒態(tài)和溶解態(tài)磷濃度計算公式的推導(dǎo)與驗證[J].長江科學(xué)院院報,2015,32(1):33-38.

LI J, JIN Z W, YANG W J. Calculation formulas for particulate and dissolved phosphorus concentrations[J]. Journal of Yangtze River Scientific Research Institute, 2015,32(1):33-38.

FormsanddistributionofnitrogenandphosphorusinFeilaixiaReservoir

XIELinhuan1,JIANGTao1,LIRui2,TANGChangyuan1,2,CHENJianyao1,CAOYingjie2,LIUXiaobo1,ZENGHongping1

(1. School of Geography and Planning, Sun Yat-sen University, Guangzhou 510275, China;(2. School of Environmental Science and Engineering, Sun Yat-sen University, Guangzhou 510275, China)

forms of nitrogen and phosphorus; spearman correlation test; N/P ratio; Feilaixia Reservoir

X524

A

0529-6579(2017)05-0112-07

10.13471/j.cnki.acta.snus.2017.05.015

2016-11-07

國家自然科學(xué)基金(41371055，41471020)；2015年廣東省水資源節(jié)約與保護專項基金

謝林環(huán)(1993年生)，女；研究方向水文與水環(huán)境；E-mail：xielh8@mail2.sysu.edu.cn；

江濤(1965年生)，女；研究方向：水文與水環(huán)境;E-mail:eesjt@mail.sysu.edu.cn