编辑推荐：
这篇研究聚焦于华北湿润与半湿润暖温带这一气候敏感过渡区，通过融合标准化潜在蒸散指数（SPEI）、标准化土壤湿度指数（SSMI）和CASA光能利用率模型（2001-2020年），系统评估了干旱时空特征及其对植被净初级生产力（NPP）的动态影响。研究揭示了该区域干旱呈现轻微缓和但短期高频高强度事件频发的复杂格局，并明确了降水是驱动植被NPP变化的主导气候因子。研究阐明了不同植被类型（以草地最敏感、林地具抗性且响应延迟）对干旱胁迫的差异化响应机制，为理解气候变化下生态脆弱区碳循环响应及制定适应性管理策略提供了关键科学依据。
1. 引言
植被是陆地生态系统和碳循环的关键组分，对全球环境变化高度敏感。净初级生产力（NPP）作为衡量生态系统生产能力和碳吸收能力的关键指标，其变化深刻反映了环境胁迫的影响。干旱作为一种长周期、大范围的气象灾害，通过引发水资源短缺、作物减产等连锁反应，威胁着生态和社会经济的稳定。在华北地区，特别是其湿润与半湿润暖温带区域，作为典型的干湿气候过渡带，生态系统对气候变化的响应尤为敏感。因此，评估该区域的干旱特征及其对NPP的影响，对于研究区域气候变化、提升生态系统适应性、促进可持续发展具有至关重要的意义。本研究旨在通过多时间尺度的干旱指数分析与NPP动态模拟，量化干旱胁迫对不同植被生产力的影响，为区域生态韧性评估和“双碳”目标实现提供定量支持。
2. 材料与方法
2.1. 研究区域
研究区域为华北的湿润与半湿润暖温带区域。其地理范围南至秦岭-淮河一线，西至乌鞘岭-祁连山，涵盖了北京、天津、山东全境，以及河北、山西、陕西、河南的大部分地区，还包括甘肃、宁夏、安徽、江苏的部分区域。该区域气候表现为年均降水量400–800毫米，夏季降雨集中，西北部以山地为主（草地或林地占优），东南部则为低地农业主产区。
2.2. 数据来源与处理
研究利用了2001年至2020年的多源数据。气象数据（月均温、月降水量）及干旱指数（Aridity Index, AI）数据集来源于青藏高原科学数据中心。总太阳辐射数据来源于欧洲中期天气预报中心（ECMWF）的ERA5-Land数据集。土壤湿度数据使用了全球每日地表土壤湿度数据集。土地利用数据基于MODIS MCD12Q1产品获取。标准化潜在蒸散指数（SPEI）数据来源于一个覆盖中国大陆的1公里分辨率数据集。此外，还使用了MOD17A3HGFv061产品作为NPP的基准数据，以及MOD13A1v061产品用于获取归一化植被指数（NDVI）。
2.3. 研究方法
- •
CASA模型模拟NPP：采用CASA（Carnegie-Ames-Stanford Approach）光能利用率模型来估算研究区域的植被NPP。模型的核心公式为NPP(x, t) = APAR(x, t) × ε(x, t)，其中APAR为植被吸收的光合有效辐射，ε为实际光能利用率，由温度和水分胁迫系数以及最大光能利用率决定。
- •
干湿变化与干旱识别：利用干旱频率指数（Drought Frequency Index, DFI）来度量研究期内干旱发生的频率。当SPEI ≤ ?1时，即认为发生了“干旱”。同时，使用标准化土壤湿度指数（SSMI）来量化土壤水分亏缺状况，作为对SPEI所表征气象干旱的补充验证。
- •
统计分析：采用Theil-Sen中值趋势分析与Mann-Kendall（M-K）检验来评估温度、降水、SPEI、SSMI和NPP的时间变化趋势。通过皮尔逊相关分析和偏相关分析来量化NPP与气候因子（温度、降水、太阳辐射）及干旱指数（SPEI， SSMI）之间的统计关联。
3. 结果与分析
3.1. 干湿变化与干旱识别
2001-2020年间，研究区气温和降水量均呈上升趋势，分别以0.02°C/年和2.65毫米/年的速率增加。干旱指数分析显示，研究区域整体干旱有轻微缓和趋势。长期尺度（SPEI-12）上干旱强度和范围相对较小，而短期尺度（SPEI-1）上则频繁发生持续时间短、强度高的干旱事件。季节性分析表明，冬季干旱范围呈增加趋势，而春、夏、秋季的干旱范围则有所下降。空间上，干旱高发区主要集中在河北南部和河南北部等地。
3.2. NPP时空变化特征
研究期内，研究区域植被NPP整体呈上升趋势，平均增长率为5.35克碳每平方米每年（5.35g C·m-2·a-1）。不同土地利用类型中，草地的NPP年平均增长率最高（8.02g C·m-2·a-1），其次为林地（4.03g C·m-2·a-1）。空间分布上，NPP呈现出从中部、中北部地区向东部和西部递减的格局，高值区集中在大行山、吕梁山和六盘山等山区，主要对应森林生态系统。
3.3. 影响NPP的因素分析
- •
气候因子的影响：降水是影响植被NPP变化的主要驱动因子，与研究区NPP的平均相关系数为0.25，且86.79%的区域呈正相关。温度对NPP的影响存在区域分异，在关中平原等地区，适宜的降水配合较高的温度促进了NPP积累。太阳辐射则表现出较弱的负相关。
- •
干旱对NPP的影响：在月尺度上，NPP与干旱指数的相关性较弱。在年尺度上，NPP与SPEI-12呈正相关的区域占80.16%，其中显著相关的区域占16.5%，高相关区主要集中在草地。在季节尺度上，夏季SPEI-3与NPP的正相关性最强（峰值0.92），而冬季则主要表现为负相关（峰值-0.91）。这表明在生长季，水分充足促进生产力；在寒冷限制期，基于蒸散-降水平衡的SPEI更能代表植被水分胁迫。
- •
不同植被类型的响应差异：草地生态系统对干旱的敏感性最高，与两个干旱指数的相关性均最强。耕地生态系统表现出中等敏感性，而林地生态系统的响应最弱，显示出一定的抗旱性和响应滞后性。在不同干旱等级下，草地NPP的下降幅度最大，林地NPP仅在重度和极端干旱条件下才表现出较明显的下降，且幅度小于草地。
4. 讨论
华北地区作为典型的气候过渡带，其生态系统和农业生产方式具有明显的过渡性特征。本研究发现，作为气候敏感过渡带的华北湿润半湿润暖温带地区，其净初级生产力对干旱事件响应迅速。研究区域呈现“暖湿化”趋势，温度与降水同步上升轻微缓解了整体干旱。长期尺度上，稳定的降水变率和水热同步使得干旱强度维持在较低水平；但短期尺度上，春季降水亏缺和季风来临前的水分短缺导致了高频高强度干旱事件的发生。
植被NPP的整体上升主要受气候驱动，降水增加增强了水分供应，适宜的温度延长了生长季并提高了光合效率。植被分异也起到重要作用，中部高NPP区以草地/林地为主，而东、西部受耕地扩张和水分限制制约，增长速度较慢。此外，人为干预（如生态恢复、灌溉）显著提升了草地的年增长率。
NPP变化受降水（主导）、温度和太阳辐射共同调控。降水与NPP呈显著正相关，尤其是在水分受限的半湿润地区。温度对NPP的影响存在地理分异，在山东等增温显著的区域，温度上升通过延长生长季等机制促进了NPP积累。太阳辐射的微弱负相关可能与辐射增强导致蒸散加剧、水分胁迫有关，但在辐射充足的背景下，植被对光变化的敏感性较低。
植被对干旱的响应机制存在差异：月尺度相关性弱，源于植被的短期缓冲能力；而季节/年尺度因累积水分整合作用，表现出显著相关性。综合使用SPEI和SSMI能更全面、准确地揭示干旱影响植被NPP的复杂机制，具有明显的协同优势。草地因其根系浅而密集，对表层水分可获得性高度敏感，因此在各类干旱事件中表现出更强的同步性。树木则对环境变化具有一定抵抗力，且响应存在一定滞后性。
5. 结论
本研究的主要结论可归纳为以下四点：
（1）华北气候敏感过渡带呈现“暖湿化”趋势，温度与降水同步上升使干旱整体略有缓和。稳定的降水变率与水热同步维持了长期较低的干旱强度，但短期尺度上频发高强度干旱。
（2）植被NPP整体呈增加趋势（55.10%的区域显著增长），空间上由中心向东、西递减。草地NPP的年均增长率最高（8.02g C·m-2a-1），其次为林地（4.03g C·m-2a-1）。
（3）降水主导植被NPP变化；温度上升促进NPP积累但存在地理分异。在华北气候敏感过渡带，充足光照导致辐射敏感性低，温度与辐射的影响弱于降水。
（4）植被响应干旱的机制各异：草地因高同步性而反应最明显，森林则表现出抗旱性与响应延迟。在时间分布上，年尺度上80.16%的区域（主要为草地）与SPEI-12指数呈正相关。季节尺度上，SPEI-3指数在夏季达到最大值0.92，冬季则为-0.91。这种相关性显著强于月尺度SPEI-1所指示的弱关联。
