基于WRF-Chem的不同气溶胶的辐射效应的研究

前言

目前我对于气溶胶辐射效应的理解就是设计敏感性实验，基础实验打开气溶胶参与辐射开关（aer_ra_feedback），其他的实验则关闭气溶胶参与辐射过程开关，也有去掉某些气溶胶的影响，如黑碳（BC）。前者好说，后者目前有两种思路，特此记录！

一、修改WRF-Chem源码

据说（我没实操过····）气溶胶在物理和化学模块中都会参与计算，而控制计算的相关文件在：

chem/module_optical_XXX.F
chem/module_aerosols_sorgam.F
phys/module_radiation_driver.F

这里面也许会有：

ext_total = ext_sulfate + ext_oc + ext_bc + ···

可以手动注释掉某个物种，然后重新编译WRF-Chem，听起来也不是很难对吧，但先不说能不能找到辐射模块文件，就重新编译WRF-Chem这一步就让我望而却步了···这操作需要懂Fortran，且在编译的过程中出现任何一个问题我可能就崩溃了，所以果断放弃！

二、修改人为排放源数据

这个一听就很棒~简单来说就是在输出人为排放源时，将某类气溶胶（如BC）的数值设置为0，然后模拟结果与基础实验做差值，这不就是BC的影响吗？但是否可能还有待研究。

三、离线辐射计算（ChatGPT-5）

• **离线辐射计算（post-processing）**是一个非常有用的替代方案：

  • 将WRF-Chem输出的三维气溶胶物种场（质量浓度、粒径信息等）导出（保存为每时刻的剖面或格点），用独立的辐射传输程序（例如 RRTMG diagnostic / libRadtran / 自己的辐射计算脚本）计算短波/长波通量，分别把不同物种的光学贡献打开/关闭，得到纯辐射通量差异。

  • 优点：不改主码、不受并行/编译影响；可以快速做很多敏感性组合；更易调参和可重复。

  • 缺点：这种方法得到的是瞬时/离线的辐射差，要把它和耦合模式中气温/动力的长时响应联系起来需要额外步骤（例如把离线辐射差作为强迫项用于单独的气候响应模拟，或估算表面/大气强迫）。

无奈~看不懂~不了解~我是🐽