基于SWAT模型的釜溪河流域非点源污染负荷模拟研究.pdf

第一作者:刘 伟,男,1 9 7 9年生,博士,副教授,主要从事环境规划与管理研究。*四川省科学技术厅“丝路经济带”生态脆弱区产业耗水效应及水资源承载力研究(N o.2 0 2 1 J D R 0 1 3 1);四川省科学技术厅四川省典型流域水环境动态管理体系研究(N o.2 0 1 7 S Z 0 1 7 4)。基于S WA T模型的釜溪河流域非点源污染负荷模拟研究*刘 伟 李逢港 李沁芮 杜 孟 汪嘉杨(成都信息工程大学资源环境学院,四川 成都6 1 0 2 2 5)摘要 非点源污染是流域污染主要来源之一。利用S WA T模型模拟分析釜溪河流域非点源氮磷污染负荷特征并识别关键源区,通过设置8种情景并结合SWA T模型评估各情景对非点源污染负荷的削减效果。结果显示:总磷、总氮及径流的率定验证结果均达到相关系数(R2)0.6、纳什系数(En s)0.5,表明S WA T模型的模拟值与观测值有较高的一致性。流域非点源污染负荷雨季高于旱季,下游大于上游,并确定8、9、1 1、1 3、1 4号子流域为关键源区。流域尺度上各情景对关键源区非点源污染负荷的削减率表现为:情景8(1 0m植被缓冲带)情景7(5m植被缓冲带)情景6(植草水道)情景5(残茬覆盖)情景3(化肥减施1 0%)情景4(免耕)情景2(化肥减施5%)。1 0m植被缓冲带削减率可达3 0%,残茬覆盖、化肥减施1 0%削减率大于1 0%,而免耕及化肥减施5%削减率大多低于1 0%。关键词 非点源污染负荷 S WA T模型 情景 模拟 釜溪河流域 D O I:1 0.1 5 9 8 5/j.c n k i.1 0 0 1-3 8 6 5.2 0 2 3.0 8.0 2 2S i m u l a t i o no fn o n-p o i n t s o u r c ep o l l u t i o n l o a d i n t h eF u x iR i v e rB a s i nb a s e do nS WA Tm o d e l L I U W e i,L IF e n g g a n g,L I Q i n r u i,DU M e n g,WANG J i a y a n g.(C o l l e g e o f R e s o u r c e s a n d E n v i r o n m e n t,C h e n g d u U n i v e r s i t y o fI n f o r m a t i o nT e c h n o l o g y,C h e n g d uS i c h u a n6 1 0 2 2 5)A b s t r a c t:N o n-p o i n ts o u r c ep o l l u t i o ni so n eo ft h e m a i np o l l u t i o ns o u r c e so fb a s i np o l l u t i o n.B a s e do nt h eS WA Tm o d e l,t h en o n-p o i n t s o u r c ep o l l u t i o n l o a dc h a r a c t e r i s t i c so f t h eF u x iR i v e rB a s i nw e r es i m u l a t e d,a n dt h ek e ys o u r c ea r e a so f t h eb a s i nw e r e i d e n t i f i e d.B ys e t t i n g8s c e n a r i o sa n dc o m b i n i n gw i t ht h eSWA T m o d e l,t h er e d u c t i o ne f f e c to f e a c hs c e n a r i oo nn o n-p o i n t s o u r c ep o l l u t i o n l o a dw a se v a l u a t e d.T h e r e s u l t s s h o w e dt h a t t h e s i m u l a t e dv a l u e so fS WA T m o d e lw e r eh i g h l yc o n s i s t e n tw i t ht h em e a s u r e do n e s.T h ec a l i b r a t i o nr e s u l t so ft o t a lp h o s p h o r u s,t o t a ln i t r o g e na n dr u no f fa c h i e v e dR20.6,n a s hc o e f f i c i e n t(En s)0.5.T h en o n-p o i n ts o u r c ep o l l u t i o nl o a di nr a i n ys e a s o nw a sh i g h e r t h a nt h a t i nd r ys e a s o n,a n d t h en o n-p o i n t s o u r c ep o l l u t i o n l o a d i n t h ed o w n s t r e a mo f t h eb a s i nw a sl a r g e r t h a nt h a t i n t h eu p s t r e a m.S u b-b a s i n s 8,9,1 1,1 3a n d1 4w e r e i d e n t i f i e da sk e y s o u r c e a r e a s.O n t h eb a s i ns c a l e,t h er e d u c t i o nr a t e so fn o n-p o i n t s o u r c ep o l l u t i o nl o a d s i nk e ys o u r c ea r e a sw e r er a n k e da sf o l l o w s:S c e n a r i o8(1 0mv e g e t a t i o nb u f f e r z o n e)S c e n a r i o7(5mv e g e t a t i o nb u f f e r z o n e)S c e n a r i o6(g r a s sp l a n t i n gw a t e r w a y)S c e n a r i o5(r e s i d u a l s t u b b l e c o v e r a g e)S c e n a r i o3(1 0%r e d u c t i o n i n f e r t i l i z e r a p p l i c a t i o n)S c e n a r i o4(n o t i l l a g e)S c e n a r i o2(5%r e d u c t i o n i n f e r t i l i z e r a p p l i c a t i o n).T h ev e g e t a t i o nb u f f e r z o n ew i t haw i d t ho f 1 0mh a dar e d u c t i o nr a t eo f 3 0%f o rn o n-p o i n t s o u r c ep o l l u t i o nl o a d s.H o w e v e r,t h er e d u c t i o nr a t eo fr e s i d u a ls t u b b l ec o v e r a g ea n d1 0%r e d u c t i o ni nf e r t i l i z e ra p p l i c a t i o n w a sg r e a t e rt h a n1 0%,w h i l et h er e d u c t i o nr a t eo fn ot i l l a g ea n d5%r e d u c t i o ni nf e r t i l i z e ra p p l i c a t i o nw a sg e n e r a l l y l e s s t h a n1 0%.K e y w o r d s:n o n-p o i n t s o u r c ep o l l u t i o nl o a d;SWA Tm o d e l;s c e n a r i o;s i m u l a t i o n;F u x iR i v e rB a s i n 水是人类赖以生存的基础自然资源,是制约社会经济发展的重要因素。随着社会经济和城镇化提速发展,水资源需求量日益增长,生态问题频频发生,水污染问题日益严重1-2。如何有效控制水污染已成为人们关注的重点3-4。非点源污染大多由氮、磷等营养物质和土壤颗粒等组成,在产生并排放后,经有效降雨和径流运载等一系列迁移转化过程进入水体5。其形成是一个复杂且随机的多因素影响过程,具有广泛性、模糊性和滞后性等特点6,使得准确量化非点源污染的“源”和“汇”存在一定难度。非点源污染关乎流域水环境与水生态安全,成为流域水环境研究中共同关注的热点7-8。6711 环境污染与防治 第4 5卷 第8期 2 0 2 3年8月 国内外研究人员针对非点源污染及防治做了大量研究。部分学者对非点源污染来源开展了研究,得出农业对非点源污染的贡献最大,主要包括农业种植、农村生活源、畜禽养殖9-1 0。SWAT模型在长时间尺度和复杂水文条件下模拟具有较大优势,被广泛应用于河流、湖泊、水库非点源污染模拟1 1-1 6。有研究表明,SWAT模型还能与管理措施联用,有效评估管理措施对非点源污染的削减效果1 7-1 9。釜溪河是沱江七大支流中长度最长、流域面积最大、人口最多的一级支流,也是四川省污染较为严重的河流之一,水体长期呈污染状态,面临严重的非点源污染问题2 0。近年来,釜溪河流域水质有所好转,但碳研所断面氨氮浓度常高于 地表水环境质量标准(G B3 8 3 82 0 0 2)类标准限值,难以稳定达标2 1。以往研究主要针对釜溪河水环境现状,分析流域水质污染规律。本研究通过构建釜溪河流域SWAT非点源污染模型,开展流域非点源污染负荷模拟研究,在此基础上设置8个情景,评估各方案对非点源污染负荷的削减效果,提出非点源氮磷污染控制对策,为釜溪河流域及其他流域开展非点源污染综合治理提供科学依据和技术支撑。1 研究区概况 釜溪河流域地处四川盆地南部,介于1 0 4 1 8 E至1 0 5 1 E、2 9 6 N至2 9 4 5 N之间,釜溪河属于沱江下游右岸一级支流。流域面积34 7 2.0 0k m2,主要包括威远河、旭水河、长滩河及镇溪河等支流,见图1。流域属于亚热带湿润气候,年平均气温1 8,年平均降水量在10 0 010 8 0mm之间,集中在69月,年平均流量4 2.4m3/s。流域碳研所断面水质难以稳定达标,水质波动大,雷公滩断面水质达标困 难。截 至2 0 2 1年,釜 溪 河 流 域 农 村 人 口 为1 5 7.4 8万人,农村生活污染源分布广泛,生活污水较难集中处置,主要通过土地漫流进入河流。同时流域耕地面积广,农业非点源污染物排放强度大2 2-2 3。2 研究方法2.1 模型构建方法2.1.1 SWAT模型简介 SWAT模型由美国农业部农业研究中心开发,是综合系统理论和技术产生的分布式模型,包括水文循环、土壤侵蚀和污染负荷子模块,可在长时间尺度上对不同管理措施影响下径流、氮磷等营养物质输移进行模拟和预测2 4。图1 釜溪河流域示意图F i g.1 S k e t c hm a po f t h eF u x iR i v e rB a s i n2.1.2 模型数据库构建 模型数据库主要由空间及非空间数据组成,空间数据用于表达位置、方向,非空间数据又叫属性数据,起补充说明作用2 5。用于SWAT模型建立、率定和验证的数据来源见表1。考虑气象数据的可得性,故模拟时段为2 0 0 82 0 1 8年,其中2 0 0 82 0 1 0年为模型预热期、2 0 1 12 0 1 5年为率定期、2 0 1 62 0 1 8年为验证期。表1 S WA T模型输入数据及来源T a b l e1 I n p u td a t aa n ds o u r c e f o rS WA T m o d e l数据分辨率来源数字高程模型(D EM)2 0m地理空间数据云土地利用类型1k m中国科学院资源环境科学与技术中心土壤类型1k m世界土壤数据库(HWS D)气象数据0.2 5 中国大气同化驱动集(CMA D SV 1.0)其他数据自贡统计年鉴、内江统计年鉴2.1.3 模型率定与验证 利用SWAT-C U P软件进行参数敏感性分析,借助不确定性适应算法(S U F I-2)中t检验值、p值对模型参数进行率定2 6。利用相关系数(R2)、纳什系数(En s)等评价指标评价模型适用性2 7。由于R2侧重于评价相关性,仅采用R2无法对系统误差进行识别;En s侧重于评价拟合度,仅采用En s会导致部分值较小的模拟误差被忽略,因此,为提高模型在流域的适用性,本研究同时采用上述两个指标进行模型7711刘 伟等 基于S WA T模型的釜溪河流域非点源污染负荷模拟研究率定和验证,当R20.6且En s0.5时,模型模拟准确性较好,模拟结果可以接受2 8。2.2 情景设置与评估 考虑釜溪河流域的实际情况与非点源污染负荷特征,将模型现状设为情景1,同时设置化肥削减、改变耕作方式、增加植被缓冲带或植草水道等情景,见表2。在SWAT模型中,通过改变参数和变量来构建不同情景,利用非点源污染负荷削减率作为评价指标,削减率指各情景下非点源污染负荷削减量与基础情景非点源污染负荷的比值2 9。3 结果与分析3.1 模型构建结果与分析3.1.1 SWAT模型构建结果 在SWAT模型中,以流域D EM为基础,将土地利用、土壤数据进行提取和重分类,得出釜溪河流域土地利用类型(分为耕地、林地、草地、水域、城乡工矿居民用地及未利用土地),其中耕地占比最大。土壤类型分为艳色高活性淋溶土、不饱和雏形土、石灰性疏松岩性土、人为堆积土、简育高活性强酸土等。以釜溪河入沱江干流处为研究区出口,将整个流域划分为1 5个子流域、3 0 9个水文响应单元。各子流域水文响应单元见表3。3.1.2 SWAT模型率定与验证 月均径流模拟值与观测值的对比结果见图2。由图2可知,多年月均径流拟合效果较好,率定期R2=0.7 7,En s=0.7 1,验证期R2=0.7 3,En s=0.6 9,符合模型的精度要求。流域总氮、总磷负荷模拟值与观测值如图3、图4所示,模拟值与观测值的相关性和拟合度较好,率定期邓关断面总氮负荷R2=0.7 2,En s=0.6 5;总磷负荷R2=0.7 5,En s=0.6 2;验证期总氮负荷R2=0.7 3,En s=0.6 8;总磷负荷R2=0.7 1,En s=0.6 3;模型模拟结果均符合精度要求。综上,采用SWAT模型对釜溪河流域非点源污染负荷模拟具有 精 确 性 和 合 理 性。模 型 敏 感 性 参 数 见表4。表2 情景及相关参数设置T a b l e2 S c e n a r i o sa n dr e l a t e dp a r a m e t e r ss e t t i n g情景类型措施设置参数调整1基础情景无新增措施无2化肥施用减少5%修改m g t文件中参数3化肥施用减少1 0%修改m g t文件中参数4非工程措施免耕m g t文件中修改耕作操作5残茬覆盖m g t文件中添加仅收割操作,径流曲线数=-2,水土保持措施因子=0.2 9,坡面漫流曼宁系数=0.36植草水道修改植草水道中参数7工程措施植被缓冲带(5m)o p s文件中植被缓冲带宽度设为5m8植被缓冲带(1 0m)o p s文件中植被缓冲带宽度设为1 0m表3 水文响应单元划分情况T a b l e3 D i v i s i o no fh y d r o l o g i c a l r e s p o n s eu n i t s子流域面积/k m2面积占比/%水文响应单元/个平均高程/m12 5 9.6 77.8 54 05 9 5.8 622 3 0.6 06.9 73 35 3 1.7 632 1 7.3 06.5 73 44 0 4.4 841 2 8.2 83.8 81 93 8 0.6 452 1 5.0 36.5 02 03 9 7.2 261 2 8.8 13.8 963 5 1.6 372 2 0.0 66.6 553 5 1.1 881 3 3.0 24.0 21 13 3 5.9 891 1.1 00.3 41 83 2 0.8 91 01 4 7.7 44.4 783 1 7.7 11 16 6.2 12.0 02 33 1 7.6 11 27 9 4.7 52 4.0 24 14 1 1.9 91 35 1 2.1 81 5.4 82 83 2 4.1 91 45.1 80.1 632 6 9.6 11 52 3 8.8 27.2 22 03 3 4.9 28711 环境污染与防治 第4 5卷 第8期 2 0 2 3年8月图2 率定期与验证期径流模拟结果F i g.2 S i m u l a t i o nr e s u l t so f r u n o f f i nc a l i b r a t i o na n dv a l i d a t i o np e r i o d 注:月份序号中的1为2 0 1 1年1月,9 6为2 0 1 8年1 2月,其余类推;图4同。图3 率定期与验证期总氮负荷模拟结果F i g.3 S i m u l a t i o nr e s u l t so f t o t a ln i t r o g e nl o a d i nc a l i b r a t i o na n dv a l i d a t i o np e r i o d图4 率定期与验证期总磷负荷模拟结果F i g.4 S i m u l a t i o nr e s u l t so f t o t a lp h o s p h o r u s l o a d i nc a l i b r a t i o na n dv a l i d a t i o np e r i o d表4 釜溪河流域S WA T模型敏感性参数T a b l e4 S e n s i t i v ep a r a m e t e r so fSWA Tm o d e l i nF u x iR i v e rB a s i n参数代码描述t检验值p值C N 2S C S径流曲线系数4.7 40E S C O土壤蒸发补偿系数3.8 90S O L_B D表层土壤湿容重2.6 00.0 1GW_D E L AY地下水滞后系数-1.9 40.0 5S O L_A L B潮湿土壤反映率-1.7 30.0 8C H_N 2主河道曼宁系数1.4 20.1 6S L S U B B S N平均坡长-1.3 60.1 8A L P HA_B NK河道调蓄系数1.2 90.1 9E P C O植被吸水补偿系数-1.2 00.2 3B I OM I X生物混合效率系数-0.9 00.3 7E R O R GN氮富集率-2.9 20.5 3N_U P D I S氮吸收分布参数-2.7 50.5 9N P E R C O氮渗透系数1.7 20.6 0S D N C O反硝化土壤含水量阈值1.3 60.9 0C D N反硝化指数速率系数-0.2 30.9 5P HO S K D土壤磷分配系数-2.0 20.4 9P P E R C O磷渗透系数1.6 80.6 2P S P磷有效性指数1.4 70.7 2E R O R G P磷富集率1.3 30.8 83.2 模拟结果与分析3.2.1 时间分布特征 根据模型运行结果统计流域内非点源污染负荷(总氮、总磷负荷)的月均值,结果见表5。总氮、总磷负荷随时间变化较为明显,雨季的总氮、总磷负荷明显高于旱季。除9月外,整体呈先上升后下降的趋势,1月至3月总氮、总磷负荷均呈较低状态;4月至6月总氮、总磷负荷出现上升,原因可能是该时间段作物处于生长期,施肥量大,且氮、磷等营养元素容易随降雨迁移入河;7月至8月作物逐渐生长,叶面积增大,氮、磷不易被降雨冲刷入河;9月由于新一轮作物耕作,施肥量有所增长,导致污染负荷上升,1 0月至1 2月污染负荷变化程度较小。表5 月均非点源污染负荷模拟结果T a b l e5 M o n t h l ya v e r a g en o n-p o i n t s o u r c ep o l l u t i o n l o a ds i m u l a t i o nr e s u l t s月份总氮负荷/t总磷负荷/t11 0 4.5 27.7 221 2 2.6 39.9 331 3 7.1 51 4.6 443 4 6.0 72 4.7 857 2 3.7 03 7.8 966 1 3.0 43 7.2 874 1 0.4 32 2.3 882 0 1.8 41 5.3 692 8 9.8 92 0.8 21 01 4 2.5 11 3.7 41 11 1 5.4 01 0.2 41 29 7.3 08.7 89711刘 伟等 基于S WA T模型的釜溪河流域非点源污染负荷模拟研究3.2.2 空间分布特征 以1 5个子流域为基本单元,结合A r c G I S软件分析2 0 1 12 0 1 8年总氮、总磷年平均负荷空间分布特征,结果如图5、图6所示。流域干流处的子流域总氮负荷较其他子流域高,主要包括1 1、1 3、1 4号子流域,该区域地势较为平缓,人口密集,且位于流域下游,所以总氮负荷较大。下游处的子流域总磷负荷高于上游子流域,1 3、1 4、1 5号子流域的总磷负荷在流域中较大,可能是因为磷肥不易被植物吸收,导致磷素在河道中累积,同时1 4号子流域为整个釜溪河流域的出水口,磷随着水流汇集在此处。图5 总氮年平均负荷分布F i g.5 D i s t r i b u t i o no f a n n u a l a v e r a g e t o t a l n i t r o g e n l o a d3.2.3 关键源区识别 采用流失强度指数法识别流域关键源区,以总氮、总磷单位面积流失强度作为评价指标,采用自然裂点分级法对各子流域氮、磷流失强度进行等级划分,分级标准见表6。依据分析结果,将流失强度高的子流域作为流域非点源污染首要治理区域。釜溪河流域总氮、总磷关键源区划分结果见图7、图8。由图7、图8可知,流失强度高的子流域分布在图6 总磷年平均负荷分布F i g.6 D i s t r i b u t i o no f a n n u a l a v e r a g e t o t a lp h o s p h o r u s l o a d威远河和干流所在区域,8、9、1 1号子流域的总氮流失强度为重度流失,3、1 3、1 4号子流域总氮为较重流失,6个子流域面积之和占流域面积的2 8.5 6%,但总氮负荷占比为7 0.8 5%。8、9、1 1、1 3、1 4、1 5号子流域的总磷流失较为严重,6个子流域占整个流域2 9.2 1%,但总磷负荷占比为7 4.2 9%。综上,研究区域内非点源污染关键源区为8、9、1 1、1 3、1 4号子流域,5个子流域面积占比仅2 1.9 9%,但输出的总氮负荷占比为6 2.3 0%,总磷负荷占比为6 0.5 4%,占比均超过一半。3.3 不同情景评估结果分析 各情景下关键源区(8、9、1 1、1 3、1 4号子流域)非点源氮、磷污染负荷的削减效果见表7。由表7可知,各情景削减率表现为情景8(1 0m植被缓冲带)情景7(5 m植被缓冲带)情景6(植草水道)情景5(残茬覆盖)情景3(化肥减施1 0%)情景4(免耕)情景2(化肥减施5%)。从单个情景来看,情景6、情景7和情景8对总氮、总磷负荷削 减效果较好,总氮负荷削减率分别为1 6.4 1%、表6 关键源区评价指标等级划分T a b l e6 C l a s s i f i c a t i o no f e v a l u a t i o n i n d i c a t o r s f o rk e ys o u r c ea r e a s项目分级标准总氮单位面积流失强度/(k g/h m2)1.5 62.9 72.9 76.3 66.3 61 2.9 71 2.9 71 8.0 71 8.0 72 5.6 8总磷单位面积流失强度/(k g/h m2)0.1 00.1 70.1 70.2 70.2 70.5 00.5 01.0 71.0 71.5 8流失强度等级轻度较轻中度较重重度0811 环境污染与防治 第4 5卷 第8期 2 0 2 3年8月图7 总氮流失强度等级F i g.7 I n t e n s i t y l e v e l o f t o t a l n i t r o g e n l o s s2 2.9 4%、3 0.3 1%,总磷负荷削减率分别为5 3.7 8%、6 1.5 8%、6 7.2 0%。其次为情景3、情景5,总氮、总磷负荷削减率均大于1 0%,情景2、情景4总氮负荷削减率在1 0%以下。从类别看,情景68(工程措施)的削减效果均优于情景25(非工程措施),这是由于工程措施可以让污染负荷迅速降低,而非工程措施通常实施起来周期较长,污染负荷削减缓慢。4 结 论 (1)通过构建釜溪河流域SWAT模型,将研究区划分为1 5个子流域和3 0 9个水文响应单元,月尺度上径流、总氮、总磷率定验证结果均达到R20.6且En s0.5,说明模型模拟值与观测值有较高一致性,模型在釜溪河流域适用性较好。图8 总磷流失强度等级F i g.8 I n t e n s i t y l e v e l o f t o t a lp h o s p h o r u s l o s s (2)对流域非点源污染负荷特征进行分析,时间上,总氮、总磷负荷集中在雨季,主要受到耕种施肥和雨水冲刷等的影响;空间上非点源污染负荷具有一定差异性,下游大于上游,主要由地势及水流汇集决定。由于8、9、1 1、1 3、1 4号子流域的总氮负荷占流域总量的6 2.3 0%,总磷负荷占比为6 0.5 4%,因此确定这5个子流域为釜溪河流域非点源污染首要治理区域。(3)对各情景下非点源污染负荷削减效果进行评估,各情景削减率表现为情景8(1 0m植被缓冲带)情景7(5 m植被缓冲带)情景6(植草水道)情景5(残茬覆盖)情景3(化肥减施1 0%)情景4(免耕)情景2(化肥减施5%)。建议可针对流域实际情况进行分区治理,采用不同措施,在非表7 关键源区各情景下的削减效果T a b l e7 R e d u c t i o ne f f e c t su n d e rd i f f e r e n t s c e n a r i o s i nk e ys o u r c ea r e a s情景总氮负荷削减量/(t/a)削减率/%总磷负荷削减量/(t/a)削减率/%情景28 5.2 53.7 82 4.0 77.7 8情景32 6 5.2 81 0.3 72 9.0 81 3.3 8情景41 2 8.3 25.4 42 5.6 81 0.2 1情景52 8 7.6 51 1.5 68 4.5 45 2.4 3情景64 0 9.6 61 6.4 18 7.3 95 3.7 8情景75 4 8.3 92 2.9 41 0 4.1 26 1.5 8情景87 4 1.7 53 0.3 11 1 3.8 66 7.2 01811刘 伟等 基于S WA T模型的釜溪河流域非点源污染负荷模拟研究点源污染负荷流失强度较大的地区可试点采用植被缓冲带,同时在农业种植较多的区域采取削减化肥使用、残茬覆盖和生态循环等自然措施,以此来削减非点源污染负荷。参考文献:1 P A S TO R A V,L UDW I G F,B I EMAN S H,e ta l.A c c o u n t i n gf o re n v i r o n m e n t a lf l o wr e q u i r e m e n t si ng l o b a lw a t e ra s s e s s-m e n t sJ.H y d r o l o g ya n dE a r t hS y s t e mS c i e n c e s,2 0 1 4,8(1 2):5 0 4 1-5 0 5 9.2 徐文静,常国梁,时迪迪,等.北运河上游非点源污染风险空间分布特征研究J.环境科学学报,2 0 2 1,4 1(1):7-1 4.3 L I U Y W,WAN G W,HU Y M.I n v e s t i g a t i n gt h ei m p a c to fs u r f a c es o i lm o i s t u r ea s s i m i l a t i o no ns t a t ea n dp a r a m e t e re s t i-m a t i o n i nS WA T m o d e lb a s e do nt h ee n s e m b l eK a l m a nf i l t e ri nu p p e r H u a iR i v e rb a s i nJ.J o u r n a lo f H y d r o l o g y a n dH y d r o m e c h a n i c s,2 0 1 7,6 5(2):1 2 3-1 3 3.4 K I MJ,P A R KB,C HO I J,e t a l.O p t i m u md e t a i l e ds t a n d a r d s t oc o n t r o l n o n-p o i n t s o u r c ep o l l u t i o np r i o r i t ym a n a g e m e n ta r e a s:c e n t e r e do nh i g h l a n da g r i c u l t u r ew a t e r s h e dJ.S u s t a i n a b i l i t y,2 0 2 1,1 3(1 7):2-1 6.5 Z HAN GT,N IJP,X I ED T.A s s e s s m e n to ft h er e l a t i o n s h i pb e t w e e nr u r a ln o n-p o i n ts o u r c ep o l l u t i o na n de c o n o m i cd e v e l-o p m e n t i nt h eT h r e eG o r g e sR e s e r v o i rA r e aJ.E n v i r o n m e n t a lS c i e n c ea n dP o l l u t i o nR e s e a r c h,2 0 1 6,2 3(8):8 1 2 5-8 1 3 2.6 裴亮,刘慧明,莫家玉,等.南水北调中线工程农业非点源污染来源及控制措施J.生态经济,2 0 1 1,2 7(6):1 6 4-1 6 7.7 L I U R M,XU F,Z HAN G P P,e ta l.I d e n t i f y i n gn o n-p o i n ts o u r c ec r i t i c a ls o u r c ea r e a sb a s e do n m u l t i-f a c t o r sa tab a s i ns c a l ew i t hS WA TJ.J o u r n a lo fH y d r o l o g y,2 0 1 6,5 3 3:3 7 9-3 8 8.8 X I AN GCY,WAN GY,L I U H W.As c i e n t o m e t r i c s r e v i e wo nn o n p o i n t s o u r c ep o l l u t i o nr e s e a r c hJ.E c o l o g i c a lE n g i n e e r i n g,2 0 1 7,9 9:4 0 0-4 0 8.9 D E NN I S L C,P E T E RJ V,K E I TH L.M o d e l i n gn o n-p o i n ts o u r c ep o l l u t i o n i nv a d o s ez o n ew i t hG I SJ.E n v i r o n m e n tS c i-e n c ea n dT e c h n o l o g y,1 9 9 7,8:2 1 5 7-2 1 7 5.1 0 张水龙,庄季屏.农业非点源污染研究现状及发展趋势J.生态学杂志,1 9 9 8,1 7(6):5 1-5 5.1 1 D I NGY,D ON GF,Z HAOJY,e ta l.N o n-p o i n ts o u r c ep o l l u-t i o ns i m u l a t i o na n db e s tm a n a g e m e n tp r a c t i c e sa n a l y s i sb a s e do nc o n t r o l u n i t s i nN o r t h e r nC h i n aJ.I n t e r n a t i o n a l J o u r n a l o fE n v i r o n m e n t a lR e s e a r c ha n dP u b l i cH e a l t h,2 0 2 0,1 7(3):8 6 8-8 8 1.1 2 张京,郑华,何梦男,等.流域水环境污染模拟及关键源区鉴别 以义乌江流域为例J.环境工程学报,2 0 2 1,1 5(4):1 1 6 7-1 1 7 7.1 3 孙浩然,边睿,李若男,等.基于SWA T模型的磷负荷削减最佳管理措施(BMP s)评估研究J.环 境科学学报,2 0 2 0,4 0(7):2 6 2 9-2 6 3 7.1 4 T AHMA S E B IN M,G R I MM K,B A Z R KA R M H,e ta l.S WA Tm o d e l i n go fn o n-p o i n ts o u r c ep o l l u t i o ni nd e p r e s s i o n-d o m i n a t e db a s i n su n d e rv a r y i n gh y d r o c l i m a t i cc o n d i t i o n sJ.I n t e r n a t i o n a l J o u r n a lo fE n v i r o n m e n t a lR e s e a r c ha n dP u b l i cH e a l t h,2 0 1 8,1