观潮 · 模式预报与气象 AI 观察

华宇全球降水反演:30 分钟一次的准全球场,CSI 0.441 压过延后 3.5 个月的 IMERG FR

2026-10-23

> 卫星降水产品一直有个换不掉的取舍:实时的快但不准,准的要等几个月。作者做的 Huayu-Global 把 GOES ABI、Meteosat SEVIRI、FY-4B AGRI 三条静止卫星数据流塞进同一套检索框架,每 30 分钟出一次 0.05° 的准全球降水场。在 5,395 个独立地面站、1,847,967 条配对的三小时记录上,它的 CSI 是 0.441,比延后 3.5 个月发布的 IMERG FR(0.418)还高,RMSE 1.961 mm 也是八类产品里最低。作者把边界划得很清楚:静止卫星只覆盖到南北纬 60°,非洲、南美和洋面的验证站稀疏,暖雨、降雪、地形雨和强对流依然是难点。

缺的是哪一块

地面雨量站、天气雷达和卫星—站点融合分析都能约束降水场,但站网密度不均、雷达有盲区、地形遮挡、洋面没有覆盖、数据共享条件也不一致,真正需要快速且跨洲一致给出降水信息的时候,这些限制最要命。

卫星产品补上了覆盖,却把选择权推给了用户:GSMaP_NOW、PERSIANN-CCS 这类实时流来得快,但探测与定量技巧都低;CMORPH、IMERG 这类经过更多处理的产品精度更好,最高质量的那一档可能要等上几个月——IMERG Final Run 是观测后 3.5 个月。原文的判断很直接:决策窗口之后送达的信息,支撑不了洪水预警和其他时效敏感的行动。

往高频静止观测上走是一条出路,因为它可以连续采样云的演变。难点在于:云顶亮温并不唯一地对应地面降水,这个关系随区域、季节、地形、风暴类型和传感器几何而变;更难的是跨平台——GOES ABI、Meteosat SEVIRI、FY-4B AGRI 在光谱覆盖、原生空间分辨率和全圆盘扫描节奏上都不一样,而一个业务化的准全球产品必须把它们无缝拼起来。

方法怎么搭的

预处理把不同平台的观测转成地理参考的多光谱云特征张量,重投影到统一的 0.05° 网格。ABI、SEVIRI、AGRI 各自做通道归一化和特征标准化,以削弱跨传感器异质性同时保留各传感器自己的信息。无效像元和低质量反演区用邻近有效像元的空间插值填充,跨越整个插值窗口的持续性缺测标记为缺失,不参与训练和评估。

训练对靠的是把静止特征张量与 IMERG FR 降水场在同一时刻配对:每 30 分钟,一个 IMERG FR 场配 2 景 AGRI 或 SEVIRI、或 3 景 ABI(按各自原生获取频率)。推理时各扇区输出转到统一网格,用重叠区感知的融合方式拼接,相邻扇区在重叠带内取平均。

三个扇区模型共用 Head—地理注意力模块(GeoAB)—Tail 的架构,但 ABI、SEVIRI、AGRI 各自一套参数。每个扇区模型含 1 个 Head Block、10 个 GeoAB 模块和 1 个 Tail Block,带残差连接;GeoAB 把局地云纹理注意力与更大范围的地理上下文结合起来,用来同时抓小尺度对流特征和大尺度降水型态。训练用 AdamW,批大小 16,100 个 epoch,4 张 NVIDIA RTX 6000 Ada(每张 48 GB),PyTorch + Python 3.12。

运行时性能:模型初始化之后,单个工作站上每个时次的数据加载加推理不到 5 分钟。

结果数字

全样本对比(1,847,967 条配对的三小时站—产品记录,评估窗口为 2022 年 7—12 月与 2025 年 1—6 月,均不在 2023—2024 训练年之内):

产品 · 分辨率 / 频次 · 延迟 · CSI ↑ · ACC ↑ · R² ↑ · CC ↑ · RMSE ↓ (mm)

PERSIANN-CCS · 0.04° / 1 h · 实时 · 0.234 · 0.853 · −0.177 · 0.334 · 2.660

GSMaP_NOW · 0.1° / 1 h · 实时 · 0.327 · 0.859 · −0.374 · 0.387 · 2.874

PERSIANN · 0.25° / 1 h · 2 天 · 0.281 · 0.855 · 0.156 · 0.421 · 2.257

IMERG ER · 0.1° / 30 min · 4 小时 · 0.404 · 0.876 · 0.195 · 0.553 · 2.203

IMERG LR · 0.1° / 30 min · 12 小时 · 0.421 · 0.887 · 0.244 · 0.590 · 2.134

CMORPH · 8 km / 30 min · 3—4 个月 · 0.401 · 0.904 · 0.255 · 0.597 · 2.118

IMERG FR · 0.1° / 30 min · 3.5 个月 · 0.418 · 0.885 · 0.307 · 0.613 · 2.042

**Huayu-Global** · **0.05° / 30 min** · **实时** · **0.441** · 0.896 · **0.362** · 0.605 · **1.961**

Huayu-Global 在 CSI、R²、RMSE 三项上都是最优;ACC 0.896 接近 CMORPH 的最高值 0.904,CC 0.605 接近 IMERG FR 的最高值 0.613。相对 GSMaP_NOW,CSI 提高 34.9%、RMSE 降低 31.8%。它是被评估的实时产品里唯一一个 CSI 高于全部基线(包括延后数月发布的产品)的。

分扇区(记录数加权):

扇区 · 记录数 · CSI · ACC · R² · RMSE (mm)

Americas(GOES ABI) · 302,755 · 0.495 · 0.814 · 0.339 · 2.452

Europe-Africa(Meteosat SEVIRI) · 135,130 · 0.398 · 0.888 · 0.247 · 1.597

Asia-Pacific(FY-4B AGRI) · 1,410,082 · 0.415 · 0.914 · 0.356 · 1.866

准全球 · 1,847,967 · 0.441 · 0.896 · 0.362 · 1.961

美洲扇区的 CSI 最高(0.495)同时 RMSE 也最高(2.452 mm),原文解释为该区域含强对流降水、事件探测强但量级误差大;欧洲—非洲扇区 CSI 最低(0.398)但 RMSE 最低(1.597 mm),属更保守的探测配更小的量级误差。

一个技术细节值得注意:Huayu-Global 是用 IMERG FR(0.1°)监督的,却在台站评估上超过它。原文的解读是,格点化的监督产品不必定义在不同空间尺度上评估时的性能上限——0.05° 网格、显式分离降雨有无与降雨量、高频云演变序列这三点让它能把雨区边界、时间和局地量级做得比 0.1° 的训练目标更锐利。另外,准全球的 R² 是在全部合并记录的范围上算的,不是各扇区 R² 的记录数加权平均,因此可能高于任一扇区值。

数据规模:原始静止观测与基准产品合计约 59.7 TB,其中 GOES-16 9.2 TB、GOES-18 11 TB、GOES-19 2.4 TB、Meteosat-10/11 14 TB、Meteosat-9 5.7 TB、FY-4B 12 TB;地面验证用 5,395 个 HadISD 与中国气象局站点。

原文图

数据处理流程
数据处理流程

Huayu-Global 的流程图:GOES、Meteosat 与 FY-4B 的扇区视图先分别转成 ABI、SEVIRI、AGRI 的云特征序列,再各自进入对应扇区的 Huayu 模型,最后把扇区输出拼成一张准全球网格降水场;整个过程只用推理时刻已经可用的卫星输入。(原文 Figure 1)

单个降水个例的六产品空间对比
单个降水个例的六产品空间对比

一个个例的准全球降水场对比,面板依次为 (a) GSMaP_NOW、(b) PERSIANN-CCS、(c) PERSIANN、(d) CMORPH、(e) IMERG FR、(f) Huayu-Global,每幅都配了四处放大区域,站点标记是放大区内的验证站位置。面板上标注的 CSI、CC、R² 与 RMSE 只对这一个上演算,不能当作正文表格里的全样本成绩。(原文 Figure 3)

亚洲验证域的站点级相关系数
亚洲验证域的站点级相关系数

亚洲验证域上以站点为单位的相关系数,站网主要集中在东亚。面板 a—f 依次为 GSMaP_NOW、PERSIANN-CCS、PERSIANN、CMORPH、IMERG FR、Huayu-Global,每个圆点是一个验证站,颜色是由 2022 年 7—12 月与 2025 年 1—6 月全部可用的三小时配对记录算出的皮尔逊相关系数,色标从 0(青色)到 1(品红)。(原文 Figure 4)

作者承认的限制

一句话判断

把三条静止卫星数据流压进同一套架构、每半小时出一次 0.05° 场,并且在不牺牲时效的前提下把 CSI 做到 0.441(高于延后 3.5 个月的 IMERG FR 0.418)、RMSE 压到 1.961 mm,这是对「实时 vs 准确」这个老取舍的实质性松动;不过覆盖止步 60°、非洲南美与洋面验证稀疏、强对流和降雪仍是硬骨头,把它当洪水预警的输入层可以,当全球降水基准还早。

在观潮中打开(含配图与交互)→