生物质炭-沼液联合施用对调控氮循环功能基因促进氮素增效的影响

唐存柳, 靳泽文, 姚光伟, 平立凤, 虞轶俊, 郑一玲, 范栩妙, 单胜道. 生物质炭-沼液联合施用对调控氮循环功能基因促进氮素增效的影响[J]. 环境化学, 2023, 42(8): 2843-2852. doi: 10.7524/j.issn.0254-6108.2022030601
引用本文: 唐存柳, 靳泽文, 姚光伟, 平立凤, 虞轶俊, 郑一玲, 范栩妙, 单胜道. 生物质炭-沼液联合施用对调控氮循环功能基因促进氮素增效的影响[J]. 环境化学, 2023, 42(8): 2843-2852. doi: 10.7524/j.issn.0254-6108.2022030601
TANG Cunliu, JIN Zewen, YAO Guangwei, PING Lifeng, YU Yijun, ZHENG Yiling, FAN Xumiao, SHAN Shengdao. Effects of combined application of biochar and biogas slurry on regulating nitrogen cycle function genes and promoting nitrogen efficiency[J]. Environmental Chemistry, 2023, 42(8): 2843-2852. doi: 10.7524/j.issn.0254-6108.2022030601
Citation: TANG Cunliu, JIN Zewen, YAO Guangwei, PING Lifeng, YU Yijun, ZHENG Yiling, FAN Xumiao, SHAN Shengdao. Effects of combined application of biochar and biogas slurry on regulating nitrogen cycle function genes and promoting nitrogen efficiency[J]. Environmental Chemistry, 2023, 42(8): 2843-2852. doi: 10.7524/j.issn.0254-6108.2022030601

生物质炭-沼液联合施用对调控氮循环功能基因促进氮素增效的影响

    通讯作者: E-mail:119044@zust.edu.cn
  • 基金项目:
    国家重点研发计划项目“数据和知识联合驱动城市易腐垃圾炭化与绿色可持续利用的关键技术及应用”(2022YFE0196000)和浙江省重点研发项目(2019C02061,2020C01017)资助.

Effects of combined application of biochar and biogas slurry on regulating nitrogen cycle function genes and promoting nitrogen efficiency

    Corresponding author: JIN Zewen, 119044@zust.edu.cn
  • Fund Project: Key Technoiogies of Carbonization and Sustainable Utilization of Urban Perishable Waste Driven by Date and Knowledge(2022YFE0196000) and Key Research and Development Program of Zhejiang Province (2019C02061,2020C01017) .
  • 摘要: 为探讨生物质炭-沼液配施条件下氮循环功能基因调控农田土壤氮素转化并影响农作物氮素吸收利用机制. 本试验以浙江省杭州市红黄壤作为研究对象,设置生物质炭和沼液两个因素,探究生物质炭-沼液配施条件下土壤基本理化性质和氮循环功能基因丰度变化情况,刻画功能基因与农田氮素利用率间的耦合关系. 结果表明,生物质炭-沼液配施可以显著降低土壤容重,提升土壤pH和土壤氮素含量,其中,高剂量生物质炭-沼液配施(C3B2)处理较单施化肥(C0B0)处理铵态氮、硝态氮、全氮含量增幅均达到显著水平(P<0.05). 与空白处理(CK)相比,生物质炭-沼液配施(C3B2)处理则显著提高了反硝化功能基因丰度,较单施化肥(C0B0)处理增幅30.98%和44.99%. 冗余分析结果显示,铵态氮、硝态氮和有机碳含量对土壤氮循环功能基因影响较为显著,结构方程模型则表明硝化作用功能基因丰度的提升对包菜氮素农学利用率呈现负相关趋势. 研究结果表明,在相同养分施用量的条件下,生物质炭-沼液配施可显著提高土壤肥力. 氮素和有机碳含量是影响功能基因丰度的关键因素,硝化作用功能基因丰度的降低可以提高农田氮素利用率. 本研究结果可以为促进农业废弃物资源化利用,提升农田土壤氮素活性提供理论依据.
  • 加载中
  • 图 1  不同生物质炭-沼液处理下土壤氮循环功能基因丰度

    Figure 1.  Abundance of functional genes of soil nitrogen cycling under different biochar and biogas slurry treatments

    图 2  不同处理包菜产量与氮素利用率

    Figure 2.  Production and nitrogen use efficiency in different treatments

    图 3  功能基因与环境因子冗余分析(RDA)

    Figure 3.  Redundancy analysis(RDA) of functional genes and environmental factors

    图 4  氮循环功能基因与氮素利用率路径分析图

    Figure 4.  Path analysis of nitrogen cycling functional genes and nitrogen use efficiency

    表 1  试验设计方案

    Table 1.  Test design scheme

    处理
    Treatment
    生物质炭a/%
    Biocher
    沼液b/%
    Biogas slurry
    处理
    Treatment
    生物质炭/%
    Biocher
    沼液/%
    Biogas slurry
    处理
    Treatment
    生物质炭/%
    Biocher
    沼液/%
    Biogas slurry
    C0B000C0B1050C0B20100
    C1B00.50C1B10.550C1B20.5100
    C2B010C2B1150C2B21100
    C3B020C3B1250C3B22100
      a代表生物质炭占小区耕层土壤总重比,b代表所施沼液全氮量占施氮总量比. 根据本项目小区面积和沼液全氮含量,0%生物质炭施用量为0 kg,0.5%生物质炭施用量为13.5 kg,1%生物质炭施用量为27 kg,2%生物质炭施用量为54 kg,0%沼液施用量为0 L,50%沼液施用量为22.7 L,100%沼液施用量为45.4 L.
      The percentage of biocher application represents the percentage of the total weight of residential tillage soil; the percentage of biogas slurry application represents the percentage of total biogas slurry applied for the total nitrogen applied. According to the residential area and total nitrogen content of biogas slurry, The application amount of 0% biochar was 0 kg, The application amount of 0.5% biochar was 13.5 kg, The application amount of 1% biochar was 27 kg, The application amount of 2% biochar was 54 kg, 0% biogas slurry application amount was 0 L, 50% biogas slurry application amount was 22.7 L and 100% biogas slurry application amount was 45.4 L.
    处理
    Treatment
    生物质炭a/%
    Biocher
    沼液b/%
    Biogas slurry
    处理
    Treatment
    生物质炭/%
    Biocher
    沼液/%
    Biogas slurry
    处理
    Treatment
    生物质炭/%
    Biocher
    沼液/%
    Biogas slurry
    C0B000C0B1050C0B20100
    C1B00.50C1B10.550C1B20.5100
    C2B010C2B1150C2B21100
    C3B020C3B1250C3B22100
      a代表生物质炭占小区耕层土壤总重比,b代表所施沼液全氮量占施氮总量比. 根据本项目小区面积和沼液全氮含量,0%生物质炭施用量为0 kg,0.5%生物质炭施用量为13.5 kg,1%生物质炭施用量为27 kg,2%生物质炭施用量为54 kg,0%沼液施用量为0 L,50%沼液施用量为22.7 L,100%沼液施用量为45.4 L.
      The percentage of biocher application represents the percentage of the total weight of residential tillage soil; the percentage of biogas slurry application represents the percentage of total biogas slurry applied for the total nitrogen applied. According to the residential area and total nitrogen content of biogas slurry, The application amount of 0% biochar was 0 kg, The application amount of 0.5% biochar was 13.5 kg, The application amount of 1% biochar was 27 kg, The application amount of 2% biochar was 54 kg, 0% biogas slurry application amount was 0 L, 50% biogas slurry application amount was 22.7 L and 100% biogas slurry application amount was 45.4 L.
    下载: 导出CSV

    表 2  功能基因PCR扩增引物

    Table 2.  Primers for PCR amplification of functional genes

    目标基因
    Target gene
    引物
    Primer
    引物序列(5'-3')
    Primer sequences
    amoA-AArch-amoAF
    Arch-amoAR
    STAATGGTCTGGCTTAGACG
    GCGGCCATCCATCTGTATGT
    amoA-BamoA-1F
    amoA-2R
    GGGGTTTCTACTGGTGGT
    CCCCTCKGSAAAGCCTTCTTC
    nxrAF1370F
    F2843R
    CAGACCGACGTGTGCGAAAG
    TCCACAAGGAACGGAAGGTC
    narGnarG 1960F
    narG 2650R
    TAYGTSGGCCARGARAA
    TTYTCRTACCABGTBGC
    nirSCd3Af
    R3cd
    GTSAACGTSAAGGARACSGG
    GASTTCGGRTGSGTCTTGA
    nirKF1aCu
    R3Cu
    ATCATGGTSCTGCCGCG
    GCCTCGATCAGRTTGTGGTT
    nosZnosZF
    nosZ1662R
    CGYTGTTCMTCGACAGCCAG
    CGSACCTTSTTGCCSTYGCG
    napAnapAV67F
    napAV67R
    TAYTTYYTNHSNAARATHATGTAYGG
    DATNGGRTGCATYTCNGCCATRTT
    目标基因
    Target gene
    引物
    Primer
    引物序列(5'-3')
    Primer sequences
    amoA-AArch-amoAF
    Arch-amoAR
    STAATGGTCTGGCTTAGACG
    GCGGCCATCCATCTGTATGT
    amoA-BamoA-1F
    amoA-2R
    GGGGTTTCTACTGGTGGT
    CCCCTCKGSAAAGCCTTCTTC
    nxrAF1370F
    F2843R
    CAGACCGACGTGTGCGAAAG
    TCCACAAGGAACGGAAGGTC
    narGnarG 1960F
    narG 2650R
    TAYGTSGGCCARGARAA
    TTYTCRTACCABGTBGC
    nirSCd3Af
    R3cd
    GTSAACGTSAAGGARACSGG
    GASTTCGGRTGSGTCTTGA
    nirKF1aCu
    R3Cu
    ATCATGGTSCTGCCGCG
    GCCTCGATCAGRTTGTGGTT
    nosZnosZF
    nosZ1662R
    CGYTGTTCMTCGACAGCCAG
    CGSACCTTSTTGCCSTYGCG
    napAnapAV67F
    napAV67R
    TAYTTYYTNHSNAARATHATGTAYGG
    DATNGGRTGCATYTCNGCCATRTT
    下载: 导出CSV

    表 3  不同生物质炭和沼液处理下土壤理化性质

    Table 3.  The physicochemical properties of the soil under different biocher and biogas slurry treatments

    处理
    Trentment
    pH容重/(g·cm-3
    SBD
    硝态氮/(mg·kg−1
    NO3--N
    铵态氮/(mg·kg−1
    NH4+-N
    全氮/(g·kg−1
    Total N
    有机碳/(g·kg−1
    Organic Carbon
    B0CK5.31±0.16c1.45±0.01a3.54±0.74c0.79±0.10d0.6±0.01c4.82±0.97d
    C05.15±0.17c1.43±0.01ab4.13±1.37c2.37±0.52c0.5±0.01c4.33±0.64d
    C16.28±2.8b1.25±0.01c8.68±1.33b3.40±0.58c1.4±0.03b11.48±1.18c
    C27.15±0.34b1.41±0.02ab14.78±2.07a7.44±0.97b1.9±0.06b15.87±2.08b
    C37.49±0.1a1.36±0.01b15.02±1.25a15.72±1.13a4.5±0.04a22.44±0.90a
    B1CK5.31±0.09c1.45±0.01a3.54±0.23c0.79±0.10d0.6±0.02c4.82±0.97d
    C05.31±0.07c1.49±0.01a4.58±0.07c3.39±0.10c0.5±0.01c4.34±0.33d
    C17.14±0.20b1.28±0.02b3.51±0.02c4.11±0.24c1.3±0.00b11.92±1.15c
    C27.35±0.06ab1.43±0.06b16.24±5.75b11.16±1.46b2.3±0.006a17.78±1.05a
    C37.39±0.23a1.42±0.06b28.34±1.19a18.66±0.61a2.2±0.016a14.22±0.59b
    B2CK5.31±0.09b1.45±0.01ab3.54±0.72d0.79±0.10d0.6±0.02d4.82±0.97d
    C05.17±0.09b1.47±0.02a4.37±0.93d4.04±0.81c0.9±0.06cd4.16±0.19d
    C17.38±0.04a1.29±0.01c8.02±0.86c4.88±0.46c1.7±0.05bc10.03±1.77c
    C27.37±0.16a1.45±0.04ab15±0.51b13.54±1.22b2.7±0.04b13.71±1.20b
    C37.46±0.05a1.36±0.02bc27.9±1.8a19.46±1.00a5.2±0.09a23.05±0.28a
      表中不同处理同列数据后不同小写字母表示各处理间差异显著(P<0.05)
      Different lowercase letters after different processing of the same column in the table indicate significant differences between treatments (P<0.05)
    处理
    Trentment
    pH容重/(g·cm-3
    SBD
    硝态氮/(mg·kg−1
    NO3--N
    铵态氮/(mg·kg−1
    NH4+-N
    全氮/(g·kg−1
    Total N
    有机碳/(g·kg−1
    Organic Carbon
    B0CK5.31±0.16c1.45±0.01a3.54±0.74c0.79±0.10d0.6±0.01c4.82±0.97d
    C05.15±0.17c1.43±0.01ab4.13±1.37c2.37±0.52c0.5±0.01c4.33±0.64d
    C16.28±2.8b1.25±0.01c8.68±1.33b3.40±0.58c1.4±0.03b11.48±1.18c
    C27.15±0.34b1.41±0.02ab14.78±2.07a7.44±0.97b1.9±0.06b15.87±2.08b
    C37.49±0.1a1.36±0.01b15.02±1.25a15.72±1.13a4.5±0.04a22.44±0.90a
    B1CK5.31±0.09c1.45±0.01a3.54±0.23c0.79±0.10d0.6±0.02c4.82±0.97d
    C05.31±0.07c1.49±0.01a4.58±0.07c3.39±0.10c0.5±0.01c4.34±0.33d
    C17.14±0.20b1.28±0.02b3.51±0.02c4.11±0.24c1.3±0.00b11.92±1.15c
    C27.35±0.06ab1.43±0.06b16.24±5.75b11.16±1.46b2.3±0.006a17.78±1.05a
    C37.39±0.23a1.42±0.06b28.34±1.19a18.66±0.61a2.2±0.016a14.22±0.59b
    B2CK5.31±0.09b1.45±0.01ab3.54±0.72d0.79±0.10d0.6±0.02d4.82±0.97d
    C05.17±0.09b1.47±0.02a4.37±0.93d4.04±0.81c0.9±0.06cd4.16±0.19d
    C17.38±0.04a1.29±0.01c8.02±0.86c4.88±0.46c1.7±0.05bc10.03±1.77c
    C27.37±0.16a1.45±0.04ab15±0.51b13.54±1.22b2.7±0.04b13.71±1.20b
    C37.46±0.05a1.36±0.02bc27.9±1.8a19.46±1.00a5.2±0.09a23.05±0.28a
      表中不同处理同列数据后不同小写字母表示各处理间差异显著(P<0.05)
      Different lowercase letters after different processing of the same column in the table indicate significant differences between treatments (P<0.05)
    下载: 导出CSV
  • [1] 董达, 王宇婕, 姜培坤, 等. 炭基肥和竹炭对土壤氮素淋失和微生物的影响 [J]. 水土保持学报, 2021, 35(2): 144-151. doi: 10.13870/j.cnki.stbcxb.2021.02.020

    DONG D, WANG Y J, JIANG P K, et al. Effects of carbon based fertilizer and bamboo charcoal on soil nitrogen loss and microorganism [J]. Journal of Soil and Water Conservation, 2021, 35(2): 144-151(in Chinese). doi: 10.13870/j.cnki.stbcxb.2021.02.020

    [2] 虞璐. 生物质炭对酸化土壤的改良效应及其对土壤硝化作用的影响[D]. 杭州: 浙江大学, 2019.

    YU L. The effect of biochar on acidic soil amelioration and soil nitrification[D]. Hangzhou: Zhejiang University, 2019(in Chinese).

    [3] 孟阳阳, 刘冰, 康建军. 不同水肥措施对甜高粱农田土壤呼吸和硝化-反硝化作用的影响 [J]. 甘肃农业大学学报, 2020, 55(6): 67-76. doi: 10.13432/j.cnki.jgsau.2020.06.009

    MENG Y Y, LIU B, KANG J J. Effects of different irrigation and fertilizer application on soil respiration and nitrification-denitrification in sweet Sorghum farmland [J]. Journal of Gansu Agricultural University, 2020, 55(6): 67-76(in Chinese). doi: 10.13432/j.cnki.jgsau.2020.06.009

    [4] 赵永鹏. 耕作和施肥对水稻土厌氧氨氧化和反硝化速率的影响[D]. 重庆: 西南大学, 2017.

    ZHAO Y P. Effects of tillage and fertilization on anammox and denitrification rates in paddy soils[D]. Chongqing: Southwest University, 2017(in Chinese).

    [5] 陈娜, 刘毅, 黎娟, 等. 长期施肥对稻田不同土层反硝化细菌丰度的影响 [J]. 中国环境科学, 2019, 39(5): 2154-2160. doi: 10.3969/j.issn.1000-6923.2019.05.044

    CHEN N, LIU Y, LI J, et al. Effects of long-term fertilization on the abundance of the key denitrifiers in profile of paddy soil profiles [J]. China Environmental Science, 2019, 39(5): 2154-2160(in Chinese). doi: 10.3969/j.issn.1000-6923.2019.05.044

    [6] 刘发波, 马笑, 张芬, 等. 硝化抑制剂对我国蔬菜生产产量、氮肥利用率和氧化亚氮减排效应的影响: Meta分析 [J]. 环境科学, 2022, 43(2): 5140-5148.

    LIU F B, MA X, ZHANG F, et al. Impact of nitrification inhibitors on vegetable production yield nitrogen fertilizer use efficiency and nitrous oxide emission reduction in China: Meta Analysis [J]. Environmental Science, 2022, 43(2): 5140-5148(in Chinese).

    [7] 王利民, 林新坚, 黄东风, 等. 红黄壤茶园不同培肥模式的土壤理化效应 [J]. 东北林业大学学报, 2012, 40(1): 54-57,65. doi: 10.3969/j.issn.1000-5382.2012.01.013

    WANG L M, LIN X J, HUANG D F, et al. Effect of different fertilization patterns on physicochemical properties of red-yellow soil in tea garden [J]. Journal of Northeast Forestry University, 2012, 40(1): 54-57,65(in Chinese). doi: 10.3969/j.issn.1000-5382.2012.01.013

    [8] 李力, 刘娅, 陆宇超, 等. 生物炭的环境效应及其应用的研究进展 [J]. 环境化学, 2011, 30(8): 1411-1421.

    LI L, LIU Y, LU Y C, et al. Review on environmental effects and applications of biochar [J]. Environmental Chemistry, 2011, 30(8): 1411-1421(in Chinese).

    [9] 简秀梅, 陈学濡, 刘富豪, 等. 不同灰分生物质炭对红壤理化特性与微生物特性的影响 [J]. 农业机械学报, 2020, 51(6): 282-291. doi: 10.6041/j.issn.1000-1298.2020.06.030

    JIAN X M, CHEN X R, LIU F H, et al. Effects of biochars with various ash contents on physiochemical properties and microbial community quantity and activity in red soils [J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(6): 282-291(in Chinese). doi: 10.6041/j.issn.1000-1298.2020.06.030

    [10] 杨润. 沼液对水稻产量、氮素利用、氨挥发的影响及氨挥发减缓措施探讨[D]. 南京: 南京农业大学, 2017.

    YANG R. Effects of biogas slurry on rice yield, nitrogen use efficiency and ammonia volatilization, and the discussion on measures for ammonia volatilization reduction[D]. Nanjing: Nanjing Agricultural University, 2017(in Chinese).

    [11] REN T T, YU X Y, LIAO J H, et al. Application of biogas slurry rather than biochar increases soil microbial functional gene signal intensity and diversity in a poplar plantation [J]. Soil Biology and Biochemistry, 2020, 146: 107825. doi: 10.1016/j.soilbio.2020.107825
    [12] 潘逸凡. 生物质炭对稻田土壤氨氧化微生物的影响研究[D]. 杭州: 浙江大学, 2014.

    PAN Y F. Influence of biochar on ammonia oxidation microbial community in paddy soil[D]. Hangzhou: Zhejiang University, 2014(in Chinese).

    [13] DEMPSTER D N, GLEESON D B, SOLAIMAN Z M, et al. Decreased soil microbial biomass and nitrogen mineralisation with Eucalyptus biochar addition to a coarse textured soil [J]. Plant and Soil, 2012, 354(1/2): 311-324.
    [14] 刘杰云, 邱虎森, 王聪, 等. 生物质炭对双季稻田土壤反硝化功能微生物的影响 [J]. 环境科学, 2019, 40(5): 2394-2403. doi: 10.13227/j.hjkx.201810182

    LIU J Y, QIU H S, WANG C, et al. Influence of biochar amendment on soil denitrifying microorganisms in double rice cropping system [J]. Environmental Science, 2019, 40(5): 2394-2403(in Chinese). doi: 10.13227/j.hjkx.201810182

    [15] 王伟, 白军红, 张玲, 等. 基于Cite Space的生物质炭对土壤氮循环影响的文献计量分析 [J]. 北京师范大学学报(自然科学版), 2021, 57(1): 76-85.

    WANG W, BAI J H, ZHANG L, et al. Biochar modulation of the soil nitrogen cycle: A Cite Space-based bibliometric study [J]. Journal of Beijing Normal University (Natural Science), 2021, 57(1): 76-85(in Chinese).

    [16] 潘逸凡, 杨敏, 董达, 等. 生物质炭对土壤氮素循环的影响及其机理研究进展 [J]. 应用生态学报, 2013, 24(9): 2666-2673. doi: 10.13287/j.1001-9332.2013.0518

    PAN Y F, YANG M, DONG D, et al. Effects of biochar on soil nitrogen cycle and related mechanisms: a review [J]. Chinese Journal of Applied Ecology, 2013, 24(9): 2666-2673(in Chinese). doi: 10.13287/j.1001-9332.2013.0518

    [17] 鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000.

    LU R K. Methods for soil agrochemical analysis [M]. Beijing: China Agricultural Science and Technology Press, 2000(in Chinese).

    [18] YANG X L, LU Y L, DING Y, et al. Optimising nitrogen fertilisation: A key to improving nitrogen-use efficiency and minimising nitrate leaching losses in an intensive wheat/maize rotation (2008-2014) [J]. Field Crops Research, 2017, 206: 1-10. doi: 10.1016/j.fcr.2017.02.016
    [19] 储成, 吴赵越, 黄欠如, 等. 有机质提升对酸性红壤氮循环功能基因及功能微生物的影响 [J]. 环境科学, 2020, 41(5): 2468-2475. doi: 10.13227/j.hjkx.201911013

    CHU C, WU Z Y, HUANG Q R, et al. Effect of organic matter promotion on nitrogen-cycling genes and functional microorganisms in acidic red soils [J]. Environmental Science, 2020, 41(5): 2468-2475(in Chinese). doi: 10.13227/j.hjkx.201911013

    [20] 董心亮, 林启美. 生物质炭对土壤物理性质影响的研究进展 [J]. 中国生态农业学报, 2018, 26(12): 1846-1854. doi: 10.13930/j.cnki.cjea.180277

    DONG X L, LIN Q M. Biochar effect on soil physical properties: A review [J]. Chinese Journal of Eco-Agriculture, 2018, 26(12): 1846-1854(in Chinese). doi: 10.13930/j.cnki.cjea.180277

    [21] ESMAEELNEJAD L, SHORAFA M, GORJI M, et al. Impacts of woody biochar particle size on porosity and hydraulic conductivity of biochar-soil mixtures: An incubation study [J]. Communications in Soil Science and Plant Analysis, 2017, 48(14): 1710-1718. doi: 10.1080/00103624.2017.1383414
    [22] 万辰, 马瑛骏, 张克强, 等. 洱海流域不同有机肥替代对土壤理化性质及油菜产量的影响 [J]. 农业环境科学学报, 2021, 40(11): 2494-2502. doi: 10.11654/jaes.2021-1026

    WAN C, MA Y J, ZHANG K Q, et al. Effects of different organic fertilizer substitutions on soil physical and chemical properties and rapeseed yield in erhai lake basin, China [J]. Journal of Agro-Environment Science, 2021, 40(11): 2494-2502(in Chinese). doi: 10.11654/jaes.2021-1026

    [23] CLOUGH T J, BERTRAM J E, RAY J L, et al. Unweathered wood biochar impact on nitrous oxide emissions from a bovine-urine-amended pasture soil [J]. Soil Science Society of America Journal, 2010, 74(3): 852-860. doi: 10.2136/sssaj2009.0185
    [24] JIN Z J, LI L Q, LIU X Y, et al. Impact of long-term fertilization on community structure of ammonia oxidizing and denitrifying bacteria based on amoa and nirk genes in a rice paddy from tai lake region, China [J]. Journal of Integrative Agriculture, 2014, 13(10): 2286-2298. doi: 10.1016/S2095-3119(14)60784-X
    [25] 周慧芳. 农田土壤厌氧氨氧化及其在氮循环中影响的研究[D]. 杭州: 浙江大学, 2018.

    ZHOU H F. Study on anammox of farmland soil and its effects on nitrogen cycling[D]. Hangzhou: Zhejiang University, 2018(In Chinese).

    [26] MURASE J, HIDA A, OGAWA K, et al. Impact of long-term fertilizer treatment on the microeukaryotic community structure of a rice field soil [J]. Soil Biology and Biochemistry, 2015, 80: 237-243. doi: 10.1016/j.soilbio.2014.10.015
    [27] 赵伟, 王宏燕, 王大庆, 等. 农肥和化肥对东北黑土土壤氮素转化作用的研究 [J]. 水土保持学报, 2009, 23(2): 99-102,108. doi: 10.3321/j.issn:1009-2242.2009.02.022

    ZHAO W, WANG H Y, WANG D Q, et al. Effects of manure and chemical fertilizers application on nitrogen transfer in black soil [J]. Journal of Soil and Water Conservation, 2009, 23(2): 99-102,108(in Chinese). doi: 10.3321/j.issn:1009-2242.2009.02.022

    [28] 马龙, 高伟, 栾好安, 等. 有机肥/秸秆替代化肥模式对设施菜田土壤氮循环功能基因丰度的影响 [J]. 植物营养与肥料学报, 2021, 27(10): 1767-1778. doi: 10.11674/zwyf.2021136

    MA L, GAO W, LUAN H A, et al. Effects of partial substitution of chemical fertilizer with manure and/or straw on the abundance of functional genes related to soil N-cycling [J]. Journal of Plant Nutrition and Fertilizers, 2021, 27(10): 1767-1778(in Chinese). doi: 10.11674/zwyf.2021136

    [29] 郭俊杰, 朱晨, 刘文波, 等. 不同施肥模式对土壤氮循环功能微生物的影响 [J]. 植物营养与肥料学报, 2021, 27(5): 751-759. doi: 10.11674/zwyf.20504

    GUO J J, ZHU C, LIU W B, et al. Effects of different fertilization managements on functional microorganisms involved in nitrogen cycle [J]. Journal of Plant Nutrition and Fertilizers, 2021, 27(5): 751-759(in Chinese). doi: 10.11674/zwyf.20504

  • 加载中
图( 4) 表( 3)
计量
  • 文章访问数:  2027
  • HTML全文浏览数:  2027
  • PDF下载数:  66
  • 施引文献:  0
出版历程
  • 收稿日期:  2022-03-06
  • 录用日期:  2022-06-05
  • 刊出日期:  2023-08-27

生物质炭-沼液联合施用对调控氮循环功能基因促进氮素增效的影响

    通讯作者: E-mail:119044@zust.edu.cn
  • 1. 浙江科技学院环境与资源学院,杭州,310023
  • 2. 浙江省废弃生物质循环利用与生态处理重点实验室,杭州,310023
  • 3. 浙江省耕地质量与肥料管理总站,杭州,310020
  • 4. 开化县新农村建设中心,开化,324300
基金项目:
国家重点研发计划项目“数据和知识联合驱动城市易腐垃圾炭化与绿色可持续利用的关键技术及应用”(2022YFE0196000)和浙江省重点研发项目(2019C02061,2020C01017)资助.

摘要: 为探讨生物质炭-沼液配施条件下氮循环功能基因调控农田土壤氮素转化并影响农作物氮素吸收利用机制. 本试验以浙江省杭州市红黄壤作为研究对象,设置生物质炭和沼液两个因素,探究生物质炭-沼液配施条件下土壤基本理化性质和氮循环功能基因丰度变化情况,刻画功能基因与农田氮素利用率间的耦合关系. 结果表明,生物质炭-沼液配施可以显著降低土壤容重,提升土壤pH和土壤氮素含量,其中,高剂量生物质炭-沼液配施(C3B2)处理较单施化肥(C0B0)处理铵态氮、硝态氮、全氮含量增幅均达到显著水平(P<0.05). 与空白处理(CK)相比,生物质炭-沼液配施(C3B2)处理则显著提高了反硝化功能基因丰度,较单施化肥(C0B0)处理增幅30.98%和44.99%. 冗余分析结果显示,铵态氮、硝态氮和有机碳含量对土壤氮循环功能基因影响较为显著,结构方程模型则表明硝化作用功能基因丰度的提升对包菜氮素农学利用率呈现负相关趋势. 研究结果表明,在相同养分施用量的条件下,生物质炭-沼液配施可显著提高土壤肥力. 氮素和有机碳含量是影响功能基因丰度的关键因素,硝化作用功能基因丰度的降低可以提高农田氮素利用率. 本研究结果可以为促进农业废弃物资源化利用,提升农田土壤氮素活性提供理论依据.

English Abstract

  • 氮素是农作物生长发育过程中必不可少的营养元素,近年来由于过量施用化肥一方面加剧了硝酸盐淋失,另一方面氮素通过反硝化作用生成N2O、N2释放到大气中,导致氮素利用率下降[1]. 氮素形态转换的本质是氮素在相关微生物功能基因调控下的转化过程,刻画氮循环功能基因与农田氮素利用率间耦合关系,明确主导氮素利用率提升的关键基因种类有利于为进一步开展土壤健康管理和促进农业增产提供科学依据. 现有研究表明,过量施加氮肥会导致土壤中积累的氨通过硝化作用在氨氧化细菌(AOB)的作用下将环境中的氨转化为亚硝酸,而亚硝酸盐氧化细菌(NOB)通过亚硝酸盐氧化还原酶(NXR)将亚硝酸盐氧化为硝酸盐,加剧了土壤酸化[2]. 同时,施肥可以改变土壤nosZ型反硝化细菌基因丰度,并且提高土壤反硝化速率[3-4]. 陈娜等[5]研究发现,增施氮肥可以提升(narG、nosZ)型反硝化功能基因丰度,从而使NO3、NO2被还原为N2O、N2释放到大气中,导致农田氮素损失. 研究表明抑制硝化作用功能基因丰度不仅可以缓解氮素以硝酸盐形式流失,还可以促进植株对氮素吸收,从而提升氮素利用率[6]. 综上,研究土壤中氮素的生物化学转化,明确施入土壤的氮素去向进而刻画氮循环功能基因丰度与氮素利用率间的耦合关系可以为减缓农田氮素流失,促进氮素增效提供理论基础.

    红黄壤主要分布在我国南方地区,土壤pH为5.25—5.44、降雨量充沛、土壤质地较黏、容重较高. 由于红黄壤本身的性质和当地气候的作用使施入土壤中的氮素通过挥发、淋失、反硝化等进入大气环境以及湖泊河流中[7]. 因此,减少氮素流失,提高农田氮素利用率已成为提高农作物增产增收的的关键所在. 生物质炭和沼液作为典型的土壤改良剂,通过增加土壤有机碳含量为氮循环微生物提供碳源,直接或间接调控氮循环过程相关功能基因丰度和关键环节,从而提高土壤氮素利用效率以及作物产量. 生物质炭是农作物秸秆或动物粪便残渣在厌氧的条件下经过高温热解所产生,具有巨大的比表面积,可以吸附土壤中的营养物质,而且作为肥料的载体,它们还影响着土壤持水性能,可以提高农作物对养分的吸收、改良土壤结构、促使土壤微生物繁殖[8-9]. 沼液作为改善土壤质量的有机物料对氮循环功能基因同样具有显著的调控作用. 由于单独施用沼液易产生氨挥发、硝酸盐淋失和N2O排放等主要环境问题,因此采用生物质炭-沼液联合施用的方式,利用生物质炭吸附和负载氮素,在降低硝酸盐淋失的同时减少N2O排放,提高土壤养分保蓄能力[10-11]. 现有研究发现生物质炭施加到土壤中可以影响硝化作用的限速步骤,促进氨氧化作用的进行,提高氨氧化菌的基因丰度,从而达到促进硝化作用的发生[12]. 而Dempster等[13]的研究则得出生物质炭的施加会对土壤氨氧化作用产生抑制的影响,这可能是由于生物质炭能释放硝化抑制剂,阻碍土壤中铵根离子向硝酸根离子的转化. 而生物质炭对土壤反硝化作用的影响至今仍有不同的观点,刘杰云等[14]研究表明,生物质炭作为土壤改良剂可以显著提高narG、nirK、nosZ基因相对丰度. 王伟等[15]研究发现,通过施加生物质炭可以改善土壤孔隙度,增加了土壤的透气状况,阻碍了厌氧条件下反硝化作用的进行. 施加沼液可以增加土壤微生物的数量,并且沼液的施加也放大了反硝化(nirK、narG)功能基因信号强度、氨化作用(ureC)功能基因信号强度[11].

    综上,前人在探讨氮循环功能基因对氮素形态调控以及影响土壤氮素循环过程做了大量研究,为明确氮循环功能基因对农田氮素利用率的影响奠定了基础. 然而,现有研究主要集中在研究功能基因丰度对氮素活化的影响,氮素增效与功能基因间的内在联系仍有待探索[16],其次,目前的研究多采用生物质炭/无机肥料或沼液/无机肥料配施的方式改良土壤,而生物质炭-沼液协同施用的研究手段仍较为少见. 因此,本研究聚焦农田氮素增效的关键问题,从氮循环功能基因角度入手,以浙江省杭州市红黄壤为试验对象,采用生物质炭-沼液联合施用的方式研究了其对土壤肥力因子、氮循环功能基因和氮素利用率的响应特征,刻画氮循环功能基因与氮素利用率间耦合关系,为改善土壤性质,提高氮素利用率以及为缓解因过量施肥而造成的农业污染问题提供理论依据.

    • 该试验地点位于浙江省杭州市浙江科技学院(30°13′N,120°1′E). 该试验土壤为红黄壤,属于亚热带季风气候,年平均气温17.8 ℃、年平均日照数1970 h、年平均降雨量为1573 mm. 基本理化性质为:土壤容重1.45 g·cm-3,pH为5.44,硝态氮含量3.54 mg·kg−1,铵态氮含量0.79 mg·kg −1,全氮含量0.6 g·kg−1.

    • 供试生物质炭原材料来自浙江省开化县顺康牧业猪粪,由浙江金华市金锅公司厌氧炭化制成猪粪生物质炭,理化性质为:硝态氮含量11.37 mg·kg−1、铵态氮含量15 mg·kg−1、全氮含量18.95 g·kg−1、有机碳含量605 g·kg−1、pH 8.5. 沼液为浙江省开化县顺康牧业猪场废水厌氧发酵制备的沼液,试验站建有专门蓄水调节池,可对沼液进行标准化处理,经标准化处理的沼液符合农田灌溉水质标准(GB5084-2005). 沼液理化性质为:硝态氮含量10 mg·L−1、铵态氮含量220 mg·L−1、全氮含量2.95 g·L−1、pH 8.5.

    • 2021年4月开始试验,采用生物质炭和沼液完全交互试验设计:生物质炭设置4个水平,分别为0%、0.5%、1%、2%,沼液设置3个水平分别为0%、50%、100%(详见表1),其中0%沼液处理采用等体积超纯水代替. 按照当地栽培方式最佳施肥量为:尿素(折合纯氮90 kg·hm−2)、氯化钾(195 kg·hm−2)、钙镁磷肥(375 kg·hm−2). 每个处理3个重复,同时设置对照(CK)处理,共计39个小区. 各小区面积3 m2(2 m×1.5 m). 试验基地种植农作物为包菜,种植采取育苗方式,再移栽至每个试验小区. 行距为40 cm,株距为20 cm. 于试验开始前两周将猪粪生物质炭一次性均匀的施于土壤表面,人工翻至15—30 cm土层,后期不再施加生物质炭. 对照组(CK)处理不施生物质炭和沼液,但仍需人工翻耕. 农作物的种植方式与当地保持一致,沼液和化肥于农作物种植一周前均匀施入其对应试验小区,并且后期每季作物种植前均需再次施用沼液和化肥,各小区每隔一个月施1次沼液,共施加3次.

    • 供试土壤于包菜成熟前两周采集,用S形采样方法采样,每个小区采集5个样本并混合均匀冷藏于带有冰块的泡沫箱内运回实验室. 2 d内除去样品中沙粒、枯叶等,把样品均匀分成3份. 一份风干,用于测土壤基本理化性质,参照《土壤农业化学分析方法》鲁如坤[17]. 一份冷藏于4 ℃冰箱测定土壤硝态氮、铵态氮,另一份储存于-80 ℃冰箱中用于提取土壤DNA. 土壤容重采用环刀法测定,土壤pH采用pH计电位法测定(水土比5:1),土壤硝态氮、铵态氮采用氯化钾浸提-流动分析仪测定,土壤全氮的测定采取浓硫酸消煮-流动分析仪测定. 每季包菜收获期各小区按照包菜颗数进行收获计产,再换算得至每公顷包菜产量. 所有植株样品经105 ℃杀青烘干,烘干样品粉碎过100目筛,经H2SO4-H2O2消解后测定全氮含量,计算氮素利用率. 氮素利用率测算公式参照Yang所述方法[18].

      包菜氮素农学利用率(ANUE,kg·kg -1) = (施氮区包菜产量-不施氮区包菜产量)/施氮量

    • 取所需土壤鲜样0.5 g,每个处理取3个平行. 使用中国Tian-Gen公司的土壤基因组DNA快速提取试剂盒TIANNAMP Soil DNA Kit提取土壤微生物总DNA,提取完成后分别利用TBS-380、NanoDrop200对土壤样品DNA进行浓度及纯度检测.

    • 采用BIO-RAD T100TM扩增仪进行,扩增反应体系为:5 μL 10×Ex Taq buffer (Mg2+ Plus),4 μL dNTP(2.5 mmol·L−1),0.5 μL上游和下游引物(20 μmol·L−1),0.25 μL Ex Taq DNA聚合酶(5 U·L−1),1 μL DNA模板,无菌去离子水补足至50 μL. 采用质粒试剂盒(Plasmid Miniprep Kit,Biomiga) 进行质粒提取,质粒DNA提取液经2%琼脂糖凝胶电泳检测、纯化后测定质粒浓度,由质粒浓度换算出质粒拷贝数. 标准曲线用无菌去离子水稀释10倍后的质粒DNA作为模板,在Bio-Rad miniopticon扩增仪上进行扩增,功能基因PCR扩增引物序列见表2[19].

    • 采用Microsoft Excel 2016、Origin 2021及SPSS 20.0软件对试验数据进行方差分析、均值的多重比较和相关性分析,显著性水平设为P=0.05,采用CANOCO 5软件 (CANOCO, Microcomputer PowerInc, Ithaca, NY, USA) 进行土壤环境因子与氮循环功能基因间的冗余分析 (redundancy analysis, RDA),采用Amos刻画功能基因与氮素利用率间的结构方程模型.

    • 本研究结果表明施加生物质炭、沼液以及生物质炭-沼液配施均可降低土壤容重,提升土壤氮素含量(表3). 单施化肥(C0B0)处理的土壤容重、pH、硝态氮、全氮、总有机碳含量与空白处理(CK)相比差异性不显著(P>0.05). 单独施用生物质炭可以显著提升土壤pH,本研究中C3B0处理土壤pH值较C0B0提升了45.43%,同时,生物质炭-沼液配施(C3B2)处理土壤pH提升幅度较大,比单施化肥处理(C0B0)提升了44.85%,其原因是由于生物质炭本身所具有的含氧官能团在土壤中发生了质子化作用导致土壤pH增加[20]. 施加生物质炭、沼液可降低土壤容重,其中C1B0、C1B1、C1B2处理土壤容重均比单施化肥(C0B0)处理显著降低了13.79%、11.72%、11.03%. 生物质炭施加量为0.5%、沼液0%替代化肥(C1B0)处理中土壤容重较单施化肥(C0B0)处理降幅最为显著(P<0.05),下降了13.79%. 由于生物质炭具有巨大的比表面积和孔隙度,施入土壤后其稀释效应增加了微生物与土壤颗粒物间的相互作用, 促进了土壤微团粒向大团粒的聚合,降低了土壤容重[21]. 土壤有机碳和氮素含量以生物质炭-沼液配施(C3B2)处理效果最佳,均较单施化肥(C0B0)处理有显著提升,并且均高于单独施生物质炭(C3B0)处理,其中C3B2处理土壤硝态氮含量较C3B0处理增加了86%,铵态氮较C3B0处理增加了23.79%. 因为生物质炭对于施入土壤中的养分具有吸附作用,沼液则富含有机质并且其主要氮素含量以无机氮为主,因此生物质炭和沼液联合施用可以提高土壤硝态氮、有机碳含量. 沼液中富含腐殖酸,其中的官能团与土壤中含氮物质形成络合物减少了土壤铵态氮的流失,并且沼液作为有机碳库,施入土壤中提升了土壤微生物的活性,从而增加了土壤速效养分含量[22].

    • 试验发现单施化肥处理可以显著提高硝化作用功能基因丰度,而施加生物质炭则能显著提高反硝化作用功能基因丰度,土壤氮循环相关基因丰度变化如图1所示.

      图1A中单施化肥(C0B0)处理较对照组(CK)硝化作用基因(amoA-1、amoA-2、amoB和nrxA)丰度提升显著,差异达显著水平(P<0.05). 而单独施加生物质炭(C3B0)处理较单施加化肥(C0B0)处理硝化作用功能基因(amoA-1、amoA-2、amoB和nrxA)显著降低,分别下降了93%、35%和16%. 说明施加生物质炭对土壤硝化作用起到了抑制作用,进而可缓解土壤中硝酸盐淋失的问题. 这与Clough的研究得到了一致的结果,其认为生物质炭可释放一种硝化抑制剂α-松萜从而抑制了土壤硝化作用[23]. 而在本试验中单施化肥处理(C0B0)与空白处理(CK)相比硝化作用功能基因丰度有显著提升,但随着生物质炭与沼液的施加导致其功能基因丰度增长率逐渐降低. 由于亚硝酸细菌的适宜生长环境为偏酸性,而生物质炭与沼液的施加提升了土壤pH,因此pH的提高抑制了硝化作用功能基因丰度的提升. 单施生物质炭(C3B0)处理较单施化肥(C0B0)处理反硝化作用基因(narG、napA、nirS-1、nirS-2和nirK-1)丰度则提升了33%、182%、189%、131%和102%(图1B). 同时,生物质炭-沼液配施(C3B2)处理与单施化肥(C0B0)处理相比,反硝化作用基因丰度提升较为显著(P<0.05). 可能是由于生物质炭、沼液的投入增加了土壤中的氮源以及土壤有机碳的含量,补充了土壤反硝化细菌所需的的底物. 也有研究发现,有机肥部分替代化肥增加了反硝化功能微生物(narG、nirK、nirS和norZ)数量与丰度,甚至形成一些较为特殊的反硝化微生物群落结构[24]. 单施化肥(C0B0)处理、单施沼液(C0B2)处理均较(CK)处理显著提升了土壤厌氧氨氧化功能基因丰度,单施生物质炭(C3B0)处理、生物质炭-沼液配施显(C3B2)处理与单施沼液(C0B0)处理相比,著降低了土壤氨氧化作用功能基因丰度(P<0.05)(图1D). 研究表明,厌氧氨氧化菌的数量与土壤尿素呈正相关,但与土壤铵态氮、硝态氮含量呈负相关趋势. 因此施肥改变了土壤肥力,提升了土壤氨氧化细菌的数量,而随着生物质炭、沼液的施加增加了土壤硝态氮与铵态氮的含量,进而阻碍了土壤厌氧氨氧化作用[25-26]. 单施加生物质炭(C3B0)处理、生物质炭-沼液配施(C3B2)处理与单施化肥(C0B0)处理相比氨化作用功能基因(ureC)均有显著增加,分别提升了36.5%、62.5%,差异均达到显著水平(P<0.05)(图1C),由于有机物料施入到土壤中后,需经过土壤中的微生物将有机氮肥转化为无机氮才能被植株所吸收,所以有机物料的施加提升了土壤氨化作用功能基因丰度[27].

    • 不同处理包菜产量与氮素利用率如图2所示,单施化肥(C0B0)处理的包菜产量较对照组(CK)提高了33.31 t·hm−2,差异达显著水平(P<0.05). 单施用生物质炭可以显著提高包菜产量,其中(C3B0)处理包菜产量为59.93 t·hm−2,较单施化肥(C0B0)处理提高了41.08%. 生物质炭-沼液配施(C3B2)处理包菜产量最高,与单施化肥(C0B0)处理相比差异显著(P<0.05),增幅为44.13%. 单独施用生物质炭可以显著提升包菜氮素利用率,其中C3B0处理包菜氮素利用率最高,较单施肥(C0B0)处理提升42.80%,差异达显著水平(P<0.05). 生物质炭-沼液配施(C3B2)处理较单施化肥(C0B0)处理包菜氮素利用率有显著提升(P<0.05)且升高了44.99%. 由于生物质炭、有机物料的投入与替代化肥提升了土壤肥力,可为作物的生长发育提供适宜的生存环境与充足的氮源,因此可对包菜的氮素利用率起到提升作用.

    • 利用氮循环功能基因与所选取的6个环境因子(pH、SBD、NH4+-N、NO3-N、TN、SOC)做冗余分析. 氮循环功能基因分别选取amoA、amoB、napA、ureC、nirK、nosZ、nirS,由图3可知,第一排序轴和第二排序轴分别解释了94.63%、4.6%. 影响土壤硝化作用功能基因丰度amoA的环境因子依次是NO3-N>NH4+-N>TN,amoB功能基因丰度与NO3-N、NH4+-N、TN呈正相关,amoA功能基因丰度与pH、NO3-N、NH4+-N、SOC呈负相关,与土壤TN含量呈正相关. 对反硝化作用功能基因nosZ、nirK丰度影响作用较为显著的环境因子是NO3-N、NH4+-N、SOC、pH、nirK、nirS功能基因丰度分别与pH、NO3-N、NH4+-N、SOC呈显著正相关关系. 在六种环境因子(pH、SBD、NH4+-N、NO3-N、TN、SOC)当中,NH4+-N(P=0.02)、NO3-N(P=0.002)、SOC(P=0.04)对土壤氮循环功能基因影响较为显著. 其分别解释了土壤中氮循环功能基因丰度变异的21.5%、66.5%、13.7%.

    • 氮循环作为土壤养分循环的重要组成部分之一,而作为硝化、反硝化又是氮循环的重要组成部分,其组成细菌的丰度和群落结构直接影响施入土壤中氮素的转化,并影响农作物对氮素的吸收利用[28]. 以氮循环功能基因和作物产量为观测变量,农田氮素利用率为潜变量构建一阶SEM(图4),模型的拟合度指数良好(CMIN/DF = 3.286<5,RMSAE = 0.078<0.08,AGFI、GFI、CFI、NFI>0.8). 由图4可以看出,土壤反硝化作用可以由narG、napA、nirS和nosZ这4个观测变量来测量,其中narG、napA、nirS,对反硝化作用影响较为显著且为正相关效应,相关系数分别为0.89、0.89、0.90,土壤硝化作用可以由amoA、amoB、nrxA和hao的4个观测变量来测量,amoA、nrxA对硝化作用的影响显著,路径系数分别为0.98、0.77. 结构方程模型结果表明,反硝化功能基因丰度的提升有利于包菜农学氮素利用率的提升(路径系数0.94),而硝化作用功能基因(amoB、hao)丰度的提升对包菜氮素农学利用率呈现负相关趋势(路径系数-0.21). 产生这一趋势的原因可能是由于施入土壤的氮肥在氨氧化细菌的参与下将氮素转化为硝态氮,导致了土壤中氮素淋失,也可能是因为铵态氮被土壤胶体吸附不易被作物吸收[29]. 同时,氨化作用、同化氮还原和厌氧氨氧化功能基因也是农田氮素利用率的正向影响因素(路径系数分别为0.10、0.17和0.22). 另外,农田氮素利用率的提升显著促进了包菜的增产(路径系数0.98).

    • 生物质炭的输入可以缓解土壤酸化问题、降低土壤容重,施加沼液可以增加土壤氮素含量提高土壤肥力. 单独施加化肥增加了土壤硝化作用(amoA、amoB)功能基因丰度,但随着施生物炭和沼液的施加会使其增长率逐渐降低. 生物质炭与沼液配施提升了土壤反硝化作用(narG、nirS、nirK)功能基因丰度. 因此,可以根据包菜在不同生理期对硝态氮和铵态氮吸收比例的不同适当减少化肥的供应量以减少因硝化作用而产生的硝酸盐淋失. 生物质炭与沼液配施提升了土壤肥力,丰富了土壤氮循环微生物群落结构,促使包菜增产并提高了包菜氮素农学利用率. 因此结合当地种植背景,生物质炭与沼液部分替代化肥模式为绿色、可持续发展的施肥管理模式.

    参考文献 (29)

目录

/

返回文章
返回