四川山区农村面源污染负荷估算与评价

刘 菊，傅 斌，王玉宽，徐 佩



四川山区农村面源污染负荷估算与评价

刘 菊1,2，傅 斌1,3，王玉宽1,3※，徐 佩1,3

（1. 中国科学院、水利部成都山地灾害与环境研究所，成都 610041；2. 中国科学院大学，北京 100049； 3.万州典型区生态环境监测重点站，重庆 404020）

对农村面源污染物负荷的估算、来源分析及控制措施制定是当前研究的热点之一。利用输出系数模型和生态系统服务与权衡综合评价模型（Integrate valuation of ecosystem services and tradeoffs tool，InVEST模型）估算宝兴县面源污染物总氮（TN）、总磷（TP）负荷，并结合水环境功能分区标准，分析面源污染物的超标情况。旨在揭示山区农村面源污染治理的空间异质性，为农村面源污染防治规划、山区水环境治理效率的提升提供科学依据。结果表明，1）2010年宝兴县面源污染物TN、TP负荷分别为1 156.44、81.41 t/a，TN、TP负荷的空间格局整体上呈现出分布不均，局部集中，靠近水体的规律；2）不同污染源对TN、TP负荷的贡献率大小顺序均为：畜禽养殖>农业用地>农村生活；3）从流域尺度看，宝兴县TN、TP负荷没有超过III类水质标准规定的污染物负荷阈值；而从30 m×30 m栅格尺度看，宝兴县TN、TP超标量分别为763.63、51.16 t/a，其超标区域主要集中在耕地、草地、居民地以及靠近水体的地方，以上超标区域即为面源污染控制和治理的关键区域。

农业；污染；模型；水环境功能分区；产水量；空间差异

0 引 言

随着城市和工业点源污染逐步得到有效的控制，山区农村面源污染成为导致水环境恶化的主要原因之一[1]。同时农村面源污染已成为制约农业可持续发展和农村水环境质量改善的重要原因[2]，因此，量化农村面源污染以确定重点控制区域，是进行非点源污染管理的主要途径，同时也是研究的热点问题之一。

输出系数模型是应用较广的估算面源污染负荷的方法，在北京[3]、三峡库区[4]等地区都得到了应用。并且其原理简单，所需参数少，又具有一定的精度[5]，被认为是较为可靠的方法，可以模拟缺乏长时间系列监测数据的大中尺度流域面源污染的年负荷量[6]。此外，面源污染输出系数的应用并不局限于输出系数模型，事实上所有的输出系数模型和许多复杂的机理模型都依赖于污染物输出系数来估算大流域污染物的产生和输移[7]。然而，输出系数模型主要考虑人为因素对农村面源污染负荷的影响，并没有考虑地形地貌、降雨量、土壤、植被等自然因素对农村面源污染输出的影响。在众多影响因素中，降雨和地形是影响面源污染的关键因素[8]。蔡明等[9]在传统的输出系数模型中引入了降雨影响系数和流域损失系数2个参数，估算渭河流域的TN负荷量，取得了较好的估算结果，但其并未考虑地形因素对面源污染输出量的影响。有学者结合DEM、土地利用和降雨资料（即考虑地形和降雨）对三峡地区面源污染负荷进行研究，结果表明改进后的模型能够使模拟误差大幅度降低[10]。有些学者尝试着将分布式水文模型以及3S技术与输出系数模型相结合，同时考虑降雨的时间和空间差异以及地形对污染物输出的影响。如龙天渝等[4]结合SLURP分布式水文模型，在输出系数模型中引入了污染负荷系数，用以表征降水、产流等将流域上产生的非点源污染物转换成流域出口污染负荷的强弱程度。耿润哲等[11]使用SWAT分布式模型模拟评价密云水库流域非点源污染负荷分布特征，结果表明降雨的时空异质性是非点源污染产生的重要影响因素。由于简单的输出系数模型忽略了影响面源污染负荷的关键因素而复杂的机理模型对基础数据要求严格，故本研究使用输出系数模型与InVEST模型估算和评价山区面源污染物负荷。InVEST模型是由美国斯坦福大学、世界自然基金会和大自然保护协会联合开发的量化和价值化生态系统服务的综合模型[12]。其中，水质净化模块不仅考虑了降雨和地形因素对面源污染物的影响，还考虑了土壤和植被等对上游污染物的截留作用。该模型结合GIS工具，能够体现不同尺度的面源污染物输出的空间差异，适合于降雨空间分布不均以及多山地区流域面源污染输出的模拟。

除了估算与分析面源污染的排放来源及负荷以外，评估面源污染对水质的危害后果及程度也得到较多的关注。因为，在有限的经费和时间内，不可能对所有的区域都加以治理，必须分清轻重缓急，首先治理危害性最大的区域[13]。面源污染对水质的危害程度不仅取决于污染物负荷的大小，还取决于受纳水体的保护目标和纳污能力。严素定等[14-15]将面源污染排放与水环境功能分区挂钩，通过面源污染物排放量与水资源量相除得到污染物排放浓度，将此浓度与水环境功能分区的水质控制类别标准值相比得到面源污染水质指数，以反映面源污染对当地水环境产生的实际影响和污染压力，更为准确地判断农村面源污染的严重程度，从而为制定有针对性的面源污染防治对策与措施提供依据。

本研究在结合水环境功能分区的基础上，估算了面源污染物的超标情况，直观地揭示山区面源污染的空间异质性。首先利用输出系数模型，结合研究区土地利用、畜禽养殖和农村生活等实际情况，计算出研究区农村面源污染物TN和TP的负荷，以此为基础，借助InVEST模型的水质净化模块，在综合考虑地形地貌、降雨量、土壤和植被等因素的影响下，进行研究区面源污染物TN和TP负荷及其空间差异情况分析，最后，结合研究区水环境功能分区的水质控制类别标准值，从污染源排放与水质评估层面分析研究区TN和TP超标量及分布。本研究旨在揭示山区面源污染治理的空间异质性，确定其重点控制和治理区域，为农村面源污染防治规划、山区环境治理效率的提升提供科学依据。

1 研究区概况

宝兴县位于四川盆地西部边缘，地理位置为102°25′～103°02′E，30°09′～30°56′N。地形以山地为主，山地面积占幅员的99.7%，是典型的山区县。宝兴县属亚热带季风性湿润气候区，受山地地形影响（海拔范围750～5 328 m），气候出现垂直带谱，呈现典型立体气候，年平均气温14.1 ℃，年平均降雨量993.7 mm。县内植被保存完整，覆盖度达80%以上，系青衣江发源地。县内有宝兴河贯穿全境，流域面积3 010 km2，宝兴县县域范围与宝兴河流域范围基本一致。宝兴县是长江上游重要的水源涵养区和洪水控制区，也是全球生物多样性的热点地区。

宝兴县隶属于四川省雅安市，全县总面积3 114 km2，下辖9个行政乡（镇）、55个行政村，有17个少数民族。2010年末，全县共有18 175户，总人口58 335人，人口密度为18.73人/km2，其中农业人口45 980人，占总人口的78.82%。2010年，宝兴县GDP为14.67亿元，农民人均纯收入达5 456元。2010年全县畜禽养殖情况为牛存栏量3.84万头，羊存栏量3.63万头，猪出栏量5.85万头，家禽出栏量13.48万只。根据全国第二次土地调查数据，该县共有农用地2 993.79 km2，占全县的96.13%。农用地中林地（2 404.14 km2）面积最大，其次为草地（498.05 km2），而园地（36.00 km2）、旱地（51.98 km2）和水田（3.62 km2），3者共占农用地的2.94%。由于宝兴县自然村落较多，农民居住相对分散，生活垃圾、人畜粪便不能得到有效处理，农药和化肥的使用也造成了一定的面源污染。《宝兴县十二五生态环保规划》明确指出在创建国家级生态县中主要制约的因素就是农村面源污染。宝兴县水系及土地利用见图1。

2 研究方法

2.1 输出系数模型

2.1.1 模型介绍

输出系数模型所需参数少，操作简便，是较为可靠的方法用来模拟缺乏长时间系列监测数据的大中尺度流域面源污染的年负荷量。1996年Johns 等在前人研究的基础上，综合考虑土地利用类型、人口和牲畜数量等因素对面源污染的影响，提出了更为细致、完备的输出系数模型计算污染物输出负荷。计算方程为[3,16]

式中L为污染物在流域的总负荷量，kg/a；为流域中的土地利用类型，共种；E为污染物在第种土地利用类型中的输出系数，kg/hm2；或第种畜禽每头排污系数，kg/a；或人口每人输出系数，kg/a；A为流域中第种土地利用类型的面积，hm2；或第种畜禽数量，头或只；或人口数量，人；为由降雨输入的污染物总量，kg/a，本研究未考虑此项。主要原因如下：首先研究区缺少相关监测和研究，难以获得雨水中TN、TP含量或者大气湿沉降速率。其次，根据宝兴县环境保护局公布信息，研究区大气中SO2、NO2和总悬浮颗粒物的日均值长期优于《环境空气质量标准》GB3095-1996中一级标准，据此推断研究区降雨污染物输入比与其他来源污染物输入少得多。

2.1.2 输出系数

输出系数的合理取值是准确估算面源污染物负荷的关键。根据公式（1）农村面源污染负荷主要从农业用地、畜禽养殖和农村生活污水及废弃物3个方面来考虑，其中，畜禽养殖主要考虑牛、猪、羊和家禽。

由于宝兴县不同土地利用的输出系数尚无相关调研工作或试验依据支撑，故参考国内相似自然条件下其他地区的研究如刘瑞民等[17]在长江上游、Fu等[18]在宝兴县、梁常德等[19]在长江干流寸滩断面至长江源头的研究结果，确定本研究区不同农业用地类型的TN、TP污染物输出系数取值（表1）。

表1 不同土地利用类型的输出系数

畜禽养殖的面源污染物输出系数参照第一次全国农业污染源普查确定的西南区养殖专业户的畜禽不同饲养阶段的排污系数结合在研究区的实地调查而确定，见表2。畜禽饲养周期参考国家环保部公布的数据以及实地调查确定：牛、羊取365 d；猪取300 d；家禽取210 d。因畜禽养殖污染物排放的具体位置难以确定，本研究按照面积进行均摊，畜禽养殖可分为圈养和放养，故结合乡镇畜禽养殖统计数据、乡镇边界以及土地利用数据，将畜禽养殖的面源污染输出分摊在居民地和可放牧的草地（该县高海拔天然牧草地因参加天保工程而被保护，禁止放牧）上。

表2 畜禽养殖的排污系数Table 2 Pollutants discharge coefficients of livestock and poultry breeding g·(头·d)-1

注：羊的排污系数按3只羊等于1只猪来折算。家禽排污系数以肉鸡为例。

Note: Pollutants discharge coefficients of three sheep are equal to the coefficients of one pig. Pollutants discharge coefficients of chicken represent the coefficients of poultry.

农村生活方面的面源污染物输出系数通过如下方法确定：根据《四川省用水定额（修订稿）》（川水发[2010]4号），农村居民与城镇居民生活用水量分别为90 L/(人·d)、100 L/(人·d)，参考袁晓燕等[20]研究，村镇生活污水中TN、TP的质量浓度分别为118.9、9.75 mg/L，选取排水系数0.60，得到研究区农村人口与城镇人口的TN输出系数分别为2.34、2.60 kg/a，TP输出系数分别为0.19、0.21 kg/a。根据村尺度人口统计数据、村界和土地利用数据将农村生活方面的面源污染输出分摊在居民地上。

2.2 InVEST模型

本研究使用的模型版本为InVEST 2.4.1，主要使用“产水”和“水质净化”2个模块。

2.2.1 产水模块

InVEST模型产水模块根据水量平衡原理，基于气候、地形和土地利用来计算流域每个栅格的径流量。产水量为区域上每个栅格单元的降雨量减去实际蒸发量，而降雨量与蒸发量之间的平衡与其他一系列的气象要素、土壤特征和地表覆盖（土地利用类型或植被覆盖类型）等密切相关。该模块已经在阿根廷南部丘布特河流域[21]，美国俄勒冈州西北部的图拉丁河流域[22]，中国三江源区[23]、都江堰市[24]、密云水库流域[25]等区域进行了应用，模块计算原理如下[26]

式中Y为第土地利用/覆被类型栅格的年产水量，mm；P为栅格单元的年均降雨量，mm；AET为第土地利用/覆被类型上栅格的实际年均蒸散发量，mm；该模块需要的数据主要包括年均降雨量、潜在蒸散量、土地利用/覆被图、土壤深度、最大根系深度、植被可利用水和蒸散系数等。其中，潜在蒸散量的获得根据InVEST模型推荐的数据难以获取的地区使用Modified-Hargreaves方法[27]。

2.2.2 水质净化模块

植物和土壤可以通过储存、转换等方式移除或减少径流中的营养盐污染物以达到净化水质的作用。InVEST模型中的水质净化模块即是用于评估生态系统中植被和土壤的水质净化功能以及污染物输出情况。计算公式如下[26]

ALV=HHSx·pol（3）

（4）

式中ALV是栅格的调整后输出值，HSS是栅格的水文敏感性得分，pol是栅格的输出系数。λ是径流指数，λ表示集水区平均径流指数。是流向栅格的所有栅格产水量之和（也包括栅格自身的产水量）。

计算出每一栅格的污染物调整输出值后，模型根据坡度决定水流的路径以及土地覆被对污染物的净化能力，计算污染物随地表径流运移过程中下游各栅格对污染物的截留量以及各栅格的污染物负荷，详细计算过程见InVEST模型用户手册[26]。

该模块需要的数据主要包括DEM、土地利用/覆被、土地利用/覆被TN、TP输出系数、土地利用/覆被污染物净化效率、产水量数据等。其中，产水量数据由产水模块的结果提供，不同土地利用/覆被类型对TN、TP的净化效率取值主要参考模型用户手册、相关文献[18]以及研究区实际情况确定（表3）。

表3 不同土地利用/覆被污染物净化效率

2.3 面源污染物超标量估算

本研究结合水环境功能区划规定的水质标准，更加直观地反映面源污染的严重程度及其空间分布。面源污染物超标量计算公式如下

式中M为第种面源污染物超标量，g；e为每个栅格上污染物的负荷，mg；y为每个栅格上的产水量，L；Cs为水质类别执行水质类别标准中第种污染物浓度阈值，mg/L。

按照《四川省水功能区划》、《成都市地面水水域环境功能类别划分管理规定》，宝兴县地表水功能区水质类别执行III类水域标准。用污染物输出浓度与《地表水环境质量标准》（GB3838-2002）III类水水质标准（TN浓度≤1.0mg/L，TP浓度≤0.2 mg/L）相比较，得到TN、TP超标区域与超标量。

2.4 数据来源

2010年宝兴县55个村人口数据来源于宝兴县统计局，各乡镇畜禽养殖情况来源于《2010年宝兴县农业年报》；宝兴县DEM（30 m×30 m）数据来自全国生态环境十年变化遥感调查项目；宝兴县土地利用图来自第二次全国土地调查项目；土壤颗粒组成（黏粒、粉粒、砂粒）、土壤深度和土壤有机质数据来自第二次全国土壤普查项目；气象数据（1981－2010年日最高温、日最低温、日降雨量、太阳辐射量）通过中国气象数据共享网下载和四川省气象站购买获得。此外还有宝兴县水文站、环境监测站实测数据用于模型校验。

3 结果与分析

3.1 产水量

图2显示了宝兴县产水量的评估结果，全县多年平均产水深为852.89 mm，多年平均产水量为26.33亿m3。与宝兴县水文站多年径流监测结果28亿m3相比较，模拟产水量的相对误差为5.96%。宝兴县产水量的空间格局是县域东北部与西南部高于县域西北部和东南部，这主要与降水的空间分布有关，同时也受到复杂地形的影响。

3.2 面源污染物TN、TP负荷

宝兴县农村面源污染TN、TP负荷计算结果如表4所示。2010年宝兴县TN和TP负荷分别为1 156.44、81.41 t。根据宝兴县多年实际监测数据（2009－2015年出界水质监测断面多年平均TN浓度为0.434 mg/L，TP浓度为0.031 mg/L，宝兴县多年平均径流量为28亿m3）计算出的研究区TN、TP负荷分别为1 215.20、86.80 t（研究区点源污染较少，故在本研究中忽略不计）。模型模拟结果与实际监测结果较为接近，模型模拟TN、TP负荷相对误差分别为4.84%、6.21%，表明本研究选取的输出系数合理可靠，符合模型的输出要求。

不同污染源类型产生的TN负荷以及对TN负荷的贡献率的大小顺序为：畜禽养殖>农业用地>农村生活。TN负荷的污染源主要是畜禽养殖，而由于全县人口较少，农村生活对TN输出贡献率仅为12.08%。污染源类型对TP负荷贡献率的影响与对TN污染负荷的影响规律相同，但畜禽养殖和农村生活对TP负荷的贡献比例有所增加。

表4 2010年宝兴县面源污染物TN、TP负荷

宝兴县面源污染物TN、TP负荷的空间格局如图3所示，从图中可以看出TN、TP负荷分布较为一致，整体上呈现出分布不均，局部集中，靠近水体的规律，与耿润哲等[11]在密云水库流域的研究结果一致。

3.3 面源污染物TN、TP超标量及其空间分布

2010年，宝兴县地表水功能区水质类别执行III类水域标准。根据III类水质标准以及宝兴县多年平均产水量（28亿m3）计算出宝兴县TN、TP负荷阈值分别为2 800、560 t。本研究估算的宝兴县TN、TP负荷分别为1 156.44、81.41 t，通过与III类水质污染物负荷阈值对比可得：从流域尺度上看，宝兴县面源污染物TN、TP负荷不存在超标情况。而从栅格尺度（30 m×30 m）上看，根据每个栅格面源污染物负荷、产水量，以及III类水质标准，使用公式（6）计算出宝兴县面源污染物超标量。宝兴县TN超标量为763.63 t，占TN负荷的66.03%；TP超标量为51.16 t、占TP负荷的62.84%。

山区面源污染物TN、TP超标量存在明显的空间差异（图4）。从栅格尺度看，宝兴县TN超标范围在0.000 36～1 772.4 kg/hm2之间，平均值为69.20 kg/hm2，标准差为103.84 kg/hm2。宝兴县TP超标范围在0.003 6～144.52 kg/hm2之间，平均值为10.21 kg/hm2，标准差为8.68 kg/hm2。宝兴县面源污染物TN、TP超标的空间格局较为相似，面源污染物超标区域主要分布在耕地、草地、居民地以及靠近水体的区域。

4 讨 论

宝兴县不同污染源类型对农村面源污染贡献率差异较大。其中，农村生活对面源污染贡献率最小，畜禽养殖对面源污染物TN、TP负荷的贡献率最大，因此为减少面源污染，应对畜禽粪尿进行重点治理，防止粪尿中的污染物进入水体。畜禽养殖污染治理包括建设沼气池、田间池、专业户粪污综合利用设施[28]。此外，农业用地产生的面源污染所占比例也较大，应优化土地利用格局，合理利用土地，提高肥料利用率，推行少耕或免耕技术。

山区农村面源污染负荷以及超标量均存在明显的空间差异（图3和图4），其中农村面源污染负荷的空间差异取决于污染源污染强度，土地利用状况、植被、河道、降雨、地形、人为管理等因素（在本文中没有考虑）和自然环境因素的空间不均匀性[29]。而面源污染超标的空间差异不仅取决于污染物负荷的大小，还取决于受纳水体的水质保护目标和纳污能力的空间差异。因此，在进行山区农村面源污染治理和预防时，若忽略其空间异质性，采用针对平原地区的治理方式，会造成治理投入针对性不强，治理效率低等情况。已有研究指出，在面积较大的流域，明确面源污染负荷较高的区域作为关键风险区是非常重要和必要的工作[30-31]。而本研究认为，在明确面源污染负荷较高的区域的基础上，再根据受纳水体的水质保护目标和纳污能力阈值识别出面源污染超标区域作为治理的优先区，能提高山区农村面源污染的治理效率，为污染控制和流域管理提供依据。

5 结 论

本研究通过输出系数模型和InVEST模型水质净化模块，在实地调研的基础上，参考已有研究成果，结合降雨、地形地貌、植被与土壤等自然因素估算宝兴县农村面源污染物TN、TP负荷，并根据水环境功能分区，分析宝兴县农村面源污染超标情况。得到结论如下：

1）宝兴县农村面源污染物TN、TP负荷分别为1 156.44、81.41 t/a，TN、TP负荷的空间格局较为一致，整体上呈现出分布不均，局部集中，靠近水体的规律。

2）不同污染源对TN和TP负荷的贡献率大小顺序均为：畜禽养殖>农业用地>农村生活。对于宝兴县而言，畜禽养殖是最主要的污染源。因此，为减少面源污染，应对畜禽粪尿进行重点治理，防止粪尿中的污染物进入水体。此外，合理利用农业用地并且减少化肥的使用仍需引起重视。

3）从流域尺度上，宝兴县TN、TP负荷均小于III类水质污染物负荷阈值，其面源污染负荷不存在超标情况。而从栅格尺度（30 m×30 m）上看，根据栅格尺度上污染物负荷、产水量以及III类水质标准计算出宝兴县TN、TP超标量分别为763.63、51.16 t/a，并且TN、TP超标量存在明显的空间差异，其超标区域主要分布在耕地、草地、居民地以及靠近水体的区域，以上超标区域即为面源污染控制和治理的关键区域。

本研究揭示山区面源污染物TN、TP负荷的空间异质性，找出面源污染超标区域，为农村面源污染防治规划、山区环境治理效率的提升提供科学依据。但本研究中输出系数取值以及InVEST模型中参数设置主要通过查阅文献与参考模型手册推荐数据获得，在今后的研究中应结合现场调研与试验，获取更为可靠的输出系数与模型参数，以提高研究准确性与合理性。

[1] 朱兆良，Norse D，孙波.中国农业面源污染控制对策[M]. 北京：中国环境科学出版社，2006.

[2] 张维理，武淑霞，冀宏杰，等.中国农业面源污染形势估计及控制对策I. 21世纪初期中国农业面源污染的形势估计[J]. 中国农业科学，2004，37(7)：1008－1017.

Zhang Weili, Wu Shuxia, Ji Hongjie, et al. Estimation of agricultural non-point source pollution in China and the alleviating strategies I. estimation of agricultural non-point source pollution in China in early 21 century[J]. Scientia Agricultura Sinica, 2004, 37(7): 1008－1017. (in Chinese with English abstract)

[3] 刘亚琼，杨玉林，李法虎.基于输出系数模型的北京地区农业面源污染负荷估算[J]. 农业工程学报，2011，27(7)：7－12.

Liu Yaqiong, Yang Yulin, Li Fahu. Estimation of pollution loads from agricultural nonpoint sources in Beijing region based on export coefficient modeling approach[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 7－12. (in Chinese with English abstract)

[4] 龙天渝，梁常德，李继承，等.基于SLURP模型和输出系数法的三峡库区非点源氮磷负荷预测[J]. 环境科学学报，2008，28(3)：574－581.

Long Tianyu, Liang Changde, Li Jicheng, et al. Forecasting the pollution load of non-point sources imported to the Three Gorges Reservoir[J]. Acta Scientiae Circumstantiae, 2008, 28(3): 574－581. (in Chinese with English abstract)

[5] Worrall F, Burt T. The impact of land-use change on water quality at the catchment scale: The use of export coefficient and structural models[J]. Journal of Hydrology, 1999, 221(1): 75－90.

[6] Ierodiaconou D, Laurenson L, Leblanc M, et al. The consequences of land use change on nutrient exports: A regional scale assessment in south-west Victoria, Australia[J]. Journal of Environmental Management, 2005, 74(4): 305－316.

[7] 薛利红，杨林章. 面源污染物输出系数模型的研究进展[J]. 生态学杂志，2009，28(4)：755－761.

Xue Lihong, Yang Linzhang. Research advances of export coefficient model for non-point source pollution[J]. Chinese Journal of Ecology, 2009, 28(4): 755－761. (in Chinese with English abstract)

[8] Noto L V, Ivanov V Y, Bras R L,et al. Effects of initialization on response of a fully-distributed hydrologic model[J]. Journal of Hydrology, 2008, 352(1): 107－125.

[9] 蔡明，李怀恩，庄咏涛，等.改进的输出系数法在流域非点源污染负荷估算中的应用[J]. 水利学报，2004(7)：40－45.

Cai Ming, Li Huaien, Zhuang Yongtao, et al. Application of modified export coefficient method in polluting load estimation of non-point source pollution[J]. Journal of Hydraulic Engineering, 2004(7): 40－45. (in Chinese with English abstract)

[10] Ding X, Shen Z, Hong Q, et al. Development and test of the export coefficient model in the upper reach of the Yangtze River[J]. Journal of Hydrology, 2010, 383(3): 233－244.

[11] 耿润哲，王晓燕，焦帅，等. 密云水库流域非点源污染负荷估算及特征分析[J]. 环境科学学报，2013，33(5)：1484－1492.

Geng Runzhe, Wang Xiaoyan, Jiao Shuai, et al. Application of improved export coefficient model in estimating non-point source nutrient load from Miyun reservoir watersheds[J]. Acta Scientiae Circumstantiae, 2013, 33(5): 1484－1492. (in Chinese with English abstract)

[12] Kareiva P, Tallis H, Ricketts T H, et al. Natural capital: theory and practice of mapping ecosystem services[M]. Oxford: Oxford University Press, 2011.

[13] Chen L, Wei G, Zhong Y, et al. Targeting priority management areas for multiple pollutants from non-point sources[J]. Journal of hazardous materials, 2014, 280: 244－251.

[14] 严素定. 黄石市农业面源污染的解析及其空间异质性研究[J]. 农业工程学报，2008，24(9)：225－228.

Yan Suding. Source apportionment and spatial heterogeneity of agricultural non-point source pollution in Huangshi, Hubei Province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(9): 225－228. (in Chinese with English abstract)

[15] 钱晓雍，沈根祥，郭春霞，等.基于水环境功能区划的农业面源污染源解析及其空间异质性[J]. 农业工程学报，2011，27(2)：103－108.

Qian Xiaoyong, Shen Genxiang, Guo Chunxia, et al. Source apportionment and spatial heterogeneity of agricultural non-point source pollution based on water environmental function zoning[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(2): 103－108. (in Chinese with English abstract)

[16] Johnes P J, O'sullivan P. The natural history of slapton ley nature reserve. XVIII, nitrogen and phosphorus losses from the catchment–An export coefficient approach[J]. Field Studies, 1989, 7(2): 285－309.

[17] 刘瑞民，杨志峰，丁晓雯，等.土地利用/覆盖变化对长江上游非点源污染影响研究[J]. 环境科学，2006，27(12)：2407－2414.

Liu Ruimin, Yang Zhifeng, Ding Xiaowen, et al. Effect of land use/cover change on pollution load of non-point source in upper reach of Yangtze river basin[J]. Environmental Science, 2006, 27(12): 2407－2414. (in Chinese with English abstract)

[18] Fu B, Wang Y, Xu P, et al. Modelling nutrient retention function of ecosystem: A case study in Baoxing County, China[J]. Procedia Environmental Sciences, 2012, 13(1): 111－121.

[19] 梁常德，龙天渝，李继承，等.三峡库区非点源氮磷负荷研究[J]. 长江流域资源与环境，2007，16(1)：26－30.

Liang Changde, Long Tianyu, Li Jicheng, et al. Importation loads of non-point source nitrogen and phosphorus in the Three Gorges Reservoir[J]. Resources and Environment in the Yangtze Basin, 2007, 16(1): 26－30. (in Chinese with English abstract)

[20] 袁晓燕，余志敏，施卫明.大清河流域典型村镇生活污水排放规律和污染负荷研究[J]. 农业环境科学学报，2010，29(8)：1547－1557.

Yuan Xiaoyan, Yu Zhimin, Shi Weiming. Domestic sewage emission dynamics and pollutant loading capacity of the Daqing River Valley: A case study on the village scale[J]. Journal of Agro-Environment Science, 2010, 29(8): 1547－1557. (in Chinese with English abstract)

[21] Pessacg N, Flaherty S, Brandizi L, et al. Getting water right: A case study in water yield modelling based on precipitation data[J]. Science of the Total Environment, 2015, 537: 225－234.

[22] Hoyer R, Chang H. Assessment of freshwater ecosystem services in the Tualatin and Yamhill basins under climate change and urbanization[J]. Applied Geography, 2014, 53: 402－416.

[23] 潘韬，吴绍洪，戴尔阜，等.基于InVEST模型的三江源区生态系统水源供给服务时空变化[J]. 应用生态学报，2013，24(1)：183－189.

Pan Tao, Wu Shaohong, Dai Erfu, et al. Spatiotemporal variation of water source supply service in Three Rivers Source Area of China based on InVEST model[J]. Chinese Journal of Applied Ecology, 2013, 24(1): 183－189. (in Chinese with English abstract)

[24] 傅斌，徐佩，王玉宽，等.都江堰市水源涵养功能空间格局[J]. 生态学报，2013，33(3)：789－797.

Fu Bin, Xu Pei, Wang Yukuan, et al. Spatial pattern of water retetnion in Dujiangyan County[J]. Acta Ecologica Sinica, 2013, 33(3): 789－797. (in Chinese with English abstract)

[25] 李屹峰，罗跃初，刘纲，等.土地利用变化对生态系统服务功能的影响：以密云水库流域为例[J]. 生态学报，2013，33(3)：726－736.

Li Yifeng, Luo Yuechu, Liu Gang, et al.Effects of land use change on ecosystem services: a case study in Miyun reservoir watershed[J]. Acta Ecologica Sinica, 2013, 33(3): 726－736. (in Chinese with English abstract)

[26] Tallis H T, Ricketts T, Guerry A D, et al. InVEST 2.4.1 User’s Guide[Z].The Natural Capital Project, Stanford, 2011.

[27] Droogers P, Allen R G. Estimating reference evapotranspiration under inaccurate data conditions[J]. Irrigation and drainage Systems, 2002, 16(1): 33－45.

[28] 马亚丽，敖天其，张洪波，等.基于输出系数模型濑溪河流域泸县段面源分析[J]. 四川农业大学学报，2013(1)：53－59.

Ma Yali, Ao Tianqi, Zhang Hongbo, et al. Non-point source analysis of Laixi River Basin in Luxian County based on output coefficient model[J]. Journal of Sichuan Agricultural University, 2013(1): 53－59. (in Chinese with English abstract)

[29] Lai Y C, Yang C P, Hsieh C Y, et al. Evaluation of non-point source pollution and river water quality using a multimedia two-model system[J]. Journal of Hydrology, 2011, 409(3): 583－595.

[30] Shen Z Y, Hong Q, Chu Z, et al. A framework for priority non-point source area identification and load estimation integrated with APPI and PLOAD model in Fujiang Watershed, China[J]. Agricultural Water Management, 2011, 98(6): 977－989.

[31] 李思思，张亮，杜耘，等.面源磷负荷改进输出系数模型及其应用[J]. 长江流域资源与环境，2014，23(9)：1330－1336.

Li Sisi, Zhang Liang, Du Yun, et al. Estimation of non-point phosphorus load using improved export coefficient model and its application in watershed management[J]. Resources and Environment in the Yangtze Basin, 2014, 23(9): 1330－1336. (in Chinese with English abstract)

Estimation and evaluation of non-point source pollutants loads in mountainous area of Sichuan

Liu Ju1,2, Fu Bin1,3, Wang Yukuan1,3※, Xu Pei1,3

(1.610041,; 2.100049,; 3.404020,)

Non-point source (NPS) pollution produced from agricultural production activity is one of the largest threats to water environment. The management and control of NPS pollutants to achieve the goals of improving surface water quality have gotten an increasing concern throughout the world. Because of the heterogeneity of key driving factors in NPS pollutions such as terrain, precipitation, and vegetation type, the management of NPS pollution in mountainous areas needs more meticulous effort than the plain areas. In this research, considering the characteristics of the spatial heterogeneity of the natural and human factors in mountainous areas, the loads of NPS pollutants, such as total nitrogen (TN) and total phosphorus (TP) are evaluated by the export coefficient model and the InVEST (integrate valuation of ecosystem services and tradeoffs tool), which is a water purification model. The excesses of NPS pollutants are analyzed according to the water environmental function zoning. The research is intended to reveal the spatial difference of NPS pollution in mountainous areas, determine the key regions for management, and provides a foundation for regional pollution control planning and environmental management efficiency promotion in mountainous areas. The method to estimate and evaluate NPS pollutants loads is developed based on the export coefficient model and InVEST water purification model, and is successfully applied to the case study of Baoxing County, showing the applicability of this method for mountainous areas. The results indicate that: 1) TN and TP loads in Baoxing County are respectively 1 156.44 and 81.41 t/a in 2010, and the spatial distribution of TN and TP loads has the characteristics of uneven distribution, local concentration, and closing to the water. 2) The contributions to pollution from various types of pollutant sources are different. The magnitude order for the contribution from various types of pollutant sources is livestock and poultry breeding > agricultural land > rural life. As above mentioned, the contribution of the livestock and poultry breeding to NPS pollution is the biggest and sustainable resource management should be prioritized in an effort to reduce the burden. Such management practices may involve building methane tanks, field ponds and the comprehensive utilization of waste facilities, which are aimed to reduce the level of NPS pollutant export into water bodies. 3) Compared with the pollutants thresholds of water quality standard at Class III at the watershed scale, the total amounts of TN and TP loads in Baoxing County do not exceed the thresholds. However, at the 30 m × 30 m grid scale, according to the pollutant load, water yield and water quality standard at Class III, TN and TP excesses in Baoxing County are calculated, and the results are respectively 763.63 and 51.16 t/a. The excesses of NPS pollutants in mountainous areas show significant spatial variation. The regions in which TN and TP loads exceed the thresholds are mainly distributed in cultivated land, grassland, residential land and the regions near the water, and these regions are the key regions for NPS pollution control and management.

agriculture; pollution; models; water environmental function zoning; water yield; spatial variation

10.11975/j.issn.1002-6819.2016.24.029

X522

A

1002-6819(2016)-24-0218-08

2016-04-25

2016-06-11

国家自然科学基金资助项目（41371539）；中国科学院“西部之光”人才培养计划项目（Y2R2180180）；中-挪生物多样性与气候变化项目四川省分项目（C/IV/S/11/242-02）

刘菊，博士生，主要从事生态补偿、流域生态与管理方面的研究。成都 中科院成都山地灾害与环境研究所610041。 Email：liuju0607@163.com

王玉宽，研究员，博士生导师，主要从事流域生态与管理、水土保持方面的研究。成都 中科院成都山地灾害与环境研究所610041。Email：wangyukuan@imde.ac.cn