喀斯特区不同植被恢复措施土壤质量评价

董 茜，王根柱，庞丹波，董 亮，张 梅，

刘玉国6，万 龙1,2，周金星1,2＊

(1. 北京林业大学水土保持学院，林业生态工程教育部工程研究中心，水土保持国家林业和草原局重点实验室，北京 100083；2. 北京林业大学云南建水荒漠生态系统国家定位观测研究站，北京 100083；3. 宁夏大学西北土地退化与生态恢复省部 共建国家重点实验室培育基地，宁夏 银川 750021；4. 宁夏大学西北退化生态系统恢复与重建教育部重点实验室，宁夏 银川 750021；5. 云南省建水县林业和草原科技推广站，云南 建水 654300；6. 中国林业科学研究院荒漠化研究所，北京 100091)

人地矛盾冲突导致超量砍伐森林，使得生态系统退化、土地荒漠化速度加快、土壤肥力降低[1-2]。土壤质量与植被类型密切相关，一方面，林木生长吸收土壤中的养分，土壤理化性质影响植物的生长；另一方面，植物凋落物经过分解将养分归还土壤，改善土壤质量[3]。不同植被恢复措施对土壤养分状况的影响程度不同[4]，土壤质量优劣综合体现于土壤生产性能、持水保水性能等方面[5]。我国滇东南喀斯特类型主要为断陷盆地，气候季节性干旱严重，土壤瘠薄，石漠化现象严重，土壤具有低磷高钙的特征。自2008 年以来，滇东南实施了大面积的石漠化综合治理工程，探究不同恢复措施下的土壤质量对评判植被恢复成功与否具有重要意义。

提取合适的评价指标是土壤质量评价过程的重要环节[6]，最小数据集法被广泛应用于土壤质量评价与检测工作[7-10]。杨文娜等[11]在贵州印江典型喀斯特区利用2 种方法对不同土地利用方式土壤进行质量评价，得出主成分分析法及模糊评价法均能较准确评价土壤质量。孙建等[12]探讨石漠化区4 种生态恢复模式对土壤质量的改善，表明金银花（Lonicera japonicaThunb.） 和 花 椒（Zanthoxylum bungeanumMaxim.）对当地植被恢复有积极作用。尽管先前的研究探索了各种评估土壤质量的指标和方法，但在特殊地质背景下的岩溶生态脆弱区土壤质量评价仍较薄弱。本文以滇东南喀斯特区石漠化治理工程中常用的5 种植被恢复措施为研究对象，深入了解不同植被恢复措施对土壤理化性质、全钙及全镁的影响；建立最小数据集土壤质量综合评价体系，并确定哪种植被恢复措施下土壤质量效果最优，为喀斯特区植被恢复提供理论指导。

1 研究区概况

研究区位于云南建水荒漠生态系统国家定位观测研究站（102°47′7″～102°57′18″ E，23°40′49″～23°44′9″ N），该地区属亚热带季风气候。年平均气温约19.8 ℃，年平均降水量约805 mm，有明显的旱季（11 月至次年4 月）和雨季（5—10月）。该地区具有典型的喀斯特地貌特征，属于喀斯特地堑盆地。土壤主要为石灰质土壤，由石灰岩基底发育而成，呈红色[13]。由于严重的人为干扰，区内土壤和植被退化，形成大面积荒地。自1996 年以来，由于退耕还林工程的实施，该地区大量的人工林用于恢复退化的生态系统。该地区亚热带原始森林大部分被人类破坏，出现了天然次生林。研究区主要造林树种为云南松（PinusyunnanensisFranch.）、 车 桑 子（Dodonaea viscosa(L.) Jacq.）、女贞（Ligustrum lucidumAit.） 、 白 枪 杆 （Fraxinus malacophyllaHemsl.） 、 栓 皮 栎 （Quercus variabilisBlume）、直杆蓝桉（Eucalyptus globulussubsp.maidenii(F. Mueller) Kirkpatrick）、 马 尾 松（Pinus massonianaLamb.）、清香木（Pistacia weinmannifoliaJ. Poisson ex Franchet）等[14]。

2 研究方法

2.1 样品采集

于2016 年6 月，分别选取不同类型林分，包括天然阔叶次生林及人工林。人工林中乔木选取了云南松及马尾松针叶林、直杆蓝桉阔叶林，同时选取了车桑子灌木林，样地基本情况见表1。根据对不同林分类型的树种恢复效果进行对比，结果更具代表性。样地大小为20 m × 20 m，每种植被类型3 个重复，共选取15 个样地进行每木检尺并采样。于2016 年8 月从每个样地中随机挖取3个土壤剖面，移除凋落物层，采集0～10、10～20、20～30 cm3层的土壤样本。将3个剖面的同一土壤层中的样品合并，并通过2 mm 与0.15 mm 孔径网筛进行筛分，进行土壤化学性质及金属元素测定；同时分3层取土壤环刀并称质量，带回实验室进行物理性质测定。

表1 不同恢复措施样地基本情况Table 1 Basic information of different recovery measures

2.2 土壤理化性质及金属元素测定

测定指标包括基础理化性质：土壤密度、总孔隙度、毛管孔隙度、非毛管孔隙度、毛管持水量、饱和含水量、pH 值、有机碳、全氮、全磷、氨态氮、硝态氮、有效磷、速效钾，其参照土壤农化分析方法进行测定[15]。研究区属于岩溶区，母岩多为石灰岩及白云岩，而石灰岩主要由CaCO3组成，白云岩主要由方解石组成，富含钙、镁元素，全钙、全镁对于植物生长发育有较大影响[16]。全钙及全镁采用HNO3-HF-HCl 微波消解法并采用电感耦合等离子体质谱仪测定[17]。

2.3 基于主成分分析法的最小数据集构建

选取不同植被恢复措施0～20 cm 表层[18]土壤指标，利用SPSS 对其理化性质及全钙、全镁进行因子分析并计算指标Norm 值[6]。Norm 值（Nik）计算公式如下：

式中：Nik是第i个指标在特征值≥1 的前k个主成分上的综合载荷[6]；uik是第i个指标在第k个主成分上的载荷[6]；λk是第k个主成分的特征值[6]。

2.4 基于聚类分析法的最小数据集构建

在R 语言中采用Ward 法（欧氏距离）进行聚类分析，将指标分为表征不同土壤质量特征的组别，组内相关性较高的可代替，结合参考文献[15]及野外调查数据，剔除冗余指标。

2.5 土壤隶属函数法评定土壤质量

利用“S”型函数、反“S”型函数及抛物线函数，进行土壤质量指数计算。

土壤有机碳、全氮、氨态氮、硝态氮、全磷、有效磷、速效钾、毛管孔隙度、总孔隙度、毛管持水量、饱和含水量在一定范围内与评价指标及土壤质量成正相关，采用“S”型函数进行计算。

土壤全钙、全镁、非毛管孔隙度在一定范围内与评价指标成负相关，采用反“S”型函数进行计算。

土壤pH 值、密度的指标值具有最佳适宜范围，偏离程度越大则对土壤功能影响越小[6]，故采用抛物线函数进行计算。

式中：u(x)为隶属函数，x为指标实际值，a为指标阈值下限，b为指标阈值上限[6]；a1指标临界值下限，a2指标临界值上限，b1最适宜值上界点，b2最适宜值下界点[6]。

采用主成分分析法中的权重进行土壤质量指数（SQI）计算：

式中：wi是第i个指标的权重，Ni为第i个指标的隶属值[19]。

3 结果分析

3.1 不同植被恢复措施土壤质量特征

图1 表明：毛管持水量、饱和含水量排序为：车桑子灌木林>马尾松针叶林>天然次生林>直杆蓝桉阔叶林>云南松针叶林。车桑子灌木林土壤饱和含水量显著高于其他植被恢复模式，其保水持水效果最佳，云南松针叶林则较差（p<0.05）。除土壤密度外，其他指标均随土层深度增加呈降低趋势。

注：不同大写字母表示不同恢复措施同一土层间差异显著（p＜0.05）；不同小写字母表示同一恢复措施不同土层间差异显著（p<0.05）。下同。Notes：Different capital letters indicate significant difference between different restoration measures and the same soil layer (p<0.05);Different lowercase letters indicate significant difference between different soil layers of the same restoration measures (p<0.05). The same below.

图2 表明：pH 值在5.43～7.1 间浮动，马尾松针叶林及车桑子灌木林土壤偏酸性。土壤有机碳、全氮、氨态氮、硝态氮、全磷整体呈随土壤深度增加而降低趋势（车桑子灌木林氨态氮除外）；有效磷及速效钾呈不规则浮动，车桑子灌木林速效钾显著高于其他恢复措施（p<0.05）。图3 表明，全钙在不同恢复措施间无显著差异，天然林及直干蓝桉阔叶林全镁含量显著高于其他恢复措施（p<0.05）。

图2 不同植被恢复措施土壤化学性质Fig. 2 Chemical properties of different vegetation restoration measures

注：不同字母表示不同恢复措施间差异显著（p＜0.05）。Notes: Different letters indicate significant difference between different recovery measures (p < 0.05).

3.2 全量数据集土壤质量评价结果

利用土壤理化性质及全钙、全镁等16 个指标对土壤质量指数进行计算，得出5 种植被恢复措施的土壤质量指数顺序为：车桑子灌木林（0.496）>马尾松针叶林（0.475）>云南松针叶林（0.460）>然次生林（0.440）>直杆蓝桉阔叶林（0.395）。

3.3 喀斯特区恢复林地最小数据集建立

3.3.1 基于主成分分析最小数据集建立 利用主成分分析法对全量数据集进行因子分析，按照累积方差贡献率≥80%的原则[20]，提取了前4 个公因子。根据因子载荷分为4 组（表2），其中，第1 组土壤有机碳与pH 值、硝态氮、有效磷、速效钾、土壤密度及全钙显著相关（表3），选取土壤有机碳进入最小数据集；总孔隙度与毛管孔隙度、毛管持水量、饱和含水量高度相关（表3），选取总孔隙度进入最小数据集。第2 组氨态氮与全磷显著相关，选取氨态氮进入最小数据集；全镁直接进入最小数据集（表3）。第3 组全氮直接进入最小数据集。第4 组非毛管孔隙度直接进入最小数据集。基于主成分分析法提取的最小数据集包括土壤有机碳、氨态氮、全镁、全氮、总孔隙度、非毛管孔隙度。根据选取出的6 个指标计算土壤质量指数，不同植被恢复措施排名为：马尾松针叶林（0.202）>车桑子灌木林（0.199）>天然次生林（0.191）>直杆蓝桉阔叶林（0.183）>云南松针叶林（0.149）。

表2 主成分因子载荷矩阵及因子权重Table 2 Principal Component Factor Load Matrix and Factor Weight of Soil Physicochemical Properties

3.3.2 基于聚类分析最小数据集建立 基于聚类分析法中Ward 法（欧氏距离）聚类，对评价指标进行分类。图4 表明：在聚合水平0～3.5 之间，指标可明显划分为4 类。第一类土壤密度与毛管孔隙度、总孔隙度、毛管持水量存在相关性，合适的土壤密度有助于植物生长，是最常用的物理指标[21]，选取土壤密度进入最小数据集。全磷及有效磷中，对于短期土壤质量有明显变化的为有效磷[22]，且其属于强变异，故选取有效磷进入最小数据集。非毛管孔隙度、全镁与组内其他因子无明显相关性，其中，非毛管孔隙度与土壤密度评价效果重叠，故将其剔除，全镁进入最小数据集。第二类速效钾与全钙均进入最小数据集。第三类全氮与氨态氮存在相关性（表3），全氮、硝态氮及氨态氮在植物生长发育过程中，起作用较大的为硝态氮及氨态氮[23]，在pH 值较低和氮限制下，植物对于氨态氮吸收更好[20]，选取氨态氮、pH 值进入最小数据集。第四类有机碳与饱和含水量高度相关（表3），碳含量处在合理范围内，可以对饱和含水量进行有效调控[24]，选取有机碳进入最小数据集。基于聚类分析法提取的最小数据集包括土壤密度、有效磷、全镁、速效钾、全钙、氨态氮、pH 值、有机碳。根据选取出的8 个指标计算的土壤质量指数排序为：车桑子灌木林（0.292）>直杆蓝桉阔叶林（0.259）>马尾松针叶林（0.226）>云南松针叶林（0.216）>天然次生林（0.211）。

3.3.3 最小数据集合理性验证 将全量数据集土壤质量指数与最小数据集土壤质量指数作散点图进行回归分析，从拟合效果看（图5），全量数据集-土壤质量指数与主成分分析-土壤质量指数的相关系数为0.843 5，全量数据集-土壤质量指数与聚类分析-土壤质量指数的相关系数为0.767 9。滇东南喀斯特区主成分分析法更适用于提取最小数据集指标进行土壤质量评价。

图4 不同植被恢复措施土壤质量评价指标聚类分析Fig. 4 Cluster analysis of soil quality evaluation indicators for different vegetation restoration models

图5 最小数据集与全量数据集的相关性Fig. 5 The correlation between the smallest data set and the full data set

表3 土壤理化性质及全钙全镁相关性矩阵Table 3 Physical and chemical properties and correlation matrix of total calcium and total magnesium

4 讨论

4.1 不同植被恢复措施土壤理化性质

森林土壤是植被类型、气候因子以及人为干扰等因素综合作用的结果。喀斯特断陷盆地区由于季节性干旱的影响，树种选择不同造成了土壤理化性质的差异。车桑子由于喜光、耐旱、耐瘠薄的特性，能在石灰岩裸露的荒山生长，是喀斯特区较干旱区域大面积恢复的植被树种。车桑子灌木林土壤密度与其他恢复措施存在显著差异，主要因为车桑子灌木林根系较浅，对于浅层土壤物理性质影响较大，表层土壤改良效果明显，土壤密度比其他植被恢复措施有显著降低[25]，但深层土壤受到根系扰动较少，不易改变其物理性质[26]，亟需发展改良深层土壤物理性质的植被恢复措施。车桑子灌木林土壤有机碳、氨态氮、硝态氮、速效钾等指标均处于优势地位，其土壤短期恢复效果较佳，故车桑子作为先锋恢复树种，对土壤改良有良好作用，但是后期车桑子林恢复措施显著降低生态系统植物多样性，需进行物种替代。化学性质是反映不同植被恢复措施改良土壤状态的重要因素[27]，分析土壤化学性质有助于了解土壤因子与群落类型之间的关系[28]。喀斯特断陷盆地不同植被恢复措施间化学性质存在显著性差异。马尾松针叶林及车桑子灌木林偏酸性，同等环境条件下针叶林土壤酸性更强，将限制细菌活动，分解缓慢，易形成粗腐殖质。不同植被恢复措施因其群落结构、凋落物种类及分解速率存在较大差异，从而导致有机质含量不同[29]。土壤有机碳表现为车桑子灌木林>马尾松针叶林>天然次生林>直杆蓝桉阔叶林>云南松针叶林，阔叶林与针叶林间无明显规律性，这与耿玉清等[30]研究不同类型对土壤有机碳影响的结果类似。从土壤全氮、全磷等养分含量看，直杆蓝桉阔叶林全氮含量位居首位，硝态氮及氨态氮含量却缺乏。不同植被恢复措施土壤化学性质仍显现出喀斯特区低磷高钙的特点，全磷及有效磷整体处于较低状态，且植被恢复后，磷元素变化幅度不大，这是因为亚热带酸性土壤中，磷易于形成不溶性磷化合物[3]，因而，断陷盆地喀斯特区磷限制性仍是植被恢复需解决的重要瓶颈问题。森林土壤全镁在不同植被恢复措施间存在显著差异，说明其含量与林分类型密切相关。全镁在天然林及直干蓝桉阔叶林中含量较高，阔叶林林下凋落物层丰富，有利于土壤阳离子交换吸附，造成土壤碱性金属相对累积[3]。

4.2 不同植被恢复措施土壤质量评价

采用全量数据集及最小数据集法，结合土壤质量指数对云南省喀斯特区不同植被恢复措施下土壤质量进行综合评价，最小数据集得出的评价指数与全量评价指数相关系数达显著相关，且基于主成分分析法的最小数据集提取方法相关系性达0.843 5。说明最小数据集评价指标体系均有较好代表性，适宜于喀斯特岩溶区土壤质量评价。喀斯特断陷盆地进行土壤质量评价时从16 个指标减少为6 个指标，分别为有机碳、全氮、氨态氮、总孔隙度、非毛管孔隙度、全镁，大幅减少了工作量，省时省力。

评价结果表明：车桑子灌木林、马尾松及云南松针叶林恢复效果较好。天然林氨态氮、全磷及有效磷含量占明显优势，其群落功能复杂，凋落物分解加快，养分循环与富集功能更强[31]，根据全国第二次土壤普查养分分级标准[6]，天然林全磷含量为一级，有效磷含量虽位于几种恢复措施首位，但仍处于第6 级，磷元素利用效率低，施用适量磷肥有助于植被更好的恢复；人工林中马尾松针叶林土壤有效磷含量较高，可以尝试与马尾松林混交措施来达到更好的恢复效果。云南松林硝态氮及速效钾含量较高，林内枯枝落叶层能有效增加地表的粗糙度，凋落物腐烂可以增加土壤氮、钾含量，起到保水保肥效果等，从而增强了植物群落改良土壤的能力[32]；但其有机碳含量处于较低水平，可能由于其枯落物含不易分解化合物，从而分解周期较长，导致云南松针叶林土壤有机碳返还量减少[33]。车桑子对于土壤短期恢复效果较好，但叶超等[34]研究发现，车桑子密度较大时，本土植物物种丰富度和多样性相对较低，群落结构简单且稳定性较差，呈现退化趋势，车桑子可作为当地恢复先锋物种，长期生态修复需考虑混交种植措施。直杆蓝桉阔叶林对于土壤恢复效果较差，其氨态氮、硝态氮、速效钾均处于最低水平。有研究发现，与邻近群落相比，直杆蓝桉阔叶林表现出对土壤养分的持续耗竭，易造成生态系统土壤养分贫瘠化[35]。侯晓龙等[36]发现，种植桉树会导致土壤表层结皮、生境条件变差，群落层次单一致使物种多样性下降。直杆蓝桉阔叶林全钙、全镁含量较高，岩溶区钙元素迁移及富集作用较强[37]，石灰岩发育的土壤上进行人工植被恢复过程中应注意考虑母岩淋溶作用。在滇东南喀斯特区优先选择马尾松针叶林进行种植，避免直杆蓝桉阔叶林单一种植。

5 结论

（1）研究区植被恢复对表层土壤改良作用较大，深层土壤受到扰动较小。马尾松针叶林及车桑子灌木林土壤偏酸性。车桑子灌木林土壤有机碳、速效钾含量显著高于其他措施，云南松针叶林土壤总孔隙度及速效钾含量低于其他植被恢复措施。

（2）滇东南喀斯特岩溶区基于主成分分析提取的最小数据集由有机碳、全氮、氨态氮、总孔隙度、非毛管孔隙度、全镁组成，与全量数据集得出的土壤质量指数相关性达0.843 5，可反应该区域土壤质量状况。

（3）根据选取出的6 个指标计算土壤质量指数，不同植被恢复措施排名为：马尾松针叶林（0.202）>车桑子灌木林（0.199）>天然次生林（0.191）>直杆蓝桉阔叶林（0.183）>云南松针叶林（0.149）。喀斯特岩溶区植被恢复措施中，马尾松及车桑子恢复效果比天然林好；直杆蓝桉及云南松恢复效果较差。喀斯特区进行人工林恢复时可优先选择马尾松针叶林进行种植；若采取车桑子长期生态修复，宜采取混交种植方式。