宽窄行垄作配置对大豆冠层结构、光合特性、根系形态及产量的影响——基于黑龙江三个地点三年田间试验

《Life》：Astrogliosis Drives Abnormal Neuronal Tau Assimilation and Increase in Post-Synaptic Levels in a Model of Alzheimer’s Disease

　　

在共同目标密度下，种植配置可改变大豆（soybean）冠层结构与产量。研究人员于2022—2024年在哈尔滨、绥化和五常对品种黑农531的四种配置进行评价，每地点-年份设三个随机完全区组，目标密度为340,000株 ha?1。数据集包含15个性状各108个小区观测值。固定环境RCBD模型纳入地点、年份、配置及其交互作用，以及嵌套于地点-年份内的区组。Tukey比较按性状和分析时期分别进行。配置影响合并籽粒产量（F3,54 = 3.516，p = 0.021）。A、B、C和常规D的平均产量分别为4020.0、4007.6、4042.0和3854.9 kg ha?1。C比D高187.1 kg ha?1（4.85%；校正p = 0.026）；A–D和B–D差异不显著。2022年，A和C高于D，而2023年或2024年年度两两产量差异均不显著。B和C的合并叶面积指数（LAI）、冠层透射率和下部叶片光合作用高于D；上部叶片光合作用、电子传递速率（ETR）和百粒重无显著差异。根系性状差异顺序为C > A > B > D。部分生理性状的显著交互作用和Fv/Fm残差离散度不等的证据限制了对合并模式的统一解释。研究结果支持在测试环境中存在配置效应，C相对D表现出合并产量优势；相关性状响应并未确立因果中介关系或超出采样条件的表现。

研究背景方面，提高大豆（soybean）生产力可能依赖于协调冠层发育、生殖生长和地下资源获取的管理措施。种植配置可在不改变群体密度的前提下改变植株空间分布，但产量对行距的响应在不同生产环境中存在差异，已有研究报道窄行大豆产量响应因地区而异，试验优势未能在生产者田间稳定重现，因此种植配置应在预期环境和管理条件下评价。种植几何响应的生理基础与冠层发育和光截获密切相关，大豆窄行的产量优势与营养生长和早期生殖生长期间较大的光截获有关，并伴随有效节数和每有效节荚数的变化；冠层内光合活性的分布也值得关注，上下部冠层叶片测定可揭示生理差异是贯穿冠层还是局限于特定叶位。地下测定提供另一视角，根系长度、直径、表面积和体积描述与资源获取相关的形态特征，但根系性状与表现的关系取决于资源条件，并非仅由根系大小决定。基于此，研究人员于2022—2024年在黑龙江省哈尔滨、绥化和五常评价春大豆品种黑农531的四种垄作配置，其中两种宽窄行四行配置和一种行内等距四行配置设于1.3 m垄上，与0.65 m垄上常规两行配置在共同目标密度340,000株 ha^?1下比较。研究目的为：评估三个地点年度籽粒产量对种植配置的响应；表征可获得的冠层、产量构成、光合、荧光和根系测定指标；探索所记录性状间的生理关联。研究检验对比性空间配置是否产生可检测的年度产量差异，并考察其性状谱能否识别值得进一步田间验证的配置。

关键技术方法方面，研究采用多点多年随机完全区组设计（RCBD），每地点-年份设三个田间区组，共108个小区；用联合固定环境RCBD模型分析地点、年份、配置及其交互作用，并以Tukey学生化极差法进行多重比较；用LI-COR LAI-2200C植物冠层分析仪测定叶面积指数（LAI）和冠层透射率；用LI-6800便携式光合系统测定上、下部叶片气体交换参数；用FluorPen FP100荧光仪测定暗适应和光适应叶绿素荧光参数并计算电子传递速率（ETR）；用Epson Expression 12000XL扫描仪结合WinRHIZO Pro 2019a分析根系形态。供试品种黑农531由黑龙江省农业科学院大豆研究所提供，试验地点为哈尔滨、绥化和五常。

研究结果方面，3.1 跨地点和年份的合并籽粒产量响应：联合分析包含108个年度小区观测值，籽粒产量在配置间（F₃,₅₄ = 3.516，p = 0.021）、地点间（F₂,₅₄ = 6.099，p = 0.004）和年份间（F₂,₅₄ = 3.971，p = 0.025）存在差异；涉及配置的产量交互作用均未达到0.05名义阈值；A、B、C、D合并均值分别为4020.0、4007.6、4042.0和3854.9 kg ha^?1，模型标准误为45.5 kg ha^?1；相对D，A、B、C分别增加165.2 kg ha^?1（4.28%）、152.7 kg ha^?1（3.96%）和187.1 kg ha^?1（4.85%）；仅C–D比较经Tukey校正后显著（同时95%置信区间16.5–357.7 kg ha^?1；校正p = 0.026），A–D和B–D未达阈值；A、B、C之间差异小且不确定，紧凑分组为A = ab、B = ab、C = a、D = b，因此显著合并对比为C–D，数值排序C > A > B > D并非完全解析的统计排序。3.2 年度产量对比和地点剖面：全部配置和地点2022、2023、2024年平均产量分别为4037.7、4014.2和3891.4 kg ha^?1；年度模型中配置效应2022年显著（F₃,₁₈ = 3.864，p = 0.027），2023年（p = 0.313）和2024年（p = 0.082）不显著；2022年A、B、C、D均值分别为4134.1、4037.8、4145.5和3833.5 kg ha^?1，A高于D 300.6 kg ha^?1（校正p = 0.044），C高于D 312.0 kg ha^?1（p = 0.035），B–D差异不显著；2023年A数值最高但四种配置共享一个字母，2024年C数值最高但所有年度两两区间均包含零；合并地点均值哈尔滨4087.6、绥化3959.1、五常3896.7 kg ha^?1；C相对D的地点-年份平均响应从哈尔滨2023年的?5.15%到哈尔滨2022年的+10.60%，C在九个环境中七个高于D。3.3 总体性状响应和环境百分比：15个合并性状中12个检测到配置效应，上部叶片光合作用（F₃,₅₄ = 0.575，p = 0.634）、ETR（F₃,₅₄ = 1.033，p = 0.385）和百粒重（F₃,₅₄ = 0.784，p = 0.508）例外；响应图谱包含九个环境内378个非产量百分比；LAI相对D的环境平均响应A、B、C分别为+1.50%、+8.19%和+8.33%，冠层透射率分别为?5.25%、+30.41%和+28.32%，下部叶片光合作用分别为+10.58%、+54.15%和+50.45%，上部叶片分别为+2.61%、?0.85%和+2.51%；B和C在每个环境中透射率和下部叶片光合作用均数值高于D，ETR响应符号不一；三种替代配置在全部九个环境中根长、直径、表面积和体积均数值高于D。3.4 冠层结构和叶片光合作用：LAI存在配置差异（F₃,₅₄ = 10.394，p < 0.001），A–D合并均值分别为4.581、4.887、4.891和4.515，B和C分别高于D 0.372和0.376（校正p < 0.001），B和C无差异，A和D形成较低组；LAI涉及配置的交互作用均不显著；冠层透射率存在配置效应（F₃,₅₄ = 265.624，p < 0.001），A–D合并分数为0.03176、0.04376、0.04295和0.03353，B和C分别高于D 0.01023和0.00941（校正p < 0.001），A低于D 0.00177（校正p = 0.010），B和C无差异，地点×配置交互作用显著（F₆,₅₄ = 13.521，p < 0.001）；下部叶片净光合速率（Pn）存在配置差异（F₃,₅₄ = 166.285，p < 0.001），A–D合并均值分别为2.647、3.700、3.619和2.403 μmol CO₂ m^?2 s^?1，A、B、C相对D分别增加0.244、1.297和1.216，B和C合并增幅为53.99%和50.59%，B和C无差异且均高于A，年份×配置和三元交互作用显著；上部叶片Pn A–D均值分别为15.888、15.392、15.880和15.535 μmol CO₂ m^?2 s^?1，所有配置共享一个合并Tukey组，配置主效应及交互作用均不显著。3.5 叶绿素荧光和电子传递：Fv/F₀存在配置差异（F₃,₅₄ = 19.771，p < 0.001），A–D合并均值为3.729、3.511、3.493和3.796，B和C分别低于D 0.285和0.303（校正p < 0.001），A–D差异不显著，A和D一组、B和C一组，年份×配置交互作用显著；Fv/Fm存在配置效应（F₃,₅₄ = 378.653，p < 0.001），合并分组顺序为D（0.791）> A（0.724）> B（0.674）> C（0.647），相对D的降低分别为8.50%、14.74%和18.17%，地点×配置、年份×配置和三元交互作用均显著；ETR A–D均值为45.818、46.849、46.825和45.355，B和C数值上分别高于D 3.30%和3.24%但差异不显著，任何年度ETR模型中配置效应均不显著，合并配置交互作用亦不显著。3.6 荚数、粒数和粒重：每株荚数存在配置差异（F₃,₅₄ = 18.980，p < 0.001），A–D合并均值为49.05、53.24、49.34和44.11，B最高且高于A和C，A和C无差异且均高于D，相对D分别增加11.21%、20.70%和11.86%；每株粒数存在配置差异（F₃,₅₄ = 16.522，p < 0.001），A–D均值为144.49、156.34、162.66和138.69，B和C形成较高组、A和D形成较低组，B和C相对D分别增加12.73%和17.28%，无配置交互作用显著，年度配置效应每年均检测到，C在三个年度组中均高于D而B仅在2024年高于D；百粒重A–D合并均值为21.413、21.048、21.357和21.356 g，所有配置共享一个字母，年度配置检验每年均不显著。3.7 根系形态和年度根系对比：四个根系性状均显示合并配置效应，根长（F₃,₅₄ = 90.910）、平均直径（F₃,₅₄ = 113.511）、表面积（F₃,₅₄ = 199.092）和体积（F₃,₅₄ = 364.231）均p < 0.001，每个性状四个合并均值均分为独立组，顺序为C > A > B > D，任何涉及配置的交互作用均不显著；根长A–D均值为368.75、344.29、401.58和313.64 cm plant^?1，相对D分别高17.57%、9.77%和28.04%；平均直径为0.297、0.268、0.333和0.245 mm，相对D分别高21.43%、9.55%和35.95%；表面积均值为50.950、44.191、59.108和36.553 cm² plant^?1，体积为2.685、2.189、3.110和1.723 cm³ plant^?1，相对D表面积增加39.39%、20.89%和61.70%，体积增加55.83%、27.00%和80.46%；年度综合配置效应在每年四个根系性状中均显著，但并非每个合并两两分离均在年度重现。3.8 主成分、相关和根系几何：主成分分析（PCA）基于36个完整地点-年份-配置聚合剖面，PC1解释44.37%、PC2解释27.75%，两成分合计72.11%，PC3增加15.51%使三成分合计87.62%；PC1对下部叶片Pn有正特征向量系数（+0.553），对Fv/F₀（?0.600）和Fv/Fm（?0.558）有负系数；PC2主要与ETR（+0.707）和上部叶片Pn（+0.635）相关；聚合Spearman矩阵显示Fv/F₀与Fv/Fm正相关（ρ = 0.649），下部叶片Pn与两个荧光比值负相关（ρ = ?0.634和?0.584）并与ETR弱正相关（ρ = 0.289），上部叶片Pn与下部叶片Pn、Fv/F₀、Fv/Fm关联微弱并与ETR弱正相关（ρ = 0.288）；根表面积与长度比A–D均值为0.13852、0.12851、0.14761和0.11673 cm，描述顺序为C > A > B > D，体积与表面积比均值为0.05281、0.04962、0.05267和0.04724 cm，A和C数值相近且A略高。3.9 残差诊断：合并产量残差无非正态证据（Shapiro–Wilk W = 0.983，p = 0.193）和离散度不等证据（中位数中心化Levene F = 1.160，p = 0.329）；全部60个性状-时期模型中Shapiro–Wilk p值均超过0.05（最小0.070），唯一低于0.05的名义Levene结果为合并Fv/Fm（F = 2.759，p = 0.046），提示该性状残差方差可能存在异质性，其合并方差分析和Tukey比较应视为暂定结论。

讨论部分总结：4.1 种植配置和产量响应尺度：小区水平分析支持采样环境内配置存在合并籽粒产量效应，经Tukey校正C高于D 4.85%，A和B未显示同样优势，数据支持C–D对比而非所有四行配置均有显著优势，A、B、C之间未检测到两两差异；共同目标群体将配置比较与有意密度处理分开，但不能确立实际成株数相等，因未记录最终成株数，C–D产量差异部分可能反映实际群体差异而非仅空间排列；垄宽、行数和间隙序列同时改变，试验对比无法将产量响应归因于任一几何组分。4.2 冠层和光合响应：B和C合并LAI和透射率增加得到小区水平比较支持，二者描述冠层不同方面；透射率比值仅量化记录位置的穿透，不量化季节截获辐射、生物量转化或辐射利用效率；下部叶片Pn响应显著而上部叶片无差异，与差异集中于采样下部冠层相符，但不能确立全冠层碳增益。4.3 生殖和根系响应：配置差异在荚数和粒数上比粒重更清晰，B荚数最高，B和C形成较高粒数组，粒重无差异；单株构成不能在没有最终成株密度和单位地面面积粒数时直接预测收获产量；所有根系维度差异顺序为C > A > B > D，处理比较表明记录采样尺度上的形态差异，但较大回收根系不直接证明更高效的养分或水分获取；尽管根系差异明显，B–C产量差异未检测到，提示根系大小本身可能不决定产量排序。4.4 荧光和探索性关联：Fv/F₀和Fv/Fm与数值ETR剖面响应不同，ETR本身无显著合并配置效应，B和C较低荧光比值不应描述为PSII效率的普遍改善；PCA和相关总结聚合生理结构，不提供处理分离的额外随机化检验，聚合、处理效应和共享环境变异可产生协变而无直接因果关系。4.5 年度推断和实践限制：2022年检测到A–D和C–D产量差异，2023和2024年未检测到两两差异，合并C–D对比涉及九个环境观测平均值，有限重复和不显著交互作用不能证明未测试环境中的一致效益；研究使用了60个独立Tukey族，未控制整个研究出现任何假阳性的概率；未记录最终成株数、各地点确切测定日期不同、缺乏阶段分辨生理记录、根系采样体积有限且无吸收测定，均限制解释；机械兼容性、可收获性和实施成本未评价，无法推断经济最优。

研究结论翻译：在共同目标群体下，配置在测试地点和年份内影响合并大豆籽粒产量和多数记录性状。C在合并产量上高于常规D 187.1 kg ha^?1（4.85%），而A–D和B–D对比经Tukey校正后不显著。2022年A和C高于D，2023和2024年无显著年度两两产量差异。B和C的合并LAI、冠层透射率和下部叶片Pn高于D，回收根系性状差异顺序为C > A > B > D。上部叶片Pn、ETR和粒重无显著合并配置效应。生理交互作用、Fv/Fm方差诊断和最终成株数缺失限制了推断。由于各配置在垄宽、行数和行内间距上同时不同，产量响应不能归因于任何单一几何因素，且最终成株数不可获得。研究结果支持在测试生产条件下进一步评价C，但未确立因果机制或超出这些条件的一般推荐。