论文部分内容阅读
摘要：为探寻较为合理且高效的石漠化植被修复的治理模式，为今后选择治理石漠化的植被组合提供参考依据。在贵州省黔东南苗族侗族自治州施秉县设置样地，选择3种治理模式：模式一为针叶林（马尾松、柏木），模式二为阔叶林（壳斗类，杂灌），模式三为灌木林（荚蒾、杂灌），分析不同治理模式下的土壤养分状况、土壤微生物群落结构、土壤微生物与环境因子相关性。结果表明，3种模式下土壤中微生物数量为细菌>真菌>放线菌;模式二治理下的土壤呈弱碱性，土壤养分含量高于其他2种模式，土壤微生物具有更高的多样性，罗尔斯通氏菌菌属（H16）、Gaiella、红游动菌属（Rhodoplanes）、芽单胞菌属（Gemmatimonas）、Haliangium、子囊菌门（Ascomycoat）、担子菌门（Basidiomycete）、接合菌门（Zygomycota）是其优势菌属，且与土壤中大多数营养物质具有显著正相关性;石漠化治理效果表现为模式二>模式一>模式三，所以模式二的石漠化治理能力更好。
关键词：石漠化;植被修复;治理模式;土壤养分;土壤微生物群落结构;土壤微生物多样性
中图分类号：S154.3 文献标志码： A
文章编号：1002-1302（2021）12-0201-08
收稿日期：2020-08-25
基金项目：贵州黎平石漠化生态系统国家定位观测研究站运行补助专项经费。
作者简介：黄 磊（1987—），男，四川凉山人，硕士，主要从事森林培育研究。E-mail：[email protected]。
通信作者：曾亚军，硕士，副研究员，主要从事森林培育理论与技术方面的研究。E-mail：[email protected]。
石漠化也称为石质荒漠化，指由于水土流失造成的地表植被减少和土壤损失、基岩裸露，最终导致土地丧失农业利用价值和生态环境退化的现象[1]。我国西南地区石漠化面积达13 000 km2，是世界上最大的喀斯特地貌分布区之一[2-3]。石漠化已成为严重阻碍我国西南地区经济增长的主要生态灾害[4]，石漠化地区的生态环境更加脆弱，人均可耕地减少;由于在该条件下植被生长困难，从而导致土壤储水功能降低，造成人畜饮水困难等问题[5]。全国范围内石漠化最为严重的是贵州省，其次是云南和广西地区，资料显示，目前贵州省还有396.33万hm2土地适宜进行生态修复，任务十分艰巨[6]。贵州省石漠化严重阻碍了当地生态环境及可持续经济发展。因此，石漠化治理在贵州地区一直都是政府关注和环境科研工作者研究的重点。
1990年以来，石漠化地区的生态退化环境一直靠植被修复的方式来治理[7]。目前，关于石漠化生态环境修复的研究主要集中在不同治理模式下植物群落组成及多样性、造林保存率、土壤肥力、土壤酶活力等方面[8-10]。近年来，研究发现土壤微生物对改善土壤质量具有重要作用[11]，土壤微生物多样性对塑造陆地生态系统的功能具有重要作用[12]。一方面，可以通过植被固定土壤来改善石漠化程度，另一方面，该过程中土壤营养物质不断积累，能促进植被生长。此外，土壤微生物对维持土壤活力具有重要意义[13]，作为土壤中物质转化的动力，土壤微生物在土壤的形成及改变土壤肥力等方面有着不可替代的作用[14]。土壤微生物对土壤环境变化十分敏感，研究显示在石漠化植被恢复过程中细菌是土壤微生物中的优势类群[15-17]，因此土壤中细菌多样性与石漠化生态系统恢复过程中土壤营养物质变化密切相关，可作为评价土壤质量变化即石漠化修复效率的指标[18]。
目前关于石漠化治理的研究主要集中在探讨一种治理模式或该治理模式下植被演替过程中土壤中微生物群落结构变化，而在不同的治理模式下土壤中的微生物群落结构及治理效果比较的研究则鲜见报道。因此，本研究拟从3种不同的治理模式入手，分别为模式一：针叶林，马尾松、柏木;模式二：阔叶林，壳斗类、杂灌;模式三：灌木林，荚蒾、杂灌。以上3种模式中的马尾松[19-20]、柏木[21]、壳斗类植物[22-23]、荚蒾[24]、杂灌植物[25]都被研究报道过是石漠化治理的先锋树种。因此，探究以上3种不同治理模式下土壤养分含量及土壤微生物群落结构的异同，旨在找到较为合理且高效的石漠化植被修复的治理模式，为今后石漠化治理工作提供一定的选择依据，提高石漠化治理效果从而达到保护生态环境的目的。
1 材料与方法
1.1 研究地区概况
研究地区位于贵州省黔东南苗族侗族自治州施秉县（108° 06′00″～108° 07′12″ E、27°12′00″～27°10′12″ N），屬亚热带湿润季风气候，年均温度 16 ℃，年均降水 1 220 mm[26]。其喀斯特石漠化面积282.95 km2，是典型的石漠化代表地区[27]。采样地点如下：编号2 （108°14′38″ 3E、26°57′31″ 3N）即治理模式一，海拔850 m，其植被类型为针叶林，主要树种为马尾松、柏木;编号12（108°14′26″ E、26°57′57″ N）即治理模式二，海拔833 m，其植被类型为阔叶林，主要树种为壳斗类、杂灌;编号26（108°15′32″ E 、26°57′46″ N）即治理模式三，海拔887 m，其植被类型为灌木林，主要树种为荚蒾、杂灌。
1.2 采样
采样时间为2018年秋季，每样地设置3个采样点，根据实际情况尽量挖深度为0～20 cm的土壤剖面，每个采样点取2份土壤样品，各500 g，1份用于测定土壤养分含量，1份用于微生物群落结构分析。
1.3 土壤养分含量测定
测定土壤pH值及有机碳、硝态氮、铵态氮、全钾、有机质、全氮、速效钾、有效磷含量。
1.4 DNA提取及高通量测序
宏基因组DNA提取，并测序分析。 1.5 数据分析
利用QIIME （V1.9.1）、Canoco 4.5、MEGA 5.0、Mothur、SPSS 25.0等软件对数据进行分析整理。
2 结果与分析
2.1 3种治理模式下各土壤养分含量
由表1可以看出，除了在模式二治理下土壤呈弱碱性外其余2种治理模式下土壤都呈酸性，这表明不同治理模式会对土壤的酸碱性产生不同的影响。在模式二治理下，土壤的全氮、铵态氮、硝态氮及有机碳含量都较其他2种治理模式高。有效磷含量在3种治理模式下并无明显差异。而全钾含量在模式一治理下的土壤中含量最高，速效钾在模式二下含量最高。
土壤有机质、氮、磷、钾等养分含量与石漠化程度密切相关，综上分析结果表明整体上在第2种治理模式下土壤具有更高的养分，其次是模式一、模式三。所以石漠化治理效果为模式二>模式一>模式三。
2.2 3种治理模式下土壤细菌群落结构
不同治理模式下土壤的营养成分即恢复效果不同，可能是由于在不同的植被治理下其土壤微生物群落结构组成差异造成的。土壤微生物是土壤的重要组成部分，在土壤中动植物残体和有机质及其有害物质的分解、生物化学循环和土壤结构的形成过程起着重要的作用。因此，本文对所采用的3种治理模式下的土壤微生物群落结构进行研究，目的在于通过解析土壤微生物群落结构为不同石漠化治理模式恢复效率，为更好地选择石漠化治理模式提供理论依据。
在石漠化治理模式植被恢复过程中开展土壤微生物区系的数量和组成研究，有助于诊断石漠化地区土壤质量和植物营养的修复状况，土壤微生物菌落统计结果显示细菌>真菌>放线菌数量。
2.2.1 土壤细菌群落结构
由物种相对丰度（图1）可以看出，在3种治理模式下的土壤微生物中优势菌群都是变形菌门（Proteobacteria，46.256%、35.727%、31.570%）、酸杆菌门（Acidobacteria，22.145%、22.988%、32.584%）以及疣微菌门（Verrucomicrobia，8.353%、14.071%、32.584%）。根据所有样品在属水平的物种注释及丰度信息，选取丰度排名前35的属，根据其在每个样品中的丰度信息进行聚类，结果显示（图2）在3种治理模式下其微生物的群落结构明显不同。在模式一植被治理下，根瘤菌属（Rhizobium，0.749%）、溶杆菌属（Lysobacter，1.344%）、马赛菌属（Massilia，0.456%）、土微菌属（Pedomicrobium，0.562%）、寡养单胞菌属（Stenotrophobacter，0.261%）以及分支杆菌属（Mycobacterium，0.367%）等是主要优势菌属。在模式二治理下，罗尔斯通氏菌菌属（H16，3.557%）、Gaiella（0.382%）、红游动菌属（Rhodoplanes，3.613%）、芽单胞菌属（Gemmatimonas，0.612%）及Haliangium（0.982%）是优势菌属。而在模式三治理土壤中根微菌属（Rhizomicrobium，0.936%）、苔藓杆菌属（Bryobacter，0.553%）及酸杆菌属（Acidibacter，0.925%）等属于优势菌属。同时，由于土壤理化性质结果显示模式二治理下土壤肥力较高，因此可推测在模式二的植被治理下其土壤中的优势微生物比其他微生物在石漠化治理过程中具有更重要的作用，为今后选择治理石漠化植被提供参考依据。
2.2.2 土壤细菌多样性
总体上来看，各个样品的覆盖率均大于98%，说明样品的采样率足够获得最丰富的OTU（表2）。选择香农指数（Shannon）、辛普森多样性指数（Simpson）和Chao1指数，对不同样品在97%一致性阈值下的α多样性进行分析。从表2可以看出，3种治理模式下土壤微生物的香农指数分别为9.20、9.29、9.33，辛普森多样性指数均为0.99，Chao1指數分别为3 707.239、3 988.579、4 128.411。表明这几种土壤样品的群落结构差异不显著，具有较高的微生物多样性。
选用Weighted Unifrac距离和Unweighted Unifrac距离这2个指标来衡量样品之间的相异系数。由图3可知，以Weighted Unifrac距离作为标准（上方），内生细菌群落相异系数最大的为模式一与模式三，为0.226。结果表明模式二与模式一、模式三之间内生菌的物种多样性差异低于模式一与模式三之间的内生菌多样性。以Unweighted Unifrac距离作为标准（下方），相异系数最大的为模式一与模式三，为0.303;其次为模式一与模式二，为0.287。结果表明3种治理模式下土壤微生物群落结构具有差异。
2.2.3 3种治理模式下土壤细菌组成与环境因子相关性
微生物菌属与环境因子相关性（图4）结果表明，全钾含量与模式一中优势菌属根瘤菌属（Rhizobium，0.828，P
关键词：石漠化;植被修复;治理模式;土壤养分;土壤微生物群落结构;土壤微生物多样性
中图分类号：S154.3 文献标志码： A
文章编号：1002-1302（2021）12-0201-08
收稿日期：2020-08-25
基金项目：贵州黎平石漠化生态系统国家定位观测研究站运行补助专项经费。
作者简介：黄 磊（1987—），男，四川凉山人，硕士，主要从事森林培育研究。E-mail：[email protected]。
通信作者：曾亚军，硕士，副研究员，主要从事森林培育理论与技术方面的研究。E-mail：[email protected]。
石漠化也称为石质荒漠化，指由于水土流失造成的地表植被减少和土壤损失、基岩裸露，最终导致土地丧失农业利用价值和生态环境退化的现象[1]。我国西南地区石漠化面积达13 000 km2，是世界上最大的喀斯特地貌分布区之一[2-3]。石漠化已成为严重阻碍我国西南地区经济增长的主要生态灾害[4]，石漠化地区的生态环境更加脆弱，人均可耕地减少;由于在该条件下植被生长困难，从而导致土壤储水功能降低，造成人畜饮水困难等问题[5]。全国范围内石漠化最为严重的是贵州省，其次是云南和广西地区，资料显示，目前贵州省还有396.33万hm2土地适宜进行生态修复，任务十分艰巨[6]。贵州省石漠化严重阻碍了当地生态环境及可持续经济发展。因此，石漠化治理在贵州地区一直都是政府关注和环境科研工作者研究的重点。
1990年以来，石漠化地区的生态退化环境一直靠植被修复的方式来治理[7]。目前，关于石漠化生态环境修复的研究主要集中在不同治理模式下植物群落组成及多样性、造林保存率、土壤肥力、土壤酶活力等方面[8-10]。近年来，研究发现土壤微生物对改善土壤质量具有重要作用[11]，土壤微生物多样性对塑造陆地生态系统的功能具有重要作用[12]。一方面，可以通过植被固定土壤来改善石漠化程度，另一方面，该过程中土壤营养物质不断积累，能促进植被生长。此外，土壤微生物对维持土壤活力具有重要意义[13]，作为土壤中物质转化的动力，土壤微生物在土壤的形成及改变土壤肥力等方面有着不可替代的作用[14]。土壤微生物对土壤环境变化十分敏感，研究显示在石漠化植被恢复过程中细菌是土壤微生物中的优势类群[15-17]，因此土壤中细菌多样性与石漠化生态系统恢复过程中土壤营养物质变化密切相关，可作为评价土壤质量变化即石漠化修复效率的指标[18]。
目前关于石漠化治理的研究主要集中在探讨一种治理模式或该治理模式下植被演替过程中土壤中微生物群落结构变化，而在不同的治理模式下土壤中的微生物群落结构及治理效果比较的研究则鲜见报道。因此，本研究拟从3种不同的治理模式入手，分别为模式一：针叶林，马尾松、柏木;模式二：阔叶林，壳斗类、杂灌;模式三：灌木林，荚蒾、杂灌。以上3种模式中的马尾松[19-20]、柏木[21]、壳斗类植物[22-23]、荚蒾[24]、杂灌植物[25]都被研究报道过是石漠化治理的先锋树种。因此，探究以上3种不同治理模式下土壤养分含量及土壤微生物群落结构的异同，旨在找到较为合理且高效的石漠化植被修复的治理模式，为今后石漠化治理工作提供一定的选择依据，提高石漠化治理效果从而达到保护生态环境的目的。
1 材料与方法
1.1 研究地区概况
研究地区位于贵州省黔东南苗族侗族自治州施秉县（108° 06′00″～108° 07′12″ E、27°12′00″～27°10′12″ N），屬亚热带湿润季风气候，年均温度 16 ℃，年均降水 1 220 mm[26]。其喀斯特石漠化面积282.95 km2，是典型的石漠化代表地区[27]。采样地点如下：编号2 （108°14′38″ 3E、26°57′31″ 3N）即治理模式一，海拔850 m，其植被类型为针叶林，主要树种为马尾松、柏木;编号12（108°14′26″ E、26°57′57″ N）即治理模式二，海拔833 m，其植被类型为阔叶林，主要树种为壳斗类、杂灌;编号26（108°15′32″ E 、26°57′46″ N）即治理模式三，海拔887 m，其植被类型为灌木林，主要树种为荚蒾、杂灌。
1.2 采样
采样时间为2018年秋季，每样地设置3个采样点，根据实际情况尽量挖深度为0～20 cm的土壤剖面，每个采样点取2份土壤样品，各500 g，1份用于测定土壤养分含量，1份用于微生物群落结构分析。
1.3 土壤养分含量测定
测定土壤pH值及有机碳、硝态氮、铵态氮、全钾、有机质、全氮、速效钾、有效磷含量。
1.4 DNA提取及高通量测序
宏基因组DNA提取，并测序分析。 1.5 数据分析
利用QIIME （V1.9.1）、Canoco 4.5、MEGA 5.0、Mothur、SPSS 25.0等软件对数据进行分析整理。
2 结果与分析
2.1 3种治理模式下各土壤养分含量
由表1可以看出，除了在模式二治理下土壤呈弱碱性外其余2种治理模式下土壤都呈酸性，这表明不同治理模式会对土壤的酸碱性产生不同的影响。在模式二治理下，土壤的全氮、铵态氮、硝态氮及有机碳含量都较其他2种治理模式高。有效磷含量在3种治理模式下并无明显差异。而全钾含量在模式一治理下的土壤中含量最高，速效钾在模式二下含量最高。
土壤有机质、氮、磷、钾等养分含量与石漠化程度密切相关，综上分析结果表明整体上在第2种治理模式下土壤具有更高的养分，其次是模式一、模式三。所以石漠化治理效果为模式二>模式一>模式三。
2.2 3种治理模式下土壤细菌群落结构
不同治理模式下土壤的营养成分即恢复效果不同，可能是由于在不同的植被治理下其土壤微生物群落结构组成差异造成的。土壤微生物是土壤的重要组成部分，在土壤中动植物残体和有机质及其有害物质的分解、生物化学循环和土壤结构的形成过程起着重要的作用。因此，本文对所采用的3种治理模式下的土壤微生物群落结构进行研究，目的在于通过解析土壤微生物群落结构为不同石漠化治理模式恢复效率，为更好地选择石漠化治理模式提供理论依据。
在石漠化治理模式植被恢复过程中开展土壤微生物区系的数量和组成研究，有助于诊断石漠化地区土壤质量和植物营养的修复状况，土壤微生物菌落统计结果显示细菌>真菌>放线菌数量。
2.2.1 土壤细菌群落结构
由物种相对丰度（图1）可以看出，在3种治理模式下的土壤微生物中优势菌群都是变形菌门（Proteobacteria，46.256%、35.727%、31.570%）、酸杆菌门（Acidobacteria，22.145%、22.988%、32.584%）以及疣微菌门（Verrucomicrobia，8.353%、14.071%、32.584%）。根据所有样品在属水平的物种注释及丰度信息，选取丰度排名前35的属，根据其在每个样品中的丰度信息进行聚类，结果显示（图2）在3种治理模式下其微生物的群落结构明显不同。在模式一植被治理下，根瘤菌属（Rhizobium，0.749%）、溶杆菌属（Lysobacter，1.344%）、马赛菌属（Massilia，0.456%）、土微菌属（Pedomicrobium，0.562%）、寡养单胞菌属（Stenotrophobacter，0.261%）以及分支杆菌属（Mycobacterium，0.367%）等是主要优势菌属。在模式二治理下，罗尔斯通氏菌菌属（H16，3.557%）、Gaiella（0.382%）、红游动菌属（Rhodoplanes，3.613%）、芽单胞菌属（Gemmatimonas，0.612%）及Haliangium（0.982%）是优势菌属。而在模式三治理土壤中根微菌属（Rhizomicrobium，0.936%）、苔藓杆菌属（Bryobacter，0.553%）及酸杆菌属（Acidibacter，0.925%）等属于优势菌属。同时，由于土壤理化性质结果显示模式二治理下土壤肥力较高，因此可推测在模式二的植被治理下其土壤中的优势微生物比其他微生物在石漠化治理过程中具有更重要的作用，为今后选择治理石漠化植被提供参考依据。
2.2.2 土壤细菌多样性
总体上来看，各个样品的覆盖率均大于98%，说明样品的采样率足够获得最丰富的OTU（表2）。选择香农指数（Shannon）、辛普森多样性指数（Simpson）和Chao1指数，对不同样品在97%一致性阈值下的α多样性进行分析。从表2可以看出，3种治理模式下土壤微生物的香农指数分别为9.20、9.29、9.33，辛普森多样性指数均为0.99，Chao1指數分别为3 707.239、3 988.579、4 128.411。表明这几种土壤样品的群落结构差异不显著，具有较高的微生物多样性。
选用Weighted Unifrac距离和Unweighted Unifrac距离这2个指标来衡量样品之间的相异系数。由图3可知，以Weighted Unifrac距离作为标准（上方），内生细菌群落相异系数最大的为模式一与模式三，为0.226。结果表明模式二与模式一、模式三之间内生菌的物种多样性差异低于模式一与模式三之间的内生菌多样性。以Unweighted Unifrac距离作为标准（下方），相异系数最大的为模式一与模式三，为0.303;其次为模式一与模式二，为0.287。结果表明3种治理模式下土壤微生物群落结构具有差异。
2.2.3 3种治理模式下土壤细菌组成与环境因子相关性
微生物菌属与环境因子相关性（图4）结果表明，全钾含量与模式一中优势菌属根瘤菌属（Rhizobium，0.828，P