WRF-Chem新手避坑指南：从零开始配置namelist.input，搞定化学和气溶胶模拟

第一次打开WRF-Chem的namelist.input文件时，面对密密麻麻的参数选项，很多初学者都会感到无从下手。作为一款强大的大气化学与气溶胶耦合模型，WRF-Chem的配置复杂度远超基础WRF，一个参数设置不当就可能导致数天的模拟计算前功尽弃。本文将带你系统梳理关键配置逻辑，避开常见陷阱，快速掌握适合自己研究需求的namelist.input配置方法。

1. 化学机制选择：chem_opt的决策逻辑

chem_opt是WRF-Chem中最重要的参数之一，它决定了模型将使用哪种化学机制和气溶胶方案。选择不当会导致模拟结果与实际情况偏差巨大，甚至直接报错终止。以下是主流选项的适用场景对比：

选择建议：

• 研究PM2.5等气溶胶污染时，至少需要选择包含MADE/SORGAM(chem_opt=2)或MOSAIC(chem_opt=4)的方案
• 模拟臭氧生成时，CBMZ机制(chem_opt=6)对VOCs的处理更精细
• 跨国界传输研究推荐MOZART机制(chem_opt=8)

注意：chem_opt一旦确定，后续的emiss_opt、bio_emiss_opt等排放参数必须与之兼容，否则会导致物种不匹配错误。

2. 排放源配置：避免输入数据与参数不匹配

排放源配置是WRF-Chem中最容易出错的环节之一，常见问题包括排放物种与化学机制不匹配、时间分辨率设置不当等。以下是关键参数配置指南：

2.1 人为排放(emiss_opt)

emiss_opt = 3 # RADM2/MADE-SORGAM机制下的典型设置 io_style_emissions = 2 # 使用时间特定的排放文件 auxinput5_interval = 3600 # 每小时更新排放数据(适用于NEI等高频排放清单)

避坑要点：

• 使用NEI排放清单时，确保emiss_opt与chem_opt匹配：
  • chem_opt=2 → emiss_opt=3
  • chem_opt=6 → emiss_opt=4
• 全球排放数据(如EDGAR)通常需要设置emiss_inpt_opt=101进行物种转换

2.2 生物排放(bio_emiss_opt)

bio_emiss_opt = 3 # 使用MEGAN在线计算生物排放 bioemdt = 60 # 生物排放更新间隔(分钟) ne_area = 41 # 必须大于实际化学物种数

常见错误：

• 忘记设置ne_area或设置值过小，导致MEGAN计算异常
• 生物排放时间步长(bioemdt)过长，影响异戊二烯等活性VOCs的日变化模拟

2.3 特殊排放源配置

# 海盐排放 seas_opt = 1 # 启用GOCART海盐排放 dust_opt = 1 # 启用沙尘排放(需提供土壤侵蚀数据) # 生物质燃烧 biomass_burn_opt = 1 # 包含生物质燃烧排放 plumerisefire_frq = 180 # 羽流上升计算频率(分钟)

3. 光解与化学过程：时间步长的艺术

化学过程的时间步长设置需要平衡计算精度和效率，不当配置要么导致计算不稳定，要么大幅延长模拟时间。

3.1 关键时间参数

chemdt = 5.0 # 化学步长(分钟)，通常为气象步长的3-6倍 photdt = 30.0 # 光解频率(分钟) phot_opt = 1 # 使用TUV光解方案

经验法则：

• chemdt建议设置为metdt(气象步长)的3-6倍
• 高分辨率模拟(<3km)需要减小chemdt(2-3分钟)
• 臭氧光化学模拟建议photdt≤30分钟

3.2 化学过程开关

gaschem_onoff = 1 # 开启气相化学 aerchem_onoff = 1 # 开启气溶胶化学 wetscav_onoff = 1 # 开启湿清除过程

调试技巧：可以临时关闭某些化学过程(gaschem_onoff=0)来隔离问题来源。

4. 初始与边界条件：避免"边缘效应"

化学初始场和边界条件的处理不当会导致模拟初期异常或边界处浓度突变。

4.1 初始场配置

chem_in_opt = 1 # 使用外部初始场 have_bcs_chem = .true. # 使用外部边界条件 auxinput12_interval = 21600 # 初始化学场输入间隔(秒) io_form_auxinput12 = 2 # WRF netCDF格式

推荐做法：

• 使用全球模型(如GEOS-Chem)输出作为初始和边界条件
• 至少设置5天的spin-up时间消除初始场影响

4.2 边界条件参数

gas_bc_opt = 1 # 使用默认气体边界轮廓 aer_bc_opt = 1 # 使用默认气溶胶边界轮廓 gas_ic_opt = 1 # 使用默认气体初始场 aer_ic_opt = 1 # 使用默认气溶胶初始场

5. 实战案例：典型城市空气污染模拟配置

以下是一个模拟华北地区PM2.5和臭氧的完整namelist.input化学部分配置：

&chem chem_opt = 4 # CBMZ/MOSAIC机制 emiss_opt = 4 # CBMZ/MOSAIC人为排放 bio_emiss_opt = 3 # MEGAN生物排放 dust_opt = 0 # 关闭沙尘排放(城市区域) seas_opt = 0 # 关闭海盐排放 chem_in_opt = 1 # 外部初始场 have_bcs_chem = .true. # 外部边界条件 chemdt = 5.0 # 化学步长5分钟 bioemdt = 30.0 # 生物排放30分钟更新 photdt = 30.0 # 光解30分钟更新 phot_opt = 1 # TUV光解方案 gaschem_onoff = 1 # 开启气相化学 aerchem_onoff = 1 # 开启气溶胶化学 wetscav_onoff = 1 # 开启湿清除 aer_ra_feedback = 1 # 气溶胶-辐射反馈 ne_area = 41 # MEGAN物种数 io_style_emissions = 2 # 时间特定排放文件 auxinput5_interval = 3600 # 人为排放每小时输入 auxinput12_interval = 21600 # 化学初始场6小时输入 /

验证技巧：

1. 先进行24小时短模拟，检查各化学物种浓度是否在合理范围内
2. 特别关注NOx/VOCs比值，异常值通常表明排放配置错误
3. 使用ncdump检查输入文件变量名与namelist设置是否匹配