华北平原夏玉米物候对齐极端降雨敏感窗口识别及可操作灌排触发指标

《Agricultural Water Management》：Identification of phenology-aligned sensitive windows to extreme-rainfall events for maize in north China Plain and actionable irrigation–drainage triggers

　　

农业水风险评估需识别能准确捕捉作物敏感期极端降雨的暴露窗口。研究人员利用华北平原2001至2015年县级数据重建物候对齐暴露窗口，定义暴露错配指数（Exposure Mismatch Index，EMI）为物候对齐窗口与静态窗口间超阈值累积降水量的方向性差值。在50 mm报告基准阈值下，静态W2（营养生长后期至生殖转换窗口）相对物候对齐窗口的漏报率与误报率分别为9.7%与10.0%，总误分类率达19.7%；在30至70 mm阈值区间总误分类率为17.3%至19.9%。分阶段模型比较显示W2在长周期灌溉覆盖层间对比最清晰：灌溉覆盖<0.90的县，调整W2响应在EMI为?120 mm至?50 mm及+40 mm至+110 mm区间呈不利关联带，而灌溉覆盖≥0.90的县该关联被削弱。全球陆面同化系统（Global Land Data Assimilation System，GLDAS）根区土壤水分数据对<0.90灌溉覆盖层的阳性带仅提供有限水文支持，未证实阴性带存在水分亏缺或直接生理介导。上述结果区分了极端降雨发生与其物候对齐产量的统计可见性，广灌溉覆盖下弱关联不应被解读为无降雨暴露或作物固有耐受。EMI可作区域暴露窗口错配筛查指标：正错配提示审查排水准备，负错配提示评估土壤水可用性，两信号均须结合预报、土壤水分、排水条件与实际灌溉信息确认。

研究背景方面，气候变化下农业韧性成为灌溉规划与粮食安全的重要目标，水文气候模式在强度、频率与时机上的改变增加了产量年际波动，但不少农业影响评估仍把气候暴露聚合在静态日历或全生长季窗口中，未追踪作物物候的年际变化。华北平原夏玉米物候存在显著时空变异，而极端降雨若发生在生殖转换期可能造成不同于营养期的损害，静态窗口会系统性错分暴露：相似雨量在营养期后果轻、在生殖期可致损。同时管理干预（如灌溉）会衰减或重塑气候—产量关联，在异质管理背景下 pooled 分析可能掩盖低灌溉覆盖区的脆弱性。因此研究人员开展此项研究，量化物候变异下的阶段相关暴露并刻画其在异质灌溉情境中的暴露—产量关联，论文发表于《Agricultural Water Management》。

关键技术方法上，研究人员构建2001至2015年华北平原县级—年面板，产量取自1 km夏玉米数据集并采用双向固定效应去县年效应得残差；物候由MODIS增强植被指数（EVI，MOD13Q1 V6.1）经插值平滑提取季节峰值P及50%振幅处的G、S锚点，将G–P与P–S各取中点划分为W1至W4四个物候对齐窗口；降雨取自CHIRPS，以各县湿日95百分位R₉₅p_i为超阈值界，分别在物候对齐与静态日历窗口累加超阈值降水；定义EMI为两者累积量之差；灌溉覆盖为历史灌溉数据集重建的耕地灌溉装备面积比，按2001至2015均值以0.90为预设描述边界分层；协变量含经纬、高程、坡度、土壤质地、有机碳、气候均值与潜在蒸散；采用广义倾向得分（generalized propensity score，GPS）加权联合广义加性模型（generalized additive model，GAM）估计灌溉 conditioned EMI—产量非线性关联；根区土壤水分取自GLDAS-2.1 Noah模型作水文兼容性差分校验。

研究结果部分，引言指出静态窗口聚合模糊了生物关键阶段风险，管理干预使气候—产量关联衰减，需物候对齐暴露与灌溉分层双视角。材料与方法小节中，研究区覆盖京津冀及山东西河南部夏玉米主产平坦温带季风区，降水集中夏季且年际变异大，灌溉覆盖空间差异显著；数据预处理完成产量、物候、降雨、环境与管理协变量及GLDAS土壤水分的同边界面积加权聚合；2.2方法论小节给出EMI构造方程（E_dyn,it=∑_d∈Wdyn,it P_itd·I(P_itd>R₉₅p_i) 与 E_sta,it 对应静态窗，EMI=E_dyn?E_sta）及GPS–GAM分层估计流程。

讨论与结论部分，研究人员发现静态W2在50 mm阈值下相对物候对齐窗口漏报9.7%、误报10.0%，总误分类19.7%，30至70 mm区间总误分类17.3%至19.9%，证明静态日历窗会系统性错分极端降雨暴露。W2在灌溉分层中对比最强：灌溉覆盖<0.90县，调整W2响应在EMI ?120至?50 mm与+40至+110 mm出现不利产量关联带；≥0.90县关联被衰减，但GLDAS根区土壤水分仅对<0.90层阳性带有限支持，未证实阴性带缺水或生理介导。由此结论为，EMI区分了极端降雨发生与其物候对齐的统计可见性，广灌溉下弱关联不等于无暴露或固有耐性；EMI可作区域筛查指标，正错配促审查排水准备、负错配促评估土壤水可用量，两信号须由预报、土壤水分、排水与实际灌溉信息确认。该研究贡献在于将预警窗误差量化为漏报误报、在县—年尺度控制环境异质性估计灌溉 conditioned EMI—产量关联、并在W2窗识别非线性关联带以支撑灌排准备的区域筛查。