褐土区农田土壤氮磷淋溶特征及其管理措施*

郭胜利，张树兰，党廷辉，郭李萍，李丽君，高鹏程，王 蕊

(1.西北农林科技大学水土保持研究所 杨凌 712100； 2.西北农林科技大学资源环境学院 杨凌 712100； 3.中国农业科学院农业环境与可持续发展研究所 北京 100081； 4.山西农业大学资源环境学院/山西省土壤环境与养分资源重点实验室 太原 030031)

氮磷化肥投入是提高粮食产量，保障粮食安全的重要措施。目前，很多作物体系氮磷养分的输入早已超过了作物需求量。据估计，作物生产体系中全球氮肥利用率已从1950年的68%降低到2010年的47%，磷肥利用率在1950—2000年间从60%降低到44%[1-2]。1940—2008年间土壤氮残留量稳定在每年50 kg·hm-2(1.1～229 kg·hm-2)[3]； 在氮肥施入量较高的区域，土壤氮含量已呈现饱和状态[4-5]。农田土壤磷的净残留量从1950年的每年1 kg 增加到2000年的每年8 Mt，1970—2010年间磷累积约210 Mt[6]。自20 世纪80年代以来，我国农田氮肥、磷肥的投入大幅度升高。据统计，2010年我国农田氮肥消费量27.9 Mt，磷肥消费量5.3 Mt[7]。目前，我国农田整体上处于富磷状态，从1980年土壤盈余4.6 kg·hm-2提高到2012年42.1 kg·hm-2[8]。盈余的氮磷养分进入水体(河流和地下水)或以气体形式排放到大气，造成水体富营养化、饮用水硝酸盐污染、生物多样性丧失、加剧气候变化等诸多负面影响[1，9-12]。如何在保障粮食安全的同时降低环境危害日益引起全球的关注[13-15]。

在欧美等发达国家，由于农业集约化发展早于中国，自20 世纪50—80年代农业面源污染问题日益严重。而自20 世纪80年代末以来欧美发达国家开始重视对农业面源污染的研究和治理。此后，欧盟国家氮磷化肥用量分别下降约30%和50%。地下水硝酸盐污染有所缓解，湖泊和近海域水体富营养化也得到一定程度改善[16]。反观中国，自1970s 以后，各大湖泊、重要水域的水体污染，特别是水体的氮磷富营养化问题急剧恶化； 由于20 世纪80年代大量化肥的投入，90年代以后我国农业面源污染日益加重，在北方集约化农区主要表现为氮素淋溶及地下水硝酸盐污染[17-18]。褐土区主要分布于关中盆地、晋南盆地和豫西，为我国传统的粮食种植区。早在1990s 褐土区农田土壤硝态氮的淋失问题就引起了学者的关注[19]，并就农业土壤中的深层硝态氮积累与施肥、作物、降水、耕作措施以及土壤性状的关系进行了系统评价[20]。与此同时，褐土地区氮肥施用对环境污染的影响[21]及地下水氮素污染逐步引起关注[22-24]。

褐土发育于富含石灰的黄土母质上，土层深厚，土壤黏化过程明显，土体构型中黏化层、淋溶淀积层在水盐运移过程中的阻控作用明显。与同纬度华北平原的潮土相比，褐土区降水偏少，残留土壤的硝态氮会产生显著积累。近年来该区域除传统粮食作物外，菜田、果园发展迅速，过量施肥和粗放灌溉现象普遍，土壤剖面中氮磷残留问题日益严重[25-28]，大量残留在土壤深层的硝态氮对地下水的潜在威胁也日益严重[29]。褐土区由于土壤碳酸盐含量高，对磷肥的固定作用强烈，这些特性在磷肥大量施用的20 世纪80年代显著影响了磷素在土壤中的吸附、转化和移动[30-34]。但随着磷肥40 余年的持续投入，褐土区农田土壤磷的饱和吸附度逐渐升高，耕层土壤中磷素移动性逐渐升高[35-36]。另外，汾渭河谷的褐土区既是农田集中区，也是大中城市、村镇的分布地。尽管已有研究发现，地下水硝酸盐与大量氮肥施用密切相关，但生活污水及动物源的硝态氮来源同样不容忽视[37]。明确地下水不同来源硝酸盐所占比重是制定其消减的基础。基于此，本研究重点分析褐土区农田土壤氮磷的盈余状况，辨识农田氮磷的淋溶特点，评价田间管理措施对氮磷淋溶的影响，对该区域绿色发展具有重要的理论价值和现实意义。

1 材料与方法

1.1 褐土区农田土壤氮、磷的盈余状况

关中盆地位于陕西省中部，西起宝鸡，东至潼关，南依秦岭，北抵北山。东西长约330 km，南北宽30～100 km，面积约2×104km2。属暖温带半干旱、半湿润气候区。平均气温12.0～13.6 ℃，年均降水量550～750 mm，多年平均蒸发量1000～1200 mm。关中盆地集中了陕西省60%以上的人口和80%的工业以及52%的耕地，为汾渭河谷典型的粮食生产基地。本研究以关中盆地为对象，从氮、磷化肥投入，土壤氮磷残留，地下水中氮磷含量及其作物产量入手，在区域尺度上深入分析褐土区农田土壤氮磷的盈余特征。

数据来源： 1)作物、土壤数据。包括关中盆地3000 多农户氮肥和磷肥投入量、作物[冬小麦(Triticum aestivum)、玉米(Zea mays)]产量调查数据(2012年)； 土壤养分调查数据，采用网格法，每1 km2设置1 个采样点，对关中盆地5 市(西安市、咸阳市、宝鸡市、渭南市与韩城市)32 个区县的农田土壤耕层(0～20 cm)采样，利用GPS 测定各采样点经纬度，利用常规分析方法测定各样点氮磷养分含量。2)地下水中氮磷含量数据。2016—2017年在陕西关中盆地采集地下水样品213 个。依据关中盆地的土地利用结构，采样区包括蔬菜种植区、粮食作物种植区、居民生活区(20～30年前为农田)，每个区设置3 个对照区(当地地下水背景，未受人类活动影响)。3)土壤深层剖面硝态氮数据。采集23 个典型关中盆地农田0～600 cm 土壤剖面样品用于测定土体中硝态氮的含量分布。

1.2 褐土区农田生态系统土壤氮磷淋失机制

1)施肥量、作物产量和氮磷淋失量耦合关系。施肥促进作物增产，但施肥过量会导致农田土壤氮磷淋失，为更好地描述施肥量、作物产量和氮磷淋失量的关系，本研究采用Johnson 等[38]1995年提出的“安全施肥区间(safety zone)”的概念。施氮量达到作物的最高产量后，土壤剖面中并没有导致矿质氮积累，但当进一步提高施氮量至一定数值时，才会导致土壤剖面中矿质氮的显著积累。从作物达到最高产量到土壤剖面中矿质氮显著积累这一施肥区间，称为“安全施肥区间”。这一概念为深入分析施肥量、作物产量和氮磷淋失间的关系提供了思路。我们依托1984年设立于陕西省渭北旱塬的长期肥料试验[39]，选取其中10 个氮磷水平处理，结合每年监测的作物产量和地上部生物量及定期开展的耕层土壤氮磷性状变化监测数据、0～300 cm 土壤剖面硝态氮含量数据(2000年)，深入分析农田生态系统施肥量、作物产量、氮磷利用效率、硝态氮淋失量或磷积累量与环境安全施肥量之间的关系。

2)土壤微生物对硝态氮的阻控作用。农田土壤硝态氮淋溶、氧化亚氮及氮气的排放和氨挥发等氮素的转化过程都与土壤微生物有关，其中由土壤微生物驱动的硝化和反硝化作用是土壤氮素循环的两个关键过程。探明土壤硝化和反硝化微生物在阻控土壤硝态氮淋溶中的潜在作用及其对土壤硝态氮淋溶的阻控效果，会更有利于生物阻控措施形成和推广。依托秸秆试验(详见1.3 节)，我们探索了春玉米种植体系中土壤微生物变化对硝态氮迁移和积累的影响。

1.3 褐土区农田土壤氮磷淋失阻控措施和影响因素

针对褐土区的农田和菜地，在陕西、山西和河北共设4 个试验。轮作(冬小麦-夏玉米)种植体系，依托西北农林科技大学“国家黄土肥力与肥料效益长期监测基地”的长期定位试验(杨凌)进行； 单作(春玉米)种植体系依托陕西长武农田生态系统国家野外科学观测研究站(长武)进行； 露地蔬菜以“春黄瓜(Cucumis sativus)/番茄(Lycopersicon esculentum)/甘蓝(Brassica oleracea)-大白菜(Brassica rapa pekinensis)”种植体系为对象，依托河北保定清苑(潮)褐土区典型露地菜地进行； 设施蔬菜以西葫芦(Cucurbita pepo)种植体系为研究对象，依托山西省农业科学院设施蔬菜地进行(表1)。

2 结果与分析

2.1 典型褐土区氮磷投入量、土壤氮磷含量与地下水氮磷含量

关中盆地大量田间调查数据显示，无论玉米种植系统还是冬小麦种植系统，随着施肥量的提高，作物并没有呈现增产效应。在玉米季，氮肥施用量大于250 kg(N)·hm−2后，产量反而呈现降低趋势； 在小麦季，超过200 kg(N)·hm−2施氮量后，产量开始降低，磷肥超过25 kg(P2O5)·hm−2水平时，产量已经无响应(图1)。这些结果表明，褐土区施肥的产量效应已经显著不同于化肥开始大量投入的20 世纪80年代初。以往的大量田间结果表明，包括褐土区农田随着氮磷化肥投入的提高，作物产量逐步提高直至达到最高产。但经过近40年的持续化肥投入后，农田土壤产量效应明显降低。农田土壤中大量的氮磷残留间接佐证了这一点。在关中盆地调查农田0～600 cm 土壤剖面中硝态氮含量分布显示，该区域农田土壤中最高硝态氮残留量约300 kg(N)·hm−2，其中75%残留在100 cm 土层以下。这些现象表明，该地区农田/蔬菜地中存在大量氮肥残留，而且大量残留的硝态氮有向土壤深层迁移的现象(图2)。相应地，关中盆地80%耕层土壤有效磷(Olsen-P) 已超过20 mg(P)·kg−1，显著不同于20 世纪80年代初的土壤缺磷状况[当时大部分耕层土壤Olsen-P 含量低于5 mg(P)·kg−1]。这些结果表明，经过40 多年的化肥投入，关中盆地农田土壤施肥过量、土壤中氮磷残留和积累现象严重。相应地，对2017年采集的213 个地下水样本氮含量的分析发现，关中地区地下水硝态氮含量超过11 mg·L−1(WHO 饮用水标准)的样点接近15%，并且存在显著的空间分布特征(图3)。超标地区主要分布在渭南北、蒲城南、临潼及渭河流域中上游一带。地下水氮素含量高值地区与采样点化肥用量高且地下水位较高有关，而地下水位相对较低或施肥量较低的其他采样点化肥的影响明显较小。与硝态氮含量分布不同，大部分取样点的磷含量较低[低于0.05 mg(P)·L−1]，但也发现蓝田县等局部的高值区域，这可能与这些区域长期有机肥或畜牧养殖废弃物的施用有关。

2.2 典型褐土区施肥量、产量与环境安全施肥量间的耦合阶段

基于15年长期定位试验的观测资料分析发现(图4)，长期不施肥下冬小麦年平均产量为1.2 t·hm−2，地上部氮素吸收量大致相当于20 kg(N)·hm−2，此时土壤剖面硝态氮残留量低于5 kg(N)·hm−2。随着氮肥投入量增加到90 kg(N)·hm−2，冬小麦产量显著提高2 倍，作物地上部氮素吸收量显著提高，约为60 kg(N)·hm−2，此时土壤剖面硝态氮残留量并没有发生显著积累。但当氮肥的投入量为135 kg(N)·hm−2时，每千克氮肥的增产量则显著降低，地上部吸氮量也未显著提高，而土壤硝态氮残留量则显著提高，为 26 kg(N)·hm−2； 若施氮量进一步提高到180 kg(N)·hm−2，每千克氮肥的增产量则显著降低，地上部的吸收量较135 kg(N)·hm−2仅提高3%，此时土壤中硝态氮残留量进一步提高，达36 kg(N)·hm−2(图4)。氮肥的投入量、作物产量/作物吸收量和土壤剖面的硝态氮残留量，三者之间耦合关系存在着3 个阶段： 1)环境友好-资源高效阶段，这一阶段土壤矿质氮残留少，氮肥生产率高； 2)环境低风险-资源有效阶段，该阶段土壤残留氮的开始随施氮量升高显著增加，氮肥生产率开始降低； 3)环境高风险-资源无效阶段，该阶段土壤残留氮随施氮量升高显著增加，同时氮肥生产率出现负增长。

与氮肥的产量和环境效应不同，在长期试验中，尽管同样可以看到产量随施磷量提高而提高，但当磷肥水平超过较低的阈值后，进一步加大磷肥投入并不显著提高作物产量，反而会导致土壤耕层有效磷随施磷量增加而积累(图4)。同样的，可以得到土壤有效磷水平、作物产量/作物磷吸收量和水溶性磷三者耦合关系存在3 个发展阶段： 1)环境友好-资源高效阶段，这一阶段土壤有效磷保持在20 mg·kg−1即可，在氮肥供应充分的条件下即可达到最高产量；2)环境低风险-资源低效阶段，该阶段土壤有效磷开始增加，但作物产量并不会显著升高； 3)环境有害-资源无效阶段，该阶段土壤有效磷超过了土壤的缓冲容量，土壤磷素开始向底层土壤迁移，同时作物产量不会升高。

2.3 秸秆阻控氮素淋失的微生物机制

与单施氮肥相比，秸秆还田增加土壤微生物生物量碳和氮，促进土壤氮素生物固定。图5 显示，相同施氮量下，在优化施肥的基础上添加秸秆处理中，表层(0～20 cm)土壤的效果更显著，0～40 cm 土层平均微生物生物量碳和氮含量较优化施肥处理分别平均增加14.1%～19.9%(P<0.05) 和 30.8%～49.6%(P<0.05)； 在常规施肥的基础上添加秸秆较单纯的常规施肥分别增加 12.7%～17.5%(P<0.05) 和23.9%～45.9%(P<0.05)。秸秆还田土壤硝化潜势较秸秆不还田处理降低幅度为15.4%～28.6%。

2.4 褐土区田间管理措施对氮磷淋失的影响

2.4.1 水肥优化措施

从表1 可知，在小麦-玉米轮作体系中，控水控肥可以维持小麦产量，与常规水肥处理相比，减量灌溉、减量施肥和优化水肥处理均未显著影响作物总产量。从表2 可知，在2016—2017年度和2017—2018年度减量灌溉、减量施肥和优化水肥处理的总氮周年淋失量均降低。氮素主要以有机氮形式淋失，其次为-N ，而-N 淋溶量很低。

在春玉米种植体系中，相较于常规肥处理，在优化施肥条件下(80%常规肥)并没有造成作物的减产，优化施肥处理3年，春玉米产量平均值较常规施氮处理增加6.1%(表1)。常规施肥和优化施肥添加秸秆均可降低0～300 cm 土壤剖面硝态氮残留，常规施肥效果更明显； 优化施肥较常规施肥可平均降低51%土壤硝态氮残留，且低于常规施肥+秸秆还田的硝态氮残留(图6)。

在华北褐土区露地菜地，优化施肥对蔬菜产量没有显著影响(表1)，但显著影响土体中氮素的淋溶量。常规施肥处理淋洗出剖面(80 cm)的总氮淋溶量占当季氮肥施用量的15.6%，优化施肥使氮淋溶量降低23.6%。与常规施肥处理相比，优化施肥+生物炭、优化施肥+减量灌溉处理总氮淋溶降低近50%，水氮协同调控对降低氮淋溶效果显著(表3)。相应地，水肥、秸秆优化管理措施也可降低200 cm 土层氮素的残留，但添加生物炭尽管可降低总氮的淋溶，但没有降低200 cm 土层中氮素的残留。

表2 2016—2018年度水肥管理对小麦-玉米轮作系统周年氮素淋失量的影响Table 2 Effects of water and nutrient management on annual nitrogen leaching loss of winter wheat and summer maize rotation in 2016-2018 kg·hm−2

在设施蔬菜种植中，各水氮处理间产量差异不显著(表1)，说明相较于常规水肥措施，减少20%施肥量和灌溉量、施用生物炭、减少20%施肥量和灌溉量同时施用生物炭不会对产量造成明显影响。从图7 可知，与常规施肥相比，降低氮肥投入可降低32%的土壤剖面硝态氮残留量； 添加生物炭可使硝态氮残留量减少49%，降低灌溉或优化施肥基础上添加生物炭可降低56%硝态氮的残留，并对残留硝态氮向下迁移起到一定的阻控效果。

2.4.2 生物炭措施

在小麦-玉米轮作体系中，与优化施肥处理相比，在优化水肥的基础上施用生物炭未能显著影响作物产量(表1)，但显著降低了总氮淋失量，也显著降低了-N 和有机氮的淋失量。华北褐土区露地菜地在减少氮肥用量20%基础上添加生物炭同样可减少氮淋失量，施用生物炭使全年总氮淋溶量比常规水肥处理降低43.0%(表3)。

2.4.3 秸秆措施

图6 显示，常规施肥基础上添加秸秆，有助于降低41%的0～200 cm 土壤剖面硝态氮残留； 即使在优化施肥基础上，添加秸秆也能降低18%的硝态氮残留(图6)。说明秸秆添加可减少土壤硝态氮残留量。此外，添加秸秆条件下，土层中硝态氮的分配比例也发生了变化，优化施肥条件添加秸秆处理0～100 cm 和100～200 cm 土层中硝态氮的分配比例1∶0.8，而单纯优化施肥处理则为1∶1.7 (图6)。这一结果表明，优化施肥条件下添加秸秆有助于阻控硝态氮的向下迁移。在灌溉条件下，露地蔬菜地上添加秸秆反而观测到硝态氮和磷素淋溶加剧的现象(表3)。这可能与土壤添加秸秆后，透水性增强有关。

3 讨论

3.1 褐土区土壤氮淋失阈值、安全施肥量及地下水氮素溯源

针对小麦-玉米种植体系的长期试验(大约 20年)，基于76 组冬小麦和73 组玉米试验资料，评价了中国北方地区适宜施肥量、土壤硝态氮的缓冲能力和环境安全施肥量的关系，发现我国北方地区适宜施肥量已经超过了该区域的环境安全施肥量[40]。值得注意的是小麦施氮量高于环境安全施氮量占总数据的64%，玉米农学施氮量高于环境安全施氮量占总数据的52%。因此，小麦-玉米轮作区目前氮肥推荐量存在硝酸盐淋失风险； 氮肥推荐不仅基于产量，也要考虑环境安全。在田间管理上，改善土壤理化性状以提高土壤矿质氮固持能力的措施都会有助于减缓或抑制硝态氮淋失。在不同种植体系中，土体中硝态氮的缓冲能力和安全施肥阈值也存在显著差异。在玉米种植体系中发现，施氮量超过 160 kg(N)·hm−2时，各土层硝态氮累积量出现跃增，而且施氮在160～200 kg(N)·hm−2内，产量、吸氮量、氮肥利用率达到或接近最高[41-42]。土壤理化性状和种植体系如何影响安全施肥量，其影响土壤硝态氮缓冲能力的机制值得进一步关注。

表3 不同管理措施下露地菜地氮肥在土壤剖面的残留和淋溶特征Table 3 Nitrogen residue and leaching under different management practices in soil profile of open vegetable field kg⋅hm−2⋅a−1

其次，汾渭河谷的褐土区也是大中城市、城镇的分布地。尽管本研究已经发现，地下水硝酸盐浓度与氮肥大量施用密切相关，但生活污水及动物源的硝态氮来源同样不容忽视。对关中盆地地下水硝酸盐含量超标的样点进行氯离子含量监测，发现关中盆地地下水硝酸盐超标地区硝酸盐与氯离子相关性达极显著水平。但各自来源的贡献需要进一步深入研究。例如，通过硝酸盐氮氧同位素组成识别硝酸盐来源和贡献率，从而为硝酸盐污染治理提供更准确、可靠的科学依据。

3.2 褐土区土壤磷淋失阈值、安全施肥量及其对环境影响的滞后效应

褐土区土壤属于石灰性土壤，pH 和碳酸盐含量较高。20 世纪80年代初期，褐土地区土壤有效磷(Olsen-P)含量是限制作物产量的重要因素[34]。因此，大量磷肥开始在该地区施用，经过30 多年持续磷肥投入，该区域的农田土壤磷素含量得到显著改善，部分区域土壤磷素显著积累。Olsen-P 含量由1980年极度亏缺(小于5 mg·kg−1)变为目前的富磷水平(大于30 mg·kg−1)，已有报道指出部分果园土壤Olsen-P含量超过100 mg·kg−1[28]。尽管石灰性土壤对磷素的固定性强、吸附能力大，但土壤中Olsen-P 含量升高同样会导致土壤吸附容量和吸附饱和度降低[35]；Khan 等[43]通过长期试验研究发现，小麦-玉米轮作体系长期偏施肥(PK)或过量施磷(MNPK)导致磷素淋洗到300 cm 或以下。对土壤全磷、Olsen-P 以及水溶性磷的动态变化进一步分析，发现当土壤全磷超过940 mg(P)·kg−1或者Olsen-P 超过37 mg·kg−1时，土壤中水溶磷显著增加。这意味着磷素淋失风险增加。此临界值低于英国洛桑试验站Broadbalk 长期定位试验结果(Olsen-P，60 mg·kg−1)[44]。

土壤残留磷对水质的污染具有滞后效应。欧美国家在20 世纪80年代开始治理土壤磷素积累对流域河流水质污染问题，但至今仍未达到预期的治理效果[46-48]。治理措施的滞后效应与土壤中残留磷素的缓慢而持久的释放有关[49]。目前农田土壤磷素积累在我国已引起广泛关注，但人们对其潜在的影响依然关注不够。

3.3 褐土区地下水硝态氮污染的综合评价

褐土区为我国传统的粮食生产基地，但20 世纪80年代以来，苹果(Malus pumila)、猕猴桃(Actinidia chinensis)等的种植日益增加，果园的单位面积氮肥的投入量是粮田的2～3 倍，多年种植的老果园中硝态氮的残留量显著高于农田。15年园龄果园土壤中100～200 cm 土层中-N 残留量达1400 kg·hm−2[26]。果园的施肥深度一般在距离地表50 cm 以内，但在果园中距离地表100 cm、200 cm 甚至400 cm 的土体中均发现了大量的硝态氮积累现象。这些报道表明，褐土区果园硝态氮淋溶的强度要大大高于相邻的农田土壤。自20 世纪80年代以来，褐土区果园面积逐年升高，2011年仅陕西省水果面积达1114.27 hm2[50]。大部分果园是由农田转化而来，果园和农田交错分布。本课题的研究主要集中于农田和蔬菜地中氮磷的淋溶和阻控，果园中氮磷的淋溶还有待进一步研究。

褐土区村庄生活污水对地下水中氮素的贡献也不容忽视。早在20 世纪70年代，科研工作者就发现，关中地区硝态氮含量较高的浅层地下水分布与居民点有极其密切的关系。在居民点内及其附近，由于生活垃圾及污水皆没有纳入城市垃圾和污水管理系统，垃圾和污水中大量的含氮物质在适宜条件下矿化而成硝态氮，之后随水下渗，进入浅层地下水。曹胜伟等[51]利用氮氧同位素技术研究发现，污水及粪便是造成南阳盆地地下水硝酸盐污染的主要原因，其平均贡献率为73%； 其次为化肥的贡献占23%。傅雪梅等[52]研究发现，河北张家口地区地下水硝态氮污染来源中粪肥及生活污水占45.37%，土壤氮来源为41.39%，降水和化肥来源占13.24%。

村庄、果园、农田镶嵌分布，不同来源的硝态氮随着水分入渗逐步向下迁移到浅层地下水[53]。判断地下环境中硝酸盐污染的来源、比例及其迁移转化过程对氮污染防治具有重要的借鉴意义。建议未来在褐土区进一步加强村庄、农田、果园等复合生态系统氮磷的淋溶研究，从而达到防治地下水污染的效果。

4 结论

经过40 余年的氮磷化肥投入，褐土区土壤氮磷的残留量显著增加，其中氮肥的残留以土体硝态氮的深层积累为特点，磷肥的残留以耕层显著积累为特征。高量施肥区域已经出现地下水氮素含量超标。与当地常规水肥投入量相比，化肥与灌溉减量、施用生物炭或秸秆都不会降低作物或蔬菜产量，但会降低氮磷淋失量； 其中化肥减量、降低灌溉可显著降低氮磷的淋失，其次是施用生物炭； 灌溉条件下，秸秆添加促进水分入渗会提高氮磷的淋失。硝态氮淋失得到阻控也与微生物碳氮的提高、土壤硝化势降低/反硝化势升高有关。此外，需要关注褐土区地下水硝酸盐污染的溯源、土壤磷素残留及其对河流污染的长期性问题、褐土区粮果生态系统中土壤氮磷淋溶的复合效应。