土壤肥力是农业生产的基础^[1]，土壤肥力的高低直接影响着作物生长，影响着农业生产的结构^[2]。土壤养分是土壤肥力的核心部分，土壤养分的演变分析有助于了解区域土壤肥力变化情况和变化规律，对稳定粮食产量具有重要意义^[3]。潮土地区地势平坦、土层深厚、生产性状良好、适种性广^[4]。潮土约占全国耕地面积的15.9%，是我国小麦和玉米等粮食作物的主产区^[5]。潮土区土壤养分含量的变化在一些地区和时间段的研究已有一些报道，如钦绳武等^[6]分析了河南省封丘县潮土区1988—1992年来土壤肥力的变化特征，发现土壤有机质、全氮、全磷、全钾、速效氮、速效磷含量均有所上升，速效钾含量有所下降；张金涛等^[7]调查了全国潮土区1987—2006年的土壤肥力状况，发现土壤有机质和全氮含量有所增加，土壤pH值降低0.28个单位；杨玉建等^[8]在山东省禹城市的土壤有机质含量研究结果表明，禹城市2003年土壤有机质含量比1980年增加了8.68 g/kg，年均提高0.38 g/kg。虽然这些研究为潮土肥力的演变提供了有价值的参考，但研究地点分散，且时间较短，不足以为国家层面的决策提供依据。所以本文在分析29年来潮土养分现状和定量演变趋势的基础上，进一步分析了土壤肥力演变的影响因素及其演变过程的贡献因子，为潮土区耕地质量管理和肥料施用提供科学依据^[9-12]。

国家级潮土长期监测点主要分布在北京 (8个)、天津 (8个)、河北 (5个)、山西 (2个)、内蒙 (2个)、江苏 (6个)、安徽 (3个)、山东 (2个)、河南 (11个)、湖北 (5个)、湖南 (2个)、新疆 (3个) 等12个省份。监测点中大部分建于1988年，后期在1997年和2004年分别新增了一些监测点。至2016年，监测时间最长的点已经有29年的历史。

潮土监测区种植的作物主要为小麦和玉米，种植制度为一年两熟制，主要以小麦−玉米、冬小麦−夏玉米、春小麦−夏玉米等轮作方式为主。监测区采用机械化耕种为主。监测点耕层 (0—20 cm) 土壤初始性质见表1。

表1(Table 1)

表1 监测点起止年限及初始耕层 (0—20 cm) 土壤基本性质 Table 1 Periods of monitoring sites and the initial basic properties of fluvo-aquic soil

表1 监测点起止年限及初始耕层 (0—20 cm) 土壤基本性质 Table 1 Periods of monitoring sites and the initial basic properties of fluvo-aquic soil 序号 No. 省 (市) Province and city 地点 Site 监测年限 Monitoring years pH 有机质 SOM (g/kg) 全氮 Total N (g/kg) 有效磷 Available P (mg/kg) 速效钾 Available K (mg/kg) 缓效钾 Slowly available K (mg/kg) 1 北京市 Beijing 顺义区 Shunyi 1998—2016 7.47 16.3 1.11 88.1 117 771 2 北京市 Beijing 大兴区 Daxing 1998—2016 8.30 18.4 1.02 545 149 678 3 北京市 Beijing 大兴区 Daxing 1998—2016 7.67 16.8 1.12 142 192 690 4 北京市 Beijing 大兴区 Daxing 2004—2016 8.08 14.0 0.85 73.3 103 660 5 北京市 Beijing 房山区 Fangshan 2004—2016 8.02 26.5 1.28 67.5 299 914 6 北京市 Beijing 房山区 Fangshan 2004—2016 7.74 14.1 0.91 10.3 107 660 7 北京市 Beijing 房山区 Fangshan 2004—2016 7.96 13.2 0.83 24.7 111 706 8 北京市 Beijing 通州区 Tongzhou 2004—2016 8.11 19.4 1.15 108 152 993 9 天津市 Tianjin 宝坻区 Baodi 2004—2016 8.18 19.5 0.94 42.2 202 759 10 天津市 Tianjin 宝坻区 Baodi 2004—2016 8.19 17.8 1.00 31.4 278 1296 11 天津市 Tianjin 汉沽 Hangu 2004—2016 8.05 26.0 1.56 162 506 1795 12 天津市 Tianjin 大港 Dagang 2004—2016 8.27 16.5 0.99 17.0 247 1261 13 天津市 Tianjin 静海区 Jinghai 1988—2016 8.16 17.4 1.05 16.3 206 925 14 天津市 Tianjin 静海区 Jinghai 1988—2016 8.12 31.6 1.75 84.3 252 967 15 天津市 Tianjin 宝坻区 Baodi 1988—2016 8.08 19.2 1.19 38.1 157 836 16 河北省 Hebei 献县 Xian county 1988—2016 8.14 14.4 0.88 13.4 157 1055 17 河北省 Hebei 盐山县 Yanshan 1988—2016 8.03 11.6 0.76 9.87 117 951 18 河北省 Hebei 盐山县 Yanshan 1998—2016 8.10 16.6 1.03 13.8 110 789 19 河北省 Hebei 深州市 Shenzhou 1998—2016 8.11 14.2 0.88 18.7 56.0 615 20 河北省 Hebei 冀州市 Jizhou 1998—2016 8.30 15.9 1.05 17.0 1156 865 21 河北省 Hebei 冀州市 Jizhou 1988—2016 8.35 15.7 0.85 21.3 102 718 22 山西省 Shanxi 盐湖区 Yanhu 1998—2016 8.62 12.8 0.88 6.33 121 944 23 山西省 Shanxi 榆次区 Yuci 1988—2016 8.26 16.4 1.05 15.3 156 922 24 内蒙古 Inner Mongolia 鄂尔多斯 Erdos 1998—2016 8.44 18.0 0.92 16.6 138 970 25 内蒙古 Inner Mongolia 鄂尔多斯 Erdos 1998—2016 8.43 8.49 0.44 14.2 63.7 450 26 江苏省 Jiangsu 姜堰区 Jiangyan 2004—2016 7.35 19.9 1.21 15.49 68.3 494 27 江苏省 Jiangsu 通州区 Tongzhou 1998—2016 7.84 20.2 1.21 9.71 106 8116 28 江苏省 Jiangsu 大丰市 Dafeng 1998—2016 8.11 17.3 1.13 13.0 138 1036 29 江苏省 Jiangsu 宿豫区 Suyu 1998—2016 7.82 21.8 1.26 21.4 81.8 576 30 江苏省 Jiangsu 铜山区 Tongshan 1998—2016 8.04 21.8 1.30 18.7 1401 919 31 安徽省 Anhui 濉溪县 Suixi 2004—2016 6.78 18.7 1.25 10.8 225 759 32 安徽省 Anhui 淮北市 Huaibei 1998—2016 7.75 14.2 1.00 13.7 144 526 33 安徽省 Anhui 颍泉区 Yinquan 1998—2016 8.08 12.2 0.85 7.45 99.8 633 34 山东省 Shandong 牡丹区 Mudan 1998—2016 8.08 13.7 1.05 22.6 184 1016 35 山东省 Shandong 牡丹区 Mudan 1998—2016 8.04 12.1 0.93 18.7 116 912 36 河南省 Henan 民权县 Mingquan 2004—2016 8.09 8.26 0.54 24.3 64.9 371 37 河南省 Henan 民权县 Minquan 2004—2016 8.13 12.9 0.86 17.6 79.3 675 38 河南省 Henan 民权县 Minquan 2004—2016 8.16 12.0 0.78 8.68 68.3 657 39 河南省 Henan 淮阳县 Huaiyang 2004—2016 7.92 16.5 1.14 10.3 1667 852 40 河南省 Henan 唐河县 Tanghe 2004—2016 6.35 14.6 0.97 30.0 130 673 41 河南省 Henan 清丰县 Qingfeng 2004—2016 8.19 14.6 0.93 11.0 86.7 688 42 河南省 Henan 滑县 Hua county 2004—2016 7.93 14.7 0.95 15.6 133 711 43 河南省 Henan 兰考县 Lankao 2004—2016 8.45 9.00 0.68 9.73 73.7 539 44 河南省 Henan 梁园区 Liangyuan 2004—2016 8.04 14.4 0.93 13.4 149 853 45 河南省 Henan 获嘉县 Huojia 1998—2016 8.21 14.8 0.90 10.7 124 944 46 河南省 Henan 尉氏县 Weishi 1998—2016 8.06 16.2 0.93 14.0 60.1 539 47 湖北省 Hubei 嘉鱼县 Jiayu 2004—2016 7.20 24.2 1.44 23.5 131 496 48 湖北省 Hubei 松滋市 Songzi 2004—2016 7.16 19.8 1.23 7.51 180 985 49 湖北省 Hubei 华容区 Huarong 2004—2016 7.49 14.1 1.04 11.7 114 565 50 湖北省 Hubei 天门市 Tianmen 1998—2016 7.28 22.1 1.41 14.0 139 701 51 湖北省 Hubei 潜江市 Qianjiang 1998—2016 7.31 19.5 1.27 14.1 87.6 479 52 湖南省 Hunan 岳阳县 Yueyang 2004—2016 5.48 27.8 1.57 14.1 173 499 53 湖南省 Hunan 鼎城区 Dingcheng 1998—2016 4.92 20.4 1.61 45.4 250 280 54 宁夏 Ningxia 灵武市 Lingwu 2004—2016 8.22 17.9 1.15 81.8 151 868 55 新疆 Xinjiang 阜康市 Fukang 1998—2016 7.74 23.7 1.17 15.9 247 1043 56 新疆 Xinjiang 阜康市 Fukang 1998—2016 7.57 26.2 1.27 18.4 257 1039 57 新疆 Xinjiang 阜康市 Fukang 1998—2016 7.45 21.5 1.12 14.8 236 1036

每个监测点设置对照 (不施肥) 和常规施肥 (农民习惯施肥) 两个处理，依照当地农民习惯进行水肥管理等农事活动，并定位记载施肥量、肥料种类、作物产量以及管理措施等信息。每季作物收获后，采集耕层 (0—20 cm) 土壤样品，测定土壤有机质、全氮、有效磷、速效钾及pH值等土壤肥力指标。土壤养分测定采用常规分析方法^[13]。

试验数据用Excel 2016整理，运用Sigmaplot 12.5和SPSS 17.0 进行相关性分析及显著性检验。为避免个别年份点位差异对土壤养分的演变规律造成影响，因此，按照建点的时间将57个监测点划分为监测初期 (1988—1997年)、监测中期 (1998—2003 年) 和监测后期 (2004—2016年) 3个阶段。本研究中所涉及的土壤养分演变趋势的数据 (盒形图) 均以每组数据的中值表示^[14]，采用SPSS17.0^[15]的方法进行主成分分析。

监测初期 (1988—1997年)、中期 (1998—2003年) 和近期 (2004—2016年) 土壤有机质、全氮、有效磷、速效钾含量、pH值的含量变化，以及三个时期平均值的比较见图1。

图1(Fig. 1)

图1 潮土土壤肥力要素演变趋势 Fig. 1 Changes of fertility factors in fluvo-aquic soil at different monitoring periods [注（Note）：盒式图下边缘线和上边缘线分别代表全部数据的 5% 和 95%，上下实心点为异常值；矩形盒上、下边缘分别代表上四分位数和下四分位数，分别代表全部数据的 75% 和 25%；盒中实线代表中值，虚线代表平均值；矩形盒上的不同字母表示不同监测时期在 5% 水平上差异显著 The lower and upper edge lines represent the 5% and 95% of the data, the solid points represent the vertical outliers. The lower and upper quartiles of the boxplots represent 25% and 75% of data, and the solid lines represent the median values and the dash lines represent the average values. Different letters above the boxplot indicate significant differences among the three mentoring periods at 0.05 level.]

潮土区土壤有机质含量总体呈上升趋势，土壤有机质含量范围在3.70～37.8 g/kg。监测初期和监测后期土壤有机质含量呈上升趋势，监测中期略有下降。监测初期土壤有机质平均含量为11.5 g/kg，监测中期为16.9 g/kg，监测后期为17.2 g/kg。监测初期到监测中期呈显著上升趋势 (P < 0.05)，提升了47.1%，从监测中期到监测后期，虽然也呈升高趋势，但未达到显著水平。

土壤全氮含量变化与有机质变化相似，整体呈现上升的趋势。监测前期、监测中期和监测后期土壤全氮含量分别为0.75、1.03和1.06 g/kg，从监测初期到监测中期土壤全氮平均含量呈显著上升趋势 (P < 0.05)，与监测初期相比显著提升37.8%。而从监测中期到监测后期增加未达到显著水平。

监测区土壤有效磷含量范围为1.99～185.80 mg/kg。三个监测阶段土壤有效磷含量均呈上升趋势，监测初期土壤有效磷含量为8.67 mg/kg，监测中期虽略有升高，但增加不显著。从监测中期到监测后期土壤有效磷继续增加，虽然监测后期与中期有效磷含量差异未达到显著水平，但比监测初期有了显著增加 (P < 0.05)，提高了226%。

土壤速效钾含量监测前期和监测后期各试验点呈上升趋势，在监测中期各监测点变异较大。监测初期到监测中期土壤速效钾含量平均值略有升高，但未达到显著水平。从监测中期到监测后期呈显著上升趋势 (P < 0.05)，监测后期土壤速效钾的平均含量比监测中期上升了30.4%。

土壤pH水平在4.00～9.40之间，各监测期间都有下降的趋势，监测初期与中期接近，从监测中期 (8.14) 到监测后期 (7.78)，土壤pH呈显著下降趋势 (P < 0.05)。

常规施肥下，小麦与玉米的产量均高于不施肥处理，而且变化趋势与不施肥处理也不同 (图2)，与监测初期相比，潮土常规施肥条件下小麦和玉米产量呈显著增加趋势，无肥区小麦产量平均值为1773 kg/hm²，常规施肥区的平均产量为6008 kg/hm²。小麦常规施肥区产量较无肥区增加238%。无肥区玉米平均产量为2393 kg/hm²，常规施肥区平均产量为7281 kg/hm²。常规施肥区玉米平均产量较无肥区增产204%，潮土区施肥措施具有显著的增产作用。

图2(Fig. 2)

图2 常规施肥与不施肥处理的小麦和玉米产量变化 Fig. 2 Yield changes of wheat and corn under conventional fertilization and no fertilization [注（Note）：矩形盒中下边缘线和上边缘线分别代表全部数据的 5% 和 95%，上下实心点为异常值。矩形盒上、下边缘分别代表上四分位数和下四分位数，分别代表全部数据的 75% 和 25%，实线代表中值，虚线代表平均值；矩形盒上的不同字母表示不同监测时期在 5% 水平上差异显著 The lower and upper edge lines represent the 5% and 95% of the data, the solid points represent the vertical outliers. The lower and upper quartiles of the boxplots represent 25% and 75% of data, the solid lines represent the median values and the dash lines represent the average values; Different letters above the boxplot indicate significant differences among the three mentoring periods at 0.05 level.]

分别统计三个监测阶段施肥条件下的作物产量 (图3)，小麦产量监测初期 (1988—1997) 平均为2902 kg/hm²，监测中期 (1998—2003) 和监测后期 (2004—2016) 小麦产量平均值与监测初期相比显著上升 (P < 0.05)，分别上升了87.6%和114%。而从监测中期到监测后期呈上升趋势，但未达到显著水平。玉米产量的变化与小麦一样，从监测初期到监测中期，上升趋势明显 ( P < 0.05)，显著提升了111%。从监测中期到监测后期增加不显著。

图3(Fig. 3)

图3 不施肥处理与常规施肥处理下小麦和玉米产量的阶段性变化 Fig. 3 Periodical change of yields of wheat and corn under non-fertilized treatments and conventional fertilized treatments [注（Note）：盒式图下边缘线和上边缘线分别代表全部数据的 5% 和 95%，上下实心点为异常值；矩形盒上、下边缘分别代表上四分位数和下四分位数，分别代表全部数据的 75% 和 25%；盒中实线代表中值，虚线代表平均值；矩形盒上的不同字母表示不同监测时期在 5% 水平上差异显著 The lower and upper edge lines represent the 5% and 95% of the data, the solid points represent the vertical outliers. The lower and upper quartiles of the boxplots represent 25% and 75% of data, and the solid lines represent the median values and the dash lines represent the average values. Different letters above the boxplot indicate significant differences among the three mentoring periods at 0.05 level.]

常规施肥与不施肥处理相比，小麦和玉米的增产量随种植年限的增加都呈升高的趋势。小麦增产量监测初期 (1988—1997年) 为1971 kg/hm²，之后逐渐升高，监测中期 (1998—2003年) 的增产量为3342 kg/hm²，监测后期 (2004—2016年) 达到最高值4514 kg/hm²，其增产率最高，达到了261%。玉米增产量在监测初期 (1988—1997年) 为1143 kg/hm²，监测中期增产量为2937 kg/hm²，监测后期 (2004—2016年) 达到最高值4246 kg/hm²，其增产率为123%。

为保证各因素间有可比性，首先对土壤pH、有机质 (SOM)、全氮 (TN)、有效磷 (AP) 和速效钾 (AK) 5个土壤肥力指标原始数据矩阵进行了标准化处理^[16]。然后，运用主成分分析法计算了5个指标对土壤肥力变化的贡献 (表2)。

表2(Table 2)

表2 主成分特征值及其在总变异方差中的占比 Table 2 Principal component eigenvalues (EV) and their proportion in total variance

表2 主成分特征值及其在总变异方差中的占比 Table 2 Principal component eigenvalues (EV) and their proportion in total variance 主成分 Principle component 指标 Index 主成分特征值 EV 占比 (%) Proportion 1 全氮TN 2.50 50.01 2 有机质SOM 1.01 20.24 3 速效钾AK 0.66 13.17 4 有效磷AP 0.60 12.02 5 pH 0.23 4.53 注（Note）：TN—Total N；SOM—Soil organic matter；AK—Arailable K；AP—Available P.

按照5个肥力指标得分系数从大到小排列为TN > SOM > AK > AP > pH，指标权重也以全氮最大，速效钾最小，第1和2主成分 (全氮和有机质含量) 的特征值分别为2.50和1.01，对总方差的贡献率分别为50.0%和20.3%，两者之和达到70.3%。换句话说，排在前面的两个主成分影响了土壤肥力全部变化的70.3%，是影响土壤肥力属性的关键因素。

主成分是原各指标的线性组合，各指标的权数为特征向量；它表示各单项指标对于主成分的重要程度并决定了该主成分的实际意义。根据主成分计算公式，可得到2个主成分与原5项指标的线性组合如下：

综合属性得分是将标准化后的数据带入函数表达式中，计算出每个主成分的得分，然后与其对应的贡献率相乘并求和，即，F = F1 × 49.31% + F2 × 20.85%。各阶段的综合肥力属性得分如图4所示。由计算结果可以看出，1988—1997年和1998—2003年的综合属性得分分别为–0.8和–0.35，而2004—2016年的综合属性得分就增加到了0.09，说明土壤肥力在监测后期得到了改善。

图4(Fig. 4)

图4 不同监测阶段土壤综合肥力得分 Fig. 4 General scores of soil fertility at different monitoring periods

综合29年来土壤养分的变化趋势来看，土壤肥力整体呈上升的趋势。除了土壤pH值外，土壤有机质、全氮、有效磷和速效钾含量基本呈上升趋势。

然而，监测的几个肥力指标的阶段性变化规律不完全一致。监测中期土壤有机质和全氮含量较监测初期有显著增加，而监测后期与监测中期没有显著变化；土壤有效磷含量则随着监测年限的延长不断增加，在监测中期，与监测初期的差异未达到显著水平，而到了监测后期，则显著高于初期；土壤速效钾含量在监测初期和中期变化不大，全氮监测后期迅速增加，与监测中期有了显著差异；土壤pH与速效钾一样，在近10年下降显著。与此相对应的是，小麦、玉米产量也与土壤有机质和全氮含量的阶段性变化出现了一致的规律性。从一个角度也证明了土壤有机质和全氮作为土壤肥力主要影响因素的作用^[17-18]。

监测潮土区土壤肥力因素的阶段性变化与29年间施肥量、肥料种类和比例的变化相关^[19-20]。在本研究的潮土区中肥料的种类主要有有机肥和化肥。与监测初期相比，中期和后期无机养分和总养分投入的量显著增加，特别是化肥投入量增加是影响土壤养分和土壤有机质的重要原因。在监测中期有机养分的施用量从监测初期的523 kg/hm²下降到了270 kg/hm²，导致土壤有机质和全氮在监测中期有所下降。长期施用有机肥能明显提高土壤中有机质含量，且有机质的增加量和年限之间具有良好的正相关性^[7]。土壤有机质含量的变化及变化量的大小与肥料类型、施用量和土壤性质有关^[21-24]。本研究中有机质含量29年来增加了10.02 g/kg。1988—2016年监测点总养分的施用量呈上升的趋势，特别是化肥。三个监测阶段总养分的施用量分别为690、813和865 kg/hm²。1998—2003年土壤有机质略有下降，可能与有机养分的施用量下降有关。

此外，潮土pH呈现下降的趋势，经过29年长期常规施肥，潮土的pH下降了3.99%。结果显示，过量施用氮肥是华北平原潮土pH值下降的重要原因^[25]。本研究结果表明，潮土pH随着时间呈下降的趋势，土壤全氮随时间的增长呈增加的趋势，可能是由于化学氮肥在土壤中转化为铵态氮，铵态氮的硝化作用引起的^[26–27]。因此，潮土培肥应该注意合理平衡增施有机肥，重视秸秆还田，适当减少氮、磷肥的施用^[28]。

作物产量受土壤条件、气候因素、施肥措施、作物品种与管理等多种因素的影响^[29]。本研究表明作物产量主要受土壤肥力和施肥的影响。长期常规施肥条件下，潮土小麦和玉米产量呈增加的趋势，通过对土壤养分含量与作物产量之间的相关性分析，结果表明作物产量与土壤养分含量呈显著的线性相关。施肥是作物获得高产的重要措施之一，也是影响土壤肥力因素变化很重要的因素^[30]。年投入总养分、有机养分、无机养分量分别为851、192和659 kg/hm²，N∶P₂O₅∶K₂O = 1∶0.23∶0.77。磷、钾的投入偏低，但土壤有效磷呈明显的上升趋势，这可能与小麦和玉米带走磷量较少有关^[31]。

从玉米和小麦产量年度变化可以看出，小麦产量的增产幅度明显大于玉米。这可能是大部分肥料都在小麦播种时施入，小麦季养分供应充足，而在玉米季有部分监测点只施入化学氮肥，不施磷钾肥^[31]。三个监测阶段N∶P₂O₅∶K₂O分别为1∶0.60∶0.56、1∶0.58∶0.36和1∶0.43∶0.33。磷和钾的量远低于氮的量。磷钾的缺乏可能是影响其产量的主要原因，对于小麦和玉米，肥料中的磷对产量的影响均较大，所以需要适当提高磷肥的用量，降低氮肥的用量。

29年来，随着化肥用量的增加，潮土区土壤有机质、有效磷、全氮和速效钾含量都呈上升趋势。监测初期与监测中期氮肥用量的大幅度增加，导致土壤有机质和全氮的显著增加，而中期、后期磷肥钾肥投入量的增加，导致后期土壤有效磷和速效钾含量的显著增加。总体而言，土壤肥力在29年来已获得显著的提高。但由于长期大量氮肥的投入，导致土壤pH在监测后期的显著下降，有可能在今后成为影响作物产量的因素。

潮土区土壤全氮和有机质是土壤肥力的主要贡献因子，因此潮土培肥的主要目标是提升全氮和土壤有机质含量。