《Journal of Hydrology: Regional Studies》：Snow cover over the Iberian mountains in km-scale global climate simulations: Evaluation and projected changes
编辑推荐：
研究区域聚焦伊比利亚半岛四大主要山脉。此类地中海山区地形复杂、气候变率大，积雪作为关键自然资源对气候变化高度脆弱。研究聚焦于传统气候模式因分辨率较粗无法准确刻画地形过程的问题，评估全球千米尺度IFS-FESOM模型模拟伊比利亚山脉季节性积雪及预估气候变化影响的
研究区域聚焦伊比利亚半岛四大主要山脉。此类地中海山区地形复杂、气候变率大，积雪作为关键自然资源对气候变化高度脆弱。研究聚焦于传统气候模式因分辨率较粗无法准确刻画地形过程的问题，评估全球千米尺度IFS-FESOM模型模拟伊比利亚山脉季节性积雪及预估气候变化影响的能力。该方法可直接从全球模拟获取高质量区域信息，无需传统区域降尺度。研究人员将模型结果与高分辨率再分析产品、卫星反演及站点观测进行对比验证，并在SSP3-7.0排放情景下预估未来积雪变化趋势。研究获得新的区域水文洞见：尽管模型存在降雪日数正偏差，但仍能较好表征积雪气候学的关键特征，包括海拔依赖性、季节循环及平均积雪面积。预估结果显示伊比利亚全境积雪季显著缩短，低至中海拔区及最南端山脉缩减幅度最大。这种变化与降雪日数减少直接相关，主要由升温驱动，降水减少亦起到重要作用。结果表明，全球千米尺度模拟可为复杂地形区提供物理一致性的局地气候预估，是传统区域化方法的简洁稳健替代方案。
该研究发表于《Journal of Hydrology: Regional Studies》，针对伊比利亚半岛山地积雪对地中海气候区水资源供给的核心作用，回应了传统粗分辨率全球气候模式（如CMIP6模式，典型分辨率约100 km）无法解析复杂地形过程、传统降尺度方法存在经验依赖性与大尺度不一致的局限。研究依托nextGEMS计划的30年千米尺度IFS-FESOM模拟（历史时段1990–2019年，SSP3-7.0预估时段2020–2049年），验证其在伊比利亚四大山脉（中央山脉、坎塔布连山脉、内华达山脉、比利牛斯山脉）的积雪模拟能力，并预估近期气候变化影响。研究的重要意义在于证明全球千米尺度模拟可跳过区域降尺度的资源与时间成本，实现区域与全球过程的双向反馈，为山地水文气候研究提供新范式。
关键技术方法包括：采用nextGEMS计划IFS-FESOM模式（大气陆地模块分辨率约9 km，海洋模块分辨率可达5 km）的历史与SSP3-7.0模拟数据；整合多源参考数据，包括ERA5-Land再分析、IPE-CSIC雪深网格数据集、CERRA与CERRA-Land区域再分析、西班牙国家气象局（AEMET）人工站点观测及MODIS卫星归一化差异雪指数（NDSI）；数据预处理统一插值至0.125°×0.125°经纬度网格，计算积雪覆盖日数（SCD，阈值1 cm）、积雪面积（SCA）、降雪日数（SfD）、霜冻日数（FD）、雨日（WD）等指标；采用多重线性回归量化温度与降水对积雪变化的贡献，通过双侧t检验与Benjamini–Hochberg程序控制错误发现率（FDR，αFDR=0.05）评估统计显著性。
研究结果分为模型评估与未来变化两部分。4.1节模型评估包含三个子部分：4.1.1平均积雪持续时间显示，IFS模拟的积雪日数海拔依赖性与CERRA-Land的S型曲线一致，优于ERA5-Land的线性偏差，但在低至中海拔存在正偏差，内华达山脉高海拔区则略偏低；偏差溯源表明正偏差源于相对湿度的负偏差（IFS较CERRA低2–4个百分点），导致相同总降水量下降雪比例偏高。4.1.2平均季节循环显示，IFS模拟的积雪面积季节峰值与参考数据集一致，但低至中海拔区面积偏大，季节曲线更宽，春季消融期更长；MODIS卫星数据验证了季节形态，但受森林覆盖与地形阴影影响存在局地低估。4.1.3分布与年际变率显示，IFS模拟的积雪面积年际变率与CERRA-Land、ERA5-Land的四分位距一致，站点尺度对比中，开阔非森林区的模拟值与观测吻合良好，森林覆盖区与低海拔区受平滑效应影响存在偏差。4.2节未来变化包含两个子部分：4.2.1积雪持续时间变化显示，2031–2049年较1991–2009年，积雪覆盖日数绝对减少量随海拔升高而增加，但相对减少率呈海拔梯度（低海拔降幅超50%）与南北梯度（内华达山脉高海拔降幅超40%，比利牛斯山脉高海拔降幅不足20%），与已有区域研究结论一致。4.2.2积雪减少的驱动因子显示，降雪日数减少是积雪缩短的直接原因，其与积雪日数变化的皮尔逊相关系数中位数为0.81；多重线性回归表明，升温导致的霜冻日数减少是主导因子，低至中海拔区贡献更显著，高海拔及南部山脉（中央山脉、内华达山脉）的雨日减少对降雪减少的贡献增强，比利牛斯山脉南坡中海拔区雨日减少亦为主要驱动。
讨论部分指出，单模式单情景模拟的结果受内部变率与降水预估不确定性的限制，但2031–2049年的升温信号足够稳健，支撑了温度驱动积雪减少的结论。研究结论总结如下：全球千米尺度IFS-FESOM模拟虽在低至中海拔存在积雪日数正偏差（源于相对湿度负偏差导致的降雪比例偏高），但仍能表征伊比利亚山脉积雪气候学的核心特征，其模拟结果处于多源参考数据集的不确定性范围内，局地站点对比亦显示兼容性。SSP3-7.0情景下，1991–2009年至2031–2049年积雪季显著缩短，低至中海拔及南部山脉受影响最严重，低海拔区积雪日数相对降幅超50%，比利牛斯山脉高海拔区降幅约20%。积雪缩短由降雪日数减少驱动，升温导致的霜冻日数减少是主要成因，降水减少对高海拔及南部山脉的贡献不可忽视。该研究证实全球千米尺度模拟可作为传统区域降尺度的简洁稳健替代方案，为地中海山区水资源适应规划提供物理一致的局地气候预估支撑。
