黄土高原是中国水土流失最严重与生态环境最脆弱的地区，同时承担着重要的农业生产任务。这一地区曾因开垦农地使生态环境遭到严重破坏^[1]，虽然退耕还林还草工程的大规模实施极大改善了区域生态环境，但也导致耕地面积减少，出现可利用耕地紧张，粮食供需不平衡的局面^[2-3]。为解决局部人—粮矛盾等突出问题，延安市治沟造地项目从2010年开始试点，2013年被列为中央支持的国家土地整治重大工程项目^[4]。治沟造地工程针对黄土丘陵沟壑区特殊地貌，通过闸沟造地、打坝修渠、垫沟覆土等措施，旨在增加区域耕地面积，提高耕地质量，保证粮食安全、解决人地矛盾^[5]。截至2018年底，已完成治沟造地33 973.33 hm²，新增耕地6 860 hm^2[6]。

土壤养分是构成土壤肥力的物质基础，其时空分布特征对土壤环境质量和作物健康状况都有着重要的影响^[7]。土地整治措施例如机械碾压、土壤深翻、表土剥离以及回填，打破了土壤的原始状态，可能对土壤的有机质及速效养分含量产生一定的影响。研究发现，经过土地整治，土壤肥力在一段时间内呈现下降趋势，土壤氮、磷、钾等养分和有机质含量减少，微量元素降低^[8-11]。自开展治沟造地工程以来，我国学者对黄土高原地区沟道土壤质量和养分的变化进行了一些研究。王鹏等^[12]以13个治沟造地项目为研究对象，分析评价了治沟造地新增耕地的土壤养分特征，表明治沟造地形成的新增耕地质量等级低，由于利用时间短，造成有机质、速效钾和有效磷含量低；付威等^[13]研究表明合适的培肥可快速提高治沟造地新构造土地的土壤有机质，改善耕层土壤物理及化学性质，提升土地生产力；陈淑梅等^[14]从陕北不同沟道土地盐碱化的研究中发现，治沟造地的土层较薄，加之水库或淤地坝的阻挡作用，使得地下水位较浅，土壤易发生盐碱化。然而，目前在沟道土地整治后较长时间尺度上土壤养分的变化还未见深入研究。因此，本文以陕西省延安市宝塔区南泥湾镇九龙泉土地整治项目为例，对不同整治年限沟道不同位置的土壤样品进行采集，研究整治前以及整治后不同年限土壤有机质、氮、磷、钾指标的变化特征，为这一地区改进土地整治方法，采取合理耕种管理措施，提高养分资源利用率提供有效建议；并通过物元模型对治沟造地工程前后土壤养分进行综合评价，为治沟造地项目效益评价提供科学依据。

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

延安市宝塔区位于北纬36°11′—37°02′，东经109°14′—110°07′，南北长96 km，东西宽76 km。属半湿润半干早大陆性季风气候，冬季寒冷干燥，降水较少；夏季炎热多雨，降雨集中，多年平均降水量为562.1 mm，降水量年内分配不均，6—9月平均降水量占全年降水量的70%左右。全年平均气温9 ℃，平均无霜期179 d。宝塔区位于黄土高原中部，属于黄土丘陵沟壑区，土壤类型为黄绵土。区内地形破碎，沟壑纵横，冲沟下切强烈，地形坡度较大，水土流失严重。

宝塔区南泥湾镇治沟造地项目是延安市治沟造地整治示范工程和样板工程。本文选取九龙泉典型治沟造地项目区作为研究对象。九龙泉沟地处汾川河流域，位于宝塔区南部的梁峁丘陵沟壑区，属河谷川、台地貌，西南—东北走向，沟道长约9.5 km，河道海拔自北向南由1 170 m降至1 093 m，平均比降0.78%，河谷宽一般在250~500 m之间，面积为397.35 hm²。沟道土地多开垦为农田，种植前为旱地，农作物一年一熟，水资源比较匮乏，无灌溉条件，粮食单产水平偏低，春玉米产量约3 750~5 700 kg/hm²。治沟造地工程于2013年初开工建设，建设规模360.91 hm²，总投资4 411.65万元。通过治沟造地工程措施，完成土地平整，建立灌溉排水、道路桥涵、生态防护工程。该区主体工程于2014年初竣工，并于当年开始耕种，种植作物主要为玉米，每年于4月底至5月初播种，10月收获，产量提高到约7 500 kg/hm²。种植年限内，施用复合肥，肥料作为基肥在播种时同时施入，尽量保证撒施均匀，同时在拔节期追施氮肥，施肥量与施肥方式与整治前相同，玉米收获后秸秆部分粉碎还田，冬季休耕，其他耕作措施与当地传统耕作方式相同。

1.2 土壤样品的采集及测定

在九龙泉沟沟道，利用GPS定位，从沟头位置开始，顺沟道向沟口每隔一段距离(800~1 300 m)设置一个点，沟头、沟道中部、沟口共设置9个固定样点(表 1)。于2013年4月项目开工前、2015年7月、2016年7月及2019年6月采样，每次按照GPS信息定位固定样点，再以每个固定样点为中心，根据地形地貌，在10 m半径内，按照“S”形随机选择5个土壤样点，采用直径40 mm的土钻采集0—30 cm土层的土壤样品，每个土样1 000 g左右，重复3次，同时注意在采样时避开施肥点。将采集的样点土壤充分混合均匀后，得到整治前(0 a)、整治后2 a(2 a)、整治后3 a(3 a)和整治后6 a(6 a)4个时期土壤样品。样品在实验室经风干、研磨、过筛后进行指标测定。土壤有机质采用K₂CrO₇-H₂SO₄氧化法测定，全氮采用半微量凯氏定氮法测定，有效磷采用NaHCO₃浸提—钼锑抗比色法测定，速效钾采用CH₃COONH₄浸提—火焰光度法测定^[15]。

1.3 数据处理

研究采用经典统计分析、地统计分析相结合的方法，采用SPSS 18.0(SPSS Inc., Chicago, USA)进行统计分析，不同整治年限的差异显著性采用单因素方差分析(ONE-WAY ANOVA，LSD)检验。

1.4 物元模型评价方法

研究应用物元模型进行土壤养分评价。由于土壤养分是一个多指标体系，明确的评价标准可能会导致信息遗漏，而基于可拓理论的物元模型显示了土壤养分的过渡状态，建立起能够变动和分解的内部结构，并计算各个指标与整体的关联度，全面地反映土壤养分的内部结构与各个指标之间的关系，由此综合判断区域土壤养分的基本情况^[16-17]。物元模型分析首先对每个评价指标进行分级，定义区间和取值范围，确定经典域和节域，然后通过单个指标的关联函数计算单个指标的土壤养分等级，再通过模型集成得到多指标的综合土壤养分水平。

2 结果与分析 2.1 土壤养分指标统计分析

就土壤养分各指标均值而言，除了速效钾以外，研究区沟道内土壤养分含量普遍较低(表 2)，根据全国第二次土壤普查土壤养分分级标准^[18]，有机质属于4—5级水平，全氮属于2—6级水平，有效磷属于3—6级水平，速效钾相对较高，属于1—3级水平。土地整治前，沟道中土壤养分变异程度较大，变异系数介于0.30~0.77之间，属于中等变异程度^[19]，其中有效磷变异程度最大，离散程度高，其次是全氮和有机质，速效钾变异程度最小。治沟造地工程实施后，沟道中的各养分变异系数整体上有变小的趋势，6 a时变异系数介于0.23~0.54之间，土壤养分各指标含量差异变小，在沟道中的分布趋于均一。

2.2 土壤养分时空变异特征 2.2.1 土壤有机质

由图 1可以看出，整治前，九龙泉沟沟头有机质含量为12.02 g/kg，沟道中部稍高于沟头，沟口有机质含量最高，为19.27 g/kg。工程实施后，沟道中的土壤有机质含量呈现出随整治年限增加而先减少后增加的趋势，但整治后6 a的土壤有机质含量平均值仍然低于整治前。沟头整治后2 a有机质含量下降为9.44 g/kg，随着整治年限的变长有机质含量逐渐增加，6 a时达到11.75 g/kg，基本恢复到整治前水平，随整治年限增长的变化并不显著。沟道中部的土壤有机质含量在2 a显著下降，下降了30%左右，3 a继续下降，6 a时稍有回升，达到11.09 g/kg，但从2 a到6 a的变化并不显著，仍然显著低于整治前。沟口土壤有机质含量的变化与沟道中部基本一致，整治后2 a显著下降，2 a以后逐渐增加，整治后6 a仍比整治前水平低15.4%。从沟道不同位置整治年限的差异显著性分析可知，治沟造地工程对沟头土壤有机质含量的影响较小，有机质含量变化不明显，对沟道中部的影响更为显著。整治后，由于表土剥离的扰动或沟坡侧面贫瘠土壤的覆盖，土壤表层有机质含量低于整治前；之后，土壤生态系统经过2~3 a的调整，通过施肥等管理措施以及作物的自肥作用，加深了土壤的腐殖质层，使表层有机质含量增加，当整治年限增至6 a时，土壤系统进一步趋于稳定，植物的生长和根系下扎也提高了下层土壤有机质的来源，促进了土壤有机质的形成，其含量进一步升高。

2.2.2 土壤全氮

整治前，全氮含量表现为：沟头＜沟道中部＜沟口，治沟造地工程后，沟道中的土壤全氮含量呈现出随整治年限先增加后减少的趋势，整治后所有年限的土壤全氮含量均高于未整治土壤。沟头整治后全氮含量逐渐增加，3 a时达到最高0.78 g/kg，随后稍有减少，3 a及6 a的含量均显著高于整治前，是整治前的2倍以上。沟道中部的土壤全氮含量在2 a显著增加，随后逐渐减少，6 a时回落到与整治前相当水平。沟口土壤全氮含量的变化趋势与沟道中部较为一致，2 a时显著增加，随后逐渐降低，其增加的比例比沟道中部更大(图 1)。

2.2.3 土壤有效磷

同全氮含量变化趋势一致，随着整治年限变化，研究区土壤有效磷含量整体上呈现出随治理年限先增加后减少的趋势，整治后土壤有效磷含量均高于未整治土壤。沟头未整治土壤有效磷含量较低，为5.01 mg/kg，整治后2 a有效磷含量逐渐上升，到3 a时达到最高值14.62 mg/kg，为整治前的近3倍，6 a时回落至6.29 mg/kg，与整治前差异不显著。沟道中部和沟口的土壤有效磷含量随整治年限变化趋势与沟头一致，其中沟道中部土壤有效磷含量变化的幅度最小(图 1)。

2.2.4 土壤速效钾

整治前土壤速效钾含量表现为：沟头＞沟道中部＞沟口，沟头最高为255.04 mg/kg。整治年限对速效钾的影响与氮、磷不同，从沟道中土壤速效钾含量的平均值可看出，整治后土壤中的速效钾含量基本上呈现下降的趋势。沟头土壤速效钾含量呈现了下降—回升—再下降的过程，整治后2 a下降较为显著，所有年限均低于整治前。沟道中部的土壤速效钾含量从整治前的184.67 mg/kg持续降低到整治后6 a的146.93 mg/kg，整个过程变化并不显著。沟口的土壤速效钾含量的变化与沟头相似，2 a显著下降，3 a稍有提高，6 a持续了下降的趋势，整治后2 a，6 a土壤速效钾含量均显著低于整治前(图 1)。速效钾含量与全氮、有效磷的变化规律不同，未见显著增加，这可能与流域内速效钾含量较高，农业施肥等措施影响不显著有关。

2.3 土壤养分综合评价

本文选取有机质、全氮、有效磷、速效钾4个指标作为评价指标，通过层次分析法确定权重^[20]，并运用物元模型，对治沟造地项目土壤养分等级进行评价研究。

2.3.1 指标权重的确定

采用AHP法对各指标进行赋权(表 3)。其中λ_max=3.98，计算出一致性指标I_c=0.002 6，查表得出平均随机性指标IR=0.9，则随机一致性比率R_c=I_c/IR=0.003＜0.1，表示结果通过一致性检验，具有满意一致性。

2.3.2 经典域、节域和待判物元的确定

依据可拓集合概念，将土壤养分分为丰富、较丰富、适量、贫乏、极贫乏5个等级进行分等评价，分别用罗马数字Ⅰ-Ⅴ表示。首先，拟定N₁，N₂，N₃，N₄，N₅这5个连续的数据集，分别对应5个评价等级。

采用全国第二次土壤普查土壤养分的分级标准^[18]，建立土壤养分评价的经典域物元R₁，R₂，R₃，R₄，R₅和节域物元R_P。由于速效钾数据中有部分超出经典域范围，故对应c₄的量值范围进行相应的变化。

$ {\mathit{\boldsymbol{R}}_{\bf{1}}} = \left[ {\begin{array}{*{20}{c}} {{N_1}}&{{c_1}}&{\langle 30, 40\rangle }\\ {}&{{c_2}}&{\langle 1.5, 2\rangle }\\ {}&{{c_3}}&{\langle 20, 40\rangle }\\ {}&{{c_4}}&{\langle 150, 260\rangle } \end{array}} \right] $

$ {\mathit{\boldsymbol{R}}_2} = \left[ {\begin{array}{*{20}{c}} {{N_2}}&{{c_1}}&{\langle 20, 30\rangle }\\ {}&{{c_2}}&{\langle 1, 1.5\rangle }\\ {}&{{c_3}}&{\langle 10, 20\rangle }\\ {}&{{c_4}}&{\langle 100, 150\rangle } \end{array}} \right] $

$ {\mathit{\boldsymbol{R}}_3} = \left[ {\begin{array}{*{20}{c}} {{N_3}}&{{c_1}}&{\langle 10, 20\rangle }\\ {}&{{c_2}}&{\langle 0.75, 1\rangle }\\ {}&{{c_3}}&{\langle 5, 10\rangle }\\ {}&{{c_4}}&{\langle 50, 100\rangle } \end{array}} \right] $

$ {\mathit{\boldsymbol{R}}_4} = \left[ {\begin{array}{*{20}{c}} {{N_4}}&{{c_1}}&{\langle 6, 10\rangle }\\ {}&{{c_2}}&{\langle 0.5, 0.75\rangle }\\ {}&{{c_3}}&{\langle 3, 5\rangle }\\ {}&{{c_4}}&{\langle 30, 50\rangle } \end{array}} \right] $

$ {\mathit{\boldsymbol{R}}_5} = \left[ {\begin{array}{*{20}{c}} {{N_5}}&{{c_1}}&{\langle 0, 6\rangle }\\ {}&{{c_2}}&{\langle 0, 0.5\rangle }\\ {}&{{c_3}}&{\langle 0, 3\rangle }\\ {}&{{c_4}}&{\langle 0, 30\rangle } \end{array}} \right] $

$ {\mathit{\boldsymbol{R}}_\mathit{\boldsymbol{p}}} = \left[ {\begin{array}{*{20}{c}} {{N_p}}&{{c_1}}&{\langle 0, 40\rangle }\\ {}&{{c_2}}&{\langle 0, 2\rangle }\\ {}&{{c_3}}&{\langle 0, 40\rangle }\\ {}&{{c_4}}&{\langle 0, 260\rangle } \end{array}} \right] $

确定待判物元。以整治前沟头为例，该点平均数据为：有机质12.02 g/kg，全氮0.34 g/kg，有效磷5.01 mg/kg，速效钾255.04 mg/kg，以此数据确定该点S₁的待判物元。

$ {\mathit{\boldsymbol{R}}_{\mathit{\boldsymbol{x}}{\bf{1}}}} = \left[ {\begin{array}{*{20}{c}} {{N_{s1}}}&{{c_1}}&{12.02}\\ {}&{{c_2}}&{0.34}\\ {}&{{c_3}}&{5.01}\\ {}&{{c_4}}&{255.04} \end{array}} \right] $

2.3.3 计算关联度，并确定质量等级

将待判物元带入关联度公式，以单指标c₁(有机质)为例(v₁=12.02)，可得到该指标对应的各评价等级关联度分别为K₁(v₁)=-0.599 2，K₂(v₁)=-0.398 8，K₃(v₁)=0.202 4，K₄(v₁)=-0.144 1，K₅(v₁)=-0.333 8，依据判断标准，判定该指标的土壤养分级别属于N₃，即适量级别。

用同样方法可以求出S₁其他各项指标对应各等级的关联度及其养分级别。根据各项指标所对应的各等级关联度及权重，计算得出所有指标对应各评价等级的综合关联度，并依据判断标准，判定各点土壤养分级别(表 4)。根据表 4，治沟造地工程前，九龙泉沟沟道中土壤养分等级均属于适量等级。工程实施后2 a，沟头的土壤养分评价等级下降到贫乏级别。但通过3 a的耕种，即恢复到适量水平。沟道中部的土壤养分等级一直维持在适量水平。沟口的土壤养分评价等级在3 a时，超过整治前水平，达到较丰富等级，6 a时再次恢复到适量水平。

3 讨论 3.1 土壤养分的时间变化

目前九龙泉沟道内，除钾素外，土壤磷素、氮素和有机质均处于缺乏状态。杨静涵、郝丽婷等人在这一地区土壤养分的研究中，均表明了相似结果^[21-22]。研究发现，经过整治的土地，由于表土剥离及生土覆盖等原因，土壤养分呈现下降趋势，土壤速效氮、速效磷、速效钾和土壤有机质显著减少。然而，经过一段时间的恢复(3~5 a)，同时加以施用化肥等措施，土壤养分各项指标均有了明显的提高，甚至高于整治前^[23-25]。本研究中，全氮、有效磷在整治后第2, 3 a就呈现出高于整治前的增长趋势，说明沟道土地整治通过合理的耕种管理措施，能较快抵消因土地整治扰动对土壤养分带来的负面影响，产生正面的影响。因为黄土高原地区土壤条件相对贫瘠，各种养分含量较低，虽然土地整治扰动对土壤养分带来了负面影响，但治沟造地工程后通过更加便利的耕作管理措施，如机器翻耕，增加了土壤孔隙，保证了耕层土壤均匀，确保降水可顺利入渗补给，且有利于养分的积累，因此施肥措施使氮磷等土壤养分含量增长较快，并在耕种一定时间后，随着作物吸收等因素，趋于稳定状态。土壤有机质含量呈现出先减少后逐渐增加的趋势。因为有机质的形成积累是个缓慢的过程，并不能靠化学施肥迅速提升，因此直到整治后6 a，也未恢复到整治前水平。相比较速效养分，其增长具有长期性和滞后性。因此在后期耕作管理中，应采取施用有机肥、秸秆还田、种植绿肥等措施，提升土壤有机质^[26]，通过优化施肥方式来保证土壤中养分的平衡，促进农业的可持续发展。速效钾的含量相对较高，可能与黄绵土中钾素含量相对较高有关，因此施肥等措施对于钾素含量的影响并不显著。

3.2 土壤养分的空间变化

因水流的运移和沉积作用，沟口通常是沟道养分的聚集区，因此沟口土壤养分含量普遍高于沟头和沟道中部，沟道内土壤养分变异程度较大，属于中等变异程度。治沟造地工程实施后，由于土地整治措施以及后期耕种管理等因素，沟道内土壤养分变异系数随整治时间有变小的趋势，土壤养分分布趋于均一，这与王瑷玲、王宏庭等的研究结果一致^[27-28]。随着整治年限继续增加，沟道不同位置的养分差异有可能因地形、水土流失等原因进一步变大。因此，可根据沟道不同位置土壤养分空间变异特征，调整沟道作物的空间布局，实施不同程度和方法的管理措施，提高养分资源的利用率，增加粮食产量。

3.3 沟道造地优势

由于退耕还林使黄土高原耕地面积减少，引起人—地矛盾问题，已经有退耕土地进行复垦现象，但退耕地复垦使土壤理化性状迅速恶化，切不可取^[29]。沟道是流域中水养的富集地，通过适当的土地整治后十分适宜农业耕作，蕴藏着丰富的潜在耕地资源。有研究表明，黄土丘陵沟壑区沟道耕地粮食产量远远大于坡耕地，有“宁种一亩沟，不种十亩坡”的说法^[30]。本研究表明，治沟造地工程后3 a土壤养分评价等级就恢复到与整治前相当水平，与黄土区坡改梯工程相比，土壤养分的恢复时间更短^[31]。同时，土壤养分含量在适量水平，整治后土壤养分水平还有继续提高的潜力。因此相比较坡面，黄土丘陵沟壑区沟道更适合建造耕地，治沟造地工程能够有效缓解这一地区粮食安全与人地矛盾问题，巩固退耕还林成果。

4 结论

(1) 除速效钾外，研究区土壤养分含量普遍较低。沟道中土壤养分变异程度较大，变异系数介于0.30~0.77之间，属中等变异程度，其中有效磷最大，其次是全氮、有机质，速效钾变异程度最小。治沟造地工程后，沟道中土壤养分变异系数有变小的趋势，土壤养分分布朝均一化方向发展。

(2) 沟道土地整治后6 a时间内，随整治年限增加，土壤有机质含量呈现出先减少后增加的趋势，6 a时含量低于整治前；土壤全氮和有效磷含量呈现先增加后减少的趋势，6 a时高于整治前；土壤速效钾含量持续降低。合理的耕种管理措施能较快抵消因土地整治扰动对土壤养分带来的负面影响，产生正面影响。后期应注重土壤有机质含量水平提升。

(3) 沟道土地整治后3 a，土壤养分评价等级恢复到整治前水平，且有继续提高的潜力，可根据沟道不同位置土壤养分变异特征，调整沟道作物的空间布局，实施不同程度和方法的管理措施，提高养分资源的利用率，增加粮食产量。