罗马尼亚区域紫花苜蓿生产系统的气候—产量关系与气候韧性

《Agriculture》：Climate–Yield Relationships and Climate Resilience of Regional Alfalfa Production Systems in Romania

　　

气候变率可通过温度、降水及胁迫性天气事件频率的变化影响饲草生产。研究人员考察了2006—2024年罗马尼亚三个发展区域——南蒙特尼亚（South-Muntenia）、东南（South-East）和东北（North-East）——的气候变率与紫花苜蓿（Medicago sativa L.）生产力。气候变量包括平均气温、最高气温、最低气温、年降水量和24 h最大降水量。采用描述性统计、基于秩的Kendall趋势分析与Sen斜率、Pearson相关分析以及含诊断检验的多元线性回归，考察了气候变量与紫花苜蓿产量之间的区域关系。研究人员还构建了一个探索性的紫花苜蓿气候韧性指数（Alfalfa Climate Resilience Index, ACRI），用于比较区域紫花苜蓿生产系统的相对气候韧性特征。该指数将两个与绩效相关的组分，即归一化平均产量和产量稳定性，与两个气候背景组分，即归一化降水量和反向归一化平均气温相结合，并评估了其对替代组分权重的敏感性。所有区域的平均年气温均显著升高，Sen斜率介于+0.071至+0.100 °C year^?1之间，而年降水量未呈现显著单调趋势。所有区域的产量斜率均为负值，但仅东北区域表现出统计学显著下降（τ = ?0.794，p < 0.001）。三个区域的平均年气温均与产量呈负相关（r = ?0.617至?0.749），而降水指标与产量均无显著双变量关系。区域多元回归模型解释了年产量观测变异的64.3%–85.1%（R² = 0.643–0.851），且气候预测因子的相对重要性和统计显著性在各区域间存在差异。ACRI将东北区域排在首位（0.760），南蒙特尼亚第二（0.469），东南第三（0.295），且这一排序在所检验的替代权重情景下保持不变。这些得分代表研究数据集中三个区域生产系统之间的相对差异，不应被解释为紫花苜蓿内在韧性的绝对度量。研究结果表明，气候—产量关系因区域而异，与本研究中所考虑的降水指标相比，温度与产量的关联更为一致。因此，ACRI被提出作为一个探索性的、依赖数据集的框架，用于数据集内区域紫花苜蓿生产系统气候韧性的比较；在更广泛应用之前，需要使用更广泛的参考数据集、统一的归一化标准以及额外的管理和环境变量进行验证。

研究背景与问题：气候变化通过升温、降水格局改变和极端天气增加对农业生产构成持续压力。紫花苜蓿（Medicago sativa L.）是重要的多年生豆科饲草，具有生物固氮、深根系统和刈后再生等特性，但仍会受到干旱和高温胁迫影响。已有研究多集中于单一胁迫因子、生理响应或具体栽培措施，将产量表现、产量稳定性与气候暴露相结合的区域尺度证据在罗马尼亚仍较缺乏。为此，研究人员利用2006—2024年的区域数据，在罗马尼亚三个发展区域评估气候—产量关系并比较相对气候韧性，论文发表于《Agriculture》。

研究概述：研究人员基于罗马尼亚国家统计局（INSSE）的县域紫花苜蓿面积和产量数据，以及罗马尼亚国家气象局（ANM）五个气象站点的气候观测数据，计算了区域产量和气候指标。通过趋势分析、相关分析、多元回归和ACRI指数，发现平均年气温在三个区域均显著上升，年降水量无显著单调趋势，东北区域紫花苜蓿产量呈显著下降；温度与产量的关联比降水指标更一致；ACRI排序为东北最高、南蒙特尼亚居中、东南最低，且该排序对权重调整稳健。该研究为区域尺度上理解苜蓿生产对气候变率的响应提供了定量框架，并强调了区域特异性适应策略的重要性。

主要技术方法：研究数据来源为INSSE县域紫花苜蓿面积和产量记录及ANM五个气象站点气候观测；以区域为单元按年汇总。方法包括描述性统计、Shapiro–Wilk正态性检验、Kendall秩趋势检验与Sen斜率、Pearson相关、多元线性回归（含残差正态性、Breusch–Pagan异方差、Durbin–Watson自相关和方差膨胀因子诊断），并构建由归一化平均产量、产量稳定性、归一化降水量和反向归一化平均气温组成的ACRI，通过改变组分权重进行敏感性分析。

研究结果：
4.1 气候变量与紫花苜蓿产量的时间趋势：Kendall分析显示，三个区域的平均年气温均显著升高，Sen斜率分别为南蒙特尼亚+0.093、东南+0.071、东北+0.100 °C year^?1；年降水量Sen斜率为负但不显著；紫花苜蓿产量Sen斜率均为负，仅东北区域显著下降（τ = ?0.794，p < 0.001，?0.571 t ha^?1 year^?1）。

4.2 气候变量与紫花苜蓿产量的关系：
4.2.1 描述性统计：南蒙特尼亚平均产量最低（13.23 ± 2.77 t ha^?1），东北最高（14.85 ± 3.23 t ha^?1）；东南平均年气温最高（12.41 ± 0.74 °C），南蒙特尼亚年降水量最高（589.69 ± 136.14 mm），东南最低（501.08 ± 118.23 mm）。

4.2.2 Pearson相关分析：平均年气温与产量在三个区域均显著负相关（r = ?0.617至?0.749）；最高气温在南蒙特尼亚和东南区域与产量显著负相关，最低气温仅在东北区域显著；年降水量和24 h最大降水量与产量均无显著相关。

4.3 多元线性回归分析：三个区域模型分别解释产量变异的64.3%、70.5%和85.1%；南蒙特尼亚仅最高气温为显著预测因子，东南仅平均气温显著，东北区域平均气温、最低气温、年降水量和24 h最大降水量均显著。

4.4 回归诊断：三个模型残差正态性和异方差性检验均不显著；东南区域年降水量和24 h最大降水量的方差膨胀因子（VIF）分别达8.05和7.13，提示存在较明显的多重共线性，该区域个别回归系数需谨慎解读。

4.5 苜蓿产量表现与稳定性：东北平均产量最高；东南变异系数（CV）最低（0.200），其次为南蒙特尼亚（0.209）和东北（0.218），三个区域年际产量变异程度相对接近。

4.6 紫花苜蓿气候韧性指数（ACRI）：等权条件下，东北ACRI最高（0.760），南蒙特尼亚第二（0.469），东南第三（0.295）。东北的归一化平均产量和反向归一化平均气温组分最高，南蒙特尼亚的归一化降水量最高，东南的归一化降水量和反向归一化平均气温均为最低；该排序不代表苜蓿内在生理韧性。

4.7 ACRI敏感性分析与区域比较：在五个权重情景下，区域排序均保持东北 > 南蒙特尼亚 > 东南，表明数据集内区域排序对权重调整具有内部稳健性，但不构成ACRI的外部验证。

讨论总结：研究人员指出，三个区域的气温显著升高与已有欧洲变暖证据一致；年降水量主要表现为年际波动，而非单调下降。温度变量比降水指标更一致地与苜蓿产量相关联，但相关不等于因果，实际产量还受蒸散、降水时程、土壤水分、品种和管理等未纳入因素的影响。ACRI应被理解为结合产量表现和气候背景的相对比较工具，而非绝对或生理意义上的韧性度量。东南区域在参考数据集内ACRI最低，提示需要关注水分利用效率和耐热耐旱品种；南蒙特尼亚需兼顾高温风险与水分管理；东北虽ACRI最高，但产量显著下降，仍需持续监测和预防性适应。研究局限包括每年仅19个观测点、产量统计口径可能为鲜重或干物质未明确、气象站点代表性有限、以及未包含土壤、品种、灌溉和管理变量等。未来应用ACRI需更广泛的数据集和统一的标准化方案。

研究结论（翻译）：在三个罗马尼亚发展区域