Noah-MP 陆面过程模型辐射输入：3种缺失数据场景下的向下短波辐射估算方案

2026/7/8 23:44:51

Noah-MP陆面模型辐射数据缺失场景的工程化解决方案

当我们在使用Noah-MP这类陆面过程模型进行地表能量平衡模拟时，向下短波辐射(DSR)作为关键输入参数，其准确性直接影响模拟结果的可信度。然而实际工作中，研究人员常会遇到观测数据不完整的情况——或是关键分量全为填充值，或是部分时段数据缺失，亦或是仅有总辐射辐照度这一单一指标。面对这些工程实践中真实存在的困境，我们需要建立系统化的数据预处理流程。

1. 辐射分量缺失问题的本质剖析

陆面过程模型中的向下短波辐射通常由直接辐射和漫射辐射组成，其波长范围集中在0.15-4微米之间。理想情况下，我们应该使用实测的直接辐射辐照度(Edirect)和漫射辐射辐照度(Ediffuse)相加得到DSR：

DSR = Edirect + Ediffuse

但现实中的数据困境往往呈现三种典型场景：

全填充值场景：直接辐射和漫射辐射字段全为无效值（如999999）
部分缺失场景：某些时段或站点缺少其中一个辐射分量
单一总辐射场景：仅能获取总辐射辐照度(Etotal)这一综合指标

这些数据质量问题会导致模型输入的不确定性增加，进而影响能量平衡计算的精度。特别是在研究地表蒸散发、土壤温度剖面等对辐射敏感的物理过程时，不准确的DSR输入可能造成模拟结果的系统性偏差。

2. 三种缺失场景的工程应对策略

2.1 全填充值场景的解决方案

当直接辐射和漫射辐射均为无效值时，可考虑以下替代方案：

总辐射反演法：利用可获取的总辐射辐照度减去估算的地表反射辐射
```
# 假设已知地表反照率(albedo)和总辐射 DSR_estimated = Etotal * (1 - albedo)
```
参数化方案：结合气象要素建立经验关系
- 使用日照时数、云量等辅助数据
- 参考FAO Penman-Monteith公式中的辐射估算方法
再分析数据替代：融合ERA5、GLDAS等再分析数据集

关键参数对比表：

方法	所需输入	优点	局限性
总辐射反演	Etotal, albedo	计算简单	依赖反照率精度
参数化方案	常规气象数据	物理基础明确	需要本地化校准
再分析数据	无	时空连续	分辨率较低

2.2 部分缺失场景的插补技术

当仅缺失部分辐射分量时，可采用更精细的插值方法：

物理模型辅助插值：

# 基于晴天辐射模型估算缺失的直接辐射 def clear_sky_direct(doy, lat, lon): # 实现晴天辐射模型计算 return Edirect_clear

统计相关法：
- 建立完整时段直接/漫射辐射的比例关系
- 使用机器学习模型预测缺失分量
时间序列分析法：
- 应用ARIMA等模型填补短时缺失
- 考虑日周期、季节周期等特征

提示：部分缺失场景下，建议优先使用基于物理的方法进行插补，以保持能量守恒特性。

2.3 单一总辐射场景的处理流程

仅有总辐射数据时，可采用分步处理策略：

质量控制阶段：
- 剔除明显异常值（如夜间非零值）
- 检测仪器漂移问题

辐射分解阶段：

# 基于Erbs模型分解总辐射 def decompose_radiation(Etotal, zenith): # 实现辐射分解算法 return Edirect, Ediffuse

结果验证阶段：
- 对比邻近站点数据
- 检查能量闭合程度

3. 工程实践中的Python实现示例

以下代码展示了基于总辐射和反照率估算DSR的完整流程：

import numpy as np import pandas as pd def estimate_dsr(dataframe, etotal_col='Etotal', albedo_col='albedo', fill_value=999999): """ 估算向下短波辐射的实用函数 参数: dataframe: 包含辐射数据的DataFrame etotal_col: 总辐射列名 albedo_col: 反照率列名 fill_value: 缺失值标记 返回: 包含估算DSR的Series """ # 处理缺失值 valid_mask = (dataframe[etotal_col] != fill_value) & \ (dataframe[albedo_col] != fill_value) # 计算DSR dsr = np.zeros(len(dataframe)) dsr[valid_mask] = dataframe[etotal_col][valid_mask] * \ (1 - dataframe[albedo_col][valid_mask]) # 标记无效值 dsr[~valid_mask] = fill_value return pd.Series(dsr, name='DSR_estimated') # 使用示例 # df['DSR'] = estimate_dsr(df, 'total_radiation', 'surface_albedo')

该实现考虑了实际工程中的常见需求：

处理原始数据中的填充值
避免负值和异常结果
保持与原始数据一致的数据类型

4. 不同方案的误差分析与优化

在实际应用中，我们需要评估各种估算方法的误差特性：

误差来源分析：
- 反照率估算的不确定性
- 辐射分解模型的假设限制
- 时空尺度不匹配问题
优化方向：
- 结合遥感数据改进反照率估计
- 使用观测数据本地化模型参数
- 实施多方案集成方法

验证方法：

def validate_estimation(observed, estimated): # 计算常用验证指标 bias = np.mean(estimated - observed) rmse = np.sqrt(np.mean((estimated - observed)**2)) r2 = np.corrcoef(observed, estimated)[0,1]**2 return {'Bias':bias, 'RMSE':rmse, 'R2':r2}

对于关键应用场景，建议实施以下质量控制步骤：