【WRF模拟技术】Linux编译排错、FNL/ERA5驱动场处理、长时序模拟配置、下垫面改造与物理参数调整

专题一：数据的空间变换与插值

1. 降尺度导论

1.1 尺度不匹配的物理本质

• GCM (>100km) 与流域/城市 (<10km) 的尺度鸿沟

• 动力(RCM) vs 统计(ESD)：物理完备性与计算成本的权衡

2. 格点转站点 (Downscaling)

2.1 最近邻法 (Nearest Neighbor)

• 缺陷：山区“阶梯效应”与海拔误差

2.2 双线性插值 (Bilinear)

• 平滑效应：为何插值削弱暴雨峰值导致灾害低估

3. 站点转格点 (Interpolation)

3.1 反距离加权 (IDW)

• 距离衰减幂指数选择与“牛眼”现象

3.2 薄板样条 (TPS/ANUSPLIN)

• 引入协变量 (DEM) 模拟“气温直减率”

4. 案例实操：

4.1数据读取与可视化

读取CMIP6 NetCDF，可视化粗网格数据

4.2插值对比实验

分别用最近邻和双线性提取某点数据，对比二者在极端降水日的误差

4.3高阶空间插值

专题二：模型评估与极端指数

1. 综合统计指标

1.1 基础指标：RMSE (对异常值敏感) vs MAE, Bias vs PBIAS

1.2 泰勒图 (Taylor Diagram)：极坐标系中同时展示 R、STD、RMSE

2. 降水与风的评估

2.1 降水探测：混淆矩阵、POD (命中率)、FAR (虚警率)

2.2 风速分布：Weibull 分布拟合

3. 极端气候指数 (ETCCDI)

3.1 核心指数：Rx5day (洪涝)、R95p (强降水占比)、CDD (干旱持续日数)

3.2 趋势分析：Mann-Kendall 检验

4. 案例实操：

4.1评估工具开发

编写evaluate_model(obs, sim) 函数

4.2泰勒图绘制

评估多个CMIP6 模式性能，筛选最优模式

4.3极端指数计算

专题三：建立物理联系(回归/CCA/类比)

1. 多元线性回归

1.1 预测因子筛选：逐步回归 (Stepwise AIC/BIC) 与 多重共线性 (VIF)

1.2 方差膨胀：解决回归导致的“极值削弱”问题 (Variance Inflation)

2. 典型相关分析 (CCA)

2.1 降维与场相关：EOF/PCA 提取环流主模态

2.2 CCA 原理：寻找大尺度场与局地要素场相关性最大的线性组合

3. 相似形法 (Analog Methods)

3.1 构造相似形：在历史库中寻找与当日环流最相似的日子 (“历史重演”)

3.2 流型识别：Lamb 环流分型在降水预测中的应用

4. 案例实操

4.1回归建模

利用ERA5 位势高度/温度建立气温回归模型，并进行方差修正

4.2 EOF 与 CCA 实战

使用eofs 库提取 500hPa 高度场模态；建立大尺度场与区域多站点的 CCA 映射

4.3相似日搜索算法

专题四：随机模型(天气发生器)

1. 降水发生器 (WGEN核心)

1.1 状态模拟 (Occurrence)：一阶/二阶马尔可夫链 (Markov Chain)，计算干湿转移概率 

1.2 强度模拟 (Intensity)：Gamma 分布或混合指数分布拟合

2. 非降水变量

2.1气温/辐射：基于降水状态的条件残差模型 (Richardson type)

3. 气候变化情景生成

3.1修改统计参数(Change Factors) 驱动发生器生成未来长序列

4. 案例实操

4.1马尔可夫链手写

计算某站30 年数据的降水转移概率矩阵

4.2蒙特卡洛模拟

基于概率矩阵，随机生成100 年日降水状态序列

4.3模型验证

专题五：偏差校正技术(QM & QDM)

1. 基础校正

1.1 线性缩放 (Linear Scaling)：仅修正均值

1.2 局地强度缩放 (LOCI)：修正 GCM 的“毛毛雨”频率偏差

2. 分位数映射 (Quantile Mapping)

2.1 CDF 匹配原理：强制模拟分布与观测一致

2.2 非参数与参数法：经验分位数表 vs Gamma 分布映射

3. 趋势保持校正 (QDM)

3.1 QM 的气候缺陷：直接 QM 会抹杀未来的极端增温/增雨信号

3.2 QDM 算法：(1)提取趋势 -> (2)校正历史 -> (3)叠加趋势

4. 案例实操

4.1基础校正练习

实现Linear Scaling，观察其无法修正方差的局限性

4.2 QM 代码实现

编写经验分位数映射代码，修正GCM 数据的 PDF

4.3 QDM 实验

专题六：机器学习降尺度

1. 特征工程 (Feature Engineering)

1.1物理因子：比湿、涡度、散度、风切变

1.2时空特征：DOY (年积日)、经纬度、海拔、地形坡向

2. 表格类模型 (Station-based)

2.1算法：Random Forest / XGBoost / LightGBM

3. 案例实操

3.1数据集构建

利用xarray 对齐 ERA5 (特征) 和 站点观测 (标签)

3.2 XGBoost 建模