基于SWAT模型的典型农业小流域氮污染时空分布特征及关键源解析*

王慧勇，遆超普，王良杰，董 磊，夏永秋，颜晓元

(1:南京林业大学南方现代林业协同创新中心，南京 210037) (2:中国科学院南京土壤研究所，土壤与农业可持续发展国家重点实验室，江苏常熟农田生态系统国家野外科学观测研究站，南京 210008)

氮是植物生长必要的营养元素，但过量的氮进入水体会造成氮污染，使得水体生态系统恶化. 例如，美国墨西哥湾流域和切萨皮克湾流域自1950年来由于氮肥的大量使用，当地的海湾水体富营养化不断加剧，破坏了其生态系统健康[1-2]. 2018年我国对1097座水库的营养状况评价结果显示, 69.6%的水库处于中营养状态，30.4%的水库处于富营养状态，主要由氮磷等营养物质过度排放所造成的[3]. 水体污染的形式分为点源污染和非点源污染，随着点源污染的有效控制，非点源污染已成为我国水体污染的主要因素，其中农业非点源氮污染问题尤为严重. 根据我国第二次全国污染源普查公报，农业活动造成的总氮排放占非点源污染排放总量的46.52%. 同时在太湖、巢湖等富营养化重点区域的研究中，由农业活动所引起的总氮污染高达60%[4-5].

模型模拟可有效应用于对非点源污染的定量研究，从而达到控制和减少非点源污染的目的，其中SWAT模型应用最为广泛[6]. SWAT模型由美国农业部农业研究中心(USDA)于1994年研发，其主要目的是在不同土地管理方式以及不同土地利用等条件下对流域内土壤的状况进行相关预测. 后经众多学者的研究改进，模型对于流域产流[7]、产沙[8]、营养物质[9]迁移的变化规律都具有较好的模拟效果，多用于大、中尺度的流域[10-11]，在小尺度流域也有较好的模拟效果[12]. 周铮等[13]运用SWAT模型对饮用水保护区北山水库小流域进行径流模拟，在实测径流数据缺乏情况下采用水库参数进行率定和验证，结果显示北山水库蓄水量的模拟值与实测值吻合性较好，表明模型适用于该区域小流域的研究. 林晓娟等[14]利用SWAT模型对岳城水库流域的污染源进行模拟，得出该流域年入库非点源污染所占比例近50%，其中总氮含量高为3025 t/a，其污染的主要来源为耕地. Bauwe等[15]基于农业小流域评估SWAT模型对径流、硝态氮负荷和作物产量三者分别模拟后的效果，三者基本符合模拟要求. 在此基础上设置相应情景模拟进行分析，结果表明在保持作物产量的情况下适当减少氮肥施用量可使得硝态氮负荷显著减少. Panagopoulos等[16]利用SWAT模型对Arachtos流域的氮负荷量进行模拟，确定了氮流失的重点区域. 韩博等[17]在蒲河流域建立非点源污染模型，分别对氮污染含量、氮污染的时空分布和关键区进行了模拟研究.

句容水库农业流域是位于江苏省句容市的典型农业小流域，流域内的句容水库为省重点中型水库，是市饮用水二级保护区和生态保护区，其水质情况直接影响当地居民的饮用水安全. 因此以句容水库农业流域为研究区域，通过构建SWAT非点源污染模型，以此分析氮污染的时空分布特点，探讨流域氮污染的关键源区，可为当地流域氮污染控制提供理论依据.

1 材料和方法

1.1 研究区概况

研究区句容水库农业流域(图1)位于江苏省句容市东北部丘陵山区的句容河上游(32°01′N, 119°13′E)，东近镇江，西临南京市，距南京大约40 km，离句容城区仅4 km，属长江水系秦淮河流域，流域集水面积45.5 km2. 研究区属北亚热带季风气候，气候湿润，多年平均气温15℃，年平均相对湿度为78%. 多年平均降雨量为1050 mm，降雨主要集中在6-9月. 流域东起句东林场，西至句容水库，地形总体上为四周高中部低，海拔在3.00～279.13 m之间，平均海拔低于50 m. 句容水库总库容2907万m3，兴利库容1202万m3，防洪库容1924万m3，水库下游为句容市华阳街道，防洪保护人口11.5万人，灌溉面积2061.3 hm2.

研究区内雨水径流自然流动或经由人工简易沟道排入就近水体，无专门修筑输水管道或渠道，且研究区为仅有一个排水出口的封闭性良好的小流域. 研究区内无大型养殖业和工业，生产活动以农业为主，主要种植方式为水稻-小麦轮作，流域内的污染主要以非点源形式存在.

图1 句容水库农业流域地理位置Fig.1 Location of Jurong Reservoir agricultural watershed

1.2 数据来源及处理

构建研究区SWAT模型所需基础数据库包括空间数据库和属性数据库，空间数据库主要包括研究区数字高程模型(DEM)、土地利用类型图、土壤类型图、河网水系图；属性数据库包括气象数据、土壤理化属性数据、污染源数据以及水文水质数据等数据. 根据SWAT模型叠加分析时的统一图层要求，研究中图层均采用Krasovsky_1940_Transverse_Mercator投影坐标系.

1.2.1 数字高程模型 DEM作为提取流域坡长、坡度等参数的必要数据，在模型中用于河网的生成和子流域的划分. DEM的分辨率对模型的模拟结果具有一定的影响，选择合适的分辨率可以提高模拟精度[18]. 经过对比，研究采用来源于中国科学院南京土壤研究所的10 m分辨率DEM，依据流域边界图进行裁剪得到句容水库农业流域DEM数据，见图2.

图2 句容水库农业流域基础数据Fig.2 Basic data of Jurong Reservoir agricultural watershed

1.2.2 土壤类型 土壤类型采用来源于中国科学院南京土壤研究所1∶5万比例尺的栅格图，结合已有数据并查阅《江苏土壤》和中国土壤数据库[19]，通过SPAW等相关软件计算出土壤属性数据，流域内土壤类型共有9种，见图2.

1.2.3 土地利用类型图和河网水系图 土地利用类型是水文循环的主要影响因素之一，是SWAT模型的必要输入数据. 土地利用类型基于2010年遥感影像进行人工目视解译所得，同时根据SWAT模型要求的土地利用分类体系，土地利用类型分为居民区、水体、林地、裸地和耕地(共5种)，详见图2.

SWAT模型应用于平原区时存在子流域划分边界不准确的情况[20]，研究采用模型中的Burn In功能结合实测河网进行修正，以改进提取的地形地貌特征[21]. 实测河网水系图通过DEM划分并结合实地勘察获得.

1.2.4 气象数据 模型构建所需的气象数据分为实测数据和天气发生器两部分，输入的实测数据包括日降水量、日最低最高气温、日平均风速、太阳辐射和相对湿度，天气发生器用于补全实测数据中的部分缺失数据. 鉴于研究区内无气象站点，采用临近站点南京气象站2000—2019年20年长时间序列逐日气象数据计算天气发生器各参数，气象站点数据来源于中国气象科学数据网. 同时采用由句容水库管理所所提供的句容水库农业流域2007-01-2012-09逐日降水和日最低、最高气温数据输入模型，其余数据由天气发生器生成.

1.2.5 水文水质数据 水文水质数据主要用于模型的校正. 水文数据主要为水库数据，其中包括水库基本参数和水库管理数据. 基本参数包括正常水位对应的水域面积和蓄水量、初始蓄水量、防洪水位对应的水域面积和蓄水量以及蒸发等相关参数。水库管理数据主要为水库入水量数据，水库入水量数据为2007-01-2012-09逐月数据. 水质数据为中间河流上水质观测点的逐月观测数据，观测时间长度为2007-01-2012-09.

1.2.6 污染源数据 根据2010年实地调研结果，获取流域内常住人口、家禽养殖、氮肥施用数据. 结合第二次全国污染普查产排污系数计算出各源头的污染量，再转化为化肥以连续施肥的方式施用到相应子流域的耕地中进行模拟计算. 通过实验观测得句容水库农业流域氮的平均沉降为39 kg /( hm2·a). 流域内水稻和小麦均在播种前施基肥，在生长发育阶段进行2次追肥，肥料种类为尿素、碳铵和复合肥. 经换算成N肥用量后，水稻和小麦基肥施肥量分别为110.05和73.37 kg/hm2，水稻在7月上旬和8月上旬进行追肥，两次施肥量都为82.54 kg/hm2；小麦在1月上旬和3月上旬的追肥量均为55.03 kg/hm2.

1.3 研究方法

1.3.1 模型建立 SAWT模型内部采用模块化结构，由水文过程子模型、土壤侵蚀子模型和污染负荷子模型3个子模型所构成[22]，可模拟流域内径流、泥沙、营养物质等各物质的输移过程. SWAT模型主要依据水量平衡方程进行，该方程是流域内所有过程的驱动力，水量平衡方程如下：

(1)

式中，SWt代表i天土壤最终含水量，mm；SW0代表土壤初始含水量，mm；t代表时间，d；Rday、Qsurf、Ea、Wseep、Qgw依次代表第i天的降水、地表径流、蒸发、土壤层的渗透和侧流、地下含水量，单位均为mm.

基于DEM和实测河网水系图模型划分子流域，Jha等[23]的研究结果表明子流域划分阈值在小于总面积的5%时，氮磷模拟能够得到较为合理的值，因此本研究基于此条件将句容水库农业流域按照最优标准共划分25个子流域. 基于划分的子流域， SWAT模型结合土地利用类型、土壤类型和坡度(阈值均为10%)进一步将子流域细分为具有相同土地利用、土壤类型、坡度的水文响应单元(HRUs)，共155个HRUs.

1.3.2 率定和验证 为使模拟的水文过程和物质输移特征更贴近流域内的实际情况，需要利用流域内实测数据进行率定和验证，对模型参数值做出合理调整. 模型的率定和验证遵循先径流后营养物质的原则，采用句容水库入库径流量对径流进行率定和验证，对氮营养物质的率定验证采用逐月实测河流水质数据，时间长度均为2007-01-2012-09. 模型设置2007年为预热期，2008-01-2010-12为率定期，验证期为2011-01-2012-09.

使用SWAT-CUP2012软件对模型进行率定和验证，软件集成SUFI2、PSO和GLUE等多种算法，研究选用较为常用的算法SUFI2. 采用决定系数(R2)和纳什系数(NSE)作为模拟结果的评判标准，在基于月尺度的模拟中当R2≥0.6和NSE≥0.5时认为模拟结果是可行的.

SWAT模型中参数众多，为减少参数率定的复杂性，常通过参数敏感性分析选出对结果影响显著的参数进行率定. 研究选取对径流和水质有重要影响的20个参数作为参数率定的重点，分别对径流和水质进行率定和验证，其率定参数见表1.

表1 率定参数敏感性及最优结果*Tab.1 Sensitivity of calibration parameters and optimal results

表2 模型评估结果Tab.2 Model evaluation results

2 结果与分析

2.1 率定和验证

图3 句容水库农业流域径流量模拟值与实测值比较Fig.3 Comparison of simulated and measured values of runoff(a), TN(c) in Jurong Reservoir agricultural watershed

图4 句容水库农业流域氮素和径流量的年内(a)和年际变化(b)Fig.4 Intra-year(a) and inter-year(b) variation of nitrogen and runoff in the Jurong Reservoir agricultural watershed

2.2 入库氮素时间分布特征

2.3 流域径流氮素负荷空间分布特征

图5 句容水库农业流域氮素单位负荷量的空间分布Fig.5 Spatial distribution of unit nitrogen load in the Jurong Reservoir agricultural watershed

图6 不同污染源对入库TN的贡献率Fig.6 Contribution of different pollution sources to TN load in the Jurong Reservoir agricultural watershed

2.4 来源分析

句容水库农业流域现状氮素主要来源于农田施肥、农村生活污水、畜禽养殖、氮沉降. 在农业非点源污染中土壤背景值对氮素污染的贡献率较大，不可忽略，因此在句容水库农业流域来源分析中考虑流域的土壤背景值. 通过采用分别模拟有、无相应污染源输入时入库水质变化情况的方法，估算出句容水库农业流域内不同污染源对入库TN的污染贡献率.

由图6可知，句容水库入库TN负荷贡献最大的为农田施肥，贡献率高达49.60%；背景值的贡献率为21.19%；氮沉降18.21%；农村生活污水次之，贡献率为8.32%；畜禽养殖最小. 因此各输入污染源对入库TN的贡献率为：农田施肥>氮沉降>农村生活污水>畜禽养殖.

3 讨论

3.1 模型率定和验证

3.2 氮污染时空分布特征

句容水库农业流域属于北亚热带季风气候，年内降水分布集中，汛期的降水可占到全年降水的60%以上，丰富的降水汇流形成径流，氮素在受到降水径流的影响同样主要集中在汛期，造成汛期是句容水库农业流域在一年内氮素污染最大的时期. 张婷等[27]针对潘家口水库流域的非点源污染研究得出，TN负荷在汛期占全年的72%，与本文结论相似.

在模拟的5年内年径流氮负荷量最大的一年为2009年，其降水量和径流量同样最大，主要原因是降水和径流是非点源污染的主要驱动力，当下垫面情况稳定时，非点源污染主要受这两者的影响[28]. 基于句容水库农业流域的模拟结果，降水量和径流量在对流域内氮素输移的影响程度上略有差别，径流氮负荷量与径流的相关性较与降水量的相关性好. 究其原因，降水虽然能够影响径流量从而影响氮素迁移，但是降水仅能够反映氮素从土壤剥离的量，无法体现输出氮的迁移量；径流不仅有降水的影响，同时包括对沿途土壤前期含水中的氮素的作用[27]. 因此在汛期减少如农田施肥等氮源输入可有效控制流域氮素污染.

3.3 氮污染关键来源

4 结论

1)从时间变化上看，入库氮素负荷量随着径流的变化而变化. 年际入库氮素负荷量总体呈现“M”型波动趋势，其中以2009年的入库量最多；年内入库氮素负荷量随降水径流的变化而呈现出集中分布的特征，主要集中在汛期. 因此合理安排施肥时间可作为减少污染的有效措施.

3)农田施肥为流域内主要的污染源，占49.60%；由于土壤的氮储存土壤背景值贡献值达21.19%. 控制农田施肥的施用量可作为治理句容水库农业流域氮素污染的有效措施，同时土地利用类型的改变也会影响氮素迁移量，这可作为在句容水库农业流域内的进一步研究方向.