以小麦“淮麦33”为试验材料，水稻秸秆与土壤均取自涟水高沟镇试验基地，秸秆不做任何操作(收割机已打碎处理)，土壤碾碎至颗粒不超过1.0 cm³。沼液取自基地周边生猪养殖场，土壤及沼液理化性质下表所示(表1)。

表1. 耕层土壤及沼液的基本理化性质

2022年10月30日进行小麦播种处理。计算单位面积水稻秸秆重量，并进行秸秆含氮量计算(水稻秸秆单位重量测算为7625 kg·hm⁻²，全氮含量为5.9 g·kg⁻¹)，最后计算秸秆含氮量占常规化肥施用量中全氮的占比(以小麦产量550 kg测算，常规施化肥为15%~15%~15% [ 1 ] 的NPK酸性复合肥，施氮量为276.45 kg·hm⁻²)，剩余施氮量由沼液或秸秆补充(水稻秸秆含氮量即为7625*5.9 = 45.0 kg·hm⁻²，其含氮量在常规施肥总施氮量276.45 kg·hm⁻²中占比为16.28%)。试验设置5个处理组：1) 对照组(不秸秆还田，不施沼液，不施化肥，记为CK)；2) 化肥处理组(100%化肥N，记为CF)；3) 秸秆全量还田 + 化肥处理(16.28%秸秆N，83.72%化肥N，记为S-CF)；4) 秸秆全量还田 + 沼液处理(16.28%秸秆N，83.72%沼液N，记为S-BS)；5) 沼液组(100%沼液N，记为BS)。其中，肥料运筹比例为基肥：追肥 = 1:1进行管理施用，每处理重复3次。试验小区面积为4 m × 7 m。小区种植管理常规进行。

1) 株高、茎基宽在小麦生长的拔节期、孕穗期和灌浆期随机选取具有代表性的植株，分别测量株高、茎基宽指标。

2) 小麦产量待小麦成熟后，各处理组随机选取1 m²进行产量测定，3次重复并取平均值，然后进行理论产量分析。

3) 土壤总DNA提取使用MoBio PowerSoil^®DNA isolation kit试剂盒提取各组处理土壤样本总DNA，用1%琼脂糖凝胶电泳检测和分光光度法(260 nm/280 nm光密度比)进行质量检测，提取得到的DNA样品检测后，于−80℃保存以备后续试验使用。

4) 特异区域扩增及测序DNA样本送至南京集思慧远生物科技有限公司，微生物多样性检测选取真菌rDNA ITS2区，扩增引物为ITS3F (5' - GCATCGATGAAGAACGCAGC-3')和ITS4R (5' - TCCTCCGCTTATTGATATGC-3') [ 11 ] 。PCR反应总体系25 µl：10 µl 2 × EasyTaq PCR SuperMix (Transgen Bio, Inc)、1 µl正向引物、1 µl反向引物、加入的DNA总量为10 ng，最后加dd H₂O补至25 µl。反应参数：94℃预变性5 min；94℃变性30 s，55℃退火30 s，72℃延伸1 min，35个循环；72℃延伸8 min。利用Illumina Novaseq PE250高通量测序平台测序。

微生物群落结构分析下机数据经过QIIME软件过滤、拼接、去除嵌处理，另外去除碱基质量分值低于20、测序长度小于150 bp、碱基模糊和引物错配序列。根据barcodes归类处理组序列信息聚类为物种分类的基本操作单元(Operational taxonomic units, OTU)，设置OTU相似性为97%，与silva真菌类数据库进行比对 [ 12 ] 。通过QIIME软件进行阿尔法(Alpha)多样性分析 [ 13 ] 。

在小麦的拔节期、孕穗期和灌浆期生长过程中，CF、S-CF、S-BS和BS处理均显著高于CK。株高方面，CF、S-CF和S-BS在不同生长时期均无显著差异；灌浆期，单独的沼液BS处理显著低于CF、S-CF和S-BS，不排除沼液后期养分欠缺的情况。茎基宽方面，在拔节期、孕穗期除CK外，其余处理组均无显著差异，但是在灌浆期，CF、S-CF与S-BS、BS存在显著差异，说明含有沼液的试验组合在小麦生长后期可能会养分供应不足，影响茎的发育(表2)。

表2. 不同处理对小麦株高、茎基宽的影响

注：表中同列数据后不同小写字母的表示差异显著(P < 0.05)，下同。

小麦理论产量方面，CF、S-CF、S-BS、BS均显著高于CK组。S-BS组的理论产量与CF、S-CF无明显差异，说明沼液和秸秆联合施用一定程度能够代替化肥的施用，这与S-BS组的株高在不同生长时期与CF、S-CF无明显差异情况一致。另外，沼液组合可以增加单位面积穗数，但是存在无效分蘖现象，故BS的千粒重较低，与化肥组合均存在显著差异。由于沼液可能导致小麦营养生长旺盛，生长过程的始穗期、齐穗期、成熟期推迟，存在贪青迟熟现象，因此BS与CK相比，熟期延迟高达7~8天(表3)。

表3. 不同处理对小麦理论产量的影响

不同试验处理真菌群落的Alpha多样性goods_coverage指标均达到了0.99，说明本次检测具有较高的微生物覆盖率。Chao1和observed_species主要是估计群落中包含物种的数目指标，具有一定的正向关联性。空白处理CK组中物种的数目最多，单独施肥的BS、CF显著减少物种的数目，说明沼液、化肥、秸秆的联合施用会降低群落中物种的数量，尤其单独施用某种肥料。Shannon指数是指信息熵，其值指数越大，表示不确定性大，说明这个群落中未知因素越多，物种多样性越高。S-CF、S-BS的联合施肥处理的shannon值与CK、CF存在显著差异。因此，联合施肥的S-CF、S-BS能够提升物种的丰富度和多样性，并与理论产量密切关联(表4)。

表4. 土壤真菌群落Alpha多样性分析

作物秸秆腐熟分解后大量的营养物质是作物果蔬生长所需肥料的重要来源。沼液可以作为一种速效液体生物肥料，与秸秆联合施用能够促进秸秆的转化利用。秸秆还田腐解处理后，会进一步增加土壤的养分含量，提高腐殖质、有机酸等，促进作物生长发育。在沼液具体运筹小麦作物方面，邵文奇等学者开展了基肥、追肥生长时期不同施肥量的研究，发现沼液的施用能够促进植株株高、茎基宽的发育。另外，在保证产量的前提下，沼液追施量增加容易导致小麦生长旺盛、贪青迟熟，这与本研究单纯施用沼液BS处理组存在生长发育的问题具有一致性 [ 14 ] 。本研究同样采取“等氮”处理的思路，潘飞飞等学者设置不同沼液代替化肥比例组合，研究结果显示沼液50%代替化肥有利于冬小麦秸秆生物量和地上部总生物量，在维持现有产量水平的条件，联合配施沼液可以节省传统化肥 [ 15 ] 。本研究中，S-BS亦与CF处理组在理论产量上无明显差异，证明了沼液可以取代部分化肥，达到“减肥”的效果。究其原理方面，王静童等学者认为短期沼液施用能够促进小麦对氮磷钾营养元素的吸收利用，缓解过酸性土壤对作物的危害，需要主要盐碱的蓄积，并提出适宜沼液的施用量 [ 16 ] 。尤其秸秆还田处理条件下，沼液中丰富的生物炭联合施用可以降低土壤反硝化(Den)和硝态氮(NO₃−-N)氨化(DNRA)速率，更有利于沼液替代化肥处理土壤DNRA过程，促进土壤N素的保存和减少土壤N₂O的损失，提高作物的生产产量 [ 17 ] 。

土壤微生物群落是土壤生态系统中重要的部分，群落的丰富度和多样性进一步反映土壤的肥力和作物生长情况，多样性指数越高表明土壤微生物群落的丰富度和多样性越高 [ 18 ] 。真菌作为土壤中数量仅次于细菌的微生物，既有益生菌，又有致病菌 [ 19 ] 。本研究中，不同处理对真菌相对丰度影响差异明显，BS、CF单独施肥处理显著减少物种的数目，S-BS和S-CF联合施肥同样减少物种的数目，说明沼液、化肥、秸秆的联合施用反而会降低群落中物种的数量，结论与刘东海等学者的研究结果类似，单独施用化肥降低土壤pH以及真菌群落的丰富度和多样性，生产上还要适当降低有机肥替代化肥比例才能保障产量 [ 9 ] 。另外，不排除沼液中不同程度的水分胁迫，改变了根际微生物群落结构及多样性，包括对不同生物酶的影响 [ 20 ] 。秸秆耦合方面，马垒等学者通过长期秸秆还田对潮土真菌群落研究发现，能够调节真菌种间互作关系，刺激特定物种集群生长，抑制病原菌活性，能够实现高产量 [ 21 ] 。吕金玲亦发现深松秸秆还田能够提高了小麦灌浆期的真菌丰富度，提高团聚体的稳定性，增加土壤中多糖等物质含量，提高小麦的产量，适合特定地区推广应用 [ 22 ] 。