WRF-Chem区域大气化学模拟：从环境搭建、排放源处理到结果验证的完整实践指南-平芜编程栈

你有没有过这样的经历：花了几天时间，好不容易在服务器上把WRF-Chem模式编译成功，结果一运行就报错，不是缺库就是路径不对，或者干脆提示“segmentation fault”。你对着满屏的日志，感觉就像在解一道没有答案的谜题。更让人头疼的是，即便模式跑起来了，生成的模拟结果和观测数据对不上，你根本不知道问题出在气象场、化学机制，还是那个看起来最不起眼但影响巨大的排放源清单上。

这几乎是每一个刚接触区域大气化学模拟的研究者或工程师都会遇到的困境。WRF-Chem，这个将气象模式WRF与大气化学过程在线耦合的利器，理论上能为我们揭示从污染物排放、传输、化学反应到最终沉降的全链条过程。然而，从“理论上的强大”到“实践中的可用”，中间隔着一道由Linux环境、依赖库、编译选项、数据制备和模式调试构成的鸿沟。很多人止步于环境配置，更多人卡在了排放源制作这个“黑箱”环节。

这篇文章不会给你一个“一键安装”的神话。相反，我会带你走完从零开始构建一个可运行、可调试、结果相对可靠的WRF-Chem模拟环境的完整流程。我的核心判断是：WRF-Chem的成功应用，关键不在于对某个复杂化学机制的深入理解，而在于能否系统化地搭建一个稳定、透明、可复现的技术栈，并把最容易出错的“数据准备”环节，尤其是排放源处理，做到足够扎实和可验证。我们将把这个过程拆解为五个环环相扣的部分，每一步都力求解释清楚“为什么这么做”以及“踩坑了怎么办”。

1. 理解核心：为什么是“在线耦合”？WRF-Chem解决了什么根本问题？

在深入命令行之前，我们必须先搞清楚WRF-Chem的独特价值。传统的大气污染模拟常常采用“离线”模式：先用WRF等气象模式模拟出气象场（如风、温、湿、压），再将结果输出为文件，最后用另一个化学传输模型（如CMAQ）读取这些气象文件来模拟污染物过程。这种方法存在一个天然的“断层”：化学过程无法实时反馈影响气象场。

WRF-Chem的“在线耦合”彻底改变了这一范式。它意味着气象模块和化学模块在同一个时间步长内协同运算。化学物种的浓度变化会影响辐射过程，从而改变气温和边界层结构；而气象场的变化又立刻作用于污染物的平流、扩散和化学反应速率。这种双向反馈对于模拟二次气溶胶生成、雾-霾相互作用、沙尘-辐射效应等过程至关重要。

举个例子，在重霾事件中，气溶胶颗粒会散射和吸收太阳辐射，导致近地面降温、大气层结稳定，进而抑制污染物的垂直扩散，形成“越脏越稳，越稳越脏”的正反馈。离线模式很难准确捕捉这种耦合效应，而WRF-Chem在架构上就为此而生。

所以，当你决定使用WRF-Chem时，你本质上是在选择处理一类气象与化学存在强相互作用的复杂环境问题。这决定了你后续所有的工作，从化学机制选择到排放源的时间分辨率要求，都需要以支持这种双向耦合为目标。

2. 奠基：构建一个专为WRF-Chem服务的Linux编译环境

WRF-Chem的编译是第一个实战关卡。它失败的原因90%不在于模式本身，而在于基础环境不纯净、不完整或不一致。

2.1 操作系统的选择与准备：稳定压倒一切

不要追求最新的Linux发行版。WRF/Chem的代码库和众多依赖库（NetCDF, HDF5, MPI等）对编译器版本非常敏感，最新系统自带的编译器可能引入兼容性问题。

推荐选择：CentOS 7/Rocky Linux 8、Ubuntu 20.04 LTS 或 22.04 LTS。这些系统有长期支持，社区资源丰富，软件库稳定。
关键一步：在安装完系统后，第一件事是更新软件源并安装基础开发工具包。
- 对于基于RHEL的系统（如Rocky Linux）：
```
sudo yum groupinstall "Development Tools" -y sudo yum install wget git cmake curl gcc-gfortran -y
```
- 对于基于Debian的系统（如Ubuntu）：
```
sudo apt update sudo apt install build-essential gfortran git wget curl cmake -y
```
这个Development Tools或build-essential套件包含了编译所需的make、gcc、g++等核心工具，避免后续出现“命令未找到”的错误。

2.2 依赖库的编译：顺序与配置决定成败

WRF-Chem依赖几个关键库，必须手动编译，且顺序有讲究。一个黄金法则是：所有库尽量使用相同的编译器（如GCC）和相同的安装前缀（/opt或/usr/local），并开启并行支持（MPI）。

我们假设将库安装在/opt目录下，需要sudo权限。也可以安装到用户目录（如~/lib），避免权限问题。

标准编译流程如下：

环境变量设置：在~/.bashrc文件中设置以下变量，然后执行source ~/.bashrc。

export DIR=/opt # 库安装目录 export CC=gcc export CXX=g++ export FC=gfortran export F77=gfortran export FFLAGS=-m64 export CFLAGS=-m64 export PATH=$DIR/bin:$PATH export LD_LIBRARY_PATH=$DIR/lib:$LD_LIBRARY_PATH export CPPFLAGS="-I$DIR/include" export LDFLAGS="-L$DIR/lib"

编译安装依赖库（按顺序）：

MPI(如OpenMPI)：WRF-Chem并行计算的基础。

wget https://download.open-mpi.org/release/open-mpi/v4.1/openmpi-4.1.5.tar.gz tar -zxvf openmpi-4.1.5.tar.gz cd openmpi-4.1.5 ./configure --prefix=$DIR make -j$(nproc) # 使用所有CPU核心加速编译 sudo make install

安装后，将MPI的路径加入环境变量：export PATH=$DIR/bin:$PATH。

NetCDF-C & NetCDF-Fortran：气象和化学数据存储的绝对标准。必须先装C库，再装Fortran库。

# 安装NetCDF-C wget https://github.com/Unidata/netcdf-c/archive/refs/tags/v4.9.2.tar.gz -O netcdf-c-4.9.2.tar.gz tar -zxvf netcdf-c-4.9.2.tar.gz cd netcdf-c-4.9.2 CPPFLAGS=-I$DIR/include LDFLAGS=-L$DIR/lib ./configure --prefix=$DIR --disable-dap make -j$(nproc) sudo make install # 安装NetCDF-Fortran wget https://github.com/Unidata/netcdf-fortran/archive/refs/tags/v4.6.1.tar.gz -O netcdf-fortran-4.6.1.tar.gz tar -zxvf netcdf-fortran-4.6.1.tar.gz cd netcdf-fortran-4.6.1 LD_LIBRARY_PATH=$DIR/lib:$LD_LIBRARY_PATH CPPFLAGS=-I$DIR/include LDFLAGS=-L$DIR/lib ./configure --prefix=$DIR make -j$(nproc) sudo make install

JASPER(用于GRIB2数据支持)：如果你需要处理来自GFS、ECMWF等机构的GRIB2格式气象初始场和边界场，这是必须的。

wget https://github.com/jasper-software/jasper/releases/download/version-4.0.0/jasper-4.0.0.tar.gz tar -zxvf jasper-4.0.0.tar.gz cd jasper-4.0.0 mkdir build && cd build cmake .. -DCMAKE_INSTALL_PREFIX=$DIR make -j$(nproc) sudo make install

2.3 验证环境：编译前的“健康检查”

在下载WRF-Chem源码之前，用一个简单的方法验证环境：

nc-config --all # 查看NetCDF配置，确认Fortran接口已启用 which mpif90 # 确认MPI Fortran编译器可用

确保这些命令都能正确输出路径和配置信息，没有报错。

3. 攻坚：WRF-Chem模式的编译与配置选择

环境就绪后，编译WRF-Chem本身反而相对直接，但配置选择至关重要。

3.1 获取源码与基础配置

从官方或GitHub获取稳定版本的WRF和WPS（WRF预处理系统）代码。通常两者版本需匹配。

wget https://github.com/wrf-model/WRF/archive/refs/tags/v4.5.tar.gz wget https://github.com/wrf-model/WPS/archive/refs/tags/v4.5.tar.gz tar -zxvf v4.5.tar.gz tar -zxvf v4.5.tar.gz cd WRF-4.5

运行配置脚本：

./configure

这时会出现一个交互菜单让你选择编译选项。这是第一个决策点：

选择编译器：选择与你编译依赖库时一致的选项，例如Linux x86_64, gfortran (gcc)配合dmpar（分布式内存并行，即MPI）。
选择嵌套：对于区域化学模拟，通常从1（基本嵌套）开始。高级嵌套选项用于更复杂的需求。
选择化学机制：这是WRF-Chem的核心。菜单会列出如RADM2、RACM、MOZART、SAPRC99、CB05等选项。你的选择应基于：
- 研究目标：模拟臭氧？选包含详细VOCs机制的（如SAPRC99）。模拟气溶胶？需包含气溶胶模块的机制（如MOZART with MOSAIC/GOCART）。
- 排放清单匹配：你准备的排放源数据必须能映射到所选化学机制的物种上。这是后续排放源制作的最大挑战。
- 计算成本：机制越复杂，物种和反应越多，计算越慢。初学者可从RADM2或CB05这类中等复杂度的机制入手。

选择后，配置脚本会生成一个configure.wrf文件。务必检查这个文件，确认其中NETCDF、MPI等路径指向你之前安装的/opt目录。

3.2 执行编译与常见问题排查

执行编译：

./compile em_real 2>&1 | tee compile.log

em_real代表真实案例（ideal案例用于测试）。tee命令将屏幕输出同时保存到compile.log文件，这是排查编译错误的生命线。

经典错误排查：

找不到 -lnetcdff 或 -lnetcdf：说明NetCDF库路径未正确链接。检查configure.wrf中的LIB_EXTERNAL部分，手动添加-L/opt/lib。
MPI相关错误：确认mpif90命令可用，且环境变量LD_LIBRARY_PATH包含了MPI库路径。
段错误 (Segmentation fault) 在编译期：可能是编译器优化过强。尝试在configure.wrf的FCOPTIM和CCOPTIM标志中降低优化等级，如将-O3改为-O2甚至-O0（仅用于调试）。
内存不足：编译WRF-Chem需要大量内存。如果服务器内存较小，可以尝试./compile -j 1 em_real单线程编译，或增加虚拟内存。

编译成功的标志是在main目录下生成wrf.exe（主程序）和real.exe（前处理程序）。同样地，进入WPS目录进行编译（需要链接JASPER库）。

4. 灵魂：排放源清单的制作、处理与质量控制

如果说模式是躯体，那么排放源就是灵魂。一个糟糕的排放清单会直接导致模拟结果失真，而这个问题在后期极难调试。

4.1 排放源数据的获取与理解

排放源通常包括人为源（工业、交通、居民、电厂等）和自然源（生物源、扬尘等）。常见数据来源有：

全球清单：如EDGAR、CEDS、MIX。分辨率较粗（如0.1度），适合作为背景或大范围模拟。
区域/国家清单：如中国的MEIC、美国的NEI。分辨率更高，物种更贴合本地实际。
卫星反演：如OMI NO₂、MOPITT CO，可用于约束或验证。

拿到数据后，首先要做的是数据对齐：

空间范围与分辨率：你的模拟域是否能覆盖清单范围？是否需要裁剪或插值？
时间分辨率：清单是年排放量、月均值，还是包含日变化和小时变化？WRF-Chem需要小时级别的排放输入。如果没有，你需要借助时间分解因子（如每月、每周、每日的权重系数）将年总量分解到每个模拟小时。
物种映射：这是最关键的一步。你的清单物种（如NOx、SO2、NMHC）必须精确地映射到你在步骤3.1中选择的化学机制所要求的物种上。例如，NMHC（非甲烷总烃）需要按照特定的VOC物种谱进行拆分。这一步通常需要编写脚本（Python/IDL/NCL）进行处理。

4.2 使用PREP-CHEM-SRC或ANTHRO_EMISS工具进行处理

原始清单数据不能直接用于WRF-Chem。需要借助工具将其转换为WRF-Chem可读的NetCDF格式，并重映射到模拟网格上。

WRF-Chem社区提供了PREP-CHEM-SRC工具包，它是一个集成的排放预处理系统。其工作流程可以概括为：

准备输入：将你的清单数据（如NetCDF或文本格式）放在指定目录，并编写一个控制文件（prep_chem_sources.inp），在其中指定：
- 模拟域信息（与namelist.wps一致）
- 使用的化学机制
- 各种排放源文件路径和类型
- 时间处理选项
运行程序：工具会根据控制文件，执行网格重映射、时间插值、物种映射，最终输出一个或多个wrfchemi_*文件（如wrfchemi_00z_d01代表0时的排放）。
质量检查：永远不要假设工具输出就是正确的。必须用ncdump、ncview或Python的xarray/netCDF4库检查输出文件：
- 排放总量是否在合理量级？（例如，中国东部NOx排放强度大概在1e-11 ~ 1e-10 kg m⁻² s⁻¹）
- 空间分布是否符合预期？（大城市、工业区是否呈高值？）
- 时间变化是否合理？（夜间排放是否降低？）

4.3 常见陷阱与调试策略

负值或异常值：检查原始清单数据是否有负值，或映射过程中出现错误计算。
空间错位：确认模拟域和清单数据的投影坐标系是否一致。PREP-CHEM-SRC通常要求清单数据是经纬度坐标。
物种缺失：运行后提示某些化学机制物种找不到排放值。这需要你返回物种映射步骤，检查是否所有必需物种都有对应的输入，或者是否需要启用工具的“缺省填充”选项（如使用背景值）。
单位不匹配：WRF-Chem内部使用kg m⁻² s⁻¹。确保你的处理工具进行了正确的单位换算。

一个实用的调试流程是：从一个极小的、你知道答案的测试案例开始。例如，创建一个只在一个网格点有固定排放量的简单清单文件，处理它，看输出文件是否在那个网格点产生了正确强度和物种的排放。这能帮你快速定位问题是出在工具配置、数据格式还是映射逻辑上。

5. 集成与实战：从首次成功运行到结果分析与问题溯源

当模式编译通过、排放源准备就绪，就可以进行第一次完整的模拟了。

5.1 配置与运行流程

气象数据准备：使用编译好的WPS程序，将GFS、ERA5等全球再分析数据插值到你的模拟网格，生成met_em.*文件。
编辑namelist：这是模式运行的“大脑”。需要仔细配置三个namelist文件：
- namelist.wps：定义模拟域（中心点、网格数、分辨率、嵌套）、地图投影、以及气象输入数据路径。
- namelist.input：核心控制文件。定义积分时间步长、物理参数化方案（微物理、积云、边界层、辐射、陆面过程）、化学选项（化学机制、排放输入频率、光解速率方案等）、输入输出设置。
- namelist.chem(WRF-Chem专用)：进一步控制化学过程，如排放缩放因子、干湿沉降选项、气溶胶方案等。
执行real.exe：将初始场、边界场和排放文件处理成WRF-Chem可读的格式。
执行wrf.exe：开始主积分计算。建议首次运行时，先模拟很短的时间（如6-12小时），并使用debug_level=300输出更详细的日志，以便快速发现问题。

5.2 结果分析与初步验证

运行成功后，你会得到一系列wrfout文件。分析的第一步不是画漂亮的图，而是做“合理性检查”：

气象场检查：温度、气压、风场、湿度是否在合理范围内？有没有出现“棋盘格”状的不稳定现象？这可能是物理方案选择不当或时间步长过长。
化学场检查：污染物浓度（如PM2.5、O₃、NO₂）的量级是否合理？例如，近地面O₃浓度通常在几十ppb量级，如果出现几千ppb，肯定是错误的。空间分布是否反映了排放源和风场的影响？
与观测对比：获取模拟区域内环境监测站的小时数据。计算简单的统计指标，如相关系数（R）、平均偏差（MB）、标准化平均偏差（NMB）。首次模拟的对比结果通常不会很好，这很正常。我们的目标是通过对比，定位问题方向：
- 所有物种普遍偏高/偏低：可能指向排放清单的整体缩放因子需要调整。
- 某种物种（如NOx）偏差很大，其他尚可：可能指向该物种的排放清单或化学机制中的相关反应存在问题。
- 日变化趋势错位：可能指向排放清单的时间分解因子不准确，或边界层日变化模拟有问题。

5.3 建立系统化的调试思维

当结果不理想时，避免盲目调整参数。建立一个自上而下的排查顺序：

第一层：输入数据
- 气象初始场/边界场是否质量良好？（用ncview快速浏览）
- 排放源文件是否被正确读取？（检查rsl.error.0000日志中关于读取wrfchemi文件的记录）
- 排放量级是否与文献或常识相符？
第二层：模式配置
- namelist.input中的物理参数化方案组合是否适用于你的区域和季节？（例如，模拟中国冬季雾霾，边界层方案和微物理方案的选择至关重要）
- 化学机制是否适合你的主要模拟物种？
- 时间步长是否满足CFL稳定性条件？（通常，水平分辨率（km）* 6 可作为最大时间步长（秒）的参考）
第三层：过程分析
- 如果浓度异常，是排放输入的问题，还是化学转化、干湿沉降、垂直混合过程的问题？可以输出模式的过程分析变量来诊断。