關(guān)鍵詞：農(nóng)業(yè)面源污染；負(fù)荷評估；水環(huán)境質(zhì)量；響應(yīng)關(guān)系

中圖分類號：X52 文獻(xiàn)標(biāo)志碼：A 文章編號：1672-2043（2024）12-2817-11 doi：10.11654/jaes.2024-0881

當(dāng)前農(nóng)業(yè)面源污染已經(jīng)取代點(diǎn)源污染成為水環(huán)境污染的最主要來源[1]，統(tǒng)計(jì)數(shù)據(jù)表明2022年我國農(nóng)業(yè)源化學(xué)需氧量（COD）、總氮和總磷的排放量分別為1785.7萬、174.4萬t和27.7萬t，分別占總排放量的68.8%、55.0% 和80.2%[2]。過多氮素、磷素等營養(yǎng)物質(zhì)進(jìn)入水體環(huán)境會(huì)誘發(fā)水體富營養(yǎng)化、破壞水環(huán)境生態(tài)，影響水資源安全利用等[3]。目前農(nóng)業(yè)面源污染管控已成為我國生態(tài)文明建設(shè)的重點(diǎn)方向，《“十四五”土壤、地下水和農(nóng)村生態(tài)環(huán)境保護(hù)規(guī)劃》指出推進(jìn)長江經(jīng)濟(jì)帶等重點(diǎn)區(qū)域農(nóng)業(yè)面源污染監(jiān)督與治理是農(nóng)業(yè)綠色發(fā)展的重要內(nèi)容。

國內(nèi)外學(xué)者圍繞農(nóng)業(yè)面源污染成因、傳輸機(jī)制、污染負(fù)荷評估、環(huán)境效應(yīng)以及防治策略等開展了一系列研究。其中，農(nóng)業(yè)面源污染與水環(huán)境質(zhì)量響應(yīng)關(guān)系的研究對于厘清農(nóng)業(yè)面源污染的環(huán)境效應(yīng)、制定科學(xué)的管控對策具有重要的理論和現(xiàn)實(shí)意義。有研究發(fā)現(xiàn)雨季洱海流域耕地、園地流失的總氮、總磷、氨氮均與周邊水質(zhì)呈正相關(guān)[4]。高福等[5]發(fā)現(xiàn)COD、總磷、氨氮等污染物與地表水水質(zhì)指數(shù)呈顯著的正相關(guān)關(guān)系，錢曉雍[6]分析了不同時(shí)間上海淀山湖區(qū)域農(nóng)業(yè)面源污染氮磷排放量與淀山湖氮磷濃度、存量的相關(guān)性，發(fā)現(xiàn)農(nóng)業(yè)面源氮磷排放量與水質(zhì)濃度僅在個(gè)別季節(jié)具有相關(guān)性，上游來水才是影響淀山湖水質(zhì)的主要因素。宋曉明等[7]研究發(fā)現(xiàn)湖南46個(gè)以農(nóng)業(yè)生產(chǎn)活動(dòng)為主的縣（市、區(qū)）單位國土面積農(nóng)業(yè)面源污染負(fù)荷與水質(zhì)類別顯著正相關(guān)。綜上可知，農(nóng)業(yè)面源污染對區(qū)域性流域水質(zhì)的影響程度因地而異，可能受到諸多環(huán)境因素的干擾。寧波是長江經(jīng)濟(jì)帶重要的工業(yè)城市，也屬典型平原河網(wǎng)地區(qū)，具有水系發(fā)達(dá)、往復(fù)流嚴(yán)重、水體自凈能力較差等特點(diǎn)。2023年寧波市地表水國控?cái)嗝嫠|(zhì)優(yōu)良率同比下降18.2個(gè)百分點(diǎn)，省控及以上斷面水質(zhì)優(yōu)良率同比下降7.4 個(gè)百分點(diǎn)[8]，水環(huán)境質(zhì)量管控形勢不容樂觀。第二次污染源普查表明長江經(jīng)濟(jì)帶農(nóng)業(yè)面源污染的COD、總氮、總磷排放量分別占該區(qū)域水污染排放總量的45.1%、48.0% 和66.6%[9]，但農(nóng)業(yè)面源污染與水質(zhì)間是否存在直接的影響尚未明確。為厘清寧波市農(nóng)業(yè)面源對地表水環(huán)境質(zhì)量的影響程度，本研究基于農(nóng)業(yè)面源污染負(fù)荷核算和水質(zhì)監(jiān)測數(shù)據(jù)，系統(tǒng)研究姚江、奉化江、甬江一級支流控制單元農(nóng)業(yè)面源污染負(fù)荷與水質(zhì)的空間分異特征，旨在解析兩者之間的響應(yīng)關(guān)系及主要污染來源，明確根據(jù)單位國土面積農(nóng)業(yè)面源污染負(fù)荷預(yù)測水環(huán)境質(zhì)量的適用范圍，以期為長江中下游經(jīng)濟(jì)帶工業(yè)發(fā)達(dá)地區(qū)農(nóng)業(yè)面源污染監(jiān)督管理提供理論支撐。

1 材料與方法

1.1 流域匯水單元?jiǎng)澐?/p>

寧波市地處浙江省東北部，地勢西南高，東北低。甬江貫穿寧波市，為東海獨(dú)流入海河流，是浙江省八大水系之一。上游源頭有姚江、奉化江兩支，在寧波市三江口匯合成甬江，并于寧波鎮(zhèn)?？诹魅霒|海。根據(jù)三江一級支流匯水范圍和行政邊界將53個(gè)鎮(zhèn)街劃定為32個(gè)控制單元，總面積1 513.35 km²。其中部分鎮(zhèn)街位于匯水范圍內(nèi)但無一級支流出水口，因此根據(jù)地理位置和水流方向?qū)⑵洳⑷胂掠捂?zhèn)街控制單元，具體包括：（1）戚家山街道并入小港街道控制單元；（2）貴駟街道、九龍湖街道、駱駝街道、招寶山街道并入蛟川街道控制單元；（3）洞橋鎮(zhèn)、古林鎮(zhèn)并入石碶街道控制單元；（4）集士港鎮(zhèn)并入高橋鎮(zhèn)控制單元；（5）外灘街道并入甬江街道控制單元；（6）三七市鎮(zhèn)并入河姆渡鎮(zhèn)控制單元；（7）白云街道、南門街道、西門街道、鼓樓街道、江廈街道、月湖街道并入文教街道；（8）百丈街道、中河街道、東勝街道、東柳街道、東郊街道并入白鶴街道控制單元。

1.2 農(nóng)業(yè)面源污染負(fù)荷計(jì)算

參考《第二次全國污染源普查產(chǎn)排污系數(shù)手冊》[10]，使用排污系數(shù)法評估農(nóng)業(yè)（包括種植業(yè)、畜禽和水產(chǎn)養(yǎng)殖）及農(nóng)村生活污水的COD、氨氮、總磷排放量。各農(nóng)業(yè)源污染物排放量計(jì)算公式如下：

1.3 水質(zhì)樣品采集與測定

分別在寧波市甬江、姚江和奉化江一級支流下游匯入干流處附近各設(shè)置一個(gè)監(jiān)測斷面，共163個(gè)監(jiān)測斷面，分布于32個(gè)控制單元。監(jiān)測頻次為每月1次。樣品采集參考《水和廢水監(jiān)測分析方法》，各監(jiān)測斷面采集0.5 m 深度處水樣500 mL，置于聚乙烯塑料瓶中，帶回實(shí)驗(yàn)室進(jìn)行測定，監(jiān)測指標(biāo)為高錳酸鹽指數(shù)、氨氮及總磷。高錳酸鹽指數(shù)采用堿性高錳酸鉀法（GB 17378.4—2007）測定，氨氮采用納氏試劑分光光度法（HJ 535—2009）測定，總磷采用鉬酸銨分光光度法（GB/T 11893—1989）測定。

1.4 基礎(chǔ)數(shù)據(jù)來源收集

以鎮(zhèn)（街道）為單位，收集控制單元內(nèi)各鎮(zhèn)（街道）的基礎(chǔ)信息，具體包括：各鎮(zhèn)（街道）國土面積、行政區(qū)劃矢量范圍、主要作物播種面積、園地面積、家禽存欄量及出欄量、生豬出欄量、牛存欄量、水產(chǎn)品產(chǎn)量（陸地人工養(yǎng)殖）、農(nóng)村常住人口與農(nóng)村生活污水處理模式（納管率），上述數(shù)據(jù)主要來源于生態(tài)環(huán)境、農(nóng)業(yè)農(nóng)村及水利相關(guān)部門。產(chǎn)排污系數(shù)數(shù)據(jù)來源于《第二次全國污染源普查產(chǎn)排污系數(shù)手冊》[10]。

1.5 數(shù)據(jù)處理與分析

由于當(dāng)前地表水水質(zhì)類別呈非線性變化，無法有效表征其與農(nóng)業(yè)面源污染負(fù)荷之間的關(guān)系，因此本研究采用污染物等差遞增的標(biāo)準(zhǔn)將水質(zhì)分為Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ和Ⅵ共6類[7]。將各控制單元斷面12個(gè)月污染物濃度進(jìn)行平均計(jì)算得出年度平均污染物濃度，再將其與各控制單元農(nóng)業(yè)面源COD、氨氮、總磷排放負(fù)荷進(jìn)行相關(guān)性分析，并利用P 值檢驗(yàn)法驗(yàn)證兩者是否顯著相關(guān)。利用Excel 2021對數(shù)據(jù)進(jìn)行初步統(tǒng)計(jì)和處理，利用ArcGIS 10.2進(jìn)行負(fù)荷與水質(zhì)空間分異表征，利用Origin 2018制圖。

2 結(jié)果與討論

2.1 農(nóng)業(yè)面源污染負(fù)荷分析

2023年農(nóng)村生活污水、種植業(yè)、水產(chǎn)養(yǎng)殖業(yè)和畜禽養(yǎng)殖業(yè)對寧波市32個(gè)控制單元農(nóng)業(yè)面源COD、氨氮、總磷污染物負(fù)荷的貢獻(xiàn)率如圖1所示。COD排放量為1 135.3 t，農(nóng)村生活污水排放貢獻(xiàn)率最高（57.3%），其次是畜禽養(yǎng)殖業(yè)（33.4%），水產(chǎn)養(yǎng)殖業(yè)貢獻(xiàn)率為9.3%。氨氮排放量為266.8 t，主要來源于農(nóng)村生活污水、種植業(yè)和畜禽養(yǎng)殖業(yè)，貢獻(xiàn)率依次為61.0%、26.1%、11.6%?？偭着欧帕繛?18.9 t，主要來自種植業(yè)，貢獻(xiàn)率達(dá)76.7%，其次是農(nóng)村生活污水，貢獻(xiàn)率為16.4%?？傮w來說，農(nóng)村生活污水和畜禽養(yǎng)殖業(yè)主要排放COD與氨氮，種植業(yè)主要排放總磷，氨氮次之。

由圖2可知，寧波市32個(gè)控制單元農(nóng)業(yè)面源污染負(fù)荷差異較大。高橋鎮(zhèn)、石碶街道以及河姆渡鎮(zhèn)控制單元的COD 污染負(fù)荷較高，分別為329.3、177.1、130.1 t。位于鄞州區(qū)的梅墟街道等7個(gè)控制單元、位于江北區(qū)的莊橋街道等4個(gè)控制單元以及位于海曙區(qū)的段塘街道等2個(gè)控制單元的COD污染負(fù)荷均為0。高橋鎮(zhèn)、橫街鎮(zhèn)以及石碶街道控制單元的氨氮排放量較高，均大于35 t；高橋鎮(zhèn)、石碶街道以及姜山鎮(zhèn)控制單元的總磷排放量較高，均大于15 t。段塘街道、新明街道、明樓街道以及福明街道控制單元的COD、氨氮、總磷排放量最低，均小于0.01 t，由于這些控制單元的城鎮(zhèn)化程度較高，均沒有規(guī)模種養(yǎng)，且農(nóng)村生活污水均已納管處理，因此污染物排放量較小。總體來說，高橋鎮(zhèn)、姜山鎮(zhèn)、橫街鎮(zhèn)、石碶街道、蛟川街道、莊市街道、方橋街道以及河姆渡鎮(zhèn)共8個(gè)控制單元的農(nóng)業(yè)面源污染較嚴(yán)重，無論COD、氨氮還是總磷排放量均位于前列。寧波市西南部是農(nóng)業(yè)面源污染較集中的區(qū)域。

圖3表示單位國土面積農(nóng)業(yè)面源污染負(fù)荷，與圖2相比部分控制單元污染負(fù)荷格局發(fā)生了變化。就單位國土面積COD 污染負(fù)荷而言，高橋鎮(zhèn)控制單元仍位列第一，但蛟川街道、橫街鎮(zhèn)及河姆渡鎮(zhèn)控制單元排名均有明顯下降。氨氮污染負(fù)荷方面，原本氨氮污染負(fù)荷排名較高的陸埠鎮(zhèn)、蛟川街道控制單元單位國土面積污染負(fù)荷排名下降顯著，其余控制單元的排名變化不明顯?？偭孜廴締挝粐撩娣e污染負(fù)荷最大的為姜山鎮(zhèn)，其次為方橋街道。相對污染總負(fù)荷，陸埠鎮(zhèn)、蛟川街道、橫街鎮(zhèn)等土地面積較大的控制單元單位國土面積污染負(fù)荷的排名呈現(xiàn)明顯下降趨勢。

隨著長江中下游地區(qū)工業(yè)和城鎮(zhèn)點(diǎn)源污染問題得到有效控制，面源污染逐漸從次要矛盾上升為主要矛盾，成為水生態(tài)環(huán)境保護(hù)工作的突出難點(diǎn)[11]。寧波市具有農(nóng)村人口稠密、生豬養(yǎng)殖密度較高、復(fù)種指數(shù)高等特點(diǎn)，因此厘清農(nóng)業(yè)面源污染來源問題是污染防治的前提。沈婷[12]研究發(fā)現(xiàn)生活污染源和畜禽養(yǎng)殖源是長江下游面源污染排放的主要來源，兩者污染物總排放占比分別為39.5%和36.5%。本研究同樣發(fā)現(xiàn)農(nóng)村生活污水和畜禽養(yǎng)殖業(yè)是農(nóng)業(yè)面源污染的主要來源。目前雖然寧波市大部分村鎮(zhèn)已通過納管進(jìn)入污水處理廠的方式處理農(nóng)村生活污水，但仍有部分行政村以農(nóng)村污水處理設(shè)施或就地消納的方式處理農(nóng)村生活污水，加之寧波市常住人口眾多，因此對農(nóng)業(yè)面源污染貢獻(xiàn)較大。另外，由于不同控制單元之間國土面積差異較大，為了剔除國土面積對總量的制約，采用單位國土面積農(nóng)業(yè)面源污染負(fù)荷可以更有效地表征地區(qū)的污染負(fù)荷強(qiáng)度[7]。本研究發(fā)現(xiàn)高橋鎮(zhèn)、石碶街道以及姜山鎮(zhèn)控制單元農(nóng)業(yè)面源污染排放量和負(fù)荷強(qiáng)度均處前列，因此可作為農(nóng)業(yè)面源污染重點(diǎn)管控對象加強(qiáng)監(jiān)督和管理。

2.2 三江一級支流水質(zhì)時(shí)空演變特征

2023年三江一級支流各水質(zhì)因子呈季節(jié)變化趨勢（圖4）。COD的平均濃度全年整體保持穩(wěn)定，5—7月有上升的趨勢，之后又逐漸下降，全年中的最大平均值出現(xiàn)在7月，最小平均值出現(xiàn)在9月。總磷的平均濃度1—4月有輕微下降的趨勢，5—7月逐步上升到全年的峰值，8月又明顯下降，之后一直保持穩(wěn)定，全年中的最大平均值出現(xiàn)在7月，最小平均值出現(xiàn)在4月。氨氮的平均濃度1—3月有下降的趨勢，4—5月有明顯上升的趨勢，6—9月有下降的趨勢，9月到達(dá)低谷后逐漸回升，全年中的最大平均值出現(xiàn)在5月，最小平均值出現(xiàn)在9月。

寧波市32個(gè)三江一級支流控制單元斷面COD、氨氮、總磷污染空間變化特征如圖5所示，其中COD平均濃度范圍為1.99～7.78 mg·L^-1，氨氮平均濃度范圍為0.08～4.50 mg·L^-1，總磷平均濃度范圍為0.10～0.46 mg·L^-1。方橋街道、姜山鎮(zhèn)、高橋鎮(zhèn)、蛟川街道控制單元范圍內(nèi)三江一級支流COD 平均濃度均超過7.0 mg·L^-1，而明樓街道、洪塘街道、段塘街道氨氮、總磷平均濃度均較高。橫街鎮(zhèn)范圍內(nèi)三江一級支流COD、氨氮平均濃度均為最低，而小港街道范圍內(nèi)總磷平均濃度最低。

如表1所示，通過對寧波市不同控制單元及不同月份各污染物間相關(guān)性進(jìn)行分析，發(fā)現(xiàn)在不同控制單元中，COD 與氨氮之間存在顯著的中低度正相關(guān)關(guān)系（r=0.372*），而COD與總磷之間不存在顯著相關(guān)關(guān)系（r=0.291）。氨氮與總磷之間的相關(guān)性則極為顯著，表現(xiàn)出高度正相關(guān)關(guān)系（r=0.920***）。在不同月份的分析中，COD 與氨氮之間表現(xiàn)出顯著的中度正相關(guān)關(guān)系（r=0.564*），COD與總磷之間則具有顯著的中度正相關(guān)關(guān)系（r=0.698**），而氨氮與總磷之間不存在顯著相關(guān)關(guān)系（r=0.390）。

整體而言，三江一級支流各監(jiān)測斷面COD、氨氮、總磷變化趨勢相對一致，在初汛期水環(huán)境污染物濃度呈上升趨勢，之后逐漸下降。由于降雨是農(nóng)業(yè)面源污染向水環(huán)境遷移的主要因素[13-14]，結(jié)合寧波市亞熱帶季風(fēng)氣候條件和多數(shù)三江一級支流建有水閘的現(xiàn)狀，推斷農(nóng)業(yè)面源污染物于枯水期在河道、溝渠淤積，初汛期時(shí)由于地表徑流的作用將淤積的污染物沖至一級支流，導(dǎo)致水環(huán)境質(zhì)量下降[15]，之后由于主汛期降雨量大，導(dǎo)致一級支流污染物濃度被稀釋而下降[16]。

通過對三江一級支流斷面水污染物間相關(guān)性進(jìn)行分析發(fā)現(xiàn)，水環(huán)境中COD、氨氮、總磷污染存在一定相關(guān)性。楊宇峰等[17]基于隨機(jī)森林模型對遼河進(jìn)行了高時(shí)間分辨率氮、磷濃度模擬與預(yù)測，發(fā)現(xiàn)斷面總氮、總磷含量的變化趨勢高度一致。本研究與前人研究結(jié)果一致，可能的原因是農(nóng)村生活污水、畜禽糞污、農(nóng)田退水等主要農(nóng)業(yè)面源污染來源均具有復(fù)合型污染特征，含有氮、磷、COD等多種污染物，在相似的污染排放驅(qū)動(dòng)因素（如降雨等）影響下，通過相似的遷移機(jī)制進(jìn)入水環(huán)境，使水環(huán)境中多種污染物呈現(xiàn)相似的變化趨勢[18-19]。

2.3 面源污染負(fù)荷與斷面水質(zhì)響應(yīng)關(guān)系分析

通過對32個(gè)控制單元的單位國土面積農(nóng)業(yè)面源污染負(fù)荷與河流水質(zhì)進(jìn)行相關(guān)性分析發(fā)現(xiàn)，兩者不存在顯著相關(guān)關(guān)系?？赡艿脑蚴峭寥缹r(nóng)業(yè)面源污染具有消納作用[20-21]，部分控制單元面源污染負(fù)荷較高但消納面積充足，因此對水環(huán)境質(zhì)量影響較小。采用單位國土面積農(nóng)業(yè)面源污染負(fù)荷進(jìn)行表征，有助于更精準(zhǔn)地評估農(nóng)業(yè)活動(dòng)對水環(huán)境的影響。

通過對寧波市三江一級支流32個(gè)控制單元污染負(fù)荷強(qiáng)度與水環(huán)境質(zhì)量進(jìn)行相關(guān)性分析發(fā)現(xiàn)，兩者之間并沒有顯著相關(guān)性。可能的原因是寧波市城鎮(zhèn)化程度整體較高，尤其是明樓街道、段塘街道、白鶴街道、文教街道等控制單元農(nóng)業(yè)面源污染負(fù)荷較低，但其工業(yè)生產(chǎn)較為發(fā)達(dá)，人口密度也相對較高，其水環(huán)境質(zhì)量仍然較差。因此整體而言寧波市32個(gè)控制單元農(nóng)業(yè)源污染物對水質(zhì)污染貢獻(xiàn)相對較少[22]。對比圖3和圖5可知，部分控制單元單位國土面積農(nóng)業(yè)面源污染負(fù)荷與各控制單元水質(zhì)呈現(xiàn)對應(yīng)關(guān)系，單位國土面積農(nóng)業(yè)面源污染負(fù)荷較高的控制單元如姜山鎮(zhèn)、高橋鎮(zhèn)、洪塘街道控制單元，其水環(huán)境質(zhì)量整體也相對較差。該類控制單元的共同特點(diǎn)是農(nóng)業(yè)生產(chǎn)發(fā)達(dá)而工業(yè)生產(chǎn)和農(nóng)村生活污水貢獻(xiàn)相對較少，因此斷面水質(zhì)對農(nóng)業(yè)面源污染負(fù)荷的響應(yīng)程度相對較高，這與湖南省農(nóng)業(yè)面源污染與農(nóng)村水環(huán)境質(zhì)量的關(guān)系相似[7]。

本研究在32個(gè)控制單元中根據(jù)污染負(fù)荷強(qiáng)度進(jìn)一步選擇15個(gè)農(nóng)業(yè)生產(chǎn)活動(dòng)強(qiáng)度較高的控制單元分析其污染負(fù)荷強(qiáng)度與水質(zhì)的響應(yīng)關(guān)系，發(fā)現(xiàn)水環(huán)境污染物濃度與污染負(fù)荷強(qiáng)度存在顯著正相關(guān)關(guān)系（圖6）。分析結(jié)果表明，COD 污染負(fù)荷強(qiáng)度與水環(huán)境COD 濃度呈顯著正相關(guān)關(guān)系（Plt;0.05），判定系數(shù)為0.27，其中10 個(gè)控制單元處于95% 置信區(qū)間。與COD 污染情況類似，農(nóng)業(yè)源氨氮負(fù)荷強(qiáng)度與水環(huán)境氨氮濃度、農(nóng)業(yè)源總磷污染負(fù)荷強(qiáng)度與總磷濃度均呈顯著正相關(guān)關(guān)系，其判定系數(shù)分別為0.30與0.35，分別有8個(gè)與7個(gè)控制單元處于95% 置信區(qū)間。研究結(jié)果說明在農(nóng)業(yè)生產(chǎn)活動(dòng)強(qiáng)度較高的控制單元范圍內(nèi)農(nóng)業(yè)面源污染是導(dǎo)致水環(huán)境污染物濃度上升的主要原因之一，水質(zhì)對污染負(fù)荷響應(yīng)程度較高。

水環(huán)境污染來源復(fù)雜多樣且時(shí)空差異較大，雖然當(dāng)前農(nóng)業(yè)面源污染成為我國水環(huán)境氮磷污染的主要來源，但農(nóng)業(yè)源污染物排放情況受農(nóng)業(yè)生產(chǎn)方式、自然條件、污染防控措施等因素影響，各種因素的綜合影響使農(nóng)業(yè)面源污染呈現(xiàn)隨機(jī)性、隱蔽性和不確定性等特點(diǎn)，因此其對水質(zhì)的影響難以評估。吳根義團(tuán)隊(duì)[7，23]研究發(fā)現(xiàn)湖南單位國土面積農(nóng)業(yè)面源污染負(fù)荷與水質(zhì)類別顯著正相關(guān)，且開發(fā)了根據(jù)單位國土面積等標(biāo)污染負(fù)荷預(yù)測水環(huán)境質(zhì)量狀況的一階指數(shù)函數(shù)增長模型。本研究發(fā)現(xiàn)寧波市單位國土面積農(nóng)業(yè)面源污染負(fù)荷與水環(huán)境污染物濃度并無顯著相關(guān)性，但農(nóng)業(yè)生產(chǎn)活動(dòng)強(qiáng)度的較高控制單元水環(huán)境質(zhì)量對農(nóng)業(yè)面源污染負(fù)荷的響應(yīng)程度較高，說明基于農(nóng)業(yè)面源污染負(fù)荷的相關(guān)水質(zhì)預(yù)測模型有其特定的適用場景，在城市化程度較高的城市適用性相對較低。胡海珍[24]基于工業(yè)源和城鎮(zhèn)生活源（污水處理廠）建立了常州市平原河網(wǎng)的水動(dòng)力-水質(zhì)模型，各污染物濃度過程線實(shí)測值與模擬值基本吻合，說明在工業(yè)廢水、城鎮(zhèn)生活污水排放量較高的地區(qū)更適合以工業(yè)源和城鎮(zhèn)生活源污染負(fù)荷為基礎(chǔ)預(yù)測水環(huán)境質(zhì)量。寧波市作為經(jīng)濟(jì)發(fā)達(dá)城市，其工業(yè)源、生活源污染整體貢獻(xiàn)較多，同時(shí)寧波市部分控制單元農(nóng)業(yè)面源污染貢獻(xiàn)比較大，農(nóng)業(yè)污染負(fù)荷與水質(zhì)響應(yīng)關(guān)系較強(qiáng)，因此需要對寧波市不同控制單元水環(huán)境污染進(jìn)行負(fù)荷評估及源解析，劃定優(yōu)先治理區(qū)域，并針對其主要污染來源制定相對應(yīng)的污染管理政策和防治措施。

3 結(jié)論

（1）寧波市農(nóng)業(yè)面源污染負(fù)荷存在顯著的空間差異，中南部控制單元化學(xué)需氧量（COD）、氨氮、總磷污染負(fù)荷均相對較高，其中姜山鎮(zhèn)、高橋鎮(zhèn)控制單元單位國土面積污染負(fù)荷最高；農(nóng)村生活污水是農(nóng)業(yè)源COD和氨氮排放的主要來源，其貢獻(xiàn)率分別為57.3%和61.0%；種植業(yè)對總磷的污染貢獻(xiàn)率達(dá)到76.7%。

（2）三江一級支流水質(zhì)呈現(xiàn)季節(jié)變化趨勢，初汛期水環(huán)境污染物濃度呈上升趨勢；姜山鎮(zhèn)、方橋街道控制單元河流COD濃度較高；白鶴街道、洪塘街道、文教街道控制單元河流氨氮、總磷含量較高。

（3）寧波市32個(gè)三江一級支流控制單元單位國土面積污染負(fù)荷與水環(huán)境質(zhì)量整體無顯著相關(guān)性；農(nóng)業(yè)生產(chǎn)活動(dòng)強(qiáng)度較高的15個(gè)控制單元，單位國土面積COD、氨氮、總磷污染負(fù)荷與水環(huán)境污染物濃度均呈顯著相關(guān)關(guān)系，其判定系數(shù)分別為0.27、0.30 及0.35，說明水質(zhì)對污染負(fù)荷響應(yīng)程度較高。