Open-Meteo API中ERA5高程数据问题的技术解析

问题背景

Open-Meteo是一个提供天气数据API的开源项目，其中包含对ERA5再分析数据集的访问接口。在最近的使用中，开发人员发现当请求参数设置为elevation=nan时，API返回的高程数据存在异常。

问题现象

当用户请求使用ERA5模型的原始网格单元高度（通过设置elevation=nan参数）时，API返回的高程值在某些情况下不正确。具体表现为：

1. 对于高山区域，API错误地返回了0米高程
2. 对于沿海区域，API同样返回0米高程，而实际应为非零值

技术原因分析

经过项目维护者的深入调查，发现该问题由多个因素导致：

1. CDS API数据获取问题：最初修复时，维护者只能从备份中恢复高程数据文件，而无法直接从CDS API获取最新数据。这个备份文件大约有两年历史，且未正确应用重力常数转换。

2. 网格调整问题：在重构CDS API新版代码时，维护者忘记将网格调整为从-180°/-90°开始的坐标系。这一错误在数周内未被发现，因为下载流程通常不涉及地表高程数据的获取。

3. 海陆网格处理：对于水域覆盖率超过50%的网格单元，系统只存储陆地网格的高程数据，而将海洋网格标记为NaN。这导致沿海和湖泊区域的网格单元高程数据丢失。

解决方案

项目维护者实施了以下修复措施：

1. 网格自动检测机制：改进了代码，使其能自动检测需要调整到-90/-180原点的网格，这降低了未来代码的复杂性。

2. 高程数据处理：重新获取并正确处理ERA5高程数据，确保正确应用重力常数转换（g=9.80665 m/s²）。

3. 未来改进计划：考虑使用负值表示海洋高程（如-1000表示0米海拔），以保留水域网格的高程信息。

技术细节

ERA5数据中的高程是通过"geopotential_at_surface"变量计算得出的。根据ECMWF的建议，正确的计算方法是将位势高度除以地球重力加速度（g=9.80665 m/s²）。在早期版本中，错误地使用了g=10 m/s²进行计算，导致了约20米的误差（在500-1200米海拔范围内）。

对用户的影响

这一修复确保了：

1. 高山区域能正确返回非零高程值
2. 温度数据不再进行错误的统计降尺度处理
3. 用户可以信任API返回的高程数据与原始ERA5数据集一致

结论

Open-Meteo团队快速响应并解决了这一数据问题，展示了开源项目对数据质量的重视。对于需要使用原始ERA5网格高程的用户，现在可以放心使用elevation=nan参数获取准确的高程数据。

对于特殊区域（如沿海和湖泊）的高程处理，项目团队将继续优化，以提供更完整的地表高程信息。