小麦，作为我国重要的粮食作物之一，是我国最重要的商品粮食和贮藏物质。长期定位试验具有时间的长期性、气候的重复性和数据的连续性等特点，既可揭示土壤肥料效益演变、肥效规律和土壤质量变化^[1-3]，又可弥补和完善许多短期试验所面临的科学问题^[4-5]。故研究长期定位下小麦不同耕作方式对土壤质量和土壤肥力累积情况的影响^[6]，能够科学地评价耕作对土壤各方面指标的长期作用，在生产上可提供决策性建议，对小麦作物产量的影响具有重要的现实意义^[7]。

国外长期定位肥料试验始于1843年Lawes爵士和Gilbert爵士建立的英国洛桑土壤肥料长期定位试验站，之后，德国、美国、丹麦和日本等国家相继建立了长期肥料试验站^[8-9]。长期肥料试验为世界农业发展、合理施用化肥、提高农作物产量和品质打下了坚实的理论和实践基础。我国长期定位试验起步较晚，虽不能与国外的经典长期定位肥料试验相比，但我国涉及的定位试验数量多、地域分布广，具有很强的针对性和目的性，在土壤碳、氮、磷、钾等养分含量方面开展了大量研究，在实践应用方面也积累了大量宝贵的经验^[10-12]，如长期施有机和无机肥能使土壤有机质、全氮以及硝态氮含量有较大幅度增加^[13-15]。土壤的综合肥力指标(integrated fertility index，IFI)在一定程度上可避免主观随意性的影响，能对大量信息进行处理，较好地综合反映土壤肥力高低^[16]，故IFI在不同区域各土壤类型和不同作物上的土壤综合肥力评价得到广泛应用，取得显著成果。然而，大量研究主要是进行不同肥料单施及配施的定位试验，涉及耕作方式、肥料及其互作的黄土高原沟壑区长期定位不同轮作系统施肥对土壤肥力的数值化综合评价报道较少。

因此，本研究以长武县为试验点，以小麦连作系统和小麦–豌豆轮作系统为研究对象，以土壤综合肥力指标为评价指标，借助长期野外定位和室内分析相结合的方法，系统分析不同小麦轮作系统下施加不同无机肥和有机肥对土壤养分的影响，综合评价在长期定位两种耕作系统不同施肥条件下土壤肥力的变化特征和肥效规律^[17]，初步探讨不同轮作系统对土壤养分和土壤肥力关系的影响，以期为改善土壤环境和提高土壤肥力提供数据支撑，为区域农业生产可持续发展提供理论依据。

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

试验区域地处陕西省咸阳市长武县，位于黄土高原的中南部，地理位置107°38′49″ ~ 107°58′02″E，34°59′09″ ~ 35°18′37″N，属于黄土高原沟壑区，塬面平坦宽阔，黄土堆积深厚，土壤类型为黏化黑垆土。试验于1984年开始，持续30余年，试验地海拔1 200 m，年均降水量578.5 mm，年均气温9.2 ℃，> 10 ℃积温3 019 ℃，无霜期172 d，属于典型的暖温带半湿润大陆性季风气候，是典型的雨养农业区。试验区域土壤肥力状况如表 1所示。

1.2 试验设计与分析方法

从长期定位试验中选取2大类，10个试验组，分别是：①轮作种植系统，分为小麦–豌豆轮作种植系统的NP(氮、磷肥配施)、NPM(氮、磷、有机肥配施)、P(单施磷肥)、CK(不施肥)；②小麦连作种植系统的NPM(氮、磷、有机肥配施)、NP(氮、磷肥配施)、P(单施磷肥)、N(单施氮肥)、M(单施有机肥)、CK(不施肥)开展研究工作，每个试验组3个重复。

定位试验30 a后，即2014年小麦收获期按小区随机在两个小麦耕作系统各处理采集耕层(0 ~ 20 cm)土壤样品，测定其有机质、全氮、全磷、碱解氮、速效磷、速效钾等含量。有机质采用重铬酸钾法测定，全氮采用凯氏定氮法测定，全磷采用H₂SO₄-HClO₄消煮–钼锑抗分光光度法测定，碱解氮采用碱解扩散法测定，速效磷采用Olsen法测定，速效钾采用火焰光度法测定，测定方法详见《土壤农化分析》^[18]。

1.3 数据处理

各项肥力指标测定数据及肥力综合评价结果采用Microsoft Excel 2007进行分析与处理，数据方差分析应用DPS 7.05软件进行方差分析和多重比较(α=0.05)，所有数据均以平均值±标准差表示。

1.4 土壤肥力综合指标值计算 1.4.1 土壤肥力指标因素的选取

在进行土壤肥力综合指数计算时，指标的选取直接影响到土壤肥力评价的真实性、合理性和科学性，故土壤肥力评价指标应全面、综合地反映土壤肥力的各方面^[19]。选取指标因素的原则是：①指标对作物生长发育和生产力影响大；②指标应以稳定性高的土壤养分含量为主；③指标间差异较大、相关性小^[20]。综合目前长期定位试验肥力评价研究^[21-22]，本文选取有机质，全氮、碱解氮、速效磷和全磷5项土壤基本养分指标计算土壤综合肥力。

1.4.2 隶属度值的确定

由于不同土壤肥力评价指标最适范围不同，量纲也存在差异，因此首先对原始数据进行标准化处理，使得处理后的数据区间处于[0.1, 1]。根据土壤养分特性和植物生长需求，参考已有研究结果可知^[23]，全氮、碱解氮、有机质、全磷、速效磷和速效钾等属性采用S型隶属函数，隶属度函数为：

$F(X) = \left\{ \begin{gathered} {\rm{ }}1.0{\rm{ \;\;\;\;\;\;\;\;\; \;\;\;\;\;\;\;\;\; }}X > {X_2} \\ 0.9(X - {X_1})/({X_2} - {X_1}) + 0.1{\rm{ \;\;\;\;\;\;\;\;\; }}{{\rm{X}}_{\rm{1}}}{\rm{ < X < }}{{\rm{X}}_{\rm{2}}} \\ {\rm{ }}0.1{\rm{ \;\;\;\;\;\;\;\;\; \;\;\;\;\;\;\;\;\; }}X < 0 \\ \end{gathered} \right.$

探讨长期施肥对不同小麦耕作系统土壤肥力状况的影响，观测各个土壤养分指标不同时期的动态变化，并对土壤养分定位观察，根据全国第二次土壤普查养分分级标准(表 2)确定隶属函数中转折点的取值X₁、X₂(表 3)。根据上述隶属度函数以及转折点的取值，计算各肥力指标的隶属度值(表 4)。隶属度值大小介于0.1 ~ 1.0，反映了其隶属的程度，最大值1.0表示土壤肥力最好，适宜作物生长；最低值0.1表示土壤肥力严重缺乏。土壤中不可能没有某种养分，为符合生产实际和消除各参数指标间的量纲差异，故将最小值定为0.1。

1.4.3 权重的确定

各项指标对土壤肥力的贡献是不同的，对各项指标应给以一定的权重。为避免人为的主观影响，应根据土壤肥力各要素之间的相关系数来确定权重系数^[24]。计算步骤为：①求单项肥力指标之间的相关系数(表 5)；②求某项肥力指标与其他肥力指标之间相关系数的平均值，并根据该平均值占所有肥力指标相关系数平均值总和的百分比(权重系数)，作为该单项肥力指标在表征土壤肥力中的贡献(表 6)，两项指标若呈负相关，计算相关系数取绝对值。

根据权重累加法，土壤的综合肥力指标(IFI)^[25-26]为：

${\rm{IFI}} = \sum\limits_i^n {Wi \times Ni} $

(1)

式中：Ni和Wi分别表示第i项隶属度值和权重系数。IFI值越大，表明土壤的综合肥力水平越高。

2 结果 2.1 小麦连作与小麦-豌豆轮作长期定位施肥对土壤养分的影响

在30 a长期定位不同施肥条件下，小麦–豌豆轮作系统土壤养分含量总体上高于小麦连作系统(表 7)。与连作系统相比较，小麦轮作系统未施肥土壤样地(CK)、磷肥添加样地(P)、氮磷配施样地(NP)和氮磷+有机肥配施样地(NPM)土壤有机质分别增加了1.1、1.13、1.25和1.22 g/kg。除NPM样地轮作系统全氮含量显著高于连作系统外，其他不同施肥样地两种耕作系统土壤全氮含量差异不显著。土壤全磷只在P处理条件下，轮作系统显著高于连作系统，增加了8%。土壤速效钾含量在NPM处理下连作系统显著高于轮作系统，其余施肥处理下两种系统差异不显著。土壤速效磷在不同施肥处理下，轮作系统均高于连作系统，以NPM增幅最明显。土壤碱解氮含量，在不同施肥条件下连作系统显著高于轮作系统。

小麦–豌豆轮作系统中，长期不同施肥对土壤养分含量变化有所差异。土壤碱解氮含量不同施肥处理表现为NPM和NP显著高于P和CK处理。土壤速效磷和全磷含量不同施肥处理之间差异显著，变化趋势表现为NPM > P > NP > CK。与CK相比较，NP和P处理下土壤速效钾含量增幅较小，而NPM处理显著增加了101%。土壤有机质和全氮含量在不同施肥处理下变化特征相似，具体表现为NPM > NP > P=CK；与CK相比较，NPM和NP有机质分别增加了73% 和27%，全氮分别增加了216% 和29%。

小麦连作系统中，NPM处理土壤养分均高于其他施肥处理与未施肥处理。土壤碱解氮含量具体表现为NPM和NP处理优于其他施肥处理和未施肥处理。土壤速效磷和土壤速效钾含量差异显著，在NPM和M处理下增量较高，速效磷含量分别增加了15.7倍和9.1倍，速效钾含量分别增加了2.1倍和1.9倍。试验结果表明，无论氮磷+有机肥配施，还是单施有机肥，都有很好的培肥土壤作用，说明有机肥在农业生产中有非常好的土壤培肥效果。土壤有机质和全氮在不同处理下变化特征相似，以NPM、NP和M处理增加最为明显，与CK相比较，土壤有机质含量分别增加79%、28% 和38%，全氮分别增加78%、29% 和34%。土壤全磷含量在NPM处理下最优，NP和P处理次之，N与M处理相当且均高于CK。

2.2 不同施肥处理和耕作方式的土壤肥力评价

由综合肥力指数(IFI)的变化特征可知，虽相同施肥条件下小麦–豌豆轮作系统和小麦连作系统土壤肥力差异不显著(除NPM处理外)，但轮作系统土壤肥力优于连作系统(图 1)。在小麦连作系统和小麦–豌豆轮作系统中，土壤肥力因施肥种类不同而有所差异，具体表现为：在小麦连作系统中，与CK处理相比，NPM、NP、P、N、M处理的IFI值分别提高186%、132%、117%、47% 和92%；小麦-豌豆轮作系统中，与CK处理相比，NPM、NP、P处理的IFI值分别提高197%、115% 和108%。

2.3 土壤养分指标与肥力综合指数相关性

土壤综合肥力指数与土壤基本理化指标的相关性分析发现(表 8)，将两种种植系统数据整合后，有机质、全磷、全氮和速效磷与IFI之间呈显著相关关系，碱解氮和速效钾与IFI不相关。在小麦连作系统中，速效磷、有机质、全氮和全磷与IFI呈显著正相关关系。在小麦–豌豆轮作系统中，碱解氮、全磷与IFI呈显著正相关关系。

3 讨论

长期定位施肥下小麦–豌豆轮作系统能够有效改善土壤肥力，提高土壤养分(如有机质、速效磷)^[27-28]，可能是因为豆科作物相对于其他作物而言，具有较强的固氮作用，能够促进土壤有机质的积累和养分的转化。土壤有机肥自身含大量腐殖质，在有机肥和豆科作物的双重作用下，土壤有机质可得到显著提高，但全氮、全磷含量增加不显著，碱解氮含量小麦连作高于小麦-豌豆轮作系统，可能是因为豆科固氮主要以有机态氮储存^[29-30]。试验结果(表 7)表明：①相同的施肥处理下，小麦–豌豆轮作系统，有机质、速效磷增幅高于小麦连作，碱解氮、全氮和全磷增幅低于小麦轮作。②小麦–豌豆轮作与小麦连作系统中，单施有机肥、氮磷配施和氮磷+有机肥配施可显著提高土壤全氮含量，是因为有机肥可以提高土壤氮贮量，改善土壤供肥能力，化肥配施土壤氮优于单施氮肥处理，全面配施提高了土壤速效氮的供应能力，从而增加土壤氮的含量^[31-32]。

从IFI的变化趋势(图 1)看，长期定位下进行相同施肥处理在小麦连作系统和小麦–豌豆轮作系统中土壤肥力表现不同。长期不施肥处理CK，由于没有养分的供应加上年年作物生长消耗，导致土壤肥力下降^[33-34]。氮磷配施和单施磷肥处理比较，IFI值差异性不显著，说明单方面施用氮磷肥或者磷肥在小麦连作与小麦-豌豆轮作系统中综合肥力不受影响。氮磷和有机肥配施处理，小麦–豌豆轮作系统IFI值显著高于小麦连作系统，说明氮磷和有机肥配施可以提高土壤肥力，究其原因：一方面氮磷和有机肥配施不仅直接增加了氮、磷和有机质等养分含量，而且显著增加有益微生物，改善作物根系环境，从而提高土壤肥力^[35]；另一方面在轮作系统下豆科作物与有机肥的双重作用提高了土壤肥力。

土壤综合肥力指数与土壤基本理化指标的相关性分析(表 8)表明：小麦两种种植系统数据综合与单一小麦连作系统中，差异显著性趋势大体一致，这可能是由于小麦连作样本较多，对IFI影响较大；有机质、全磷、全氮和速效磷含量与IFI差异显著正相关，说明其对黏化黑垆土的土壤肥力贡献较大，与有关学者研究紫色土结论一致，是表征土壤肥力的重要指标^[36]。小麦-豌豆轮作系统中，碱解氮与IFI显著相关的原因可能与固氮作物具有固氮能力有关^[37]。以上试验数据表明，有机质、全磷、全氮和速效磷含量对黏化黑垆土土壤肥力的贡献较大，应用相关性分析法计算的综合肥力指数能较好地反映土壤肥力状况。

4 结论

1) 小麦-豌豆轮作系统和小麦连作系统不同施肥处理比较，有机质、全氮、速效磷含量均有所增加，土壤全磷差异不显著，土壤速效钾和碱解氮含量均有所降低。小麦-豌豆轮作系统中，不同肥料配施与CK比较，碱解氮、速效磷、有机质、全氮和全磷含量均不同程度增加，其中氮磷肥和有机肥配施处理增幅明显。小麦连作系统中，不同肥料配施与CK相比较，土壤有机质、全氮、全磷、碱解氮、速效磷和速效钾含量增加。

2) 经过30 a长期不同施肥，不同处理间土壤IFI表现为：小麦-豌豆轮作+NPM > 小麦连作+NPM > 小麦-豌豆轮作+NP、小麦连作+NP > 小麦-豌豆轮作+P、小麦连作+P > 小麦连作+M > 小麦-豌豆轮作+不施肥 > 小麦连作+不施肥。长期氮磷肥和有机肥配施相较单施氮磷肥、有机肥的土壤综合肥力要高，小麦-豌豆轮作系统中施用氮磷肥和有机肥土壤综合肥力优于小麦轮作系统，对于区域农业实践具有指导意义。

3) 不同耕作系统对土壤养分和土壤肥力的关系有一定的影响。小麦连作系统有机质、全磷、全氮和速效磷与土壤综合肥力指数IFI之间呈显著相关，碱解氮和速效钾与IFI不相关。小麦-豌豆轮作系统中，碱解氮、全磷和IFI呈显著正相关。