中国农业区秸秆还田对作物产量与土壤属性的综合效应：基于全国监测网络的定量评估

《Land》：Benefits in Crop Yield and Soil Properties of Straw Return Across China’s Agricultural Regions

　　

摘要专业翻译 已发表的综合研究（syntheses）表明，秸秆还田（straw return）总体上有利于作物生产和土壤质量。然而，仍需直接证据来确定这些效益是否具有一致性，以及其在不同区域和种植制度间的变异规律。本研究分析纳入了一个覆盖中国六大主要农业区的全国性田间试验网络，共24个监测站点、379组处理-对照记录。具体而言，共汇编了2800个结局特异性效应量（outcome-specific effect sizes），每个效应量均以秸秆还田处理与其配对秸秆移除对照的比值表示，涵盖作物产量及九项土壤指标。与秸秆移除相比，秸秆还田与作物产量提高3.6%（95%置信区间，0.9%至6.5%）相关；汇总估计值在以下指标中亦较高：土壤有机质（soil organic matter）提高11.5%、全氮（total nitrogen）提高5.7%、速效磷（available phosphorus）提高9.0%、速效钾（available potassium）提高14.7%、体积土壤含水量（volumetric soil water content）提高6.1%、水稳性团聚体（water-stable aggregates）提高17.8%、耕作层厚度（plough-layer thickness）提高16.4%。作物产量、速效磷、速效钾、水稳性团聚体和耕作层厚度在不同农业区域间检测到异质性；由于在该网络中农业区域与种植制度高度混杂，这些差异是区域关联性的而非可归因于地理本身，且若干区域分层仅由两至三个监测站点代表。在种植制度间，速效钾、水稳性团聚体和耕作层厚度检测到异质性，但作物产量未检测到异质性。较短持续期类别具有较大的速效磷响应，而较长持续期类别具有较大的团聚体和耕作层响应，但未显示额外的平均产量效益。所观察到的响应特征提示，中国南方稻田的磷钾平衡与土壤结构监测应作为优先事项。类似地，东北地区的残茬处置（residue placement）、华北地区基于土壤测试的钾素管理，以及西北和西南地区的水肥协调管理，均被指示为监测优先事项而非已验证的处方。总体而言，秸秆还田与较高的作物产量及多项土壤属性值相关，尽管响应特征在不同区域间存在差异。这些模式表明，若要实现秸秆还田在作物产量和土壤质量方面的全国性效益，需要经本地验证的区域特异性策略。

论文解读

一、研究背景与问题

秸秆还田作为中国农田生态系统中的重要残茬管理措施，通过将生物质和养分在田间作物系统内循环利用，对维持土壤肥力和促进农业可持续发展具有关键作用。已有综合研究普遍报道秸秆还田后作物产量有所提高，但增产幅度因作物类型、管理措施和持续年限而异。然而，在全国尺度上，秸秆还田对作物生产和土壤属性的影响尚未被充分理解。核心实践问题不仅在于秸秆还田平均而言是否有益，更在于其效应在不同生产条件下如何变化。已有定量证据大多来源于独立的田间试验及其后续综合研究，但不同研究在指标选择、采样方案、测定方法和观测周期上的差异，使得将点位尺度的发现整合为全国性一致性评估存在较大困难。基于模型的评估虽可扩展空间覆盖范围，但其估计对数据输入、参数化和分析假设较为敏感。在大规模、长期、统一指标框架下的田间观测数据亦十分稀缺。因此，建立协调统一的监测网络，为比较不同区域间作物和土壤响应特征提供基于田间观测的互补性视角，具有重要的科学和实践意义。

二、研究设计与核心发现

研究人员利用一个前瞻性协调的全国监测网络数据集，该网络于2021年建立，共包含47个监测站点，其中24个站点报告了符合配对秸秆还田与秸秆移除对照设计要求的作物产量和土壤指标数据。最终数据集包含来自24个站点的379组处理-对照记录，覆盖中国东北、华北、长江中下游、华南、西北和西南六大农业区，代表水稻基系统、玉米基系统、小麦-玉米轮作、小麦单作和棉花系统五类主要种植制度。与以往基于独立发表实验的荟萃分析不同，该数据集具有三个独特特征：其一，记录无论观测效应的方向和统计可检测性均纳入数据集，不存在发表偏倚；其二，所有处理-对照对均满足严格单因素准则，即处理与对照仅在是否还田上存在差异，土壤类型、气候和施肥等点位属性在比较的两臂间共享；其三，十项结局在同一站点、采用统一测定方法同步测定，可在站点内描述多指标响应特征。

研究得出以下主要结论：网络尺度上，秸秆还田与作物产量提高3.6%相关，同时土壤有机质提高11.5%、全氮提高5.7%、速效磷提高9.0%、速效钾提高14.7%、体积土壤含水量提高6.1%、水稳性团聚体提高17.8%、耕作层厚度提高16.4%。容重平均响应为?2.7%，土壤pH无检测到的平均响应。区域关联性异质性在作物产量、速效磷、速效钾、水稳性团聚体和耕作层厚度五项结局中被检测到；种植制度间异质性在速效钾和两项结构指标中被检测到，但作物产量未检测到差异。较短持续期类别的速效磷响应较大，而较长持续期类别在团聚体和耕作层响应上较大，但未显示额外的平均产量效益。该研究的重要意义在于，通过多站点统一指标框架的定量评估，揭示了秸秆还田的效益是多维度的——不仅限于产量，还涵盖土壤养分、水分和结构等多个方面，且这些响应在不同区域和种植制度间存在显著差异，为制定区域适应性秸秆还田策略提供了基于多站点观测的证据基础。

三、关键技术方法

研究人员采用的主要技术方法包括：第一，基于前瞻性协调监测网络的数据采集框架，通过年度技术规程统一各站点的指标定义、采样深度和测定方法，确保跨站点可比性；第二，效应量计算采用自然对数响应比（ln RR），以处理均值与对照均值之比的对数表示，消除不同指标单位差异，使各结局可在共同的无量纲尺度上比较；第三，抽样方差根据处理和对照的标准差及样本量计算，缺失离散度信息时按均值的10%比例插补；第四，总体效应采用限制最大似然（REML）估计的随机效应模型进行综合，权重为抽样方差与比较间方差之和的倒数；第五，亚组分析以农业区域、种植制度和还田持续期为分组变量，采用REML综合检验评估组间异质性。数据来源为全国监测网络中的24个监测站点。

四、研究结果

网络覆盖与总体响应特征。 379组记录覆盖全部六个农业区和五类种植制度，但覆盖不均衡，长江中下游和东北贡献216组记录，西北和西南分别仅34组和32组。网络尺度上，秸秆还田将产量提高3.6%（p=0.009），土壤有机质、全氮、速效磷和速效钾均显著提高，速效钾以14.7%的响应幅度居养分指标之首；体积土壤含水量、水稳性团聚体和耕作层厚度亦显著提高，两项结构指标为十项结局中汇总响应最大者。容重和pH无显著平均响应。

区域关联性异质性。 五项结局检测到区域间异质性。作物产量响应在区域间存在差异（p=0.007），长江中下游和华南为正估计，其余区域置信区间包含零。速效磷在区域间变异显著（p<0.001），长江中下游、华南和西南为正估计。速效钾显示最强的区域关联性异质性（p<0.001），长江中下游、华南和西北为正估计，而华北基于两个站点为负估计。水稳性团聚体在区域间有差异（p=0.004），长江中下游和华北为正估计。耕作层厚度亦存在区域差异（p<0.001），长江中下游和东北为正估计。土壤有机质、全氮、体积土壤含水量、容重和pH未检测到区域间异质性。

种植制度响应特征。 种植制度间异质性在速效钾和两项结构指标中被检测到，但作物产量未检测到差异（p=0.203）。水稻基系统内产量估计为6.3%，玉米基系统为?0.5%，小麦-玉米轮作为?0.6%，但系统内估计未确立系统间产量差异。速效钾在水稻基系统为24.5%，玉米基系统为7.0%，小麦-玉米轮作为?3.1%（p<0.001）。水稳性团聚体响应在三个可估计系统间有差异（p=0.016），分别为23.1%、?0.5%和14.9%。耕作层厚度亦在系统间有差异（p<0.001），对应估计为29.0%、9.4%和4.4%。土壤有机质、全氮、速效磷、体积土壤含水量、容重和pH未检测到系统间异质性。

持续期类别差异。 较短持续期类别的速效磷响应为16.3%，较长持续期类别无可检测响应（?3.6%，p<0.001）。水稳性团聚体在两个持续期类别均较高，但较长持续期类别的汇总响应更大：9.0%对32.7%（p=0.023）。耕作层厚度亦在类别间有差异，估计分别为5.1%和38.3%（p<0.001）。作物产量、土壤有机质、全氮、速效钾、体积土壤含水量、容重和pH未检测到持续期间异质性。

作物情境间的产量响应变异。 探索性分析提示11种响应作物情境间存在产量响应异质性（Q=18.64，df=10，p=0.045），早稻、晚稻和油菜的正估计置信区间排除零，其余情境区间包含零。该结果描述作物情境间的变异而非作物物种间的两两差异。

五、讨论与结论

讨论部分指出，网络尺度上的响应特征将3.6%的产量增加与土壤肥力、水分和结构的更广泛变化相结合。产量估计小于既往文献报道值，可能反映了本网络所代表的特定站点和管理对比组合，且网络估计不设发表偏倚过滤。速效钾具有最大的养分响应，秸秆钾素释放迅速，应纳入田间钾素预算。速效磷的增加幅度与一至数年秸秆磷素还田量相当，可通过难溶性磷库的再分配实现。区域关联性响应特征提示南方稻田应优先关注磷钾平衡与结构监测，东北地区应关注残茬处置和土壤条件，华北小麦-玉米轮作应保留基于土壤测试的钾素管理，西北地区应将残茬保留与水分管理相协调，西南地区应将秸秆与轮作养分管理相结合。持续期结果提示较长期还田在结构指标上的响应更大，但速效磷响应较小，不支持渐进累积假设。研究局限性包括：站点记录数不均衡、管理信息不完整、区域与种植制度高度混杂、部分区域仅由两至三个站点代表、持续期类别为不同记录组而非同一田块的重复观测。未来监测应扩展和平衡站点覆盖、完善管理和环境数据、尽可能纳入同一田块的重复测定。

研究结论为：在24站点网络中，秸秆还田与3.6%的平均产量增加及七项土壤养分、水分和结构指标的较高汇总值相关。区域关联性异质性具有结局特异性：长江中下游（5站点）和华南（3站点）呈现最清晰的正向产量和磷钾信号；东北站点在土壤有机质和耕作层上的信号比产量更清晰；华北两个站点的速效钾估计为负；西北三个站点的主要正向信号为速效钾；西南两个站点呈现正向速效磷响应。由于区域与种植制度在该网络中混杂，这些对比不能归因于地理本身，且由一至两个站点代表的分层估计受站点限制。主要种植制度未分化产量响应，较长持续期类别在代表记录中未显示额外平均产量效益。管理应保持区域适应但范围宽泛：南方水稻系统应强调磷钾平衡和结构监测；东北系统应强调残茬处置和土壤条件；集约化小麦-玉米轮作应保留土壤测试钾素管理；水分受限系统应将还田与土壤水分保持相协调；多样化轮作应将秸秆与作物特异性养分管理相结合。在所有情境中，产量应与速效磷、速效钾及至少一项结构指标共同评估。因此，需要更具区域适应性的秸秆还田措施，以实现秸秆还田的农艺效益最大化。