煤矸石充填复垦土壤细菌群落变化

焦 赫1,2,李新举1,2

(1.山东农业大学 资源与环境学院,山东 泰安 271018; 2.土肥资源高效利用国家工程实验室,山东 泰安 271018)

摘 要:煤炭开采造成大面积矿区土地塌陷,导致当地生态系统的脆弱性显著增加。煤矸石充填复垦是修复矿区塌陷地的常用方法,土壤细菌群落对于反映复垦土壤质量和生态健康状况具有重要的指示作用。为了了解矿区塌陷地煤矸石充填复垦对土壤细菌群落的影响,采用野外采样和高通量lllumina Miseq 16S rRNA测序技术,以周边普通耕地作为对照组(CK),分析了中国东部典型高潜水位矿区采煤塌陷地煤矸石充填复垦耕地(G)不同土层深度(0~20,20~40,40~60 cm)下土壤理化性质、土壤细菌群落多样性和结构组成的变化,以及土壤细菌群落与土壤理化性质的关系。结果表明:① 煤矸石充填复垦显著改变了研究区土壤理化性质,不同土层深度下土壤性质变化趋势差异明显;② 煤矸石充填复垦地的土壤细菌群落多样性与丰富度均低于普通耕地,且与土层深度呈显著负相关(P<0.05),40~60 cm层因煤矸石充填造成细菌群落组成差异大于其他土层。另外,变形菌门、放线菌门和酸杆菌门是所有处理中的优势菌门;③ 典型相关分析(CCA)和相关性分析(Pearson)结果表明,pH、含水量(MoiC)、有机碳(SOC)、速效磷(AP)和速效钾(AK)是影响土壤细菌群落组成的主要土壤理化性质,优势菌门仅与pH,AP,MoiC等部分理化性质显著相关。

关键词:复垦土壤;煤矸石充填;细菌群落;高通量测序;采煤塌陷区

中图分类号:TD88

文献标志码:A

文章编号:0253-9993(2021)10-3332-10

移动阅读

收稿日期:2020-06-28

修回日期:2020-11-06

责任编辑:黄小雨

DOI:10.13225/j.cnki.jccs.2020.1104

基金项目:国家自然科学基金资助项目(41771324);山东省煤田地质局科研专项奖励基金资助项目(鲁煤地科字(2019)8号)

作者简介:焦 赫(1995—),女,山东潍坊人,硕士研究生。E-mail:986843108@qq.com

通讯作者:李新举(1965—),男,山东济宁人,教授。Tel:0538-8246203,E-mail:lxj0911@126.com

引用格式:焦赫,李新举.煤矸石充填复垦土壤细菌群落变化[J].煤炭学报,2021,46(10):3332-3341.

JIAO He,LI Xinju.Variation in the soil bacterial community of reclaimed land filled with coal gangue[J].Journal of China Coal Society,2021,46(10):3332-3341.

Variation in the soil bacterial community of reclaimed land filled with coal gangue

JIAO He1,2,LI Xinju1,2

(1.College of Resources and Environment,Shandong Agricultural University,Taian 271018,China;2.National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources,Taian 271018,China)

Abstract:Coal mining often causes large areas of subsidence land collapse in mining areas, leading to a significant increase in the vulnerability of the local ecosystem. Reclamation with coal gangue filling is commonly adopted to repair the subsidence of mining areas, and soil bacterial community plays an important indicative role in reflecting the quality and ecological health of reclaimed soil. In order to understand the impact of reclamation of subsidence area with coal gangue filling on soil bacterial community, takeing surrounding normal cultivated land (CK) as controla comparison, field sampling and high throughput Illumina Miseq 16S rRNA gene sequencing were used to analyze the soil physical and chemical properties, the variations of soil bacterial community composition and diversity, and the relationship between them in different soil depths (0-20, 20-40 and 40-60 cm) of reclaimed land filled with coal gangue (G). The results showed that:① Ccoal gangue reclamation significantly changed the soil physical and chemical properties in the study area, the changing trend of soil properties under different soil depths was obviously different.② The diversity and abundance of soil bacterial community in the reclaimed land filled with coal gangue were lower than those in normal cultivated lands, which also showed a significant negative correlation with soil depth (P<0.05). The composition difference of bacterial community in the 40-60 cm soil layer due to coal gangue filling was greater than that in other soil layers. Besides, Proteobacteria, Actinobacteria, and Acidobacteria were the dominant phyla in all treatments.③ Results from canonical correlation analysis (CCA) and Pearson correlation analysis indicated that pH, Moisture Content (MoiC), Soil Organic Carbon (SOC), Available Phosphorus (AP), and Available Potassium (AK) were the main soil physical and chemical properties in shaping soil bacterial communities, and the dominant bacteria phyla were only significantly associated with some soil properties such as pH, AP and MoiC.

Key words:reclaimed soil;coal gangue filling;bacterial community;high-throughput sequencing;coal mining subsidence area

高潜水位煤矿区煤炭开采常造成严重的地面沉降,多年来一直采用回填复垦技术来减轻沉降影响[1-2]。煤矸石因为价格低廉且长期堆放会对地表结构和生态环境造成严重破坏,而被广泛用作矿区沉陷土地充填材料[3-4]。然而煤矸石作为固体废弃物,养分含量低,且含有较多重金属元素,其充填后的耕地质量较低,易造成严重的土壤、水和空气污染等问题[5]。因此,煤矸石充填复垦地土壤质量的研究是国内外相关专家关注的焦点[6-7]。目前,对采煤塌陷地煤矸石充填复垦地的土壤质量研究主要集中于植物保护、重金属污染、土壤酶活性等方面[8-10],在煤矸石充填复垦对土壤微生物生态影响的研究较缺乏。

土壤微生物参与了土壤生态系统中的物质循环、信息传递以及能量转换,在养分的转化循环、维持生态系统稳定性和植被演替恢复中起着决定性作用[11-14]。细菌作为土壤地下和地上系统交流的重要媒介,参与了土壤生态系统中大多数养分的循环转化,是推动植物凋落物和土壤有机质分解的重要驱动力[15]。其中,土壤细菌群落被认为是生物地球化学循环、气候调节或植物生产力等土壤功能的关键参与者[16],对土壤物理化学性质、植被类型和水文变化等环境因素变化敏感[17-18]。因此,土壤细菌群落组成和结构变化常被用来反映土壤环境质量的变化。

研究表明,土壤细菌群落组成和结构的变化是退化土壤环境恢复效果的重要反映,包括湿地生态恢复、退化森林治理、盐碱地修复[19-21]。此外,土壤细菌群落在衡量矿区复垦土壤质量和生态效益方面具有重要作用。王瑾[22]研究发现土壤细菌群落对不同复垦时间、不同破坏程度、不同复垦区域的采煤沉陷复垦土壤响应不尽相同。KNELLER等[23]通过研究采矿区不同复垦技术下土壤细菌等微生物对于土壤外源有机质添加后的反映,发现细菌群落在评估矿山土壤质量恢复效果的过程中起着关键指示作用。郭晓明等[24]研究矿区沉陷复垦地土壤细菌数量及酶活性空间分布特征发现,与非沉陷地相比,沉陷区对土壤细菌群落的生长存在一定抑制作用,采煤沉陷导致耕地土壤细菌数量空间异质性显著增加。李兆龙等[25]在铜尾矿应用土壤细菌解释了土壤微生物群落结构对重金属污染的响应,证明细菌等微生物对矿区退化土壤进行生态修复评估的重要价值。LI等[26]研究了黄土高原西北部露天煤矿复垦区土壤微生物丰度和多样性,说明细菌等微生物群落演替对于不同复垦方案的响应敏感程度存在显著差异。由此可见,土壤细菌及其群落在矿区采煤沉陷复垦领域的传统研究较为广泛。然而在高潜水位矿区,煤矸石充填复垦作为主要的复垦方式之一,目前利用现代高通量测序方法研究煤矸石充填后土壤细菌群落的变化是有限的。

笔者利用高通量lllumina Miseq 16S rRNA测序方法,探讨了中国东部典型高潜水位矿区采煤塌陷地煤矸石充填复垦耕地不同土层深度下土壤细菌群落多样性和结构组成的变化,并分析了特定优势细菌门与土壤理化性质的关系。主要目的是确定:① 煤矸石充填复垦后土壤性质变化;② 煤矸石充填复垦后土壤细菌群落多样性和结构组成的差异变化;③ 可能影响细菌群落结构和多样性的潜在关键土壤特性。

1 材料与方法

1.1 研究区概况

研究区位于山东济宁市东滩矿区,地理坐标为35.41°N~116.96°E,总面积约为1 616.56 km2(图1)。地形以平原为主,作物种类丰富,气候为四季分明的温带大陆季风性气候,雨量集中,水热同步。根据美国农业部土壤质地划分和美国农业部土壤系统分类的数据[27],研究区域的土壤类型主要为湿润淋溶土。土壤质地为砂壤土,砾石含量约7.3%。

图1 研究区位置
Fig.1 Location of study area

东滩矿区因常年煤炭开采造成了大面积采煤塌陷地,自2001年开始进行塌陷地的充填复垦。该区充填复垦方式根据充填物的不同主要有煤矸石充填、粉煤灰充填和客土充填等方式。研究区主要采用的典型复垦工艺是煤矸石充填复垦。具体做法是在煤矿区的土地塌陷平稳后,前期先进行表土剥离,然后利用煤矸石进行充填,最后再回填表土层的土壤。复垦规划设计标高为47.0 m,未复垦土地平均标高约为39.8 m,平均充填矸石厚度约为6.2 m,矸石充填后覆土厚度在60~100 cm[28]。充填煤矸石来源于研究矿区煤矸石山和部分选煤厂矸石,煤矸石以碳酸盐类和碳质砂岩为主,其风化物以砾石为主,颗粒间差异较大。由于早年复垦过程中多次受到机械压实,土壤较为紧实。复垦后的利用方式主要为耕地、种植小麦、玉米等作物。

1.2 样品采集

2018年8月,按照室内设计路线,在试验区设置3个16 m×16 m的样地(3个重复),在每个样地内,采用直径5 cm的土壤螺旋钻,用对角线取样法和环刀法分3层(0~20,20~40,40~60 cm)于5个点均匀采样,并以周边正常耕地作为对照。每层土样采集5份并均匀混合,质量约2 kg,且每层取1 kg环刀样,用于土壤含水量和容重的测定,所有样品均装入无菌保鲜袋放置铺满冰袋的保温箱中,立即运回实验室。除去环刀样,将其他土样中明显的植物残渣、石子砖块去除后,按照四分法,将每层所取5份土样均匀混合,每个样品的一半经风干研磨并过2 mm的筛子,储存在室温下用于化学和物理性质分析;另一半立即放置-80 ℃下冷冻用于微生物分析。

1.3 土样理化性质

把土样置于铝制铁盒中,在105 ℃下烘干48 h,进行土样含水量(MoiC)测量。容重(BD)采用环刀法测定。采用高温外加热重铬酸钾氧化-容重法测定土壤有机碳(SOC)含量。用pH计(PHSJ-3F)和TDR(TRIME-IPH)测定土壤中的pH和电导率(EC)。分别用凯氏法测定碱解氮(AN)、奥尔森法测定速效磷(AP)、用盐酸铵萃取后用原子吸收光谱法测定速效钾(AK)。对上述指标均进行了3次重复测试,并根据相关性做进一步分析。

1.4 土壤细菌16S rRNA基因序列分析

1.4.1 DNA抽提和PCR扩增

根据制造商说明,使用E.Z.N.A.® soil试剂盒(Omega Bio-tek,Norcross,GA,U.S.)说明书从每种土壤样品中称取0.5 g进行总DNA抽提。DNA浓度和纯度利用NanoDrop 2000进行检测,利用1%的琼脂糖凝胶电泳检测DNA提取质量;用引物515F(5′-GTGYCAGCMGCCGCGGTAA-3′)和806R(5′-GGACTACNVGGGTWTCTAAT-3′)对16S rRNA基因V5-V6可变区进行PCR扩增。扩增程序为:94 ℃预变性10 min,35个循环(94 ℃,变性40 s,55 ℃ 退火30 s,72 ℃延伸90 s),最后72 ℃延伸10 min (PCR采用TransGen AP221-02:TransStart Fastpfu DNA Polymerase;PCR仪:ABI GeneAmp® 9700型)。扩增体系为24 μL,18.35 μL 水,2.5 μL 10×buffer,2 μL 2.5 mmol/L dNTPs,各0.5 μL引物(515F和806R),0.15 μL Taq聚合酶,1 μL所提取DNA。

1.4.2 高通量测序

用2%琼脂糖凝胶电泳对3个重复PCR产物进行检测并汇集,然后利用AxyPrep DNA Gel Extraction Kit(Axygen Biosciences,Union City,CA,USA)进行纯化,并利用QuantiFluorTM-ST(Promega,USA)进行检测定量。最后,根据Illumina MiSeq平台(Illumina,San Diego,USA)标准操作规程将纯化后的扩增片段构建PE 2×300的文库。通过高通量测序对所有纯化和混合样本测序读数进行分析。最终,实验数据利用Illumina公司的Miseq PE300平台进行测序(上海美吉生物医药科技有限公司)。

1.5 数据处理和分析

采用SPSS 22.0软件对土壤性质数据进行单因素方差分析(ANOVA),在P<0.05时观察到显著性。使用R vegan软件包计算了α多样性指数,包括丰富度(Ace和Chao指数)、微生物群落多样性(Shannon指数)。用β多样性分析评价物种复杂度的差异,利用R vegan软件包做基于Bray-Curtis差异矩阵的主坐标分析(PCoA)以比较和可视化土壤样本之间的相似性。为探讨细菌群落与土壤性质之间的关系,利用CANOCO 5.0对门水平下前10个优势菌门与土壤理化性质进行典型相关分析(CCA),以Pearson相关分析法辅助分析。

2 结果与分析

2.1 土壤理化性质

所考虑的土壤理化性质(AN,AP,AK,SOC,pH,EC,BD,MoiC)在不同处理间差异显著(表1)。G的pH值分布在8.23~8.68,3层pH均显著高于CK(P<0.05),且在深度变化上呈现逐渐降低的变化趋势。G的EC在不同土层深度下均大于CK,随土层深度增加而增大。G和CK的EC在DL(Deep Layer)均高于SL(Surface Layer)和ML(Middle Layer)层。G的AN和AK变化趋势一致,ML和DL土层质量分数(均值为72.85 mg/kg和75.82 mg/kg)远低于表层(120.67 mg/kg 和108.83 mg/kg),G的AP质量分数随土层深度增加呈现下降趋势。与CK相比,G的AN质量分数(120.67 mg/kg)在SL层略高于CK(120.57 mg/kg)。AK,AP在不同土层深度下质量分数均明显低于CK。G的SOC质量分数随土层深度增加呈现逐渐下降趋势,在SL层与CK质量分数接近,但其他层远低于CK。G的BD在各层均高于CK,且ML质量分数高于其他土层。整体上G的含水量低于CK。

表1 复垦地和对照地的土壤性质

Table 1 Soil properties of the reclaimed and control sites

样品编号深度/mAN/(mg·kg-1)AP/(mg·kg-1)AK/(mg·kg-1)SOC/(g·kg-1)pHEC/(μS·cm-1)BD/(g·cm-3)MoiC/%G-SL0~20120.67±13.32a8.03±2.32b108.83±4.54a7.64±1.18a8.68±0.24a368.00±29.31a1.65±0.09a0.18±0.97aG-ML20~4069.63±5.09b2.71±0.26b72.53±1.50b4.93±0.36bc8.60±0.14ab378.33±33.02a1.76±0.08a0.13±1.04cG-DL40~6076.07±4.26ab2.60±0.27b79.11±4.48b4.07±0.75c8.23±0.30bc403.33±46.52a1.63±0.06a0.10±0.67dCK-SL0~20120.57±38.17a15.10±5.89a113.67±19.22a6.37±1.14ab7.98±0.24c361.33±43.25a1.40±0.04b0.18±0.97aCK-ML20~4092.40±17.50ab5.30±2.46b90.33±7.57b5.66±0.62bc8.07±0.09c358.67±35.53a1.39±0.06b0.18±0.20abCK-DL40~60114.53±42.58ab4.43±2.15b85.67±10.69b6.04±1.49ab7.51±0.21d374.00±35.73a1.41±0.11b0.17±0.50b

注:不同上标字母表示处理间差异显著(P<0.05),下同。

2.2 土壤细菌群落组成特征

从18个样品整合的数据集中保留了1 607 790个高质量的序列(细菌群落平均每个样本89 322个),共检测到11 430个OTUs,有效序列长度均>200 bp。G和CK的ACE,Chao和Shannon指数在丰富度和多样性方面存在显著差异(表2)。所有样本的覆盖率都在96%以上,说明测序读数足以进行分析。通过对Chao,ACE和Shannon指数在不同土层深度下的计算,结果表明,在SL层,G的丰富度和多样性略高于CK,但在DL层,G的土壤细菌ACE指数、Chao指数和Shannon指数的最低值分别为1 062,1 090 和4.38,远低于CK。 ANOVA分析表明,与CK相比,Ace,Chao,Shannon指数在不同土层深度差异显著(P<0.05),在DL层,G的Ace,Chao,Shannon指数分别较SL层降低74.44%~80.92%,73.81%~79.90% 和29.93%~41.00%。因此,说明整体上煤矸石充填后的耕地的细菌群落多样性和丰富度低于CK,且与土层深度呈显著负相关。PCoA发现,G和CK在不同土层深度下土壤细菌群落组成存在显著差异(图2)。

表2 不同处理下观测到的操作分类单元、覆盖率、丰富度和多样性的估计值

Table 2 Estimated number of observed OTUs,coverage,richness and diversity under different treatments

相品编号深度/mOTUsACEChaoShannonCoverageG-SL0~203 434±171.55a4 578±163.47ab4 536±222.36ab6.73±0.12a96.98±0.28bcG-ML20~403 088±64.01b4 069±621.26bc4 068±591.53bc6.47±0.13ab97.63±0.89bG-DL40~60935±102.51d1 062±148.66d1 090±128.29d4.38±0.35d99.56±0.21aCK-SL0~202 849±128.49c3 841±320.31c3 872±281.65c6.29±0.11c97.60±0.30bCK-ML20~403 028±135.90ab4 249±153.18ab4 242±119.69ab6.23±0.09c97.23±0.09bcCK-DL40~603 360±99.57a4 817±127.79a4 789±174.12a6.48±0.10ab96.81±0.09c

图2 主坐标分析
Fig.2 Principal coordinates analysis

在门的水平上,样品的总体细菌群落组成相似,但观察到一些样品的不同比例,如图3所示,图3中,不同小写字母为处理间差异显著(P<0.05)。在所有样品中,变形菌门是最丰富的种类,平均相对丰度为32.43%。放线菌门是第2丰富的门,平均相对丰度为17.60%。酸杆菌门的平均相对丰度为14.46%。其他优势菌门平均相对丰度值按降序排列依次为:绿弯菌门(5.90%)、拟杆菌门(4.75%)、芽单胞菌门(4.53%)、厚壁菌门(3.39%)、浮霉菌门(2.77%)、罗卡氏菌门(2.68%)、疣微菌门(1.49%)、球状细菌(1.30%)。与CK相比,G的变形菌门和酸杆菌门平均相对丰度值随土层深度增加呈下降趋势,变形菌门相对丰度下降了14.13%~24.25%,酸杆菌门相对丰度下降了91.16%~97.41%,而放线菌门呈明显上升趋势,相对丰度上升了113.85%~167.17%。优势菌门相对丰度在CK各土层中变化没有明显规律。另外,与自身相比,G的优势菌门在DL层相对丰度远低于其他土层。

图3 细菌门的相对丰度
Fig.3 Relative abundance of bacterial phyla

2.3 土壤理化性质与土壤细菌群落的相关性

采用CCA方法分别研究了G和CK土壤细菌群落(前10个优势门的相对丰度)与土壤理化性质(AN,AP,AK,SOC,pH,EC,BD,MoiC)之间的关系(图4)。研究发现,该区土壤细菌群落的总体结构与pH,MoiC,SOC,AP,AK显著相关。另外,笔者还计算了G和CK前10个优势菌门与土壤理化性质之间的皮尔逊系数(Pearson)见表3,4。研究发现,相对丰度较高的优势菌群(变形菌门、放线菌门和酸杆菌门)仅与土壤某些理化性质显著相关。在G中,MoiC与变形菌门、酸杆菌门、绿弯菌门等优势菌门相对丰度均呈显著正相关(P<0.05)。pH与拟杆菌门的相对丰度呈极显著负相关(P<0.01),与芽单胞菌门呈显著正相关(P<0.05)。在CK中,相对丰度较高的前3类优势菌仅有变形菌门相对丰度与AP呈显著正相关(P<0.05)。此外,放线菌门、酸杆菌门、绿弯菌门、厚壁菌门与大多数土壤理化性质无显著相关。

表3 煤矸石充填复垦地土壤细菌群落与土壤理化性质的相关系数

Table 3 Correlation coefficients between soil bacterial community and soil physical and chemical properties of reclaimed land filled with coal gangue

菌群pHrPANrPAPrPAKrPBDrPSOCrPECrPMoiCrPProteobacteria0.5460.1280.5000.1700.6290.0700.5120.159-0.2290.5530.4480.227-0.0830.8320.672∗0.048Actinobacteria-0.1470.706-0.5250.146-0.4850.186-0.5150.1560.1530.694-0.5610.116-0.0240.951-0.5790.102Acidobacteria-0.5400.1340.6320.0680.5610.1160.5530.1230.2090.5890.755∗0.019-0.4480.2270.749∗0.011Chloroflexi-0.5980.0890.2860.4560.4320.2460.2900.4490.4900.1800.5030.168-0.4480.2260.675∗0.046Bacteroidetes0.853∗∗0.003-0.1720.659-0.3190.402-0.1260.747-0.5200.151-0.3410.3690.5670.111-0.5440.130Gemmatimonadetes-0.718∗0.0290.4030.2820.5600.1170.4010.2850.2380.5370.6470.060-0.4770.1940.680∗0.044Firmicutes-0.0040.991-0.5450.129-0.5310.141-0.5500.1250.2340.545-0.5930.0920.4210.260-0.5580.119Planctomycetes0.5910.0930.825∗∗0.0060.840∗∗0.0050.830∗∗0.0060.0350.9280.876∗∗0.002-0.2070.5930.943∗∗0Rokubacteria0.2230.565-0.3290.387-0.2290.553-0.4020.2840.4490.225-0.1550.690-0.791∗0.011-0.0150.969Verrucomicrobia0.5450.1290.667∗0.0500.761∗0.0170.712∗0.0310.1710.6590.6200.075-0.1940.6180.889∗∗0.001

注:*表示在0.05的水平上呈显著相关;**表示在0.01的水平上呈极显著相关,下同。

图4 典型相关分析
Fig.4 Canonical correspondence analysis

3 讨 论

3.1 土壤理化性质变化

G与CK存在显著差异,不同土层深度下土壤理化性质变化趋势差异很大。G的AN,AP,AK,SOC在SL层含量高于其他土层,这可能与人工施肥、表土回填有关[29-30]。但从ML层开始,特别是AP,SOC含量随土层深度增加呈现明显下降趋势。这是由于煤矸石养分含量低,充填效果差,导致土壤养分下渗速率变低[31-32]。与CK相比,G的EC在不同土层深度下均大于对照地,且随土层深度增加而增大,这也与煤矸石充填有关,由于高地下潜水位从底部煤矸石充填区带来可溶性物质,并逐年随水分蒸发向上层迁移,导致与煤矸石接触层的EC明显高于表层,与前人研究结果[28]类似。研究区复垦地由于早年受机械压实影响严重,较多次数的碾压破坏了土壤结构[33],使得土壤紧实,容重增大,深层次土壤群落恢复速率慢,从而影响土壤养分的充分转化,同时也导致复垦地的MoiC低于普通耕地[34-35]。这些发现表明,充填复垦质量影响土壤养分质量的提高,通过提高煤矸石充填技术和加强土地耕作深翻,将会在稳定土壤生态系统和控制土地退化等方面发挥重要作用。

表4 对照地土壤细菌群落与土壤理化性质的相关系数

Table 4 Correlation coefficients between soil bacterial community and soil physical and chemical properties of contrast land

菌群pHrPANrPAPrPAKrPBDrPSOCrPECrPMoiCrPProteobacteria0.0540.8910.0690.8590.772∗0.0150.6110.080-0.3290.387-0.1560.689-0.2030.601-0.0500.898Actinobacteria-0.0670.863-0.1970.6120.1870.629-0.0150.969-0.5480.127-0.0360.927-0.0310.938-0.3930.295Acidobacteria-0.0640.8710.3710.326-0.4200.261-0.4750.1960.4330.2450.3840.3080.2320.5470.3130.412Chloroflexi-0.1560.688-0.2340.574-0.5910.094-0.3630.3370.2360.5410.3370.3750.2410.532-0.2130.583Bacteroidetes0.0320.9360.2170.6590.946∗∗0.0000.6220.074-0.3520.3530.2750.4730.0610.8750.0870.823Gemmatimonadetes-0.5430.131-0.5280.144-0.5170.155-0.2370.540-0.2150.579-0.3410.369-0.5540.122-0.6400.064Firmicutes0.4820.250-0.2560.505-0.2540.5090.0080.9380.2570.505-0.2630.4940.2440.5260.3060.423Planctomycetes0.1740.6540.4160.2660.5690.1100.4080.2750.0950.8070.5410.1330.3970.2910.5720.108Rokubacteria-0.4050.280-0.0510.897-0.816∗∗0.007-0.707∗0.0330.3440.364-0.3520.3530.5100.160-0.3720.324Verrucomicrobia0.3370.3760.2390.5360.833∗0.0050.725∗0.027-0.0840.8300.3850.3060.3950.2920.5150.156

3.2 土壤细菌群落组成变化

根据表1的结果,煤矸石充填复垦对OTUs、细菌群落丰富度(ACE和Chao1)和多样性(Shannon)指数有显著影响(P<0.05)。在本研究中,煤矸石充填复垦后的总OTUs和细菌多样性指数(包括ACE,Chao和Shannon)在SL层均增加,但在DL层出现显著降低。SL层中由于较高的灌水量和施肥量,促进了土壤细菌的生长。此外作物种植也会为土壤细菌提供更多生长营养成分,与前人研究[29-30,36]一致。但由于煤矸石充填质量差,营养成分缺失,在DL层复垦地的土壤多样性和丰富度骤减。另外,整体上G的土壤细菌群落多样性是低于CK的,究其原因可能是:一方面,由于煤矸石充填造成土壤肥力下降,影响了有机质和矿物质的分解,抑制了细菌生长和繁殖[31-32];另一方面,由于研究区属于高潜水位采煤塌陷地,复垦时间也较长,研究区存在一定盐渍化威胁,且土壤普遍呈碱性,在一定程度上抑制了细菌活性,导致复垦地土壤细菌群落多样性较低[37]

变形菌门、放线菌门和酸杆菌门是研究区内最具优势的门,这与前人研究一致[38-39],说明这些门并不是某一区域特有,不同区域相对丰度不同。在G中观察到放线菌随土层深度增加呈现上升趋势。LI等[40]在研究矿区重金属污染环境时得出放线菌适合在低营养环境中生长。而酸杆菌的相对丰度随土层深度下降而降低,在G的DL层观察到的相对丰度平均值仅有0.7%。这与煤矸石充填质量差,且研究区土壤偏碱性,土壤结构被破坏,进而影响了酸杆菌的生长环境有关[41]。另外,由于煤矸石填充所造成的土层扰动,也令酸杆菌在各个土层的分布远不及普通耕地平稳[30]。在G和CK的各层均观察到较高的变形菌门相对丰度。据研究,变形菌在草地[37]、森林[38]和牧场土壤[42]的细菌群落中占主导地位,近年来,变形菌在矿区中也被判断为优势菌[38-39]。此外,由于煤矸石充填影响土壤结构稳定和营养输送,复垦土壤中细菌间的竞争压力增大,所以G的优势菌门在DL层相对丰度远低于其他土层。

3.3 土壤理化性质与土壤细菌群落关系变化

CCA(图4)和Pearson相关系数(表3,4)显示了土壤细菌群落与土壤理化性质之间的关系。在本研究中,pH是影响复垦地与对照地土壤细菌群落的主要因素。因为高潜水位特有的水土条件导致土壤存在盐渍化风险,尤其是研究区的复垦地中充填的煤矸石中又含有大量的碳酸盐矿物,碳酸盐矿物溶解后溶液导致土壤呈弱碱性,pH过高,将破坏土壤结构,导致土壤通气状况变差,间接影响土壤养分转化及其有效性[43]。通过图4可以看出,在煤矸石充填复垦地中,MoiC和SOC是影响细菌群落结构的主要理化性质。已有研究证明煤矸石结构性差,大孔隙多,毛细孔隙少,土壤保水能力差[44]。而煤矸石充填复垦土壤会影响含水量的下渗速率,抑制细菌等微生物活跃度,进而影响细菌群落相对丰度稳定性[32]。研究发现,煤矸石充填复垦地的SOC是细菌群落能量和营养物的稳定来源,细菌分解代谢速率与SOC含量相关。尹宁宁等[45]在研究矿区SOC与土壤微生物功能多样性之间关系时也有同样的结论。另外,AP和AK是影响对照地土壤细菌群落结构的主要理化性质。长期施肥效应对细菌群落的影响不容忽视,AP的施用能增加土壤细菌多样性以及各细菌物种在群落中分布的均匀性[46]。而钾肥的生物可用性一直被认为是调节农田生态系统产量的关键因素,PENG等[47]的研究中发现AK是影响油菜籽根际微生物群落的重要因子。也有研究认为AK含量与农田黑土细菌各细菌群落生物量呈高度正相关[48]

根据Pearson相关分析结果,发现变形菌门与AP存在显著相关性(P<0.05),说明变形菌门具有嗜营养的特点[10]。芽单胞菌门具有抗极端环境条件和清除土壤中重金属的能力,能适应碱性土壤环境,与复垦地pH呈正相关[30]。而拟杆菌门因为不能适应碱性环境与pH呈负相关[49]。另外,MoiC与复垦地中变形菌门、酸杆菌门、绿弯菌门等优势菌门的相对丰度均呈显著正相关,这与煤矸石填充土壤保水能力差,影响细菌群落生长繁殖有关[44]。此外,优势菌门仅与部分理化性质显著相关(P<0.05),放线菌门、酸杆菌门、绿弯菌门、厚壁菌门对大多数土壤理化性质的变化无明显响应。

4 结 论

(1)由于受煤矸石自身特性及高潜水位特有的水土条件影响,煤矸石充填复垦地的土壤细菌群落多样性与丰富度低于普通耕地,且随土层深度增加呈显著下降趋势,40~60 cm土层细菌群落多样性与丰富度明显低于其他土层。

(2)煤矸石充填复垦地土壤细菌群落组成与普通耕地有很高的相似性,在不同土层深度下,土壤细菌群落结构的变化主要体现在优势菌群落比例的变化。在所有的处理中,优势菌门是变形菌门、放线菌门和酸杆菌门。与对照地相比,复垦地表层的优势菌相对丰度较高,但经过煤矸石充填复垦工程扰动后,土壤生态环境受到严重影响,尤其在40~60 cm层,细菌等微生物难以在煤矸石层正常生长,细菌群落组成差异变大。

(3)煤矸石充填复垦显著改变了研究区土壤理化性质,不同土层深度下土壤理化性质变化趋势差异明显。相关性分析表明,pH,MoiC,SOC,AP和AK是影响土壤细菌群落组成的主要土壤理化性质。但部分优势菌门对理化性质的响应并不明显,可能与土壤微生物群落的高度复杂性及该区域条件的特殊性等多方面因素相关。

参考文献(References):

[1] ZHEN Qing,ZHENG Jiyong,ZHANG Xingchang,et al.Changes of solute transport characteristics in soil profile after mining at an opencast coal mine site on the Loess Plateau,China[J].Science of the Total Environment,2019,669:142-153.

[2] 胡振琪,李勇,陈洋,等.2019年土地科学研究重点进展评述及2020年展望——土地工程与信息技术分报告[J].中国土地科学,2020,34(2):93-102.

HU Zhenqi,LI Yong,CHEN Yang,et al.Progress review on land science research in 2019 and prospects for 2020:The sub-report on land engineering and land information technology[J].China Land Science,2020,34(2):93-102.

[3] MA Dan,DUAN Hongyu,LIU Jiangfeng,et al.The role of gangue on the mitigation of mining-induced hazards and environmental pollution:An experimental investigation[J].Science of the Total Environment,2019,664:436-448.

[4] 胡振琪,肖武,赵艳玲.再论煤矿区生态环境“边采边复”[J].煤炭学报,2020,45(1):351-359

HU Zhenqi,XIAO Wu,ZHAO Yanling.Re-discussion on coal mine eco-environment concurrent mining and reclamation[J].Journal of China Coal Society,2020,45(1):351-359.

[5] JIANG Xue,LU Wenxi,ZHAO Haiqing,et al.Potential ecological risk assessment and prediction of soil heavy-metal pollution around coal gangue dump[J].Natural Hazards and Earth System Sciences,2014,14(6):1599-1610.

[6] 白中科,周伟,王金满,等.试论国土空间整体保护、系统修复与综合治理[J].中国土地科学,2019,33(2):1-11.

BAI Zhongke,ZHOU Wei,WANG Jinman,et al.Overall protection,systematic restoration and comprehensive management of land space[J].China Land Science,2019,33(2):1-11.

[7] WANG Jinman,LI Xinfeng,BAI Zhongke,et al.The effects of coal ga-ngue and fly ash on the hydraulic properties and water content distribution in reconstructed soil profiles of coal-mined land with a high groundwater table[J].Hydrological Processes,2017,31(3):687-697.

[8] 张变华,靳东升,张强,等.不同植物种植对矿区复垦土壤微生物多样性的影响[J].农业资源与环境学报,2019,36(3):355-360.

ZHANG Bianhua,JIN Dongsheng,ZHANG Qiang,et al.Effects of different plant cultivations on the microbiological diversity of recultivated soil in mining area[J].Journal of Agricultural Resources and Environment,2019,36(3):355-360.

[9] BLAGODATSKAYA E,KUZYAKOV Y.Active microorganisms in soil:Critical review of estimation criteria andapproaches[J].Soil Biology and Biochemistry,2013,67:192-211.

[10] 李金融,侯湖平,王琛,等.基于高通量测序的复垦土壤细菌多样性研究[J].环境科学与技术,2018,41(12):148-157.

LI Jinrong,HOU Huping,WANG Chen,et al.Soil bacteria diversity of reclaimed soil based on high throughput sequencing[J]. Environmental Science and Technology,2018,41(12):148-157.

[11] 毕银丽,申慧慧.西部采煤沉陷地微生物复垦植被种群自我演变规律[J].煤炭学报,2019,44(1):307-315.

BI Yinli,SHEN Huihui.Effect of micro-reclamation on different planted forest on the vegetation self-succession in the western mining subsidence area[J].Journal of China Coal Society,2019,44(1):307-315.

[12] JACOB T O,OGHENEKUME G E,ASHTON K C.Bacterial degradation of coal discard and geologically weathered coal[J]. International Journal of Coal Science & Technology,2020,7(2):405-416.

[13] SEKHOHOLA L M,COWAN A K.Biological conversion of low-grade coal discard to a humic substance-enriched soil-like material[J].International Journal of Coal Science & Technology,2017,4(2):183-190.

[14] GOLEBIEWSKI M,DEJA-SIKORA E,CICHOSZ M,et al.16S rDNA pyrosequencing analysis of bacterial community in heavy metals polluted soils[J].Microbial Ecology,2014,67(3):635-647.

[15] PRAMANIK K,MITRA S,SARKAR A,et al.Alleviation of phytotoxic effects of cadmium on rice seedlings by cadmium resistant PGPR strain Enterobacter aerogenes MCC 3092[J].Journal of Hazardous Materials,2018,351:317-329.

[16] YOUNG J M,SKVORTSOV T,KELLEHER B P,et al.Effect of soil horizon stratigraphy on the microbial ecology of alpine paleosols[J].Science of the Total Environment,2018,657:1183-1193.

[17] GUO Yanqing,CHEN Xiaotian,WU Yuanyuan,et al.Natural revegetation of a semiarid habitat alters taxonomic and functional diversity of soil microbial communities[J].Science of the Total Environment,2018,635:598-606.

[18] 林耀奔,叶艳妹,吴次芳,等.水田土壤细菌群落对不同重金属污染水平的响应分析——以A县为例[J].环境科学学报,2020,40(1):224-233.

LIN Yaoben,YE Yanmei,WU Cifang,et al.Response analysis of soil bacterial community to different heavy metal pollution levels in paddy fields:A case study of A county[J].Acta Scientiae Circumstantiae,2020,40(1):224-233.

[19] GUO Yundan,LU Yuzhong,SONG Yuying,et al.Concentration and characteristics of dissolved carbon in the Sanjiang Plain influenced by long-term land reclamation from marsh[J].Science of the Total Environment,2014,466-467:777-787.

[20] CHEN Zhan,MALTZ M R,CAO Jixin,et al.Elevated O3 alters soil bacterial and fungal communities and the dynamics of carbon and nitrogen[J].Science of the Total Environment,2019,677:272-280.

[21] 李明,毕江涛,王静.宁夏不同地区盐碱化土壤细菌群落多样性分布特征及其影响因子[J].生态学报,2020,40(4):1316-1330.

LI Ming,BI Jiangtao,WANG Jing.Bacterial community structure and key influence factors in saline soil of different sites in Ningxia[J].Acta Ecologica Sinica,2020,40(4):1316-1330.

[22] 王瑾.西部煤矿区开采扰动对根际微生态影响及微生物复垦效应[D].北京:中国矿业大学(北京),2015:49-52.

WANG Jin.Impacts of coal mining disturbance on rhizosphere microecological and microbial reclamation effect in western area of China[D].Beijing:China University of Mining and Technology-Beijing,2015:49-52.

[23] KNELLER T,HARRIS R J,BATEMAN A,et al.Native-plant amendments and topsoil addition enhance soil function in post-mining arid grasslands[J].Science of the Total Environment,2018,621:744-752.

[24] 郭晓明,赵同谦.采煤沉陷区耕地土壤微生物数量及酶活性的空间特征[J].环境工程学报,2010,4(12):2837-2842.

GUO Xiaoming,ZHAO Tongqian.Spatial characters of soil microbial quantities and enzyme activity in coal mining subsided area[J]. Chinese Journal of Environmental Engineering,2010,4(12):2837-2842.

[25] 李兆龙,梁红,刘文,等.稀土矿区生态修复过程中的土壤改良及细菌群落变化[J].仲恺农业工程学院学报,2013,26(1):9-13.

LI Zhaolong,LIANG Hong,LIU Wen,et al.Soil amelioration and bacterial colony changes during ecological restoration at rare earth mining land[J].Journal of ZhongKai University of Agriculture and Technology,2013,26(1):9-13.

[26] LI Junjian,ZHOU Xiaomei,YAN Junxia,et al.Effects of regenerating vegetation on soil enzyme activity and microbial structure in reclaimed soils on a surface coal mine site[J].Applied Soil Ecology,2015,87:56-62.

[27] 吴克宁,赵瑞.土壤质地分类及其在我国应用探讨[J].土壤学报,2019,56(1):227-241.

WU Kening,ZHAO Rui.Soil texture classification and its application in China[J].Acta Pedologica Sinica,2019,56(1):227-241.

[28] 王坤.高潜水位采煤塌陷区充填复垦土壤碳动态研究[D].徐州:中国矿业大学,2016:14-15.

WANG Kun.Study on carbon dynamics of reclaimed soil in coal mining subsidence area with high groundwater level[D]. Xuzhou:China University of Mining and Technology,2016:14-15.

[29] HAMONTS K,BISSETT A,MACDONALD B C T,et al.Effects of ecological restoration on soil microbial diversity in a temperate grassy woodland[J].Applied Soil Ecology,2017,117/118:117-128.

[30] HOU Huping,WANG Chen,DING Zhongyi,et al.Variation in the soil microbial community of reclaimed land over different reclamation periods[J].Sustainability,2018,10(7):2286.

[31] 刘曙光,徐良骥,余礼仁.煤矸石充填复垦地土壤理化性质的时空变化特征[J].环境工程,2018,36(2):173-177.

LIU Shuguang,XU Liangji,YU Liren.Spatiotemporal variation of physical and chemical properties of reclaimed soil in coal fangue filling land[J].Environmental Engineering,2018,36(2):173-177.

[32] 郑永红,张治国,姚多喜,等.煤矸石充填复垦对土壤特性影响研究[J].安徽理工大学学报(自然科学版),2013,33(4):7-11.

ZENG Yonghong,ZHANG Zhiguo,YAO Duoxi,et al.Study on influence of gangue on reclaimed soil properties[J].Journal of Anhui University of Science and Technology(Natural Science),2013,33(4):7-11.

[33] 鲁叶江.压实处理对矸石充填复垦土壤水分的影响[J].煤炭科学技术,2010,38(11):125-128.

LU Yejiang.Compaction treatment affected to moisture content of coal reject backfilling reclaimed land[J].Coal Science and Technology,2010,38(11):125-128.

[34] 焦晓燕,王立革,卢朝东,等.采煤塌陷地复垦方式对土壤理化特性影响研究[J].水土保持学报,2009,23(4):123-125,145.

JIAO Xiaoyan,WANG Lige,LU Chaodong,et al.Effects of two reclamation methodologies of coal mining subsidence on soil physical and chemical properties[J].Journal of Soil and Water Conservation,2009,23(4):123-125,145.

[35] 徐良骥,黄璨,章如芹,等.煤矸石充填复垦地理化特性与重金属分布特征[J].农业工程学报,2014,30(5):211-219.

XU Liangji,HUANG Can,ZHANG Ruqin,et al.Physical and chemical properties and distribution characteristics of heavy metals in reclaimed land filled with coal gangue[J].Transactions of the Chinese Society of Agricultural Engineering,2014,30(5):211-219.

[36] 杨亚东,王志敏,曾昭海.长期施肥和灌溉对土壤细菌数量、多样性和群落结构的影响[J].中国农业科学,2018,51(2):290-301.

YANG Yadong,WANG Zhimin,ZENG Zhaohai.Effects of long-term different fertilization and irrigation managements on soil bacterial abundance,diversity and composition[J].Scientia Agricultura Sinica,2018,51(2):290-301.

[37] CHENG Zhibo,CHEN Yun,ZHANG Fenghua.Effect of reclamation of abandoned salinized farmland on soil bacterial communities in arid northwest China[J].Science of the Total Environment,2018,630:799-808.

[38] ZENG Quanchao,AN Shaoshan,LIU Yang.Soil bacterial community response to vegetation succession after fencing in the grassland of China[J].Science of the Total Environment,2017,609:2-10.

[39] WANG Ziting,LU Liu,CHEN Qing,et al.Conservation tillage increases soil bacterial diversity in the dryland of northern China[J]. Agronomy for Sustainable Development,2016,36(2):28.

[40] LI Yuanyuan,CHEN Longqian,WEN Hongyu.Changes in the composition and diversity of bacterial communities 13 years after soil reclamation of abandoned mine land in eastern China[J]. Ecological Research,2015,30:357-366.

[41] 徐嘉兴,李钢,陈国良,等.矿区土地生态质量评价及动态变化[J].煤炭学报,2013,38(S1):180-186.

XU Jiaixng,LI Gang,CHEN Guoliang,et al.Dynamic changes and evaluation of land ecological quality in coal mining area[J]. Journal of China Coal Society,2013,38(S1):180-186.

[42] 刘君,王宁,崔岱宗,等.小兴安岭大亮子河国家森林公园不同生境下土壤细菌多样性和群落结构[J].生物多样性,2019,27(8):911-918.

LIU Jun,WANG Ning,CUI Daizong,et al.Community structure and diversity of soil bacteria in different habitats of Da Liangzihe National Forest Park in the Lesser Khinggan Mountains[J]. Biodiversity Science,2019,27(8):911-918.

[43] 谷得明.煤矸石堆存对地表与浅层地下水环境的影响研究[D].淮南:安徽理工大学,2015:37-43.

GU Deming.Research of coal gangue storage impact on surface and shallow groundwater environment—A case of Panji mining area in Huainan[D].Huainan:Anhui University of Sience and Technology,2015:37-43.

[44] 黄凯.煤矸石充填复垦土壤水分含量与作物有效性的研究[J].安徽农业科学,2014,42(6):1598-1599,1602.

HUANG Kai.Moisture content of reclaimed soil filled with coal gaugue and its validation of crops[J].Journal of Anhui Agricultural Science,2014,42(6):1598-1599,1602.

[45] 尹宁宁,王丽萍.复垦矿区土壤有机碳和微生物活性变化及其解析[J].环境科技,2014,27(6):5-13.

YIN Ningning,WANG Liping.The changes and analysis of organic carbon and microbial activity in reclaimed soil in mine area[J]. Environmental Science and Technology,2014,27(6):5-13.

[46] MOROR H,KUNITO T,SAITO T,et al.Soil microorganisms are less susceptible than crop plants to potassium deficiency[J]. Archives of Agronomy and Soil Science,2014,60:1807-1813.

[47] PENG Caijuan,LAI Songsong,LUO Xuesong,et al.Effects of long term rice straw application on the microbial communities of rapeseed rhizosphere in a paddy-upland rotation system[J].Science of the Total Environment,2016,557:231-239.

[48] 白震,张明,张旭东,等.长期磷肥处理农田黑土细菌群落特征[J].应用生态学报,2008,19(7):1567-1573.

BAI Zhen,ZHANG Ming,ZHANG Xudong,et al.Characteristics of bacterial community in an arable mollisol under long-term phosphorus fertilization[J].Chinese Journal of Applied Ecology,2008,19(7):1567-1573.

[49] 侯建伟,邢存芳,邓晓梅,等.pH对花椒根区土壤细菌群落结构的影响[J].西北农林科技大学学报(自然科学报),2020,48(5):115-122.

HOU Jianwei,XING Cunfang,DENG Xiaomei,et al.Effect of pH on soil bacterial community structure in root zone of prickly ash[J]. Journal of Northwest A&F University(Natural Science Edition),2020,48(5):115-122.