Connecticut Regional Haze State Implementation Plan November 2009



### ATTACHMENT T

NYSDEC Technical Support Document TSD-1d1 TSD-1d



## 8hr Ozone Modeling using the SMOKE/CMAQ system

Bureau of Air Quality Analysis and Research

Division of Air Resources

New York State Department of Environmental Conserva tion

Albany, NY 12233

February 1, 20062

Air Quality Modeling Domain

The modeling domain utilized in this application re presented a sub-set of the inter-RPO’s continental modeling domain that covered the entire 48-state region with emphasis on the Ozone Transport Region. The OTC modeling domain at 12km horizontal mesh is displayed in Figure 1 is part of the 36km continent al domain that is designed to provide boundary conditions (BCs). The particulars of the t wo modeling domains are:

The 36km domain covered the continental US by a 14 9 by 129 mesh in the east-west and north-south directions, respectively. The domain is based on Lambert Conformal Projection with the center at (97oW 40oN) and paral lels at 33oN and 45oN. As evident from Figure 1, the 12km domain utilized in this ana lysis covers most areas of the eastern US and has 172 by 172 mesh in the horizontal. Both domains utilize 22 layers in the vertical extending to about 16km with 16 layers pla ced within the lower 3km.

Photochemical Modeling -- CMAQ

The CMAQ (version 4.5.1) with CB4 chemistry, aeroso l module for PM 2.5 and RADM cloud scheme was utilized in this study. Photochemi cal modeling was performed with the CCTM software that is part of the CMAQ modeling pac kage. Version 4.5.1 of this modeling software was obtained from the CMAS modeli ng center at http://www.cmascenter.org . The following module options were used in compili ng the
CCTM executable:

• Horizontal advection: yamo

• Vertical advection: yamo

• Horizontal diffusion: multiscale

• Vertical diffusion: eddy

• Plume-in-Grid: non operational

• Gas phase chemical mechanism: CB-4

• Chemical solver: EBI

• Aerosol module: aero3

• Process analysis: non operational

The following computational choices were made durin g compilation:

• Compiler version: PGI 6.0

• Fortran compiler flags:-Mfixed -Mextend -Bstatic -O 2 -module ${MODLOC} -I.

• C compiler flags : -v -O2 -I${MPICH}/include

• IOAPI library: version 3.0

• NETCDF library: version 3.6.0

• Parallel processing library version: mpich 1.2.6

• Static compilation on 32-bit system The following choices were made for running the exe cutable:3

• Number of processors: 8

• Domain decomposition for parallel processing: 4 col umns, 2 rows

• Number of species written to the layer-1 hourly-ave rage concentration output (ACONC) file: 39 (O3, NO, CO, NO2, HNO3, N2O5, HONO , PNA, PAN, NTR, NH3, SO2, FORM, ALD2, PAR, OLE, ETH, TOL, XYL, ISOP , ASO4I, ASO4J, ANO3I, ANO3J, ANH4I, ANH4J, AORGAI, AORGAJ, AORGPAI , AORGPAJ, AORGBI, AORGBJ, AECI, AECJ, A25I, A25J, ACORS, ASEA S, ASOIL)

• Each daily simulation was performed for 24 hours st arting at 05:00 GMT (00:00 EST)

The following postprocessing steps were performed u sing utility tools from the “ioapi” software package obtained from http://www.baronams.com/products/ioapi/AA.html#tool s :

• Extract and combine the following species for each hour for the first 16 model layers from the full 3-D instantaneous concentratio n output file: O3, CO, NO, NO2, NOY_1 ( =NO + NO2 + PAN + HNO3 ), NOY_2 ( =NO + NO2 + PAN +
HNO3 + HONO + N2O5 + NO3 + PNA + NTR ), HOX ( =OH + HO2 ), VOC ( =2*ALD2 + 2*ETH + FORM + 5*ISOP + 2*OLE + PAR + 7*T OL + 8*XYL ), ISOP, PM2.5 ( =ASO4I + ASO4J + ANO3I + ANO3J + ANH4I + ANH4J +
AORGAI + AORGAJ + 1.167*AORGPAI + 1.167*AORGPAJ + A ORGBI +
AORGBJ + AECI + AECJ + A25I + A25J ), PM_SULF ( =ASO4I + ASO4J ), PM_NITR ( =ANO3I + ANO3J ), PM_AMM ( =ANH4I + ANH4J ), PM_ORG_SA ( =AORGAI + AORGAJ

), PM_ORG_PA ( =1.167*AORGPAI +
1.167*AORGPAJ

), PM_ORG_SB( =AORGBI + AORGBJ ), PM_ORG_TOT ( =AORGAI + AORGAJ + 1.167*AORGPAI + 1.167*AORGPAJ + AORGBI +

AORGBJ ), PM_EC ( =AECI + AECJ ), PM_OTH ( =A25I + A25J ), PM_COARS ( =ACORS + ASEAS + ASOIL ), SO2, HNO3, NH3, H2O2

• Extract all species for all model layers for the la st hour of each daily instantaneous concentration output file to enable “ hot” restarts of modeling
simulations

• Create daily files of hourly running-average 8-hr o zone concentrations with time stamps assigned to the first hour of the averaging interval

The following files are archived on LTO2 computer t apes (each tape holds approximately 200 Gb of data) for each day:

• Aerosol/visibility file

• Layer-1 hourly-average concentration output file (c ontains 39 species)

• Dry deposition file

• Wet deposition file

• Extracted 16-layer species file

• Restart file (last hour of full 3-D instantaneous c oncentration file)

• Hourly 8-hr concentration file4

Photolysis Rates

One of the inputs to CMAQ is the photolysis rates. In this study, photolysis rate lookup tables were generated for each day of 2002 with the JPROC software that is part of the CMAQ modeling package. This software was obtained f rom the CMAS modeling center

at http://www.cmascenter.org . Rather than using climatological ozone column dat a, daily ozone column measurements from the NASA Earthprobe TOMS instrument were
downloaded from ftp://toms.gsfc.nasa.gov/pub/eptoms/data/ozone/Y200 2/ and used as input to the JPROC processor. It should be noted th at TOMS data were missing for the time period from August 3 – 11, 2002. The missing p eriod was filled as follows-- TOMS data file for August 2 was used as JPROC input for August 3 rd through August 7 th , and the TOMS data file for August 12 th was used as JPROC input for August 8 th through August 11 th

.

Boundary Conditions (BCs)

The boundary conditions for the 12km grid were extr acted from the 36km CMAQ simulation. The 36km simulation utilized boundary c onditions that were based on a oneway nest approach to GEOS-CHEM global model outputs (Moon and Byun 2004, Baker 2005). As stated above, the intent of the 36km CMA Q simulation was to provide the
BCs for the 12km model that would be more reflectiv e of the emissions and meteorology rather than to use either clean or arbitrary pollut ant fields. Also, in this study the CMAQ simulations utilized a 15-day ramp-up period, there by minimizing the propagation of the boundary fields into the areas of concern. A report on the setup and application of the 36km CMAQ and the extraction of the BCs is availabl e from NYSDEC.

Meteorological data

The meteorological data for this study was based on MM5 modeling (see Meteorological Modeling, 2007). The MM5 fields are then processed by MCIP version 3.0, a utility available as part of the CCTM software from CMAS Mo deling Center (see http://www.cmascenter.org ) to provide CMAQ model-ready inputs.

Emissions

The emissions data for 2002 were generated by indiv idual states within the OTR and were assembled and processed through the Mid Atlant ic Northeast Visibility Union (MANE-VU), a Regional Planning Organization (RPO). These emissions were then processed by NYSDEC using SMOKE processor to provid e CMAQ compatible inputs (Anthro-Emissions 2006). The 2002 emissions for the non-OTR areas within the modeling domain were obtained from the correspondin g RPOs and were processed using SMOKE, in a manner similar to that of the OTR.emiss ions. Details of this processing are outlined in the report (Pechan 2007), and the hourl y biogenic emissions (Bio-Emissions, 2006)5

CMAQ simulations

CMAQ simulations were performed using the one-way n esting approach in which we perform the continental CMAQ simulation at 36km gri d spacing. For this simulation we utilized clean initial conditions with boundary con ditions extracted from the simulation of GEOS-CHEM global chemical model. The interface prog ram used in this application was developed by University of Huston (Moon and Byun 20 04), which was applied to obtain hourly 36km boundary concentrations from GEOS-CHEM outputs. The CMAQ 36km simulation was initiated from December 15, 2001 wit h the first 15 days as spin up period and terminated on December 31, 2002. The simulation utilized the 2002 emissions data available from the RPOs and 2002 MM5 meteorological fields developed by the University of Maryland (TSD-1a). The hourly boundar y fields for the 12km CMAQ domain were obtained by application of BCON program to the 3-D concentration fields generated by the 36km CMAQ simulation.

The 12km simulations for both base and future year were assigned the boundary conditions based on the 36km CMAQ simulation and cl ean initial conditions. The simulation period covered was from April 15 through September 30, with the first 15 days of April set as ramp-up or spin-up period and that only data from May 1 through September 30 were used in the analysis. Details on CMAQ setup and run scripts are
available from NYSDEC.

References

Baker, K.: (2005) http://www.ladco.org/tech/photo/present/ozone.pdf

Moon, N. and D. Byun: (2004) A simple user’s guide for “geos2cmaq”code: Linking CMAQ with GEOS-CHEM. Version 1.0. Institute for Mul tidimensional Air quality Studies (IMAQS), University of Houston, Houston TX.

Meteorological Modeling: (2007) Meteorological Mode ling using Penn State/NCAR 5 th
Generation Mesoscale Model (MM5). TSD-1a

Pechan: (2006) Technical Support document for 2002 MANE-VU SIP Modeling inventories, version 3. Prepared by E. H. Pechan & Associates, Inc. 3622 Lyckan Parkway, Suite 2005, Durham, NC 27707.

Bio-Emissions: (2006) Processing of Biogenic Emissi ons for OTC/MANE-VU Modeling. TSD-1b

Anthro-Emissions: (2006) Emission Processing for th e Revised 2002 OTC Regional and
Urban 12 km Base Case Simulations. TSD-1c