印度尼西亚北苏门答腊安息香树脂农林复合系统的可持续集约化与生产动态

《Trees, Forests and People》：Sustainable Intensification and Production Dynamics of Styrax Resin Agroforestry in North Sumatra, Indonesia

　　

北苏门答腊的安息香树脂农林复合系统提供了重要的经济和生计效益，但关于长期生产动态和生产力模式的定量证据仍然有限。本研究考察了2010年至2021年印度尼西亚北苏门答腊六个安息香生产县在生产、生产力、增长和生产波动性方面的变化。分析结合时间序列评估、生产力测量、变异系数、复合年增长率（CAGR）和三阶多项式回归，以表征时空生产动态。结果表明，安息香树脂产量的增加主要与单产提高相关，而非种植面积扩张。尽管种植面积下降，总产量自2014年起显著增加，表明现有农林复合系统的生产效率提高和集约化增强。Toba县和Pakpak Bharat县表现出相对稳定且高效的增长，而Central Tapanuli县和Humbang Hasundutan县则表现出更大的生产波动性。North Tapanuli县和Dairi县表现出相对有限或停滞的生产趋势。三阶多项式模型对总体生产轨迹拟合良好（R² = 0.98），捕捉了研究期间产量的非线性变化。这些发现凸显了以生产力为导向的可持续集约化在加强安息香树脂农林复合系统方面的潜力，且主要不依赖耕地扩张。与生产力相关的潜在农林复合系统属性包括生物物理方面，如冠层和根系结构、有机凋落物循环和微气候调节，以及土地管理方面，如混作、修剪和间伐及土壤维护。然而，这些属性并未在本研究中直接测量。该研究受限于对二手数据的依赖，缺乏原始田间观测或直接树脂产量测量。未来研究应整合一手数据收集以验证这些模式，增强经验稳健性，并阐明安息香生产力、恢复力和可持续性结果的潜在机制。

安息香树脂（Styrax resin）农林复合系统在北苏门答腊具有重要的经济与生计价值，但长期定量研究不足。已有研究多关注其生态特征、栽培收获实践、树脂加工及社会经济意义，缺乏对产量长期变化和区域生产绩效差异的系统评估。仅凭产量增加并不能判断系统是否可持续或具有韧性，因为增长可能来自耕地扩张、单产提高、生产结构或经营方式变化。因此，研究人员利用2010–2021年北苏门答腊六个安息香生产县的县级数据，从生产趋势、生产力、增长率、波动性和区域差异五个维度开展纵向定量评估，以区分产量增长究竟受耕地扩张驱动还是受生产效率提升驱动，并为该地区传统树脂农林复合系统的可持续发展提供证据。

研究人员发现，安息香产量在该时期内显著增加，而种植面积总体下降，说明产量增长主要来自单产提高而非面积扩张。区域比较显示，Toba县和Pakpak Bharat县增长相对稳定高效，Central Tapanuli县和Humbang Hasundutan县波动较大，North Tapanuli县和Dairi县增长有限或停滞。三阶多项式模型（R²=0.98）表明总体产量轨迹呈非线性变化。这些结果突出表明，以生产力为导向的可持续集约化可以在不主要依赖耕地扩张的情况下加强安息香树脂农林复合系统。这项研究的重要意义在于提供了县级层面的长期定量证据，将产量增长与生产稳定性结合，构建了区域基准框架，为后续实地研究和差异化政策制定奠定基础。

研究数据来自印度尼西亚中央统计局（BPS）发布的六个安息香生产县（即Toba县、Pakpak Bharat县、Central Tapanuli县、Humbang Hasundutan县、North Tapanuli县和Dairi县）2010–2021年度产量与种植面积官方统计。研究人员以单位面积产量作为生产力指标；使用变异系数（Coefficient of Variation, CV）衡量生产波动性；采用皮尔逊相关系数（Pearson Correlation Coefficient, r）分析面积与产量的关联强度；利用复合年增长率（Compound Annual Growth Rate, CAGR）量化年平均增长；并拟合二阶和三阶多项式回归模型捕捉非线性趋势，以决定系数（R²）评估拟合优度。同时结合P值分析和95%置信区间（Confidence Interval, CI）检验统计显著性。这些方法共同构成对产量动态、稳定性及区域差异的定量分析框架。

3.1 产量与种植面积趋势：时间序列分析显示，种植面积在2010年最高后持续下降，尤其2012年明显减少；但总产量从2013年起加速上升，2021年达最高值。面积与产量反向变动，表明单位面积产出效率提高。

3.2 相关性分析：土地面积与安息香产量（2010–2021）：皮尔逊相关显示所有县面积与产量均呈正相关且显著（p<0.01），但强度不同。Central Tapanuli相关性最强，Toba和Pakpak Bharat中等，North Tapanuli和Dairi较弱；敏感性分析提示种植技术、气候和采脂强度等非面积因素同样影响产量。

3.3 县域间的区域差异：比较分析显示六个县的产量轨迹差异明显。Central Tapanuli面积较大但产量轨迹独特；North Tapanuli变化有限；Dairi产量低且动态弱；Toba、Humbang Hasundutan和Pakpak Bharat则表现出显著产量变化。

3.4 按县划分的安息香产量时间模式：各县时间模式各异。Central Tapanuli总体上升但年际变化大；North Tapanuli缓慢增长但幅度小；Dairi稳定；Humbang Hasundutan上升伴随明显年际波动；Toba和Pakpak Bharat波动显著。

3.5 种植面积与产量之间的关系：相关分析表明面积与产量的关联强度在不同县不一致，某些县弱相关说明不能仅用面积解释产量变化。

3.6 产量趋势的多项式建模：三阶多项式回归对总体产量拟合优度极高（R²=0.98），捕捉了非线性变化；但该模型仅描述时间模式，不识别具体驱动因素。

3.7 波动性与CAGR矩阵分析：将CV与CAGR结合，可将六县分为三类：高增长高波动（Central Tapanuli、Humbang Hasundutan）、高增长中等波动（Toba、Pakpak Bharat）、低增长低波动（North Tapanuli、Dairi）。

3.8 基准分析：基于增长、波动、趋势及面积-产量关系，进一步确认上述三类生产表现形态，为区域比较和后续优先序提供参照。

3.9 农林复合可持续性：产量增加伴随面积下降，表明可从土地效率提升实现增长，符合可持续集约化方向；但研究未直接测量土壤质量、生物多样性、碳汇等可持续性指标，需谨慎解读。

讨论部分指出，产量与面积背离表明生产增长可能来自土地利用效率改进而非耕地扩张。区域异质性说明农林复合发展受当地生态、社会经济、制度和市场条件影响，统一策略不适用。Central Tapanuli和Humbang Hasundutan的高增长伴随高波动，需要重视稳定性；North Tapanuli和Dairi的低增长低波动可能反映发展受限，但无法从数据直接归因。波动性-CAGR矩阵将区域分为“复兴区”“稳定区”和“复制区”，Toba和Pakpak Bharat可作为相对稳定增长参考。三阶多项式结果强调非线性轨迹，但并非因果或预测模型。整体而言，将增长与稳定性整合的定量基准有助于筛选和优先化区域，并支撑因地制宜的农林复合发展策略。

本研究表明，2010–2021年间，尽管北苏门答腊的安息香树脂种植面积总体下降，但其产量有所增加。种植面积与产量之间的背离表明，产量变化与土地面积变化并不成比例，说明单位种植面积产量发挥了越来越重要的作用。然而，现有县级数据无法将这一模式直接归因于特定的管理、技术、社会经济或环境因素。分析还揭示了安息香生产动态中显著的区域异质性。Toba县和Pakpak Bharat县表现为相对较高的产量增长并伴有中等波动，而Central Tapanuli县和Humbang Hasundutan县则表现为高增长伴随较大的生产波动性。North Tapanuli县和Dairi县表现为相对较低的增长和较低的波动性。这些截然不同的轨迹表明，安息香农林复合的发展不能仅用产量规模来充分代表，区域评估应同时考虑增长和稳定性。本研究建立的波动性-CAGR基准框架提供了一种定量方法，可根据增长和稳定性比较区域生产轨迹，可作为筛选和基准化工具，用于识别生产特征不同的县，并为进一步调查和针对具体情境的农林复合规划提供信息。三阶多项式分析还表明，安息香总量生产呈现非线性时间轨迹，但该模型应被解释为对时间模式的描述，而非因果模型或独立验证的预测模型。未来研究应以县域和农户层面的调查补充当前的定量分析，探讨产量和波动性差异背后的机制，重点关注社区管理、本土与传统文化生态知识、采收和采后实践、市场条件、制度支持、环境变异以及农户决策。将这些社会经济、文化和生态维度与纵向生产数据相结合，将为理解北苏门答腊安息香树脂农林复合系统的韧性和长期可持续性提供更全面的基础。