Page 63 - 2022年第53卷第4期

的等效水高变化 ［44］ ，数据序列最早始于 2002 年 4 月。本文使用 CSR 和 JPL 2 家研究机构公布的 2003—
2016 年 RL06 Mascon 数据产品，为了减少由不同 GRACE 后处理策略带来的不确定性，采用这 2 种数据
计算结果的算数平均值作为研究区域的陆地水储量距平 ［45］ 。由于 GRACE 卫星在其运行过程中因自身轨
道调整等因素导致部分数据缺失，采用 K 近邻算法（K Neareast Neighbor，KNN）插补缺失值 ［46］ 。
3.1.3 陆面水文模型数据 本文采用 ERA-Interim ［47］ 、GLDAS2 ［48］ 的天然蒸散发、天然径流、降水数
据。ERA-Interim 数据为逐日时间步长，空间分辨率为 0.25°。GLDAS2 采用目前发布的最新版本 GL⁃
DAS2.1，包含 Noah、VIC、CLSM 模式，三个模式由于缺少人类活动的参数仅模拟自然条件下的蒸散
发 ［14］ 。其中，Noah 数据空间分辨率为 0.25°，VIC 和 CLSM 数据空间分辨率为 1.0°，为保持一致性，用
双线性插值法将 VIC 和 CLSM 数据从 1.0°插值到 0.25°。
3.2 水文模型适用性分析

3.2.1 降水 图 3 为 2003—2016 年 4 个陆面水文模型模拟的降水、径流与实测降水、R nat 的年内变化
特征。受东南季风和西南季风的影响，珠江流域降水具有明显的季节性，4 种模拟数据均呈现出夏季
高、冬季低的特点，其中，CLSM、Noah 的降水强迫数据是相同的，CLSM、VIC 的季节变化特征与
CMFD 降水最为接近。在珠江流域，4 个模型模拟的降水与 CMFD 降水的 r 均在 0.9 以上，Noah 的精度
最高，与 CMFD 的 r 为 0.98，ERA-Interim 与 CMFD 的 r 最低，为 0.91。从西江、北江、东江三个子流域
来看，Noah、CLSM、VIC的 RMSE（17.63 ~ 31.30 mm/月）小于 ERA-Interim（47.05 ~ 65.37 mm/月）。在珠
江流域多年平均降水量上，CMFD 为 1418.38 mm/a，Noah 为 1439.82 mm/a，CLSM 和 VIC 为 1440.52
mm/a，ERA-Interim 为 1707.74 mm/a，ERA-Interim 高估了实际降水量。降水作为水文模型最重要的

400 400
300 300

200 200
100 100

0 0
1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12

400 400

300 300
200 200

100 100
0 0
1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12

400 400

300 300
200 200

100 100
0 0
1 2 3 4 5 6 7 8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12

图 3 四种陆面水文模型模拟的多年平均月降水与月径流
— 437 —