Linking GEOS-Chem to CMAQ
This page contains information about how to link the GEOS-Chem with the CMAQ regional air quality model.
- 1 Overview
- 2 pncglobal2cmaq
- 3 GEOSCHEM4CMAQ
- 4 Source code
- 5 Validation
- 6 References
- 7 Known issues
- 8 CMAQv5.0
Global models are increasingly used as time-variable boundary conditions to regional models. Previously, representative concentrations (from models or observations) were used to construct time-invariant boundaries. This prevents important variability in the contribution of grid-internal and grid-external concentrations. By linking to global models, regional models conceptually improve the grid-external sources and provide their own estimates of grid-internal sources. To connect regional and global models, there are conceptually 4 steps:
- Horizontal Interpolation - Since most regional models use a curvilinear grid, this is the simplest step.
- Vertical Interpolation - Many regional models use a pure sigma-pressure coordinate, which is different from GEOS-Chem's hybrid eta-sigma.
- Tracer translation - GEOS-Chem uses explicit and lumped chemical species (gas and particle), which require one-to-one, many-to-one, and one-to-many mapping to regional model tracers.
- Unit conversions - Translated tracers must be converted to regional model units to ensure compatibility.
- Meta Data - Regional model meta-data must be provided to ensure compatibility.
Any utility that accomplishes these steps can be used.
pncglobal2cmaq is part of the PseudoNetCDF software package (http://github.com/barronh/pseudonetcdf). This is compatible with chemistry and SOA aerosols from GEOS-Chem v8, v9 and v10. By default, it is configured with v9. It can easily be adapted to different versions by editing the mappings.json file (see below for more details).
- Horizontal Interpolation - nearest neighbor with automatic conversion of CMAQ coordinates
- Pressure Interpolation - pressure vertical interpolation (recommended); pure GEOS-Chem eta to pure CMAQ sigma is an option if met files are unavailable.
- Tracer translation - Gas-phase translations for CB05 and SAPRC07 are provided; Aerosol translations to CMAQ AE5 and AE6 are provided; custom translations are also easy
- Unit conversion - dynamically converts units
- Meta Data - adds CMAQ meta-data to outputs automatically.
A typical process has 3 steps and 2 pre-requisites and is described in the figure below. The pre-requisites are at least one CMAQ 3-D meteorology input for the inner domain (METCRO3D; only required for initial conditions) and one for the boundary (METBDY3D).
Run GEOS-Chem with ND49 enabled for the domain of interest
This section provides a minimal description of the GEOS-Chem simulation options necessary for pncglobal2cmaq.
- Enable ND49 for all layers (safer than trying to subset) with hourly (if not hourly, use the --tstep=HHMMSS option later).
- Define the domain to be inclusive of your domain (more is better).
- Output all tracers and, with v8, output pure NO, NO2 and O3.
- Output PSURF and TMPU -- these will be used for unit conversion and vertical interpolation. This can be problematic with SOA and versions v9-02 and higher where there are more than 90 tracers. If you run these diagnostics separately, they can be combined later (not described here).
Refer to the diagnostics appendix of your GEOS-Chem version to ensure that correct tracers are being output. I recommend running one day and then reviewing the output tracers using pncdump as described below. Careful review will prevent long delays later.
pncdump --header nd49/ts20100227.bpch | less
Process one or more days ND49 file
Prepare by creating a mapping definitions file
pncglobal2cmaq.py --template=cb05tucl_ae6_aq /path/to/anyfile > mappings.json
- "--template" can be set to cb05tucl_ae5_aq, cb05tucl_ae6_aq, saprc07_ae5_aq, saprc07_ae6_aq
- mappings.json can be modified as desired.
pncglobal2cmaq.py nd49/ts20100227.bpch --METBDY3D met/METBDY3D_20100227 --METCRO3D met/METCRO3D_20100227 --mapping=mappings.json --outifolder=icon/ --outbfolder=bcon/
- Folders nd49, met, bcon, and icon are not required -- adjust paths to be consistent with your setup.
- "mappings.json" is the path that you created in the last step.
- "--METCRO3D" is only necessary to create an initial condition file and can be omitted for any
- "--METBDY3D" provides temperature and pressure for vertical interpolation and unit conversion.
Combine for CMAQ Simulation
CMAQ boundary files must have 25 instantaneous times per CMAQ simulation day. GEOS-Chem ND49 files typically have 24 times. Combine multiple days of output to ensure coverage (too much is okay) of all instantaneous times in CMAQ simulation. It is typical to use monthly BCON files for long simulations.
Combine Using PseudoNetCDF
pncgen --stack=TSTEP bcon/ts20100227.bpch.BCON.nc bcon/ts20100228.bpch.BCON.nc [...] bcon/ts2010041.bpch.BCON.nc bcon/ts201003.BCON.nc
Combine Using NCO
ncrcat bcon/ts20100227.bpch.BCON.nc bcon/ts20100228.bpch.BCON.nc [...] bcon/ts2010041.bpch.BCON.nc bcon/ts201003.BCON.nc
Output files can be visualized using the PseudoNetCDF pncboundaries.py script.
With fewer options:
pncboundaries.py bcon/ts20060801.bpch.BCON.nc -v O3 -r TSTEP,mean bcon.ts201003.BCON.nc
With more options:
pncboundaries.py bcon/ts20060801.bpch.BCON.nc --states=True --norm="BoundaryNorm(np.logspace(0.025, 1, 9), 256) -v O3 -r TSTEP,mean bcon.ts201003.BCON.nc
Below is a comparison of the CMAQ v5.0.2 standard BCON from "profiles" compared to the GEOS-Chem derived BCON. Note in GEOS-Chem, increasing ozone with altitude is seen as expected. PROFILE-based BCON has much lower ozone at high altitude because the data was collated before CMAQ top pressures were typically as high as they are now. Using the old profiles leads to exceptionally low ozone concentrations entering the domain in the upper boundaries.
--Barron H. 2016-09-08 13 (EST)
Tao Zeng wrote:
- The code to convert GEOS-Chem simulations as the boundary and initial conditions for CMAQ CB05 is ready. I have tested it successfully on GEOS-Chem v08-01-04 using pgf90 compiler. Some highlights are:
- (1) 2 output options:
- (a) in CMAQ domain after regridding. The outputs are CMAQ ready BC/IC files with NetCDF format. The MCIP met file (GRIDCRO3D and GRIDCRO2D) and default CMAQ IC/BC are required. IOAPI and NetCDF libs are also needed.
- (b) in G-C subdomain without regridding. It is in binary format.
- (2) Conversion of GOS-Chem species to CB05 tracers
- Tracer mapping table are designed to convert GEOS-Chem species to CMAQ tracers. Rob Pinder from EPA helped me to review it.
- CSPEC array is used to get direct NO and NO2 concentrations. Some other species in CSPEC but not in STT are also included. Several inactive species (ACTA, EOH, HCOOH, and MOH) are turned on into active species to get better mapping of CB05 tracers. This resulting differences are <1% in average over the domain after 1 month integration. And the max difference is ~2% for most tracers. But for NH3 and NIT, the max differences are up to 50% when they are at very low level.
- (3) To avoid excess downward transport from stratosphere in CMAQ, tracer concentrations at the tropopause in GEOS-Chem are used to for the tropopause and the layers above.
- Any suggestions are welcome. The whole package includes a new fortran module file and some lines adding into 3 files. If anybody is interested, please let me know.
--Bob Y. 10:19, 10 March 2010 (EST)
Authors and collaborators
- Tao Zeng (Georgia EPD)
|Georgia EPD||Tao Zeng||...|
|EPA||Havala Pye||Boundary conditions with a focus on aerosols|
|Add yours here...||...||...|
Tao Zeng has provided the GEOS-Chem to CMAQ source code linkage. It is available for download from:
ftp ftp.as.harvard.edu get pub/geos-chem/downloads/geoschem4cmaq.20100316.tgz
To extract the archive, type:
tar xvzf geoschem4cmaq.20100316.tgz
You will find a README.txt file in the archive with more information.
--Bob Y. 15:27, 16 March 2010 (EDT)
SOA1, SOA2, SOA3, and SOA4 in G-C are outputed but are not converted to CB05 tracers. In CMAQ v4.7, 19 SOA species are defined. The CMAQ release notes do not fully explain them. So I will update this package when possible.