《Forests》：Ecofunctional Factors Associated with Resin Yield of Bursera bipinnata in Agroforestry Systems of Tropical Deciduous Forest
编辑推荐：
来自热带落叶林（tropical deciduous forest，TDF）的本地物种双羽裂橄榄属（Bursera bipinnata）所产出的柯巴树脂（copal），是墨西哥农村社区的重要收入来源。研究人员采用主成分分析（principal componen
来自热带落叶林（tropical deciduous forest，TDF）的本地物种双羽裂橄榄属（Bursera bipinnata）所产出的柯巴树脂（copal），是墨西哥农村社区的重要收入来源。研究人员采用主成分分析（principal component analysis，PCA）、层次聚类分析和多元线性模型，评估了生态功能梯度与树脂产量之间的关联。在墨西哥莫雷洛斯州特帕尔辛戈市农林复合系统（agroforestry systems，AFS）的27株树中，评估了生态、土壤、测树和经营管理变量。Horn平行分析支持保留3个主成分，解释了63.66%的总方差。多元线性模型显著。PC1与对数转换后的树脂产量呈正相关，其特征是总树高（TH）、总真菌（TF）、土壤有机质（OM）、放线菌（AC）、有效磷（P）和全氮（TOTN）为负载荷，海拔（ALT）为正载荷。相反，PC3与对数转换后的树脂产量呈负相关，其特征是土壤容重（BD）和镁（MG）为负载荷，土壤含水量（MC）为正载荷。某站点树木的对数转换树脂产量低于参照站点，而另一站点与参照站点无显著差异。尽管聚类组观测平均树脂产量不同，但差异无统计学意义。总体结果表明，在抽样树木和研究站点内，树脂产量与整合生态、土壤、测树和经营属性的多元梯度之间存在探索性关联。
研究背景方面，柯巴树脂（copal）是墨西哥农村社区具有重要文化和经济价值的非木材林产品（non-timber forest product，NTFP），由双羽裂橄榄属（Bursera bipinnata）等物种在热带落叶林（tropical deciduous forest，TDF）及传统农林复合系统（agroforestry systems，AFS）中产出。目前在莫雷洛斯、普埃布拉和格雷罗等地，柯巴树脂采收多在无森林经营方案下开展，缺乏可持续采收技术和恢复休养期，可能危及产量与社区生计。关于树脂产量的生态、土壤、测树和经营管理驱动因素研究有限，因此研究人员开展该项研究，以识别与双羽裂橄榄属（Bursera bipinnata）树脂产量相关的生态功能因子，为可持续经营提供科学基础。论文发表于《Forests》。
关键技术方法方面，研究人员在墨西哥莫雷洛斯州特帕尔辛戈市Los Sauces集体农地设置3个采样点，每点从15株候选树中依分类、产脂史和卫生状况选9株，共27株活跃采脂树。采用主成分分析（principal component analysis，PCA）前以Z分数标准化，用Horn平行分析保留成分；层次聚类基于PC得分、欧氏距离与Ward.d2法；多元线性模型以log(产量+1)为响应变量，PC1至PC3为连续预测因子，采样点为分类固定效应，并进行残差诊断。
研究结果如下。
- 2.
Materials and Methods（材料与方法）：研究人员通过站点生态与土壤变量、单株测树与经营变量及树脂产量测定，建立27株树数据集，明确树脂产量为传统刻伤采收季总质量，并说明抽样为非概率目的抽样，推断限于抽样树木与站点。
3.1. Dataset and Variable Selection（数据集与变量筛选）：描述统计显示海拔均一，坡度、采收年限、树脂产量变异大；树脂产量均值37.56±44.76 g，分布右偏含1个零值。
3.2. Correlation Analysis（相关分析）：经Benjamini–Hochberg错误发现率（false discovery rate，FDR）校正后，有机质（OM）与全氮（TOTN）最强正相关，放线菌（AC）与总树高（TH）正相关，有效磷（P）与钾（K）、真菌（TF）正相关。
3.3. Sampling Adequacy Diagnostics for Principal Component Analysis（PCA抽样适切性诊断）：完整变量KMO（Kaiser–Meyer–Olkin）为0.342，最终保留13变量后KMO为0.612，Bartlett球度检验显著，支持后续探索性PCA。
3.4. Principal Component Analysis（PCA）：Horn平行分析保留PC1、PC2、PC3，解释63.66%方差。PC1整合树高、真菌、海拔、有机质、放线菌、磷、氮；PC2整合干径、全氮、有机质、采脂年限与磷反向；PC3为容重、含水量、镁的土壤物理水分梯度。自助重抽样显示主要载荷方向较一致但稳定性因变量而异。
3.5. Hierarchical Clustering and Multiple Linear Model（层次聚类与多元线性模型）：平均轮廓宽和Gap统计量未支持稳定聚类数，三簇Jaccard相似低，Kruskal–Wallis检验组间产量无显著差异（p=0.0533）。多元线性模型显著，R2=0.435，PC1正关联、PC3负关联树脂产量，PC2不显著；某站点产量低于参照站点S1，另一站点与S1无显著差异。
讨论部分总结：相关分析仅少数配对关系经FDR后显著，反映土壤肥力、微生物与测树属性共变。PCA三类梯度分别表征结构—地形—土壤生物、干径—有机质—采脂史、物理水分特征，但均为共变模式而非因果或生理机制。聚类似是而非，模型显示PC1和PC3与log树脂产量关联，站点差异可能来自未测环境或经营异质性。局限性包括目的抽样偏差、单株采收努力未量化、样本量小致PCA仅作探索。
结论部分翻译：双羽裂橄榄属（Bursera bipinnata）的树脂产量与抽样树木和站点内整合生态、土壤、测树和经营属性的多元梯度相关。PCA保留3个主成分，解释63.66%总多元变异。多元线性模型显示PC1与对数转换树脂产量正相关，PC3与响应变量负相关，PC2无统计学显著性。采样站点也与树脂产量相关，其中某站点数值低于参照站点S1，而另一站点与S1无显著差异。层次聚类用于探索三组分划，但验证标准显示该解支持与稳定性有限，探索组间树脂产量差异无统计学意义。总体结果为树脂产量与生态、土壤、树体结构和经营属性组合存在关联提供探索性证据，但未确立因果或生理机制。所识别的多元模式可为未来研究及双羽裂橄榄属（Bursera bipinnata）可持续采收中站点、土壤、树木和经营属性监测提供依据。
