随着空气污染和气候变化问题的日益严重,准确模拟大气化学过程成为研究和政策制定的重要支撑工具。WRF-Chem(Weather Research and Forecasting model coupled with Chemistry)是一款集成大气化学和气象模式的数值模拟系统,被广泛用于空气质量预测、气溶胶对气候的影响研究等领域。
本文将详细介绍 WRF-Chem 的基本概念、安装配置、关键物理过程、应用案例以及一些使用技巧,帮助读者更好地理解并应用 WRF-Chem 进行大气化学模拟。
1. 什么是 WRF-Chem?
WRF-Chem 是 WRF(Weather Research and Forecasting)模式的一个扩展版本,专门用于模拟大气化学成分的演变,包括气态污染物、气溶胶和化学过程。WRF-Chem 允许用户在气象模式运行时同时计算空气污染物的输送、化学反应和沉降等过程。
1.1 WRF-Chem 的特点
- 全耦合系统:WRF-Chem 将大气动力学、物理过程和化学过程紧密耦合,使得化学过程能够实时影响气象条件,例如气溶胶对辐射的影响。
- 灵活的化学机制:支持多种化学反应机制,如 RADM2, CBMZ, MOZART, RACM,以及气溶胶模型(如 MADE/SORGAM, MOSAIC)。
- 支持外部排放数据:可以引入全球或区域的排放清单,如 EDGAR, MEIC, NEI,或用户自定义的排放数据。
- 广泛的应用领域:包括空气质量模拟、气溶胶气候效应研究、火灾烟雾影响、城市污染等。
2. WRF-Chem 的安装与配置
由于 WRF-Chem 依赖多个外部库(如 NetCDF, MPICH, HDF5 等),安装相对复杂,建议在 Linux 系统(如 Ubuntu 或 CentOS)下进行。
2.1 依赖环境
在安装 WRF-Chem 之前,需要准备以下软件:
- Fortran & C 编译器(如 GNU 或 Intel 编译器)
- NetCDF(C & Fortran)
- MPICH 或 OpenMPI
- HDF5
- Perl、NCL、Python(可选用于前后处理)
2.2 下载与编译
- 获取 WRF-Chem 源码:
git clone https://github.com/wrf-model/WRF.git cd WRF
- 选择 WRF-Chem 版本(推荐使用最新稳定版)。
- 运行
configure
选择计算架构:./configure
- 编译 WRF:
./compile em_real >& log.compile
- 编译成功后,
WRF/main/
目录下应生成 wrf.exe 和 real.exe。
2.3 测试 WRF-Chem
编译后可以运行 test/em_chem_inp
进行基本测试,以确保 WRF-Chem 能够正确运行。
3. WRF-Chem 的核心物理过程
WRF-Chem 耦合了多个大气化学和物理过程,包括:
- 化学反应机制(Chemical Mechanisms):控制气态污染物的演变,如 CBMZ, RADM2。
- 气溶胶模块(Aerosol Schemes):计算气溶胶生成、成长、沉降的过程,如 MADE/SORGAM, MOSAIC。
- 气象与化学相互作用(Feedbacks):如气溶胶-辐射相互作用(direct effect)和气溶胶-云微物理相互作用(indirect effect)。
- 排放源(Emissions):包括人为排放(交通、工业)和自然排放(生物排放、火灾烟雾)。
4. WRF-Chem 的应用案例
4.1 城市空气质量预测
许多研究机构使用 WRF-Chem 进行 PM2.5、O3 和其他污染物的预测。例如:
- 模拟北京及华北地区的空气污染过程
- 评估治理政策对空气质量的影响
4.2 气溶胶对气候的影响
WRF-Chem 允许研究气溶胶如何改变云和辐射,比如:
- 分析气溶胶导致的温度变化
- 研究火灾烟雾对气候的长期影响
4.3 火灾烟雾传输
WRF-Chem 还可用于火灾烟雾的传输模拟,例如:
- 模拟加州山火烟雾的扩散
- 研究火灾烟雾对人体健康的影响
5. WRF-Chem 的数据前处理
为了成功运行 WRF-Chem,需要准备多个输入数据,如气象数据、排放数据等。
5.1 气象数据
通常使用 WPS(WRF Preprocessing System) 处理全球再分析数据(如 ERA5, GFS)。
5.2 排放数据
常用排放数据源:
- 全球排放数据库(如 EDGAR, HTAP)
- 地区排放清单(如 NEI, MEIC)
- 在线排放模型(如 MEGAN, GFED 火灾排放)
5.3 化学初始条件
使用mozbc工具为 WRF-Chem 生成化学初始和边界条件,可从MOZART或CAM-Chem获取数据。
6. WRF-Chem 运行与后处理
6.1 运行 WRF-Chem
- 运行
real.exe
生成初始气象条件:./real.exe
- 运行
wrf.exe
进行模拟:mpirun -np 8 ./wrf.exe
6.2 结果可视化
- NCL:
ncview
或wrf_python
- Python: 使用
wrf-python
进行可视化 - Paraview/VisIt: 用于3D可视化分析
7. WRF-Chem 模拟优化建议
- 调整物理方案:根据区域特征选择最适合的参数化方案。
- 优化排放数据:确保使用高分辨率、最新的排放清单。
- 增加观测数据同化:结合卫星数据或地面观测改进模拟精度。
- 使用更高分辨率:对于城市污染模拟,可使用 1 km 或更高的分辨率。