补充灌溉而无养分替代条件下小粒咖啡（Coffea arabica）幼苗早期生长受限：基于FarmBot平台的试验研究

《Sci》：Lower Early Growth of Coffea arabica Seedlings Under Supplemental Irrigation Without Nutrient Replacement: A FarmBot Trial

　　

机器人平台可在咖啡早期试验中标准化处理施予，但关于水—养分联合响应的证据仍然有限。研究人员评估了补充灌溉与施肥对小粒咖啡（Coffea arabica） cv. Castillo幼苗的主效应及交互效应，开展了一项为期180天的户外FarmBot试验。66株幼苗被分配至对照组（A区，n = 18）、仅灌溉组（B区，n = 18）、仅施肥组（C区，n = 15）以及灌溉加施肥组（D区，n = 15）。监测项目包括株高和叶直径，以及土壤湿度和温度。9株幼苗未存活至第6个月评估；在剩余57株中，另有3个株高值和4个叶直径值因与活体植株尺寸不符而被编码为零并视为缺失；分析样本为株高n = 54、叶直径n = 53。第6个月时，B区均值最低（31.0 ± 5.3 cm；30.9 ± 12.1 mm），而C区和D区数值均值最大。单因素方差分析检测到处理间差异：株高F(3, 50) = 14.81, p < 0.001；叶直径F(3, 49) = 20.26, p < 0.001。Tukey检验显示，B区与A、C、D区均存在差异，而A、C、D区之间差异不显著。多元方差分析检测到灌溉、施肥及交互效应；然而，回归模型仅分别解释了株高和叶直径方差的36%和43%。在这些特定场地条件下，无养分替代的补充灌溉与较差的幼苗生长相关。研究结果支持在特定场地条件下将补充灌溉与养分状况和排水状况一并评估。

**解读文章：基于FarmBot的补充灌溉与施肥对咖啡幼苗早期生长影响的试验研究**

**1. 研究背景与问题**
精准农业通过集成传感器、自动化及空间解析数据，旨在改善投入管理并减少操作者相关变异。然而，在拉丁美洲生产系统中，其应用仍受成本、基础设施和技术能力制约。对于咖啡而言，幼苗早期生长对水分供应和养分供给高度敏感：干旱会改变生长、光合作用和胁迫生理，而氮素响应则取决于土壤水分状况。现有咖啡水肥研究常因处理时间、空间坐标及重复测量依赖人工操作而难以比较。FarmBot作为一种可编程门架系统，具有固定植株坐标、可重复供水和时间关联的传感器及图像记录功能，与人工操作相比，能够在环境条件下实现标准化处理施予和可追溯的重复观测。因此，有必要在标准化小区平台上评估补充灌溉与施肥及其交互作用对咖啡幼苗形态指标的影响。

**2. 研究目的与设计**
研究人员提出以下问题：补充灌溉、施肥及其交互作用是否改变两个预设形态学结果——株高和叶直径——在小粒咖啡 cv. Castillo生长前180天内的表现。假设施肥可调节补充灌溉的响应，使联合响应不能仅由单一投入效应解释。为此，采用FarmBot实施的2×2析因设计，估算水分和养分的单独主效应及其交互作用。试验地点位于哥伦比亚桑坦德省El Páramo的一个18 m²露天实验区，66株幼苗按A（对照，n=18）、B（仅灌溉，n=18）、C（仅施肥，n=15）、D（灌溉+施肥，n=15）四个处理区分配。B区每两天接受约96 mL补充水分；C区每月每株施用10 g磷酸二铵（DAP，18–46–0）和12 g尿素（46–0–0），总施肥量为每株22 g，持续6个月；D区同时接受两种程序化输入。由于试验地无遮挡，程序化水分输入被定义为补充灌溉。第6个月时，9株死亡，另有部分零编码无效值被剔除，最终分析样本为株高n=54，叶直径n=53。

**3. 主要研究方法**
- **机器人平台系统**：使用FarmBot Genesis XL V1.7（FarmBot Inc., San Luis Obispo, CA, USA），具备可编程门架、固定植株坐标、可重复供水及时间关联的传感器和图像记录功能，减少操作者相关变异。
- **统计分析方法**：采用单因素方差分析（ANOVA）比较第6个月处理均值；析因方差分析和多元方差分析（MANOVA）评估灌溉、施肥及交互效应；Tukey HSD进行成对比较；探索性回归分析生长与环境因素关系；残差正态性用Shapiro–Wilk检验，方差齐性用Levene检验。
- **数据筛选规则**：零编码的株高或叶直径值被视为无效占位符而非生物学零值，因存活幼苗不可能有零尺寸；该规则独立于处理分配，统一应用于所有区域。

**4. 研究结果**
- **描述性统计与生长曲线**：第6个月，D区和C区数值均值最大（D区：株高48.7 ± 9.4 cm，叶直径68.8 ± 11.8 mm；C区：株高47.2 ± 8.2 cm，叶直径64.9 ± 15.5 mm），A区居中（43.7 ± 8.7 cm；55.7 ± 15.1 mm），B区最低（31.0 ± 5.3 cm；30.9 ± 12.1 mm）。生长曲线显示，处理差异随时间增大；A、C、D区株高轨迹相似，B区自第2个月起生长较低，至第6个月差距显著扩大。
- **单因素方差分析**：第6个月，株高F(3, 50) = 14.81，p < 0.001；叶直径F(3, 49) = 20.26，p < 0.001，表明处理间存在显著差异。残差诊断支持模型假设（株高：Shapiro–Wilk W = 0.979，p = 0.44；Levene F(3, 50) = 1.95，p = 0.13；叶直径：W = 0.983，p = 0.66；Levene F(3, 49) = 0.58，p = 0.63）。
- **Tukey HSD成对比较**：B区与A、C、D区在株高和叶直径上均显著不同（调整后p < 0.001）；A、C、D区之间无显著差异（调整后p在0.100至0.968之间），且这些对比的置信区间均包含零，因此数值排序不能解释为施肥或联合处理优于对照。
- **交互作用图与MANOVA**：交互作用图显示非平行轮廓，提示灌溉×施肥交互；MANOVA检测到灌溉（Wilks’ λ = 0.740，F = 9.15，p < 0.001）、施肥（Wilks’ λ = 0.665，F = 13.09，p < 0.001）及交互效应（Wilks’ λ = 0.855，F = 4.42，p = 0.017）均显著。施肥的Wilks’ λ小于灌溉，表明肥料因子在数据中产生更大的多元分离。
- **回归分析**：回归模型解释了株高方差的36%（F(2, 51) = 14.22，p < 0.001）和叶直径方差的43%（F(2, 50) = 18.75，p < 0.001），解释力有限；模型系数显示灌溉项为负、施肥项为正，与区域均值一致，但被视为探索性描述而非确认性估计。

**5. 讨论与结论**
B区的低值表示在试验条件下特定处理关联，可能的解释包括无养分替代下的养分稀释或淋失、过量水分下根区通气性降低以及幼嫩咖啡植株对施用时间的敏感性。但研究人员强调这些都是需后续实验检验的合理推测，而非已证实的因果机制。实际意义方面，补充灌溉应结合养分状况、排水和微气候一起评估，而非作为孤立的日历输入。试验局限性包括：单一地点、单一季节、一个品种、区域大小不等；场地无遮蔽且场地表征仅限于平台自身的湿度和温度序列，降雨、排水、土壤类型和初始养分状况未纳入协变量；肥料配方固定且仅一种施用量；分析为组级设计，损失和零编码分布不均匀。未来研究应跨土壤、降雨、季节和基因型重复析因设计，校准水肥梯度，并记录气孔导度、SPAD叶绿素指数、叶水势、根冠比、根生物量、组织养分和排水量，并采用混合效应模型、随机森林和梯度提升树进行预测。结论表明，水分可用性和施肥在小粒咖啡幼苗期早期生长中显著交互；未补充养分的补充灌溉与第6个月较低的株高和叶直径相关，而对照、仅施肥和灌溉+施肥三种处理统计上相当。析因分析检测到灌溉×施肥交互，支持补充水分对早期生长的响应随养分状况而变化。FarmBot提供了可重复的处理施予和空间关联监测，但更广泛的验证需在不同条件下进行。