-
土壤是植物生长代谢必不可少的物质基础,其性质与结构受气候、母质类型、地形、生物、人类活动等因素的影响[1]。土壤肥力是植物与生长环境相协调的产物,由土壤各基本属性共同反应而成,能供给植物生长所需的各种养分和水分[2-3],植物生长发育以及果实产量和品质直接受土壤肥力高低影响[4],因此对土壤肥力的空间异质性和分布规律进行合理分析和客观评价,才能科学合理地施肥,进一步提高果实的产量及品质。
大别山山核桃(Carya dabieshanensis M. C Liu & Z. J. Li)属胡桃科(Juglandaceae)山核桃属( Carya Nutt.),是我国特有干果和木本油料树种[5],特产于皖豫鄂三省交界的大别山山区(29~31° N,115~116° E),主要包括安徽金寨县、霍山县,湖北罗田县和河南商城县等地。前人对大别山山核桃的分类学、坚果表型和油脂组分、遗传多样性进行系统研究[6-9],但对土壤养分的研究较少,缺乏合理有效的施肥方案。研究表明土壤肥力指标能够反映土壤养分状况、植物吸收养分能力以及环境对植物的影响,苗庆选[10]等分析了不同生态区核桃林地土壤全氮、有机质、碱解氮、有效磷和速效钾含量等肥力特征的差异;杨慧思[11]等探讨了氮、磷、钾等影响山核桃林地土壤肥力的因子;袁紫倩[12]等对山核桃林地土壤有机质、碱解氮、有效磷、速效钾等进行主成分分析,探明了影响山核桃品质的土壤主要养分元素。本试验以大别山山核桃主产区金寨县、霍山县、罗田县的10个乡镇共88个土壤样品为研究对象,通过测量土壤pH值、有机质、全氮、全磷、全钾、碱解氮、有效磷和速效钾8个土壤肥力指标,解析土壤肥力状况,利用模糊数学隶属函数来计算土壤肥力综合指数(IFI),综合客观评价土壤肥力,以期为大别山山核桃林地土壤高效管理和科学施肥提供依据。
-
研究区域包括安徽金寨县、霍山县以及湖北省罗田县,属于北亚热带湿润季风气候,金寨县年平均气温15.6 ℃,极端最高气温41.7 ℃,极端最低气温−13.9 ℃,无霜期平均228 d,年平均降雨量1 400 mm;霍山县年平均气温15.0 ℃,极端最高气温43.0 ℃,极端最低气温−17.0 ℃,无霜期平均220 d,年平均降雨量1 366 mm;罗田县年平均气温16.4 ℃,极端最高气温41.6 ℃,极端最低气温−14.6 ℃,无霜期平均240 d,年平均降雨量1 330 mm。
-
土壤样品采集位置见表1,分别为安徽省金寨县的吴家店镇、果子园乡、古碑镇、关庙乡、沙河乡、燕子河镇、长岭乡、天堂寨镇,霍山县的太阳乡和湖北省罗田县的九资河镇。根据大别山山核桃投产林的分布情况确定采样林分,每个林分根据面积确定采样数量,在林分中按“S”形路径布点,采样深度为0~30 cm,采用四分法去除作物残体和杂物后留取1 kg土样,并记录采样位置、日期等信息。土壤样品带回实验室风干、磨碎、过筛后用塑料封口袋密封后置于4 ℃冰箱避光保存备用。
编号
Number乡镇
Towns纬度(N)
Latitude/(°)经度(E)
Longitude/(°)海拔
Altitude/m取样数/(份)
Sample numberJGM 安徽省金寨县关庙乡 31°30′19.3"~31°28′23.6" 115°25′20.8"~115°22′53.3" 528~827 11 JWJD 安徽省金寨县吴家店镇 31°09′17.3"~31°13′59.2" 115°35′00.1"~115°38′42.8" 654~806 9 JSH 安徽省金寨县沙河乡 31°20′40.7"~31°22'11.0" 115°23′06.0"~115°28'51.7" 450~825 18 JTTZ 安徽省金寨县天堂寨镇 31°15′13.0" 115°50′53.6" 644~691 6 JGZY 安徽省金寨县果子园乡 31°19'49.6" 115°39'17.0" 658~754 4 JGB 安徽省金寨县古碑镇 31°19'54.0"~31°19'57.0" 115°46'35.0"~115°47'58.5" 583~920 11 JYZH 安徽省金寨县燕子河镇 31°16'14.9"~31°21'15.2" 115°49'43.5"~116°02'31.4" 559~720 13 JCL 安徽省金寨县长岭乡 31°08'59.6" 115°55'45.9" 722~752 4 HTY 安徽省霍山县太阳乡 31°07′38.0" 116°10′27.2" 664~870 5 LJZH 湖北省罗田县九资河镇 31°10′12.7" 115°36′10.7" 887~929 7 Table 1. Soil sampling of forest land in main production area of C. dabieshanensis
-
土壤pH值采用电位法;有机质含量采用外加热法;全氮含量采用半微量凯氏法;土壤碱解氮含量采用碱解扩散法;全磷和有效磷含量采用钼锑抗比色法;土壤全钾和速效钾含量采用火焰光度法[13]。
-
大别山山核桃林地土壤各肥力指标评价分级标准根据第二次全国土壤普查并依据有关资料制定(表2),同时参考前人的研究成果[14-16],构建林地土壤综合肥力评价指标体系。
指标
Target等级 Grade 强酸性/极丰富
Strong acidity/
Extremely rich酸性/丰富
Acidity/
Rich弱酸性/中等
Faintly acid/
Moderate中性/缺乏
Neutral/
Poor弱碱性/很缺乏
Weakly alkaline/
Very poor碱性/极缺乏
Basicity/
Extremely poorpH <4.5 4.5~5.5 5.5~6.5 6.5~7.5 7.5~8.5 >8.5 有机质 Organic matter/(g·kg−1) >40 30~40 20~30 10~20 6~10 <6 全氮 Total N/(g·kg−1) >2 1.5~2 1~1.5 0.75~1 0.5~0.75 <0.5 全磷 Total P/(g·kg−1) >1 0.8~1 0.6~0.8 0.4~0.6 0.2~0.4 <0.2 全钾 Total K/(g·kg−1) >25 20~25 15~20 10~15 5~10 <5 碱解氮 Available N/(mg·kg−1) >150 120~150 90~120 60~90 30~60 <30 有效磷 Available P/(mg·kg−1) >40 20~40 10~20 5~10 3~5 <3 速效钾 Available K/(mg·kg−1) >200 150~200 100~150 50~100 30~50 <30 Table 2. Classification standards of soil nutrient status
-
不同肥力指标的权重Wi通过pearson相关分析法确定,公式如下:
$ Wi=\frac{r{i}_{均}}{{\displaystyle \sum r{i}_{均}}} $ ,其中,ri均——第i个土壤肥力指标与其他指标间相关系数的均值。 -
结合实际土壤状况,确定大别山山核桃林地土壤肥力各评价指标的隶属度函数拐点的取值(表3)。并计算出各项肥力指标的隶属度值N。pH值采用峰值型隶属度函数,有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾采用戒上型隶属度函数。
转折点
Turning pointpH 有机质
Organic matter/
(g·kg−1)全氮
Total N/
(g·kg−1)全磷
Total P/
(g·kg−1)全钾
Total K/
(g·kg−1)碱解氮
Available N/
(mg·kg−1)有效磷
Available P/
(mg·kg−1)速效钾
Available K/
(mg·kg−1)X1 4.5 6 0.5 0.2 5 30 3 30 X2 5.5 40 2 1 25 150 40 200 X3 6.5 X4 7.5 Table 3. Value of inflection point of membership function curve
峰值型隶属度函数:
戒上型隶属度函数:
其中,x1、x2、x3和x4分别为相应函数的拐点取值。
-
土壤肥力综合指数:IFI=
$\sum\limits_{{\text{i}} = 1}^{\text{n}} {{{NiWi}}}$ ,其中,Wi——第i个土壤肥力指标的权重,Ni——第i个土壤肥力指标的隶属度。结合大别山山核桃林地土壤的实际情况,本研究采用等间距法将林区土壤养分状况分为5个土壤肥力等级:IFI>0.80,Ⅰ级,土壤肥力水平高;0.6<IFI<0.8,Ⅱ级,土壤肥力水平较高;0.4<IFI<0.6,Ⅲ级,土壤肥力水平中等;0.2<IFI<0.4,Ⅳ级,土壤肥力水平较低;IFI<0.2,Ⅴ级,土壤肥力水平低。IFI取值在0到1之间,值越接近1,代表土壤肥力越好。
-
使用Microsoft Excel 2019和 SPSS 20.0软件系统对样品数据进行统计分析。
-
大别山山核桃林地土壤pH值介于4.34~6.55,平均值为5.19,属酸性土(表2),变异系数为7.32%(表4)。不同乡镇间土壤pH差异不显著,除JGZY(4.87)和JGB(4.99)外,其余乡镇均大于5.0。由图1可知,全区强酸性、酸性、弱酸性和中性土壤分别占3.41%、71.60%、23.86%和1.14%,无弱碱性和碱性土壤。
乡镇
Town土壤肥力指标
Soil fertility indexpH 有机质
Organic matter/
(g·kg−1)全氮
Total N/
(g·kg−1)全磷
Total P/
(g·kg−1)全钾
Total K/
(g·kg−1)碱解氮
Available N/
(mg·kg−1)有效磷
Available P/
(mg·kg−1)速效钾
Available K/
(mg·kg−1)JGM 5.18 ± 0.13 a 58.20 ± 6.97 ab 1.22 ± 0.16 ab 0.61 ± 0.09 b 9.04 ± 1.00 a 206.17 ± 20.43 abc 19.23 ± 5.65 ab 92.16 ± 5.82 a JWJD 5.25 ± 0.11 a 58.73 ± 6.15 ab 0.76 ± 0.09 ab 1.09 ± 0.16 ab 8.34 ± 0.45 a 172.38 ± 12.33 bc 21.95 ± 3.21 ab 68.86 ± 6.59 abc JSH 5.40 ± 0.10 a 67.47 ± 6.89 a 0.86 ± 0.08 ab 0.87 ± 0.13 ab 7.54 ± 0.42 a 226.70 ± 25.78 ab 20.11 ± 2.61 ab 82.69 ± 7.99 ab JTTZ 5.35 ± 0.12 a 56.67 ± 5.69 ab 1.18 ± 0.19 ab 1.40 ± 0.05 a 8.32 ± 0.42 a 154.70 ± 21.34 bc 36.76 ± 7.31 a 68.56 ± 6.59 abc JGZY 4.87 ± 0.24 a 44.36 ± 12.93 ab 1.24 ± 0.23 ab 1.26 ± 0.24 ab 6.91 ± 0.79 a 136.66 ± 54.60 bc 20.56 ± 6.99 ab 47.80 ± 12.18 c JGB 4.99 ± 0.09 a 42.14 ± 5.60 ab 1.30 ± 0.11 a 0.93 ± 0.15 ab 9.37 ± 0.96 a 155.98 ± 11.23 bc 21.11 ± 5.01 ab 60.88 ± 6.22 abc JYZH 5.14 ± 0.12 a 34.33 ± 3.49 b 0.89 ± 0.08 ab 1.35 ± 0.18 a 6.71 ± 0.68 a 116.37 ± 11.27 c 22.96 ± 3.82 ab 55.36 ± 6.01 bc JCL 5.13 ± 0.02 a 64.56 ± 3.30 a 0.86 ± 0.09 ab 0.83 ± 0.17 ab 8.52 ± 1.43 a 274.96 ± 31.41 a 9.33 ± 3.41 b 85.24 ± 4.52 ab HTY 5.22 ± 0.13 a 50.75 ± 10.37 ab 1.21 ± 0.29 ab 0.82 ± 0.28 ab 10.07 ± 1.85 a 180.60 ± 37.35 abc 17.79 ± 7.72 ab 84.71 ± 15.11 ab LJZH 5.06 ± 0.11 a 50.22 ± 8.26 ab 0.98 ± 0.10 ab 0.96 ± 0.22 ab 8.93 ± 1.31 a 168.91 ± 34.59 bc 20.58 ± 7.36 ab 59.34 ± 5.52 bc 林区范围
Range of
forest area4.34~6.55 10.88~121.34 0.37~2.16 0.19~2.46 2.59~16.55 51.18~552.95 3.03~62.20 25.87~160.52 平均值
Mean5.19 53.11 1.02 0.99 8.24 179.41 21.30 71.52 标准差
SD0.38 23.30 0.42 0.54 2.66 82.20 14.91 26.70 变异系数
CV/%7.32 43.87 40.78 54.37 32.34 45.82 70.00 37.34 级别
Grade2级(酸性) 1级(极丰富) 3级(中等) 2级(丰富) 5级(很缺乏) 1级(极丰富) 2级(丰富) 4级(缺乏) 注:不同小写字母表示显著性差异(P<0.05)
Note: Different lowercase letters indicate significant differences (P<0.05)Table 4. pH、organic matter and available nutrient contents of soil of C. dabieshanensis forest
-
土壤有机质含量介于10.88~121.34 g·kg−1,平均值为53.11 g·kg−1,有机质含量丰富,变异系数为43.87%,属于中等变异。不同乡镇间土壤有机质含量存在显著差异(P<0.05),其中以JSH土壤有机质含量最高,达到67.47 g·kg−1,其次为JCL(64.56 g·kg−1)、JWJD(58.73 g·kg−1)和JGM(58.20 g·kg−1),JYZH最低(仅为34.33 g·kg−1)。全区有机质主要集中于极丰富水平占比为73.86%,丰富、中等、缺乏的土壤分别占9.09%、13.64%和3.41%,无有机质极缺乏、很缺乏土壤(图1)。
-
土壤全氮的含量介于0.37~2.16 g·kg−1,均值1.02 g·kg−1,变异系数40.78%,属于中等变异。不同乡镇间土壤全氮含量差异不显著(P>0.05),其中以JGB含量最高,为1.30 g·kg−1,其次是JGZY(1.24 g·kg−1)和JGM(1.22 g·kg−1),JWJD最低(仅为0.76 g·kg−1)。全区全氮极丰富和丰富水平较少,占比为4.55%和9.09%,中等、缺乏、很缺乏、极缺乏的土壤分别占30.68%、28.41%、21.59%和5.68%。土壤碱解氮含量变化范围51.18~552.95 mg·kg−1,平均值为179.41 mg·kg−1,碱解氮含量总体较高,变异系数为45.82%,变异水平中等。不同乡镇间土壤碱解氮含量存在显著差异(P<0.05),其中以JCL含量最高,为274.96 mg·kg−1,其次为JSH(226.70 mg·kg−1)和JGM(206.17 mg·kg−1),均大于200 mg·kg−1;JYZH最低,仅为JCL的42.32%。全区土壤碱解氮含量主要分布于极丰富和丰富水平,分别占60.23%和19.32%,无极度缺乏的土壤。
-
土壤全磷含量介于0.19~2.46 g·kg−1,均值为0.99 g·kg−1,变异系数54.37%,属于中等变异。不同乡镇间土壤全磷含量存在显著差异(P<0.05),其中以JTTZ含量最高,为1.40 g·kg−1,其次是JYZH(1.35 g·kg−1)和JGZY(1.26 g·kg−1),JGM最低,为0.61 g·kg−1。全区全磷主要分布于极丰富水平,占47.73%,其余分别占6.82%(丰富)、14.77%(中等)、19.32%(缺乏)、10.23%(很缺乏)和1.14%(极缺乏)。土壤有效磷含量介于3.03~62.20 mg·kg−1之间,平均为21.30 mg·kg−1,变异系数为70.00%,变异幅度较大。不同乡镇间土壤有效磷含量存在显著差异(P<0.05),其中以JTTZ最高,达到36.76 mg·kg−1,JWJD、JSH、JGZY、JGB、JYZH、LJZH之间含量相近,JCL含量最低,为9.33 mg·kg−1,仅为JTTZ的25.38%。林区土壤有效磷含量主要集中于丰富和中等水平,均占29.55%,无极度缺乏的土壤。
-
土壤全钾含量介于2.59~16.55 g·kg−1,平均为8.24 g·kg−1,变异系数为32.34%,变异程度中等。不同乡镇间土壤全钾含量差异不显著(P>0.05),其中以HTY最大,达到10.07 g·kg−1,其次是JGB(9.37 g·kg−1)和JGM(9.04 g·kg−1),JGZY和JYZH最低,分别为6.91 g·kg−1和6.71 g·kg−1。林区土壤全钾含量主要分布于很缺乏水平,占比为69.32%,全区无丰富和极丰富水平。土壤速效钾含量变化范围25.87~160.52 mg·kg−1之间,平均为71.52 mg· kg−1,变异系数为37.34%,与土壤全钾变异程度一致。不同乡镇间土壤速效钾含量存在显著差异(P<0.05),其中以JGM最高,达到92.16 mg·kg−1;JGZY最低,仅为JGM的51.87%。林区土壤速效钾含量相对较匮乏,分别占59.09%(缺乏),26.14%(很缺乏),3.41%(极缺乏),10.23%(中等),丰富仅占1.13%。
-
隶属度越接近于理想值1,表示土壤肥力越高。由雷达图(图2)可知,有机质的平均隶属度值最高,达0.91;碱解氮次之,为0.88,全钾的平均隶属度值最低,仅为0.25,其次为速效钾,为0.31。根据土壤各养分指标权重(表5)可知,有机质权重值最大,为0.29;其次为速效钾和有效磷,分别为0.22和0.20,其中权重低于0.10的有pH、全磷和全钾,说明土壤肥力的提高更依赖于有机质、碱解氮、有效磷和速效钾含量的增加。
pH 有机质
Organic
matter全氮
Total N全磷
Total P全钾
Total K碱解氮
Available N有效磷
Available P速效钾
Available K相关系数均值
Mean value of
correlation
coefficient权重
weightpH 1 −0.07 −0.13 0.02 0.10 −0.17 0.14 0.13 0.008 1 0.016 9 有机质 Organic matter 1 0.27* −0.07 −0.06 0.65** −0.06 0.36** 0.138 0 0.287 2 全氮 Total N 1 0.04 0.08 0.13 0.10 −0.02 0.061 3 0.127 5 全磷 Total P 1 −0.17 −0.34** 0.63** −0.19 0.006 4 0.013 4 全钾 Total K 1 −0.02 0.12 −0.03 0.000 4 0.000 9 碱解氮 Available N 1 −0.25** 0.48** 0.063 9 0.132 9 有效磷 Available P 1 0.02 0.096 0 0.199 8 速效钾 Available K 1 0.106 4 0.221 5 注:**表示相关性达极显著水平(P<0.01);*表示相关性达显著水平(P<0.05)
Notes: ** indicates a highly significant level of correlation (P<0.01);* indicates a significant level of correlation (P<0.05)Table 5. Correlation coefficient and weight of soil fertility indicators
-
相关性分析表明(表5),有机质与全氮呈显著正相关(P<0.05),与碱解氮、速效钾呈极显著正相关(P<0.01),全磷与碱解氮呈极显著负相关(P<0.01),与有效磷呈极显著正相关(P<0.01),碱解氮与有效磷呈极显著负相关(P<0.01),与速效钾呈极显著正相关(P<0.01),其他指标间相关性均不显著,其中土壤的pH具有很好的独立性。
-
选取土壤pH值、有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾8个指标对大别山山核桃林地土壤进行肥力综合评价(表6),10个乡镇的土壤综合肥力指数(IFI)变幅范围是0.52~0.72,平均值为0.61,各乡镇IFI由大到小依次为:JTTZ(0.719,Ⅱ级)>JGM(0.675,Ⅱ级)>JWJD(0.653,Ⅱ级)>JSH(0.648,Ⅱ级)>JGB(0.612,Ⅱ级)>JCL(0.569,Ⅲ级)>JGZY(0.566,Ⅲ级)>HTY(0.555,Ⅲ级)>JYZH(0.553,Ⅲ级)>LJZH(0.523,Ⅲ级)。大别山山核桃天然分布区内50.00%的土壤肥力较高,其余处于中等水平,大别山山核桃土壤肥力综合评价等级为Ⅱ级,土壤综合肥力较丰富。
土壤综合指数
IFI肥力等级
Fertility grade排序
SortJGM 0.675 Ⅱ 2 JWJD 0.653 Ⅱ 3 JSH 0.648 Ⅱ 4 JTTZ 0.719 Ⅱ 1 JGZY 0.566 Ⅲ 7 JGB 0.612 Ⅱ 5 JYZH 0.553 Ⅲ 9 JCL 0.569 Ⅲ 6 HTY 0.555 Ⅲ 8 LJZH 0.523 Ⅲ 10 均值 Mean 0.607 Ⅱ Table 6. Comprehensive index and evaluation grade of soil fertility
-
通过主成分分析对土壤各肥力指标与IFI之间的载荷情况进行分析(表7),根据特征值大于1共提取到4个主成分,方差贡献率分别为28.20%、23.27%、13.43%和12.49%,累计方差贡献率为77.39%,基本能反应土壤肥力的基本信息。IFI与有机质、碱解氮在第一主成分(PC1)上占主要决定因素,载荷值均大于0.70;第二主成分(PC2)中主要决定因子为全磷和有效磷,分别为0.82和0.90;第三主成分(PC3)pH的载荷值为0.79,表明土壤酸碱度对肥力的影响;全钾在第四主成分(PC4)中具有较大载荷,为0.85。
主成分
Principal componentPC1 PC2 PC3 PC4 pH −0.04 0.25 0.79 −0.05 有机质 Organic matter 0.83 −0.10 −0.09 −0.10 全氮 Total N 0.42 0.24 −0.46 0.51 全磷 Total P −0.11 0.82 −0.18 −0.29 全钾 Total K −0.02 0.03 0.36 0.85 碱解氮 Available N 0.77 −0.44 −0.03 −0.05 有效磷 Available P 0.10 0.90 0.09 0.02 速效钾 Available K 0.65 −0.13 0.45 −0.23 土壤综合指数 IFI 0.79 0.53 0.02 0.05 特征值 Characteristic value 2.55 2.10 1.21 1.12 贡献率/% Contribution rate 28.20 23.27 13.43 12.49 Table 7. Results of principal component analysis of soil fertility indicators
Soil Nutrient Characteristics and Fertility Evaluation of Carya dabieshanensis Forest Land
- Received Date: 2023-03-24
- Accepted Date: 2023-07-24
- Available Online: 2023-10-20
Abstract: