《Agriculture》：Rain-Fed and Irrigated Vineyard Yields in Spain: Regional Differences and Associations with Annual Hydrological Conditions
编辑推荐：
区域产量比较可以补充葡萄园试验，但无法识别灌溉的因果效应。研究人员将西班牙2019–2023年期间14个自治区70个自治区-年份配对组的农业与水文统计数据进行关联分析。采用等权重配对差异法，以及具有制度特异性自治区和年份效应的联合回归分析，并使用自治区层面的野
区域产量比较可以补充葡萄园试验，但无法识别灌溉的因果效应。研究人员将西班牙2019–2023年期间14个自治区70个自治区-年份配对组的农业与水文统计数据进行关联分析。采用等权重配对差异法，以及具有制度特异性自治区和年份效应的联合回归分析，并使用自治区层面的野生聚类自助法检验。置信区间采用近似CR1-t推断。平均灌溉产量超过雨养产量1677 kg/ha（95% CI 785–2568；自助法p = 0.0013）。年降水量/参考蒸散量（P/ETo）每增加0.1个单位，雨养和灌溉斜率估计值分别为307和?73 kg/ha。两者直接差异为?380 kg/ha（95% CI ?1223至462；自助法p = 0.6306），排除加利西亚2023年后变为104（自助法p = 0.3782）。正的均值产量差异在缩尾处理和逐一排除自治区分析中仍然存在。不同水文斜率的证据仍然不精确且对影响性观测敏感。年度区域P/ETo不能衡量葡萄藤水分胁迫，且缺乏灌溉量、管理和质量数据限制了解释。研究结果确立了描述性的区域产量差异，但未证明生理缓冲效应、灌溉因果效益或更高的产量稳定性。
一、研究背景与问题提出
葡萄栽培是地中海地区最具经济、社会和文化重要性的农业活动之一。西班牙在全球葡萄种植领域占据突出地位：2024年葡萄园面积约93万公顷，居全球首位，约占全球葡萄园面积的13%。然而，地中海葡萄栽培日益面临环境脆弱性挑战——气候变化通过升温、降水格局改变、大气蒸发需求增加、干旱和极端高温等途径，显著影响葡萄发育、产量、浆果成分和葡萄酒特性。水资源可用性因此成为地中海葡萄栽培未来可持续性的核心决定因素。
现有研究存在明显的证据缺口：大量实验研究聚焦于亏缺灌溉、调节性亏缺灌溉、局部根系干燥等策略，但实验通常涉及单个葡萄园、品种或相对有限的地理区域，其结果未必能在更广泛的区域尺度上、在异质性生产系统中得到验证。西班牙提供了独特的研究场景：伊比利亚半岛涵盖从湿润西北部到干旱地中海和内陆地区的对比性气候条件，且各自治区同时存在雨养和灌溉两种生产系统。因此，区域层面的核心问题不仅是灌溉葡萄园是否产量更高，更在于两种系统的产量与自然水资源可用性之间的关联是否存在差异。然而，区域比较需要谨慎：同一自治区内的灌溉与雨养葡萄园可能在品种、土壤和管理上存在差异，年度区域指标无法解析葡萄藤水分状况或生长季胁迫的时间特征，这些生态学关联不能确立灌溉本身产生生理缓冲效应的因果结论。
二、研究设计与核心发现
该研究是一项生态学研究，利用西班牙农业、渔业和食品部的作物面积与产量调查（ESYRCE）数据，以及西班牙国家气象局（AEMET）的水量平衡统计数据，对2019–2023年间14个自治区70个自治区-年份配对组进行分析。研究目标包括：量化平均配对产量差异；评估各制度下年度水资源可用性与产量的自治区内关联；直接比较这些关联。主要发现为：平均灌溉产量超过雨养产量1676.6 kg/ha（95% CI 784.9至2568.2；自助法p = 0.0013），且这一正向差异在排除加利西亚2023年、缩尾处理和逐一排除自治区分析中均持续存在。在联合制度特异性模型中，雨养产量与年P/ETo的关联为306.8 kg/ha每0.1单位（95% CI 21.5至592.1；自助法p = 0.0312），而灌溉关联为?73.4 kg/ha每0.1单位（95% CI ?698.6至551.8；自助法p = 0.8074），两者直接对比为?380.2 kg/ha每0.1单位（95% CI ?1222.6至462.2；自助法p = 0.6306），未提供不同斜率的清晰证据。排除加利西亚2023年后，雨养估计降至188.2 kg/ha，灌溉估计升至292.2 kg/ha，直接对比变为104.0 kg/ha（自助法p = 0.3782），斜率对比的符号发生改变。
三、主要技术方法
研究人员采用的主要方法包括：等权重配对差异分析，通过截距回归估计平均产量差异并按自治区聚类；联合固定效应模型，允许自治区效应和年份效应按制度（雨养/灌溉）不同，以直接比较两种制度的斜率；少聚类推断采用空假设强制的野生聚类自助法t检验（WCR11），使用Rademacher权重在自治区层面进行，对14个自治区枚举全部16384种符号模式；置信区间采用近似CR1-t推断；敏感性分析包括排除加利西亚2023年、逐一排除每个自治区、对核心变量进行1st/99th百分位缩尾处理，以及用降水（P）、参考蒸散量（ETo）、P减ETo、PADMAX、AD25和AD75等替代指标分别建模。数据来源为ESYRCE区域汇总工作簿和AEMET水量平衡统计。
四、研究结果分节解读
3.1 面板特征与配对差异： 合格样本包含70对配对，14个自治区各5年。雨养平均产量为5045.5 kg/ha，灌溉平均产量为6722.0 kg/ha。灌溉产量在59对（84.3%）中超过雨养产量。平均配对差异为1676.6 kg/ha（95% CI 784.9至2568.2；自助法p = 0.0013），排除加利西亚2023年后为1826.8 kg/ha（95% CI 996.4至2657.2；自助法p = 0.0012）。年度非加权平均差异在整个观测期内均为正。
3.2 水文斜率与直接比较： 联合制度特异性模型估计雨养关联为306.8 kg/ha每0.1单位P/ETo（95% CI 21.5至592.1；自助法p = 0.0312），灌溉关联为?73.4 kg/ha每0.1单位（95% CI ?698.6至551.8；自助法p = 0.8074），灌溉区间宽泛且包含正负关联。直接灌溉减雨养斜率对比为?380.2 kg/ha每0.1单位（95% CI ?1222.6至462.2；自助法p = 0.6306），未提供不同斜率的清晰证据。排除加利西亚2023年后，雨养估计变为188.2 kg/ha（95% CI ?7.3至383.6；自助法p = 0.0548），灌溉估计变为292.2 kg/ha（95% CI ?31.0至615.3；自助法p = 0.0868），直接对比变为104.0 kg/ha（95% CI ?179.0至387.0；自助法p = 0.3782），对比不精确且改变符号。
3.3 敏感性与诊断发现： 1st/99th百分位缩尾处理后，平均差异为1691.7 kg/ha，斜率对比为?299.3 kg/ha每0.1 P/ETo（95% CI ?1020.2至421.7；自助法p = 0.5142）。逐一排除自治区分析中，平均差异范围从1426.9到1866.2 kg/ha，斜率对比范围从?545.7到58.9 kg/ha每0.1 P/ETo，符号取决于保留的自治区。排除全部加利西亚后，对比为58.9 kg/ha每0.1 P/ETo。替代指标分析中，雨养降水斜率为310.5 kg/ha每100 mm（95% CI ?32.8至653.8；自助法p = 0.0824），气候平衡斜率为249.7 kg/ha每100 mm（95% CI ?55.8至555.3；自助法p = 0.1006），未建立制度间一致的差异。
3.4 与原模型比较： 重建的原始共同效应交互为?191.9 kg/ha每0.1 P/ETo（WCR11自助法p = 0.0430），表明原始强显著性在改变模型设定前就依赖于推断程序。允许自治区和年份效应按制度不同后，交互估计变为?380.2 kg/ha（WCR11 p = 0.6306），这是修订结论的两个不同原因。
五、讨论与结论
区域层面的正向平均产量差异约为1.7 t/ha，在各项敏感性分析中持续存在，但该差异不估计灌溉相同葡萄园的产量增益，因为自治区和年份内配对无法消除两种生产系统在品种、立地和管理上的差异。雨养产量与年P/ETo的正向关联在主要自治区内分析中约为307 kg/ha每0.1单位，排除加利西亚2023年后降至约188 kg/ha，其强度和推断分类对样本和程序敏感。灌溉斜率估计不精确，直接灌溉减雨养斜率对比在排除加利西亚后改变符号，数据不能确立不同水文斜率或灌溉缓冲区域水分胁迫的结论。
研究结论为：在合格的西班牙自治区-年份面板中，平均灌溉产量超过雨养产量约1.7 t/ha，正向非加权差异在指定敏感性分析中持续存在。年P/ETo在主要自治区内分析与雨养产量正相关；排除加利西亚2023年后，估计减弱，近似区间跨越零，自助法p值接近阈值。雨养与灌溉斜率的直接比较不精确且在敏感性分析中改变符号，数据不能确立灌溉下较低水文敏感性的结论。这些区域关联不识别因果灌溉响应或更高产量稳定性。亏缺灌溉或其他节水策略的决策需要关于水资源利用、管理和葡萄品质的位点特异性证据，这些考虑也与企业社会责任相关，特别是通过葡萄栽培价值链中负责任的水资源管理。
