-
森林土壤是维持森林生态系统的结构组成和动态平衡的重要因子,它是植物生存的基质,也是物质循环、能量流动、水分平衡等众多生态活动的载体[1]。土壤肥力是土壤理化性质和生物性质的综合反映,影响植被的发生、发育和演替进程,同时植物群落结构及其产生的凋落物以及营造的生态环境和生物环境都会影响生态系统中土壤养分的分布和循环[2-3]。植被的凋落物和枯木的理化性质决定了养分的降解速度和生物有效性,枯枝落叶分解完成后,根系会重新分配养分,从而决定养分积累的垂直和水平变化[4],各个土壤类型在不同的植被类型中均表现出不同的土壤理化性质,并随着植被的演替而发生变化[5-6]。因此,森林土壤的肥力因植被的类型、结构和演替年限的不同而不尽相同[7]。
近年来,人工林所引起的土壤养分循环受阻和生态系统不稳定等诸多问题已引起广泛关注,土壤肥力已成为人工林恢复与重建工程的关键所在[8]。国内外众多学者对不同气候带的不同植被类型人工林土壤的理化性质和土壤质量的研究表明,植被恢复后,不同林分类型土壤养分含量差异显著,土壤改良效果阔叶林优于针叶林,有林地优于荒草地[9-10];不同植被恢复模式下,乡土树种的土壤质量改善能力优于外来树种,造林时应遵循适地适树原则,以植被自然正向演替改良土壤质量[11-12];人工林土壤养分的吸收量和实际归还量不均衡是导致土壤肥力下降重要原因,其中,针叶人工林加剧了土壤地力退化[13-14]。尽管有关植被恢复对土壤理化性质的报道较多,但主要集中在不同的海拔、地形、成土母质、林分起源等因素对森林土壤理化性质的影响,而不同的立地条件可能造成森林土壤的理化性质具有较强的异质性[15]。然而,在南亚热带区域,针对在相同的立地条件下,不同树种在同一时间营造的人工纯林对土壤理化性质的影响却鲜见报道。
为此,本文选取柳州市沙塘林场马尾松采伐迹地上1987年营造的木荷(Schima superba Gardn. et Champ.)、杉木(Cunninghamia lanceolata (Lamb.) Hook.)、蓝果树(Nyssa sinensis Oliv.)、米老排(Mytilaria laosensis Lec.)4种人工纯林为研究对象,以未造林灌草地为对照,旨在分析不同树种组成的人工纯林在经过33年生长后土壤的理化性质和化学计量学特征,揭示在相同立地条件下不同人工林恢复对土壤理化性质的影响规律,并对不同林分的土壤质量进行评价,以期为柳州丘陵地区人工林营造及可持续经营提供科学依据。
-
试验地位于广西柳州市沙塘林场(109°21′~109°28′ E,24°27′~24°34′ N),属南亚热带季风气候区,光照充足,雨热同季,年均气温20.1 ℃,年均日照时间1 634.9 h,年均降水量1 200~1 500 mm,月年均蒸发量1 599.8 mm,年均相对湿度78%,年均无霜期长达357 d。林场海拔区间150~230 m,主要地貌类型为低山丘陵,土壤为砂岩发育成的红壤,土壤深厚且层次不明显。造林地之前为马尾松林分,马尾松皆伐后,在其采伐迹地上营造人工林,初植密度2 500 株·hm−2,现存人工林1 000.5 hm2,其中,木荷1.3 hm2,杉木2.1 hm2,蓝果树2.0 hm2,米老排4.7 hm2,4种林分均为1987年造林;灌草地则是在马尾松采伐迹地上次生演替而来,长年无人为干扰。
-
2020年8月中旬,在马尾松采伐迹地的连续山地选取木荷、杉木、蓝果树、米老排4种林分,对照灌草地遵循就近原则选取,林分平均海拔160 m。在各个林分中至少间隔20 m设置3个代表性样地(20 m × 20 m)。再在每个样地内选取3个土壤取样点,除去表面植被和凋落物,用100 cm3的土壤环刀(土钻:型号Edelman,产地Germany)分别采集不同层次(0~15、15~30、30~50、50~70、70~100 cm)土样,用于测定土壤水分物理性质;同时取500 g左右土样装入样品袋,剔除较大的植物残体和石渣,风干后研磨用于测定土壤化学性质。对样方内胸径≥1 cm的乔木每木检尺,记录种名、胸径、树高,同时在每个样方内随机布设2个2 m × 2 m小样方进行灌草调查,记录灌草的种名、株数、盖度、高度。各林分基本情况见表1,其中,主要灌草植被有罗伞(Brassaiopsis glomerulata (Bl.) Regel)、淡竹叶(Lophatherum gracile Brongn.)、狗脊(Cibotium barometz (L.) J.Sm)、弓果黍(Cyrtococcum patens (L.) A. Camus)、玉叶金花(Mussaenda Pubescens Ait.f.)、半边旗(Pteris semipinnata L.)、小蜡(Ligustrum sinense Lour.)、粗叶榕(Ficus hirta Vahl)等。
林分类型
Stand types坡度
Slope/(°)林分密度
Density/(株·hm−2)郁闭度
Canopy/%胸径 DBH/cm 平均树高
Mean height/m灌草种类
Species/种灌草盖度
Coverage/%Dmin Dmax $\overline {D}$ 木荷 Schima superba 7 459 87 16.8 42.4 29.4 18.7 27 17.3 杉木 Cunninghamia lanceolata 8 1592 93 3.0 31.9 20.2 15.6 19 60.3 蓝果树 Nyssa sinensis 11 792 82 4.4 36.7 20.4 16.4 40 64.7 米老排 Mytilaria laosensis 14 1292 91 3.0 35.3 16.7 14.1 35 11.9 灌草地 Irrigate grassland 13 — — — — — — 52 75.6 Table 1. Basic information of sample plot
-
土壤物理性质测定:土壤密度、最大持水量、毛管孔隙度、总孔隙度参照中华人民共和国林业行业标准-森林土壤分析方法[16]。
土壤化学性质测定[17]:pH值用土水比1∶2.5 pH计法,土壤有机碳用重铬酸钾容量外加热法,全氮含量用半微量凯氏定氮法,全磷含量用王水酸熔—钼锑抗比色法,全钾含量用王水酸熔—火焰光度计法,碱解氮含量用碱解扩散法,有效磷含量用盐酸—氟化铵法,速效钾含量用醋酸铵浸提—火焰光度法,阳离子交换量用乙酸铵静置交换法,土壤水溶性硫酸根离子用分光光度法。
-
用模糊数学原理对数据进行标准化处理,实现对各土壤性质的量纲归一化,使得各个指标具有可比性。当土壤指标对土壤功能呈正相关时,将最大持水量、毛管孔隙度、总孔隙度、土壤有机碳、全N、全P、全K、碱解N、有效P、速效K、阳离子交换量、硫酸根离子分为“S”形,当土壤指标对土壤功能有一个最佳的适宜范围时,将pH值、土壤密度分为抛物线形[18-19]。其公式分别为:
式(1)中:Xi为各个评价指标的隶属度值;X表示各个评价指标原始数值;Xmax和Xmin分别表示各个评价指标的最大值和最小值。
然后,进行KMO和Bartlett球形度检验,检验后的数据经过主成分分析得到成分矩阵、成分得分系数矩阵、特征值、贡献率和公因子方差。提取特征值>1的主成分,根据成分得分系数构建数学模型,如公式(2)。最后根据各主成分的方差贡献率为权重,加权求和。
式(2)中:F为土壤质量综合指数,n为主成分数量;ai为各个因子的得分系数。
-
不同林分类型的土壤理化性质采用Excle 2016进行数据统计,采用SPSS 26进行显著性检验、KMO和Bartlett球形度检验以及主成分分析,用Origin 2018作图。
-
由图1(a)可知:各个林分的土壤密度随着土层加深而增加。5种林地的平均土壤密度从大到小依次为杉木林(1.39 g·cm−3)>蓝果树林(1.37 g·cm−3)>灌草地(1.35 g·cm−3)>木荷林(1.28 g·cm−3)>米老排林(1.25 g·cm−3),但差异不显著(p>0.05)。杉木林每层土壤密度均高于其他林分的同层土壤,说明杉木林林地可能因枯落物少、须根多粗根少等原因,对土壤物理性质(密度)改善不如阔叶树种。
土壤孔隙度衡量土壤透气和水分渗透程度,影响土壤与大气之间的水气交换和植物对土壤中水分和养分的吸收。由图1(c、d)可知:同种林分随着土层深度的增加,毛管孔隙度和总孔隙度都整体上呈现逐渐递减的趋势。不同林分间毛管孔隙度和总孔隙度差别明显,阔叶林的土壤孔隙度大于针叶林,可能是阔叶林的凋落物量大,凋落物的分解有利于土壤动物和微生物的活动,并且阔叶林的根系发达,使得土壤颗粒疏松,透气性和透水性增加;而针叶林土壤透气保水性相对较差,这也证明了杉木林土壤密度大,含水率低。灌草地土壤孔隙度与米老排林相近,表明灌木和草本能够改良土壤孔隙度。
-
土壤的pH值是由土壤母质、气候和林分经营措施综合控制的,与微生物的活动、土壤元素的转移有关,对植物的生长发育有重要影响。图2(a)表明:不同森林类型土壤pH值的平均变化范围为3.94~5.12,土壤呈强酸性或酸性,pH值均随土壤深度的增加而升高。各林分的土壤pH值无明显差异,其中,米老排林的平均pH值最高,为4.77,木荷林的平均pH值最低,为4.09。
土壤有机碳作为林木营养的来源,直接影响和改变土壤的理化性质和微生物特性。图2(b)表明:不同林分的土壤有机碳具有明显差异,且随土壤深度的增加而逐渐减少,0~15 cm层土壤有机碳含量明显高于其它土层,表聚效应明显,其中,木荷林0~15 cm层的土壤有机碳含量高,达18.36 g·kg−1。平均土壤有机碳含量依次为米老排林(9.70 g·kg−1)>木荷林(9.04 g·kg−1)>蓝果树林(5.51 g·kg−1)>灌草地(4.15 g·kg−1)>杉木林(3.26 g·kg−1)。与对照灌草地相比,营造米老排、木荷、蓝果树人工林能分别提高133.7%、117.8%、32.8%的土壤有机碳,而营造杉木人工林对土壤有机碳的贡献还不如灌草地,说明阔叶树对土壤有机碳积累的效应强于针叶树。
-
全N是土壤N的储备指标,说明了土壤供应N元素的能力,碱解N能较好的反应出近期土壤N的供应状况和N的释放速率。图3(a、b)表明:全N含量的范围为0.30~1.43 g·kg−1,碱解N含量的范围为17.17~163.09 mg·kg−1,土壤N素含量差异显著(p<0.05),N素含量阔叶林>针叶林>灌草地,同种林分中随土层深度的增加而减少。米老排林、木荷林、蓝果树林全N含量和碱解N含量在0~15 cm层中差别明显,全N含量分别为1.43、1.33、0.99 g·kg−1;碱解N含量分别为163.09、122.50、90.94 mg·kg−1,但随着土层深度逐渐增加,3种林分的全N和碱解N在土壤70~100 cm层的差别不明显,表明土壤N素异质化的过程受到植物群落组成的制约,米老排林分的N素积累优于其他林分。与对照灌草地相比,营造米老排、木荷、蓝果树、杉木人工林能分别提高117.7%、132.3%、95.4%、33.1%的土壤全N含量。
-
全P是衡量土壤中各种形态P素的总和,有效P则衡量土壤供应P素水平的高低。图4(a)表明:同种林分各土层的全P含量差异不显著(p>0.05),各林分平均全P含量依次为蓝果树林(0.27 g·kg−1)>米老排林(0.24 g·kg−1)>木荷林(0.22 g·kg−1)>灌草地(0.21 g·kg−1)>杉木林(0.19 g·kg−1)。与对照灌草地相比,营造米老排、木荷、蓝果树人工林能分别提高12.9%、3.5%、25.5%(依据原始数据计算)的土壤全P含量,杉木林土壤的全P含量在所有林分中最低,杉木林生长可能消耗更多的P元素。图4(b)表明:5种林分土壤平均有效P的含量为0.15~0.46 mg·kg−1,相同林分不同土层间有效P含量存在显著差异(p<0.05),总体上随着土壤深度的增加而减少,表层有效P含量明显大于深层,各个林分0~30 cm层有效P含量急剧下降,随后下降趋势减缓,木荷林和杉木林0~15 cm层的有效P含量约为深层的2倍。有林地的有效P含量均高于灌草地,其中,蓝果树林的平均有效P的含量最高,为0.26 mg·kg−1,而灌草地的平均有效P含量最低,为0.19 mg·kg−1。
-
土壤K素与植物生存息息相关,速效K含量虽占比极少,但能真实反映土壤中K含量的供应情况。图5(a、b)表明:土壤K素含量在不同林分间存在显著差异(p<0.05),同种林分间差异不显著(p>0.05)。全K含量随着土层深度的增加呈上升趋势,5种林分土壤0~15 cm层的全K含量均低于土壤深层。灌草地全K含量最高,均值为23.81 g·kg−1,蓝果树林和米老排林全K含量相近,均值分别为17.21、19.03 g·kg−1,木荷林和杉木林土壤的全K含量较少,均值分别为10.63、9.83 g·kg−1。米老排林、蓝果树林、灌草地林的速效K含量在0~50 cm层呈下降趋势,之后随土层深度的增加缓慢上升,木荷林和杉木林总体上呈现下降的趋势。与对照灌草地相比,营造米老排、木荷、蓝果树、杉木人工林使土壤全K含量分别减少了20.1%、55.3%、27.7%、58.7%,速效K(除米老排林外)也有相似规律,这种变化规律暗示了人工林对K元素的需求量较大,尤其是木荷林和杉木林,因此,在人工林经营过程中需要对不同树种人工林进行针对性的K元素调控管理。
-
阳离子交换量是指土壤胶体所能吸附各种阳离子的总量,是土壤肥力形成的基础,直接决定了植物产量和品质。图6(a)表明:5种林分的阳离子交换量变幅为6.15~11.31 cmol·kg−1,平均阳离子交换量含量依次为蓝果树林(10.71 cmol·kg−1)>米老排林(10.56 cmol·kg−1)>灌草地(9.60 cmol·kg−1)>木荷林(7.61 cmol·kg−1)>杉木林(6.68 cmol·kg−1),蓝果树林、米老排林、灌草地的阳离子交换量平均含量相近,木荷林、杉木林的阳离子交换量含量明显低于其它林分,各林分不同土层的阳离子交换量无明显差异。根据土壤阳离子交换量>20 cmol·kg−1为保肥能力强、10~20 cmol·kg−1为中等保肥能力、<10 cmol·kg−1为弱保肥能力可知[20],蓝果树、米老排保肥能力中等,灌草地、木荷林、杉木林保肥能力弱。
硫酸盐在自然界中广泛存在,水溶性SO42−的含量会影响土壤含盐量,从而影响林分的生长。由图6(b)可知:平均水溶性SO42−含量依次为蓝果树林(52.49 mg·kg−1)>木荷林(48.71 mg·kg−1)>杉木林(47.31 mg·kg−1)>米老排林(40.19 mg·kg−1)>灌草地(37.57 mg·kg−1)。有林地的水溶性SO42−含量在0~30 cm层相近,都呈上升的趋势,并且50~100 cm层的SO42−含量(除米老排林50~100 cm层外)明显大于表层,而灌草地0~50 cm层显著大于50~100 cm层的SO42−含量,说明环境对不同森林浅层土壤的SO42−含量的影响相同,林分生长吸附SO42−后,经过渗透作用富集在深层土壤中。
-
由表2可知:5种林分的土壤除C:N和C:P比值外,其它化学计量比值存在显著差异。C:N比值的变化范围为9.100~10.359,低于中国土壤C:N的平均范围(10~12)[21];C:P比值的变化范围为25.061~31.245,低于中国土壤C:P比均值(61)[22];N:P比值的变化范围为1.917~4.378,低于中国土壤N:P比均值(5.2)[22]。5种林分相比,木荷林(除C:N外)土壤元素比值均最高,而阔叶林C:N、C:P、N:P比值整体上高于针叶林,针叶林高于灌草地。
比值
Ratios木荷
Schima superba杉木
Cunninghamia lanceolata蓝果树
Nyssa sinensis米老排
Mytilaria laosensis灌草地
Irrigate grasslandC:N 9.100±3.068 a 9.103±2.386 a 9.511±1.860 a 10.359±2.378 a 9.123±2.971 a C:P 30.00±7.024 a 26.835±11.090 a 29.413±9.655 a 31.245±7.965 a 25.061±12.322 a N:P 4.378±1.289 a 2.903±0.633 b 3.057±0.599 b 3.737±1.853 a 1.917±0.609 c C:K 0.933±0.664 a 0.550±0.318 b 0.504±0.220 b 0.532±0.328 b 0.175±0.131 c N:K 0.095±0.039 a 0.058±0.019 b 0.051±0.015 b 0.049±0.023 b 0.017±0.007 c P:K 0.021±0.004 a 0.020±0.003 a 0.017±0.004 b 0.013±0.002 c 0.009±0.001 d 注:数据为“平均值 ± 标准差”;不同小写字母表示同一元素比值在不同林分间差异显著(p < 0.05)。
Notes: the data is "mean ± standard deviation"; Different lowercase letters indicate that the ratio of the same element has significant difference among different stands (p < 0.05).Table 2. Soil stoichiometric characteristics of different artificial forests
-
主成分分析法在可以损失很少的信息的前提下,将不同林分类型的多个土壤理化性质指标降维简化为几个重要指标,能够很好的处理各个指标的多重相关性而又保留独立性[23-24]。通过KMO和Bartlett球形度检验,KMO值等于0.650,说明各个变量间存在较高的相关性,巴特利特球形度检验的近似卡方为631.475,Sig值为0.000,说明各个指标数据能够进行因子降维分析,两项结果检验都表明该数据适用主成分分析。表3表明:公因子方差最小为0.683,说明各个变量都能被很好的表达,按照特征值>1的原则,提取了3个主成分因子,其特征值分别为6.257、4.146、1.880,累积方差累积为87.733%,说明提取的3个主成分可以解释各个林分87.733%土壤指标的综合信息,其中,速效K、阳离子交换量、最大持水量、毛管孔隙度、总孔隙度在第1主成分发挥重要作用;土壤有机碳、全N、全K、碱解N、有效磷P在第2主成分发挥重要作用;全P在第3主成分发挥重要作用。
土壤属性
Soil properties主成分 Principal component 公因子方差
Common factor variancePCA1 PCA2 PCA3 土壤酸碱度 pH value 0.519 −0.677 0.131 0.745 土壤有机碳 Soil organic content 0.635 0.759 0.021 0.979 全氮 Total nitrogen 0.519 0.800 0.172 0.939 全磷 Total phosphorus 0.492 0.044 0.829 0.932 全钾 Total potassium 0.530 −0.804 0.120 0.941 碱解氮 Alkali hydrolyzed nitrogen 0.628 0.717 0.092 0.918 有效磷 Available phosphorus 0.370 0.790 −0.038 0.763 速效钾 Available potassium 0.783 −0.446 0.243 0.871 阳离子交换量 Cation exchange capacity 0.731 −0.298 0.552 0.927 硫酸根离子 Sulfate ion −0.440 0.270 0.645 0.683 土壤密度 Bulk density −0.743 −0.416 0.275 0.802 最大持水量 Maximum water capacity 0.917 −0.044 −0.305 0.937 毛管孔隙度 Capillary porosity 0.859 −0.318 −0.292 0.923 总孔隙度 Total porosity 0.888 −0.220 −0.294 0.922 特征值 Eigenvalue 6.257 4.146 1.880 — 方差贡献率 Variance contribution/% 44.691 29.612 13.431 — 累积方差贡献率 Accumulateed contribution of variance/% 44.691 74.302 87.733 — Table 3. Principal component factor load matrix, factor weight and common factor variance of soil physical and chemical properties
根据3个主成分得分系数建立主成分模型表达式[25],将隶属度值代入表达式,从而求出各主成分的得分(表4),然后按照各主成分的方差百分比占总累积方差的百分率,得到各主成分的权重,最终求得土壤质量综合评价函数,计算公式为:F=0.509 × F1 + 0.338 × F2 + 0.153 × F3。5种林分的土壤质量综合得分排名为:米老排林(0.393)>蓝果树林(0.363)>木荷林(0.353)>灌草地(0.265)>杉木林(0.244)。
林分类型
Frosttypes主成分 Principal component 综合排名
Comprehensive rankF1 F2 F3 F 木荷 Schima superba 0.474 0.053 0.614 0.353 3 杉木 Cunninghamia lanceolata 0.303 −0.002 0.588 0.244 5 蓝果树 Nyssa sinensis 0.472 0.020 0.760 0.363 2 米老排 Mytilaria laosensis 0.685 −0.045 0.388 0.393 1 灌草地 Irrigate grassland 0.646 −0.286 0.212 0.265 4 Table 4. Comprehensive evaluation results of soil quality of different stand types
Effects of Different Plantations on Soil Physical and Chemical Properties and Soil Quality Evaluation in South Subtropical Zone
- Received Date: 2021-12-11
- Accepted Date: 2022-02-23
- Available Online: 2022-06-20
Abstract: