ERA5-Land雪相关变量的区域频率分析、偏差校正和降尺度

《Weather and Climate Extremes》：Regional frequency analysis, bias correction and downscaling of ERA5-Land snow related variables

　　

本文提出了对ERA5-Land雪深（snow depth, ds）和三日雪深增量（increase of snow depth over three consecutive days, DH3gg）进行区域频率分析（regional frequency analysis）、偏差校正（bias correction）和降尺度（downscaling），旨在量化无测站地点的设计荷载（design loads）并准确评估雪崩危险性（avalanche hazard）。局部统计分析通常通过提取雪深时间序列的极值进行，但这些序列往往不完整，尤其是在极端降雪事件期间。此外，大多数情况下，测站并不位于雪崩通常发生的海拔高度。相比之下，ERA5-Land提供了更广泛的空间和时间覆盖，尽管其受到偏差（biases）影响，导致极值被低估或高估。为解决这些局限性，研究人员对观测和ERA5-Land雪相关变量均提出了区域频率分析。采用的方法是指数值区域法（index value regional method），该方法根据常见的降雪和雪崩发生模式识别出同质区域（homogeneous areas）。然后，基于观测和ERA5-Land区域化变量的累积分布函数（cumulative distribution functions）之间的比较，应用了偏差校正和降尺度方法，以意大利中部地区作为案例研究。从该比较中得出的校正值（calibration values）随后用于调整ERA5-Land分位数，从而提供了一种在局部尺度外推设计荷载的简易程序。

**论文解读：ERA5-Land雪相关变量的区域频率分析、偏差校正与降尺度**

**研究背景与问题**
雪崩危险性评估在工程和土地规划中至关重要，尤其在山区城市化扩张和旅游发展的背景下。然而，雪深观测数据通常存在空间覆盖不足、时间序列不完整等问题，尤其在极端降雪事件期间，且测站多位于低于雪崩释放区的海拔。ERA5-Land全球再分析数据集提供自1950年以来的连续时间序列，具有0.1°×0.1°的空间分辨率，能覆盖无观测区域，并在雪相关变量模拟上优于其他网格产品。但ERA5-Land受偏差影响，可能高估或低估极端值，原因包括空间分辨率有限导致地形平滑、陆面模型参数化简化以及观测约束不足。为克服这些限制，研究人员提出对观测和ERA5-Land雪深及三日雪深增量进行区域频率分析，并应用偏差校正与降尺度，以提供无测站区域的设计荷载。该研究发表在《Weather and Climate Extremes》。

**主要技术方法**
研究人员采用指数值区域法进行区域频率分析，基于L-moment计算异质性度量（heterogeneity measure H）和不一致性度量（discordancy measure D_i），并利用广义极值分布（Generalized Extreme Value, GEV）拟合区域样本。偏差校正通过分位数映射（quantile mapping）实现，使用简单缩放因子C（校正系数）调整ERA5-Land无量纲分位数。降尺度通过建立指数值与海拔的线性回归关系，将区域分位数转换为局部设计荷载。案例研究中，观测数据来源于MeteoMont（意大利中部亚平宁山脉的39个手动测站），ERA5-Land节点经筛选后保留51个。

**研究结果**
- **区域识别与异质性检验**：根据雪崩发生日期相似性，将研究区划分为6个同质区域（A1至A6）。对每个区域，MeteoMont和ERA5-Land的ds和DH3gg样本均通过异质性检验（H<0），表明区域同质。
- **偏差分析**：通过比较MeteoMont与ERA5-Land的无量纲累积分布函数（CDF），计算偏差度量B(TR)。对ds，区域A3、A4、A5的偏差较大（B(300)=1.20、0.63、1.06）；对DH3gg，区域A2和A5偏差较大（B(300)=-0.78、0.55）。偏差校正后，所有区域的B(TR)显著降低，且RMSE、MAE下降，R²上升。
- **局部设计荷载外推**：利用校正后的ERA5-Land分位数和MeteoMont指数值-海拔回归方程，计算各区域在30、100、300年重现期下的局部ds和DH3gg设计荷载。以Rigopiano（1930 m a.s.l.）为例，300年重现期ds为3.30 m，DH3gg为1.75 m，释放深度约2 m，与文献一致。
- **空间分布映射**：生成300年重现期下ds和DH3gg的局部设计荷载空间分布图，无需空间插值，可直接用于指定位置。

**总结与讨论**
研究人员指出，ERA5-Land能够捕捉观测数据CDF的形状，部分区域无需偏差校正即可直接使用。区域频率分析克服了局部统计中样本量不足的问题，尤其适用于高重现期（如300年）。偏差校正通过简单缩放有效减小了ERA5-Land与观测之间的差异，且校正效果随海拔升高而增强。降尺度方法利用观测指数值-海拔关系，避免了ERA5-Land网格平均海拔带来的偏差。该程序为专业从业者和公共管理机构提供了快速估算雪崩设计荷载的工具，简化了从区域到局部的外推过程，并为雪崩动力学模拟和防御结构设计提供了现实边界条件。研究结论表明，ERA5-Land结合区域频率分析、偏差校正和降尺度，是评估无测站区域雪崩风险的有效方法。