洱海流域不同施肥模式对稻田氮磷径流流失的影响

姜海斌，沈仕洲，谷艳茹，于双，邹洪涛，Arif Husain，张磊，张克强*

（1.农业农村部环境保护科研监测所，天津300191；2.沈阳农业大学土地与环境学院，沈阳110866；3.国家农业环境大理观测实验站，云南 大理671004；4.大理市农业环境保护监测站，云南 大理671000）

洱海是云南省第二大高原湖泊，水域面积252 km2，是大理市主要饮用水源地，也是苍山洱海国家级自然保护区和风景名胜区的核心，被称为大理人民的“母亲湖”[1]。自20世纪90年代以来，洱海水质呈下降趋势（从Ⅱ类水质下降至Ⅲ类水质），处于中度富营养化阶段[2]。研究表明，目前农业面源污染已成为洱海流域富营养化的主要污染源之一[3]，而水稻又是流域主要的种植作物，占农作物播种总面积34.80%，因水稻种植期降雨量较大[4]，径流是农田氮磷流失的重要途径，因此研究水稻种植过程中氮磷径流流失是削减农业面源污染的重要方式之一。多项研究表明，施肥是影响稻田氮磷流失的主要原因之一。姚金铃等[5]对洱海流域稻田氮磷径流流失的研究表明，常规施肥条件下，稻田氮、磷径流流失量分别为6.80 kg·hm-2和1.30 kg·hm-2，占施肥总量的8.92%和16.27%。Lee等[6]研究表明，化肥减量施用可有效减少氮磷流失量，提高化肥利用率。蔡佳佩等[7]研究表明，化肥减量20%配施1 500 kg·hm-2有机肥能够有效降低施肥后7 d内田面水中总氮和总磷含量，同时降低氮磷径流流失风险，有效增加水稻籽粒和秸秆产量。石敦杰等[8]的研究结果表明，控释氮肥发挥了肥效长的优势，能够满足水稻中后期的生长需求，使水稻产量提高了8.30%，而田面水总氮含量比常规施肥降低了89.11%。前人在洱海流域开展的氮磷流失研究，多是测定稻田田面水氮磷浓度，分析氮磷流失风险，并未设置径流池直接测定氮磷径流流失量，且对不同施肥模式下稻田的氮磷径流流失规律研究较少。为此，本研究设置化肥减施、有机肥替代和缓控释肥等施肥处理，通过田间径流小区试验，研究洱海流域不同施肥模式对稻田氮磷径流流失的影响，并结合水稻产量，探究有效减少稻田氮磷径流流失的施肥模式，为洱海流域水稻种植氮磷减排提供科学数据支撑。

1 材料与方法

1.1 试验地概况

试验地位于云南省大理市喜洲镇，农业农村部环境保护科研监测所大理综合实验站内（25°53′34″N，100°10′27″E）。该地区属于典型低纬高原中亚热带西南季风气候类型，海拔1 980 m，气候温和，光照充足，年平均气温14.6℃，多年平均降雨量为1 048 mm[9]。供试土壤为水稻土，0～20 cm土壤基本理化性质为：有机质57.3 g·kg-1，pH 7.1，全氮3.3 g·kg-1，全磷0.9 g·kg-1，硝态氮21.64 mg·kg-1，铵态氮14.22 mg·kg-1，有效磷35.3 mg·kg-1。

1.2 试验设计

试验共设8个处理：不施肥（CK）；常规施肥（CF）；常规施肥减量20%（T1）；有机肥等氮替代T1（T2）；有机肥等磷替代T1，尿素补齐氮肥（T3）；考虑有机肥矿化率25%，以氮计，有机肥替代T1（T4）；考虑有机肥矿化率25%，以磷计，有机肥替代T1，尿素补齐氮肥（T5）；缓控释肥（T6）。各处理施肥量见表1，每个处理重复3次，共24个小区，随机区组设计。每个试验小区面积为30 m2（6 m×5 m），小区间用水泥砂浆筑埂，埂宽0.24 m、埂高0.20 m，地表面以下筑入1 m，径流口距田面0.15 m。每个径流小区在排水端设径流收集池，径流池容积为4.5 m3（宽×长×深=1 m×4.5 m×1 m），池表面覆盖彩钢瓦密封，防止雨水进入。

1.3 田间管理

供试水稻品种为云粳25，种植时间为2019年5月30日—10月11日。水稻种植施肥管理按当地种植管理方式，其中常规施肥氮肥为尿素，磷肥为过磷酸钙，钾肥为硫酸钾；缓控释肥氮磷钾比例为23∶10∶12；有机肥为牛粪商品有机肥，含有机质14.5%，氮（N）2.3%，磷（P2O5）2.4%，钾（K2O）5.7%。有机肥作基肥，在翻耕前一次性施入，翻耕深度约20 cm，缓控释肥和磷钾肥在分蘖期一次性施入，单施化肥处理（CF、T1）和有机无机配施处理（T3、T5）的尿素分两次施入，在分蘖期施入70%，在穗肥期施入30%。所有施肥处理肥料均人工撒施。2019年5月30日施入有机肥并且翻耕土地，6月3日插秧，株距9 cm，行距23 cm，每蔸2～3株，6月12日施分蘖肥，8月2日施穗肥，10月11日水稻收获。稻田田面水高度维持在10 cm左右，水面低至2～3 cm时补充灌溉水，每周约2次，在水稻生长中期人工除草1次。

1.4 样品采集和测定

稻田产生径流时采集径流水样，并测量径流池水面高度。采集前先混匀水样，再将水样装于100 mL塑料瓶带回实验室测定，不能及时测定时将水样放于-4℃冰箱冷冻保存。采集完毕后排空径流池，关闭上部彩钢瓦以防止降雨直接进入径流池。分别测定水样中总氮、铵态氮、硝态氮、总磷、可溶磷浓度：总氮采用碱性过硫酸钾消解紫外分光光度法，铵态氮采用纳氏试剂分光光度法，硝态氮采用紫外分光光度法，总磷采用钼酸铵分光光度法，可溶性磷采用钼锑抗比色法[10]。水稻成熟后分别收获每个小区的水稻籽粒和秸秆测产。

表1不同处理施肥类型和施肥量（kg·hm-2）Table 1 Types and amount of fertilization application of different treaments（kg·hm-2）

1.5 数据处理与分析

使用Microsoft Excel 2010进行数据整理，Origin 2018绘图，采用SPSS19单因素方差分析进行显著性检验（P＜0.05表示差异显著），LSD和Duncan法进行均值比较。

2 结果与分析

2.1 水稻季降雨和径流特征

水稻季降雨量和径流量如图1所示。水稻季降雨量变化范围为0.20～43.20 mm·d-1，降雨量总和为407.20 mm，集中发生在7月，其占整个水稻季降雨量的49.40%。生长期内共采集径流水样6次，径流变化范围为25.40～195.00 m3·hm-2，径流总量为476.50 m3·hm-2，其中4次径流集中在7月发生，占总径流量的62.40%。径流量除直接受降雨量影响以外，还受灌溉和水稻生长情况以及降雨强度等因素影响，这些导致径流量与降雨量不呈线性关系。

2.2 不同施肥模式稻田径流氮素浓度

不同施肥模式稻田径流总氮浓度如图2所示。整个水稻季径流总氮浓度范围为0.27～10.39 mg·L-1，平均浓度为3.69 mg·L-1。施分蘖肥后第6 d总氮浓度最高，为1.61～10.39 mg·L-1。7月2日除CF处理外，其他各处理总氮浓度均低于3.72 mg·L-1。7月5日降雨量较大，各施肥处理径流总氮浓度都有升高且较为接近，为6.85～7.98 mg·L-1。7月27日所有施肥处理的总氮浓度降至最低，平均为0.27 mg·L-1。8月24日由于穗肥的施入，各施肥处理浓度略有升高。在水稻不同阶段施肥后，各施肥处理径流总氮浓度均表现为分蘖肥期＞穗肥期。整个水稻生长期，化肥处理（CF和T1）径流总氮浓度高于施用有机肥替代化肥（T2～T5）和缓控释肥（T6）处理，这与宁建凤等[11]的研究结果一致。

不同施肥模式稻田径流铵态氮浓度如图3所示。整个水稻季径流铵态氮浓度范围为0.17～6.19 mg·L-1，平均浓度为0.80 mg·L-1。6月18日CF处理铵态氮浓度高至6.19 mg·L-1。从7月2日开始径流铵态氮浓度处于较低水平，且趋于稳定，在0.13～1.35 mg·L-1之间。分蘖肥的施入提高了稻田径流铵态氮浓度，6月18日采集的稻田径流铵态氮浓度均高于其余5次采集的浓度，后期径流铵态氮浓度低于1.35 mg·L-1，由于穗肥的施入量较小并且施入时间距离发生径流时间间隔较长，因此稻田径流铵态氮浓度没有显著升高，各施肥处理铵态氮浓度峰值均表现为分蘖肥＞穗肥。在整个水稻生长阶段，CK处理铵态氮浓度均低于0.35 mg·L-1。

不同施肥模式径流硝态氮浓度如图4所示。水稻径流硝态氮浓度受施肥和降雨影响较为显著，整个水稻季径流硝态氮浓度范围为0.21～3.00 mg·L-1，平均浓度为1.33 mg·L-1。施入分蘖肥后，第一次采集的稻田径流各施肥处理硝态氮浓度都处于相对较高水平，高于1.35 mg·L-1。7月2日各施肥处理硝态氮浓度显著下降，7月5日受强降雨影响，各施肥处理硝态氮浓度都略有回升，提高了径流池水样硝态氮浓度，这和陈永高等[12]的研究结果一致。补施穗肥氮素后，8月24日各施肥处理硝态氮浓度呈现升高趋势。T4处理初始硝态氮浓度仅为1.76 mg·L-1，在8月24日其硝态氮浓度达到3.00 mg·L-1。

2.3 不同施肥模式径流磷素浓度

不同施肥模式径流总磷浓度如图5所示。稻田径流总磷浓度在0.03～0.91 mg·L-1之间，平均浓度为0.30 mg·L-1，在整个水稻生长期内，径流总磷浓度趋于稳定。T6处理总磷浓度高于其他施肥处理，平均浓度为0.71 mg·L-1，其中7月24日和7月27日采集的径流样总磷浓度为0.52 mg·L-1和0.91 mg·L-1。T4处理前3次采样总磷浓度较低，在0.23～0.28 mg·L-1之间，从7月24日开始升高，7月27日达到峰值，为0.61 mg·L-1。

不同施肥模式径流可溶磷浓度如图6所示。稻田径流可溶磷的浓度在0.04～0.45 mg·L-1之间，平均浓度为0.34 mg·L-1，各施肥处理可溶磷浓度普遍低于0.50 mg·L-1，较为稳定。6月18日CF处理可溶磷浓度高于其他施肥处理，7月24日和27日采集的水样可溶磷浓度显著低于T4和T6处理。T4的可溶磷浓度较高，在0.17～0.35 mg·L-1之间，前5次采样中可溶磷浓度缓慢升高，在7月27日达到峰值0.35 mg·L-1，随后可溶磷浓度开始降低。T6的可溶磷浓度同样有升高的趋势，在7月24日达到峰值0.45 mg·L-1。

2.4 不同施肥模式对稻田径流氮磷流失量的影响

稻田径流氮素流失量如表2所示。稻田径流总氮流失量呈T1＞CF＞T6＞T3＞T4＞T5＞CK＞T2，CF和T1处理的总氮流失量分别为32.26 kg·hm-2和32.30 kg·hm-2，都处于较高水平，两者之间差异不显著。T6处理总氮流失量为27.18 kg·hm-2，比CF的总氮流失量降低了15.7%。施用有机肥的T3、T4、T5和T2处理总氮流失量分别比CF低10.19、12.47、12.62 kg·hm-2和13.03 kg·hm-2，即分别降低31.6%、38.6%、39.1%和40.4%。稻田径流铵态氮流失量在2.56～8.55 kg·hm-2之间，占稻田径流总氮的13.2%～30.6%，呈CF＞T1＞T4＞T2＞T3＞T5＞CK＞T6，CF处理的铵态氮流失量最大，为8.55 kg·hm-2。T1处理可降低稻田径流铵态氮流失量，较CF处理降低了2.03 kg·hm-2。施用有机肥的T4、T2、T3、T5处理铵态氮流失量分别较CF处理降低2.49、2.98、4.00、5.17 kg·hm-2，即分别降低29.1%、34.9%、46.8%和60.5%。T6处理对稻田径流铵态氮流失量的削减效果最好，铵态氮浓度低至2.16 kg·hm-2。稻田径流硝态氮占总氮流失量的22.5%～37.2%，稻田径流硝态氮的流失量呈CF＞T1＞T4＞T5＞T3＞T2＞T6＞CK的趋势，CF处理的硝态氮流失量最大，达到9.12 kg·hm-2。T1处理有效降低了硝态氮浓度，比CF低了14.7%。T6处理也降低了硝态氮的流失量，降幅为32.9%。施用有机肥处理的T4、T5、T3和T2处理都不同程度降低了硝态氮浓度，降幅为19.3%～30.2%。

表2稻田径流氮磷流失量（kg·hm-2）Table 2 Nand Ploss in runoff water fromrice fields（kg·hm-2）

稻田径流磷素流失量如表2所示。稻田径流总磷流失量和可溶磷流失量均呈T4＞T6＞T1＞CF＞T2＞CK＞T5＞T3。与CF处理相比，T1、T4和T6处理增加了总磷流失量的16.3%、66.3%和45.3%，而T3和T5降低了总磷流失量的23.3%和17.4%，T2无显著差异。稻田径流可溶磷流失量占总磷流失量的83.6%～98.8%，磷素流失形态以可溶磷为主。

2.5 不同施肥模式对水稻产量的影响

不同施肥模式水稻籽粒和秸秆产量见表3。稻田籽粒产量呈T6＞CF＞T1＞T5＞T3＞T4＞T2＞CK。T6处理的籽粒产量达到11.24 t·hm-2，与CF处理无显著差异，T1处理的籽粒产量为10.00 t·hm-2，比CF处理下降9.5%，施用有机肥处理（T2、T3、T4和T5）的籽粒产量比CF处理下降11.8%～42.9%，其中T2籽粒产量仅为6.31 t·hm-2。但和等氮（磷）投入的化肥减量模式T1处理相比，施用有机肥模式的T3和T5处理的籽粒产量无显著差异。CK水稻籽粒产量最低，仅为5.91 t·hm-2。稻田秸秆产量与稻田籽粒产量较为相似，呈T6＞CF＞T1＞T4＞T3＞T5＞CK＞T2。T6、T1和T4处理秸秆产量与CF处理无显著差异，施用有机肥的T2、T3和T5处理均不同程度降低了稻田秸秆的产量，降幅为15.5%～36%。CK水稻秸秆产量最低，仅为5.50 t·hm-2。

3 讨论

3.1 不同施肥模式对氮素径流流失的影响

本研究中，稻田径流总氮浓度受施肥和降雨影响较大，整个水稻季径流总氮浓度范围为0.27～10.39 mg·L-1。6月18日总氮浓度在1.61～10.39 mg·L-1，其原因为6月12日施入分蘖肥所致，化肥施入稻田后除一部分被作物吸收外，大部分会存在于田面水中，从而使此次采集的径流水样浓度最高。7月5日降雨量较大，各施肥处理下径流总氮浓度都有升高且较为接近，处于6.85～7.98 mg·L-1之间。强降雨对稻田土壤的冲击力较大，导致农田土壤随径流水流失的风险加大，原本存在于土壤中的氮素也会溶于水，提高了径流总氮浓度。6月18日采集的稻田径流铵态氮浓度在0.49～6.19 mg·L-1，高于其余5次稻田径流铵态氮的浓度，其余径流铵态氮浓度均低于1.35 mg·L-1，说明铵态氮在水稻生长前期明显高于水稻生长后期，这主要由于施肥主要在水稻生长前期，且水稻生长量偏小，对肥料养分需求小。水稻径流硝态氮浓度受施肥和降雨影响较为显著，整个水稻季径流硝态氮浓度变化范围为0.21～3.00 mg·L-1。6月18日各施肥处理硝态氮浓度均值为1.92 mg·L-1，8月24日各施肥处理硝态氮浓度均值为2.48 mg·L-1，高于其他4次稻田径流硝态氮浓度，这可能是因为前期水稻秧苗较小对氮素利用较少，且灌溉水、雨水和施入穗肥补充了氮素所致。

表3水稻籽粒和秸秆产量（t·hm-2）Table 3 Grain and straw yield of rice（t·hm-2）

不同施肥处理下稻田径流总氮流失量为19.23～32.30 kg·hm-2之间，稻田径流总氮流失量呈T1＞CF＞T6＞T3＞T4＞T5＞CK＞T2。CF处理施用的尿素为无机养分，在土壤中释放较快，而水稻在生长期内需要持续供应养分，所以导致部分未被利用的养分随径流流失[13]。如表2所示，CF处理的总氮流失量为32.26 kg·hm-2，T1和T6处理的总氮流失量和CF处理无显著差异，但显著高于有机肥处理T2～T5和CK处理（P＜0.05）。与肖雪玉等[14]的研究结果一致，T6处理的总氮流失量、铵态氮流失量和硝态氮流失量分别为27.18、2.16 kg·hm-2和6.12 kg·hm-2，均低于CF处理的氮素流失量。与吴美玲等[15]研究结果相似，施用有机肥可以显著降低水稻径流总氮浓度，本研究T2、T3、T4和T5处理总氮流失量为19.23、22.07、19.79和19.64 kg·hm-2，分别较CF处理降低40.4%、31.6%、38.6%、39.1%，施用有机肥同样降低了铵态氮和硝态氮流失量。张玉平等[16]研究结果也表明，在等氮、磷养分施用条件下，有机无机肥配施有利于降低地表氮磷径流损失。

3.2 不同施肥模式对磷素径流流失的影响

本研究中，稻田径流总磷和可溶磷浓度较为稳定，稻田径流总磷浓度在0.03～0.91 mg·L-1之间，平均浓度为0.30 mg·L-1。T6处理总磷浓度在水稻生长前期和水稻生长后期均高于其他施肥处理，前期可能与缓控释肥溶解缓慢，部分未溶解的缓控释肥直接通过径流损失，导致径流总磷浓度升高。后期缓控释肥养分缓慢释放提高了田面水的浓度，进而提高了径流总磷浓度。T4处理总磷浓度在7月24日开始升高，这与有机肥使用量大和有机肥养分释放缓慢有关。稻田径流可溶磷的浓度在0.04～0.45 mg·L-1之间，平均浓度为0.34 mg·L-1。水稻生长中期各施肥处理径流可溶磷浓度高于水稻生长前期，这可能是稻田长期在淹水条件下受降雨扰动，释放出了一部分被固持在土壤中的磷素，从而使径流可溶磷浓度升高。后期各施肥处理径流水样可溶磷浓度开始降低，这与作物中后期需磷量升高有关。不同处理间可溶磷浓度变化规律与总磷一致，T4和T6处理在7月24日和7月27日磷素浓度明显高于其他施肥模式磷素浓度。

不同施肥模式下稻田径流总磷流失量为0.66～1.43 kg·hm-2，稻田径流总磷流失量呈T4＞T6＞T1＞T2＞CF＞CK＞T5＞T3。与CF处理相比，T1、T4和T6处理增加了总磷流失量的16.3%、66.3%和45.3%，而T3和T5处理降低了总磷流失量的23.3%和17.4%，T2处理无显著差异。T4处理有机肥投入量高达26 000 kg·hm-2，导致该处理径流总磷流失量为1.43 kg·hm-2[17-18]。T6处理水稻径流总磷流失量较大，为1.25 kg·hm-2，这可能是由于高强度降雨对土壤的冲刷强度较大，雨滴溅蚀破坏土壤团聚体结构，分散表层土壤，土壤颗粒被径流水携带进入径流池，径流中泥沙颗粒含量较高，泥沙中的缓控释肥导致产流中磷流失量升高[19]。有机肥处理中的T3和T5处理的总磷流失量为0.66 kg·hm-2和0.71 kg·hm-2。T3处理投入的磷素全部来自有机肥，有机肥中的磷素以有机态磷为主，不能直接被作物吸收利用，有机态磷需要一定时间水解矿化成无机态磷，同时T3处理有机肥投入量为2 609 kg·hm-2，低于其他3个有机肥处理施肥量，所以T3总磷流失量处于最低水平[20]。T5处理的有机肥投入量高于T2处理，有机肥提供所需磷素，用尿素补充所缺氮素，导致T5处理的有机肥投入量的总磷流失量低于T2处理，化肥配施有机肥可明显降低磷素流失量，有利于水稻高产稳产[21-22]。T4处理的总磷流失量有机肥施用量为26 000 kg·hm-2，总磷流失量为1.43 kg·hm-2，比CF处理增加66.3%，说明施用有机肥量增至一定程度时，反而会增加总磷流失量[23]。水稻径流可溶磷流失量趋势与总磷流失量类似，总体呈T4＞T6＞T1＞CF＞T2＞T5＞T3＞CK，可溶磷流失量占总磷流失量83.6%～98.8%，说明磷素的流失以可溶态为主。

3.3 不同施肥模式对水稻产量的影响

水稻籽粒产量呈T6＞CF＞T1＞T5＞T3＞T4＞T2＞CK，T6处理籽粒产量达到11.24 t·hm-2，与CF处理的籽粒产量无显著差异。缓控释肥作为新型肥料，能很好地满足作物不同生长时期的养分需求，通过提高作物生长中后期的供肥能力，促进速效养分对作物生长的供应，从而提高了水稻籽粒和秸秆的产量[24]。T1处理籽粒产量为10.00 t·hm-2，和CF处理相比未出现显著减产。与T1处理相比，施用等磷投入的有机肥处理T3和T5在施用有机肥的基础上增施尿素提供氮源，在有机无机配施的基础上能保证作物的产量，3个处理间籽粒产量无显著差异，这和蔡媛媛等[25]的研究结果一致。有机肥施用处理T4和T2其籽粒产量均有所下降，其中T2降低最多，仅为6.31 t·hm-2，其原因是仅投入有机肥料，有机肥释放养分较为缓慢，不能满足水稻生育前期营养需求，进而造成水稻减产[26]。CK处理水稻籽粒产量最低，仅为5.91 t·hm-2，说明施肥可有效提高水稻籽粒产量。水稻秸秆产量趋势与水稻籽粒产量较为相似，呈T6＞CF＞T1＞T4＞T3＞T5＞CK＞T2。可以发现，T6处理同时提高了水稻籽粒和秸秆的产量，有机肥施用的T5、T3和T2处理秸秆产量同样下降，而T4处理秸秆产量大于其他有机肥处理，可能原因是T4处理有机肥投入较高，造成水稻贪青晚熟，秸秆产量较高，从而使水稻籽粒产量降低[27]。

4 结论

（1）与常规施肥相比，化肥减量处理的总氮径流流失量无显著差异，有机肥替代处理和缓控释肥处理均降低总氮径流流失量；与常规施肥相比，有机肥替代处理T2总磷径流流失量无显著差异，高量施用有机肥和施用缓控释肥处理均增加总磷径流流失量。

（2）与常规施肥相比，施用缓控释肥处理的籽粒产量无显著差异，化肥减量和有机肥替代处理的籽粒产量分别降低9.5%和11.8%～42.9%，与减量施肥相比，施用有机肥处理T3和T5的籽粒产量无显著差异；与常规施肥相比，化肥减量、施用有机肥处理T5和施用缓控释肥处理的秸秆产量无显著差异，其他施肥处理的秸秆产量降低15.5%～36.0%。

（3）综合考虑不同施肥模式的稻田径流流失量和水稻产量等因素，化肥减量、有机无机配施模式和施用缓控释肥是可供选择的环境友好型施肥模式。