Evaluation of Soil Fertility after Planting Ginseng in Yanbian Prefecture

通过对延边州走访调查，最后选定延吉市、珲春市、敦化市、汪清县、安图县5个地区采集土壤样品，共选取28个采样点来代表延边州参后地土壤状况，各地区土壤类型均为砂壤土，各采样点的详细信息见表 1. 土壤样品为延边州3市2县进行农田栽参(种植的人参均为2年生进行移栽2年)起参后的土壤样品. 采集样品时，利用五点取样法，将上层土壤拨去，采样深度在20 cm左右，混合5点的样品最后留1 kg土壤. 回实验室将土壤样品在室内自然风干，粉碎过筛备用.

土壤pH值和电导率采用5∶1的水土比，震荡后静止用pH计和电导率仪测定；有机质采用外加热常规方法测定；全氮用凯氏法消化，连续流动分析仪测定；铵态氮、硝态氮用氯化钙浸提，连续流动分析仪测定；全磷、全钾用高氯酸-硫酸酸溶法消化，连续流动分析仪测定；速效磷用碳酸氢钠浸提，连续流动分析仪测定；速效钾用乙酸铵浸提，火焰光度法测定. 具体方法参照鲍士旦编著的《土壤农化分析》^[23].

应用Office Excel 2010对原始数据进行记录并整理，采用SPSS 22.0对数据进行方差分析及主成分分析，并用Excel对数据进行隶属函数得分计算，最后进行图形绘制.

为研究延边州土壤肥力状况，对28个采样点的数据进行描述性统计分析，结果见表 2. 结果显示，延边州农田栽参参后地土壤均呈酸性，其中最小值为4.53，酸化现象比较严重. 参后地土壤中的电导率较低，一般适宜人参生长的电导率在100~200 μs/cm之间，而各地区中电导率的最大值为118.67 μs/cm，平均值为69.17 μs/cm，超过100 μs/cm的土壤样品占21%. 各土壤样品中有机质质量分数较高，平均值达到52.11 g/kg，低于30 g/kg的样品仅占7.1%. 全氮质量分数平均值为3.65 g/kg，总体属于氮丰富土壤，其中全氮质量分数极差达6.56 g/kg；硝态氮和铵态氮质量分数的变化幅度为8.41~50.50 mg/kg和10.31~22.78 mg/kg. 在全磷质量分数中来看，变化幅度为0.31~1.78 g/kg；速效磷质量分数平均值为54.40 mg/kg，质量分数较低，其中最小值为15.03 mg/kg，个别地区速效磷处于缺乏状态. 全钾质量分数平均值为5.33 g/kg，低于适宜作物生长的平均水平，低于5.00 g/kg的样品占43%；速效钾质量分数的变化幅度为94.59~378.36 mg/kg，其中低于200 mg/kg的样品占54%，总体来看速效钾质量分数较低. pH值的变异系数(CV)为8%，属于低等变异水平(CV＜10%)，其余指标属于中等变异水平.

28个采样点为3市2县的农田栽参参后地土壤，不同地区的土壤状况有所不同(表 3). 结果显示，5个地区的土壤状况均呈酸性，除安图地区的pH值为最低(5.42)外，其余地区间差异无统计学意义(p＜0.05). 汪清地区的电导率最高为80.09 μs/cm，珲春地区的电导率最低为47.66 μs/cm，珲春地区的电导率与其余各地区间差异有统计学意义，其余地区间差异无统计学意义. 有机质质量分数的由高到低依次为敦化、汪清、安图、延吉、珲春，且珲春地区的有机质质量分数显著低于其余各地区. 土壤全氮质量分数由高到低依次为敦化、汪清、安图、延吉、珲春，其中珲春地区的全氮质量分数最小为2.30 g/kg，与延吉差异无统计学意义，但与其余地区差异有统计学意义；延吉地区的硝态氮质量分数最高为29.14 mg/kg，与汪清差异有统计学意义，其中汪清最低为19.50 mg/kg；土壤铵态氮质量分数中，安图质量分数最高为15.17 mg/kg，且显著高于延吉和汪清地区，与敦化和珲春间差异无统计学意义. 土壤全磷质量分数由高到低依次为汪清、延吉、敦化、安图、珲春；土壤速效磷质量分数由高到低依次为汪清、敦化、延吉、安图、珲春. 土壤全钾质量分数中汪清最高为5.60 g/kg，显著高于珲春地区，珲春地区的全钾质量分数最小为4.25 g/kg；在速效钾质量分数中汪清和敦化相差不大，分别为241.54 mg/kg和241.58 mg/kg，其次是延吉地区为201.27 mg/kg，最后为珲春(179.9 mg/kg)和安图(169.86 mg/kg)地区.

土壤肥力评价有很多种方法，主成分分析法也是其中的一种. 主成分分析法是把很多指标简化为几个综合评价的指标，具体可以代表所有的指标成分，从而更加清晰地找到众多成分中重要的部分，从而便于对数据进行分析.

利用SPSS 22.0对28个样品的10个土壤养分指标进行了主成分分析，结果见表 4. 将各项指标的数据经标准化处理后，计算相关矩阵的特征值和所对应的特征向量. 结果表明，前4项综合指标的特征值均大于1，此时累计贡献率为69.581%，可代表原来所有指标的绝大部分信息. 主成分1(F₁)中主要综合了pH(-0.667)、电导率(0.838)、全氮(0.572)、硝态氮(0.583)和全磷(0.771)，代表原始数据信息量的29.720%. 主成分2(F₂)中主要综合了有机质(-0.578)、速效磷(0.540)和全钾(0.625)，代表原始数据信息量的15.092%. 主成分3(F₃)中主要综合了速效钾(0.449)，代表原始数据信息量的12.996%. 主成分4中主要综合了铵态氮(0.565)，代表原始数据信息量的11.772%.

结果见表 5. 根据各主成分的因子负荷和标准化后的养分质量分数数据来计算各主成分的得分，主成分得分(F)计算公式为

式中，FAC为标准化数据，λ为特征值的算数平方根.

综合主成分得分(F_综合)是计算每个主成分得分与其对应贡献率之乘积的总和，即

利用SPSS 22.0对5个地区的10个土壤养分指标进行主成分分析，结果见表 6. 将各项指标的数据经标准化处理后，计算相关矩阵的特征值和所对应的特征向量. 结果表明，前3项综合指标的特征值均大于1，此时累计贡献率为90.387%，可代表原来所有指标的绝大部分信息. 主成分1(F₁)中主要综合了电导率(0.917)、有机质(0.878)、全氮(0.868)、全磷(0.849)、速效磷(0.936)、全钾(0.960)和速效钾(0.797)，代表原始数据信息量的55.821%. 主成分2(F₁)中主要综合了pH值(-0.770)和铵态氮(0.889)，代表原始数据信息量的21.951%. 主成分3(F₁)中主要综合了硝态氮(0.953)，代表原始数据信息量的12.615%.

结果见表 7. 各主成分得分计算公式同2.3.1.1.2.

综合主成分得分(F_综合)计算公式为

根据主成分的综合得分显示，5个地区的土壤肥力状况排名由高到低依次为敦化(1.147)、汪清(0.962)、延吉(0.327)、安图(0.163)、珲春(-2.599).

模糊综合评价法是模糊数学中的隶属度理论将定性评价转化为定量评价，从而对土壤的养分状况做出综合评价的一种方法. 评价时首先根据土壤养分的各个指标进行隶属度值计算，然后根据养分指标的贡献水平确定权重，最后求出土壤养分综合指数^[24].

根据前人研究结果^[25]，结合延边州参后地土壤状况，本研究测定的10种养分指标中，pH值采用抛物线型函数计算隶属度值，见公式(1)，转折点取值为x₁=4.5，x₂=5.5，x₃=6.5，x₄=7.5；其余的养分指标采用S型函数计算隶属度值，见公式(2).

其余养分指标在S型函数中转折点的取值见表 8. 根据建立的隶属度函数公式(1)和公式(2)计算得到各土壤样品中每种养分指标的隶属度值.

每种养分指标对土壤肥力的贡献率不同，因此赋予每种养分指标一定的权重值是十分必要的. 本研究中采用相关系数法来确定每种养分指标的权重系数：首先计算出各养分指标间的相关系数(表 9)，然后计算某种养分指标与其他养分指标之间相关系数绝对值的平均值，该平均值占所有养分指标相关系数绝对值的平均值的总和的百分比，也就是权重系数. 各养分指标的权重系数见表 10.

将每种养分指标的隶属度值乘以权重系数后进行累加，可得到每个土壤样品(图 1)和每个地区(图 2)基于模糊综合评价法的肥力综合指数(INI)，有学者将其分级：高(INI≥0.8)、较高(0.6≤INI＜0.8)、中(0.4≤INI＜0.6)、较低(0.2≤INI＜0.4)、低(INI＜0.2). 按照分级方法，本研究中28个样点中3个为高、14个为较高、10个为中、1个为较低. 因此，按照模糊综合评价法，延边州农田栽参参后地土壤的肥力水平在中等以上. 按照地区由高到低依次为敦化(0.747)、汪清(0.731)、延吉(0.678)、安图(0.649)、珲春(0.404).

将两种评价方法结果进行相关性分析，表明两种方法的评价结果具有极显著的正相关关系(图 3)，说明两种方法可以相互佐证.

人参的生长对土壤的环境要求极高，农田栽参同样面临老参地的问题，延边州的老参地土壤养分上呈现出一定的空间变异性，因此调查分析农田土壤的空间变异性十分重要. 本文通过对延边州3市2县的农田栽参参后地土壤的养分状况进行调查，发现pH值呈低等变异水平，电导率、有机质、全氮、铵态氮、硝态氮、全磷、速效磷、全钾和速效钾均呈中等变异水平. 黄安等^[26]研究发现土壤养分呈现中等强度空间变异，种植不同作物时土壤养分状况同样呈中等强度空间变异水平^[1-2]. 各地区的土壤酸化现象严重，随着人参种植年限的增加，参后地逐渐酸化^[27]，今后可以添加碱性物质来调节土壤的pH值. 相较于农作物种植后的土壤，酸化现象严重，参后地的磷钾肥质量分数较农田土壤多，主要是因为种植人参施入了多量的有机肥，而农田地多为化肥为主^[28].

有众多学者运用主成分分析法对土壤的肥力状况进行评价^[29]，此方法可以将因子标准化，具有更好的准确性. 本文利用主成分分析法选取土壤的10个基本理化性质对延边州28个参后地土壤样品进行土壤评价，其中敦化市的24号样品排名为第一，其次是敦化市的20号样品. 对3市2县的综合得分进行排名，由高到低依次为敦化(1.147)、汪清(0.962)、延吉(0.327)、安图(0.163)、珲春(-2.599). 陈欢等^[30]利用主成分分析法对沙浆黑土的肥力状况进行了评价，提取出2个主成分，反映了原信息量的88.22%，并探究不同施肥效果的排序. 陈玉芹等^[31]同样利用主成分分析和聚类分析法对德宏橡胶林进行了土壤肥力评价，发现要侧重施用磷肥.

模糊综合评价法也是较为常用的一种评价方法，有众多学者采用此方法对土壤进行了肥力评价^{[19, 32-34]}. 如赵蛟等^[35]采用模糊数学评价法对土壤的9种理化指标进行了分析，从而评价竹林的土壤肥力状况，李自林等^[36]采用模糊数学评价法对植烟土壤肥力进行了评价，韦建玉等^[37]采用模糊数学和偏相关分析法对各植烟区的肥力状况进行了评价. 本文通过模糊综合评价法发现：28个样点土壤等级中3个为高、14个为较高、10个为中、1个为较低，3市2县由高到低排名为敦化、汪清、延吉、安图、珲春，同主成分分析方法结论一致. 珲春的得分为0.404，接近较低水平，所以对于后期种植作物时建议施入有机肥料对土壤进行改良，提高土壤肥力.

通过两种方法进行的相关性分析，结果呈现极显著的正相关关系，所以两种方法可以相互佐证，表明两种方法均具有可行性. 农田栽参土壤状况差，参后地更为严重，土壤养分的缺失会造成作物减产，并且短期内无法种植人参，所以如何改善参后地土壤状况、有效进行轮作、缩短连作年限，是未来人参发展的重要研究方向.

延边州pH值为5.78，电导率为69.17 μs/cm、有机质、全氮、硝态氮、铵态氮、全磷、速效磷、全钾、速效钾的平均质量分数分别为52.11 g/kg，3.65 g/kg，24.90 mg/kg，14.12 mg/kg，0.94 g/kg，54.40 mg/kg，5.33 g/kg，211.28 mg/kg. 除pH值为低等变异水平外(CV=8%)，其余指标均呈中等变异水平.

通过主成分分析法提取出3个主成分代表所有信息，反映了原信息量的90.387%. 3市2县的土壤肥力综合得分由高到低依次为敦化(1.147)、汪清(0.962)、延吉(0.327)、安图(0.163)、珲春(-2.599).

通过模糊综合评价法对3市2县的28个样点进行排名，由高到低依次为敦化(0.747)、汪清(0.731)、延吉(0.678)、安图(0.649)、珲春(0.404).

两种方法的评价结果呈极显著正相关关系，可以相互佐证，表明延边州农田参参后地土壤状况为中等.