陳 丹，張 冰，曾逸凡，張海平（同濟(jì)大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 200092）

基于SWAT模型的青山湖流域氮污染時(shí)空分布特征研究

陳 丹，張 冰，曾逸凡，張海平*（同濟(jì)大學(xué)環(huán)境科學(xué)與工程學(xué)院，上海 200092）

基于SWAT模型構(gòu)建了青山湖流域的氮負(fù)荷模型，利用2005～2011年的徑流觀測(cè)資料和2011年7月15日～8月4日的水質(zhì)監(jiān)測(cè)資料對(duì)模型進(jìn)行率定和驗(yàn)證.青山湖流域氮污染的主要來(lái)源包括農(nóng)田化肥、農(nóng)村生活、畜禽養(yǎng)殖、城鎮(zhèn)點(diǎn)源及濕沉降，其中對(duì)流域出口TN濃度貢獻(xiàn)最大的來(lái)源是化肥施用，占34.3%.對(duì)流域現(xiàn)狀氮負(fù)荷的時(shí)空分布特征進(jìn)行分析，結(jié)果表明:NO3--N流失量的時(shí)空分布特征明顯，6～9月的流失量最大，占62.5%，空間上中部流失量比西北部自然林區(qū)域大，坡面流是NO3--N主要流失途徑；有機(jī)氮流失量較大的區(qū)域主要分布在中部農(nóng)田比例較高的子流域.研究結(jié)果可為流域的氮污染控制和研究提供科學(xué)依據(jù).

SWAT模型；青山湖流域；氮負(fù)荷；時(shí)空分布特征

氮是生物的基本營(yíng)養(yǎng)元素，同時(shí)也是引起水體富營(yíng)養(yǎng)化和水質(zhì)變差的根源因子.氮素污染已引起全球范圍內(nèi)的廣泛關(guān)注， 很多學(xué)者正致力于區(qū)域和全球范圍的氮循環(huán)和氮平衡問(wèn)題［1-4］，其中氮素流失是研究流域氮循環(huán)與氮平衡的基本環(huán)節(jié).資料顯示，歐洲由于農(nóng)業(yè)活動(dòng)造成的氮污染達(dá)到入海氮通量的60%［5］，在丹麥，270條河流94%的氮負(fù)荷由非點(diǎn)源污染引起［6］，荷蘭由農(nóng)業(yè)非點(diǎn)源提供的總氮（TN）也占到水環(huán)境污染總量的60%之多［7］.在我國(guó)富營(yíng)養(yǎng)化日趨嚴(yán)重的巢湖、滇池和太湖地區(qū)，由農(nóng)業(yè)活動(dòng)所造成的總氮污染貢獻(xiàn)率已達(dá)到60%～70%［5］.開展對(duì)流域氮污染的研究及流域內(nèi)氮污染特征分析，可為流域非點(diǎn)源污染治理提供理論依據(jù).

SWAT為美國(guó)農(nóng)業(yè)科學(xué)研究所（ARS）開發(fā)的一個(gè)流域尺度模型，目前在國(guó)內(nèi)外農(nóng)業(yè)非點(diǎn)源污染研究中得到了廣泛的應(yīng)用［8-11］，模型已被應(yīng)用于流域徑流模擬、流量預(yù)測(cè)和非點(diǎn)源污染評(píng)價(jià)等方面，并取得了顯著研究成果.Grizzetti等［12］利用SWAT模型研究了芬蘭Vantaanioki流域的氮素循環(huán)和轉(zhuǎn)化.Panagopoulos等［13］基于GIS利用SWAT模型模擬了希臘Arachtos中等尺度流域的氮流失情況，并指出了氮負(fù)荷關(guān)鍵源區(qū)，為降低流域氮污染提供了參考.Bouraou等［14］用SWAT來(lái)確定控制營(yíng)養(yǎng)物流失的主要過(guò)程和途徑. Daniel等［15］利用SWAT模擬氮肥施用和排水溝渠對(duì)氮流失的影響.SWAT在中國(guó)已有十多年的運(yùn)用實(shí)踐，一些學(xué)者在丹江口水庫(kù)［16］、北京密云水庫(kù)［17］、三峽庫(kù)區(qū)大寧河流域［18］、山美水庫(kù)［19］的非點(diǎn)源污染模擬均取得了較好的效果.

流域氮素流失量與降雨徑流過(guò)程、土壤性質(zhì)、植被、管理措施、農(nóng)業(yè)生產(chǎn)方式（如化肥類型、施肥方式）等因素密切相關(guān).這些影響因素的時(shí)間尺度差異較大，從年周期到一次暴雨過(guò)程不等，因此在應(yīng)用SWAT模型模擬流域氮素輸出時(shí)，應(yīng)盡可能采用高頻率的水文水質(zhì)監(jiān)測(cè)資料率定模型參數(shù)，以準(zhǔn)確反映流域各種物理、化學(xué)、生物過(guò)程及人類活動(dòng)對(duì)氮素輸出的影響，提高模擬精度.目前，國(guó)內(nèi)外已有研究多以月尺度或年尺度的實(shí)測(cè)資料率定SWAT模型.本研究嘗試在青山湖流域應(yīng)用日徑流、水質(zhì)監(jiān)測(cè)資料率定SWAT模型；以此為基礎(chǔ)，分析流域氮素污染的主要來(lái)源，探討流域氮負(fù)荷的時(shí)空分布特征，以期為青山湖流域氮污染控制提供依據(jù)，并為其他流域的氮污染研究提供參考.

1 研究區(qū)域概況

青山湖位于浙江省臨安市境內(nèi)（圖1），流域面積654km2，主要入湖河流有南苕溪、錦溪、橫溪、靈溪和雙林溪等.青山湖流域?qū)儆谥衼啛釒駶?rùn)型季風(fēng)氣候區(qū)，多年平均氣溫15.9℃，多年平均降水1613.9mm，降水日158d.降水量隨海拔高程遞減，山區(qū)大于平原，西部大于東部；降水量年際變化較大，差異在1.71～2.25:1（最大與最小年降水量比）.年內(nèi)降水量分布規(guī)律明顯，5～7月雨量充沛，其中梅雨季節(jié)常持續(xù)20d左右；7～8月間，處于副熱帶高氣壓控制下，天氣晴朗；11月至次年的3月為枯水季節(jié)，雨量稀少.目前青山湖流域面臨嚴(yán)重的氮污染問(wèn)題，對(duì)太湖氮素入湖量的貢獻(xiàn)較大，迫切需要治理.

圖1 青山湖及水文水質(zhì)監(jiān)測(cè)點(diǎn)位置Fig.1 Location of Qingshan Lake and hydrological and water quality stations

2 研究方法

2.1 SWAT模型簡(jiǎn)介

SWAT模型是一個(gè)半分布式的水文模型，主要包括水文過(guò)程子模型、土壤侵蝕子模型和污染負(fù)荷子模型［20］.在運(yùn)行模型之前，要首先進(jìn)行亞流域的劃分和最小水文響應(yīng)單元（HRU）的劃分.本研究中，整個(gè)流域劃分為46個(gè)亞流域、189個(gè)HRU，進(jìn)行分布式模擬計(jì)算.

2.2 模型率定與驗(yàn)證

圖2 橋東水文站實(shí)測(cè)和模擬月均徑流量比較Fig.2 Observed and simulated monthly runoff at Qiaodong Station

為合理模擬流域的水文過(guò)程和氮負(fù)荷特征，不僅需要掌握流域內(nèi)實(shí)際情況，而且需要通過(guò)模型率定和驗(yàn)證，對(duì)模型參數(shù)值進(jìn)行合理調(diào)整，使模型能夠真實(shí)反映流域特征.2011年7～8月對(duì)研究流域?qū)嵤┝舜笠?guī)模調(diào)研，掌握人口、畜禽養(yǎng)殖、化肥施用（量及時(shí)間）、工業(yè)污染等信息，用于模型輸入.選取2003～2004年為模型運(yùn)行的預(yù)熱期，南苕溪橋東水文站2008～2011年日觀測(cè)徑流數(shù)據(jù)用于水文模型的率定，南苕溪臨安水文站的2005～2007年日觀測(cè)徑流數(shù)據(jù)作為水文模型的驗(yàn)證數(shù)據(jù).2011年7月15日～8月4日的實(shí)測(cè)日氮負(fù)荷通量（實(shí)測(cè)河道斷面流量與濃度的乘積）作為水質(zhì)模型的率定數(shù)據(jù).流域水文水質(zhì)監(jiān)測(cè)點(diǎn)位置見圖1.

采用確定性系數(shù)R2、Nash-Suttcliffe系數(shù)Ens作為模型精度評(píng)價(jià)指標(biāo).一般認(rèn)為， R2≥0.6、Ens≥0.5時(shí)，模型的模擬結(jié)果是可以接受的.評(píng)價(jià)結(jié)果表明，率定期月徑流模擬的R2、Ens分別為0.87和0.78，日徑流模擬的R2、Ens值分別為0.82 和0.71；驗(yàn)證期月徑流模擬的R2、Ens分別為0.86 和0.69，日徑流模擬的R2、Ens值分別為0.77和0.63. 圖2～圖3分別為率定期月尺度及日尺度橋東水文站實(shí)測(cè)與模擬徑流量的比較.對(duì)長(zhǎng)達(dá)6年的歷史日徑流觀測(cè)資料的高精度重演表明，所構(gòu)建的模型較為準(zhǔn)確地反映了流域水文特征，可合理描述流域的地表徑流過(guò)程.對(duì)氮素負(fù)荷通量的模擬結(jié)果（圖4）顯示，各監(jiān)測(cè)點(diǎn)硝態(tài)氮（NO3-N）通量模擬的R2、Ens值分別在0.81、0.55以上，TN通量模擬的R2、Ens值分別在0.78、0.53以上，表明模型能較好模擬各亞流域的氮素產(chǎn)生和運(yùn)移過(guò)程.

圖3 橋東水文站實(shí)測(cè)和模擬日徑流量比較Fig.3 Observed and simulated daily runoff at Qiaodong Station

圖4 各監(jiān)測(cè)點(diǎn)N-N率定結(jié)果Fig.4 Calibration results for N-N at all stations

3 流域氮污染來(lái)源分析

青山湖流域現(xiàn)狀氮素來(lái)源主要包括農(nóng)田化肥、農(nóng)村生活、畜禽養(yǎng)殖、城鎮(zhèn)點(diǎn)源（城鎮(zhèn)生活污水和工業(yè)廢水）及濕沉降.各項(xiàng)污染源對(duì)流域現(xiàn)狀氮素的年貢獻(xiàn)量分別為:農(nóng)田化肥施用4173t、農(nóng)村生活與畜禽養(yǎng)殖為1140t、城鎮(zhèn)點(diǎn)源量243t以及來(lái)源于濕沉降的870t.

表1 流域氮源輸入量及輸出貢獻(xiàn)率Table 1 Input amount of nitrogen source and output contribution rate in Qingshan Lake Basin

應(yīng)用青山湖流域SWAT模型，分析了各氮源對(duì)流域氮素流失的貢獻(xiàn)，見表1.有機(jī)氮流失多來(lái)源于植物的殘留，因此在有機(jī)氮輸出中，流域背景值（主要來(lái)自植物殘留）貢獻(xiàn)最高，為24t.農(nóng)田化肥施用造成的N-N輸出量最高，達(dá)341t，其次是來(lái)自濕沉降，第三是農(nóng)村生活與畜禽養(yǎng)殖.與其他氮源相比，濕沉降的N-N更容易流失，流失率為24.3%，這是由于部分NN隨降雨直接沉降于溪流、湖泊等水體表面，或者被坡面徑流攜帶進(jìn)入溪流. N-N的輸出量絕大部分來(lái)源于城鎮(zhèn)點(diǎn)源，為147t.從對(duì)流域出口TN輸出的貢獻(xiàn)上分析，農(nóng)田化肥的貢獻(xiàn)最高，對(duì)TN濃度貢獻(xiàn)為1.07mg/L，占34.3%，表明削減農(nóng)田化肥用量是流域氮素污染治理的一個(gè)主要方向.城鎮(zhèn)點(diǎn)源的輸入量雖然不高，但由于其氮素沒(méi)有經(jīng)過(guò)土壤系統(tǒng)削減，直接排入水體，削減率很低，造成在流域輸出量很高，表明控制點(diǎn)源量是流域氮素污染治理的一種有效方案.

4 流域氮素流失時(shí)間分布特征

4.1 年際變化特征

降雨是氮素流失的驅(qū)動(dòng)力，對(duì)氮素流失影響顯著.表2為流域水文及氮素流失年際變化趨勢(shì). 2005～2011年降雨量1250～1450mm.徑流深年均770mm，徑流系數(shù)為0.57. NN流失量年均1170t，呈現(xiàn)出逐年減少趨勢(shì)，有機(jī)氮（ON）流失量年均426t，最高發(fā)生在2006年，達(dá)616t.從年尺度上分析，NN流失量與降雨量流失年相關(guān)系數(shù)僅為0.40，有機(jī)氮流失量與降雨量相關(guān)系數(shù)為-0.46，這表明氮素流失量與降雨量之間的分布無(wú)明顯規(guī)律，兩者的相關(guān)性不大.

表2 流域水文及氮素流失年際變化Table 2 Yearly variation of precipitation-runoff and nitrogen loss in Qingshan Lake Basin

4.2 月際變化特征

圖5 流域氮素流失的月變化趨勢(shì)Fig.5 Seasonal trend of nitrogen loss in the basin

2005～2011年流域降雨的月分布特征明顯，降雨量主要集中在6～8月份，約占年降雨總量的43.6%.圖5為氮素流失的月變化趨勢(shì). N-N流失月變化特征明顯，集中在6～9月份，與降雨分布基本一致，占總流失量的62.5%.通過(guò)分析模型結(jié)果發(fā)現(xiàn)，影響NN月變化特征的因素除降雨外，施肥時(shí)機(jī)也十分重要，7～9月份農(nóng)田施肥量較大，當(dāng)遇到強(qiáng)度較大的降雨時(shí)，氮素很容易從土壤系統(tǒng)中大量流失.因此，優(yōu)化施肥時(shí)機(jī)和方式可有效降低流域氮素流失量.

5 流域氮素流失空間分布特征

除降雨徑流外，氮素流失受下墊面的影響也較大.圖6為流域內(nèi)氮素流失量的空間分布，顯示各區(qū)域具有較大的差異性.2005～2011年，青山湖流域單位N-N流失量為6～50kg/hm2，流失量較大區(qū)域主要集中在流域中部農(nóng)田比例較高的子流域.流域西北部雖然坡度較大，但由于土地利用以自然林為主，氮素輸入量低， N-N流失大多低于10kg/hm2. N-N為3種不同形式徑流攜帶，坡面流流失量約占總流失量的47.2%，并且與后者空間分布基本一致. NN側(cè)向流流失量占總流失量的34.6%，在流域中部較高，流失量為15～20kg/hm2，其余區(qū)域流失量較低，為2～10kg/hm2. NN基流流失量占總流失量的18.2%，在0～10kg/hm2之間，空間分布差異不大，與基流深分布一致，這可能是由于N-N容易在深層土壤中殘留累積［21］，使區(qū)域間濃度差異變小.ON流失量較大的區(qū)域分布在流域中部農(nóng)田比例較高的子流域，年均流失負(fù)荷達(dá)到24.7kg/hm2，其余流域的ON流失量多在0～17kg/hm2之間.

圖6 流域氮素流失單位負(fù)荷量的空間分布Fig.6 Spatial distribution of the unit load of nitrogen loss in Qingshan Lake Basin

6 結(jié)論

6.1 構(gòu)建適合青山湖流域氮素動(dòng)態(tài)模擬的SWAT模型，利用實(shí)測(cè)水文（2005～2011年日徑流數(shù)據(jù)）及氮素負(fù)荷通量（2011年7月15日～8月4日水文水質(zhì)同步監(jiān)測(cè)數(shù)據(jù)）資料對(duì)模型進(jìn)行率定和驗(yàn)證.月徑流量的R2和Ens值分別達(dá)到0.87和0.78，日徑流量的和Ens值分別為0.82和0.71，8個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)TN和NO3--N日負(fù)荷通量的R2值均大于0.81、Ens值大于0.53，表明模擬精度較高，所建立的模型能合理描述流域氮素輸出負(fù)荷對(duì)地表徑流過(guò)程的動(dòng)態(tài)響應(yīng).

6.2 在ON輸出中，流域植物殘留導(dǎo)致的流域背景值貢獻(xiàn)最大，為24t；農(nóng)田化肥施用是NO3--N輸出的最大貢獻(xiàn)者，達(dá)341t； NO3--N的輸出量絕大部分來(lái)源于城鎮(zhèn)點(diǎn)源，為147t.對(duì)流域出口TN濃度貢獻(xiàn)最大的氮素來(lái)源是農(nóng)田化肥，其貢獻(xiàn)為1.07mg/L，占34.3%，表明控制農(nóng)田化肥施用是流域氮素污染治理的一個(gè)主要方向.

［1］ShresthaJ，Niklaus P A， Pasquale N， et al. Flood pulses control soilnitrogencyclingin a dynamic river floodplain ［J］. Geoderma，2013，228:14-24.

［2］Bouwman L， Goldewijk K K， Van Der Hoek K W， et al. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900-2050period ［J］. Proceedings of The National Academy of Sciences of The United States of America， 2013，110（52）:20882-20887.

［3］Rabouille S， Van de Waal D B， Matthijs H C P， et al. Nitrogen fixation and respiratory electron transport in the cyanobacteriumCyanothece under different light/dark cycles ［J］. FemsMicrobiology Ecology， 2014，87（3）:630-638.

［4］ElbertW， Weber B， Burrows S， et al. Contribution of cryptogamic covers to the global cycles of carbon and nitrogen ［J］. Nature Geoscience， 2012，5（7）:459-462.

［5］韓 博，周麗麗，范昊明.蒲河流域氮污染負(fù)荷模擬及時(shí)空分布［J］. 中國(guó)水土保持科學(xué)， 2013，11（2）:90-96.

［6］Kronvang B， Graesboll P， Larsen S E. Diffuse nutrient losses in Denmark ［J］. Water science and Technology， 1996，33（4）:81-88.

［7］Boers P C M. Nutrient emissions from agriculture in Netherlands: causes and remedies ［J］. Water Science and Technology， 1996， 33（4）:183-190.

［8］Wang X， Shang S， Yang W， et al. Simulation of land use-soil interactive effects on water and sediment yields at watershedscale［J］. Ecological Engineering， 2010，36:328-344.

［9］Lam Q D， Schmalz B， Fohrer N. Modelling point anddiffuse source pollution of nitrate in a rural lowlandcatchment using the SWAT model ［J］. Agricultural Water Management， 2010，97:317-325.

［10］Park M J， Park J Y， Shin H J， et al. Projection of futureclimate change impacts on nonpoint source pollutionloads for a forest dominant dam watershed by reflectingfuture vegetation canopy in a soil and water assessmenttoolmodel ［J］. Water Science and Technology， 2010，61:1975-1986.

［11］范麗麗，沈珍瑤，劉瑞民.基于SWAT模型的大寧河流域非點(diǎn)源污染空間特征研究 ［J］. 水土保持通報(bào)， 2008，28（4）:133-137.

［12］Grizzetti B， Bouraoui F， Granlund K. Modeling diffuse emission and retention of nutrients in the Vantaanjoki watershed （Finland） using the SWAT model ［J］. Ecological Modeling， 2003，169:25-38.

［13］Panagopoulos Y， Makropoulos C， Baltas E. SWAT parameterization for the identification of critical diffuse pollution source areas under data limitations ［J］. Ecological Modelling， 2011，222（19）:3500-3512.

［14］Bouraoui F， Grizzetti B. An integrated modeling framework to estimate the fate of nutrients: Application to the Loire （France） ［J］. Ecological Modeling， 2008，212（3）:450-459.

［15］Daniel N M， Prasanna H G， Jeffrey G A， et al. Modeling the impact of nitrogen fertilizer application and tile drain configuration on nitrate leaching using SWAT ［J］. Agricultural Water Management， 2013，130:36-43.

［16］王 曉，郝芳華，張 璇.丹江口水庫(kù)流域非點(diǎn)源污染的最佳管理措施優(yōu)選 ［J］. 中國(guó)環(huán)境科學(xué)， 2013，33（7）:1335-1343.

［17］王曉燕，秦福來(lái)，歐 洋，等.基于SWAT模型的流域非點(diǎn)源污染模擬-以密云水庫(kù)北部流域?yàn)槔?［J］. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2008，27（3）:1098-1105.

［18］Nash J E， Sutcliffe J V. River flow forecasting throughconceptual models. Part I-A discussion of principles ［J］. Journal of Hydrology， 1970，10:282-290.

［19］劉梅冰，陳冬平，陳興偉.山美水庫(kù)流域水量水質(zhì)模擬的SWAT與CE-QUAL-W2聯(lián)合模型 ［J］. 應(yīng)用生態(tài)學(xué)報(bào)， 2013，24（12）: 3574-3580.

［20］唐 寅.運(yùn)用SWAT模型研究小流域氣候及土地利用變化的水文響應(yīng) ［D］. 北京:北京林業(yè)大學(xué)， 2011.

［21］王激清，韓寶文，劉社平.施氮量和耕作方式對(duì)春玉米產(chǎn)量和土體硝態(tài)氮累積的影響 ［J］. 干旱地區(qū)農(nóng)業(yè)研究， 2011，29（2）: 129-135.

致謝：本論文的前期調(diào)研由同濟(jì)大學(xué)環(huán)境科學(xué)專業(yè)11屆及10屆環(huán)科師生協(xié)助完成，在此表示感謝.

Analysis of temporal and spatial characteristics of nitrogen pollution in Qingshan Lake basin based on SWAT model.

CHEN Dan， ZHANG Bing， ZENG Yi-fan， ZHANG Hai-ping*（College of Environmental Science and Engineering，Tongji University， Shanghai 200092， China）. China Environmental Science， 2015，35（4）：1216～1222

A nitrogen load model has been constructed for Qingshan Lake Basin based on SWAT model. The model was calibrated and verified using the runoff measurements during 2005to 2011and water quality measurements during July 15to August 4， 2011. Main nitrogen sources include fertilizer use， domestic sewage， livestock， urban point source and wet deposition， among which fertilizer use is the biggest contributor to TN concentration at the outlet of the basin， amounting to 34.3% of the total. Temporal and spatial characteristics of the nitrogen load distribution have been analyzed using the model. The results show that seasonal characteristics of NO3--N loss are distinct， with 62.5% of the yearly total loss concentrated on the period of June to September； highest in the middle region and lowest in the northern west region （covered mostly by natural forest）； overland flow is the main loss route for NO3--N； the highest loss of organic nitrogen is found in the middle region of the basin where farmland is the dominant land-use type. The results of the present study can contribute to the management and study of nitrogen pollution in basins.

SWAT model；Qingshan Lake Basin；nitrogen load；temporal and spatial characteristics

X524，X171

A

1000-6923（2015）04-1216-07

陳 丹（1990-），女，河南新安人，碩士，主要從事生態(tài)水力學(xué)研究.

2014-09-20

國(guó)家水體污染控制與治理科技重大專項(xiàng)課題（2008ZX07101-006-07）

* 責(zé)任作者， 教授， hpzhang@#edu.cn