《Agronomy》：Data-Driven Framework Integrating Database Analysis and Bayesian Segmented Quantile Regression to Determine Critical Nutrient Levels in Soil and Leaves of Peach Trees
编辑推荐：
施肥推荐常基于少数地区、少量品种的小规模校准试验提出，导致养分临界水平（CL）与充足范围（SR）未必适用，且CL与SR是否因桃树产量与品质指标（如果肉硬度、总可溶性固形物（TSS）浓度）而异尚不清楚。研究人员将数据库与贝叶斯分段分位数回归（BSQR）结合，以估
施肥推荐常基于少数地区、少量品种的小规模校准试验提出，导致养分临界水平（CL）与充足范围（SR）未必适用，且CL与SR是否因桃树产量与品质指标（如果肉硬度、总可溶性固形物（TSS）浓度）而异尚不清楚。研究人员将数据库与贝叶斯分段分位数回归（BSQR）结合，以估算土壤与叶片养分的CL与SR。研究旨在针对巴西南部亚热带气候下桃树，融合数据库与BSQR模型，提出土壤及叶片养分CL与SR。数据集含巴西南部南里奥格兰德州208个观测，涉及‘Maciel’与‘Chimarrita’桃（Prunus persica）品种。基于BSQR建立平台回归模型，量化果实产量、品质与土壤/叶片养分浓度的关联。结果显示：叶片氮（N）、钾（K）针对果实品质（果肉硬度、TSS）的CL与SR低于针对产量者；研究人员首次发现叶片N、K与产量相关的CL与SR在品种间存在差异，不同土壤类型与气候站点间叶片N的CL与SR亦不同，低有机质砂土站点叶片N（产量相关）CL更高。该研究对巴西南部桃产业有推进作用，所提“数据库+BSQR”范式对全球果园具应用价值。
研究背景方面，巴西为拉美第三大桃生产国，但平均单产仅约欧美主产国一半，原因在于区域养分推荐依赖少数校准试验，数据库覆盖品种、砧木、产区有限，且CL与SR是否因产量与品质指标（果肉硬度、TSS）而分异、是否随品种与地点变化均属空白。Boundary Line方法以群体最优表现为边界线，常配线性平台（linear-plateau）分段回归定CL；BSQR可在宽环境变异下给出CL与SR的后验估计，故研究人员以“数据库+BSQR”框架填补上述缺口。研究目的为评估桃品种与地点对产量、硬度、TSS方差的影响，并借BSQR推求土壤与叶片养分CL与SR；假设为品质相关N、K的CL与SR低于产量相关者，且CL与SR随品种与站点变化。
关键技术方法上，研究人员构建观察性数据库（非控制试验）208株桃树观测，含‘Maciel’（106）与‘Chimarrita’（102），来自南里奥格兰德州Bento Gon?alves、Pelotas、Porto Alegre三产区（2009–2010季），砧木含Aldrighi等五种且跨区分布以减混杂；土壤取0–20 cm冠层投影样，Mehlich-1提P、K，KCl提Ca、Mg，SMP法测交换酸；叶片花后13周取中冠层百叶，凯氏定N，ICP-OES定P、K、Ca、Mg；果实计产并测硬度与TSS。统计上先描述统计与ANOVA（仅描变异），PCA看多变量结构，主体用BSQR：相对产量（%）对养分浓度做分段平台贝叶斯推断，MCMC/Gibbs四链各1万自适应+1万抽样，Rhat近1收敛，CL取平台拐点，SR取90%后验可信区间；并按全库与“品种/站点”分层估CL/SR，Student t检验比CL差异。
研究结果部分，3.1探索性数据库分析显示单株产量受“品种×站点”交互显著，Pelotas最高（Chimarrita 13.01、Maciel 14.46 kg/株），Porto Alegre内Maciel（7.26）高于Chimarrita（2.84）；TSS仅受站点影响（Porto Alegre最高），硬度受交互影响（Porto Alegre的Chimarrita最高）。PCA（PC1 48%+PC2 23%=71%）示产量与冷量、叶N/P正关，与辐射、降雨、盐基饱和度负关；Bento Gon?alves组受黏土、有机质、盐基饱和度驱动，Pelotas+Porto Alegre组受叶N、产量、TSS、硬度驱动，砧木无分离效应。3.2叶片养分临界浓度中，通用法（n=208）得叶N产量CL 26.0 g kg?1、SR 24.7–28.4，叶K产量CL 28.4（Maciel）/25.9（Chimarrita）g kg?1且品种间SR异；站点上Pelotas+Porto Alegre叶N产量CL 34.5 g kg?1显著高于Bento Gon?alves 24.4 g kg?1，叶K站点间CL无差。通用叶P CL 2.4（SR 1.9–2.9）、Ca CL 13.1（SR 12.7–13.5）、Mg CL 3.1（SR 1.0–4.7）g kg?1。针对硬度：叶N CL 26.0（SR 24.7–28.4）、P CL 3.1（SR 2.3–4.0）、K CL 21.1（SR 15.6–23.1）、Mg CL 2.2（SR 0.5–4.1）g kg?1，Ca无CL仅SR 4–19.3。针对TSS：叶N CL 27.3（SR 25.4–30.2）、P CL 1.5（SR 0.8–1.7）、K CL 21.2（SR 19.9–25.2）、Ca CL 16.5（SR 14.2–19.1）、Mg CL 2.0（SR 0.8–2.5）g kg?1。总体叶N/K/Mg产量CL高于品质（硬度、TSS）CL；与南里奥格兰德/圣卡塔琳娜州CQFS-RS/SC推荐比，BSQR的N/Ca/Mg充足范围偏低、K临界水平偏高。3.3土壤养分临界水平中，两品种合并估得土P（Mehlich-1）产量CL 23.8 mg dm?3（SR 21.4–26.2）、K CL 190 mg dm?3（SR 187–192）、Ca CL 5.8 cmolc kg?1（SR 4.5–6.8）、Mg CL 1.2 cmolc kg?1（SR幅度宽，峰值0–2）、盐基饱和度CL 65%（SR 62–67%）；土壤养分与硬度/TSS无可用BSQR边界故未估CL。
讨论部分总结：观察库虽为2009–2010年数据，但环境变异宽，足以支撑养分–表现基本关系阈值推断，现代 orchard需未来新库验证。品质相关N、K的CL低于产量相关，反映目标差异与气候/品种/物候影响；高叶N反伴低TSS可能源于营养生长稀释糖分（葡萄等前例支持，但本研究未直测生理故为解释而非结论）。叶K高伴品质优合K促糖转运渗透调节之说。气候中冷量正关产量、负关硬度/TSS，降雨高（Bento Gon?alves）站点产量/品质偏低或因淋溶与病害。品种间叶N CL差异或关Chimarrita长势旺、叶库大；Maciel叶K CL高或关其高产致果采携出K多。Pelotas+Porto Alegre叶N CL高于Bento Gon?alves，或关砂土低有机质需更高叶N保产。土P CL 23.8 mg dm?3较CQFS高约15%，K CL 190 mg dm?3落其CECpH7 15.1–30.0 cmolc kg?1档内，盐基饱和度CL 65%可作参。BSQR自宽观察库所得值异于小样本旧标准，应视作巴西南部商业桃园区域校准参。
结论部分翻译：桃产量与品质响应受品种与站点共同影响；单株产量与果肉硬度受“品种×站点”交互显著，TSS主要受站点条件影响。大观察库+BSQR可有效定桃园养分CL与SR，覆盖叶N/P/K/Ca/Mg（产量与品质）及土P/K/Ca/Mg/盐基饱和度（产量）。假设获证——品质（硬度、TSS）相关叶N、K的CL与SR一般低于产量相关者，故营养靶标应随“增产”或“提质”而调。另一假设亦获证：叶N、K产量CL与SR在‘Chimarrita’与‘Maciel’间及对比土壤气候站点间分异，低有机质砂土站点叶N产量CL更高。研究强调商业桃园养分标准与施肥推荐须含品种特异与站点特异条件。
（全文登于《Agronomy》）
