环境科学学报  2015, Vol. 35 Issue (6): 1914-1921
小麦-玉米轮作体系农田氮素淋失特征及氮素表观平衡    [PDF全文]
丁燕1, 2, 杨宪龙1, 2, 同延安1, 2 , 韩仲宇1, 2, 谌琛1, 2, 唐希望1, 2    
1. 西北农林科技大学资源环境学院, 杨凌 712100;
2. 农业部西北植物营养与农业环境重点实验室, 杨凌 712100
摘要:连续6年采用渗漏计法研究了不同施氮处理下陕西关中小麦-玉米轮作区农田土壤90 cm深度处氮素(N)淋失特征和土壤-作物体系氮素表观平衡状况.结果表明:该地区农田氮素淋溶主要发生在降雨量较多的玉米季,且集中在8月和9月.监测期内,TN和NO3--N年平均流失量分别为2.72~23.07 kg · hm-2和1.53~18.72 kg · hm-2,年流失率分别为0.65%~3.44%和0.82%~3.32%,且年总氮、硝态氮流失量均随年施氮量增加呈指数增加.氮素淋失形态中,NO3--N比例较高,可占总氮淋失量的56.00%~81.00%,且随着氮肥用量的降低,其占总氮淋失量的比例也随之减小.可见,施氮量的大小在一定程度上会影响淋失液中各形态氮的比例.氮素表观平衡结果显示,随着施氮量提高,氮素在土壤中的残留和表观氮盈余均呈现指数增加趋势.长期施氮条件下,土壤-作物体系氮素表观损失率的幅度为32.60%~55.20%,土壤表观残留率为-0.17%~8.20%.多年监测结果表明,优化施氮模式下,作物不仅可以获得较高的产量和氮肥利用率,农田氮素淋失量也大幅降低,在节约肥料资源的同时减轻了潜在的环境风险.
关键词施氮量    NO3--N    总氮    淋溶    表观平衡    
Characteristics of N leaching and apparent N budget in cultivated lands under a Winter wheat-summer maize rotation system
DING Yan1, 2, YANG Xianlong1, 2, TONG Yan'an1, 2 , HAN Zhongyu1, 2, CHEN Chen1, 2, TANG Xiwang1, 2    
1. College of Natural Resources and Environment, Northwest A&F University, Yangling 712100;
2. Key Laboratory of Plant Nutrition and the Agri-environmental in Northwest China Ministry of Agriculture, Yangling 712100
Abstract: A six-year field experiment, using the Lysimeter method, was conducted in Guanzhong Plain of Shaanxi Province to investigate the effects of nitrogen (N) fertilizer application amount on N loss characteristics and N budget in the soil profile at 90 cm depth under winter wheat-summer maize cropping system. Results showed that N leaching mainly occurred in August and September, the rainy period of maize growing season. In the monitoring period, the average annual loss of total N and nitrate (NO3--N) were 2.72~23.07 kg · hm-2 and 1.53~18.72 kg · hm-2, respectively, and the average annual leaching rate from fertilizer for total N and NO3--N were 0.65%~3.44% and 0.82%~3.32%, respectively. N leached from the soil were mostly in the form of NO3--N comprising with 56.00%~81.00% of total N, whereas ratio of NO3--N decreased with low N fertilizers. Thus, different N fertilizer amounts affected the forms of nitrogen loss from the soil. The nitrogen budget result illustrated that the nitrogen application significantly enhanced N surplus and nitrate accumulation in the 90 cm depth of soil profile. The annual residual NO3--N and N surplus had a positive relationship with the annual N fertilizer inputs. Continuous N application resulted in N loss rate of 32.60%~55.20% and annual residual NO3--N in the 90 cm depth of soil of-0.17%~8.20%. According to the multi-year monitoring results, OPT treatment showed less nitrogen leaching amounts at high yield and better N use efficiency, reduced potential environmental risks and saved N fertilizer resources.
Key words: nitrogen application amount    NO3--N    total nitrogen    leaching    N budget    
1 引言(Introduction)

充足的氮供应是作物获得高产的必要条件,然而,自20世纪80年代以来,中国农田生态系统在长期不合理施肥下氮素总体上处于盈余状态,而且呈现持续增长趋势(Zhu and Chen, 2002).大量研究结果表明,氮肥过量施用不仅不能提高作物产量,反而使得NO3--N在土壤中大量积累淋失,从而对地下水产生造成风险(贾良良等,2001张明等,2011).农田氮素淋失作为氮肥输入土壤-作物体系后一个重要的损失途径,一直以来受到环境科学工作者的广泛关注(刘宏斌等,2004王西娜等,2006).张维理等(1995)对我国北方69个地点的地下水和饮用水硝酸盐含量的调查结果表明,由于氮素肥料的大量施用,半数以上的地下水样中硝酸盐含量超标.众所周知,土壤发生淋溶需要两个条件:N在土壤中的富集和土壤水迁移(张国梁和章申,1998李世清和李生秀,2000王朝辉等,2006巨晓棠等,2003aLiang and Mackenzie, 1994),而降雨作为不易控制因素,所以研究氮肥施入土壤后的转化和去向,确定氮肥施用量成为现在的研究热点(刘新宇等,2010巨晓棠等,2002).

冬小麦-夏玉米轮作是陕西关中平原的主要种植模式,在粮食生产中发挥重要作用.研究报道,陕西省部分农田存在严重施氮过量问题,使得土壤中硝酸盐大量累积,不仅增加了氮素淋溶的风险,且表观损失量大幅提高(贾良良等,2001同延安等,2004).目前,关于农田氮素淋失国内外开展了大量的研究,然而大部分研究较多关注于NO3--N,对其他氮素淋失形态(如有机氮等)关注较少.事实上,有研究显示,有机氮也是氮素从土壤生态系统流失的主要形态之一(Hagedorn et al., 2001Perakis and Hedin, 20012002).故在陕西关中平原地区利用渗漏计法连续6年展开长期定位试验,研究小麦-玉米轮作下N淋失特征和氮素表观平衡,为陕西关中地区合理施肥提供理论依据.

2 材料与方法(Material and methods) 2.1 试验地概况

试验于2007年10月—2013年10月在陕西省杨凌区大寨乡国家黄土肥力与肥料效益监测基地进行,地处陕西关中平原西部.试验地所在经纬度34°4′N、108°2′E,海拔534 m,属于暖温带半湿润偏旱气候,年平均气温13 ℃,≥10 ℃,积温4196.2 ℃.年平均降水量550~600 mm,主要集中在7—10月份,年均蒸发量993 mm,无霜期184~216 d.试验土壤属褐土类、塿土亚类、红油土属、厚层红油土种,土壤质地为粉砂粘壤土,地下水位28 m,深井抽水灌溉.实行冬小麦-夏玉米轮作的一年两熟制.播前0~100 cm土壤理化性状见表 1,生育期内降雨见图 1.

表 1 试验土壤的基本理化性状 Table 1 Physical and chemical properties of tested soil

图 1 2007年10月—2013年10月小麦-玉米生育期内降雨分布 Fig.1 Rainfall distribution during the growth period of wheat and maize in 2007—2013
2.2 实验设计

试验设4个处理,年施氮量分别为:CK不施肥(纯N:0 kg · hm-2)、OPT-N减量施氮(纯N:165 kg · hm-2)、OPT优化施肥(纯N:330 kg · hm-2)、OPT+N增量施氮(纯N:495 kg · hm-2)(施肥方案见表 2).其中优化施肥是根据资料调研、专家推荐关中地区施肥量上设定的,每个处理重复3次,共计12个小区,小区面积为5 m×6 m=30 m2,随机区组排列.种植方式为冬小麦-夏玉米轮作.在每个小区按照每20 cm一层挖至 90 cm土体深处安装 2 个规格为40 cm×50 cm的PVC管材淋溶盘后分层回填,通过侧向管连接至集液管,以收集淋溶水样,通过真空泵抽取淋溶液测定(图 2),供试小麦品种为小偃22,玉米品种为郑单958.冬小麦于前一年10月撒施播种,次年6月收获,其中50%的氮肥和100%的磷肥在播种前撒施翻耕入土,另50%的氮肥在返青后拔节前追施.夏玉米于6月穴施播种,10月收获,其中50%的氮肥和100%的磷钾肥在五叶期穴施,另50%的氮肥在大喇叭口期追施.氮肥用尿素(含N 46%),磷肥用普通过磷酸钙(含P2O512%),钾肥用氯化钾(含K2O 60%),除施肥操作外,灌溉及其他田间管理措施均同当地农民习惯保持一致.

表 2 试验施肥方案 Table 2 Fertilizer application rate in the experiment

图 2 田间原位淋溶装置示意图 Fig.2 Sketch map of field device in situ leaching
2.3 样品采集与分析 2.3.1 植株样品采集与分析

在小麦收获时,各小区分别采8行样品,风干以计算其干重(样品分籽粒和秸秆两部分分析),然后以采样面积折算生物量.玉米收获时,每个小区全部收割(籽粒和秸秆分开),分别称量各小区秸秆和籽粒鲜重,然后各小区取20棒玉米籽粒和部分秸秆,带回实验室风干计算其干重,以小区总鲜重折算其生物量.四分法取部分样品烘干,烘干样品粉碎后,用浓H2SO4-H2O2消解,然后用流动注射分析仪测定籽粒和秸秆全氮.

2.3.2 土壤样品采集与分析

玉米收获后,用土钻取0~90 cm土样,每个小区避开淋溶盘取5个样点混合均匀,每20 cm一层,最后80~90 cm分作一层,新鲜土样带回实验室后立即用CaCl2浸提后用流动注射分析仪测定NO3--N.

2.3.3 灌溉与干湿沉降样品采集与分析

每次灌溉当天收集灌溉水样品,记录灌溉量后带回实验室低温保存待测.干湿沉降样品采集利用APS-2A型降水降尘自动采样器收集,量取降雨量后低温保存待测,水样NO3--N、NH+4-N用clever chem200流动注射分析仪测定.

2.3.4 淋溶水样采集与分析

每次降雨灌溉4~5 d后收集淋溶水样,量取体积后取100 mL样品,带回实验室放入冰箱冷冻保存待测,NO3--N测定方法同灌溉水.总氮(TN)用碱性过硫酸钾氧化-紫外分光光度计比色法测定.

2.4 计算方法及公式

采用Excel软件对试验数据进行统计分析,用PASW Statistics(SPSS)18.0进行Duncan新复极差法多重比较和利用Sigma plot作图,分析过程中所用公式如下:

TN(NO3--N)年表观淋洗量(kg · hm-2)=施氮处理TN(NO3--N)淋溶量-空白处理TN(NO3--N)淋溶量 (1)

TN(NO3--N)年淋失率=TN(NO3--N)年表观淋洗量×100%/年施氮量 (2)

NO3--N累积量(kg · hm-2)=土层厚度×土壤容重×NO3--N含量/10 (3)

土壤氮素矿化量(kg · hm-2)=不施氮肥区地上部分氮积累量+不施氮肥区土壤残留无机氮量-不施氮肥区土壤起始无机氮量-灌溉水无机氮量-干湿沉降无机氮量 (4)

氮素表观损失量(kg · hm-2)=氮输入量-作物吸收氮量-土壤残留无机氮量 (5)

氮盈余量(kg · hm-2)=残留无机氮量+表观损失氮量 (6)

氮肥表观利用率=(施氮区植株吸氮量-对照区植株吸氮量)×100%/施氮量(7)

氮肥表观残留率=(施氮区土壤残留Nmin-对照区土壤残留Nmin)×100%/施氮量 (8)

氮肥表观损失率=100%-氮肥表观利用率-氮肥表观残留率 (9)

3 结果与分析(Results and analysis) 3.1 2008—2013年小麦-玉米农田90 cm处渗漏液量变化

试验从2007年10月进行到2013年10月,各处理小麦-玉米生育期内在土层90 cm处平均渗漏液量如图 3所示,监测周期内,共发生了17次淋溶,年产流频率在2到5次之间.2008年2次、2009年前两次,以及2010年7月、2012年8月共6次产流均由灌溉引起外,11次都是由大的降雨或者连续降雨造成,当中只有1次发生在5月,其余10次均发生在8—9月份.监测周期内2011年9月发生的超过300 mm强降雨造成两次产流量之和超过80 mm,但由于夏季温度较高,玉米叶蒸腾作用强,导致强降雨之后单次产流并未达到平常降雨产流的4~5倍.另外,发生连续降雨也容易产生产流,如2010年7月到9月和2012年的8月到9月,土壤中若含水率一直处在较高的状态,其后的较小的降雨都将产生淋溶.另外,统计结果表明,降雨量、灌溉量与产流并无显著相关性,说明,在降雨量高的周期内,产流量同时也受到土壤中干湿状况、作物蒸腾、作物生长状况等因素的影响,这也解释了强降雨不一定导致高产流量,反而较小的降雨也能产生产流的原因.

图 3 2008年到2013年间小麦-玉米轮作农田渗漏液的动态变化 Fig.3 Dynamic change of leachate in wheat-maize rotation field during the year 2008—2013
3.2 不同施氮对N淋失形态的影响

2008—2013年间监测时间内共发生的17次淋溶(2008年4月6日渗漏液未测总氮浓度),其NO3--N、TN的浓度及淋失量动态变化见图 4(由于在干旱区域硝化作用强烈,且容易吸附,NH+4-N浓度偏低,本文中不考虑此形态的流失).从渗漏液中两种形态氮浓度和流失量动态变化上看,浓度和流失量变化趋势具有高度一致性,NO3--N与TN的流失也同样如此,NO3--N、TN浓度范围分别为0~103.5 mg · L-1和0~120.3 mg · L-1.NO3--N、TN流失量范围分别为0~18.01 kg · hm-2和0~18.41 kg · hm-2,其都随着施氮量的增加而呈现增加趋势,CK和OPT-N施肥处理前期流失TN和NO3--N浓度较高,后期逐渐平稳,在其他处理出现峰值时仍然很低,而OPT+N施肥处理较其他施肥处理显著增加了总氮和硝态氮浓度以及流失量.

图 4 008—2013年监测期内TN、NO3--N浓度及流失量的动态变化 Fig.4 Dynamics of total N,nitrate concentration and leaching amount during the year of 2008—2013
3.3 不同施氮对TN和NO3--N的年流失量和流失率的影响

不同施氮处理对TN和NO3--N的年流失量和流失率如表 3所示.各年间TN和NO3--N各处理间流失量和流失率变化较大,其TN从2009年到2013年平均流失量为2.72~23.07 kg · hm-2,平均淋失率为0.65%~3.44%.NO3--N在2008年到2013年平均流失量为1.53~18.72 kg · hm-2,平均淋失率为0.82%~3.32%,且TN和NO3--N淋失量和淋失率均随着施氮量的提高而呈现增加趋势.统计结果表明TN和NO3--N年流失量与施氮量存在指数相关,其TN年流失量y和NO3--N年流失量z与施氮量x的回归关系分别是y=2.266e0.0044x(R2=0.76* *);z=1.598e0.0048x(R2=0.70* *),可见随着施氮量的提高,其淋失TN和NO3--N的比例也随之提高,大大增加了其淋失量.2009—2013年的监测期内,NO3--N是土壤中氮素淋失的主要形态,其占TN流失量的比例为56.00%~81.00%,除去占有很少比例的NH+4-N外,其余是有机氮,高忠霞等(2010)利用大型回填土渗漏池研究得出淋溶水样中除无机氮外基本以可溶性有机氮为主.长期监测结果得出NO3--N占TN流失量比例从小到大排列顺序为CK<OPT-N<OPT<OPT+N,不施氮处理最小为56%,表明不同施氮处理会影响氮素淋失形态.

表 3 不同施氮对TN、NO3--N年流失量和流失率的影响 Table 3 Effects of different N application amounts on the amount and ratio of total nitrigon and nitrate nitrigon leaching
3.4 不同施氮对N来源与去向的研究

从2008年到2013年,每一年的氮素输入和输出如表 4所示,计算氮素平衡时,一般按照0~100 cm作为作物根际吸收氮肥主要吸收范围(巨晓棠等,2002),但由于淋溶盘设置在90 cm深度,故按照0~90cm土层计算.从结果可以看出,氮的来源包括肥料氮、灌溉和干湿沉降携入氮、矿化氮以及土壤中起始NO3--N(塿土主要以NO3--N为主,铵态氮含量极低,不计算在内),氮矿化不考虑氮的激发效应,不同处理视作相同.氮的输出主要包括作物携出氮、土壤残留氮和其他形式损失,在肥料氮的损失途径中,主要有氨挥发,硝化-反硝化和淋洗到根层以下或发生地表径流损失.在关中平原小麦-玉米轮作地,基本不存在地表径流,且由施肥引起的氮表观淋失只占氮表观损失的2.40%~6.00%,所以气态损失有可能是该地区氮主要损失途径,梁东丽等(2002)在塿土上的研究也表明,在玉米生长季节反硝化损失最高.从氮的输出结果看出,OPT+N处理吸氮量比OPT处理增幅小,但土壤残留氮、表观无机氮淋洗量、表观氮损失量和氮盈余量在OPT+N下明显提高.从结果可以得到土壤中残留Nmin(y)和盈余氮(z)均随施氮量(x)增加呈现增加,其回归关系分别为y=29.45e0.0038x(R2=0.80* *)和z=35.75e0.0054x(R2=0.88* *),采用F值检验法进行显著性分析,得到两者均达到1%极显著性水平,可见利用这样的关系估算本地区土壤中氮残留和氮盈余是可行的.从当年看,除去损失的氮素,残留土壤氮也是盈余氮的主要去向,且土壤中残留与作物生长有密切关系,当吸氮量较高时,土壤中残余氮随之降低.氮表观损失不仅与施氮量有关,也与当年作物生长状况有很大关系,如2010年、2011年和2012年吸氮能力强,氮表观损失普遍有所降低.反之,2008年、2009年和2013年则比较高.而在高施氮下,高吸氮量对氮表观损失的降低并不十分明显.

表 4 2008—2013年氮的来源与去向 Table 4 Input and output of nitrogen in 2008—2013
3.5 小麦-玉米轮作体系中的N素平衡

将6年内所有氮素输入和输出作为一个整体,计算氮平衡,如表 5所示是2008—2013年小麦-玉米轮作系统中氮输入和输出平衡,除了对照处理外,氮肥是氮的主要来源,占52.00%~77.00%,降雨和灌溉携入氮、矿化氮和起始无机氮之和为903.4 kg · hm-2,占氮来源的23.00%~48.00%.OPT+N施肥较OPT施肥,吸氮量增幅基本不变,而土壤残留、表观NO3--N淋失和表观损失大幅度增加,增加幅度分别为128.00%、122.00%和95.00%.氮肥表观利用率由53.50%降低到36.60%,表观残留率由4.00%升高到8.20%,表观损失率则由42.50%升高到55.20%.6年试验后OPT+N处理下其盈余氮约占总投入氮50.00%,加之继续如此高的肥料投入,最终对环境造成损害.OPT-N处理与OPT处理相比,土壤中氮残留和表观氮损失有很大幅度的降低,但与此同时作物的吸氮量也大幅度下降,且经过6年OPT-N施肥后,土壤供氮由起始46.2 kg · hm-2降低到32.1 kg · hm-2,继续减氮施肥,土壤必定会供氮不足,影响产量.而OPT施肥处理在吸氮量上与OPT+N施肥处理大致相同,但在降低氮盈余上减少了约50.00%.

表 5 2008—2013年小麦-玉米轮作系统氮素表观平衡 Table 5 Nitrogen budgets of the wheat-maize rotation system during the year 2008—2013
4 讨论(Disscusion)

经过陕西关中平原的长达6年渗漏池试验,土层较为接近原状土,其结果代表了该地区淋溶的一般性规律.可以发现:降雨和灌溉都有可能导致淋溶的发生,小麦-玉米轮作农田上每年发生的频率在2~5次左右,基本都集中在降雨量多的玉米季,即8、9月份,小麦季基本不发生淋溶.而在渗漏液中,以NO3--N为主的无机氮占TN的56.00%~81.00%,除去低含量的NH+4-N,其余以可溶性有机氮为主的有机氮,同样是不可忽视的N淋失形态,李宗新等(2007)也提出小麦-玉米轮作年周期内农田土壤发生淋溶的氮素以NO3--N为主,可溶性有机氮占30%左右,这是由于NO3--N不易被土壤颗粒和土壤胶体吸附(余贵芬等,1999),易于遭雨水或灌溉水淋洗迅速渗漏.且低氮处理下,其NO3--N含量占TN含量低,不同施肥处理会影响氮素淋失形态.N的流失伴随着水分移动,所以降雨致使淋溶的发生(Moreno et al., 1996),但渗漏液量和引起产流的降雨灌溉量并未表现出显著相关,这是由于大田条件下淋溶影响因素较多,不仅包括降雨,同时作物生长状况、土壤中干湿状况,蒸腾作用等同样影响水分下渗,这也造成了不同年份陕西关中区域淋溶差异.

从6年长期试验看出渗漏液中TN及NO3--N的浓度和淋失量均随着施氮量的增加而呈现增加趋势,尤其在OPT+N处理下,NO3--N的浓度高达103.5 mg · L-1,是国家标准地下水NO3--N质量标准(最高限值为20 mg · L-1)(GB/T14848-93)的5倍,长此以往必定对地下水饮用水造成健康风险.NO3--N年平均淋失量为1.53~18.72 kg · hm-2,由施肥引起的肥料NO3--N流失量为1.35~16.43 kg · hm-2,相当于年流失尿素2.93~35.72 kg · hm-2,与前人在塿土上试验结果类似(高忠霞等,2010).

巨晓棠等(2003b)对北京郊区氮去向研究表明,氮肥利用率随施氮量的升高而降低,而损失率相应的增加,长期施氮下,氮肥表观利用率为36.60%~68.10%,氮肥表观损失率为32.70%~57.10%,残留率在-1.70%~8.20%,只占很小一部分,说明时间跨度越长,在土壤中间的残留并不按比例增加,且有可能被后茬作物所利用,所以以损失方式脱离作物根系是N素的主要去向.OPT+N处理(年施氮量为495 kg · hm-2)经过了6年定位试验,其氮肥表观损失率大于氮肥表观利用率,将造成大量的氮损失于土壤-作物这一系统,对环境造成严重氮负荷.

关于淋失,一般认为NO3--N移动到根层(作物根区层的土层深度为100 cm)以外便称作淋失,本文因淋溶盘设定位置将之定义在0~90 cm.也有学者(王西娜等,2006)将根际范围定义为0~200 cm,但刘学军等(2004)通过减氮施肥试验得出100~200 cm层次中累积的NO3--N却不因减量施N而下降,说明这一土层的硝态氮难以被作物吸收利用,即使根系吸收氮肥区域会在90 cm之下,但总体上还是非常有限的,最终还是会脱离土体进入到水环境.对于当季作物,氮的主要损失去向是残留在土体中(杨宪龙等,2013巨晓棠等,2002),但对于长期施氮,由于残留在土壤中的氮仍可被吸收利用,所以盈余氮最终便是以淋失或反硝化等其他形式损失到根系范围之外.本试验在大田开展了6年,得出了N素淋失特征和表观平衡特征,但在区域尺度上准确估算地下水中农田TN和硝态氮排放量仍然需要进一步研究.

5 结论(Conclusions)

以小麦-玉米轮作为主要种植模式的陕西关中平原,淋溶主要发生在降雨量较大的玉米季,以8、9月份居多,TN、NO3--N浓度和流失量均随着施氮量的提高而增加,且NO3--N为主要淋失形态,占TN的56.00%~81.00%.高施氮量会增加N的淋失和土壤中NO3--N的残留,不仅给农田带来土壤环境氮素负荷,而且对地下水带来健康风险.多年监测结果表明,3种施氮量下,优化施肥较OPT-N处理不仅能保证作物吸氮量和土壤供氮能力,获得较高的产量,且与OPT+N处理相比能够提高氮肥利用率,降低表观损失.综合来看,OPT是节约肥料资源和减轻潜在的环境风险的最佳施氮量.

参考文献
[1] 高忠 霞, 杨学云, 周建斌, 等. 2010. 小麦-玉米轮作期间不同施肥处理氮素的淋溶形态及数量[J]. 农业环境科学学报, 29(8): 1624-1632
[2] Hagedorn F, Bucher J B, Schleppi P. 2001. Contrasting dynamics of dissolved inorganic and organic nitrogen in soil and surface waters of forested catchments with Gleysols [J]. Geoderma, 100(1/2): 173-192
[3] 贾良良, 陈新平, 张福锁, 等. 2001. 北京市冬小麦氮肥适宜用量评价方法的研究[J]. 中国农业大学学报, 6(3): 67-73
[4] 巨晓棠, 刘学军, 邹国元, 等. 2002. 冬小麦/夏玉米轮作体系中氮素的损失途径分析[J]. 中国农业科学, 35(12): 1493-1499
[5] 巨晓棠, 刘学军, 张福锁. 2003a. 冬小麦/夏玉米轮作中NO\|3-N在土壤剖面的累积及移动[J]. 土壤学报, 40(4): 538-546
[6] 巨晓棠, 潘家荣, 刘学军, 等. 2003b. 北京郊区冬小麦/夏玉米轮作体系中氮肥去向研究[J]. 植物营养与肥料学报, 9(3): 264-270
[7] 李世清, 李生秀. 2000. 半干旱地区农田生态系统中硝态氮的淋失[J]. 应用生态学报, 11(2): 240-242
[8] 李宗新, 董树亭, 王空军, 等. 2007. 不同肥料运筹对夏玉米田间土壤氮素淋溶与挥发影响的原位研究[J]. 植物营养与肥料学报, 13(6): 998-1005
[9] Liang B C, Mackenzie A F. 1994. Changes of soil nitrate-nitrogen and denitrification as affected by nitrogen fertilizer on two Quebec soils[J]. Journal of Environmental Quality, 23(3): 521-525
[10] 梁东丽, 同延安, Emterdy O, 等. 2002. 黄土性土壤剖面中N2O排放的研究初报[J]. 土壤学报, 39(6): 802-809
[11] 刘宏斌, 李志宏, 张维理, 等. 2004. 露地栽培条件下大白菜氮肥利用率与硝态氮淋溶损失研究[J]. 植物营养与肥料学报, 10(3): 286-291
[12] 刘学军, 巨晓棠, 张福锁. 2004. 减量施氮对冬小麦-夏玉米种植体系中氮利用与平衡的影响[J]. 应用生态学报, 15(3): 458-462
[13] 刘新宇, 巨晓棠, 张丽娟, 等. 2010. 不同施氮水平对冬小麦季化肥氮去向及土壤氮素平衡的影响[J]. 植物营养与肥料学报, 16(2): 296-303
[14] Moreno F, Cayuela J A, Fernandez J E, et al. 1996. Water balance and nitrate leaching in an irrigated maize crop in SW Spain[J]. Agricultural Water Management, 32(1): 71-83
[15] Perakis S S, Hedin L O. 2001. Fluxes and fates of nitrogen in soil of an unpolluted old-growth temperate forest, southern Chile[J]. Ecology, 82(8): 2245-2260
[16] Perakis S S, Hedin L O. 2002. Nitrogen loss from unpolluted South American forests mainly via dissolved organic compounds [J]. Nature, 415(6870): 416-419
[17] 同延安, Emteryd O, 张树兰, 等. 2004. 陕西省氮肥过量施用现状评价[J]. 中国农业科学, 37(8): 1239-1244
[18] 王西娜, 王朝辉, 李生秀. 2006. 黄土高原旱地冬小麦/夏玉米轮作体系土壤的氮素平衡[J]. 植物营养与肥料学报, 38(6): 759-764
[19] 王朝辉, 李生秀, 王西娜, 等. 2006. 旱地土壤硝态氮残留淋溶及影响因素研究[J]. 土壤, 38(6): 676-681
[20] 杨宪龙, 路永莉, 同延安, 等. 2013. 施氮和秸秆还田对小麦-玉米轮作农田硝态氮淋溶的影响[J]. 土壤学报, 50(3): 564-573
[21] 余贵芬, 吴泓涛, 魏永胜, 等. 1999. 氮在紫色土中的移动和水稻氮素利用率的研究[J]. 植物营养与肥料学报, 5(4): 316-320
[22] 张国梁, 章申. 1998. 农田氮素淋失研究进展[J]. 土壤, 30(6): 291-297
[23] 张明, 同延安, 郭俊炜, 等. 2011. 陕西关中小麦/玉米轮作区氮肥用量及施氮现状评估[J]. 西北农林科技大学学报(自然科学版), 39(4): 152-158
[24] 张维理, 田哲旭, 张宁, 等. 1995. 我国北方农用氮肥造成地下水硝酸盐污染的调查[J]. 植物营养与肥料学报, 1(2): 80-87
[25] Zhu Z L, Chen D L. 2002. Nitrogen fertilizer use in China Contributions to food production, impacts on the environment and best management strategies[J]. Nutrient Cycling in Agroecosystems, 63(2/3): 117-127