不同植物恢复措施对采石矿废弃地土壤物理性质的改良效果及评价

闫烨琛，赵廷宁，张 艳，王美琪，胡 平，杨 凝，陈 童

（1.北京林业大学 水土保持学院,北京100083；2.北京林业大学 水土保持国家林业和草原局重点实验室，北京100083；3.北京林业大学北京市水土保持工程技术研究中心，北京100083）

随着社会经济的发展和工业化程度的提高，人们对砂石料需求迅猛增长，采石场数量也呈激增之势，露天采石场的特点是规模小、数量多、距城市近、对环境影响的范围大。众多采石场正面临着十分严重的生态环境退化问题，如环境恶化、水土流失、生物多样性锐减等［1］；而且在开山采石的过程中造成山体破碎，岩石裸露，边坡稳定失衡，容易产生滑坡、泥石流等地质灾害，对周边居民正常生活产生严重威胁，对社会经济健康持续发展产生不良影响［2］。面对这一状况，急需对矿山废弃地进行治理改良。中国近年来高度重视生态文明建设，倡导绿色发展，对废弃采石场的生态重建越来越重视［3］。国内外对于废弃矿山的生态恢复研究主要集中在恢复策略［4］、 工程措施［5］、 恢复演替［6］、 土壤改良［7］、 绿化基质［8］等方面，但针对采石场这种特殊的矿种类型相关恢复治理措施对土壤物理性质改良效果研究很少，并且少有提出针对采石矿生态重建的高效、适宜措施，针对不同植物措施恢复治理对土壤物理环境变化研究不够重视，因此不同植物措施对采石场土壤物理性质的改良效果需做进一步研究。本研究以北京市周边废弃采石场生态重建过程中不同植被恢复措施为研究对象，通过调查分析恢复治理中不同植物恢复措施对土壤物理性质的影响，探讨不同植物恢复措施与采石矿土壤物理性质改良之间的关系，旨在为采石矿这类特殊的矿区废弃地重构选择适宜的治理措施及对土壤物理性质快速改良提供科学依据。

1 材料与研究方法

1.1 研究区概况

治理区位于北京市房山区河北镇境内，海拔约为200～380 m，属暖温带半湿润半干旱大陆季风气候区。春季干旱多风，夏季炎热多雨，秋季天高气爽，冬季寒冷干燥。年平均气温10.8℃（山区），1月平均气温-4.6℃，7月平均气温24.3℃，年温差30.0℃左右。该地区春秋季节结霜频繁，年平均无霜期约160.0～180.0 d，冬季地面下最大冻土深度56 cm，年平均结冰天数124 d。多年平均降水量为650～700 mm。降水多集中在夏季（6-8月），占全年降水量的80%以上，且汛期多暴雨。采石场于2015年由国土资源部投资对其实施土地复垦工程，首先进行场地平整、压实，然后进行客土回填。客土来源于附近荒地，土壤物理性质基本一致，种植侧柏Platycladus orientalis，油松Pinus tabulaeformis，黄栌Cotinus coggygria和白杨Populus tomentosa等；研究区常见的草本植物有猪毛蒿Artemisia scoparia，狗尾草Setaria viridis，艾草Artemisia argyi，马唐Digitaria sanguinalis，青蒿Artemisia apiacea，地黄Rehmannia glutinosa，曼陀罗Datura stramonium，黄花蒿Artemisia annua，灰菜Chenopodium album等14种植物，其中猪毛蒿、狗尾草、马唐、青蒿、黄花蒿是调查样地各群落的优势种。

1.2 研究方法

1.2.1 试验地选定 试验于2018年7-8月在北京市房山区河北镇露天采石矿生态修复治理区进行实地调查取样，治理区的坡度为0°。调查前根据设计资料及现场踏查，最终选取5种植物（白杨、侧柏、油松、黄栌、油松+黄栌）修复措施（表1），种植每种植物时每穴施加客土。每种措施设置20 m×20 m的标准样地。（1）采样方法。对标准地内的乔木进行每木检尺测量，测量树高、胸径、冠幅、郁闭度；在20 m×20 m的标准样地内设置5个1 m×1 m的草本样方，草本植物调查其种类、数量、平均高和盖度等基本情况；在样地内标准木的树穴内，用100 cm3体积的环刀取得原状土，按0～20和20～40 cm分层进行取样，每层重复3个，并用塑封袋分层采取一定的土样，自然风干。（2）草本植物相关指标计算方法。①丰富度指数［9］：Margalef指数R的计算公式为R=（S-1）/lnN。其中：S为群落中的总种数；N为群落中的个体总数。②多样性指数［10］：Simpson指数D的计算公式为Shannon-Wiener指数H的计算公式为其中:Pi为第i种个体数占群落中所有种个体数的比例；S为群落中总种数。③均匀度指数［10］：Pielou均匀度指数J的计算公式为J=H/Hmax=H/lnS。其中:H为实际观察的种类多样性，Hmax为最大的种类多样性；S为群落中的总种数。（3）土壤物理性质测定指标与方法。采用环刀法测定土壤容重及土壤水分物理性质（土壤含水量、饱和持水量、毛管持水量、田间持水量、土壤孔隙度等）［11］；用筛分法测土壤粒径，将自然风干后的土壤研碎，分别过10.00，2.00，1.00和0.25 mm筛孔，每层进行称量，计算土壤粒径分布占比。

表1 研究区样地基本情况Table 1 Basic situation of the studied sites

1.2.2 数据处理与评价方法 （1）采用Excel 2016和Origin 2018进行数据处理及图表绘制，使用 SPSS 18.0统计软件进行单因素方差分析和主成分分析。（2）土壤物理性质主要采用土壤含水量、土壤容重、饱和持水量、毛管持水量、田间持水量、毛管孔隙、非毛管孔隙7个指标作为评价指标。采用主成分分析法建立土壤物理性质的综合得分评价体系。计算出主成分得分后，再通过综合得分表达式得出综合得分。将综合得分按照高低进行排序，通过综合排序结果衡量其土壤理化性质情况。

2 结果与分析

2.1 草本植物物种多样性分析

如表2所示：侧柏和油松+黄栌措施对应草本层的D，H和J的差异不显著，表明这2种措施修复下的草本层生物多样性相似，且这2种措施治理后林下草本生物多样性相对于其他3种措施的值较大，说明这2种措施治理后草本恢复效果更好。白杨措施修复下的草本层R，D和H相对于侧柏、黄栌和油松+黄栌措施显著偏小（P＜0.05），说明白杨措施修复治理后的林下草本植被单一，群落稳定性较差。

表2 不同植物配置下草本植物多样性比较Table 2 Comparison of herbaceous plant diversity under different plant configurations of the quarry

2.2 土壤物理性质分析

2.2.1 不同植物恢复措施修复下的土壤粒径分布特征 如图1所示：在0～40 cm土层，分布在2.00～10.00 mm的土壤粒径占比显著大于在大于10.00和1.00～2.00 mm（P＜0.05）的土壤粒径占比；修复后的土壤粒径主要分布在0.25～10.00 mm。在0～20 cm土层，分布在0.25～10.0 mm的土壤粒径占比从高到低依次为侧柏、黄栌、白杨、油松+黄栌、油松；在20～40 cm土层，分布在0.25～10.00 mm土壤粒径占比从高到低依次为白杨、侧柏、油松、黄栌、油松+黄栌。可以看出，在0～20 cm土层，侧柏措施修复下的0.25～10.00 mm土壤粒径含量最高，占百分比为89%，白杨措施在20～40 cm土层0.25～10.00 mm土壤粒径含量最高，占百分比为80%；且在0～20 cm土层0.25～10.00 mm土壤粒径含量高于20～40 cm。

2.2.2 不同植物恢复措施修复下的土壤容重特征 土壤容重可以判断土壤的松紧程度，疏松、有团粒结构的土壤容重小，反之，紧实板结的容重大［12］。如图2所示：在0～20和20～40 cm土壤中，侧柏措施修复作用下的土壤容重都显著小于（P＜0.05）其他几种措施修复下的土壤容重；侧柏、黄栌、油松+黄栌措施在0～20 cm的土壤容重都大于20～40 cm的土壤容重。这主要是由于长期植被恢复过程中，主要通过浇灌措施来补充植被所需水分，致使地表土壤板结。在20～40 cm土层白杨、油松2种措施修复作用下的土壤容重大于其他几种措施修复的土壤容重，其值都大于1.5 g·cm-3，说明这2种措施修复下的土壤在20～40 cm深处土壤孔隙数量少，土壤的水分、空气、热量状况较差。

图1 不同植物恢复措施治理后的土壤粒径分布Figure 1 Distribution of soil particle size after different vegetation restoration measures

2.2.3 不同植物恢复措施下的土壤孔隙度特征 如图3所示：土壤孔隙度包括土壤毛管孔隙和土壤非毛管孔隙。土壤孔隙是储存土壤水分的场所，也是土壤微生物等迁移的通道，其大小取决于土壤质地、土壤结构等条件。不同治理措施的总孔隙度在土层深度0～20 cm内变化范围为42.75%～52.10%；在20～40 cm为40.65%～53.60%。侧柏措施治理恢复的矿区土壤在0～20 cm及20～40 cm的土壤空隙度均大于其他几种措施，说明侧柏措施治理恢复的土壤储存水分的能力更好，能更好地为植被生长提供水分条件。同时可以看出：5种不同治理措施下在0～20 cm土层的土壤孔隙度从大到小依次为侧柏、油松、白杨、油松+黄栌、黄栌，在20～40 cm土层土壤孔隙度从大到小依次为侧柏、油松、油松+黄栌、白杨、黄栌。

图2 不同植物恢复措施治理后的土壤容重Figure 2 Soil bulk densities after different vegetation restoration measures

图3 不同植物恢复措施治理后的土壤空隙度Figure 3 Soil porosities after different vegetation restoration measures

2.2.4 不同植物恢复措施修复下的土壤田间持水量的特征 如图4所示：5种不同植物修复措施下，土壤在0～20 cm土层的田间持水量均比20～40 cm土层的田间持水量高，且5种不同治理措施下的土壤田间持水量在0～20 cm和20～40 cm土层的大小规律呈现一致性，从大到小依次为侧柏、黄栌、油松+黄栌、白杨、油松。说明矿区复垦地土壤保水、供水能力改良效果在浅层优于深层。土壤田间持水量的分布特征在不同措施间存在一定的差异，侧柏治理措施恢复下的矿区各层土壤田间持水量都大于其他4种措施，油松治理措施治理后的土壤田间持水量显著小于其他几种措施。

图4 不同植物恢复措施治理后的土壤田间持水量Figure 4 Soil field water-holding capacities after different vegetation restoration measures

2.3 不同植物恢复措施的综合评价

运用主成分分析法，针对采石场5种不同植物措施修复改良土壤物理性质效果进行综合分析［13］。根据特征值和相应的方差贡献率，将7个指标分为2个主成分，其中第1主成分以土壤含水量、饱和持水量、毛管持水量、田间持水量、毛管孔隙、非毛管孔隙为主的指标贡献较大，第2主成分中土壤容重的贡献较大。最终确定主成分得分和综合得分表达式如下：F1=-0.47x1-1.45x2+0.92x3+0.52x4+0.47x5+1.33x6-0.897x7;F2=1.95x1+0.11x2-0.33x3-0.31x4-0.39x5+0.28x6-1.32x7。其中：F1为第1主成分；F2为第2主成分；x1为标准化的土壤含水量；x2为标准化土壤容重；x3为标准化饱和持水量；x4为标准化毛管持水量；x5为标准化田间持水量；x6为标准化毛管孔隙；x7为标准化非毛管孔隙。 综合得分=87.18×F1/100+8.85×F2/100。

根据综合得分信息可知：侧柏、油松+黄栌、油松最后综合得分分别为11.29，11.19和11.07，因此，施行不同治理措施后，土壤物理性质质量从高到低依次为侧柏、油松+黄栌、油松、黄栌、白杨（表3）。可以发现侧柏、油松+黄栌、油松这3种植物措施的最后得分相差不大，即在采石矿的生态恢复中优先依次推荐。

表3 主成分分析评价Table 3 Evaluation of principal component analysis

3 讨论

土壤物理性质是反映土壤潜在生产力的重要指标，它能影响土壤水分含量、热特性、通气性、土体稳定性等，进而影响作物的生长发育［14］。土壤颗粒是土壤最基本的组成单元，土壤颗粒的粒径分布对土壤孔隙、土壤水分和土壤团聚体都有着巨大的影响［15-16］。本研究结果表明：5种不同植被措施恢复下的土壤粒径主要分布在0.25～10.00 mm，土壤有机质是形成土壤颗粒的重要胶结物质，在土壤大粒径的形成过程中发挥着重要作用［17］。由于表层土壤受植被的枯枝落叶和土壤微生物分解作用，使得表层土壤有机质高于深层土壤，因此出现表层土壤的大粒径含量高于深层。土壤容重是反映土壤紧密程度的一个重要指标，直接影响着土壤的孔隙状况［18］。研究发现：白杨、油松2种植被措施恢复治理后的土壤容重在20～40 cm土层值都大于1.5 g·cm-3，这表明这2种植被措施对矿区恢复治理后的土壤容重普遍偏大，空隙度较小，储存土壤水分能力较弱，不利于植物根系生长，这与张志玲等［19］的研究基本一致。这主要是因为白杨、油松2种措施对应的树种在幼林龄期以主根生长为主，侧根较少，对深层土壤的修复效果欠佳；侧柏、黄栌、油松+黄栌措施恢复治理后的表层土壤容重大于深层土壤容重，这主要是因为矿区植被的水分主要靠后期养护灌溉，由于灌溉水的瞬时量大，水压力强，对表层土壤造成严重扰动和结构重组，土壤发生板结、致密，这也就出现了本次研究出现的结果。本研究得出不同措施土壤田间持水量范围为18.18%～26.89%，5种治理措施在深层土壤中的田间持水量都低于表层。这一研究规律与陈晓燕等［20］研究结果一致。由于深层土壤已经与该地原有的岩石风化物充分混合，其机械组成已经发生了变化，使其呈现出较低的田间持水量指标；另外，由于植物生长，其枯枝落叶对表层土壤的改良效果也优于深层，所以呈现出上层土壤的理化性质较好的趋势。利用主成分分析法可知：5种不同恢复措施治理后，土壤物理性质质量从高到低依次为侧柏、油松+黄栌、油松、黄栌、白杨。这与不同措施恢复治理后的矿区土壤物理性质单因素分析结果一致，也与5种不同植物措施恢复治理后的林下草本多样性指标调查分析结果相对应。

4 小结

不同植物恢复措施对采石矿生态修复后，土壤粒径主要集中分布在0.25～10.00 mm；侧柏相对于其他4种植物恢复措施降低土壤容重，增大土壤孔隙度，提高土壤的持水能力更优；5种不同植物恢复措施恢复治理后的表层土壤田间持水量高于深层；根据5种不同植被恢复措施下土壤土壤容重、田间持水量、土壤孔隙度等指标，运用主成分分析法得出矿区土壤物理性质质量从高到低依次为侧柏、油松+黄栌、油松、黄栌、毛白杨。