环境科学  2020, Vol. 41 Issue (8): 3646-3656   PDF    
三峡库区规模化顺坡沟垄果园氮、磷输出过程及流失负荷
严坤1,2, 王玉宽1,3, 刘勤1,3,4, 徐佩1,3, 闫洋洋1,3     
1. 中国科学院、水利部成都山地灾害与环境研究所, 成都 610041;
2. 中国科学院大学, 北京 100049;
3. 水利部三峡工程运维系统万州典型区水土保持监测站, 万州 404020;
4. 成都理工大学地球科学学院, 成都 610059
摘要: 通过对三峡库区顺坡垄沟构型的规模化柑橘园集水区次降雨过程径流氮、磷进行动态监测,分析典型降雨事件氮和磷流失负荷,并探讨了规模化柑橘种植对土壤氮磷流失及入库河流水环境的影响.结果表明:①顺坡沟垄柑橘园集水区径流氮和磷年流失负荷分别为13.43 kg·(hm2·a)-1和1.26 kg·(hm2·a)-1,春季施肥及强降雨冲刷是集水区污染物高负荷的主要原因;②集水区全年总氮(TN)和总磷(TP)的EMC为8.49 mg·L-1和0.87 mg·L-1,超过发生水体富营养化含量标准;③春季施肥后的2场典型降雨中,长历时暴雨径流硝态氮(NN)和溶解态磷(DP)负荷为4.94 kg·hm-2和0.28 kg·hm-2,分别占TN和TP流失负荷的92.90%和64.69%;短历时大雨径流NN和DP负荷为0.52 kg·hm-2和0.05 kg·hm-2,分别占TN和TP的65.92%和74.88%,溶解态氮和磷是顺坡沟垄坡面果园径流氮和磷流失的主要途径;④集水区氮磷流失表现出显著的"初始冲刷效应",初期20%的地表径流流失了58.0%的TN,57.0%的DN,58.5%的NN,79.0%的AN,62.0%的TP,63.5%的DP和60.0%的PP,控制初期地表径流对降低入库径流养分具有重要作用.
关键词: 三峡库区      顺坡沟垄      集水区      面源污染      规模化果园     
Dynamic Process of Nitrogen and Phosphorus Export and Loss Load in an Intensive Orchard with Ridge and Furrow Plantation in the Three Gorges Reservoir Area
YAN Kun1,2 , WANG Yu-kuan1,3 , LIU Qin1,3,4 , XU Pei1,3 , YAN Yang-yang1,3     
1. Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China;
3. Wanzhou Key Regional Soil and Water Conservation Monitoring Station of Three Gorges Project Operation and Maintenance System, Wanzhou 404020, China;
4. Earth Sciences Faculty, Chengdu University of Technology, Chengdu 610059, China
Abstract: To comprehend the runoff load of nitrogen (N) and phosphorus (P) and the impact on the receiving river in an agricultural area with an intensive orchard plantation and a longitudinal ridge and furrow morphology in the Three Gorges Reservoir Area, the runoff and N and P concentrations were dynamically monitored in a typical citrus orchard catchment in Wanzhou Country, Chongqing, China. The results showed that the nutrient concentration in runoff water from the intensive citrus planting catchment was very high. The average annual event mean concentrations (EMC) were 9.31 mg·L-1 for total nitrogen (TN), 8.11 mg·L-1 for dissolved nitrogen (DN), 5.66 mg·L-1 for nitrate nitrogen (NN), 0.51 mg·L-1 for ammonium nitrogen (AN), 0.87 mg·L-1 for total phosphorus, 0.56 mg·L-1 for solved phosphorus (DP), and 0.32 mg·L-1 for particulate phosphorus (DP). In addition, the annual loss loads were 13.43, 12.20, 8.77, 0.75, 1.26, 0.84, and 0.42 kg·(hm2·a)-1 for TN, DN, NN, AN, TP, DP, and PP, respectively. The annual average concentrations of TN and TP were 8.49 mg·L-1 and 0.87 mg·L-1, respectively, which exceeded the category V values of the surface water quality standards (GB3838-2002) by 4.25 times and 2.2 times, respectively, and also exceeded the internationally recognized thresholds for the eutrophication of waterbodies. The TN and TP loss load from storm runoff was one of the main reasons for the degradation of the river water quality, thus suggesting the need to treat surface runoff and control runoff nutrient losses, especially during the first storm events after fertilization. During two typical long-duration springtime rainfall events after fertilization, the loads of nitrate nitrogen (NN) and dissolved phosphorus (DP) were 4.94 kg·hm-2 and 0.28 kg·hm-2, respectively, which accounted for 92.90% and 64.69% of the total annual TN and TP loss loads, respectively. The loads of NN and DP in a short-duration high-intensity rainfall event were 0.52 and 0.05 kg·hm-2 respectively, which accounted for 65.92% and 74.88% of the total annual TN and TP loss loads, respectively. The DN and DP were the main forms of nitrogen and phosphorus losses from the intensive citrus orchard with a longitudinal ridge and furrow morphology. Meanwhile, the catchment showed a significant first-flush phenomenon during a typical rainfall event, with a total of 58.0%, 57.0%, 58.5%, 79.0%, 62.0%, 63.5%, and 60.0% of the mass of TN, DN, NN, AN, TP, DP, and PP in the initial 20% of the runoff, respectively. Hence, controlling the surface runoff at the early runoff stage plays an important role in reducing nutrient losses.
Key words: Three Gorges Reservoir Area      longitudinal ridge and furrow morphology      catchment      non-point source pollution      intensive orchard     

三峡库区是我国重要的柑橘生产适宜区, 柑橘业对库区农村经济发展和移民安稳致富发挥了重要作用[1].近年来, 快速城镇化导致耕地撂荒或闲置, 土地流转和三峡库区生态屏障建设加快了柑橘种植, 逐渐形成了1 667 km2柑橘生态屏障带, 对库周土地利用格局造成重大改变.由于大部分规模化柑橘园分布在紫色土丘陵坡地上, 紫色土土层瘠薄且“上覆土壤、下伏岩石”的岩土二元结构, 土地生产力普遍较低, 规模化柑橘园在整地时沿顺坡高厢聚土起垄, 在坡面顺坡方向形成“宽垄-窄沟”相间的坡面构型, 达到快速排水与增厚土层, 形成柑橘适宜的生长条件[2].坡面是地表过程发生和发展的基础单元[3], 也是径流形成和养分迁移主要的界面, 坡面形态能显著影响地表径流和壤中流以及伴随发生的溶质运移[4~6].三峡库区中部规模化柑橘园整地形成的顺坡垄沟不同于传统意义上自然坡面、人工改造坡面(如梯田、鱼鳞坑和水平阶等)和农田耕作微地形, 具有独特的坡面水文过程和养分输移机制.目前虽有众多水文模型与养分迁移模型, 但大多数养分迁移机制模型假设条件为理想坡面, 难以模拟复杂微地形下养分流失负荷, 并不清楚“宽垄-窄沟”坡面构型的集中连片柑橘园养分流失特征及规律.

在三峡库区柑橘种植面积占农地面积的35%以上[7], 因其化肥施用量大, 氮磷流失负荷远高于传统粮食作物, 是区域重要的面源污染源[8].目前, 在三峡库区面源污染方面研究多集中于传统农作物(如玉米和红薯等)坡耕地氮素迁移特征和过程[9~11], 而柑橘面源污染主要通过地块尺度平顺径流小区监测, 分析不同农业措施(如覆盖、生草和植物篱等模式)径流、泥沙与养分流失特征[12~15], 缺乏对径流、养分流失等复杂过程的研究, 难以制定基于过程的柑橘园污染物控制措施[16].本文以三峡库区中部广泛分布的规模化柑橘园为研究对象, 连续监测顺坡垄沟坡面构型的柑橘集水区径流、氮及磷养分输出形态、浓度和负荷变化, 分析该区域规模化果园普遍采用的整地方式对水质的影响, 明确春季施肥后典型降雨事件对非点源污染贡献, 以期为三峡库区规模化果园面源污染负荷定量评估与氮磷面源污染防控提供科学依据.

1 材料与方法 1.1 研究区概况

研究区位于三峡库区中部重庆市万州区龙沙镇龙安村, 地理位置(108°10E, 30°25′N, 图 1), 属于三峡库区入库支流瀼渡河流域.

图 1 研究区位置及规模化柑橘集水区监测点示意 Fig. 1 Location of the orchard catchment

该流域所在的区域属亚热带季风湿润气候, 年均气温18.0℃, 年均降雨量1 229 mm, 年均相对湿度82%, 年均日照1 258 h, 适宜的柑橘生长气候使该区成为三峡库区重要的柑橘种植区, 近5年万州区大力发展晚熟柑橘[品种为塔罗科血橙(Tarocco)]产业, 新建、改造和更新品种面积超过80 km2, 目前规模化柑橘园主要采用顺坡沟垄方式整地(垄沟规格为:垄宽2.0~2.5 m, 垄高0.4~0.5 m, 沟宽0.5~0.8 m, 沟深0.3~0.5 m), 并在垄上按4 m株距栽种柑橘(图 2).

图 2 顺坡沟垄整地方式及汇流示意 Fig. 2 Longitudinal ridge-furrow morphology and flow concentration sketch

该区域土壤类型为紫色土, 呈微碱性或中性, 土壤理化性质如下:pH 7.23±0.29、有机质(15.2±4.45)g·kg-1、速效氮(82.03±17.5) mg·kg-1、速效磷(85.3±2.26)mg·kg-1和速效钾(285.78±20.0)mg·kg-1, 降雨后规模化柑橘园地表径流通过顺坡垄沟和排水沟道汇入瀼渡河并排入水库.

1.2 研究方法 1.2.1 降雨量与径流量监测

在0.6 hm2的柑橘坡面集水区下侧修建V型测流堰, 通过完整的坡面沟渠系统, 将地表径流汇入堰槽进行实时监测.利用自计水位仪(HOBO U20), 每5 min记录一次过堰水位, 并根据“水位-流量”关系计算径流量, 同时利用测流堰站房上安装自动降雨监测设备(HOBO RG3-M), 连续记录降雨量及降雨历时.

1.2.2 水样采集及分析

利用全自动水质采样器(ISCO 6712)采集径流水样.当降雨达到一定条件时(连续4 h累积降雨量超过10 mm), 采样器自动触发并按照设定的采样条件采集瞬时水沙样品, 前30 min时每10 min采样一次, 30~90 min时每20 min采样一次, 90~540 min时每30 min采样一次, 540~720 min时每60 min采样一次, 之后每120 min采样一次, 每瓶采集1 L, 共采集24个过程样.取样结束后, 立即将水样转移到聚乙烯瓶并放入4℃冰箱保存, 48 h内完成分析.分析指标包括总氮(TN)、溶解性总氮(DN)、硝态氮(NN)、氨氮(AN)、总磷(TP)和溶解性总磷(DP).其中, TN和TP用摇匀后的水样直接测定; DN、AN、NN和DP用经0.45 μm滤膜过滤后的水样测定.TN和TP采用碱性过硫酸钾消解-紫外分光光度法测定; NN、AN、DN和DP采用流动注射分析仪(AutoAnalyzer3)测定, PP用TP与DP差值表示.

1.2.3 数据处理

(1) 次降雨事件地表径流平均浓度EMC

降雨径流污染物浓度随降雨随机性也表现出相应的波动特征[17], 目前EMC是评价单次降雨径流对水质影响的重要指标, 计算公式为:

(1)

式中, EMC为次降雨径流事件流量加权平均浓度, mg·L-1; ctt时刻污染物浓度, mg·L-1; qtt时刻径流流量, L·s-1; T为单次降雨径流结束时间, min.

(2) 无因次累积负荷-累积径流量变化曲线M-V

单次降雨事件期间, 普遍使用累积污染物负荷-累积径流量变化曲线M-V能反映初期冲刷效应(污染物负荷增加率高于径流量增加率)[18].该方法在水平轴上绘制累积径流比Vt, 在垂直轴上绘制累积污染物负荷比Mt, 其中:

(2)
(3)

式中, t为单次降雨径流某一采样时间, min.

(3) MFFn

MFFn值是另一种定量分析污染物“初期冲刷效应”标准[19], 表示累积径流量占总径流量的n%时, 累积污染物负荷百分比与累计径流量百分比之比, 其计算公式为:

(4)

式中, n是从径流产生到t时刻累积径流百分比, 0~100%; M是次降雨事件污染物流失总量, g; V是次降雨事件径流总量, L.

本实验数据采用SigmaPlot 12.0和SPSS 19.0进行绘图和统计分析.

2 结果与分析 2.1 集水区降雨与径流特征

规模化柑橘种植区降雨主要分布在4~10月, 累积降雨量904.5 mm, 占全年降雨量的81%.由图 3可知, 4~5月降雨较为集中, 夏季存在季节性干旱.整个观测期集水区径流量为620.55 m3, 最大次降雨地表径流量为173.08 m3(2018年5月5日), 春季6场降雨径流量占全年的71.0%.对集水区降雨与径流进行相关分析发现, 次降雨径流量与雨量存在显著相关性(r=0.687, P < 0.01).

图 3 2018年集水区降雨与径流关系 Fig. 3 Rainfall and runoff volumes for the entire rainfall-runoff events during 2018

2.2 EMC与年流失负荷

表 1可知, 监测期间14场产流径流水中TN、DN和NN的EMC平均值为8.49、7.38和4.98 mg·L-1, 变异系数89.27%、95.61%和96.34%; AN的EMC平均值为0.45 mg·L-1, 变异系数275.83%.TP、DP、PP的EMC平均值为0.87、0.59和0.28 mg·L-1, 变异系数75.99%、72.71%和92.41%.较高变异系数体现了不同降雨事件径流氮、磷EMC存在较大差异, 特别是AN具有较高的变异系数, 施肥后首场降雨AN的EMC甚至达到最小次降雨径流EMC的240倍.

表 1 EMC与年降雨径流氮磷流失负荷(n=14) Table 1 EMC and annual N and P loss load of rainfall-runoff(n=14)

整个监测期TN、DN、NN、AN、TP、DP和PP的流失负荷分别为13.43、12.20、8.77、0.75、1.26、0.84和0.42 kg·hm-2, 略低于曹瑞霞等[20]对三峡库区新政小流域氮磷流失负荷估算值, 该研究得到的氮磷流失负荷是根据观测到的两场降雨径流平均氮磷浓度与这两场降雨占全年降雨量的比例估算而得, 监测的两场降雨均为强降雨事件, 使用其径流氮磷浓度作为全年氮磷的平均浓度, 可能会高估该流域氮磷流失负荷.本研究还发现, 虽然TN和DN的浓度虽随降雨事件变化, 但NN始终是其主要流失形态, 监测期内超过60%的降雨事件中NN/DN比值超过1/2, 且NN与DN之间呈显著的线性正相关(NN=0.751 6×DN-2.882, r=0.899 7, P<0.01), 而AN/DN比值较小(0.16%~18.69%), 且与降雨无显著相关关系(P>0.05).

2.3 氮磷流失过程

4~11月(除10月外)集水区均产生地表径流和径流养分流失, 整个观测期单场降雨事件TN、DN、NN、TP、DP和PP最大流失负荷为5.32、5.09、4.95、0.44、0.28和0.15 kg·hm-2(5月5日), 单场降雨事件AN最大流失负荷为0.69 kg·hm-2(4月21日, 施肥后第4 d, 施肥后的首场降雨).通过对各养分浓度和距施肥天数进行相关分析发现, AN流失浓度与距施肥天数呈显著负相关(P < 0.01), 施肥是影响顺坡沟垄柑橘园AN流失最显著的影响因子[21].

利用养分流失累积曲线描述集水区养分年内流失过程(图 4), 氮磷养分流失具有明显的阶段性特征, 大致可以分为5个时期.其中, 第I阶段(春季施肥之前)、第Ⅳ阶段(6月22日~9月20日)氮磷流失较为平稳; 第Ⅲ阶段(5月5日~6月22日)、第5阶段(9月20日之后)氮磷流失有小幅增加的趋势; 第Ⅱ阶段(春季施肥后~5月5日)也是春季施肥后的前两场产流降雨(两场降雨分别达到大雨、暴雨标准), 是全年径流养分快速流失期, 此阶段贡献了全年TN、DN、NN、AN、TP、DP和PP总流失量的70.4%、72.1%、68.9%、94.1%、67.1%、64.1%和73.0%, 顺坡沟垄柑橘园土壤养分流失是降雨与施肥耦合作用的结果.

图 4 集水区养分流失年累积曲线 Fig. 4 Cumulative curve of N and P loss in the orchard catchment

2.4 典型降雨径流养分流失特征 2.4.1 典型降雨径流特征

选择两场典型降雨事件(5月5日和5月17日)进行分析, 其降雨特征见表 2.两场降雨分别代表长历时暴雨和短历时大雨事件, 其中5月5日降雨为春季施肥后首场暴雨, 也是全年降雨量最大的次降雨事件.

表 2 典型降雨事件特征 Table 2 Characteristics of typical rainfall events

长历时暴雨事件(5月5日)降雨具有波动性, 前期径流随雨强增加迅速增大并出现第一次径流峰值, 而后伴随雨强降低径流有所减少, 当雨强再次增大后出现第二次径流峰值; 而短历时大雨事件(5月17日)降雨量与径流量趋势基本一致, 且具有明显的径流单峰峰值(图 5).两场降雨事件共产生径流213 m3, 约占全年径流量的1/3, 其中仅5月5日暴雨事件径流量占全年径流量的27.89%.受前期降雨量的影响[23], 5月5日降雨初损值为6 mm, 5月17日为11.4 mm.径流峰值较降雨峰值表现出滞后性, 其中5月17日大雨事件径流峰值出现在雨峰后6 min, 峰值流量达19.54 L·s-1.

图 5 两场典型降雨事件降雨-径流关系 Fig. 5 Rainfall-runoff relationship for the two typical rainfall events

此外还发现, 两场降雨事件均出现径流快速回落的现象, 这可能与顺坡沟垄坡面构型和集水区面积较小有关[24].在土层浅薄、基岩埋深较浅的坡面上, 土壤层自身涵水能力较低[25], 而对于沟垄坡面而言, 高强度降雨或前期连续降雨能够促进垄面土壤迅速饱和, 并伴随优先流的产生, 在基岩处形成侧向流, 回归进入垄沟地表径流[26], 沟垄构型加快径流的汇集, 同时固定的顺坡流路促进了径流排导[27], 因此该坡面构型下地表径流对降雨的响应快于紫色土区传统自然坡面和小流域[28].

2.4.2 典型降雨径流氮磷浓度

长历时暴雨径流样中TN、DN和NN浓度分别为(17.64±5.18)、(16.72±4.89)和(16.05±5.02) mg·L-1, AN浓度为(0.06±0.13) mg·L-1, TP、DP和PP浓度分别为(1.37±0.26)、(0.89±0.17)和(0.48±0.12)mg·L-1.通过配对样本t检验发现, 长历时暴雨径流TN、DN和NN之间无显著性差异(df=24, P>0.05), 且表现较好的一致性; 短历时大雨径流样中TN、DN、NN浓度分别为(18.79±6.13)、(16.47±5.42)和(10.62±2.93)mg·L-1, AN浓度为(0.02±0.01)mg·L-1, TP、DP和PP浓度分别为(0.91±0.23)、(0.65±0.18)和(0.26±0.06) mg·L-1, 见图 6.

图 6 典型降雨事件径流养分各形态氮磷浓度变化 Fig. 6 Changes of runoff N and P concentrations during two typical rainfall events

结合图 5的降雨径流关系发现, 长历时暴雨径流中各形态氮浓度峰值出现的时间早于径流峰值, 而各形态磷浓度峰值略晚于径流峰值.图 6可知, 暴雨径流前期(径流产生后400 min内)TN、DN和NN浓度峰值为22.57、20.68和20.41 mg·L-1, 浓度峰值出现后, 随着径流延长, TN、DN和NN浓度趋于相对稳定; AN在径流形成时浓度最高, 随着降雨推移而迅速降低且保持较低的平稳状态, 这与唐柄哲[29]的研究结果相似.本研究的典型暴雨事件距前期春季施肥时间较短, 且API值为29.29 mm, 土壤处于偏湿润状态, 王辉等[30]的研究发现土壤前期较高含水量能够快速溶解表层NN, 导致NN在径流前期快速流失, 而本研究NN与TN存在显著的线性关系(TN=0.956 7×NN-0.872 3, r=0.990 6), 各形态氮素流失浓度在短时间迅速增加达到峰值, 后期由于降雨量变化、淋溶和侧流补给, 浓度总体处于平稳波动, 但始终维持较高的浓度水平.径流磷各形态浓度随时间变化趋势较为接近, DP占径流TP的质量分数为65%, 溶解态磷是顺坡垄沟坡面集水区磷素随径流流失的主要形态.

短历时大雨径流前期(径流产生后150 min内)各形态氮浓度差异不显著(df=9, P>0.05), 当后期径流量小于1 L·s-1后, 主要是壤中流出流形成的地表径流, 以NN为主, 此阶段与长历时暴雨事件径流后期地表径流TN、DN平均浓度无显著性差异(df=15, P>0.05);各形态磷峰值出现的时间略晚于径流峰值, 这与于兴修等[25]在沂蒙山区小流域典型次降雨各形态磷素浓度峰值早于径流峰值的研究结果有所差异.

两场典型降雨事件径流氮磷流失负荷较高, 长历时暴雨事件地表径流TN和TP流失负荷为5.32 kg·hm-2和0.44 kg·hm-2, NN和DP流失负荷为4.95 kg·hm-2和0.28 kg·hm-2, 分别占TN和TP的95.69%和64.49%;短历时大雨事件地表径流TN和TP流失负荷为0.78 kg·hm-2和0.06 kg·hm-2, NN和DP流失负荷为0.52 kg·hm-2和0.05 kg·hm-2, 分别占TN和TP的88.08%和74.88%. NN和DP是顺坡沟垄柑橘集水区径流氮磷的主要流失形态, 特别在长历时暴雨事件中氮素几乎全部以NN形态流失, 氮素流失主要形态与张林等[31]在库首兰陵溪小流域观测结果基本一致.

2.4.3 初期冲刷效应

并不是任何降雨事件都能够出现“初期冲刷效应”[32], 降雨特征、下垫面特征和流域面积等均可能对污染物“初期冲刷效应”产生影响[33, 34].大雨和暴雨事件各形态氮磷的M-V曲线均位于1:1对角线的上方(图 7), 各形态氮流失主要集中在径流初期, 这与杨小林等[24]的研究结果有所不同, 其发现在紫色土丘陵区农林复合小流域径流中后期TN、DN和NN流失负荷较高, 在相对不透水比例较高的集镇小流域, TN和AN主要集中在径流初期阶段; 张林等[31]的研究发现在三峡库区农林复合小流域AN流失主要集中在暴雨径流前期, 而TN和NN在径流后期集中流失, 本研究对象柑橘园地表径流具有快速回落的特征, 加上表层土壤氮素随土壤水迅速通过垄沟沟壁侧渗出流补给地表径流, 能够实现氮素特别是NN的快速流失.各形态磷流失也主要集中在径流初期, 这与徐亚娟等[28]的研究发现在紫色土丘陵区小流域大雨及暴雨径流过程中, TP初期径流流失负荷较高, “初期冲刷效应”明显相似.

图 7 典型降雨事件氮磷污染物M-V曲线 Fig. 7 The M-V curves of different pollutants during two typical rainfall events

使用式(4)计算出不同形态氮磷污染物MFF20和MFF30(图 8), 可见各形态氮磷污染物的MFFn均值都大于1.0, 且MFF20值均大于MFF30值, 其中TN、DN、NN、AN、TP、DP和PP的MFF20比MFF30分别高22.5%、20.4%、20.2%、35.4%、26.5%、28.7%和29.5%, 说明在顺坡沟垄集水区前20%地表径流比前30%地表径流更容易发生“初期冲刷效应”. MFF20的计算结果可知, 前期20%的地表径流输出了58.0%的TN、57.0%的DN、58.5%的NN、79.0%的AN、62.0%的TP、63.5%的DP和60.0%的PP.

图 8 不同形态氮磷污染物MFF20与MFF30 Fig. 8 MFF20 and MFF30 for the different N and P pollutants indices

3 讨论

径流氮磷流失过程比较复杂, 降雨条件、农事活动、下垫面特征、地表覆盖和土壤理化性质等都会影响地表径流氮磷流失过程、特征[35, 36].三峡库区中部春季易出现集中性强降雨, 当与施肥同期时, 将可能造成土壤养分大量流失, 成为径流氮磷高输出期, 流失负荷甚至超过全年的70%, 蒋锐等[37]的研究发现紫色土丘陵区农耕地施肥后暴雨径流流失的TN负荷占年流失负荷的58%.施肥后首场暴雨径流氮磷平均浓度较高, 且氮磷流失以溶解态为主, 分别占总氮和总磷的80%和60%以上, 与顺坡种植能加快地表产汇流过程和径流冲刷能力[38], 促进颗粒态氮磷流失迁移的结论[39]有所不同, 本研究区规模化果园经营者因劳动力投入成本高, 在节省开支情况下一般仅对树盘周围进行人工除草作业, 其他部位杂草往往难以及时清除, 柑橘园常年保持较高的地表覆盖, 增加了地表粗糙度, 减少坡面侵蚀和降低地表径流颗粒态氮和磷的流失, 2018年的观测发现本集水区全年仅在5月5日暴雨事件中测流堰有悬移质泥沙输出(输沙模数仅为0.79 t·km-2), 也印证了地表覆盖在降低径流冲刷中的作用.

强降雨时地表径流水平运动强于垂直淋溶, 对地表具有一定的冲刷效应, 相对容易出现“初期冲刷效应”[28, 40].虽然强降雨对地表污染物不仅具有冲刷作用, 同样也具有稀释和溶解作用, 但在较小集水面上, 雨强越大, 污染物“初期冲刷效应”也越显著[41].本研究两场降雨虽然长历时暴雨降雨量大, 然而雨强(I60=11.2 mm·h-1)却远低于短历时大雨(I60=32.2 mm·h-1), 短历时大雨“初期冲刷效应”更为明显.不同形态氮磷污染物表现出不同“初期冲刷效应”[42], 在顺坡沟垄柑橘集水区, AN是“初期冲刷效应”最强的污染物, 其次为各形态磷和氮.Zeng等[43]的研究发现城市不透水下垫面AN、TP是径流初期阶段最主要的污染物, 但在川中丘陵紫色土小流域径流中后期氮素通过壤中流迁移集中汇流, TN、DN和NN呈现“末期冲刷效应”[24, 44], 这可能与降雨径流过程有关, 顺坡沟垄集水区强降雨径流过程曲线与城市不透水下垫面较为相似, 与具有“平缓退水”的紫色土小流域具有明显的差异.

柑橘园顺坡沟垄坡面构型以及粗放经营管理增加了面源氮磷流失浓度, 地表径流TN和TP年EMC均值高达8.49 mg·L-1和0.87 mg·L-1, 高于华玲玲等[45]对三峡库区农林复合为主的古夫河小流域监测到的氮磷年均排放浓度(总氮年平均浓度2.34 mg·L-1, 总磷年平均浓度0.056 mg·L-1), 也超过了《地表水环境质量标准》(GB 3838-2002)Ⅴ类标准的4.25倍和2.2倍, 且远高于水体富营养化的临界浓度(总氮为0.2 mg·L-1, 总磷为0.02 mg·L-1)[46].特别是春季施肥后的首场暴雨径流中溶解态氮和磷EMC甚至高达17.57 mg·L-1和0.97 mg·L-1, 这些溶解态氮和磷在适宜水流流态、光照和水温等条件下能够促进藻类生长繁殖[47, 48], 从而加剧了三峡库区入库支流水体的富营养化风险.

4 结论

(1) 规模化柑橘园顺坡沟垄的坡面构型加快了地表产汇流, 地表径流控制着集水区水文过程, 这与紫色土区传统自然坡面或农林复合小流域高比例透水下垫面以壤中流为主的汇流过程存在较大的差异.通过果园林下生草, 增加沟壁和沟道粗糙度, 建立截水沟等措施降低产汇流速度.

(2) 柑橘集水区年径流氮和磷流失负荷高达13.43 kg·(hm2·a)-1和1.26 kg·(hm2·a)-1, 远高于库首柑橘径流小区总氮和总磷流失负荷; 分别为1.68 kg·(hm2·a)-1和0.08 kg·(hm2·a)-1][49], 规模化顺坡沟垄柑橘园正成为三峡库区中部重要的面源污染输出地类, 将成为该区域种植业面源污染防控重点.

(3) 目前规模化柑橘园因管理粗放, 顺坡沟垄柑橘集水区地表径流TN和TP年EMC分别为8.49 mg·L-1和0.87 mg·L-1, 特别是春季施肥后的强降雨地表径流氮磷流失浓度达到极值, 是规模化柑橘园面源氮磷控制的关键时期.可以通过调整施肥时间、优化水肥管理(如采用水肥一体化)等措施降低规模化柑橘园面源氮磷浓度, 减轻对三峡库区水环境的影响.

(4) 顺坡宽沟窄垄柑橘园集水区具有明显的“初期冲刷效应”, 暴雨和短历时大雨期间前期20%的地表径流输出了超过50%的污染物负荷.因此, 在高强度降水期特别是施肥后的强降雨期, 可以利用规模化果园坡面配套的灌溉水塘或新建小型人工湿地等措施, 截留初期径流, 对于降低入库面源氮、磷污染物负荷具有重要的防控效果.

参考文献
[1] 徐江. 三峡库区柑桔产业发展的几点经验及启示[J]. 中国南方果树, 2010, 39(6): 78-80.
[2] 朱波, 况福虹, 高美荣, 等. 土层厚度对紫色土坡地生产力的影响[J]. 山地学报, 2009, 27(6): 735-739.
Zhu B, Kuang F H, Gao M R, et al. Effects of soil thickness on productivity of sloping cropland of purple soil[J]. Journal of Mountain Science, 2009, 27(6): 735-739.
[3] 樊登星, 余新晓, 贾国栋, 等. 北京山区灌草坡面水土流失特征及其影响因素[J]. 中国水土保持科学, 2014, 12(2): 24-28.
Fan D X, Yu X X, Jia G D, et al. Characteristics of soil and water loss and its influencing factors on slope scale in rocky mountain area of Beijing[J]. Science of Soil and Water Conservation, 2014, 12(2): 24-28.
[4] 王辉, 王全九, 邵明安. 人工降雨条件下黄土坡面养分随径流迁移试验[J]. 农业工程学报, 2006, 22(6): 39-44.
Wang H, Wang Q J, Shao M A. Laboratory experiments of soil nutrient transfer in the loess slope with surface runoff during simulated rainfall[J]. Transactions of the CSAE, 2006, 22(6): 39-44.
[5] 王全九, 杨婷, 刘艳丽, 等. 土壤养分随地表径流流失机理与控制措施研究进展[J]. 农业机械学报, 2016, 47(6): 67-82.
Wang Q J, Yang T, Liu Y L, et al. Review of soil nutrient transport in runoff and its controlling measures[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(6): 67-82.
[6] Dong W C, Cui S, Fu Q, et al. Modelling soil solute release and transport in run-off on a loessial slope with and without surface stones[J]. Hydrological Processes, 2018, 32(10): 1391-1400.
[7] 中华人民共和国环境保护部.长江三峡工程生态与环境监测公报2017[R].北京: 中华人民共和国环境保护部, 2017.
[8] 朱波, 汪涛, 王建超, 等. 三峡库区典型小流域非点源氮磷污染的来源与负荷[J]. 中国水土保持, 2010(10): 34-36.
Zhu B, Wang T, Wang J C, et al. Source and load of non-point source nitrous and phosphorus pollution of typical small watersheds in the Three Gorges Reservoir Area[J]. Soil and Water Conservation in China, 2010(10): 34-36.
[9] 刘方谊, 夏颖, 黄敏, 等. 湖北省三峡库区不同种植模式下农田地表径流氮磷流失特征[J]. 农业资源与环境学报, 2018, 35(6): 550-558.
Liu F Y, Xia Y, Huang M, et al. Characteristics of nitrogen and phosphorus losses from farmlands through surface runoff under different cropping patterns in the Three Gorges Reservoir area of Hubei Province[J]. Journal of Agricultural Resources and Environment, 2018, 35(6): 550-558.
[10] 陈玲, 刘德富, 宋林旭, 等. 不同雨强下黄棕壤坡耕地径流养分输出机制研究[J]. 环境科学, 2013, 34(6): 2151-2158.
Chen L, Liu D F, Song L X, et al. Characteristics of nutrient loss by runoff in sloping arable land of yellow-brown under different rainfall intensities[J]. Environmental Science, 2013, 34(6): 2151-2158.
[11] 何丙辉, 黄巍, 郭甜, 等. 不同施肥处理对紫色土坡耕地地表径流磷素流失的影响[J]. 重庆师范大学学报(自然科学版), 2012, 29(3): 42-47, 134-141.
He B H, Huang W, Guo T, et al. Characteristics of N and P loss in the soil of purple sloping farmland at different fertilization levels[J]. Journal of Chongqing Normal University (Natural Science Edition), 2012, 29(3): 42-47, 134-141.
[12] 江才伦, 彭良志, 曹立, 等. 三峡库区紫色土坡地柑橘园不同耕作方式的水土流失研究[J]. 水土保持学报, 2011, 25(4): 26-31, 35.
Wang C L, Peng L Z, Cao L, et al. Correlation between farming methods and soil erosion in the purple soil slope of citrus orchard in the Three Gorges Reservoir Area[J]. Journal of Soil and Water Conservation, 2011, 25(4): 26-31, 35.
[13] 廖晓勇, 罗承德, 陈治谏, 等. 三峡库区坡地果园间植草篱的水土保持效应[J]. 长江流域资源与环境, 2008, 17(1): 152-156.
Liao X Y, Luo C D, Chen Z J, et al. Functions of soil and water conservation by grass hedgerow intercropping of slope orchard in Three Gorges Reservoir Area[J]. Resources and Environment in the Yangtze Basin, 2008, 17(1): 152-156.
[14] 栾好安.三峡库区橘园绿肥的生态效应及其对柑橘产量和品质的影响[D].武汉: 华中农业大学, 2015.
Luan H A. Effects of planting green manure on ecological environment, yield and quality of citrus in citrus orchard at Three Gorges Reservoir[D]. Wuhan: Huazhong Agricultural University, 2015.
[15] 刘月娇, 倪九派, 张洋, 等. 三峡库区新建柑橘园间作的截流保肥效果分析[J]. 水土保持学报, 2015, 29(1): 226-230.
Liu Y J, Nie J P, Zhang Y, et al. Analysis of effects on water interception and nutrient conservation from citrus intercropping in hilly area of Chongqing city[J]. Journal of Soil and Water Conservation, 2015, 29(1): 226-230.
[16] Wickham J D, Jones K B, Wade T G, et al. Uncertainty in scaling nutrient export coefficients[A]. In: Wu J G, Jones K B, Li H, et al (Eds.). Scaling And Uncertainty Analysis in Ecology[M]. Dordrecht: Springer, 2006.
[17] 汪楚乔, 陈柔君, 吴磊, 等. 宜兴典型村落不同下垫面降雨径流污染物排放特征[J]. 生态与农村环境学报, 2016, 32(4): 632-638.
Wang C Q, Chen R J, Wu L, et al. Release of pollutants with rainfall-triggered runoff from different underlying surfaces in villages typical of Yixing[J]. Journal of Ecology and Rural Environment, 2016, 32(4): 632-638.
[18] Hathaway J M, Hunt W F. Evaluation of frst fush for indicator bacteria and total suspended solids in urban stormwater runoff[J]. Water Air and Soil Pollution, 2011, 217(1-4): 135-147.
[19] Ma Z B, Ni H G, Zeng H, et al. Function formula for first flush analysis in mixed watersheds:A comparison of power and polynomial methods[J]. Journal of Hydrology, 2011, 402(3-4): 333-339.
[20] 曹瑞霞, 刘京, 邓开开, 等. 三峡库区典型紫色土小流域径流及氮磷流失特征[J]. 环境科学, 2019, 40(12): 5330-5339.
Cao R X, Liu J, Deng K K, et al. Characteristics of nitrogen and phosphorus losses and runoff in a typical purple soil watershed in the Three Gorges Reservoir Area[J]. Environmental Science, 2019, 40(12): 5330-5339.
[21] Cookson W R, Rowarth J S, Cameron K C. The fate of autumn-, late winter- and spring-applied nitrogen fertilizer in a perennial ryegrass (Lolium perenne L.) seed crop on a silt loam soil in Canterbury, New Zealand[J]. Agriculture, Ecosystems & Environment, 2001, 84(1): 67-77.
[22] Perrone J, Enright P, Papineau F, et al. Overview of canadian watershed monitoring programs for improved water quality[J]. Canadian Water Resources Journal, 1998, 23(2): 121-134.
[23] Appels W M, Bogaart P W, van der Zee S E A T M. Surface runoff in flat terrain:how field topography and runoff generating processes control hydrological connectivity[J]. Journal of Hydrology, 2016, 534: 493-504.
[24] 杨小林, 朱波, 董玉龙, 等. 紫色土丘陵区小流域非点源氮迁移特征研究[J]. 水利学报, 2013, 44(3): 276-283.
Yang X L, Zhu B, Dong Y L, et al. Transport processes of diffuse nitrogen in typical catchments in the hilly area of purple soil[J]. Journal of Hydraulic Engineering, 2013, 44(3): 276-283.
[25] 于兴修, 李振炜, 刘前进, 等. 沂蒙山区典型小流域降雨径流的磷素输出特征[J]. 环境科学, 2012, 33(8): 2644-2651.
Yu X X, Li Z W, Liu Q J, et al. Output characteristics of rainfall runoff phosphorus pollution from a typical small watershed in Yimeng mountainous area[J]. Chinese Journal of Environmental Science, 2012, 33(8): 2644-2651.
[26] 边锋, 郑粉莉, 徐锡蒙, 等. 东北黑土区顺坡垄作和无垄作坡面侵蚀过程对比[J]. 水土保持通报, 2016, 36(1): 11-16.
Bian F, Zheng F L, Xu X M, et al. Comparison of soil erosion process between longitudinal ridge slope and non-ridge slope in mollisol region of northeast China[J]. Bulletin of Soil and Water Conservation, 2016, 36(1): 11-16.
[27] 王宇, 韩兴, 赵占军, 等. 垄沟秸秆覆盖对黑土顺坡耕地氮、磷养分阻控效果[J]. 水土保持学报, 2016, 30(1): 137-140.
Wang Y, Han X, Zhao Z J, et al. Controlling effects of furrow straw mulching on N and P nutrient loss on black soil slope farmland[J]. Journal of Soil and Water Conservation, 2016, 30(1): 137-140.
[28] 徐亚娟, 高扬, 朱宁华, 等. 紫色土流域次降雨条件下碳、磷非点源输出过程及其流失负荷[J]. 生态学报, 2014, 34(17): 5021-5029.
Xu Y J, Gao Y, Zhu N H, et al. The dynamic process of carbon and phosphorus losses by runoff and associated loss load in a farmland-dominated watershed of purple soil region[J]. Acta Ecologica Sinica, 2014, 34(17): 5021-5029.
[29] 唐柄哲.不同施肥类型及耕作方式下紫色土坡耕地径流氮、磷流失研究[D].重庆: 西南大学, 2016.
Tang B Z. Study on nitrogen and phosphorus loss process of purple soil under different fertility and tillage[D]. Chongqing: Southwest University, 2016.
[30] 王辉, 王全九, 邵明安. 前期土壤含水量对黄土坡面氮磷流失的影响及最优含水量的确定[J]. 环境科学学报, 2008, 28(8): 1571-1578.
Wang H, Wang Q J, Shao M A. Effect of antecedent soil moisture on Nitrate-N and Phosphorus loss from loess slope-land and determination of optimal moisture values[J]. Acta Scientiae Circumstantiae, 2008, 28(8): 1571-1578.
[31] 张林, 黄志霖, 肖文发, 等. 三峡库区兰陵溪小流域径流氮磷输出及其降雨径流过程特征[J]. 环境科学, 2018, 39(2): 792-799.
Zhang L, Huang Z L, Xiao W F, et al. Characteristics of nitrogen and phosphorus output in runoff and rainfall runoff in Lanlingxi Watershed, Three Gorges Reservoir Area[J]. Environmental Science, 2018, 39(2): 792-799.
[32] Hall M J, Ellis J B. Water quality problems of urban areas[J]. Geojournal, 1985, 11(3): 265-275.
[33] Taebi A, Droste R L. First flush pollution load of urban stormwater runoff[J]. Journal of Environmental Engineering and Science, 2004, 3(4): 301-309.
[34] Lang M, Li P, Yan X Y. Runoff concentration and load of nitrogen and phosphorus from a residential area in an intensive agricultural watershed[J]. Science of the Total Environment, 2013, 458-460: 238-245.
[35] 冯源嵩, 林陶, 杨庆媛. 百花湖周边城市近郊小流域氮、磷输出时空特征[J]. 环境科学, 2014, 35(12): 4537-4543.
Feng Y S, Lin T, Yang Q Y. Temporal and spatial characteristic of nitrogen and phosphorus output in the suburb watershed around the Baihua lake[J]. Environmental Science, 2014, 35(12): 4537-4543.
[36] 陈成龙, 高明, 倪九派, 等. 三峡库区小流域稻田空间格局对氮磷流失影响[J]. 环境科学, 2017, 38(5): 1889-1897.
Chen C L, Gao M, Ni J P, et al. Influence of spatial pattern of paddy field on the losses of nitrogen and phosphorus in Three Gorges Reservoir Area[J]. Environmental Science, 2017, 38(5): 1889-1897.
[37] 蒋锐, 朱波, 唐家良, 等. 紫色丘陵区小流域典型降雨径流氮磷流失特征[J]. 农业环境科学学报, 2008, 27(4): 1353-1358.
Jiang R, Zhu B, Tang J L, et al. Characteristics of nitrogen and phosphorus losses in typical rainfall-runoff events in a small watershed in hilly area of purple soil[J]. Journal of Agro-Environment Science, 2008, 27(4): 1353-1358.
[38] 任雨之, 郑江坤, 付滟, 等. 不同耕种模式下遂宁组紫色土坡耕地产流产沙特征[J]. 水土保持学报, 2019, 33(2): 30-38.
Ren Y Z, Zheng J K, Fu Y, et al. Characteristics of runoff and sediment yield in purple soil sloping farmland under different tillage patterns in Suining formation[J]. Journal of Soil and Water Conservation, 2019, 33(2): 30-38.
[39] 刘宗岸, 杨京平, 杨正超, 等. 苕溪流域茶园不同种植模式下地表径流氮磷流失特征[J]. 水土保持学报, 2012, 26(2): 29-32, 44.
Liu Z A, Yang J P, Yang Z C, et al. Characteristics of nitrogen and phosphorus losses with surface runoff in tea field of Tiaoxi watershed under different planting patterns[J]. Journal of Soil and Water Conservation, 2012, 26(2): 29-32, 44.
[40] 李振炜, 吴迪, 于兴修, 等. 沂蒙山区典型小流域次暴雨的氮磷输出特征[J]. 安徽农业科学, 2012, 40(12): 7269-7270, 7312.
Li Z W, Wu D, Yu X X, et al. Output characteristics of non-point nitrogen and phosphorus from a typica1 small watershed in Yimeng mountainous area during a rainstorm[J]. Journal of Anhui Agricultural Sciences, 2012, 40(12): 7269-7270, 7312.
[41] Lee J H, Bang K W. Characterization of urban stormwater runoff[J]. Water Research, 2000, 34(6): 1773-1780.
[42] Stutter M I, Langan S J, Cooper R J. Spatial contributions of diffuse inputs and within-channel processes to the form of stream water phosphorus over storm events[J]. Journal of Hydrology, 2008, 350(3-4): 203-214.
[43] Zeng J J, Huang G R, Luo H W, et al. First flush of non-point source pollution and hydrological effects of LID in a Guangzhou community[J]. Scientific Reports, 2019, 9: 13865.
[44] 汪涛, 朱波, 罗专溪, 等. 紫色土坡耕地硝酸盐流失过程与特征研究[J]. 土壤学报, 2010, 47(5): 962-970.
Wang T, Zhu B, Luo Z X, et al. Nitrate loss from sloping cropland of purple soil[J]. Acta Pedologica Sinica, 2010, 47(5): 962-970.
[45] 华玲玲, 李文超, 翟丽梅, 等. 三峡库区古夫河小流域氮磷排放特征[J]. 环境科学, 2017, 38(1): 138-146.
Hua L L, Li W C, Zhai L M, et al. Characteristics of nitrogen and phosphorus emissions in the Gufu river small watershed of the three Georges reservoir area[J]. Environmental Science, 2017, 38(1): 138-146.
[46] 罗固源, 卜发平, 许晓毅, 等. 三峡库区临江河回水区总氮和总磷的动态特征[J]. 土木建筑与环境工程, 2009, 31(5): 106-111.
Luo G Y, Bu F P, Xu X Y, et al. Dynamics of TN and TP in the backwater region of the Linjiang River in the Three Gorge Reservoir[J]. Journal of Chongqing Jianzhu University, 2009, 31(5): 106-111.
[47] 雷沛, 张洪, 单保庆. 丹江口水库典型入库支流氮磷动态特征研究[J]. 环境科学, 2012, 33(9): 3038-3045.
Lei P, Zhang H, Shan B Q. Dynamic characteristics of nitrogen and phosphorus in the representative input tributaries of Danjiangkou Reservoir[J]. Environmental Science, 2012, 33(9): 3038-3045.
[48] 高永霞, 朱广伟, 贺冉冉, 等. 天目湖水质演变及富营养化状况研究[J]. 环境科学, 2009, 30(3): 673-679.
Gao Y X, Zhu G W, He R R, et al. Variation of water quality and trophic state of lake Tianmu, China[J]. Environmental Science, 2009, 30(3): 673-679.
[49] 曾立雄.三峡库区兰陵溪小流域养分的分布、迁移与控制研究[D].北京: 中国林业科学研究院, 2010.
Zeng L X. Distribution, migration and control of nutrient in Lanlingxi watershed, Three Gorges Reservoir Area[D]. Beijing: Chinese Academy of Forestry, 2010.