GEOS-Chem碳模拟中分子数量变化的数值分析
项目地址: https://gitcode.com/gh_mirrors/ge/geos-chem 背景介绍
在使用GEOS-Chem 14.2.2版本进行碳模拟时，研究人员发现了一个有趣的现象：当关闭化学反应速率和平流过程后，虽然物种的混合比保持恒定，但网格点内的分子数量却出现了变化。这一现象最初是在添加新的NOx物种时发现的，随后在CO和CO2的测试中也得到了验证。
问题现象
在关闭以下所有过程的情况下：
- 化学反应速率（k_Trop/k_Strat设为0）
- 平流传输
- 排放源
理论上，每个网格点内的分子数量应该保持不变。然而实际模拟结果显示：
- 物种混合比（SpeciesConcVV）保持恒定
- 通过空气质量计算或直接从SpeciesConcMND*AIRVOL获得的分子数量却发生了变化
问题诊断
经过深入分析，发现导致这一现象的主要原因是忽略了边界层混合（PBL mixing）过程的影响。边界层混合会改变混合比，同时也会改变空气质量量。当关闭PBL混合后，分子数量的波动显著减小。
数值分析
即使关闭了所有主要物理化学过程，分子数量仍会出现微小波动。这种波动属于数值噪声范畴，其特点包括：
- 波动幅度通常比数值本身小5-6个数量级
- 是由数值计算中的舍入误差和截断误差引起的
- 在长时间积分过程中难以完全避免
技术建议
对于需要进行类似敏感性测试的研究人员，建议：
- 确保关闭所有可能影响物种浓度的过程，包括：
- 化学反应
- 平流传输
- 对流过程
- 边界层混合
- 使用瞬时输出（instantaneous output）而非时间平均输出
- 设置高时间分辨率（如每小时）的输出频率
- 重点关注SpeciesConcMND诊断量
结论
在GEOS-Chem模型中，当关闭所有物理化学过程后观察到的分子数量微小变化属于正常数值现象。这种波动不会影响模型的整体性能和研究结论的有效性。研究人员在进行类似测试时，应全面考虑所有可能影响物种浓度的过程，并正确理解数值模拟中固有的微小波动现象。
创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
