高田田，謝暉，萬能勝，熊竹陽，胡正華，賴錫軍

（1.南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，應(yīng)用氣象學(xué)院，南京 210044；2.中國科學(xué)院南京地理與湖泊研究所，中國科學(xué)院流域地理學(xué)重點實驗室，南京 210008；3.安徽省巢湖管理局湖泊生態(tài)環(huán)境研究院，合肥 238000）

氮磷是湖泊生態(tài)系統(tǒng)所必需的生源要素，但由于人類活動加劇導(dǎo)致的過量氮磷入湖極易誘發(fā)水體富營養(yǎng)化和藻華暴發(fā)[1]。隨著點源污染得到有效管控，面源污染開始成為氮磷入湖的重要驅(qū)動過程與來源，闡明入湖流域面源污染機(jī)理與調(diào)控機(jī)制是湖泊污染治理的關(guān)鍵瓶頸問題，而面源污染模擬則是解決這一瓶頸問題的重要環(huán)節(jié)。面源污染的隨機(jī)性、分散性、廣泛性等復(fù)雜特征加大了污染輸出核算的難度，僅用監(jiān)測手段難以完全反映污染輸出在流域尺度上的時空異質(zhì)性?，F(xiàn)階段常用的面源污染模型有流域水文模型和統(tǒng)計型模型。其中常用的流域水文模型包括SWAT 模型、HSPF 模型、AnnAGNPS 模型、HYPE 模型等[2]，模型從最初的流域水循環(huán)過程模擬，到逐漸耦合泥沙和營養(yǎng)鹽等模塊，可對流域子系統(tǒng)內(nèi)污染物遷移轉(zhuǎn)化的物理過程進(jìn)行定量化描述，但此類模型通常包含較多概念性參數(shù)，需要翔實的實測數(shù)據(jù)對其進(jìn)行參數(shù)率定以及驗證，參數(shù)識別復(fù)雜[3]。統(tǒng)計型模型的代表是輸出系數(shù)模型[4]，其方法成熟且應(yīng)用簡單，因此被國內(nèi)外學(xué)者廣泛應(yīng)用。輸出系數(shù)模型是基于不同來源污染物的輸出強度對流域內(nèi)輸出總量的多年均化進(jìn)行估算，考慮到降雨作為面源污染的直接驅(qū)動力以及地形對面源污染傳輸路徑的重要影響，有研究在輸出系數(shù)模型中引入了降雨和地形因子[5-8]，這為流域面源污染的時空差異模擬提供了更為客觀的方法。

巢湖是長江中下游大型淺水湖泊，作為我國水污染防治的重點水體，其富營養(yǎng)化問題一直備受關(guān)注，已有研究表明面源污染是巢湖入湖氮磷污染的主要貢獻(xiàn)源[9-10]。巢湖入湖流域多為農(nóng)村流域，農(nóng)村流域的面源污染來源廣泛、成因復(fù)雜、核算難度大。2018—2020 年巢湖流域水雨情變化顯著，2019 年是巢湖流域的極枯年，2020 年巢湖流域則經(jīng)歷了百年一遇的特大洪水，這為分析流域面源污染特征在典型水雨情條件下的差異提供了難得的現(xiàn)實場景。烔煬河作為直接入湖的重要河流之一，其流域自然地理和社會經(jīng)濟(jì)特征在巢湖農(nóng)村流域中極具代表性。鑒于此，本文以烔煬河流域為研究區(qū)域，模擬近3 a 不同典型水雨情條件下面源污染氮磷輸出，分析其時空變化特征，解析流域分區(qū)、分時的污染組成結(jié)構(gòu)與來源，以期為削減入湖面源污染和保障巢湖水環(huán)境健康提供科學(xué)支撐。

1 材料與方法

1.1 研究區(qū)概況

烔煬河流域是巢湖北岸一支入湖農(nóng)村小流域（圖1a），流域總面積為89.3 km2。流域內(nèi)的烔河（自西北向東南）和煬河（自西向東）兩條主干河流在入湖口上游交匯為烔煬河后流入巢湖。根據(jù)巢湖小流域劃分及編碼技術(shù)規(guī)定，將烔煬河流域劃分為24 個子流域（圖1b）。流域總體地勢為西北高、東南低，上游以丘陵崗地為主，中下游的鎮(zhèn)區(qū)-2 號子流域為高密度人口聚集地，下游烔煬河干流-1 號子流域為平原圩區(qū)，其地形地貌在巢湖流域中具有代表性。流域內(nèi)土地利用類型以林地和農(nóng)田為主，種植的主要農(nóng)作物有水稻、油菜和小麥等，建設(shè)用地主要位于烔煬鎮(zhèn)鎮(zhèn)區(qū)。烔煬河流域?qū)俦眮啛釒夂騾^(qū)，年均氣溫約16 ℃，多年平均降水量約1 030 mm，雨量充沛，汛期時降雨事件頻發(fā)，流域面源污染問題突出[11]。根據(jù)前期調(diào)研，該流域面源氮磷污染輸出主要來自于種植業(yè)化肥流失、農(nóng)村生活污水和畜禽養(yǎng)殖等，具有巢湖農(nóng)村流域面源污染的統(tǒng)一特征，因此在該流域研究面源氮磷的污染特征和組成結(jié)構(gòu)，對治理農(nóng)村流域面源氮磷污染和防治巢湖富營養(yǎng)化具有重要意義。

圖1 烔煬河流域地理位置及子流域和水系Figure 1 Location，sub-catchments，and rivers in the Tongyang River catchment

1.2 數(shù)據(jù)來源和模擬方法

1.2.1 數(shù)據(jù)來源

（1）氣象水文。烔煬河流域內(nèi)建有烔煬雨量站，具有連續(xù)小時尺度的降雨量數(shù)據(jù)。筆者前期針對烔煬河流域已搭建了HSPF 水文模型，烔河干流處設(shè)有一處矩形薄壁堰，基于每日河流流量監(jiān)測值已進(jìn)行模型參數(shù)率定（模擬值和監(jiān)測值的納什系數(shù)為0.7）[11]，因此可獲得流域內(nèi)各河流的模擬流量數(shù)據(jù)，用于計算該流域入湖污染通量。

（2）下墊面。烔煬河流域數(shù)字高程模型（Digital Elevation Model，DEM）數(shù)據(jù)（圖2a）和土地利用類型圖（圖2b）由當(dāng)?shù)毓芾聿块T提供，精度均為1∶10 000，其中土地利用類型劃分為水田、旱田、林地、草地、城鎮(zhèn)和水域。

圖2 烔煬河流域數(shù)字高程模型和土地利用類型Figure 2 DEM and land use type of the Tongyang River catchment

（3）污染源。烔煬河流域內(nèi)點源污染來自于烔煬鎮(zhèn)污水處理廠，尾水直接排入烔河干流（3 號河流）尾端，通過每日的尾水水質(zhì)和水量監(jiān)測數(shù)據(jù)可計算流域內(nèi)點源氮磷污染輸出量。流域面源污染源主要考慮農(nóng)村生活污水、畜禽養(yǎng)殖和土地利用。農(nóng)村人口數(shù)量根據(jù)自然村位置以及常住人口調(diào)研獲取，根據(jù)現(xiàn)場走訪調(diào)研和烔煬鎮(zhèn)污染源調(diào)查數(shù)據(jù)得知，研究區(qū)內(nèi)共有13家畜禽養(yǎng)殖場，養(yǎng)殖種類為羊、肉雞、蛋雞和鵝。

（4）水質(zhì)監(jiān)測。烔煬河水質(zhì)自動監(jiān)測站根據(jù)《地表水自動監(jiān)測技術(shù)規(guī)范（試行）》（HJ 915—2017）建設(shè)并運維，高頻觀測數(shù)據(jù)被生態(tài)環(huán)境部門用于考核烔煬河流域入湖斷面水質(zhì)達(dá)標(biāo)情況。該自動監(jiān)測站位于入湖河段中下游（117°38′10.26″E，31°39′40.20″N），距離入湖口約700 m，可對總氮、氨氮、總磷、化學(xué)需氧量、pH 值、水溫等水質(zhì)指標(biāo)進(jìn)行實時監(jiān)測，監(jiān)測頻率為每4 h 一次。根據(jù)入湖河流的流量和總氮、總磷濃度即可計算烔煬河流域氮、磷污染物總量，由于4 h頻次的氮、磷監(jiān)測濃度和日均模擬流量數(shù)據(jù)存在尺度差異，因此對每日濃度數(shù)據(jù)取平均值，再結(jié)合日均流量數(shù)據(jù)進(jìn)行入河污染負(fù)荷的計算。

1.2.2 面源氮磷污染輸出模擬

（1）改進(jìn)的輸出系數(shù)模型。輸出系數(shù)模型是用于模擬流域面源污染輸出的常用方法，其基本原理是基于不同污染源的源強系數(shù)，結(jié)合土地利用面積和污染源單位數(shù)量等進(jìn)行計算。傳統(tǒng)的輸出系數(shù)模型未能將降雨和地形的影響考慮在內(nèi)，但這二者是決定流域面源污染輸出的關(guān)鍵因素。本研究采用改進(jìn)的輸出系數(shù)模型，對烔煬河流域面源氮磷輸出進(jìn)行模擬，具體公式如下：

式中：Li為第i種污染物的輸出量，kg；α為降雨因子；β為地形因子；Ei,j為第i種污染物第j個污染源的輸出系數(shù)，kg·hm-2·a-1、kg·人-1·a-1或 kg·只-1·a-1；Aj為第j種污染源的土地利用面積、人口數(shù)量或畜禽養(yǎng)殖數(shù)量等。

（2）降雨因子取值。本研究采用降雨因子表征降雨年際差異對流域面源污染輸出的影響。前期的相關(guān)研究表明面源污染氮磷負(fù)荷與降雨量存在一定正相關(guān)性[12]，則第i種污染物第t年的降雨因子（αi,t）可表示為：

式中：Li,t為第i種污染物第t年的面源污染負(fù)荷為第i種污染物的多年平均面源污染負(fù)荷，kg；rt為第t年的降雨量為多年平均降雨量，mm；fi表示第i種污染物面源污染負(fù)荷與降雨量之間的函數(shù)關(guān)系。本研究根據(jù)入湖口斷面處氮磷監(jiān)測濃度與流量計算流域氮磷污染負(fù)荷總量，再減去點源負(fù)荷即可計算烔煬河流域面源氮磷負(fù)荷，并采用2018—2020 年烔煬雨量站月降雨量與月負(fù)荷進(jìn)行擬合，方程如下：

式中：LTN和LTP分別為總氮和總磷的月負(fù)荷，kg；r為月降雨量，mm。

（3）地形因子取值。本研究采用地形因子表征子流域坡度空間差異對流域面源污染輸出的影響。前期研究表明，流域面源污染輸出與坡度存在冪函數(shù)關(guān)系[12]，因此地形因子可根據(jù)下式進(jìn)行計算：

式中：Lj為第j個子流域的面源污染輸出量，kg；θj為第j個子流域的坡度，（°）為研究區(qū)平均坡度，（°）；c和d為常量。根據(jù)相關(guān)文獻(xiàn)d取值為0.610 4[13]，基于烔煬河流域的DEM 數(shù)據(jù)計算出該流域平均坡度為3.34°，由此可根據(jù)各子流域坡度計算地形因子。

（4）輸出系數(shù)取值。本研究涉及到的污染源包括生活污水、畜禽養(yǎng)殖以及土地利用下的地表徑流，為保證輸出系數(shù)的取值具有相對的準(zhǔn)確性，本文選擇巢湖流域以及長江中下游流域的文獻(xiàn)作為相關(guān)系數(shù)的取值依據(jù)。畜禽養(yǎng)殖的氮磷污染輸出系數(shù)由飼養(yǎng)周期、糞尿排泄系數(shù)、糞尿流失系數(shù)以及污染物平均質(zhì)量比[14-16]綜合計算獲得，生活污水和各類土地利用的氮磷輸出系數(shù)參考相關(guān)文獻(xiàn)[17-23]取值。烔煬河流域污染源的具體輸出系數(shù)匯總?cè)绫?所示。

表1 烔煬河流域污染源輸出系數(shù)Table 1 Export coefficients of different pollution sources in the Tongyang River catchment

2 結(jié)果與討論

2.1 面源氮磷輸出時間變化

傳統(tǒng)輸出系數(shù)模型未能考慮降雨年際差異對流域面源污染輸出的影響，估算結(jié)果僅為多年平均值。本研究針對烔煬河流域近年來雨情年際大幅變化這一特征，在輸出系數(shù)模型中引入降雨因子對該流域2018—2020 年的面源污染氮磷輸出進(jìn)行重新模擬，結(jié)果如表2 所示。近3 a 的烔煬河流域面源污染總氮和總磷輸出平均值分別為（261.33±173.87）t·a-1和（11.03±9.98）t·a-1，年際差異顯著，這與降雨量的強變化幅度密切相關(guān)。巢湖流域2019 年為大旱年，烔煬河流域全年降雨量不足700 mm，遠(yuǎn)低于多年平均水平。而2020 年巢湖流域在梅雨季節(jié)暴雨事件頻發(fā)且強度大，累積降雨量居近30 a 之首，形成大面積洪澇災(zāi)害[24]；烔煬河流域在此期間降雨量高達(dá)832.3 mm，全年降雨量為1 424.7 mm，約是2019年年降雨量的兩倍。由于面源污染與降雨量的非線性關(guān)系[公式（3）和公式（4）]，2020 年面源污染總氮輸出總量高達(dá)453.53 t·a-1，為 2019 年的 3.9 倍，2020 年面源總磷輸出總量為 22.20 t·a-1，是 2019 年的 7.4 倍，降雨是造成氮磷輸出顯著差異的主要原因，國內(nèi)其他相關(guān)研究也有類似報道。王宏等[25]對2012—2013 年間沱江農(nóng)業(yè)小流域的氮磷排放特征展開研究，結(jié)果表明，2012 年7—9 月的高降雨量造成該年份氮磷流失量遠(yuǎn)高于2013 年。張林等[26]通過研究蘭陵溪小流域在雨季的氮磷輸出發(fā)現(xiàn)，豐水期降雨是土壤中氮磷污染輸出的主要驅(qū)動力，期間由降雨徑流攜帶流失的總氮和總磷負(fù)荷分別占雨季氮、磷流失量的88%和90%。降雨量對面源總磷輸出的影響相對于總氮更為明顯，磷相比于氮更容易吸附在顆粒物上，降雨-徑流和土壤侵蝕過程增強了磷的排放與傳輸[27-29]，因此總磷輸出對2020 年高降雨量的響應(yīng)比總氮更敏感。王趙飛等[30]的研究表明，烔煬河流域中農(nóng)田顆粒態(tài)磷的流失較高，而流域農(nóng)田面積占總面積的50%以上，因此應(yīng)密切關(guān)注強降雨期間的農(nóng)田總磷流失。

表2 烔煬河流域面源氮磷污染源輸出總量Table 2 Export of nitrogen and phosphorus pollution in the Tongyang River catchment

2018—2020 年流域面源污染入河量由流域入湖污染總量與點源尾水入河量相減獲得，流域面源污染入河量與面源污染輸出量之比即為入河系數(shù)[31]。2018—2020 年流域面源總氮和總磷的入河系數(shù)分別為0.08~0.15 和0.19~0.26，與國內(nèi)外相關(guān)研究中的入河系數(shù)取值基本處于同一水平[32-34]，這也進(jìn)一步驗證了本研究采用的改進(jìn)輸出系數(shù)方法及系數(shù)取值在巢湖典型農(nóng)村流域的適用性。

2.2 面源氮磷輸出空間變化

本研究考慮了生活污水、畜禽養(yǎng)殖和土地利用為流域面源污染源，各污染源的氮磷輸出在空間上的差異明顯。如圖3a 和圖3b 所示，鎮(zhèn)區(qū)-2 號子流域和煬河源頭-24 號子流域具有生活污水總氮和總磷的高污染輸出。鎮(zhèn)區(qū)-2號子流域內(nèi)統(tǒng)計常住人口為4 500余人，居子流域之首，煬河源頭-24 號子流域位居其次。但由于煬河源頭-24 號子流域地處上游，坡度較高（坡度因子為1.21），鎮(zhèn)區(qū)-2 號子流域處于流域中下游，地勢平緩（坡度因子為0.67），地形因子的引入使得煬河源頭-24 號子流域的生活污水總氮和總磷年均輸出最高，分別達(dá)到677.87 kg·a-1和113.73 kg·a-1。烔河左支-6 號子流域內(nèi)雖自然村常住人口僅400 余人，但該片區(qū)內(nèi)包含烔煬中學(xué)和烔煬敬老院兩處高密度人口聚集點，據(jù)統(tǒng)計烔煬中學(xué)現(xiàn)有師生1 500余人，敬老院常住人口110 余人，這使得烔河左支-6號子流域的生活污水氮磷輸出僅次于鎮(zhèn)區(qū)-2 號子流域和煬河源頭-24號子流域。

烔煬河流域內(nèi)畜禽養(yǎng)殖帶來的面源總氮和總磷年均輸出為 7.73 t·a-1和 1.87 t·a-1，其中超過 64.31%的總氮和超過55.25%的總磷輸出來自于烔河源頭-11 號子流域（圖3c 和圖3d）。烔煬河流域現(xiàn)有13 家畜禽養(yǎng)殖場，總體數(shù)量上以家禽為主，烔河源頭-11號子流域內(nèi)有一家擁有600 頭羊的養(yǎng)殖廠。由于羊的糞尿排泄系數(shù)遠(yuǎn)高于家禽，使得該分區(qū)內(nèi)的畜禽養(yǎng)殖氮磷污染非常突出。在煬河左干-17 號子流域中有3 家畜禽養(yǎng)殖場，其氮磷污染輸出僅次于烔河源頭-11號子流域。

烔煬河子流域在2018—2020 年的土地利用氮磷污染輸出受到地形因子、土地利用面積及其輸出系數(shù)的綜合影響。如圖3e和圖3f所示，煬河源頭-24號子流域的土地利用氮磷污染輸出年均值最高（總氮32.08 t·a-1、總磷 1.17 t·a-1）。該分區(qū)總面積為 10.83 km2，為所有子流域之首，且其中農(nóng)田面積高達(dá)4.91 km2，農(nóng)田的氮磷輸出系數(shù)受化肥施用影響而相對較高。此外，該分區(qū)地處上游低山丘陵區(qū)，地形因子也高達(dá)1.21。綜合以上因素，煬河源頭-24 號子流域的土地利用氮磷污染輸出最高。此外，烔河源頭-11 號子流域和烔河左支-10 號子流域由于同樣原因也造成氮磷污染輸出較高。地形坡度對氮磷污染輸出的影響可通過烔河左支-6 號子流域的模擬結(jié)果體現(xiàn)，該分區(qū)面積第二，農(nóng)田面積最大，但由于子流域平均坡度小，地形因子僅為0.78，使得土地利用氮磷污染輸出低于10號和11號子流域。

圖3 烔煬河子流域不同污染源的氮磷輸出空間分布Figure 3 Spatial distribution of the nitrogen and phosphorus export from different sources in the Tongyang River catchment

綜合生活污水、畜禽養(yǎng)殖和土地利用的模擬結(jié)果（圖4），面源總氮污染年均輸出最大的子流域為煬河源頭-24 號子流域[（32.76±21.80）t·a-1]，總磷污染年均輸出最大的子流域為烔河源頭-11 號子流域[（1.81±1.64）t·a-1]，二者分別是煬河和烔河的源頭子流域。以單位面積污染強度來看，煬河源頭-24 號子流域的輸出強度略有下降，而烔河源頭-11 號子流域具有最高的總氮年均污染強度（4.48 t·km-2·a-1）和總磷年均污染強度（0.29 t·km-2·a-1）。對總氮而言，煬河左支-20號子流域和煬河左支-22號子流域雖面積較小，但污染強度僅次于烔河源頭-11 號子流域，主要是由于其包含農(nóng)田、水域以及羊糞尿等高輸出系數(shù)的污染來源。對總磷而言，煬河左干-17 號子流域的排放強度高達(dá)0.22 t·km-2·a-1，僅次于烔河源頭-11號子流域。

圖4 烔煬河子流域面源氮磷污染輸出量與強度的空間分布Figure 4 Spatial distribution of nitrogen and phosphorus export and export intensity in the Tongyang River catchment

2.3 面源氮磷輸出來源分析

2018—2020 年烔煬河流域面源總氮和總磷的年均輸出量分別為 261.33 t·a-1和 11.03 t·a-1，其組成結(jié)構(gòu)（圖5）中農(nóng)田占據(jù)了絕大部分比例，特別是總氮，農(nóng)田的輸出量占比高達(dá)81%。烔煬河流域現(xiàn)階段土地開發(fā)強度較低，農(nóng)田依然超過總面積50%，且存在尿素、復(fù)合肥等化肥的過量施用以及利用率低等問題。據(jù)當(dāng)?shù)剞r(nóng)戶問卷調(diào)查統(tǒng)計，烔煬河流域的化肥折純量施用強度為307 kg·hm-2，高于發(fā)達(dá)國家為防止化肥施用造成環(huán)境危害設(shè)定的上限值（225 kg·hm-2）[35]。這種過量施肥的現(xiàn)象在巢湖農(nóng)村流域普遍存在。前人研究表明，巢湖流域的種植業(yè)施肥強度居全國前列，是長江中下游化肥流失的重污染區(qū)域[36]。在這樣的背景之下，合肥市政府明確了“巢湖流域一級保護(hù)區(qū)水稻生產(chǎn)化肥2021 年底實現(xiàn)零增長、2022 年底實現(xiàn)零使用”的目標(biāo)任務(wù)。這一舉措勢必將從源頭上控制營養(yǎng)鹽向土壤中的直接輸入，從而減少面源氮磷污染輸出。水域是烔煬河流域氮磷輸出貢獻(xiàn)的第二來源，坑塘為該流域水域的主要組成部分?？犹潦黔h(huán)巢湖地帶的特色水適應(yīng)性景觀單元，在烔煬河流域中分布廣而密[37]，據(jù)統(tǒng)計坑塘約占該流域總面積的9%，且其中大部分被開發(fā)利用為養(yǎng)殖池塘。在現(xiàn)階段高投入、高產(chǎn)出的養(yǎng)殖模式下，餌料和飼料通常被過量投放，從而極易造成氮磷的流失。在長江流域重點水域十年禁漁的政策背景下，以內(nèi)陸?zhàn)B殖為導(dǎo)向的坑塘開發(fā)利用模式所帶來的氮磷污染新形勢需要被重點關(guān)注。畜禽養(yǎng)殖對面源總磷輸出的貢獻(xiàn)達(dá)到了17%，是僅次于農(nóng)田的第二大污染源。據(jù)調(diào)研，烔煬河流域現(xiàn)存養(yǎng)殖場雖規(guī)模較小，但糞污仍以直接還田這種“自循環(huán)”處理方式為主，應(yīng)加強好氧堆肥和厭氧治沼等資源化和無害化處理，從而減少畜禽養(yǎng)殖氮磷污染輸出。

圖5 烔煬河流域面源污染氮磷輸出組成結(jié)構(gòu)Figure 5 Composition of nonpoint source nitrogen and phosphorus export in the Tongyang River catchment

子流域中污染源的組成結(jié)構(gòu)（圖6）呈現(xiàn)空間異質(zhì)性。從面源總氮污染輸出來看，超過90%的子流域中土地利用的貢獻(xiàn)高于75%。畜禽養(yǎng)殖的總氮污染貢獻(xiàn)比例在烔河源頭-11 號子流域和煬河左干-17號子流域中較為突出，分別達(dá)到了17.84% 和13.96%。生活污水總氮污染輸出貢獻(xiàn)僅在鎮(zhèn)區(qū)-2 號子流域中高于8%。各子流域中面源總磷污染輸出的組成比例與總氮具有一定差異，總磷污染強度最高的烔河源頭-11 號子流域中畜禽養(yǎng)殖的貢獻(xiàn)達(dá)到了56.96%，在煬河左干-17 號子流域中的貢獻(xiàn)超過了50%，在煬河左支-16 號子流域也達(dá)到了40.13%。生活污水總磷污染貢獻(xiàn)同樣表現(xiàn)為鎮(zhèn)區(qū)-2 號子流域最大，高達(dá)28.14%，因此應(yīng)切實推進(jìn)該區(qū)域人口聚集區(qū)的污水管網(wǎng)改造，加強生活污水的收集納管，整治雨污混接、錯接問題。

圖6 烔煬河子流域面源氮磷污染輸出組成結(jié)構(gòu)Figure 6 Composition of nonpoint source nitrogen and phosphorus export in each sub-catchment of the Tongyang River catchment

3 結(jié)論

（1）采用改進(jìn)的輸出系數(shù)模型模擬了烔煬河流域2018—2020 年面源污染氮磷輸出，降雨因子和地形因子的引入細(xì)化了模擬結(jié)果，結(jié)果表明近3 a 烔煬河流域面源總氮和總磷污染輸出年均值分別為（261.33±173.87）t·a-1和（11.03±9.98）t·a-1，2019 年和2020 年水雨情條件的大幅變化造成污染輸出的年際差異顯著。農(nóng)田是研究區(qū)面源污染的首要來源，其對總氮和總磷的污染輸出貢獻(xiàn)分別高達(dá)81%和52%。

（2）子流域?qū)用嫦旅嬖吹孜廴据敵隹偭?、強度和組成結(jié)構(gòu)呈現(xiàn)空間異質(zhì)性。源頭子流域的面源氮磷污染輸出總量居首，其中烔河源頭子流域具有最高的畜禽養(yǎng)殖氮磷污染輸出，煬河源頭子流域具有最高的生活污水及土地利用氮磷污染輸出。

（3）烔煬河流域土地利用中農(nóng)田占比高，農(nóng)田化肥過量施用以及利用率低的問題，是巢湖其他農(nóng)村流域普遍存在的現(xiàn)象。切實推廣“化肥使用零增長、零使用”等政策將從源頭上減少污染輸出。農(nóng)業(yè)面源污染防治更應(yīng)把握重點分區(qū)，基于污染來源因地制宜地配置管控措施，針對烔河源頭-11 號和煬河左干-17號子流域應(yīng)加強畜禽糞污的無害化和資源化處理，針對鎮(zhèn)區(qū)和烔煬中學(xué)等人口密集區(qū)應(yīng)加強生活污水收集治理。