环境科学  2021, Vol. 42 Issue (4): 2066-2079   PDF    
微生物群落驱动AM真菌、生物炭及联合改良沙化土壤作用潜力
张哲超, 杨久扬, 郝百惠, 郝利君, 罗俊清, 李雪, 刁风伟, 张璟霞, 郭伟     
内蒙古大学生态与环境学院, 蒙古高原生态学与资源利用教育部重点实验室, 内蒙古草地生态学重点实验室, 呼和浩特 010021
摘要: 沙化土壤作为土地荒漠化的重要过渡形式,实现其有效恢复对减缓土地荒漠化进程意义重大.本研究显示,丛枝菌根(arbuscular mycorrhizal,AM)真菌和生物炭已应用于沙化土壤改良过程,但二者联合对沙化土壤改良影响研究较少;此外,细菌和真菌群落在沙化土壤改良过程中的作用尚不清楚.采用温室盆栽试验的方法,分别设置对照处理(CK)、单独接种AM真菌处理(RI)、单独施加生物炭处理(BC)和二者联合改良处理(BC_RI),研究不同改良方式对小果白刺(Nitrariasi birica Pall.)菌根侵染率、生物量、矿质营养元素(N、P、K、Ca和Mg)含量及土壤有机碳、营养元素(全N、全P和全K)、水稳性团聚体组成影响.采用高通量测序技术,考察细菌和真菌群落在沙化土壤改良过程中的作用,结合多元分析手段,探究不同改良方式改良作用机制,旨在为合理有效改良沙化土壤提供基础数据和理论依据.结果表明,接种处理(RI和BC_RI)小果白刺根系均观察到明显的菌根侵染现象,但RI和BC_RI处理间菌根侵染率无显著性差异.与CK相比,RI处理显著增加了小果白刺地上部生物量和地上部N、K、Ca和Mg含量,BC和BC_RI处理显著增加了小果白刺地上部、根部生物量及N、P、K、Ca和Mg含量;BC_RI处理与RI和BC相比,根部生物量及P、K、Ca和Mg含量显著增加.与CK相比,BC和BC_RI处理显著增加了土壤有机碳含量,RI处理使得土壤全N含量显著增加了152.54%,BC处理使得土壤全P和全K含量分别显著降低了12.5%和18.18%.各处理0.25~0.05 mm粒径土壤团聚体比例最高,BC_RI处理能够显著促进大粒径(>0.25 mm)土壤团聚体形成.与CK相比,RI和BC_RI处理显著降低细菌、真菌群落Sobs指数和Shannon指数;各处理细菌及真菌菌门组成及丰度存在显著差异性.RDA及网络分析结果显示,AM真菌、生物炭及二者联合改良方式对土壤基质环境及土壤微生物群落结构影响显著,不同改良方式下微生物分子生态网络关系发生显著变化,不同改良处理中核心物种组成具有差异性;BC和BC_RI较RI处理,网络连接密集程度更高且核心物种组成更丰富;生物炭与AM真菌联合,弱化了Rhizophagus intraradices的核心作用,并增强其他微生物(特别是细菌菌种)的核心地位.SEM结果显示,AM真菌和生物炭施用通过直接影响土壤细菌和真菌群落结构,进而实现对植物生长及土壤性质变化的影响,微生物群落结构差异(特别是微生物间的互作关系变化)是导致土壤改良效果差异的主要驱动力.综上所述,不同改良方式对沙化土壤改良效果影响具有差异性,微生物群落在土壤改良过程中具有关键影响作用,AM真菌和生物炭联合对加速沙化土壤生态恢复具有潜在优势和应用价值.
关键词: 微生物群落      丛枝菌根真菌      生物炭      联合改良      沙化土壤      改良机制     
Potential of Arbuscular Mycorrhizal Fungi, Biochar, and Combined Amendment on Sandy Soil Improvement Driven by Microbial Community
ZHANG Zhe-chao , YANG Jiu-yang , HAO Bai-hui , HAO Li-jun , LUO Jun-qing , LI Xue , DIAO Feng-wei , ZHANG Jing-xia , GUO Wei     
Inner Mongolia Key Laboratory of Environmental Pollution Control and Waste Resource Recycle, Inner Mongolia Key Laboratory of Grassland Ecology, School of Ecology and Environment, Inner Mongolia University, Hohhot 010021, China
Abstract: Sandy soils are considered as a significant transition phase to desertification. The effective recovery of sandy soils is of great significance to mitigate the desertification process. Some studies have shown that arbuscular mycorrhizal (AM) fungi and biochar improved the sandy soil, but there have been very few studies regarding the combined effects of AM fungi and biochar amendments on sandy soil improvement. Additionally, the roles of the bacterial and fungal community during the process of sandy soil improvement remain unclear. A greenhouse pot experiment with four treatments, including a control (CK, no amendment), single AM fungi-assisted amendment (RI), single biochar amendment (BC), and combined amendment (BC_RI, biochar plus AM fungi), was set up. This study investigated the effects of different amendment methods on the Nitrariasi birica mycorrhizal colonization, biomass, nutrient (N, P, K, Ca, and Mg) content, soil organic carbon, soil nutrient (TN, TP, and TK) content, and soil water-stable aggregate composition. High throughput sequencing technology was used to investigate the roles of the bacterial and fungal communities during the process of sandy soil improvement. Combined with multiple analysis methods, the improvement mechanisms of different amendment methods were explored. The aim was to provide basic data and scientific basics for reasonably and effectively improving sandy soils. The results indicated that a significant mycorrhiza colonization was observed in the inoculation (RI and BC_RI) treatments, but there was no substantial difference in the mycorrhiza colonization with the RI and BC_RI. Compared with the CK, the shoot biomass and shoot element (N, K, Ca, and Mg) contents were significantly increased in the RI, and the shoot element (N, P, K, Ca, and Mg) contents were significantly increased in the BC and BC_RI; compared with the RI and BC, the root biomass and the root element (P, K, Ca, and Mg) contents were significantly increased in the BC_RI. Compared with the CK, the soil organic carbon contents were significantly increased in the BC and BC_RI, the soil TN contents were significantly increased by 152.54%, and the soil TP and TK contents were significantly decreased by 12.5% and 18.8%, respectively. The proportion of soil aggregates with particle sizes of 0.25-0.05 mm was the highest in each treatment, and the large particle size (>0.25 mm) soil aggregate was significantly increased in the BC_RI. Compared with the CK, the Sobs and Shannon indices of the bacterial/fungal community were significantly decreased in the RI and BC_RI. There was a difference in the microbial community compositions and abundance in the various treatments. The results of the RDA and network analysis were as follows: the effects of AM fungi, biochar, and combined amendment on the soil environment and microbial community structure were significant; in the different amendment treatments, the relationship of the microbial molecular ecological network was significantly changed, and the composition of the core species varied; compared with the RI, there was a higher network connection degree and a richer core species composition in the BC and BC_RI; moreover, the essential role of Rhizophagus intraradices was weaken and the core roles of the other microorganisms (especially bacterial species) were enhanced under the combined effects of biochar and AM fungi. The SEM results demonstrated that the application of AM fungi and biochar could directly affect the bacteria/fungi community structure, and further affect the plant growth and soil properties. The differences in the microbial community structure (especially the change in the microbial interaction) were the key driving factors that led to the difference in the soil improvement effectiveness. In summary, the effects of the different amendment methods on the improvement effectiveness of sandy soils varied. The microbial community played key roles in the process of sandy soil improvement, and there were potential advantages and applications in accelerating the ecological restoration of sandy soils under the combined AM fungi and biochar amendment.
Key words: microbial community      arbuscular mycorrhizal (AM) fungi      biochar      combined amendment      sandy soil      improvement mechanism     

随着人口增长、人类活动和全球气候变化等因素影响, 加剧了世界各地土地的退化和沙化, 土地荒漠化已成为目前我国乃至世界公认的最严重的生态环境问题之一[1].中国是世界上受土地退化危害最严重的国家之一, 我国荒漠化(沙化)土地分布较广, 主要集中于北方干旱、半干旱地区和部分半湿润地区[2].截至目前我国荒漠化土地面积已达2 620 000 km2, 沙化土地面积已达173 000 km2[3], 其中土地沙化作为土地荒漠化的重要过渡形式, 较荒漠化土地更易被恢复[4].因此, 实现沙化土地的有效恢复, 对于减缓土地荒漠化进程、改善退化土壤危害具有重要意义.

受沙化土壤理化性质差(如养分匮乏、保水保肥性能差等)和土壤生物群落复杂程度低等因素限制, 沙化土壤的自然恢复或植物群落修复进程缓慢[5].因此, 有必要通过施加外源调节剂或改良剂加速沙化土壤的生态恢复[5].生物炭由于其具有富碳、疏松多孔、比表面积大、持水性高和吸附能力强等优点, 常被作为理想的土壤改良剂应用于退化土壤改良过程[6, 7].有研究表明, 生物炭施用能够显著影响土壤理化性质(如改善土壤容重、土壤水分条件、土壤酸碱度和提高土壤养分等)和土壤生物学特性(如改变土壤微生物丰度、微生物群落结构和微生物功能等)[8~11], 但由于受生物炭材料本身性质、热解过程和土壤性质等不同的影响[12], 其对于沙化土壤理化性质改良效果具有明显的差异性[13~19].此外, 目前关于生物炭施用对沙化土壤微生物群落影响研究还较少.然而, 土壤微生物群落对于调控土壤物质循环和改善土壤质量具有关键性作用.因此, 除土壤基本理化性质变化外, 有必要进一步关注到土壤微生物群落结构变化, 从微生物作用角度完善生物炭施用对沙化土壤改良的作用机制.丛枝菌根(arbuscular mycorrhizal, AM)真菌对于极端环境中植物定植、植物养分吸收利用、土壤抗侵蚀能力和土壤碳储存等方面存在着有益影响, 因此也常被与植物联合, 用于退化土壤改良[20, 21].已有研究显示, AM真菌对促进退化土壤(如沙化土壤、矿区土壤、盐碱土壤)植物生长、养分吸收[22~24]以及污染土壤、盐碱土壤的生态恢复发挥着积极作用[25, 26], 然而对沙化土壤性质(基本理化性质、土壤微生物群落)改良及生态恢复方面关注较少.此外, 有研究发现土壤中生物炭容易导致植物和菌根真菌之间响应关系的重大变化, 如通过影响土壤理化性质改变土壤养分有效性, 进而影响AM真菌定植[27]; 亦或是增强AM真菌协助宿主植物抵抗病原侵染的能力[28].外源生物炭施加也会对菌根真菌或外生菌根生长产生增效作用, 进而对土壤质量及植物生长产生促进作用[29~31].基于此, 构建植物-微生物-土壤改良剂之间的协同过程, 对于加速土壤恢复具有巨大潜能.因此, 深入研究植物-AM真菌-生物炭修复沙化土壤作用机制, 对日后建立沙化土壤的高效修复技术, 加速退化土壤生态恢复具有十分重要的科学意义和实际应用价值.

土壤微生物作为土壤基质与植物间物质、能量和信息传递的桥梁, 对土壤性质改变、土壤养分循环及植物生长等方面具有重要影响[32, 33].尽管, 目前关于AM真菌和生物炭单独改良退化土壤方面的作用已有相关研究, 但二者联合的改良作用研究较少[34], 且很少有研究关注不同微生物群落(特别是真菌群落)在沙化土壤改良中的潜在作用.基于此, 本研究采用温室盆栽试验, 通过设置AM真菌和生物炭单独或联合处理的方式, 探究不同改良方式对植物生长、沙化土壤性质及微生物群落结构影响; 从土壤微生物作用角度探究AM真菌、生物炭施用及二者交互作用对沙化土壤改良的影响作用机制; 并通过多元分析手段, 识别引起改良效果差异的主要驱动力, 厘清微生物(AM真菌、细菌和真菌)、土壤及植物互作关系, 以期为合理有效改良沙化土壤提供基础数据和理论依据.

1 材料与方法 1.1 培养基质

供试土壤取自内蒙古呼和浩特西部土默特左旗小营子村(N40°30′48.60″, E111°45′59.69″)附近, 收集表层(0~20 cm)无污染风沙土, 去除表面杂质(石块、枯枝落叶等), 并自然风干、研磨过筛(2.0 mm), 土壤基本性质见表 1, 其中沙化土壤水稳性团聚2~0.25、0.25~0.05、0.05~0.01和0.01~0.005 mm粒径所占比例分别为11.37%、77.54%、9.41%和1.12%.在盆栽试验前对土壤进行高温灭菌处理(121℃, 2 h), 灭菌后自然风干保存备用.

表 1 供试沙化土壤基本理化性质 Table 1 Physical and chemical properties of the sandy soil

1.2 供试植物、菌种与生物炭改良剂

小果白刺(Nitrariasi birica Pall.)是荒漠生态系统的优势物种之一, 是一种典型的沙漠灌丛植物, 具有独特的防风固沙、耐干旱、抗风蚀和盐碱等功能[35]; 此外, 由于其根系发达, 在植物固沙和土壤改良等应用方面一直备受关注[36].基于此, 本研究选用小果白刺为供试植物, 小果白刺种子购于中国宁夏远胜陆洋草业生态工程有限公司.在进行温室盆栽试验前, 用10%体积比过氧化氢(H2O2)浸泡15 min, 以对种子表面进行消毒.之后在25℃潮湿滤纸上催芽萌发48~72 h, 直到露出胚根后进行播种.本研究供试的AM真菌为根内球囊霉[Rhizophagus intraradices, RI, BGC BJ09, 1511C0001BGCAM0042, 370个孢子·(10 g)-1].该菌剂购买于北京农林科学院植物营养与资源研究所, 以玉米(Zea mays L.)和白三叶草(Trifoliumrepens L.)混合栽种形式对菌株进行扩繁培养3个月, 获得内含真菌孢子、菌丝以及宿主植物根段等繁殖体的根际砂土混合物[37].基于前人研究发现秸秆类生物炭对沙土改良的有益作用[38, 39].因此, 本研究选用秸秆类生物炭为土壤改良剂, 生物炭购于江苏省溧阳市德胜活性炭厂, 其生产基本技术指标见表 2.在盆栽试验前对生物炭进行高温灭菌处理(121℃, 2 h), 灭菌后自然风干保存备用.

表 2 秸秆生物炭生产技术指标 Table 2 Technical indicators of the straw biochar

1.3 温室盆栽试验设计

本试验采用生物学盆栽法, 处理分为:对照处理(CK)、单独接种RI处理(RI)、单独施加生物炭处理(BC)和接种RI与生物炭联合改良处理(BC_RI)4个处理, 每个处理设置6盆重复(共计24盆).以圆形花盆(上口径25 cm×下口径15 cm×高20 cm)为培养容器, 内衬无菌塑料自封袋.每盆放入4 kg灭菌沙化土壤, 依据已有研究生物炭添加最佳施用[40]及预试验结果, BC处理生物炭添加量选择20 g·kg-1, AM真菌接种处理采用混接法进行, 其中接菌处理(即RI和BC_RI)每盆每千克土壤加入50 g菌剂[41], 不接菌处理(即CK和BC)加入同等质量的灭菌菌剂和10 mL的菌剂过滤液(10 μm滤膜过滤), 所有添加物均保证与土壤基质充分混匀.将每盆播种15粒均匀预发芽的小果白刺种子, 15 d后间苗, 每盆保留6株长势相近的植株.所有处理随机摆放于内蒙古大学的温室中, 试验期间采用自然采光(日间温度为20~30℃, 夜间温度为10~20℃).通过称重法定期添加去离子水, 维持土壤基质含水量为田间最大持水量的80%[37], 自出苗之日起植物生长150 d后收获.

1.4 样品制备与基本理化指标分析测定

盆栽试验结束后, 分别收获植物地上部分以及地下部分(根系)样品, 用蒸馏水彻底冲洗去除附着土壤, 70℃烘干至恒重, 称重记录地上部分以及地下部分植物样品干重.在收获时, 随机选取部分新鲜根段(0.7 g)保存于乙醇(50%), 用于测定菌根侵染率[42].植物样品烘干称重后粉碎过筛(< 150 μm), 室温保存用于后续植物相关指标测定.称取研磨过筛植物样品(2~3 mg), 利用元素分析仪(Vario MACRO, CHNOS Elemental Analyzer, Elementar Co., Germany)测定植物样品C和N含量.称取0.5 g研磨过筛植物样品, 加入5 mL优级纯HNO3于120℃开放式消化96 h, 利用电感耦合等离子体发射光谱仪(ICP-OES, Optima 7000DV, PerkinElmer USA)测定消煮液中P、K、Ca和Mg含量.

收获植物样品后, 将整盆土壤混合.从每个处理6盆重复中随机挑选3盆, 收集3盆重复的新鲜混合土壤样品.其中一部分土壤样品使用无菌离心管(50 mL)收集, 并立即转移到-80℃冰箱中储存, 用于后续的分子生物学分析; 剩余土壤样品风干、研磨、过筛后置于室温下保存, 以供后续土壤理化性质分析.利用HCl(1 mol·L-1)对土壤进行预处理以去除土壤中的无机碳组分, 将预处理后的土壤样品105℃烘干至恒重, 称取土壤样品(120 mg), 利用元素分析仪测定土壤样品总有机碳及全N含量[43].称取风干过筛后的土壤样品(1 g), 加入等体积微量金属级(trace metal grade)氯化氢(HCl)和硝酸(HNO3)对样品进行消解[44], 使用ICP-OES分析样品中全P和全K含量.称取风干过筛土壤样品, 分别按照固液比(W/V)1/2.5和1/5比例加入去离子水, 振荡混匀静置, 利用pH计(Seven Compact S210)和电导率仪(Seven Compact S230)分别测定土壤样品pH值及电导率(EC)[45].参照中国农业行业标准(NY/T 1121.20-2008), 采用湿筛法对土壤水溶性团聚体组分比例进行分离测定.沙化土壤粒径组成采用比重法进行测定[46].

1.5 分子生物学分析

采用FastDNA® SPIN试剂盒(MP Biomedicals, CA, USA)对新鲜土壤样品(0.5 g)中DNA进行提取, 通过微量紫外分光光度计(Nano Drop® ND-2000, Wilmington, DE, USA)测定DNA浓度和纯度.利用PCR热循环仪(ABI GeneAmp® 9700, CA, USA), 选用特异性功能基因引物(表 3)分别对细菌16S rRNA基因V3~V4区和真菌ITS rRNA基因ITS区域进行PCR扩增; 再用2%琼脂糖凝胶电泳对PCR扩增产物进行定量检测, 用凝胶回收试剂盒(Axygen Biosciences, Union City, CA, USA)对目标片段进行切胶回收.按照前述电泳定量结果, 通过荧光计(QuantusTM Fluorometer, Promega, USA)对PCR扩增回收产物进行荧光定量.按照上海美吉生物技术有限公司(Majorbio Bio-Pharm, 中国上海)标准, 将纯化后的扩增子以等分子量混合, 借助Illumina MiSeq平台, 分别对土壤细菌、真菌群落进行高通量测序分析.

表 3 PCR引物序列信息 Table 3 PCR primer sequence

1.6 数据处理分析

所有植物及土壤基本理化指标试验数据通过Excel计算平均值和标准误差, 并采用SPSS 17.0软件对数据进行统计分析:采用单因素方差分析(ANOVA)分析不同处理差异性, 并通过Duncan多重复极差检验差异显著性(P < 0.05).通过R语言(ggplot2包)对植物及土壤基本理化指标数据进行可视化.采用双因子方差分析AM真菌、生物炭施用和二者交互作用的显著性影响.

采用Mothur软件对原始16S rRNA或ITS rRNA基因测序序列进行去重复、质量过滤和优化.利用UPARSE算法(UPARSE 7.1 http://drive5.com/uparse/)将优化后的原始数据进行拼接, 将相似性大于97%的优化序列划分为一个操作分类单元(OTU), 并与16S rRNA数据库(Silva:Release132 http://www.arb-silva.de)和ITS rRNA数据库(Unite:Release 7.2 http://unite.ut.ee/index.php)进行比对.通过Ribosomal Database Project(RDP http://rdp.cme.msu.edu/)数据库对每个OTU代表序列进行识别和分类.为保证分析结果的有效性, 对所有样本的子样本序列按照最小样本序列数抽平进行后续数据分析[47].最终, 各处理每个样品得到475 044条(细菌)和63 754条(真菌)优化后的有效序列, 各样本分别包含1 637个有效的细菌OTUs和775个真菌OUTs.

利用Mothur(http://www.mothur.org/wiki)计算微生物多样性指数(Shannon)及群落丰富度指数(Sobs); 利用Qiime 1.7.0(http://qiime.org/index.html)进行主坐标分析(Principal coordinate analysis, PCoA), 以判定各处理微生物群落分布及结构组成差异; 利用Circos-0.67-7(http://circos.ca/)对样本与物种对应关系进行可视化, 并通过R语言(stats包)和Python(scipy包)计算, 对多组样本中物种进行克氏秩和检验(Kruskal-Wallis H test), 识别不同处理中具有显著性差异的物种; 通过计算RI、BC和BC_RI处理中物种间Spearman相关性系数, 利用Cytoscape 2.6.0软件分别构建不同改良方式下细菌、真菌物种网络关系图, 以探究识别不同改良方式下微生物群落作用模式; 采用Canoco5(version 5.0)软件, 利用冗余分析(redundancy analysis, RDA)探讨AM真菌、生物炭及联合施用对土壤结构、土壤性质及微生物群落结构差异的影响; 基于结构方程模型(structural equation modeling, SEM)构建AM真菌、生物炭、土壤性质、土壤结构以及微生物群落间作用关系, 利用SPSS套件Amos(version 24.0)进行SEM分析, 以探究AM真菌、生物炭施用对沙化土壤改良作用机制.

2 结果与分析 2.1 菌根侵染率、小果白刺生物量及矿质营养元素含量分析

各处理菌根侵染率、小果白刺生物量及营养元素(N、P、K、Ca和Mg)含量如图 1所示.不接种处理(CK和BC)小果白刺根系中未观察到菌根真菌的侵染; 接种R. intraradices处理中, RI和BC_RI处理间菌根侵染率无显著性差异[图 1(a)].双因子分析结果表明, 接种处理对菌根侵染率有极显著影响(P < 0.001), 而生物炭施加及二者交互作用对菌根侵染率无显著影响.

CK: 对照处理; RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理; 柱状图上不同小写字母表示同一类指标不同处理间在P < 0.05显著水平差异 图 1 各处理小果白刺根系菌根侵染率、生物量及矿质营养元素(N、P、K、Ca和Mg)含量 Fig. 1 Mycorrhizal colonization, biomass, and mineral nutrient element (N, P, K, Ca, and Mg) content in the different treatments

与CK相比, RI、BC和BC_RI处理使得小果白刺地上部干重分别显著增加了113.60%、264.08%和285.44%(P < 0.05), 使得小果白刺地下部干重分别增加了126.32%、568.42%和884.21%; 与RI相比, BC_RI处理使得小果白刺地上部与根部干重均显著增加(P < 0.05); 与BC相比, BC_RI处理仅使得小果白刺根部干重显著增加[P < 0.05, 图 1(b)1(c)].双因子分析结果表明, 接种处理对小果白刺地上部及根部干重有显著影响(P < 0.05), 生物炭施加对小果白刺地上部及根部干重均有极显著影响(P < 0.001), 而在二者交互作用无显著影响.

N、P、K、Ca和Mg是植物生长所必需的营养元素.与CK相比, RI处理使得小果白刺地上部的N、K、Ca和Mg元素含量显著增加(P < 0.05), 而对根部元素含量无显著影响, BC处理使得地上部的N、P、K、Ca和Mg含量以及根部的N、K、Ca和Mg含量均显著增加(P < 0.05), BC_RI处理使得地上部及根部的N、P、K、Ca和Mg含量显著增加(P < 0.05); 与RI相比, BC_RI处理使得地上部P和K含量以及根部N、P、K、Ca和Mg含量显著增加(P < 0.05); 与BC相比, BC_RI处理使得根部P、K、Ca和Mg含量显著增加[P < 0.05, 图 1(d)1(e)].双因子分析结果表明, 接种处理对小果白刺地上部N和Ca含量有显著影响(P < 0.05), 对地上部Mg及根部P、K和Mg含量有非常显著影响(P < 0.01), 对根部Ca含量有极显著影响(P < 0.001); 生物炭施加对小果白刺地上部N和Ca含量有显著影响(P < 0.05), 对根部N含量有非常显著影响(P < 0.01), 对地上部P和K含量及根部P、K、Ca和Mg含量有极显著影响(P < 0.001); 二者交互作用, 对小果白刺地上部N含量及根部P和Mg含量存在显著影响(P < 0.05), 对根部Ca含量有非常显著影响(P < 0.01).

2.2 土壤有机碳、营养元素含量和水稳性团聚体差异分析

不同处理土壤有机碳含量如图 2(a)所示.与CK相比, BC和BC_RI处理使得土壤有机碳含量分别显著增加了58.00%和63.00%(P < 0.05), 而RI处理对土壤有机碳含量无显著影响; 与RI相比, BC_RI处理使得土壤有机碳含量显著增加了45.08%(P < 0.05); 与BC相比, BC_RI处理相对土壤有机碳含量无显著影响[图 2(a)].双因子分析结果表明, 仅生物炭施加对土壤有机碳含量有极显著影响(P < 0.001), 而接种处理及二者交互作用对其无显著影响.

(a)土壤有机碳含量; (b)土壤营养元素含量; (c)土壤水稳性团聚体组成; CK: 对照处理; RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理; (a)和(b)柱状图上不同小写字母表示同一类指标不同处理间在P < 0.05显著水平差异; (c)柱状图内不同小写字母表示2~0.25 mm粒径水稳性团聚体组成不同处理间在P < 0.05显著水平差异, 不同大写字母表示0.25~0.05 mm粒径水稳性团聚体组成不同处理间在P < 0.05显著水平差异; 图例内星号表示该粒径土壤水稳性团聚体组成差异显著性, *表示P < 0.05 图 2 各处理土壤有机碳、土壤营养元素含量及土壤水稳性团聚体组成 Fig. 2 Soil organic carbon, soil nutrient element content, and soil water-stable aggregate composition in the different treatments

不同处理土壤营养元素全N、全P和全K含量如图 2(b)所示.与CK相比, RI处理使得土壤全N含量显著增加了152.54%(P < 0.05), BC处理使得土壤全P和全K含量分别显著降低了12.5%和18.18%(P < 0.05), BC_RI处理对土壤全N、全P和全K含量无显著影响; 与RI相比, BC_RI处理使得土壤全N含量显著降低了51.68%(P < 0.05); 与BC相比, BC_RI处理使得土壤全P含量显著增加了10.71%[P < 0.05, 图 2(b)].双因子分析结果表明, 接种处理对土壤N含量有非常显著影响(P < 0.01); 生物炭施用对土壤N、P和K含量有显著影响(P < 0.05); 二者交互作用对土壤N和P有非常显著影响(P < 0.01), 对土壤K有显著影响(P < 0.05).

各处理水稳性团聚体组成中, 0.25~0.05 mm粒径土壤团聚体比例最高, 约占总体水稳性团聚体的71%~79%.与CK相比, RI和BC处理对土壤水稳性团聚体组成变化无显著影响, BC_RI处理使得2~0.25 mm粒径及0.25~0.05 mm粒径土壤团聚体比例发生显著变化(P < 0.05), 而对其他粒径团聚体比例无显著影响; 与RI相比, BC_RI处理对土壤水稳性团聚体组成变化无显著影响; 与BC相比, BC_RI处理显著促进了土壤大粒径团聚体(>0.25 mm)的形成(P < 0.05), 显著降低了0.25~0.05 mm粒径土壤团聚体比例(P < 0.05), 对其他粒径土壤团聚体比例无显著影响[图 2(c)].双因子分析结果表明, 接种处理对2~0.25 mm粒径和0.25~0.05 mm粒径土壤团聚体比例变化存在显著(P < 0.05)和非常显著(P < 0.01)影响; 而生物炭施用和二者交互作用对土壤水稳性团聚体组成变化无显著影响.

2.3 微生物群落组成及多样性差异分析

各处理门(phylum)水平下物种群落组成如图 3所示.细菌群落中, 变形菌门(Proteobacteria)、酸根菌门(Acidobacteria)、拟杆菌门(Bacteroidetes)、绿弯菌门(Chloroflexi)、放线菌门(Actinobacteria)和厚壁菌门(Firmicutes)相对丰度总和占细菌菌群的80%~86%, 其中变形菌门(Proteobacteria)丰度最高, 酸根菌门(Acidobacteria)、拟杆菌门(Bacteroidetes)、绿弯菌门(Chloroflexi)、放线菌门(Actinobacteria)和厚壁菌门(Firmicutes)丰度次之[图 3(a)]; 真菌菌群中, 除未分类真菌菌门(unclassified_k_Fungi)外, 子囊菌门(Ascomycota)在各处理中相对丰度最高, 壶菌门(Chytridiomycota)相对丰度次之[图 3(b)].对比发现, 尽管各处理菌门组成具有相似性, 但不同菌门丰度具有明显差异性, 细菌菌门中Bacteroidetes、Actinobacteria、Firmicutes、Saccharibacteria、WWE3、Nitrospirae、Chlamydiae、Latescibacteria、Peregrinibacteria和WS2在不同处理间存在显著差异[P < 0.05, 图 3(a)], 真菌菌门Basidiomycota在不同处理间存在显著差异[P < 0.05, 图 3(b)].

CK: 对照处理; RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理 图 3 各处理土壤细菌及真菌群落组成(门水平) Fig. 3 Bacterial and fungal community composition (phylum level) in the different treatments

微生物α-多样性指数直观地反映出微生物群落多样性和丰富度差异, 不同改良方式对土壤细菌、真菌群落多样性(Shannon指数)和丰富度(Sobs指数)影响如图 4所示.与CK相比, RI和BC_RI处理使得土壤细菌Shannon指数和Sobs指数及真菌Sobs指数显著降低(P < 0.05), BC处理与CK相比无显著变化; 与RI相比, BC_RI使得土壤真菌Shannon指数显著增加(P < 0.05); 与BC相比, BC_RI处理使得土壤细菌和真菌α-多样性指数均显著降低(P < 0.05).

(a)细菌群落α-多样性指数; (b)真菌群落α-多样性指数; CK: 对照处理; RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理; 星号表示差异显著性, *表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 4 各处理细菌及真菌α-多样性指数 Fig. 4 The α-diversity indices of the bacterial and fungal community in the different treatments

2.4 不同改良方式对土壤基质性质及微生物群落结构差异影响

为探究不同改良方式下土壤基质环境(物理结构和性质)及微生物群落结构组成差异, 本研究基于RDA分析, 揭示不同改良方式对各处理土壤基质环境及微生物群落结构差异的影响, 其中土壤物理结构和性质分别通过土壤水稳性团聚体组成和土壤营养元素含量进行表征, 微生物群落结构基于OTU水平进行分析, 分别包含细菌群落及真菌群落.图 5中不同颜色图例代表不同处理样本, 红色箭头代表AM真菌、生物炭及二者联合改良方式, RDA1轴和RDA2轴对结果的解释分别为34.13%~55.68%和6.41%~19.57%, 两轴解释度总和较高, 表明不同改良方式对土壤基质环境影响显著.对比发现, 各处理土壤结构、性质差异明显, AM真菌、生物炭(BC)及二者交互(Interaction)均与不同处理土壤结构和性质呈显著相关关系[P < 0.05, 图 5(a)5(b)].此外, 细菌群落和真菌群落分布分散, 各处理细菌及真菌群落结构差异显著; 同样不同改良方式均与各处理细菌、真菌群落结构呈显著相关关系[P < 0.05, 图 5(c)5(d)].

CK: 对照处理; RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理; AMF: AM真菌; BC: 生物炭; Interaction: AM真菌与生物炭联合改良 图 5 不同改良方式对土壤基质环境(物理结构和性质)及土壤微生物群落结构影响(RDA) Fig. 5 Effects of different amendment treatments on the soil environment (physical structure and properties) and soil microbial community structure (RDA)

2.5 不同改良方式微生物分子生态网络关系分布

分子生态网络关系图显示改良剂施用处理条件下各物种之间的相互关系, 网络图仅显示具有显著相关关系(|r|>0.6, P < 0.001)的物种联系.网络关系图中(图 6图 7), 不同节点代表不同物种, 不同节点颜色代表该物种所属不同菌门, 红色连线表示正相关, 绿色连线表示负相关, 线的多少表示节点之间联系的密切程度.阿拉伯数字表示细菌及真菌群落网络关系中的模块数, 根据模块中物种连接度高低顺序, 将模块依次进行排序, 其中模块1为细菌和真菌群落中物种连接度最高的模块, 该模块被定义为网络关系的核心模块, 该模块中的物种为核心物种[48].

RI: 单独接种RI处理; BC: 单独施加生物炭处理; BC_RI: 接种RI与生物炭联合改良处理; 节点表示不同物种(不同节点颜色代表该物种所属不同菌门), 其中红色名称物种为该网络关系中连接度最高物种(定义为核心物种)、未显示物种名称的节点为未分类(unclassified)菌种; 连线表示物种间存在极显著(|r|>0.6, P < 0.001)相关关系(红色连线代表正相关, 绿色连线表示负相关, 线的多少表示节点之间联系的密切程度的高低); 阿拉伯数字表示不同模块(按照连接度高低顺序依次排序), 下同 图 6 细菌群落单因素相关性网络分析 Fig. 6 Network analysis applied to the bacterial community

图 7 真菌群落单因素相关性网络分析 Fig. 7 Network analysis applied to the fungal community

图 6所示, 不同改良处理条件下, RI处理细菌群落网络关系包含2类模块, BC和BC-RI处理均包含3类模块.对比发现, 不同改良处理中, 物种间相关性程度不同, RI、BC和BC_RI处理细菌群落中正相关关系分别为183、216和262个, 负相关关系分别为198、212和235个, 其中BC和BC_RI处理相关性联系高于RI处理.此外, 不同处理同一模块物种组成不同, 其中核心物种组成具有差异性, RI处理中细菌核心模块主要包含Acidobacteria_bacterium_CB_286306等19种核心物种, BC处理主要包含Domibacillus_robiginosus等23种细菌核心物种, BC_RI处理主要包含Archangium_gephyra等28种细菌核心物种.

图 7所示, 不同改良处理条件下, 真菌群落网络关系模块数、物种间相关性程度及核心物种组成也存在差异性.RI处理真菌群落网络关系包含5类模块, 其中分别包含134和92个正相关和负相关关系, 核心模块主要包含Rhizophagus_intraradices等14种真菌核心物种; BC处理真菌网络关系包含3类模块, 分别存在156和136个正相关和负相关关系, 主要包含Acremonium_tubakiiGibberella_baccata等32种真菌核心物种; BC_RI处理真菌群落网络关系包含4类模块, 其中正相关和负相关关系分别为126和94个, 核心模块中主要包含Gibellulopsis_nigrescens等13种真菌核心物种, 其中Rhizophagus_intraradices并未出现在BC_RI真菌网络关系中.

2.6 AM真菌、生物炭对沙化土壤改良作用机制分析

本研究结果显示, AM真菌和生物炭施用能够影响土壤基质环境(土壤性质、结构和微生物群落结构组成), 不同改良方式对沙化土壤中植物生长具有正向促进作用, 但促进效果具有差异性.考虑到SEM是一种先验方法, 可以用直观的图形表示生态环境中复杂关系网络[49].为探明影响改良效果差异的主要驱动力, 厘清土壤微生物群落及土壤基本性质在改良过程中的重要角色, 本研究通过SEM构建AM真菌、生物炭、土壤性质、土壤结构以及微生物群落关系, 结果如图 8所示.AM真菌分别与细菌群落、真菌群落呈极显著相关关系(r=0.98, P < 0.001; r=0.69, P < 0.001); 生物炭分别与细菌群落、真菌群落和土壤性质呈极显著相关关系(r=-0.17, P < 0.001; r=-0.56, P < 0.001; r=0.93, P < 0.001), 与植物生物量呈显著相关关系(r=0.60, P < 0.05); 此外, 细菌群落、真菌群落分别与植物生物量(r=0.73, P < 0.01; r=-0.44, P < 0.05)、土壤性质(r=-0.88, P < 0.001; r=0.56, P < 0.01)存在显著相关关系.

红色实线箭头表示呈极显著相关关系, 黑色实线箭头表示呈显著相关关系, 虚线表示相关关系不显著; 线条粗细表示相关系数的强弱, 线条越粗, 相关性越强; 箭头数字表示r值; 星号表示差异显著性, *表示P < 0.05, **表示P < 0.01, ***表示P < 0.001 图 8 SEM揭示沙化土壤改良作用机制 Fig. 8 SEM revealed the mechanism of sandy soil improvement

3 讨论

植物生物量能够直接表征植物生长状况, 并间接反映土壤质量与肥力状况的好坏.因此, 植物生长状况对于直观地反映不同改良方式对沙化土壤改良的实际效果具有指导意义.本研究结果发现, AM真菌、生物炭施用以及联合改良对小果白刺生长均有促进作用, 但不同改良方式对植物生长促进效果具有差异性.生物炭施用(BC)及联合改良处理(BC_RI)比AM真菌单独接种(RI)对小果白刺生长促进效果显著, 其中联合改良方式对小果白刺根部生长较单独改良方式促进效果显著.植物生长受土壤基质环境的直接影响, 土壤基质环境变化可能是引起改良效果的差异的主要影响因子.许多研究已证实生物炭施用对提高土壤养分(特别是土壤有机碳)具有重要影响[50~52], 本研究结果同样发现生物炭施用对土壤基质环境有机碳含量具有显著提升效果.由于受沙化土壤养分匮乏等条件限制, 生物炭施用能够为土壤补充直接的养分含量, 这对于植物生长促进效果可能比AM真菌作用更直接.此外, AM真菌、生物炭施用能够影响土壤微生物群落结构[53~55]; 生物炭施用同样也会改变植物菌根真菌的信号传导过程[56], 进而影响菌根真菌协助植物生长的作用效果.生物炭施用后, 直接影响到植物根系与AM真菌定植关系, 因此在短期改良过程中, 联合改良对植物根部生长影响作用更明显.

目前大部分研究已发现AM真菌、生物炭施用对微生物群落组成结构的重要影响[57~59], 但却忽略了不同微生物群落(细菌和真菌群落)在调控土壤基质和改良土壤性质方面的关键作用.本研究通过高通量测序分析, 探究不同改良方式对细菌、真菌群落影响.微生物多样性研究结果显示, AM真菌接种及联合改良能够显著降低微生物群落多样性及丰富度, 这可能主要归因于AM真菌菌根对营养物质竞争的优势作用[60].由于菌根快速生长, 降低了其他微生物营养供给需求, 使得微生物生长受限, 从而降低了土壤微生物多样性.尽管有研究证实生物炭的多孔结构能够为微生物提供一定的生存空间[61], 但由于受到沙化土壤本身条件限制及AM真菌菌根等因素影响, 生物炭对微生物作用优势并不能在短期内得以发挥.因此, 联合改良处理中微生物多样性也显著降低.除微生物多样性差异外, 不同改良方式下土壤物理结构、营养环境及土壤微生物群落结构差异显著, 结合RDA结果发现, AM真菌、生物炭及二者联合改良对土壤基质环境变化具有显著影响, 其中二者联合作用条件下微生物群落结构均与其他改良方式存在明显差异, 表明二者联合可能更易导致微生物群落间互作关系的显著变化.有研究已发现, 生物炭对土壤中菌根群落变化具有显著影响[56], 在菌根群落与生物炭的共同影响下, 使得土壤微生物群落关系发生显著变化, 进而导致微生物群落结构的显著差异.网络分析结果进一步证实了, 不同改良方式对微生物群落互作关系的显著影响.结果表明, 不同改良方式下微生物分子生态网络关系能够发生显著变化, 不同处理网络关系中核心物种组成也发生明显改变; 生物炭施加处理(BC和BC_RI)较单独接种AM真菌改良处理(RI), 网络连接密集程度更高且核心物种组成更丰富; 生物炭与AM真菌联合, 弱化了Rhizophagus intraradices的核心作用, 并增强其他微生物(特别是细菌菌种)的核心地位.这些结果表明, 不同改良方式对微生物菌群互作关系具有重要影响, 生物炭施加可能能够平衡微生物群落对生态位和营养等的合理利用, 这对于有效调控微生物群落间协同分工意义重大.AM真菌单独作用时, 由于AM真菌菌根的优势作用, 使得其他微生物生长受限; 当施加生物炭后, 一方面由于生物炭本身多孔结构, 能够为其他微生物在其孔隙形成活性生物膜提供合适的生态位[61]; 另一方面, 生物炭能够通过吸附植物根系分泌物等方式改变根系周围环境, 进而影响AM真菌菌根生长[62].生物炭施加能够适当减弱AM真菌对土壤营养、环境等的争夺优势, 从而为其他更多的有益微生物生长提供基础.微生物群落协同联系对土壤生态恢复具有重要影响[63], 尽管AM真菌、生物炭和二者联合对沙化土壤改良均有正向促进作用, 但生物炭施加更有利于增强微生物群落互作关系, 因此, 生物炭对促进沙化土壤改良效果影响更显著.此外, AM真菌与生物炭联合, 更有利于平衡微生物群落间营养供给需求, 这可能对实现沙化土壤生态恢复具有潜在优势.

改良剂施用对土壤微环境变化有着直接影响, 为探明微生物群落、土壤基质环境在沙化土壤改良过程中的影响作用, 本研究通过SEM直观地揭示不同改良方式对沙化土壤改良的作用机制.结果表明, AM真菌和生物炭施用通过直接影响土壤细菌和真菌群落结构, 进而实现对植物生长及土壤性质变化的影响; 生物炭施用除能够直接影响土壤微生物群落外, 也会直接作用于土壤以及植物本身; 因此, 生物炭施用较AM真菌对沙化土壤改良优势更显著.综上结果发现, 微生物群落在土壤改良过程中起着关键性作用, 微生物群落结构差异(特别是微生物间的互作关系变化)是导致土壤改良效果差异的主要驱动力.由于生物炭施用能够瞬时弥补土壤养分匮乏的限制, 因此就短期改良效果来看, 生物炭较AM真菌对沙化土壤改良优势更显著.尽管AM真菌与生物炭联合施用对沙化土壤短期改良效果并不一定完全优于单独施加生物炭改良作用效果, 但本研究结果已显示出联合改良对植物根部生长的显著促进作用, 说明AM真菌-生物炭-植物联合修复对加速沙化土壤改良的潜在优势.然而想要构建完善高效的植物-微生物-改良剂联合改良体系, 除了需要考虑土壤基质、改良剂本身性质影响外, 还需要关注微生物群落, 特别是有益微生物(如菌根有益细菌、真菌等)在其中发挥的关键作用.日后可以考虑寻找更多的有益微生物, 用以加强联合改良的协同过程, 并针对不同土壤基质, 通过野外实地试验探讨其在退化土壤生态恢复中的实际作用效果, 为后续理解并提高土壤改良修复效率提供有力依据.

4 结论

(1) AM真菌、生物炭及联合改良方式均对沙化土壤中小果白刺地上部生物量及矿质营养元素吸收具有显著促进作用; 生物炭和联合改良处理促进了小果白刺根部矿质营养元素吸收和生长, 联合改良处理效果要优于生物炭处理.

(2) 不同改良方式使得沙化土壤营养环境、物理结构具有显著差异性, AM真菌、生物炭及联合改良方式对土壤基质环境变化存在显著影响.

(3) 不同改良方式使得沙化土壤细菌、真菌群落α-多样性、组成及结构存在显著差异性, AM真菌、生物炭及二者联合对土壤微生物群落结构影响显著.

(4) 不同改良方式下微生物分子生态网络关系发生显著变化, 不同处理网络关系中核心物种组成具有差异性; 生物炭和联合改良处理较单独接种AM真菌改良处理, 网络连接密集程度更高且核心物种组成更丰富; 生物炭与AM真菌联合, 弱化了Rhizophagus intraradices的核心作用, 并增强其他微生物(特别是细菌菌种)的核心地位.

(5) AM真菌、生物炭施用通过直接影响土壤细菌、真菌群落结构, 进而实现对植物生长及土壤性质变化的影响; 微生物群落在土壤改良过程中具有关键性作用, 微生物群落结构差异(特别是微生物间的互作关系变化)是导致土壤改良效果差异的主要驱动力.

参考文献
[1] Verón S R, Paruelo J M, Oesterheld M. Assessing desertification[J]. Journal of Arid Environments, 2006, 66(4): 751-763. DOI:10.1016/j.jaridenv.2006.01.021
[2] 韩洋. 荒漠(沙)化土地治理现状及对策研究[J]. 山西农经, 2019(23): 82, 84.
[3] 刘佳, 刘远妹, 杜忠. 中国退化沙化草地治理研究进展[J]. 安徽农学通报, 2018, 24(21): 161-164.
Liu J, Liu Y M, Du Z. Review on the research of degraded and desertification grassland restoration in China[J]. Anhui Agricultural Science Bulletin, 2018, 24(21): 161-164. DOI:10.3969/j.issn.1007-7731.2018.21.061
[4] Sial T A, Liu J, Zhao Y, et al. Co-application of milk tea waste and NPK fertilizers to improve sandy soil biochemical properties and wheat growth[J]. Molecules, 2019, 24(3). DOI:10.3390/molecules24030423
[5] 王桂君, Ohsowoski B, Klironomos J, 等. 菌根菌剂及土壤改良剂对退化生态系统的修复潜能分析[J]. 生态经济, 2014, 30(7): 179-184.
Wang G J, Ohsowoski B, Klironomos J, et al. The possibility of fungi inoculums and soil amendment for restoring degraded ecosystems[J]. Ecological Economy, 2014, 30(7): 179-184. DOI:10.3969/j.issn.1671-4407.2014.07.043
[6] 丁黎, 吴胜春, 单胜道, 等. 沙化土壤施用羊粪生物质炭及堆肥的养分渗流特征[J]. 科技通报, 2019, 35(3): 61-69, 74.
Ding L, Wu S C, Shan S D, et al. The nutrients leaching characteristics of desertification soil by the application of sheep droppings biochar and fertilizer[J]. Bulletin of Science and Technology, 2019, 35(3): 61-69, 74.
[7] Zhang L Y, Jing Y M, Chen G H, et al. Improvement of physical and hydraulic properties of desert soil with amendment of different biochars[J]. Journal of Soils and Sediments, 2019, 19(7): 2984-2996. DOI:10.1007/s11368-019-02293-8
[8] 魏秋兰. 浅谈生物炭对土壤改良影响的研究进展[J]. 南方农业, 2019, 13(30): 189-190.
[9] 张千丰, 王光华. 生物炭理化性质及对土壤改良效果的研究进展[J]. 土壤与作物, 2012, 1(4): 219-226.
Zhang Q F, Wang G H. Research progress of physiochemical properties of biochar and its effects as soil amendments[J]. Soil and Crop, 2012, 1(4): 219-226.
[10] Tan Z X, Lin C S K, Ji X Y, et al. Returning biochar to fields: a review[J]. Applied Soil Ecology, 2017, 116: 1-11. DOI:10.1016/j.apsoil.2017.03.017
[11] Ogura T, Date Y, Masukujane M, et al. Improvement of physical, chemical and biological properties of aridisol from Botswana by the incorporation of torrefied biomass[J]. Scientific Reports, 2016, 6. DOI:10.1038/srep28011
[12] Blanco-Canqui H. Biochar and soil physical properties[J]. Soil Science Society of America Journal, 2017, 81(4): 687-711. DOI:10.2136/sssaj2017.01.0017
[13] Wang H, Huang Y F, Shen C Y, et al. Co-transport of pesticide acetamiprid and silica nanoparticles in biochar-amended sand porous media[J]. Journal of Environmental Quality, 2016, 45(5): 1749-1759. DOI:10.2134/jeq2016.02.0073
[14] Trifunovic B, Gonzales H B, Ravi S, et al. Dynamic effects of biochar concentration and particle size on hydraulic properties of sand[J]. Land Degradation & Development, 2018, 29(4): 884-893.
[15] Głąb T, Palmowska J, Zaleski T, et al. Effect of biochar application on soil hydrological properties and physical quality of sandy soil[J]. Geoderma, 2016, 281: 11-20. DOI:10.1016/j.geoderma.2016.06.028
[16] Khan W U D, Ramzani P M A, Anjum S, et al. Potential of miscanthus biochar to improve sandy soil health, in situ nickel immobilization in soil and nutritional quality of spinach[J]. Chemosphere, 2017, 185: 1144-1156. DOI:10.1016/j.chemosphere.2017.07.097
[17] Carvalho M L, De Moraes M T, Cerri C E P, et al. Biochar amendment enhances water retention in a tropical sandy soil[J]. Agriculture, 2020, 10(3). DOI:10.3390/agriculture10030062
[18] Hardie M, Clothier B, Bound S, et al. Does biochar influence soil physical properties and soil water availability?[J]. Plant and Soil, 2014, 376(1-2): 347-361. DOI:10.1007/s11104-013-1980-x
[19] Jeffery S, Meinders M B J, Stoof C R, et al. Biochar application does not improve the soil hydrological function of a sandy soil[J]. Geoderma, 2015, 251-252: 47-54. DOI:10.1016/j.geoderma.2015.03.022
[20] Rillig M C, Mardatin N F, Leifheit E F, et al. Mycelium of arbuscular mycorrhizal fungi increases soil water repellency and is sufficient to maintain water-stable soil aggregates[J]. Soil Biology and Biochemistry, 2010, 42(7): 1189-1191. DOI:10.1016/j.soilbio.2010.03.027
[21] Singh S, Kapoor K K. Inoculation with phosphate-solubilizing microorganisms and a vesicular-arbuscular mycorrhizal fungus improves dry matter yield and nutrient uptake by wheat grown in a sandy soil[J]. Biology and Fertility of Soils, 1999, 28(2): 139-144.
[22] Medina A, Azcón R. Effectiveness of the application of arbuscular mycorrhiza fungi and organic amendments to improve soil quality and plant performance under stress conditions[J]. Journal of Soil Science and Plant Nutrition, 2010, 10(3): 354-372.
[23] Acosta-Martínez V, Cotton J. Lasting effects of soil health improvements with management changes in cotton-based cropping systems in a sandy soil[J]. Biology and Fertility of Soils, 2017, 53(5): 533-546. DOI:10.1007/s00374-017-1192-2
[24] Song Z H, Bi Y L, Zhang J A, et al. Arbuscular mycorrhizal fungi promote the growth of plants in the mining associated clay[J]. Scientific Reports, 2020, 10(1). DOI:10.1038/s41598-020-59447-9
[25] Li S P, Bi Y L, Kong W P, et al. Effects of arbuscular mycorrhizal fungi on ecological restoration in coal mining areas[J]. Russian Journal of Ecology, 2015, 46(5): 431-437. DOI:10.1134/S1067413615050173
[26] Zhang Y F, Wang P, Yang Y F, et al. Arbuscular mycorrhizal fungi improve reestablishment of Leymus chinensis in bare saline-alkaline soil: implication on vegetation restoration of extremely degraded land[J]. Journal of Arid Environments, 2011, 75(9): 773-778. DOI:10.1016/j.jaridenv.2011.04.008
[27] DeLuca T H, MacKenzie M D, Gundale M J, et al. Wildfire-produced charcoal directly influences nitrogen cycling in ponderosa pine forests[J]. Soil Science Society of America Journal, 2006, 70(2): 448-453. DOI:10.2136/sssaj2005.0096
[28] Matsubara Y, Hasegawa N, Fukui H. Incidence of fusarium root rot in asparagus seedlings infected with arbuscular mycorrhizal fungus as affected by several soil amendments[J]. Engei Gakkai Zasshi, 2002, 71(3): 370-374. DOI:10.2503/jjshs.71.370
[29] Verma B, Reddy M S. Biochar augmentation improves ectomycorrhizal colonisation, plant growth and soil fertility[J]. Soil Research, 2020, 58(7): 673-682. DOI:10.1071/SR20067
[30] Ishii T, Kadoya K. Effects of charcoal as a soil conditioner on citrus growth and vesicular-arbuscular mycorrhizal development[J]. Journal of the Japanese Society for Horticultural Science, 1994, 63(3): 529-535. DOI:10.2503/jjshs.63.529
[31] 潘丽冰, 徐凡珍, 沙丽清. 生物炭对土壤理化性质及橡胶幼苗生物量的影响[J]. 山地学报, 2015, 33(4): 449-456.
[32] Zhang T A, Chen H Y H, Ruan H H. Global negative effects of nitrogen deposition on soil microbes[J]. The ISME Journal, 2018, 12(7): 1817-1825. DOI:10.1038/s41396-018-0096-y
[33] Liu Y, Guo K Y, Zhao Y Z, et al. Change in composition and function of microbial communities in an acid bamboo (Phyllostachys praecox) plantation soil with the addition of three different biochars[J]. Forest Ecology and Management, 2020, 473. DOI:10.1016/j.foreco.2020.118336
[34] Yang Q, Ravnskov S, Neumann Andersen M. Nutrient uptake and growth of potato: Arbuscular mycorrhiza symbiosis interacts with quality and quantity of amended biochars[J]. Journal of Plant Nutrition and Soil Science, 2020, 183(2): 220-232. DOI:10.1002/jpln.201900205
[35] Zhou H, Zhao W Z, Zheng X J, et al. Root distribution of Nitraria sibirica with seasonally varying water sources in a desert habitat[J]. Journal of Plant Research, 2015, 128(4): 613-622. DOI:10.1007/s10265-015-0728-5
[36] 刘利平, 斯琴毕力格, 徐志明, 等. 阿拉善荒漠区白刺属果实营养成分及其利用分析[J]. 内蒙古林业科技, 2016, 42(3): 29-31, 42.
Liu L P, Siqinbilige, Xu Z M, et al. Nutritional components of fruits of Nitraria L. and its utilization in Alashan desert area[J]. Journal of Inner Mongolia Forestry Science & Technology, 2016, 42(3): 29-31, 42.
[37] 王芳, 郭伟, 马朋坤, 等. 丛枝菌根真菌对铈污染土壤上玉米生长和铈吸收的影响[J]. 环境科学, 2016, 37(1): 309-316.
Wang F, Guo W, Ma P K, et al. Effects of arbuscular mycorrhizal fungi on the growth and Ce uptake of maize growth in Ce-contaminated soils[J]. Environmental Science, 2016, 37(1): 309-316.
[38] 陈萌萌, 唐东山, 张晓文, 等. 生物炭对荒漠土壤的改良作用研究[J]. 山东化工, 2019, 48(7): 217-221, 224.
Chen M M, Tang D S, Zhang X W, et al. Research on the improvement of biochar on desert Soil[J]. Shandong Chemical Industry, 2019, 48(7): 217-221, 224.
[39] 王道涵, 邢聪慧, 梅傲雪, 等. 生物炭改良风沙土对土壤团聚体稳定性的影响[J]. 辽宁工程技术大学学报(自然科学版), 2018, 37(5): 816-823.
Wang D H, Xing C H, Mei A X, et al. Effect of biochar on the stability of aggregates in Aeolian soil[J]. Journal of Liaoning Technical University (Natural Science), 2018, 37(5): 816-823.
[40] 王桂君, 许振文, 路倩倩. 生物炭对沙化土壤理化性质及作物幼苗的影响[J]. 江苏农业科学, 2017, 45(11): 246-248.
[41] 肖雪毅, 陈保冬, 朱永官. 丛枝菌根真菌对铜尾矿上植物生长和矿质营养的影响[J]. 环境科学学报, 2006, 26(2): 312-317.
Xiao X Y, Chen B D, Zhu Y G. The influences of arbuscular mycorrhizal fungi on growth and mineral nutrition of plants grown in copper mine tailing[J]. Acta Scientiae Circumstantiae, 2006, 26(2): 312-317.
[42] 常青, 郭伟, 潘亮, 等. 镧-铅复合污染下AM真菌对玉米生长和镧、铅吸收的影响[J]. 环境科学, 2017, 38(9): 3915-3926.
Chang Q, Guo W, Pan L, et al. Effects of arbuscular mycorrhizal fungi on the growth and uptake of La and Pb by maize grown in La and Pb-contaminated soil[J]. Environmental Science, 2017, 38(9): 3915-3926.
[43] 王巧环, 任玉芬, 孟龄, 等. 元素分析仪同时测定土壤中全氮和有机碳[J]. 分析试验室, 2013, 32(10): 41-45.
[44] Pisani O, Liebert D, Bosch D D, et al. Element losses from fields in conventional and conservation tillage in the Atlantic Coastal Plain, Georgia, United States[J]. Journal of Soil and Water Conservation, 2020, 75(3): 376-386.
[45] Mohamed B, Mounia K, Aziz A, et al. Sewage sludge used as organic manure in Moroccan sunflower culture: effects on certain soil properties, growth and yield components[J]. Science of the Total Environment, 2018, 627: 681-688.
[46] 吕贻忠, 李保国. 土壤学[M]. 北京: 中国农业出版社, 2006.
Lü Y Z, Li B G. Soil science[M]. Beijing: China Agriculture Press, 2006.
[47] Fang D X, Zhao G, Xu X Y, et al. Microbial community structures and functions of wastewater treatment systems in plateau and cold regions[J]. Bioresource Technology, 2018, 249: 684-693.
[48] 马静, 卢永强, 张琦, 等. 黄土高原采煤沉陷对土壤微生物群落的影响[J]. 土壤学报, 2020.
Ma J, Lu Y Q, Zhang Q, et al. Effects of coal mining subsidence on soil microbial community in the Loess Plateau[J]. Acta Pedologica Sinica, 2020. DOI:10.11766/trxb202003160122
[49] Wasilah S, Fahmyddin T. The advancement of the built environment research through employment of structural equation modeling (SEM)[J]. IOP Conference Series: Earth and Environmental Science, 2018, 126. DOI:10.1088/1755-1315/126/1/012001
[50] Van Zwieten L, Kimber S, Morris S, et al. Effects of biochar from slow pyrolysis of papermill waste on agronomic performance and soil fertility[J]. Plant and Soil, 2010, 327(1-2): 235-246.
[51] 周志红, 李心清, 邢英, 等. 生物炭对土壤氮素淋失的抑制作用[J]. 地球与环境, 2011, 39(2): 278-284.
[52] Lehmann J, Kleber M. The contentious nature of soil organic matter[J]. Nature, 2015, 528(7580): 60-68.
[53] Fan X X, Chang W, Sui X, et al. Changes in rhizobacterial community mediating atrazine dissipation by arbuscular mycorrhiza[J]. Chemosphere, 2020, 256. DOI:10.1016/j.chemosphere.2020.127046
[54] Lopes L D, Fernandes M F. Changes in microbial community structure and physiological profile in a kaolinitic tropical soil under different conservation agricultural practices[J]. Applied Soil Ecology, 2020, 152. DOI:10.1016/j.apsoil.2020.103545
[55] Sarfraz R, Yang W H, Wang S S, et al. Short term effects of biochar with different particle sizes on phosphorous availability and microbial communities[J]. Chemosphere, 2020, 256. DOI:10.1016/j.chemosphere.2020.126862
[56] Warnock D D, Lehmann J, Kuyper T W, et al. Mycorrhizal responses to biochar in soil-concepts and mechanisms[J]. Plant and Soil, 2007, 300(1-2): 9-20.
[57] Zhu X M, Chen B L, Zhu L Z, et al. Effects and mechanisms of biochar-microbe interactions in soil improvement and pollution remediation: a review[J]. Environmental Pollution, 2017, 227: 98-115.
[58] Kolton M, Graber E R, Tsehansky L, et al. Biochar-stimulated plant performance is strongly linked to microbial diversity and metabolic potential in the rhizosphere[J]. New Phytologist, 2017, 213(3): 1393-1404.
[59] Nuccio E E, Hodge A, Pett-Ridge J, et al. An arbuscular mycorrhizal fungus significantly modifies the soil bacterial community and nitrogen cycling during litter decomposition[J]. Environmental Microbiology, 2013, 15(6): 1870-1881.
[60] Jacoby R, Peukert M, Succurro A, et al. The role of soil microorganisms in plant mineral nutrition-current knowledge and future directions[J]. Frontiers in Plant Science, 2017, 8. DOI:10.3389/fpls.2017.01617
[61] 崔丙健, 崔二苹, 胡超, 等. 生物质炭施用对再生水灌溉空心菜根际微生物群落结构及多样性的影响[J]. 环境科学, 2020, 41(12).
Cui B J, Cui E P, Hu C, et al. Effects of selected biochars application on the microbial community structures and diversities in rhizosphere of water Spinach (Ipomoea aquatica Forssk.) irrigated with reclaimed water[J]. Environmental Science, 2020, 41(12). DOI:10.13227/j.hjkx.202006087
[62] Gu Y A, Hou Y G, Huang D P, et al. Application of biochar reduces Ralstonia solanacearum infection via effects on pathogen chemotaxis, swarming motility, and root exudate adsorption[J]. Plant and Soil, 2017, 415(1-2): 269-281. DOI:10.1007/s11104-016-3159-8
[63] De Vries F T, Griffiths R I, Bailey M, et al. Soil bacterial networks are less stable under drought than fungal networks[J]. Nature Communications, 2018, 9. DOI:10.1038/s41467-018-05516-7