刘宽梅,周秋文

1 贵州师范大学地理与环境科学学院, 贵阳 550001 2 兰州大学资源环境学院, 兰州 730000

土壤是人类和陆地生物赖以生存的物质基础。在成土过程中,生物是其中最活跃的要素,生物中的绿色植物,通过光合作用制造活体有机质,借助微生物作用,让土壤产生了肥力特性,有效推动了土壤的形成和演化[1]。在不同的森林植被群落条件下,土壤发育的理化性质会有所差别,而不同的土壤理化性质对土壤侵蚀的产生与发生强度都有重要影响[2]。土壤侵蚀会带来土壤肥力下降、土地生产力降低、土地荒漠化等危害[3],而土壤抗蚀性是评估土壤抵抗侵蚀营力破坏的性能[4]。

国内外学者对土壤抗蚀性的研究从用单一指标评价演变成多指标综合评价,还从土壤理化性质角度、统计学角度、侵蚀动力学角度、人工模拟法等来探讨土壤抗蚀性能[5-6]。在土壤抗蚀性研究的发展历程中,研究成果越来越全面和丰富,从不同区域、不同土地利用方式、不同植被覆盖、不同岩石风化物、不同土壤类型、不同管理措施、不同侵蚀模式等多方面来研究土壤抗蚀性[7- 11]。国内对土壤抗蚀性的研究从1960年初开始,研究内容主要集中在分散率、分散系数、侵蚀率和分散度等度量指标[12]。目前我国的土壤抗蚀性研究成果较丰富,尤其对黄土高原地区土壤抗蚀性研究比较全面[13- 15]。

贵州喀斯特地区植被覆盖率低、降水量大、土层薄、土壤蓄水性差,极易发生土壤侵蚀[16],植被要素与土壤侵蚀的关系十分密切[17],因此在贵州喀斯特地区开展不同植被群落下的土壤抗蚀性研究内容有重要意义。目前对于贵州喀斯特地区土壤抗蚀性能的研究主要有,赵洋毅等探讨了花溪地区4种岩性发育土壤的抗蚀性,并比较了4种岩性上发育的自然植被下土壤的抗蚀性差异[18]。吴鹏等研究了黔中杠寨小流域5种主要人工林地的水土保持机能[19]。虽然目前针对喀斯特地区不同林地类型的土壤抗蚀性开展了一定研究工作,也取得了一些显著的进展[20- 23]。但是喀斯特地区地表要素空间异质性极大,气候、岩性等条件的差异将直接或间接影响土壤抗蚀性。鉴于不同环境条件下喀斯特林地土壤抗蚀性存在明显差异,有限的研究不足以说明喀斯特林地土壤抗蚀性的普遍规律,仍需要开展一定的样地研究。

本研究以贵州师范大学花溪校区内生态实验站内阔叶林、针叶林和灌丛为研究对象,选择土壤有机质含量、水稳性团聚体含量、团聚度和分散系数等与土壤抗蚀性有关的指标,对比讨论各植被群落下土壤抗蚀性能的差别,并采用相关分析法,分析该区域中各土壤抗蚀性指标间的相关性及其相互影响下的土壤抗蚀性能,望为贵州喀斯特地区制定合理的水土保持方案提供一定的理论依据。

1 研究方法

1.1 研究区概况

试验样地位于贵州省花溪区,地理坐标为106°37′58″—106°38′26″E,26°22′51″—26°23′20″N,属中亚热带湿润季风气候区,平均海拔1200m,年平均气温14.9℃,无霜期平均246d,无霜期长,全年降雨多集中于4—8月份,年均降雨量约为1229mm,雨量充沛。境内地貌破碎,以山地和丘陵为主,为典型的亚热带喀斯特区域。该研究区内有针叶林、阔叶林和灌丛三种植被型组,乔木平均胸径为15cm,郁闭度为0.85,土壤以黄壤和石灰土为主。

1.2 样地设置

为了减少地域差异(如降水量、温度和气候等因素)对结果产生的影响,不同植被类型的样地均设置在花溪大学城贵州师范大学生态实验站内。样地基本情况见表1。

表1 样地基本概况

图1 采样实地照片Fig.1 Photographs of the samples

1.3 样品的采集

每种植被群落样地中选取3个点,每个样点重复采样三次,在0—10cm和10—20cm土层内,用环刀取原状土,用于测定土壤容重、含水量等性质；在每种样地中以S形布点法分别取0—10cm土层和10—20cm土层的混合土样,然后用四分法将土样分取到大约1kg带回实验室风干备用,用于土壤抗蚀指数、土壤机械组成、土壤有机质和土壤团聚体等指标的测定。

1.4 样品的测定

1.4.1分析方法

土壤含水量测定采用酒精烘烤法；土壤容重测定采用环刀法；土壤有机质测定采用高温外热重铬酸钾氧化-容量法；土壤团聚体测定采用人工筛分法(干筛法和湿筛法)；土壤微团聚体及土壤机械组成的测定采用比重法；土壤抗蚀指数的测定采用土粒浸水实验,其过程为：先将土壤样品风干、筛分；再将筛分后3—6mm粒径的土壤颗粒数100粒,分四次放入盛水容器中的1mm土壤筛上进行实验,水温以实验地水源为准,水恰好浸没土粒,同时记录每1分钟土粒完全散开而崩塌的个数,持续记录10分钟,然后按公式计算抗蚀指数。

1.4.2评价指标

土壤抗蚀性是一个由多个指标共同作用的综合因子,单一指标只能揭示土壤对侵蚀营力的相对敏感水平[24]。因此,只能在一定的控制条件下通过实际测定土壤某些参数作为土壤抗蚀性指标, 继而达到评估土壤抗蚀性的目的。本实验采用以下常用指标综合评价土壤抗蚀性：

(1)

式中,S为土壤抗蚀指数(%)；W为总土粒数(个)；V为崩塌土粒数(个)。

(2)

式中,O为土壤有机质含量(%)；V1为滴定空白时所用FeSO4的量(mL)；V2为滴定土样时所用FeSO4的量(mL)；5.00为所用K2Cr2O7的量(mL);0.8000为1/6K2Cr2O7标准溶液的浓度(mol/L);0.003为1/4 毫摩尔碳的克数(g)；1.724为碳转化为有机质的系数；1.1为校正系数；WT为烘干土重(g)。

(3)

式中,Ad为>0.25mm干筛团聚体含量(%)；md为干筛各级团粒质量(g)；Md为干筛试验土样总质量(g)。

(4)

式中,Aw为>0.25mm水稳性团聚体含量(%)；mw为湿筛各级团粒质量(g)；Mw湿筛试验土样总质量(g)。

(5)

式中St为结构破坏率,Ad,Aw含义同上。

R=a-b

(6)

式中,R为土壤团聚状况(%)；a为>0.05mm微团聚体分析值(%)；b为>0.05mm机械组成分析值(%)。

(7)

式中,R1为团聚度(%),R,a含义同上。

(8)

式中,K1为分散率(%)；a1为<0.05mm微团聚体分析值(%)；b1为<0.05mm机械组成分析值(%)。

(9)

式中,K2为分散系数(%)；a2为<0.001mm微团聚体分析值(%)；b2为0.001mm机械组成分析值(%)。

1.5 数据处理方法

一般的数据整理、各个指标的具体计算、统计图表作图均在Excel中进行,土壤抗蚀指标间的相关分析采用SPSS软件进行处理。

2 结果与分析

2.1 土壤抗蚀指数

由图2可知,三种植被群落下的土壤抗蚀性有所差别,但三种植被类型下的土壤抗蚀指数都比较高,反映出植被覆盖对缓解土壤侵蚀的重要作用。在0—10cm范围内,土壤抗蚀性强弱表现为：阔叶林>灌丛>针叶林。阔叶林表层植物种类多样且覆盖度高,枯枝落叶等物质较多,腐殖质含量高,同时能够充分利用外界环境条件和树种间相互促进的有利作用产生含量较高的有机胶体,从而促成土壤团粒结构的构成,使土层更疏松,透水性增强[25-26]。而灌丛和针叶林表层下枯枝落叶种类和数量较阔叶林少,在生物作用下产生的土壤团聚力较弱,所以土壤抗蚀能力弱于阔叶林。在10—20cm范围内,土壤抗蚀性强弱表现为：灌丛>针叶林>阔叶林。阔叶林内木本植物群落下的腐殖质会从表层土壤向下层急剧减少[1],影响下层土壤胶体的形成,进而影响土壤团粒结构的形成。根据图3所示,针叶林表层的容重在三者中是最低的,说明针叶林表层土壤土质疏松,有机质通过淋溶等作用向下淀积到下层内,为形成较好的土壤结构提供条件,因此下层土壤抗蚀性能有所改善。灌丛下层土壤样品质地粘重,有机质含量高,不易被侵蚀[27]。因此针叶林下层的土壤抗蚀性能高于阔叶林,但灌丛下层土壤的抗蚀性较另外两者高。

2.2 有机质含量

如图4所示,在0—20cm土层范围内,三种植被群落下土壤有机质含量随土层深度增加都呈递减趋势。在阔叶林下0—10cm土层中有机质含量最高,但在10—20cm土层中下降趋势明显,相比下降36.13%。这主要是因为阔叶林中植物种类多样,表层生物量流动大,物质循环较快,枯落物容易分解转化,致形成含量较高的土壤活性有机质[25],但喀斯特地区土层较薄,下层生物活动量锐减,枯落物在表层大部分就已经被分解,土壤肥力较表层差,所以产生的有机质下降趋势明显。草本植物以枯残根系进入土体上部,腐殖质自表层向下不会锐减,所以灌丛在该范围内有机质下降趋势比较平缓,说明该植被群落下受有机质影响的土壤抗蚀性比较稳定。针叶林表层容重值小,说明土壤土质疏松,有机质通过淋溶等作用向下淀积,加上10—20cm土层原有的有机质,所以在10—20cm土层中针叶林有机质含量比阔叶林和针叶林高。

图2 不同植被群落的土壤抗蚀指数 Fig.2 Soil erosion resistance inhex in different vegetation communities 超出图例符号的垂直棒代表平均值的标准误差

图3 不同植被群落下土壤容重Fig.3 Soil bulk density under different vegetation communities

图4 不同植被群落下土壤有机质含量 Fig.4 Soil organic matter content in different vegetation communities

2.3 土壤团聚体

2.3.1干筛团聚体含量

由表2可知,三种植被类型下的土壤团聚体结构组成多以>3mm以上为主。在0—10cm土层内,阔叶林土壤中>3mm的土壤团聚体含量较针叶林和灌丛高29.96%和9.75%,灌丛土壤下>3mm粒级的土壤团聚体含量比纯针叶林高出5.41%,针叶林土壤下的>3mm的土壤团聚体含量所占比例最低。在10—20cm土层内,灌丛土壤下>3mm粒级的土壤团聚体含量明显高于其他两种植被类型,但针叶林土壤中>3mm的土壤团聚体含量略高于阔叶林。>0.25mm干筛团聚体含量的排序为：在土壤表层内,阔叶林(94.93%)>灌丛(94.75%)>针叶林(89.34%)；在10—20cm土层内,灌丛(94.04%)>针叶林(94.01%)>阔叶林(92.42%)。这些在一定程度上反映出阔叶林对改善表层土壤结构的优势,对增强喀斯特地区的土壤抗蚀性十分有用,另一方面也反映出喀斯特地区的纯针叶林表层土壤团聚体容易受外界条件影响而分散,阔叶林10—20cm土层内土壤团聚性能较差,表明喀斯特地区生态环境的脆弱性,以及保护植被的重要性。

2.3.2水稳性团聚体含量和结构破坏率

由表3可知,三种植被类型土壤中>0.25mm水稳性团聚体的含量为：0—10cm土层内,阔叶林(70.06%)>灌丛(65.96%)>针叶林(60.83%)；10—20cm土层内,灌丛(66.76%)>针叶林(64.25%)>阔叶林(10.21%)。阔叶林表层枯枝落叶物多,腐殖质的含量高,植被类型多样,生物活动丰富,有机质含量丰富[25],所以阔叶林表层土壤的水稳性团聚体含量占优势。相比之下,针叶林内植被类型单一,且上层土壤疏松,有机质等营养物质通过淋溶作用向下淀积,所以针叶林土壤下的水稳性团聚体在10—20cm土层比较占优势。灌丛下土壤的有机质虽然不多,但其本身黏性很强,所以该种植被类型下的土壤水稳性团聚体含量并不低。通常情况下,用土壤水稳性团聚体含量和有机质含量作为评价土壤的抗蚀性指标时,两者之间的变化规律应基本一致[28]。但是土壤团聚体的形成是一个非常复杂的过程,受物理、化学、生物作用的共同影响[29],胶结物质在团聚体形成的过程中也起着非常重要的作用,但是胶结物质可以是无机的、有机的或有机无机结合的,无机胶结物质主要受成土母质影响,有机胶结物质的形成主要受微生物数量、活性及其代谢产物和植物根系分泌物,有机质的输入等因素控制[30],所以该样地内的土壤水稳性团聚体变化与土壤有机质变化不太一致很大程度上是受这些因素影响。结构破坏率的大小表示土壤抗蚀能力的强弱,此样地中结构破坏率的变化规律与>0.25mm水稳性团聚体含量的变化规律一致。

表2 不同植被群落下土壤团聚体(干筛)组成/%

表3 不同植被群落下水稳性团聚体组成(/%)和结构破坏率(/%)

2.4 土壤机械组成及微团聚体对抗蚀性的影响

2.4.1土壤机械组成分析

由表4可知,样地内土壤机械组成差异比较大,0—10cm土层内：阔叶林下土壤以物理性粘粒(<0.01mm)和中细粉粒(0.01—0.001mm)为主,其含量分别占52.41%和31.01%；针叶林以沙粒(1—0.05mm)和物理性黏粒(<0.01mm)为主,其含量分别占29.41%和48.61%；灌丛以沙粒(1—0.05mm)和物理性黏粒(<0.01mm)为主,其含量分别占47.41%和43.61%。在10—20 cm土层内：针叶林和灌丛土壤的机械组成分布特征与表层一致,但是阔叶林土壤以沙粒(1—0.05mm)和物理性黏粒(<0.01mm)为主。表层土壤内针叶林和灌丛土壤的沙粒比重高,主要是由于土体中夹着很多岩石碎屑,这会导致土壤蓄水蓄肥性不好,所以在这三种植被类型下,阔叶林表层的蓄水保肥能力是最好的,自然植被生长优势较大,其土壤抗蚀性能也是最强的。

表4 不同植被群落下土壤机械组成状况

2.4.2以微团聚体含量为基础的土壤抗蚀性指标

由表5可知,不同植被类型0—10cm土壤的团聚状况排序为：阔叶林(34.09%)>灌丛(24.21%)>针叶林(14.94%)；10—20cm土壤的团聚状况排序为：灌丛(39.61%)>针叶林(27.50%)>阔叶林(22.09%)。团聚度的变化特征与团聚状况一致。从团聚状况可以看出阔叶林土壤表层结构明显优于另外两种,团聚度较针叶林和灌丛土壤表层高24.39%和23.61%,但下层结构差,一旦表层遭到破坏,对下层的影响很大。如表5所示,0—10cm土层内,阔叶林土壤的分散率低于针叶林和灌丛土壤分散率,但灌丛土壤分散率较针叶林下土壤分散率低23.67%,表明在阔叶林下表层土具有较强的抗蚀性,灌草地次之,针叶林最弱；在10—20cm土层内,土壤抗蚀性强弱为：灌丛(24.71%)>针叶林(58.84%)>阔叶林(68.17%)。

表5 不同植被群落下土壤团聚状况、团聚度、分散率、分散系数

Table 5 Soil agglomeration status, soil agglomeration degree, soil dispersion rate and soil dispersion coefficient in different vegetation communities

植被类型Vegetation types土壤层次Soil layers/cm团聚状况Soil agglomeration status%团聚度Soil agglomeration degree/%分散率Soil dispersion rate /%分散系数Soil dispersion coefficient /%阔叶林0—1034.09 58.10 54.80 51.28 Broad-leaved forest10—2022.09 41.94 68.17 90.91 针叶林0—1014.94 33.71 78.84 55.56 Coniferous forest10—2027.50 45.32 58.84 79.37 灌丛Bushwood0—1024.21 34.49 55.17 55.56 10—2039.61 45.53 24.71 64.52

2.5 抗蚀指标间的相关性分析

对土壤有机质含量、团聚体含量、水稳性团聚体含量、结构破坏率、团聚状况、团聚度、分散率、分散系数等8个因子进行主成分分析发现：有机质含量在土壤表层和10—20cm土层的贡献率最大,分别为75.073%和86.024%,因此用团聚体含量和有机质含量之间的相关性分析对各抗蚀指标间的关系进行检验。由表6可以得出,在0—10cm土层范围内,水稳性团聚体较干筛团聚体占优势,有机质含量与水稳性团聚体呈显著正相关关系,相关系数为0.977,土壤中水稳性团聚体含量会随着有机质含量的增加而增加, 水稳性团聚体在水中浸湿后稳定性比较高,不容易解体, 在喀斯特区,降雨量大且土层薄,土壤易遭水流侵蚀, 通过改善植被状况来增加土壤有机质含量, 从而增加土壤水稳性团聚体含量来达到提高土壤抗蚀性的目的[1]。结构破坏率、分散系数和有机质含量、水稳性团聚体的相关系数分别为-0.994、-0.931和-0.947、-0.832,呈负相关关系；团聚状况和团聚度和有机质含量、水稳性团聚体的相关系数分别为0.991、0.941和0.997、0.847,呈正相关关系,表明结构破坏率、团聚状况亦可用于评价该区表土土壤抗蚀性能。

表6 0—10cm土层土壤抗蚀指标的相关性分析

**. 在 0.01 水平上显著相关；*. 在 0.05 水平上显著相关

在10—20cm土层范围内,干筛团聚体较水稳性团聚体占优势,降水时靠土壤的下渗能力可以减少水蚀力。在10—20cm土层中,大团聚体越多,土壤间的空隙越大越有利于渗水。由表7可知,有机质含量与干筛团聚体的相关系数呈正相关关系,系数为0.94。在10—20cm土层内,团聚度来表征土壤抗蚀性能比较明显,其与有机质和干筛团聚体含量的相关系数分别为0.927和0.999,呈正相关关系。

表7 10—20cm土层土壤抗蚀指标的相关性分析

*. 在 0.05 水平上显著相关

3 讨论

3.1 不同区域林下土壤抗蚀性的差异

土壤抗蚀性影响因素较多,不同区域在地形、岩性、气候、植被等条件均有差异,导致不同区域林下土壤抗蚀性可能存在差别。由于土壤抗蚀性这一概念在国内应用较多,国外应用较少,因此国外的研究中很少涉及土壤抗蚀性,国内开展的不同植被类型土壤抗蚀性的研究也很少出现在国外期刊上。Zheng等人的研究表明,杉木林和桉树林下土壤的抗蚀性高于茶园、枇杷园和撂荒地[31]。刘旦旦等在黄土残塬沟壑区的研究结果表明,林地土壤抗蚀性最高,荒草地次之,农耕地土壤抗蚀性最差[32]。黄土丘陵沟壑区的研究结果表明,撂荒后植被从一年生草本群落阶段演替到多年生蒿禾类草本群落阶段的过程中,土壤的级配状况在不断改良,颗粒分形维数增大,团粒结构有所改善,土壤团粒结构的分形维数与结构体破坏率均不断减小,稳定性不断增强,从而提高了土壤的抗蚀性[33]。上述研究结果都表明,林地下的土壤抗蚀性高于灌木、草地,或者高级植被群落下土壤的抗蚀性高于低级群落。本研究中将土壤抗蚀性分为两层分析,考虑到土壤侵蚀首先发生在表层,因此将本研究中0—10cm层的研究结果与上述研究结果进行比较。结果表明,土壤表层的抗蚀性为阔叶林>灌丛>针叶林,除针叶林外,本研究的结果与上述其他研究结果是基本一致的。

李阳芳等在云南元阳梯田区的研究结果表明,荒草地的抗蚀性最好,灌木林地次之,有林地,坡耕地最差；>3mm的土壤团聚体含量、3—2mm的土壤结构水稳度、<0.05mm的土壤分散率、土壤的团聚状况、植物根系含量5个指标可以较全面地反映出元阳梯田区土壤的抗蚀性能[34]。吴胡强等在安徽省岳西县的研究表明,土壤抗蚀性强弱表现为马尾松林大于桑树林[35]。南方红壤区的研究表明,土壤抗蚀性强弱的排序为灌木林>松树林>材用竹林>经济林>茶园,结构性颗粒指数、土壤团聚度和土壤结构体破坏率等指标能够很好地衡量土壤抗蚀性强弱[36]。上述研究结果表明,排除人类直接干扰的用地类型外,总体上荒草地土壤抗蚀性最好,灌木林地次之,林地抗蚀性最差。这些结果与前述的研究结果以及与本研究的结果是相反的。这进一步说明土壤抗蚀性的影响因素较多,植被类型并不是决定土壤抗蚀性强弱的唯一因素,需要综合考虑成土母质、气候、物种等因素。

3.2 喀斯特地区林下土壤抗蚀性的差异

如前所述,土壤抗蚀性的影响因素众多,即使是同为喀斯特地区,不同地区、不同自然条件下开展的研究结果仍可能存在一定差异。陈佳等在桂西北喀斯特地区的研究表明,土壤抗蚀性强弱顺序为:原生林>次生林>撂荒地>坡耕地>人工林[21]。蒋丽伟等在花江峡谷示范区的研究表明,土壤抗蚀性在退耕香椿林>退耕撂荒地>退耕油桐林>退耕车桑子地>退耕花椒地>耕地[20]。李会等在普定陈旗流域的研究结果表明,不同土地利用方式土壤的抗蚀性大小顺序为,灌草地>稀疏灌丛地>火烧迹地>复合植被>幼林地>坡耕地[23]。王佩将等在贵州省大方县桶井石漠化综合治理示范区的研究结果表明,土壤抗蚀性大小排序为灌草丛> 乔灌过渡林> 灌木林> 乔木疏林> 人工林> 草坡> 坡耕地。喀斯特山区,坡耕地土壤抗蚀性能最差; 植被恢复过程中,土壤的抗蚀性能先逐渐变好,后逐渐变差,转折点在灌草阶段[22]。上述研究虽然同在喀斯特地区开展,但是研究结果有的表明植被演替阶段越高级,土壤抗蚀性越强；有的研究结果表明,土壤抗蚀性强弱并非与植被演替阶段同向发展,而是存在转折。

总体而言,本研究的结果表明在上层土壤中阔叶林的土壤抗蚀性大于灌木林,而针叶林的土壤抗蚀性最低。上述研究表明,喀斯特地区的土壤抗蚀性存在明显的空间异质性,受自然条件变化的影响很大。中国西南喀斯特地区覆盖范围较广,降水、气温等气候因素区域差异大,会直接或间接地影响土壤属性,从而影响土壤抗蚀性；岩性条件的差异也是造成喀斯特地区土壤抗蚀性差异的重要原因。喀斯特地区主要岩石类型有石灰岩、白云岩,而不同地层的石灰岩、白云岩其纯度、物质组成、夹非碳酸盐岩的状况都有差异,这将直接影响成土母质,从而影响土壤抗蚀性。此外,中国西南喀斯特地区地表破碎、土层薄,土壤属性在微观尺度上也存在显著差异,加上不同研究树种、树龄选择的差异,都可能造成同样是在中国西南喀斯特地区开展的研究,但不同植被类型下土壤可蚀性的差异较大。上述情况进一步说明,为弄清喀斯特地区不同植被下土壤可蚀性的特征,已有研究仍然显得不足,需要在不同的地貌、岩性、气候条件、植被类型区域开展更多的相关研究。

4 结论

通过对实验样地中三种植被群落下0—10cm和10—20cm土层土壤的抗蚀性指标进行评价,发现除有机质指标外,其余指标的变化特征均一致。土壤0—10cm层的抗蚀性：阔叶林>灌丛>针叶林；土壤10—20cm层的抗蚀性：灌丛>针叶林>阔叶林。通过综合两个土层土壤抗蚀性指标的情况,可认为有机质含量、水稳性团聚体、团聚度能较好地表征喀斯特林地土壤抗蚀性强弱。鉴于阔叶林的表层土壤抗蚀性高于灌丛和针叶林,而表层土壤是最先受到土壤侵蚀影响的层次,在喀斯特地区的植被恢复、生态重建、水土保持等工作中,如果以防治土壤侵蚀为目的,应考虑种植阔叶树种,避免针叶树种。此外,相关研究结果对于指导水土保持实践的区域差异与局限性也应予以重视。