张富林 吴茂前 夏 颖 翟丽梅 段小丽 范先鹏† 熊桂云刘冬碧 高 立
(1 湖北省农业科学院植保土肥研究所,湖北省农业面源污染防治工程技术研究中心,农业农村部潜江农业环境与耕地保育科学观测实验站,农业农村部废弃物肥料化利用重点实验室,农业环境治理湖北省工程研究中心,武汉 430064)
(2 中国农业科学院农业资源与农业区划研究所,北京 100081)
(3 浠水县农业环境保护站,湖北浠水 438200)
改革开放以来,我国农业取得了举世瞩目的成就,以不足世界 9%的耕地养活了22%的人口,粮食产量实现“十二连增”,但同时,我国肥料利用率较低,肥料养分通过径流、渗漏和挥发等途径损失导致的农业面源污染问题日益突出[1-2]。根据2007年第一次全国污染源普查结果,农业面源化学需氧量、总氮、总磷年排放量已达1 320 万t、270.5 万t和28.5 万t,分别占全国排放总量的43.7%、57.2%和67.4%。农业面源污染已成为水体污染、湖泊富营养化的主要原因[3]。
江汉平原是我国水稻主要产区之一,该区域水稻生产对保障我国粮食安全发挥着重要作用。但因肥、水管理粗放,降雨频繁,该区域水稻种植引发的氮磷面源污染问题不容小视。特别是近年来,随着种植方式的改变,直播稻和机插秧稻已成为该区域主要的水稻种植方式。在这些种植方式下,人为排放刚施肥泡田水的现象非常普遍,这极大地加重了该区域的氮磷面源污染负荷。又由于该区域水网密集,而且农田沟渠水网与河湖水网紧密相连,稻田面源污染会直接威胁区域农业用水安全和饮用水安全。因此,防控稻田氮磷面源污染已是该区域目前亟待开展的重要任务。
田面水中氮磷是稻田面源污染的直接来源,掌握田面水氮磷动态特征是防控稻田氮磷面源污染的重要前提条件之一。目前,已有一些关于水稻田面水氮磷动态变化的研究,但多集中于太湖、洱海等地区[4-9],鲜有关于江汉平原水稻田的研究报道。而且已有研究显示,不同地区稻田田面水中氮磷动态变化特征不同[10]。在洱海北部地区的研究表明,施肥后 9 d 内是控制稻田氮素损失的关键时期,控制磷素损失的关键时期则是在施肥后的两周内[9]。对陕南汉江、丹江流域的研究表明[10],控制氮损失的关键时期是施肥后5 d 内,施磷后9 d内是控制田面水磷流失的关键时期。因此,有必要对江汉平原地区水稻田面水氮磷动态特征进行研究。本研究通过大田试验,研究了江汉平原稻田田面水中氮磷形态及浓度动态变化特征,以期为指导当地施肥、防控稻田氮磷面源污染和氨挥发提供科学依据。
本试验于江汉平原腹地的农业农村部潜江农业环境与耕地保育科学观测实验站内(湖北省潜江市浩口镇柳洲村)进行。试验点所在区域属亚热带季风气候区,年均降水量为1 100 mm,年均气温16.1 ℃,全年无霜期约250 d,土壤为潴育型水稻土。试验初期土壤耕层基本理化性质:容重1.07 g∙kg-3,有机质22.8 g∙kg-1,全氮1.87 g∙kg-1,全磷0.78 g∙kg-1,有效磷14.68 mg∙kg-1,硝态氮14.14 mg∙kg-1,铵态氮4.59 mg∙kg-1,pH 6.90。
氮用量试验:设不施氮(N0)、施纯氮157.5 kg∙hm-2(N1)、施纯氮210 kg∙hm-2(N2)、施纯氮262.5 kg∙hm-2(N3)、施纯氮315 kg∙hm-2(N4)、施纯氮420 kg∙hm-2(N5)6个处理,每个处理三次重复,共计18个小区,随机区组排列。各处理均配施等量的磷、钾肥,施用量为P2O575 kg∙hm-2,K2O 90 kg∙hm-2。
磷用量试验:设不施磷(P0)、施P2O537.5 kg∙hm-2(P1)、施P2O575 kg∙hm-2(P2)、施P2O5112.5 kg∙hm-2(P3)、施P2O5150 kg∙hm-2(P4)、施P2O5300 kg∙hm-2(P5)6个处理,每个处理三次重复,共计18个小区,随机区组排列。各处理均配施等量的氮、钾肥,施用量为N 210 kg∙hm-2,K2O 90 kg∙hm-2。
小区面积均为40 m2。试验所用氮肥为尿素(含N 463 g∙kg-1),分基肥(40%)、分蘖肥(30%)和穗肥(30%)施用;磷肥为过磷酸钙(含P2O5120 g∙kg-1),全部作为基肥施用;钾肥为氯化钾(含K2O 600 g∙kg-1),分基肥(50%)和穗肥(50%)施用。底肥全小区撒施,施后再耖田一次;分蘖期和穗肥全小区撒施。水稻品种为广两优476,按16.7 cm × 26.7 cm的密度人工栽秧,每蔸2株。试验于4月28日播种,5月30日移栽,9月17日收获。试验过程中水分管理按当地习惯进行,插秧时不排放泡田水,群体80%够苗后自然断水晒田,收获前一周自然断水搁田。
每次施氮肥后,每天上午8:00~10:00采集田面水样,每个小区采集8个样点,混合后测定总氮(TN)、可溶性总氮(DTN)、颗粒态氮(PN)、可溶性有机氮(DON)、和施用磷肥后,每天上午8:00~10:00采集田面水样,每个小区采集8个样点,混合后测定总磷(TP)、可溶性总磷(DTP)、颗粒态磷(PP)。TN采用碱性过硫酸钾消煮—紫外分光光度法测定;DTN先用0.45 μm滤膜过滤,再用碱性过硫酸钾消煮—紫外分光光度法测定;采用紫外分光光度法测定;采用靛酚兰比色法测定;PN用TN和DTN差减计算获得,DON用DTN与差减计算获得。TP采用碱性过硫酸钾消煮—钼蓝比色法测定;DTP先用0.45 μm滤膜过滤,再用碱性过硫酸钾消煮—钼蓝比色法测定;PP用TP 与DTP差减计算获得。
数据分析和绘图采用Microsoft Office Excel 2010、SPSS11.5和R 3.0.1进行。差异显著分析用最小显著差异(LSD)法进行。
对稻田田面水中氮素存在形态分析(表1)可知,施氮肥8 d内(分蘖肥为7 d)田面水中的氮素以DTN为主,即使不施氮肥条件下DTN占TN的比例较低,但也高达62.4%,而其他施氮肥处理DTN占TN的比例均在88%以上,且随施氮量的增加而线性增加。田面水中PN比例较低,虽然不施氮肥的比例较高,但也仅为37.6%,而其他施氮肥处理的均低于12%,且随施氮量的增加而逐渐降低。几种可溶态氮素(和DON)的存在比例在施氮与不施氮条件下表现不同:在不施氮肥条件下,以为主,其占TN的比例为35.4%,DON和的比例分别仅为11.9%和15.1%;但在施氮条件下,则以DON和为主,二者占TN的比例分别在31.6%和44.7%以上,且随施氮量的增加而线性增加,而的比例均低于12%,且随施氮量的增加而线性降低。
进一步分析主要形态氮素浓度与施氮量关系发现,TN、DTN、和DON的浓度均随施氮量的增加而增大,而且当施氮量超过一定量后,其增幅会明显提高(图1)。对于TN,当施氮量超过287.8 kg∙hm-2后,其浓度增幅会增加2.4倍;对于DTN,施氮量超过289.9 kg∙hm-2后,其浓度增幅会增加2.5倍;对于和DON,当施氮量分别超过231.5 kg∙hm-2和336.7 kg∙hm-2时,其浓度也会明显提高(图1)。
表1 不同施氮水平下田面水中氮素形态构成Table 1 Fractionation of N in surface water relative to N application rate
图1 施氮量与田面水中不同形态氮素浓度的关系Fig.1 Relationships between N application rate and concentrations of different forms of N in surface water
施氮肥后各主要形态氮素的动态变化见图2。从图中可以看出,各施氮处理的TN和DTN浓度变化动态基本一致,无论施基肥、分蘖肥,还是施穗肥,二者均在施肥后1 d立即达到峰值,而后逐渐降低,在施基肥和分蘖肥后5 d降低至与不施氮肥基本接近,在施穗肥后2 d与不施氮肥处理基本接近。与TN和DTN相同,各施氮处理DON浓度在每次施肥后1 d就迅速达到峰值,但降低速度较TN和DTN更快,在基肥和分蘖肥后3 d降低至与不施氮处理基本接近。各施氮处理NH4+-N浓度动态变化特点与TN和DTN不同,其浓度在施基肥和分蘖肥后2 d才达到峰值,而后逐渐降低,在施肥后5 d基本与不施氮处理接近,在穗肥后1 d立即达到峰值,而后在第2天迅速降低至与不施氮基本接近。此外,从施肥期来看,施基肥和分蘖肥后各施氮处理的各形态氮素浓度均明显高于施穗肥后的浓度。比如N2处理,其TN浓度在施基肥和分蘖肥后1 d分别高达73.7 mg∙L-1和62.6 mg∙L-1,而穗肥后仅为12.0 mg∙L-1,仅为基肥和分蘖肥的16%和19%。
对磷素的存在形态分析(表2)可知,施氮肥后1~8 d内稻田田面水中磷素以PP为主,不同施磷水平下PP占TP的比例为76%~93%。田面水中DTP的比例较低,各施磷肥处理均低于24%。此外,DTP占TP的比例会随着施磷量的增加而线性增加,而PP占TP的比例则线性降低。
分析施磷肥后1~8 d内各形态磷素的浓度(图3)后发现,无论是何种形态磷,其浓度均随施磷量的增加而增加,二者之间存在极显著的线性相关关系,相关系数均大于0.97。但进一步分析发现,TP、PP浓度随施磷量增加的增幅较大,而DTP的增幅较小。
图2 不同施氮水平下田面水中不同形态氮素浓度的动态变化Fig.2 Dynamic changes of different forms of N in surface water relative to N application rate
表2 不同施磷水平下田面水中磷素形态构成Table 2 Fractionation of P in surface water relative to P application rate
施磷肥后TP、DTP和PP浓度的变化动态基本一致,均在施肥后1 d达到峰值,并在施肥后3 d内急剧降低,各处理TP、DTP、PP的平均降幅分别高达80.0%、83.6%、79.0%,而后缓慢降低,在施肥5 d后保持平稳(图4)。
图3 施磷量与田面水中不同形态磷素浓度的关系Fig.3 Relationships between P application rate and concentrations of different forms of P in surface water
图4 不同施磷水平下田面水中各形态磷素浓度的动态变化Fig.4 Dynamic changes of different forms of P in surface water relative to P application rate
目前,关于施用尿素后稻田田面水中氮素主要以何种形态存在的研究结果不尽相同。有研究表明,田面水氮素主要以为主,而且占TN的比例会随施氮量的增加而增加[10-11]。但也有研究[9]表明,施尿素后田面水中氮以DTN为主,DTN中又以DON和为主。本研究表明,施用尿素后,稻田田面水中的氮素以DTN为主,不同施氮水平下DTN占TN的比例均在88%以上,且DON和是稻田氮素主要存在形态,和DON占TN的比例分别在31.6%和44.7%以上。稻田氮素存在形态的差异,可能主要与尿素在稻田水—土系统中的水解速度有关。水解速度快,表现出以为主的特点,水解速度慢,会有大量的尿素存在于田面水中,则会表现出DON和共同为主的特点。本研究中,尿素施用后并不会立即水解为,在施肥后3 d内,田面水中DON的浓度很高(图1),表明仍有大量的尿素溶于田面水中,因而,表现出田面水氮素以DON和共同为主的特点。此外,本研究中各形态氮素占TN的比例与施氮量关系的结果与施泽升等[9]的类似,均是DTN、DON、占TN的比例随施氮量的增加而明显增大,而和PN的比例随施氮量的增加而逐渐降低。
关于施磷肥后田面水中磷的主要存在形态,已有研究的结果也不一致。有研究表明,施磷肥后,DTP是田面水中磷素的主要成分,占TP的比例介于50.0%~79.0%之间[9-10]。但也有研究表明,PP是田面水中磷素的主要存在形态,DTP/TP基本在0.5以下[12]。本研究表明,施磷肥后,田面水中磷素主要以PP为主,不同施磷水平下PP占TP的比例介于76%~93%之间,而DTP的比例均低于24%。磷肥存在形态的差别可能主要与施磷肥时对土壤的扰动程度有关,扰动较大,田面水中土壤颗粒较多,被吸附固定的磷素就多,因而PP浓度较高。本研究的磷肥均作为基肥一次使用,在施基肥后要耙田,对土壤的扰动比较剧烈,因而表现出磷素以PP为主的特点。关于各形态磷占TP比例与施磷量的关系,已有的研究结果比较一致,均为DTP占TP 的比例随施磷量增加而增加,PP占TP的比例则随施磷量的增加呈不断降低趋势[9-10]。本研究也获得了类似的结果。DTP和PP占TP的比例与施磷量的关系不一致,可能主要与PP容易沉淀有关。
诸多研究表明,氮肥施用量是影响田面水中氮素浓度的主要因子,而且氮肥用量与田面水中不同形态氮素浓度之间均有显著的线性相关关系[4,9-11]。本研究中,TN、DTN、和DON浓度虽然均随施氮量的增加而不断增大,但当施氮超过一定量后,这些氮素浓度随施氮量增加的增幅会明显变大(图1)。田面水中氮素浓度之所以出现这种跃变特点,可能与施氮量以及土壤对氮素的吸附和固定有关。当施氮量较低时,施用的氮素一部分会被土壤吸附固定,一部分溶于田面水中,呈现出田面水中氮素浓度随施氮量增加而增大的特点;但当施氮量较高,超过土壤对氮素的最大吸附固定量后,溶于田面水中的氮素就会明显增加,从而呈现出氮素浓度跃增的特点。田面水中氮素浓度跃变的特点,在以往研究中有一定的反映,如已有研究表明,稻田氨挥发量会随施氮量增加而增加,但当超过某个施氮水平后,氨挥发量会跃增[13-15]。稻田氨挥发量与田面水中铵态氮浓度有密切的相关关系[16-18],因此,氨挥发量跃增间接反映田面水中铵态氮浓度也会发生跃增。有研究表明,磷肥用量与田面水中不同形态磷素浓度之间均有显著的线性相关关系[4,9-11]。本研究结果也表明,TP、DTP和PP的浓度会随施磷量的增加而线性增加(图3)。
上述研究结果表明,尽可能降低氮、磷肥施用量,是降低江汉平原地区稻田氮、磷浓度,控制氮、磷径流和气态损失的重要措施,而且为了防止氮素损失的大幅增加,氮肥施用量应尽可能控制在231.5 kg∙hm-2以内。
关于施尿素后田面水中TN浓度动态变化趋势,以往的研究结果基本一致[4,7,11-12],即无论基肥还是追肥,均在施肥后1 d就达到峰值,而后逐渐降低。但不同的研究,TN浓度降低速度不同。金洁等[11]在浙江杭州余杭区及吴俊等[12]在苕溪流域的研究结果表明,无论基施还是追施尿素,田面水中TN浓度1周以后明显降低,控制氮素田面流失主要时期为施肥后1周内。而施泽升等[9]在洱海北部地区的研究则显示,基施和追施尿素9 d后TN才降低至较低水平,施肥后 9 d 内是控制氮素损失的关键时期。本研究中,TN也是在施尿素后1 d就达到最大值,但TN降低的速度更快,基肥和分蘖肥后5 d就降低至与不施氮肥趋同,施穗肥后2 d就降至与不施氮肥趋同(图2)。有研究表明,DTN浓度在基肥后1~3 d 内达到峰值,在穗肥后1 d就立即达到峰值,但均在9 d 后明显降低。本研究中DTN浓度的变化规律与TN相同,均在施尿素后1 d达到峰值,基肥和分蘖肥后5 d、穗肥后2 d就降至与不施氮肥趋同(图2)。本研究中TN和DTN降低速度快可能主要与本地区氨挥发较快有关。从本研究TN和DTN的动态变化规律来看,在江汉平原地区,施基肥和分蘖肥后5 d内、施穗肥后2 d内是控制稻田氮素损失的关键期。
氨挥发是稻田氮肥损失的主要途径之一[19-21],氨挥发损失的氮量占施氮量的8%~39%[22-23]。而且田面水中的浓度是氨挥发的决定因素之一。鉴于此,以往的研究对稻田田面水中动态变化的关注较多。诸多研究表明,施尿素后,田面水中浓度呈现先升高后迅速降低的特点[4,9,12]。但也有研究表明,无论基施还是追施尿素,施肥后1 d田面水中就达到峰值[10-11]。此外,关于浓度降低速度的研究结果也不尽相同。有研究表明,田面水中浓度在施尿素9 d后明显降低[9],有的则显示,施肥后4 d内就迅速降低[10]。本研究中,浓度在基肥和分蘖肥后呈先增加后降低的特点,施肥2 d后达到峰值,5 d后降低至与不施肥趋同,但穗肥后1 d就立即达到最大,2 d后迅速降低至很低水平,而且穗肥后浓度明显低于基肥和分蘖肥的浓度(图2)。可见,施基肥和分蘖肥后5 d内、施穗肥后2 d内,特别是施基肥和分蘖肥后5 d内是控制江汉平原地区稻田氮素氨挥发损失的关键期。而且,适当将基肥和分蘖肥的氮肥后移至穗肥是控制该区域氮素氨挥发损失的重要措施。
本研究中,分蘖肥和穗肥的施氮量相同,但分蘖肥后田面水中氮素浓度明显高于穗肥后的,这可能主要与施肥前土壤氮素水平的差异有关。施分蘖肥仅距施基肥8 d时间,在这8 d时间内,水稻还处于返青阶段,对氮素的吸收量很小,有较多的基肥氮残留于土壤中,这使得施分蘖肥时土壤氮素含量较高,分蘖肥氮被土壤吸附的量会比较少,因而,分蘖肥氮在田面水中存在的量较多,田面水氮素浓度较高。而施穗肥时距分蘖肥和基肥已有近两个月的时间,水稻生长也经历了从返青至分蘖盛期2个生育时期,水稻会吸收大量的氮素,在施穗肥时土壤中氮素含量相对较低,因而,穗肥氮就会有较多部分迁移至土壤中,田面水中氮素浓度就较低。
从以往研究结果来看,施磷肥后田面水中磷素的动态变化趋势比较一致,均是施磷肥后1 d各形态磷的浓度就达到峰值,而后逐渐降低[4,9-12]。但田面水中磷素降低的速度差别较大。有研究表明,磷素易被固定和沉淀,4~5 d内就迅速降至稳定[4]。但也有研究显示,施磷肥12 d后,各施磷处理的田面水磷素浓度才与不施磷的基本接近,施磷后两周内是磷流失高风险期[9]。有的研究甚至认为稻季土壤磷素流失的最大风险时期约在水稻移栽后一个月内[24]。本研究中,TP、DTP和PP浓度变化趋势与以往研究一致,但各形态磷的降速较快,施磷肥3 d后,TP、DTP和PP的平均降幅分别高达80.0%、83.6%和79.0%(图4)。从本研究结果来看,施磷肥后1~3 d是控制江汉平原稻田磷素流失的关键时期。
江汉平原是我国水稻主产区,且直播稻和机插秧稻已成为该区域最主要的两种稻作方式。这两种稻作方式普遍存在施基肥后1~2 d排放刚施肥泡田水的情况,由于人为排放刚施肥泡田水,再加上水稻种植中重施基肥,轻施或不施穗肥,因而,江汉平原水稻种植中氮、磷流失和氮素氨挥发损失均较严重。从本研究的结果来看,要防控该区域稻田氮、磷流失,应通过节水泡田等措施,力争做到不排泡田水,如果要排,也应尽可能在旋田施肥5 d后再排,严格禁止旋田施肥后3 d内排泡田水;在施分蘖肥和穗肥时,应密切关注天气状况,力求避免施肥后5 d内有大雨或者暴雨发生。此外,要防控该区域水稻氮素氨挥发损失,应重点控制施基肥和分蘖肥后5 d内的氨挥发损失,而且应尽可能将基肥和分蘖肥的氮肥后移至穗肥施用。
在江汉平原地区,施尿素后,水稻田面水中氮素以DTN为主,占TN 88%以上。在可溶态氮素中,又以DON和为主,二者占TN 76.3%以上。施磷肥后,田面水中磷素以PP为主,占TP 76%以上。减少氮、磷肥用量可降低稻田氮、磷损失,且氮肥施用量应尽可能控制在231.5 kg∙hm-2以内。施基肥和分蘖肥后5 d内、施穗肥后2 d内是控制江汉平原地区稻田氮素损失的关键期,施磷肥后3 d内是控制磷素流失的关键期。