••

•
•

•
•

•
•••

•
•••

•
••

•
•
••••

•
•

•
•
•
••••••••

•
•

•
•
•
••

• 

CALPUFF BART MODELING PROTOCOL • REFINED ANALYSIS 

LAFARGE BUILDING MATERIALS INC. 

TULSA CEMENT PLANT 

Prepared by: 

Vern Choquette • Principal Consultant 

Don Whitney, PE • Managing Consultant 

Eugene Chen, PE • Senior Consultant 

Jeremy Townley • Consultant 

TRINITY CONSULTANTS 

120 East Sheridan 

Suite 205 

Oklahoma City, OK 73104 

REeEl\J E0 

NOV n2 701'1-\ .I 

A\R QlJJ\LlTY -­

(405) 228-3292 

November 2, 2006 

Project 063701.0044 

n. TritlrityA
 '-LJnsu tants 

•

•

•

• • TABLE OF CONTENTS •

•

•

•

•

•

•

•

•

•

•

•

•
•

•

•

•
•

•

•

•

•

•

•

•

•

•

•
•

•

•
•

•

•

•

•

• Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 •

• 

1. EXECUTIVE SUMMARY •••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• 1-1 


1.1 BACKGROUND.................................•.......•.....•....•.....................•............................ 1-1 


1.2 OBJECTIVE ............................................................................................................ 1-1 


1.3 LOCATION OF SOURCE AND RELEVANT CLASS I AREAS ....................................... 1-2 


2. CALPUFF MODELING SYSTEM .•••••••••..••.•....••.•...•.•••..••.•...•••.•.•.••..••.••••.•••••.•..••. 2-1 


2.1 MODEL VERSIONS ............ .................................................................................... 2-1 


2.2 MODELING DOMAIN ............................................................................................. 2-2 


3. CALMET ..••••••••••••••••••••••••••••••.••.•••...••••••.••••.•••••••••••••.••••••••••••••••..•••••••••••••.•••••••• 3-1 


3.1 GEOPHYSICAL DATA............................................................................................. 3-1 


3.1.1 TERRAIN DATA................................................................................................... 3-1 


3.1.2 LAND USE DATA ................................................................................................ 3-1 


3.2 METEOROLOGICAL DATA ..................................................................................... 3-1 


3.2.1 MESOSCALE MODEL DATA ..........................................................................;..... 3-1 


3.2.2 SURFACE DATA .................................................................................................. 3-1 


3.2.3 UPPER AIR DATA................................................................................................ 3-2 


3.2.4 PRECIPITATION DATA......................................................................................... 3-2 


3.3 CALMET PROCESSING CONTROL ........................................................................ 3-3 


4. CALPlTFF ........................................................................................................... 4-1 


4.1 BART-ELIGIBLE SOURCES ................................................................................... 4-1 


4.2 SOURCE EMISSIONS DATA .................................................................................... 4-1 


4.2.1 PMIO EMISSIONS ................................................................................................. 4-2 


4.2.2 STACK PARAMETERS .......................................................................................... 4-3 


4.3 RECEPTOR LOCATIONS.......................................................................................... 4-3 


4.4 BACKGROUND OZONE .......................................................................................... 4-6 


4.5 BACKGROUND AMMONIA ..................................................................................... 4-6 


4.6 CALPUFF MODEL CONTROL P ARAMETERS ........................................................ 4-7 


5. CALPOST........................................................................................................... 5-1 


5.1 LIGHT EXTINCTION ALGORI11IM........................................................................... 5-1 


5.2 CALPOST PROCESSING MElHOD......................................................................... 5-1 


5.3 NATURAL BACKGROUND ...................................................................................... 5-2 


5.4 EVALUATING BART-EXEMPTION......................................................................... 5-2 


5.5 SUMMARY OF CALPOST CONTROL PARAMETERS............................................... 5-2 


APPENDIX A - EMISSION DATA 

APPENDIX B - DATA FILE LIST 
• • • • 

•

• 


LIST OF TABLES 


TABLE 2-1. CALPUFF MODELING SYSTEM VERSION NUMBERS ........................................................ 2-1 


TABLE 2-2. VERTICAL GRID STRUCTURE .......................................................................................... 2-3 


•

•

•

•

•

•

•

• TABLE 3-1. SUMMARY oFCALMET INPUTS .................................................................................... 3-4 


• TABLE 4-1. CONSTRUCTION DATE SUMMARy................................................................................... 4-1


• TABLE 4-2. ESTIMATED 24-HOUR AVERAGE EMISSION RATES ......................................................... 4-2 


TABLE 4-3. SPECIATED PM to EMISSIONS .......................................................................................... 4-3 


TABLE 4-4. POINT SOURCE STACK PARAMETERS.............................................................................. 4-3 


•

•

•

•

•

• 


TABLE 4-4. SUMMARY OF CALPUFF INPUTS ................................................................................... 4-7 


TABLE 5-1. MONTHLY HUMIDITY FACTORS " ..•...;.;.;........;.;....;.....;..,....:,.......~; ..:.;....;........;;;;;....:,::.... 5-1
 •

•

• 


TABLE 5-2. DEFAULT AVERAGE ANNUAL NATURAL BACKGROUND LEVELS .................................. 5-2 


TABLE 5-3. SUMMARY oFCALPOST INPUTS ................................................................................... 5-2
 •

•

• 

•

•
•

•

•
•

•

•

•

•

•

• 

Lafarge Building Materials 11 Trinity Consultants 

Refmed BART Modeling Protocol 063701.0044 

•

•

•

• 
•

•

•

•

• 

LIST OF FIGURES 

• FIGURE 1-1. LOCATION OF TuLSA PLANT AND NEARBY CLASS I AREAS ......................................... 1-3 


•

•

• 


FIGURE 2-1. PROPOSED MODELING DOMAIN .................................................................................... 2-2
 •

•

• FIGURE 3-1. PLOT OF MODELING DOMAIN USGS TERRAIN DATA ................................................... 3-1 


FIGURE 3-2. PLOT OF SURFACE STATIONS......................................................................................... 3-2
 •

•

• FIGURE 3-3. PLOT OF PRECIPITATION STATIONS ............................................................................... 3-3 


• FIGURE 4-1. CANEY CREEK WILDERNESS RECEPTOR LOCAnONS.................................................... 4-4 


FIGURE 4-2. HERCULES-GLADES WILDERNESS RECEPTOR LOCATIONS ........................................... 4-4 


• 


•

•

• FIGURE 4-3. UPPER BUFFALO WILDERNESS RECEPTOR LOCAnONS ................................................ 4-5 


FIGURE 4-4. WICHITA MOUNTAINS WILDERNESS RECEPTOR LOCAnONS ........................................ 4-6
 •
•

•

•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

• 

•

•

• Lafarge Building Materials ill Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

• 
• • • • • • • 

• • • • 

• • • • • 

1. EXECUTIVE SUMMARY 

•

•

• 

1.1 BACKGROUND 

On July 1, 1999, the U.S. Environmental Protection Agency (EPA) published the fmal Regional Haze •

•

• 

Rule (RHR). The objective of the RHR is to improve visibility in 156 specific areas across the United 

States, known as Class I areas. The Clean Air Act defines Class I areas as certain national parks (over 

6,000 acres), wilderness areas (over 5,000 acres), national memorial parks (over 5,000 acres), and 

international parks that were in existence on August 7, 1977. 

On July 6, 2005, the EPA published amendments to its 1999 RHR, often called the BART rule, which 

included guidance for making source-specific Best Available Retrofit Technology (BART) 

determinations. The BART rule defmes BART-eligible sources as sources that meet the following 

criteria: 

•

•

•

•

•

• (1) Have potential emissions ofat least 250 tons per year of a visibility-impairing pollutant,

•

• 

(2) Were "in existence" between August 7, 1962 and August 7, 1977, and 

(3) Are listed as one ofthe 26 listed source categories in the guidance. 

•

•

• 

A BART-eligible source is not automatically subject to BART. Rather, BART-eligible sources are 

subject-to-BART if the sources are "reasonably anticipated to cause or contribute to visibility 

impairment in any federal mandatory Class I area." EPA has determined that sources are reasonably 

anticipated to cause or contribute to visibility impairment if the visibility impacts from a source are 

greater than 0.5 deciviews (dv) when compared against a natural background. •

•

• Air quaJity modeling is the tool that is used to determine a source's visibility impacts. States have the 

authority to exempt certain BART-eligible sources from installing BART controls ifthe results of the 

dispersion modeling demonstrate that the source cannot reasonably be anticipated to cause or •

contribute to visibility impairment in a Class I area. Further, states also have the authority to define •

•

• 


the modeling procedures for conducting modeling related to making BART determinations. 


To promote consistency between states in the development ofBART modeling protocols and to •

• 

harmonize the approaches between adjacent RPOs, the Central States Regional Air Planning 

(CENRAP) organization developed BARTModeling Guidelines (December 15, 2005). The intent of 

the guidelines is to assist CENRAP states and source operators in the development of statewide and 

source-specific modeling protocols. 

·1.2 OBJECTIVE 

•

• 

The objective ofthis document is to provide a report summarizing the modeling methods and 

procedures that will be followed by Trinity when conducting refined modeling for Lafarge Building 

Materials' Inc. (Lafarge) Portland Cement Plant in Tulsa, OK. 

1-1
 Lafarge Building Materials Trinity Consultants, 


Refmed BART Modeling Protocol 063701.0044 

•

•
• 
 •

1.3 LOCATION OF SOURCE AND RELEVANT CLASS I AREAS
 •

•

• Lafarge's Tulsa Plant is located due east ofTulsa, OK and operates as a portland cement 

manufacturing plant. It operates two dry-process kilns that primarily consume a combination of coal 

and coke. Lafarge has identified the sources listed in Appendix A at the Tulsa Plant as sources that 

meet the three criteria for being BART-eligible sources. These sources include: 

.. Unit 221, Raw Material Grinding #2 


.. Unit 222, Auxiliary Raw Material Grinding #2 


.. Unit 223, Auxiliary Raw Material Grinding #2 


.. Unit 404, Kiln #2 Feed System 


.. Unit 502, Kiln #2 


.. Unit 706, Finish Grinding #2 


.. Unit 705, Finish Grinding #2 System 


.. Unit 704, Auxiliary Finish Grinding #2 


.. Unit 602, Clinker Cooler #2 


.. Unit 1101, Coal Unloading from Railcar 


.. Unit 1102, Coal Silo (drag conveyor and bucket elevators) 


.. Unit 1103, Coal Silo (crusher) 


A detailed summary of BART-eligibility is included in Section 4.1. Figure 1-1 provides a map ofthe 

location of the Tulsa Plant in relation to its neighboring Class I areas. 

As described in CENRAP's BART Modeling Guidelines, Class I areas within 300 km ofthe plant 

will be included in this refilled analysis. The following Class I areas fell within this 300 km radius: 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• Class I Area Name Distance from source (km)

•

•

• 

Caney Creek Wilderness 247 

Hercules-Glades Wilderness 258 

Upper Buffalo Wilderness 213 

Wichita Mountains National Wildlife Refuge 300 •

•

•

•

•
•

•

•

•

•

•

•
•

•

•

• 

1-2
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

•••••••••••••••••••••••••••••••••••••••••••• 

FIGURE 1-1. LOCATION OF TULSA PLANT AND NEARBY CLASS I AREAS 

Note: Tulsa Plant, Class I area, and domain boundary locations presented above are approximations for visual purposes and are not exact. 

1-3 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • 

• • • • • • 

2. CALPUFF MODELING SYSTEM 

•

•

•

•

•

• 

•

•

•

•

•

•

• The main components of the CALPUFF modeling system are CALMET, CALPUFF, and CALPOST. 

CALMET is the meteorological model that generates hourly three-dimensional meteorological fields 

such as wind and temperature. CALPUFF simulates the non-steady state transport, dispersion, and 

chemical transformation of air pollutants emitted from a source in "puffs". CALPUFF calculates 

hourly concentrations of visibility affecting pollutants at each specified receptor in a modeling 

domain. CALPOST is the post-processor for CALPUFF that computes visibility impacts from a 

source based on the visibility affecting pollutant concentrations that were output by CALPUFF. 

•

• 2.1 MODEL VERSIONS 

•

•

• 

Earth Tech, Inc. is the primary developer of the CALPUFF modeling system and all related programs. 

The versions of the CALPUFF modeling system programs that Trinity will use for Lafarge's 

modeling are listed in Table 2-1. Table 2-1 also compares the program versions that Trinity used 

with the program versions recommended by CENRAP. Note that some of the program versions are 

not the same as the program versions recommended by CENRAP. Certain programs will be utilized 

in order to avoid data compatibility issues with previous versions of the CALPUFF system. Trinity 

proposes the use of the following software version numbers listed below. Model selection was 

primarily driven by a desire to utilize the most recent release of CALMET. CALPlJFF and 

CALPOST executable versions were then selected as appropriate for compatibility. 

•

• 

TABLE 2-1. CALPUFF MODELING SYSTEM VERSION NUMBERS 

•

•

•

• 

Program CENRAP Suggested 

Lafarge Ref"med 

Analysis 

Version Level Version Level Reason for Difference

TERREL 3.311 030709 3.311 030709 

CTGCOMP 2.42 030709 2.22 030528 Version recommended is not 

available 

CTGPROC 2.42 030709 2.42 030709 

MAKEGEO 2.22 030709 2.22 030709 

CALMM5 2.4 050413 2.4 050413 

CALMET 5.53a 040716 6.211 060414 Most recent release of 

CALMET 

CALPUFF 5.711a 040716 6.112 060412 Use version compatible with 

CALMET v6.211 

POSTUTIL 1.4 040818 1.4 040818 

CALPOST 5.51 030709 6.131 060410 - Use version compatible with 

CALPUFF 6.112 •

•

•

•
•

•

•

•

• 

2-1
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • • • • • • • 

2.2 MODELING DOMAIN 

Trinity will utilize a modeling domain that extends at least 50 km in all directions beyond the Tulsa 

Plant and each ofthe four surrounding Class I areas. The map projection for the modeling domain is 

Lambert Conformal Conic (LCC) and the coordinate system will be North American Datum 1927 

(NAD27), which is an LCC projection. The meteorological grid spacing is 3 km. 

The southwest comer of the modeling domain will be Latitude 34.00'N, Longitude 100.5°W which is 

assigned as the 0, 0 reference point for the domain. The northeast comer of the modeling domain is 

approximately Latitude 37.25°N, Longitude 92°W. At a grid spacing of 3 km, the number of X grid 

cells is 256 and the number ofY grid cells will be 155. 

Figure 2-1 provides a plot ofthe modeling domain with respect to the sources and Class I areas. 

500 

450 

400 

~ 350 

! 

~ 300 

'c 

u" 

~ 250 

§" 

;; 200 Cll 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• .c 

E 

.::J 150 

100 

50 

o 

FIGURE 2-1. PROPOSED MODELING DOMAIN 

r------------------------------­

Hercules-GIlldes • 

Okiahom~City_OK 

NDrm~tOK 

WIChita_Mountaln.l' ED 


Lawton_OK 


~-----------------------

-50 0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 

Lam bert Conformal Conic (X) (km 1 

For the purposes ofthis analysis, the CALMET and CALPUFF doma41s will be and are singularly 

referred to herein as the "domain." 

2-2
 Lafarge Building Materials Trinity Consultants 


Refmed BART Modeling Protocol 063701.0044 

• • • • • • • • • • • • 

• 


The horizontal domain is comprised of grid cells, each containing a central grid point at which 

meteorological and computational parameters are calculated at each time step. For this analysis, grid 

spacing intervals of 3 Ian are selected to resolve terrain features within the domain. Table 2-2 

summarizes the vertical grid structure selected for both analyses. The cell face height of each cell 

indicates its vertical extent. The vertical domain is also composed of terrain-following grid cells, the 

number and size of which are chosen so as to constrain the boundary layer in which dispersion and 

chemical transformations take place. 

TABLE 2-2. VERTICAL GRID STRUCTURE 

Cell Face Height 

Vertical Grid Cell (meters) 

1 20 

2 40 

3 80 

4 160 

5 320 

6 640 

7 1,200 

8 2,000 

9 3,000 

10 4,000 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

2-3 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

•

•

• 
3. CALMET 


•

•

•

•

•

•

•

• Trinity will conduct a three-year CALMET analysis incorporating both mesoscale model and

•

•

•

• 

observation meteorological data to generate three years of CALMET outputs to be used as inputs to 

CALPUFF. The CALMET model requires the input of geophysical data, meteorological data, and 

model parameter settings. }he CALMET modeling procedures that Trinity will use generally follows 

the recommendations in CENRAP's protocol. It should be noted that Trinity will use surface and 

precipitation data but not upper air observations in this CALMET analysis. As a result, some of 

CENRAP's recommendations may not apply. 

3.1 GEOPHYSICAL DATA 

CALMET requires geophysical data to characterize the terrain and land use parameters that 

•

•

•

•

•

• potentially affect dispersion. Terrain features affect flows and create turbulence in the atmosphere 

and are potentially subjected to higher concentrations of elevated puffs. Different land uses exhibit 

•

•

• 

variable characteristics such as surface roughness, albedo, Bowen ratio, and leaf-area index that also 

effect turbulence and dispersion. 

3.1.1 TERRAIN DATA

•

• Trinity will use terrain data obtained from the United States Geological Survey (USGS) in 

I-degree (I :250,000 scale or approximately 90 meter resolution) digital format. A list of

•

• 

the USGS terrain files to be used in the analyses is provided in Appendix B. A plot of the 

land elevation for the modeling domain based on the referenced files is provided in Figure 

3-1. •

•

•
•

•

•
•

•

•

•

•

•

•

•

• 

• 

3-1 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Protocol 063701.0044 

•

•

• 
•••••••••••••••••••••••••••••••••••••••••••• 

FIGURE 3-1. PLOT OF MODELING DOMAIN USGS TERRAIN DATA 

tl 

~ ~ 35 II 

"d 

'i.....l 

34.5 

34 

315 

Longitude (degrees W) 

BOO 


700 


600 


500 


400 


300 


200 


100 


o 

Meters (Above MSL) 

3-1
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 06370l.0044 

• • • • 

• • • • • • 

•
•

•

•

•

•

• Trinity will use the TERREL program to generate grid-cell elevation averages from the 

USGS terrain data across the modeling domain.

•

• 3.1.2 LAND USE DATA 

Trinity will use USGS Composite Theme Grid (CTG) format Land Use and Land Cover •

• (LULC) data files at 1:250,000 resolutions. Where 1:250,000 land use data is not 

available, Trinity will use USGS data at 1: 100,000 resolutions. A list ofthe USGS land 

use files to be used is provided in Appendix B. 

• 

•

•

• Trinity will use the CTGPROC program to generate land use for each grid cell from the 

LULC data across the modeling domain. Since the USGS CTG format LULC data files 

must be compressed prior to use in the CTGPROC utility processor, Trinity will compress 

the files using the program CTGCOMP. 

3.2 METEOROLOGICAL DATA 

• CALMET will be used to assimilate data for 2001, 2002, and 2003 using mesoscale model output and 

National Weather Service (NWS) surface station observations, upper air station observations, and 

precipitation station observations to develop the meteorological field. 

3.2.1 MESOSCALE MODEL DATA 

• 


•

•

•

• Hourly mesoscale data will be used to supplement the hourly surface, upper air, and 


• 


precipitation observation data. The mesoscale data will be used to defme the initial guess 


field for the CALMET simulations. Trinity will use the following 5th generation mesoscale 


• 


model meteorological data (or MM5 data): 


• • 2001 MM5 data (12 km resolution), 

• 2002 MM5 data (36 km resolution), and • • 2003 MM5 data (36 km resolution) • 3.2.2 SURFACE DATA 

The use ofmultiple stations for meteorological observations in CALMET/CALPUFF· 

•

• 

provides a substantial enhancement over the steady-state treatment ofobservations from a 

single meteorological station. Parameters affecting turbulent dispersion that are observed 

hourly at surface stations include wind speed and direction, temperature, cloud cover and. 

ceiling, relative humidity, and precipitation type. Surface data will be extracted from 

among the surface stations listed in Appendix B and are illustrated in Figure 3-2. These 

sources were selected from the available data inventory to optimize spatial coverage and 

representation ofthe domain. Raw observations from these stations were obtained from 

the National Climatic Data Center (NCDC), quality assured and were merged using the 

•

•

• 

SMERGE pre-processor to create a single assimilated data file of surface observations for 

each year analyzed. 

•

•

• 

3-1 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 
•

•

• 


FIGURE 3-2. PLOT OF SURFACE STATIONS 

i<SLG

• KFW. .. 

KCNU 

•

•

•

• 	 500 •

I<I'PF KSGF r- - - - -KP2ii" -110 - - -KJUi - &- 1 ~ I<t:FV 

• i<l..BL ~ 

I<BVO 400 I<PNC • 	 Hercules-G 

•

lades 

E 	 KFI.P ~G KHRO 

•

•

•

•

•

• 

2:S 

i<SWO ~ 

u '(5 300 


u 


ro 
 KOKIfI1K I<CSM .. KOUN 0 • ~ 

.. 0 KH8R 200 u • I -e I<M1I1IT KHOT (]) Wichita_Mountainal I<LAW 	 .. .. .0 ..

•

•

• 	

E 

('iJ KFDR KADM 

•

•

•

• 

c 

E 

.. .. -l 

100 I<5PS .. KF3I I<I'RX. .. .. KTXK .. IELD.. L ___________ ~_~ _______________ l

. , , , , , I , , , , ! , , , , I , , , , I , ,', , I , , , , , , I , , , 

a 100 200 300 400 500 600 700 800 

•

•

• 

o 

Lambert Conformal Conic (Xl (km) •

•

• 3.2.3 UPPER AIR DATA 

Upper air observations will not be used in this CALMET analysis (i.e., NOOBS 1). •

•

•

• 

3.2.4 PRECIPITATION DATA 

The effects of chemical transformations and deposition processes on ambient pollutant •

• 

concentrations will be considered in this analysis. Therefore, it will be necessary to 

include observations of precipitation in the CALMET analysis. The precipitation stations 

from which Trinity will extract data are listed in Appendix B. The locations ofthe • precipitation stations with respect to the modeling domain are shown in Figure 3-3. These •

• 	

stations were selected from the available data inventory to optimize spatial coverage and 

representation ofthe domain. Data from the stations were processed for use in CALMET

•

• 	

using Version 5.31, Level 030528 of EPA's PMERGE program. 

•

•

•

•

• 	 3-2 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 •

• 
• • • • 

• • • • • 

•
•

•

• 	 FIGURE 3-3. PLOT OF PRECIPITATION STATIONS 

latA1 STO~\lA LE~O 500 

•

•

•

•

•

• 400 ........ 

E 

£ 

•

•

•

•

• .~ 

s:: 300' 0 

(.) 

Ci 

E 

.E 

'­

•

• s:: 

0 200 

• 

(.) 

1:: 

II) 

.c 

E

•

• m 

....J 

100 

o 

-e!ikUT- LOCt:;..,..,t - MNfSHF 

r,1'il.m'iING 

IND~ JO~UN 

• GRhT 

WlYNOK 

G~Y 

t:J:CI ~O

CHI~KA I 

VVichita Mountainsl 

I 

DUNCAN •

~T 

Ct..4KS 

FOREMII. • MILLIJIO

~O 

TEXARK • . LE\M •S\I 

I CHltoR 

I 

I 

L 

ARDttOFl 

cot4JdiHbtoONT ­ J,ocf<ss OelrON FRI~CO 

o 	 100 200 300 400 500 600 700 800 

Lambert Conformal Conic (X) (km) 

3.3 CALMET PROCESSING CONTROL 

•

•

•

•

•

• 

Table 3-1 provides a listing ofthe CALMET parameters that Trinity will use in Lafarge's modeling 

analysis. In addition to the parameters that will be used for Lafarge's modeling, the table also lists 

CENRAP's recommended parameters for comparison. In cases where a parameter to be used is 

different than what CENRAP recommended, a short explanation as to the difference is proved. Note 

that most ofthe differences from CENRAP's recommended parameters are due to the inclusion of 

observation data into the modeling analysis, since CENRAP's parameters are based on a noobservation analysis. 

•

•

•

•

•

•

• 	 3-3 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 •

• 
• • 

• • • • • • • • • • • 

•

•

•

• TABLE 3-1. SUMMARY OF CALMET INPUTS •

•

•
•

•
•

•

•

•
•

•
•

•

•
•

• 


•

•
•

• 


•

•

•

•

•

•

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

NUSTA Number of upper air 

data sites 

0 0 Only surface and 

precipitation 

observations will be 

used 

NOWSTA Number of overwater 

data sites 

0 0 

IBYR Starting year 2001 l~ate Years 2001, 2002, 

ar 2003

IBMO Starting month 1 Appropriate 

month 

Due to file size, 

analysis will be 

completed by month 

IBHR Starting hour 1 I 

IBTZ Base time zone 6 6 

IRLG Length of run 6 Varies with 

month 

Due to file size, 

analysis will be 

completed by month 

IRTYPE Run type (I for 

CALPUFF) 

I 1 

LCALGRD Compute CALFRID 

data fields (T run 

CALGRID) 

F F 

ITEST Stop run after SETUP 

to do input QA (2 

run) 

2 2 

PMAP Map proiection LCC LCC 

RLATO Latitude (decimal 

degrees) of projection 

origin 

40N 33.25N Appropriate for 

domain 

RLONO Longitude (decimal 

degrees) ofprojection 

origin 

97W 100.5W Appropriate for 

domain 

XLATl Latitude of 1st 

standard parallel 

33N 33N Appropriate for 

domain 

XLATI Latitude of 2nd 

standard parallel 

45N 40N Appropriate for 

domain 

DATUM Datum region for 

output coordinates 

WGS-G WGS-G Selected datum to 

match datum of 

land use data 

NX Number of X grid 

cells in 

meteorological grid 

300 256 Appropriate for 

domain 

3-4
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • 

• • • • • • • • • • • 

•
•

•

..


•

•

•
•

•

•

•

• 


•

•

•

•
•

•

•

• 


•

• 


CALMET Variable Desc 

Value Included:ln 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

NY Number ofY grid 

cells in 

meteorological grid 

192 155 Appropriate for 

domain 

DGRIDKM Grid spacing (lan) 6.0 3.0 Refined grid 

spacing 

XORIGKM Ref. coordinate of 

SW comer ofgrid 

cell 

-1008 ° Appropriate for 

domain 

YORIGKM Ref. coordinate of 

SW comer ofgrid 

cell 

0.0 ° Appropriate for 

domain 

NZ Number of vertical 

layers 

10 10 

ZFACE Vertical cell face 

heights (NZ + 1 

values) 

0,20,40,80, 160, 

320,640,1200,2000, 

3000,4000 

0, 20,40, 80, 

160,320,640, 

1200,2000, 

3000,4000 

LSAVE Save met. data fields 

in an unformatted 

file? 

T T 

IFORMO Type of unformatted 

output file (I for 

CALPUFF) 

1 1 

LPRINT Print met. fields F F 

IPRINF Print intervals 1 I 

IUVOUT(NZ) Specify layers ofu,v 

wind components to 

print 

NZ*O NZ*O 

IWOUT(NZ) Specify layers ofw 

wind component to 

print 

NZ*O NZ*O 

ITOUT(NZ) Specify layers of3D 

temperature field to 

print 

NZ*O NZ*O 

LDB Print met data and 

variables 

F F 

NNl First time step for 

debug data to be 

printed 

1 1 

NN2 Last time step for 

debug data to be 

printed 

1 2 Will generate debug 

data for a total of 2 

time steps 

IOUTD Control variable for 0 

• ° •

•

• 

3-5 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

• 
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

• 


Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

•

• 

• 

•

• 

CALMET Variable Description 

writing test/debug 

wind fields 

NZPRN2 Number of levels 

starting at surface to 

print 

IPRO Print interpol 

wind compo 

IPRI Print terrain adjusted 

surface wind 

components 

IPR2 Print initial 

divergence fields 

IPR3 Print fmal w' ~ ~ 

and direction 

IPR4 Print fmal di9 

fields 

IPR5 Print winds after I 

kinematic effects 

IPR6 Print winds after 

Froude number 

adjustment 

IPR7 Print winds after 

slope flows are added 

IPR8 Print fmal wind field 

components 

NOOBS No observation mode 

(2 = No surface, 

overwater, or upper 

air observations; use 

MM5 for surface, 

overwater, and upper 

air data) 

NSSTA Number of 

meteorological 

stations in 

SURF.DAT file 

NPSTA Number of 

precipitation stations 

in PRECIP.DAT file 

ICLOUD Gridded cloud fields 

(0 = no, 3 Gridded 

cloud cover from 

prognostic relative 

Value Included iIi 

CENRAP Protocol 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

3 

3-6 

Value Trinity 

Will Use 

I 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I 

47 

l35 

3 

Notes 

Default 

Will use surface and 

precipitation 

observations 

Number of stations 

Number of stations 

• 
• • • 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

•
•

•

•

• 


•

•
•

•

•

•

•

• 

CAL 

• 

scription 

Value Included in 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

humidity) 

IFORMS Format of surface 

data (2 :; formatted) 

2 2 

IFORMP Format of 

precipitation data (2 

formatted) 

2 2 

IFORMC Format of cloud data 

(2 formatted) 

2 I NtA - No cloud data 

used in model 

IWFCOD Generate winds by 

diagnostic wind 

module? (1 yes) 

1 1 

IFRAD] .Adjust winds using 

Froude number 

effects? (1 == yes) 

I 1 

IKINE Adjust winds using 

kinematic effects? (0 

= no) 

0 1 Will compute 

kinematic effects in 

this analysis 

IOBR Use O'Brien 

procedure for vertical 

winds? (0 = no) 

0 0 

rSLOPE Compute slope 

flows? (1 = yes) 

1 I 

IEXTRP Extrapolate surface 

winds to upper layers 

(-1 = no extrapolation 

and ignore layer 1 of 

upper air station data) 

-1 -1 

rCALM Extrapolate surface 

winds even if calm? 

(0 = no) 

0 0 

BIAS Layer dependent 

biases weighting aloft 

measurements 

0, 0, 0, 0, 0, 0, 0, 0, 0, 

° 

0,0,0,0,0,0, 

0,0,0,0 

RMIN2 Minimum vertical 

extrapolation distance 

Distance (lan) around 

an upper air site 

where vertical 

extrapolation is 

excluded (set to -1 if 

IEXTRP=±4) 

-1 -1 

IPROG 

• 

Using prognostic or 14 14 

3-7 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

•

•

•

e: 
•
•

•

•

•

•

•
•

•
•

•

•
•

•
•

•
•
•

•
•
•

•

•

•

•

•
•

•

•

•
•

•

•
•
•
•

•
•

•
•

•

•

•

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

MM-FDDA data? (14 

=Use winds from 

MM5.DAT as initial 

guess wfud field) 

ISTEPPG Timestep (hours) of 

the MM5 data 

I 1 

LVARY Use varying radius of 

influence to develop 

surface winds? 

T T Use FLAG default 

RMAXl Maximum radius of 

influence over land in 

surface layer (km) 

30 100 As included in 

IWAQM Phase 2 

case study 

RMAX2 Maximum radius of 

influence over land 

aloft (km) 

30 500 As included in 

IWAQM Phase 2 

case study 

RMAX3 Maximum radius of 

influence over water 

(km) 

50 36 Saine as MM5 data 

spacing 

RMIN Minimum radius of 

influence used 

anywhere (km) 

0.1 0.1 

TERRAD Radius of influence 

ofterrain features 

(km) 

12 12 

RI Weighting of first 

guess surface field 

(km) 

1 50 As included in 

IWAQM Phase 2 

case study 

R2 Weighting of first 

guess aloft field (km) 

1 100 As included in 

IWAQM phase 2 

case study 

RPROG MM5 windfield 

weighting parameter 

(k:tD.) 

0 0 

DlVLIM Maximum acceptable 

divergence 

5.E-6 5.E-6 

NITER Max number of 

passes in divergence 

minimization 

50 50 

NSMTH Number of passes 

through smoothing 

filter in each layer of 

CALMET(NZ 

values) 

2,4,4,4,4,4,4 2,4,4,4,4,4, 

4,4,4,4 

3-8
 Lafarge Building Materials Trinity Consultants 


Refmed BART Modeling Protocol 063701.0044 

•

•

•

•

•

•

•

•

•
•

•

•
•
•

•

•
•

•

•
•

•

•

•

•
•

•

•
•

•

•
•

•

•

•

•
•

•
•

•

•

•

•

•

• 

CALMET Variable Description 

Value Included iIi 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

NlTR2 Max number of 

stations used in each 

layer for the 

interpolation of data 

to a grid point (NZ 

values) 

5,5,5,5,5,5,5,5,5, 

5 

5,5,5, 5,5,5, 

5,5,5,5 

CRlTFM Critical Froude 

number 

1.0 1.0 

ALPHA Kinematic effects 

parameter 

0.1 0.1 

FEXTR2 Scaling factor for 

extrapolating surface 

winds aloft 

NZ*O.O NZ*O.O 

NBAR Number of terrain 

barriers 

0 0 

IDIOTPI Compute temperature 

from observations (0 

= true) 

0 0 

ISURFT Surface station to use 

for surface 

temperature (between 

1 and NSSTA) 

4 4 

IDIOPT2 Domain-averaged 

wind component 

switch 

0 0 

IUPT Station for lapse rates 

(between 1 and 

NUSTA) 

2 2 

ZUPT Depth through which 

lapse rate is 

calculated 

200 200 

IDIOPT3 Domain averaged 

wind component 

switch 

0 0 

IUPWND Upper air station for 

domain winds 

-1 -1 

ZUPWND Bottom and top of 

layer through which 

the domain scale 

winds are computed 

1., 1000. 1., 1000. 

IDIOPT4 Observed surface 

wind component 

switch 

0 0 

3-9
 Lafarge Building Materials Trinity Consultants 


Refilled BART Modeling Protocol 063701.0044 

• • • • • • •••••

• •••••••••••••

•
•••

•
••••••••

• 

Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

• •••

• 

CALMET Variable Description 

IDIOPT5 Observed aloft wind 

component switch 

LLBREZE Use lake breen > 

module? 

NBOX Number oflake 

breeze regions 

NLB Number ofstations in 

the region 

METBXID(NLB) Station IDs in the 

region 

CONSTB Neutral stability 

mixing height 

coefficient 

CONSTE Convective stability 

mixing height 

coefficient 

CONSTN Stable stability 

mixing height 

coefficient 

CONSTW Overwater mixing t=coefficient 

- FCORIOL lute value of 

olis parameter 

lAVEZI Conduct spatial 

averaging? (1 = yes) 

MNMDAV Max search radius in 

averaging process 

(number ofgrid cells) 

HAFANG Half-angle ofupwind 

looking cone for 

averaging (degrees) 

ILEVZI Layers ofwind use in 

upwind averaging 

(between 1 and NZ) 

DPTMIN Minimum potential 

temperature lapse rate 

in the stable layer 

above the current 

convective mixing 

height 

DZZI Depth oflayer above 

current convective 

mixing height 

Value Included in 

CENRAP Protocol 

0 

F 

0 

0 

0 

1.41 

0.15 

2400 

0.16 

I.E-4 

1 

10 

30 

1 

0.001 

200 

3-10 

Value Trinity 

Will Use 

0 

F 

0 

0 

0 

1.41 

0.15 

2400 

0.16 

1.E-4 

1 

10 

30 

1 

0.001 

200 

Notes 
• • • • • •••••••••••••••••••••••••••••

• ••••

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

through which lapse 

rate is computed (m) 

ZIMIN Minimum overland 

mixing height (m) 

50 50 

ZIMAX Maximum overland 

mixing height (m) 

3000 3000 

ZIMINW Minimum overwater 

mixing height (m) 

50 50 

ZIMAXW Maximum overwater 

missing height (m) 

3000 3000 

ITPROG 3D temperature from 

observations or from 

MM5? 

2 I Will use surface 

observations and 

precipitation data 

IRAD Type of interpolation 

(1 = lIr) 

I I 

TRADKM Temperature 

interpolation radius 

of influence (lan) 

36 36 

NUMTS Max number of 

stations for 

temperature 

interpolations 

5 5 

IAVET Spatially average 

temperature? (1 = 

yes) 

I 1 

TGDEFB Temperature gradient 

below the mixing 

height over water 

(KIm) 

-.0098 -0.0098 

TGDEFA Temperature gradient 

above the mixing 

height over water 

(KIm) 

-.0045 -0.0045 

JWATl Beginning land use 

categories over water 

55 55 

JWAT2 Ending land use 

categories for water 

55 55 

NFLAGP Precipitation 

interpolation flag (2 = 

1I?) 

2 2 

SIGMAP 

••

• 

Radius ofinfluence 50 50 

3-11 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

• 
•

•
•

•
••

•
•
•
•
•
•
•
•

•
•
•••••••

•
•••

•
•

•
••

•
•••••

•
••••

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Will Use Notes 

for precipitation 

interpolation (km) 

CUTP Minimum 

precipitation rate cut 

off(mmIhr) 

0.01 0.01 

3-12
 Lafarge Building Materials Trinity Consultants 


Refmed BART Modeling Protocol 063701.0044 

4. CALPUFF 


Trinity will conduct a three-year CALPUFF analysis. TheCALPUFF model requires the input of 

meteorological data output by CALMET, source emissions data, receptor data, ozone and ammonia 

data, and model parameter settings. 

4.1 BART-ELIGmLE SOURCES 

As a Portland Cement plant, the Tulsa Plant is considered BART-eligible. In order to determine what 

specific units constitute the "BART-eligible" source, the construction dates ofthe each of the plant's 

units must be examined. A list ofthose units that were "in existence" within the 15-year BART 

window (7 August 1962-1977) is provided below. 

TABLE 4-1. CONSTRUCTION DATE SUMMARY 

Unit No. PointID Description Construction Date 

221 EP04 #2 discharge chute, separator, bucket elev. 1964 


222 EP05 Raw mill pump hopper 1964 


223 EP05 Raw mill pump hopper 1964 


404 EP13 Kiln #2 feed system 1964 


502 EPt5 Kiln No.2, screw conveyors, & bucket elev. 1964 


706 EP24 #2 bucket elev. and separator 1964 


705 EP26 Finish mill #2 1964 


704 EP25 Finish mill #2 feed belt 1964 


602 EP17 Clinker cooler #2 1964 


1101 

R! 

EP27 CoaVcoak unloading 

EP28 CoaVcoak drag conveyor, bucket elev., crusher 

1974 


1974 


EP28 CoaVcoak drag conveyor, bucket elev., crusher 1974 


~:j\'\=l-l r-fJh oJ'?SOI 

4.2 SOURCE EMISSIONS DATA 

Trinity will include emissions of S02, NO", and PM lOin the model. For source-specific analyses, 

CENRAP has suggested several sources of emission data in the following priority:i 

.... Continuous Emission Monitoring (CEM) data, for period 2001-03 


.... CEM data for 2002 


.... Actual 2002 emission rates in state emission inventories 


.... Emission factors from AP-42 source profiles 


.... Permit allowable emissions, or 


.... 

Refined BART Modeling Protocol 063701.0044 


Potential to emit 


I CENRAP BART Modeling Guidelines. p. 7-/. 

4-1
 Lafarge Building Materials Trinity Consultants 

i

••
 •

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 4.2.1 PM10 EMISSIONS 

PMlO emissions are not directly entered into the model but are instead speciated out into its smaller 

constituents. The PMlO emissions from the kiln are speciated according to the Federal Land 

Managers' guidance to include the fol1owing: 2

,3 

• Coarse particulate matter (PMc, considered PM2•5- 1O) 

• Fine particulate matter (PMF, considered PM<2.5) 

• Sulfates (S04) 

• Secondary organic aerosols (SOA) 

• Elemental carbon (EC) 

Based upon FLM guidance, PMIO emissions from natural gas combustion consists of primarily coarse 

•

•

•

•

•

•

•

•

•

•

•

• 


(25%) and fine (75%) particulate matter and emissions from fabric filters consists ofprimarily coarse 


(70%) and fme (30%). A summary ofspeciated PMlO e!Uissions is presented in Table 4-3. Full 

details ofPMIO speciation can be found in Appendix A. •

•

• 


A detailed summary of emission rates and source data is included in Appendix A. Table 4-2 provides 

a brief summary ofthese emission rates. Emissions were calculated using emissions factors that were 

based on either AP-42 emission factors or baghouse manufacturer emission data. Based upon 

reported annual hours of operation from 2001-2005, it is reasonable to assume that each of these units 

operated continuously over a single 24-hour period. As a result, emissions are calculated based upon 

a full 24 hours of operation. In addition, emissions were based on the maximum achievable 

production rate (or heat duty in the case of combustion units) over a one hour period. 

TABLE 4-2. ESTIMATED 24-HOUR AVERAGE EMISSION RATES 

Unit ID NOx PM10 S02 

Ibs/hr Ibs/hr Ibs/hr 

EP04 0.089 3.857 0.001 

EP05 -­ 0.429 -­

EP13 -­ 0.557 -­

EP15 576.00 17.143 960.00 

EP24 -­ 3.402 -­ ···EP26 -­ 0.429 -­ EP25 -­ 0.429 -­ EPl7 -­ 4.607 -­ EP27 -­ 0.375 -­

EP28 -­ 0.321 -­

2 http://www2.nature.nps.gov/airlPermits/ect/ectCementKiln.cfin 


3 As recommended by Visibility Improvement State and Tribal Association of the Southeast (VISTAS). 


• 

4-2 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

•

•

• 

• 
TABLE 4-3. SPECIATED PM10 EMISSIONS 

•


•

•

•

•

•

•

•

•

•

•

•

•

•


•

•

•


•

•


•


•

•

•

•

•

•

•


•

•

•

•

• 


The base elevation of each the units is based upon visual inspection of USGS topographic maps. •

Source locations were based upon site inspections performed by Lafarge personnel. •

•

•

• 

4.3 RECEPTOR LOCATIONS 

The National Park Service (NPS) has electronic files for each Class I area available on their website 

•

•

• 

containing the locations and elevations of discrete Class I area receptors. The receptor files for Caney 

Creek Wilderness, Hercules-Glades Wilderness, Upper Buffalo Wilderness, and Wichita Mountains 

Wildlife Refuge will be downloaded from the NPS website, converted into the LCC NAD27 

•

•

•

• 

PM (coarse) PM (f"me) 804 80A EC 

Ibs/hr Ibs/hr Ibs/hr Ibs/hr Ibs/hr 

Unit ID 

EP04 0.964 2.893 -- -- -­

EP05 0.300 0.129 -- -­

EP13 0.390 0.167 -- -- -­ EP15 1.149 1.320 .. 12.91 1.766 -­

EP24 2.381 1.021 -- -- -­

EP26 0.300 0.129 -- -- -­

EP25 0.300 0.129 -- -- -­ EP17 3.225 1.382 -- -- -­

EP27 0.263 0.113 -- -- -­ EP28 0.225 0.096 -- -- -­

4.2.2 STACKPARAMETERS 

Table 4-4 provides a summary of source exhaust parameters. 

TABLE 4-4. POINT SOURCE STACK PARAMETERS 

Emission Latitnde Longitnde Base Exit Stack Stack Exit 

unit (degrees) (degrees) Elevation Velocity Height Diameter Temperatnre 

(m) (ft/s) (feet) (feet) (Fahrenheit) 

EP04 36.194777 95.811055 216 
 11.94 100 4 180 


EP05 36.19475 95.810722 216 
 8.29 55 1.6 70 


EP13 36.195472 95.811916 216 
 9.55 97 1.7 200 


EP15 36.195833 95.811777 216 
 17.83 60 6.9 500 


EP24 36.19425 95.811222 216 
 10.53 100 4 200 


EP26 36.194166 95.811166 216 
 33.16 64 0.8 70 


EP25 36.194166 95.811166 216 
 33.16 55 0.8 200 


EP17 36.19375 95.811805 216 
 3.56 53 8.0 200 


EP27 36.194027 95.811027 216 
 4.64 7 2.0 70


EP28 36.194138 95.811333 216 
 4.91 135 1.8 70 


4-3
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • 

• • • • • • • 

• • • • • • • • 

• • 

• • • 

• • • • • 

• • 

• • 

• • 

• • 

• 

• • 

• • 

• • 

• • 

• 

• • 

• 

• 

• • 

• • 

• • 

• • 

~~~.ili:~6rup AugU51 22,2003 

• .. 

.. • 

• 

• • 

• • 

• • 

•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• • •

•

•

•

•

•

•

•

• 

•

•

• 4-4 Lafarge Building Materials 

Refined BART Modeling Protocol •

• 

projection, and incorporated into the CALPUFF model. Receptor locations for these Class I areas are 

shown in Figures 4-1 to 4-4. 

FIGURE 4-1. CANEY CREEK WILDERNESS RECEPTOR LOCATIONS 

Caney Creek Wilderne~s 

80 Rt:ceptaTs 

1.J SFS Boundary Source: 

FS National Co.erag& tile: 

.. .. .. .. 

• • • • • • • • • .. • • 

• .. • .. • • • .. • • • • 

.. • .. • • • • 

~ • 

• 

• 

• 

.. 

• 

• 

FIGURE 4-2. HERCULES-GLADES WILDERNESS RECEPTOR LOCATIONS 

• 

• 

• 

Hercules-Glades Wilderness 

80 RSCIlplors 

USFS Boundary Source: 

FS National Coverage file: 

NRIS - ,ALp group 

Corvallis, OR 

Aug usl 22.2003 

Trinity Consultants 

063701.0044 
• • • • • • • • • • • • • • • • • • 

• • 

• • • • 

• • • • 

• • • 

• • • 

• • • 

• • • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• • 

• 

• • • 

• • • 

• 

• 

• 

•

•

• •

FIGURE 4-3. UPPER BUFFALO WILDERNESS RECEPTOR LOCATIONS • 

• •• • 

• • 

Uppe r Buffalo Wilderne ss 

72 Receptors 

USFS Boundary Source: 

FS National Coverage file: 

NRIS· ALP group 

Corvallis, OR 

August 22,2003 

• 

•

•

•

•

•

•

•

•

•

•
•

•

• 


• 

•

•
•

•

•

•

•

• 

4-5
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

" '", •

• 

FIGURE 4-4. WICHITA MOUNTAINS WILDERNESS RECEPTOR LOCATIONS 

•

•

•

•

• Wichita Mountains Wilderness ""thin Wichita Mountains NWR 

59 Receptors 

FWS Wilderness Boundary Source: 

• 

Region2GIS 

500 Gold A.B., SW 

PO 80.1306 

Albuquerque. Nt/! 87103 

FWS NWR Boundary Source: •

•

•

•

•

•

•

•

•

•

•

•

•

• •

4.4 BACKGROUND OZONE • 

• • •• 

• • • 

hltp :lIwww.fws.govldatalr2gislboundary.hlml 

SeplBm ber 1e, 2003 

• Background ozone concentrations are required in order to model the photochemical conversion of •

• 

S02 and NOli. to sulfates (S04) and nitrates (N03)' CALPUFF can use either a single background 

value representative of an area or hourly ozone data from one or more ozone monitoring stations.

•

•

•

• 

CENRAP recommends either developing background ozone estimates from ambient monitors located 

within the particular domain being modeled or developing background ozone estimates from 

CENRAP's most recent CMAQ or CAMx simulation for the 2002 base year. Trinity proposes the 

use of a conservative value of 40 ppb (as recommended in CENRAP protocol Appendix) in lieu of the 

CENRAP recommendation.

•

• 4.5 BACKGROUND AMMONIA 

• 

Background ammonia concentrations are required to model the formation of ammonium sulfates and •

• ammonium nitrates. CENRAP recommends developing background ammonia estimates from 

CENRAP's most recent CMAQ or CAMx simulation for the 2002 base year. Since CMAQ/CAMx 

modeled and observed monthly averaged ammonia concentrations exhibit wide spatial variability, 

•

•

• CENRAP recommends obtaining separate monthly-averaged ammonia concentrations from CMAQ 

or CAMx for the CENRAP north, central and south modeling domains, respectively. These would 

then be used as input to CALPUFF. Since the data from CENRAP's CMAQ and CAMx simulations 

is not readily available, Trinity proposes the use of a conservative monthly background concentration 

• 

4-6 Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

•

•

•

• 
..,


•

•

•

•

•

•

•

•

•

•

•

•

•
•
•

•
•

•

•

•

•
•

•

•

•

•

•

•

•

•

•

•
•

•

•
•

•
•

•
•

•

•

• 

of 3 ppb. This background concentration is the value included in CENRAP's protocol as a default 

background value for the CENRAP region. 

4.6 CALPUFF MODEL CONTROL PARAMETERS 

Table 4-4 provides a listing of the CALPUFF parameters that Trinity will use in Lafarge's modeling 

analysis. In addition to the parameters used for Lafarge's modeling, the table also lists CENRAP's 

recommended parameters for comparison. In cases where a parameter is different than what 

CENRAP recommended, a short explanation as to the difference is provided. 

TABLE 4-4. SUMMARY OF CALPUFF INpUTS 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

METRUN All model periods 

in met files will be 

run 

0 0 

IBYR Starting year 2001 Appropriate met year Years 2001, 2002, 

2003 

IBMO Starting month 1 =:i 

IBDY Starting day I 1 

IBHR Starting hour 1 1 

XBTZ Base time zone (6 

= CST) 

6 6 

IRLG Length of run 8760 8760 

NSPEC Number of 

MESOPUFF II 

chemical species 

10 10 

NSE Number of 

chemical species 

to be emitted 

8 6 S02, S04, NOx, 

PMC, PMF, SOA 

lTEST Program is 

executed after 

SETUP phase 

2 . 2 

MRESTART Do not read or 

write a restart file 

during run 

0 0 

NRESPD File written only at 

last period 

0 0 

METFM CALMETbinary 

file 

(CALMET.MET) 

I 1 

Averaging time in 

minutes 

60 60 

PGTIME PG Averaging 60 60 

4-7
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • 

•
•

•

•

•

•

•
•

•

•

•
•

•

•

•
•

•
•

•

•

•

•

•

•

•

•
•

•
•

•
•

•

•

•

•

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes. 

time in minutes 

MGAUSS Gaussian 

distribution used 

in near field 

I I 

MCTADJ Partial plUme path 

terrain adjustment 

3 3 

MCTSG Sub-grid-scale 

complex terrain 

not modeled 

0 0 

MSLUG Near-field puffs 

not modeled as 

elongated 

0 0 

MTRANS Transitional plume 

rise modeled 

1 1 

MTIP tip 

ash used 

1 1 

MSHEAR (0, 1) Vertical 

wind shear (not 

modeled, 

modeled) 

0 0 

MSPLIT Puffs are not split 0 1 Allowed for puff 

splitting due to 

significant sourcereceptor distances 

MCHEM MESOPUFFII 

chemical 

parameterization 

scheme 

1 1 

MAQCHEM Aqueous phase 

transfonnation not 

modeled 

0 0 

MWET Wet removal 

modeled 

I 1 

MDRY Dry deposition 

modeled 

I 1 

MDISP PG dispersion 

coefficients 

3 3 

MTURBVW Use both cry and 

crw from 

PROFILE.DAT to 

compute cry and crz 

(n/a) 

3 3 

MDISP2 PG dispersion 3 3

•

•

•

• 

4-8 Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

• 
• • • • • • • • • • 

• • • • • 

•

• 


•

•

•

•

•
•

•

•

•

• 


•

•

•

•

•

•

•

•

•
•

•

•
• •

Lafarge Building Materials Trinity Consultants •

•

• 

CALPUFF Value Included in 

Variable Description CENRAP Protocol 

coefficients 

MROUGH PG O'y and O'z not 0 

adjusted for 

roughness 

MPARTL No partial plume I 

penetration of 

elevated inversion 

MTINV Strength of 0 

temperature 

inversion 

computed from 

default gradients 

MPDF PDF not used for 0 

dispersion under 

convective 

conditions 

MSGTIBL Sub-grid TIBL 0 

module not used 

for shoreline 

MBCON Boundary 0 

concentration 

conditions not 

modeled 

MFOG Do not configure 0 

for FOG model 

output 

MREG Technical options 1 

must conform to 

USEPALong 

Range Transport 

(LRT) guidance 

CSPEC CENRAP 

Output 

Group 

Species Modeled Emitted 

S02 1 1 

S04 1 0 

NOX 1 1 

HN03 1 0 

N03 1 0 

NH3 1 1 

PMC 1 1 

PMF 1 1 

Dry 

Deposition 

1 

2 

1 

1 

2 

1 

2 

2 

4-9 

Value Trinity Will 

Use Notes 

0 

I 

0 

0 

0 

0 

0 

1 

Lafarge 

Output 

Group 

Species Modeled 

S02 1 

S04 1 

NOX 1 

HN03 1 

N03 1 

-

NH3 1 

PMC 1 

PMF 1 

Dry 

Emitted Deposition 

1 

1 

1 

0 

0 

0 

1 

1 

1 

2 

1 

I 

2 

1 

2 

2" 

Refined BART Modeling Protocol 063701.0044 
•

•
•
•

•

•

•

•

•
•

•

•
•
•

•
•

•
•

•

•
•

•
•

•

•
•

•
•
•
•

•
•
•

•

•

•
•
•
•
•

•

•

•

• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

EC I I 2 EC 1 0 2 

SOA I 1 2 SOA 1 1 2 

PMAP Map projection 

Number of X grid 

cells in 

meteorological 

grid 

UTM LCC 

NX 66 256 Based on refmed grid 

parameters as 

specified in 

CALMETruns 

NY Number of Y grid 

cells in 

meteorological 

grid 

66 155 Based on refmed grid 

parameters as 

specified in 

CALMETruns 

NZ Number of vertical 

layers in 

meteorological 

grid 

10 10 

DGRIDKM Grid spacing (km) 6 3.0 

0,20,40,80, 160, 

320,640, 1200, 

4000 

Refmed grid size 

ZFACE Cell face heights 

in meteorological 

grid(m) 

0, 20, 40, 80, 16O, 

320,640,1200,2000, 

3000,4000 

XORIGKM Reference X 

coordinate for SW 

comer of grid cell 

ofmeteorological 

grid (km) 

5 ° Appropriate for 

domain (refmed grid) 

YORIGKM Reference Y 

coordinate for SW 

comer of grid cell 

ofmeteorological 

grid (km) 

3327 ° Appropriate for 

domain (refmed grid) 

IUTMZN UTM zone of 

coordinates 

(NAD83) (use 

only if 

PMAP=UIM) 

12 14 N/A; LCC 

coordinates used 

IBCOMP X index oflower 

left comer of the 

computational grid 

I I 

JBCOMP Y index of lower 

left comer ofthe 

computational 

grids 

I I 

IECOMP X index of upper 

right comer of the 

66 256 Appropriate for 

domain (refined grid) 

4-10
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • 

• • • • • • • • • • • • • 

• • • • 

•

•

•
•

•

•

•

•
•

• 

•

•

•

•

•

•

•

•

• 

•

•

• 

CALPUFF 

Variable 

• 

Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

computational grid 

JECOMP Y index of upper 

right comer of the 

computational grid 

66 155 Appropriate for 

domain (refmed grid) 

LSAMP Sampling grid is 

not used 

F F 

IBSAMP X index of lower 

left comer of 

sampling grid 

1 1 

JBSAMP Y index of lower 

left comer of 

sampling grid 

1 1 

IESAMP X index ofupper 

right comer of 

samplingDid 

66 256 Appropriate for 

domain (refined grid) 

JESAMP Y index ofupper 

right comer of 

sampling grid 

66 155 ':Appropriate for 

domain (refmed grid) 

MESHDN Nesting factor of 

sampling grid 

1 1 

ICON Output file 

CONC.DAT 

containing 

concentrations is 

created 

1 1 

IDRY Output file 

DFLX.DAT 

containing dry 

fluxes is created 

1 1 

IWET 

. 

Output file 

WFLX.DAT 

containing wet 

fluxes is created 

1 1 

!VIS Output file 

containing relative 

humidity data is 

created 

1 1 

LCOMPRS Perform data 

compression in 

output file 

T T 

IMFLX Do not calculate 

mass fluzes across 

specific 

0 0 

4-11 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Protocol 063701.0044 
•
•

•

•

•
e 

e.

•

•

•
•

•

•

•
 e

•

•

•
e

•

•

e


•

•

•

•

•

•

•
•
•

•

•

•

•

•

•
•

• 


CALPUFF Value Included in 

Variable Description CENRAP Protocol 

boundaries 

IMBAL Mass balances for 0 

each species not 

reported hourly 

ICPRT Print concentration 1 

, fields to output list 

file 

IDPRT Do not print dry 0 

flux fields to 

output list file 

IWPRT Do not print wet 0 

flux fields to 

output list file 

ICFRQ Concentration 1 

fields are printed 

to output list file 

every hour 

IDFRQ Dry flux fields are 1 

printed to output 

list file every 1 

hour 

IWFRQ Wet flux fields are 1 

printed to output 

list file every 1 

hour 

IPTRU Units for line 3 

printer output are 

in glm3 for 

concentration and 

g/m2/s for 

deposition 

IMESG Messages tracking 2 

the progress of run 

written to screen 

Logical value for F 

debug output 

IP First uffto track 1 

NP puffs 1 

NNI Meteorological 1 

period to start 

output 

NN2 Meteorological 10 

Value Trinity Will 

Use 

0 

1 

0 

0 

I 

I 

1 

3 

2 

F 

1 

1 

1 

10 

Notes 

•

•

•

• 

4-12 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Protocol 063701.0044 

• 
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

• 


CALPUFF 

Variable 

• 

Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

period to end 

output 

NHILL Number of terrain 

features 

0 0 

NCTREC Number of special 

complex terrain 

receptors 

0 0 

MHILL Input terrain and 

receptor data for 

CTSG hills input 

in CTDM format 

2 2 

XHILL2M Conversion factor 

for changing 

horizontal 

dimensions to 

meters 

1 1 

ZHILL2M Conversion factor 

for changing 

vertical 

dimensions to 

meters 

1 1 

XCTDMKM X origin of CTDM 

system relative to 

CALPUFF 

coordinate system 

(Ian) 

0 0 

YCTDMKM Y origin ofCTDM 

system relative to 

CALPUFF 

coordinate system 

(Ian) 

0 0 

S02 Diffusivity 0.1509 0.1509 

Alpha star 1000 1000 

Reactivity 8 8 

Mesophyll 

resistance 

0 0 

Henry's Law 

coefficient 

0.04 0.04 

NOX Diffusivity 0.1656 0.1656 

Alpha star 1 1 

Reactivity 8 8 

•

•

• 

Mesophyll 

resistance 

5 5 

4-13 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

• 
• • • 

• • • • • • • • • • • • • • • 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 


•
•

•

•

.' 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

Henry's Law 

coefficient 

3.5 3.5 

HN03 Diffusivity 0.1628 0.1628 

Alpha star I I 

Reactivity 18 18 

Mesophyll 

resistance 

0 0 

Henry's Law 

coefficient 

8.e-8 8.e-8 

S04-2 Geomatric mass 

mean diameter of 

S04-2 (!-lm) 

0.48 0.48 

N03­ Geometric mass 

mean diameter of 

N03- (].lm) 

0.48 0.48 

PMC Geometric mass 

mean diameter of 

PMC (!-lm) 

6 6 

PMF Geometric mass 

mean diameter of 

PMF (!-lm) 

0.48 0.48 

EC Geometric mass 

mean diameter of 

EC (!-lID) 

0.48 0.48 

SOA Geometric mass 

mean diameter of 

SOA (!-lm) 

0.48 0.48 

RCUTR Reference cuticle 

resistance (s/cm) 

30 30 

RGR Reference ground 

resistance (s/cm) 

10 10 

REACTR Reference 

pollutant reactivity 

8 8 

NINT Number ofparticle 

size intervals for 

effective particle 

deposition velocity 

9 9 

IVEG Vegetation in nonirregated areas is 

active and 

unstressed 

I I 

S02 

•

• 

Scavenging 

coefficient for 

3.2IE-05 3.2IE-05 

4-14 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 
CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

liquid precipitation 

(S·I) 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

0 0 

S04-2 Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

I.OE-04 I.OE-04 

Scavenging 

coefficient for 

frozen 

precipitation (S'I) 

3.0E-05 3.0E-05 

RN03 Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

6.0E-05 6.0E-05 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

0 0 

N03­ Scavenging 

coefficient for 

liquid precipitation 

. (S·I) 

1.0E-04 I.OE-04 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

3.0E-05 3.0E-05 

NH3 Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

8.0E-05 NA 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

0 NA 

PMC Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

1.0E-4 1.0E-4 

Scavenging 3.0E-05 3.0E-05 

4-15
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

:.




•


,.• 


Ie.' 

~.

•

•

o 

.

e 

'- ' 

•

• 

•• 4-16

•

• 

' 

Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

., .... '. ""',"':.-' 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

coefficient for 

frozen 

precipitation (S-I) 

PMF Scavenging 

coefficient for 

liquid precipitation 

(S-I) 

1.OE-05 l.OE-05 

Scavenging 

coefficient for 

frozen 

precipitation (S-I) 

3.0E-05 3.0E-05 

EC Scavenging 

coefficient for 

liquid precipitation 

(S-I) 

1.OE-04 l.OE-04 

,Scavenging 

coefficient for 

frozen 

precipitation (s-I) 

3.0E-05 3.0E-05 

OC Scavenging 

coefficient for 

liquid precipitation 

(S-I) 

l.OE-04 l.OE-04 

" 

Scavenging 

coefficient for 

frozen 

precipitation (S-I) 

3.0E-05 3.0E-05 

MOZ Read ozone 

background 

concentrations 

from ozone.dat file 

(measured values) 

1 1 

BCK03 Background ozone 

concentration 

(ppb) 

12*40 12*40 

BCKNH3 Background 

ammonia 

concentration 

(ppb) 

12*3 12*3 

RNITE1 Nighttime S02 

loss rate is %/hour 

0.2 0.2 

RNITE2 Nighttime NOX 

loss rate is %/hour 

2 2 
•• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

RNITE3 Nighttime HN03 

loss rate is %lhour 

2 2 

MH202 Background H202 

concentrations 

1 ° Need to choose 

°in 

order to use monthly 

background value 

BCKH202 Background 

monthly H202 

concentrations 

1 12*1 

BCKPMF Fine particulate 

concentration for 

SOA option 

(llg/m3 ) 

1 1 

OFRAC Organic fraction of 

fme particulate for 

SOA option 

.2 0.15,0.15,0.2,0.2,0.2, 

0.2,0.2,0.2,0.2,0.2, 

0.2,0.15 

Irrelevant, since 

MCHEM not equal to 

4 

VCNX VOC/NOX ratio 

! for SOA option 

50 50 

SYDEP Horizontal size of 

a puff in meters 

beyond which the 

time dependant 

dispersion 

equation ofHeffter 

is used 

550 550 

MHFTSZ Do not use Heffier 

formulas for sigma 

z 

° 0 

JSUP Stability class used 

to determine 

dispersion rates for 

puffs above 

boundary layer 

5 5 

CONK] Vertical dispersion 

constant for stable 

conditions 

0.01 0.01 

CONK2 Vertical dispersion 

constant for 

neutral/stable 

conditions 

0.1 0.1 

TBD Use ISC transition 

point for 

determining the 

transition point 

0.5 0.5 

•

•

•

•

•

•

-.• 

4-17 Lafarge Building Materials Trinity Consultants • Refmed BART Modeling Protocol 063701.0044 • 
• • • • • • • • • • •••••••

•
••

•
••

•
••

•
••••••

• ••••

• 

CALPUFF 

Variable 

• 

Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

between the 

Schulman-Scire to 

Huber-Snyder 

Building 

Downwash 

scheme 

IURBI Lower range of 

land use categories 

for which urban 

dispersion is 

assumed 

10 10 

IURB2 Upper range of 

land use categories 

for which urban 

dispersion is 

assumed, 

19 19 

ILANDUlN Land use category 

for modeling 

domain 

* * 

XLAIIN Leaf area index for 

modeling domain 

* * 

ZOIN Roughness length 

in meters for 

modeling domain 

* * 

ELEVIN Elevation above 

sea level 

* * 

XLATIN North latitude of 

station in degrees 

-

XLONIN South latitude of 

station in degrees 

-

-

meter 

in meters 

I 10 10 

ISIGMAV Sigma-v is read 

for lateral 

turbulence data 

1 1 

IMIXCTDM Predicted mixing 

heights are used 

0 0 

XMXLEN Maximum length 

of emitted slug in 

meteorological 

grid units 

1 1 

XSAMLEN 

•

• 

Maximum travel 

distance ofslug or 

10 10 

4-18 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

• 
• • • • 

• • • • ••

•
•••••

•
•

• •••••••••••••••••••••••

• 

CALPUFF 

Variable 

• 

Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

puffin 

meteorological 

grid units during 

one sampling unit 

MXNEW Maximum number 

of puffs or slugs 

released from one 

source during one 

time step 

60 60 

MXSAM Maximum number 

ofsampling steps 

during one time 

step for a puffor 

slug 

60 60 

NCOUNT Number of 

iterations used 

when computing 

the transport wind 

for a sampling step 

that includes 

,transitional plume 

rise 

2 2 

SYMIN Minimum sigma y 

in meters for a 

new puff or slug 

I I 

SZMIN Minimum sigma z 

in meters for a 

new puff or slug 

I I 

SVMIN Minimum lateral 

turbulence 

velocities (mls) 

0.5 0.5 

SWMIN Minimum vertical 

turbulence 

velocities (mls) 

0.20,0.12, 0.08, 0.06, 

0.03,0.016 

0.20,0.12, 0.08, 

0.06,0.03,0.016 

WSCALM Minimum noncalm wind speeds 

(mls) 

0.5 0.5 

XMAXZI Maximum mixing 

height (m) 

3000 3000 

XMINZI Minimum mixing 

height (m) 

20 20 

SL2PF Maxmium cry/puff 

length 

10 10 

4-19 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 
~­

•

•

•

•

•

•

.,

•

•

•

•

• 

Lafarge Building Materials Trinity Consultants 

Refilled BART Modeling Protocol 063701.0044 

CALPUFF 

Variable Description 

PLXO Wind speed 

power-law 

exponents 

WSCAT Upper bounds of 

1st 5 wind speed 

classes 

PGGO Potential temp 

gradients PG E & N 

CDIV Divergence 

criterion for dw/dz 

(lis) 

PPC Plume path 

coefficients (only 

if MCTADJ = 3) 

NSPLIT Number of puffs 

when split 

IRESPLIT Hours when puff 

is eligible to split 

ZISPLIT Previous hours 

minimum mixing 

height,m 

ROLDMAX Previous max 

Olixing. 

height/current 

height ratio, must 

be less than this 

value to allow puff 

to split 

NSPLITH Number ofpuffs 

resulting from a 

split 

SYSPLITH Minimum sigma-y 

ofpuff before it 

may split 

SHSPLITH Minimum puff 

elongation rate 

from wind shear 

before puffmay 

split 

CNSPLITH Minimum species 

concentration 

before a pluff may 

Value Included in 

CENRAP Protocol 

0.07,0.07, 0.1 0,0.15, 

0.35,0.55 

1.54,3.09,5.14,8.23, 

10.80 

0.020, 0.035 

0.01 

0.5,0.5,0.5,0.5, 

0.35,0.35 

3 

1900 

100 

0.25 

5 

1.0 

2.0 

I.OE-07 

4-20 

Value Trinity Will 

Use 

0.07,0.07,0.10, 

0.15,0.35, 0.55 

1.54,3.09,5.14, 

8.23,10.80 

0.020, 0.035 

0.01 

0.5, 0.5, 0.5, 0.5, 

0.35,0.35 

3 

Hour 17 

100 

0.25 

5 

1.0 

2.0 

I.OE-07 

Notes 

Should be by hour of 

day-1900 is hour 17 

, 

., .~ .
CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity Will 

Use Notes 

split 

EPSSLUG Criterion for 

SLUG sampling 

1.0E-04 1.0E-04 

EPSAREA Criterion for area 

Trajectory step 

length for 

numerical site 

algorithm 

1.0E-06 1.0E-06 

DSRlSE 1.0 1.0 

NPTl Number of point 

sources with 

constant stack 

parameters or 

variable emission 

rate scale factors 

Varies by scenario 10 Appropriate for 

Lafarge 

IPTU Units for point 

source emission 

rates are gls 

1 3 Used different units 

(3 = lblhr) 

NSPTl Number of sourcespecies 

combinations with 

variable emissions 

scaling factors 

- -

,. 

NPT2 Number ofpoint 

sources with 

variable emission 

parameters 

provided in 

extemal file 

- -

MISC Other point source 

inputs include 

stack height, stack 

diameter, exit 

temperature, exit 

velocity, 

downwash flag 

and emissions by 

species 

- -

NARI Number of 

polygon area 

sources 

Varies by scenario 0 None modeled 

IARU Units for area 

source emission 

1 1 

4-21 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • • ••••••••

•
••••••••••••••••••••••••••

• 


CALPUFF 

Variable ~ 

rates are g/m2/s 

NSARI Number of source 

- species 

combinations with 

variable emissions 

scaling factors 

NAR2 Number of 

- buoyant polygon 

area sources with 

variable location 

and emission 

parameters 

NLN2 Number of 

- buoyant line 

sources with 

variable location 

and emission 

parameters 

NLINES Number of 

- buoyant line 

sources 

ILNU Units for line 

source emission 

rates in 't!s 

NSLNI Number ofsource­

- species 

combinations with 

variable emissions 

scaling factors 

MXNSEG Maximum number 

- of segments used 

to model each line 

NLRISE Number of 

-

distance at which 

transitional rise is 

computed 

XL Average line 

- source length (m) 

HBL Average height of 

-

line source height 

(m) 

WBL Average building 

- width(m) 

--

0 --------

Value Trinity Will 

Use Notes 

None modeled 

4-22
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

CALPUFF 

Variable Description 

WML Average line 

source width (m) 

DXL Average 

separation 

between buildings 

(m) 

FPRlMEL Average buoyancy 

parameter (m4/s3) 

NVLl Number of volume 

sources 

IVLU Units for volume 

source emission 

rates in 

grams/second 

NSVLl Number of sourcespecies 

combinations with 

variable emissions 

scaling factors 

IGRDVL Gridded volume 

source data is not 

used 

VEFFHT Effective height of 

emissions (m) 

VSIGYI Initial sigma-y 

value 

VSIGZI Initial sigma-z 

value 

NREC Number of nongridded receptors 

Value Included in 

CENRAP Protocol 

-

-

-

-

-

-

-

-

-

-

5630 

Value Trinity Will 

Use 

-

-

-

0 

-

-

-

-

-

-

271 

Notes 

None modeled 

-

Receptor data 

provided in 

Appendix 

..•••/. 

~. 

4-23 '. Lafarge Building Materials Trinity Consultants •

•

•

• 

Refined BART Modeling Protocol 063701.0044 
•
•

•
• 


5. CALPOST 


• 

•

•

•

•

• Trinity will conduct a three-year CALPOST analysis to determine the change in light extinction. 

caused by Lafarge's BART-eligible sources when compared to a natural background. The CALPOST 

model requires the input of concentration data output by CALPUFF. 

5.1 LIGHT EXTINCTION ALGORITHM 

Trinity will utilize EPA's currently approved algorithm for reconstructing light extinction (as opposed 

to the new equation for reconstructing light extinction recommended by the IMPROVE Steering 

Committee). The light extinction equation is provided below. 

bext = 3*f(RH)*[(NH4)2S04] + 3*f(RH)*[NH4N03] + 4*[OC] + l*[PMrJ 

+ 0.6 *[PMc ] + lO*[EC] + bRay 

•

•

•

•

•

•

•

•

•

• 


The algorithm will be used to calculate the daily light extinction attributable to Lafarge's BART­


eligible sources and light extinction attributable to a natural background. The change in deciviews 

based on the source and background light extinctions was evaluated using the equation below. 

/), dv = 10* In bext, background +bext, source 1 

[ bext, background 

•

•

•

•

• 5.2 CALPOST PROCESSING METHOD 

Trinity will use CALPOST Method 6, which calculates hourly light extinction impacts for the source •

• and background using monthly average relative humidity adjustment factors. Trinity will use 

monthly Class I area-specific relative humidity adjustment factors based on the centroid ofthe Class I 

area as included in Table A-3 of EPA's Guidance for Estimating Natural Visibility Conditions Under 

the Regional Haze Program. The factors for the relevant Class I areas in the analysis are provided in 

Table 5-1. 

TABLE 5-1. MONTHLY HUMIDITY FACTORS 

Class I Area Jan Feb Mar Apr May Jun Jut Aug Sep Oct Nov Dec 

Caney Creek 3.4 3.1 2.9 3.0 3.6 3.6 3.4 3.4 3.6 3.5 3.4 3.5 

Hercules-Glades 3.2 2.9 2.7 2.7 3.3 3.3 3.3 3.3 3.4 3.1 3.1 3.3 

Upper Buffalo 3.3 3.0 2.7 2.8 3.4 3.4 3.4 3.4 3.6 3.3 3.2 3.3 

Wichita Mountains 2.7 2.6 2.4 2.4 3.0 2.7 2.3 2.5 2.9 2.6 2.7 2.8 

•

•

•

•

•

•

•

•

•

•

•
•

•

•

• 

5-1 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

• 
•

•

•

•

•

•

•

•

• 


5.3 NATURALBACKGROUND 

For Caney Creek and the Wichita Mountains, Trinity will use default west natural background values. 

Hercules-Glades and Upper Buffalo will be processed using default east natural background values.4 

The annual average concentrations are included in Table 2-1 ofEPA's Guidancefor Estimating 

Natural Visibility Conditions Under the Regional Haze Program and are provided in Table 5-2. 

TABLE 5-2. DEFAULT AVERAGE ANNUAL NATURAL BACKGROUND LEVELS 

Component 

Average Annual 

Natural Background 

West (llg/m3 

) 

Average Annual 

Natural Background 

East (llg/m3 

) 

Ammonium Sulfate 0.12 0.23 

Ammonium Nitrate 0.10 0.10 

Organic Carbon Mass OA7 lAO 

Elemental Carbon 0.02 0.02 

Soil 0.50 0.50 

Coarse Mass 3.0 3.0 

5.4 EVALUATING BART-EXEMPTION 

Trinity will compare the 2001 through 2003 daily ~dv values output by CALPOST to a contribution 

threshold of 0.5 ~dv. If the 98th percentile ~dv values (the 7th highest value, 22nd highest over three 

years) output by CALPOST is less 0.5 ~dv, then the source is exempt from BART requirements. If 

the 98th percentile value of the ~dv values output by CALPOST is greater than 0.5 ~dv, then the 

source is considered subject to BART regulations and must start on a BART-determination process. 

5.5 SUMMARY OF CALPOST CONTROL,PARAMETERS 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

Table 5-3 provides a listing ofthe CALPOST parameters that Trinity will use in Lafarge's modeling 

analysis. In addition to the parameters used for Lafarge's modeling, the table also lists CENRAP's 

recommended parameters for comparison. In cases where a parameter to be used is different than 

what CENRAP recommended, a short explanation as to the difference is provided. 

•

• 

TABLE 5-3. SUMMARY OF CALPOST INPUTS 

CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity will 

use Notes 

ISYR Starting year 2001 Appropriate met 

year 

Years 200 I, 2002, 

2003 

ISMO Starting month I 1 

ISDY Starting day I 1 

4 Per comments received from EPA Region 7 to Trinity Consultants' Kansas City office, May 22, 2006. 

5-2 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

•

•

•

•

•

•

•

•

•

•

• 
• 

•
•

•

•

•

•
•

•

•

• 5-3 Lafarge Building Materials Trinity Consultants 

CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity will 

use Notes 

ISHR Starting hour 0 0 All CALPUFF 

periods will be 

included 

NPER Number ofperiods to 

process 

8760 8760 

NREP Process every hour of 

data? 1 =yes 

1 1 

ASPEC Process species for 

visibility 

VISIB VISIB 

ILAYER Layer/deposition code; 

1 for CALPUFF 

concentrations 

1 1 

A Scaling factor, slope 0 0 

B Scaling factor, 

intercept 

0 0 

LBACK Add hourly 

background 

concentrations of 

fluxes? 

F F 

LG Process gridded 

receptors? 

F F 

LD Process discrete 

receptors? 

T T 

LCT Process complex 

terrain receptors? 

F F 

LDRING Report receptor ring 

results? 

F F 

NDRECP Select all discrete 

receptors 

-1 Varies As appropriate for 

Class I area being 

analyzed 

IBGRID X index of LL comer 

of receptor grid 

-1 -1 

JBGRID Y index of LL comer 

ofreceptor grid 

-1 -1 

IEGRID X index ofUR comer 

ofreceptor grid 

-1 -1 

JEGRID Y index of UR corner 

ofreceptor grid 

-1 -1 

NGONOFF Number of gridded 

receptor rows 

0 0 

NGXRECP Exclude specific 

gridded receptors, Yes 

=0 

0 0 

• Refmed BART Modeling Protocol 063701.0044

•

• 

• • • • • • 

• • • • • 

•
•

•
•

•

•

•

•

•

• 


•

•
•

•
•

•

•

•

•
•

•
•

•

•
•

•
•

• 


•

•

•

•

• 

CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity will 

use Notes 

RHMAX Maximum RH% used 

in particle growth 

curve 

95 

• 

95 

LVS04 Compute light 

extinction for sulfate? 

T T 

LYN03 Compute light 

extinction for nitrate? 

T T 

LVOC Compute light 

extinction for organic 

carbon? 

T T 

LVMPC Compute light 

extinction for coarse 

particles? 

T T 

LVMPF Compute light 

extinction for fine 

particles? 

T T 

LVEC Compute light 

extinction for 

elemental carbon? 

T T 

LVBK Include background in 

extinction calculation? 

T T 

SPECPMC Coarse particulate 

species 

PMC PMC 

SPECPMF Fine particulate 

species 

Extinction efficiency 

for coarse particulates 

PMIO PMF Notation difference 

EEPMC 0.6 0.6 

EEPMF Extinction efficiency 

for fine particles? 

1.0 1.0 

EEPMCBCK Extinction efficiency 

for coarse part. 

Background 

0.6 0.6 

EES04 Extinction efficiency 

for ammonium sulfate 

3.0 3.0 

EEN03 Extinction efficiency 

for ammonium nitrate 

3.0 3.0 

EEOC Extinction efficiency 

for organic carbon 

4.0 4.0 

EESOIL Extinction efficiency 

for soil 

1.0 1.0 

, 

EEEC Extinction efficiency 

for elemental carbon 

10.0 I 10.0 

5-4
 Lafarge Building Materials Trinity Consultants 


Refmed BART Modeling Protocol 063701.0044 

-


CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity will 

use Notes 

MVISBK Method 6 for 

background light 

extinction 

6 6 

BEXTBTBK Background extinction 

for MVISBK=1 

12 Not Used Not necessary since 

MVISBK=6 

RHFRAC % ofparticles affected 

byRH 

10 Not Used Not necessary since 

MVISBK=6 

RHFAC Extinction coefficients 

for modeled and 

background 

hygroscopic species 

computed using EPA 

(2003) monthly RH 

adjustment factors 

Depends on Class I 

Area 

For Caney Creek: 

3.4,3.1,2.9,3.0, 

3.6,3.6,3.4,3.4, 

3.6,3.5,3.4,3.5 

Hercules-Glades: 

3.2,2.9,2.7,2.7, 

3.3,3.3,3.3,3.3, 

3.4,3.1, 3.1, 3.3 

Upper Buffalo: 

3.3,3.0,2.7,2.8, 

3.4,3.4,3.4,3.4, 

3.6,3.3,3.2,3.3 

Wichita Mountains: 

2.7, 2.6, 2.4, 2.4, 

3.0,2.7,2.3,2.5, 

2.9,2.6,2.7,2.8 

As appropriate for 

Class I area 

BKS04 Background sulfate 

extinction coeff­ east 

0.23 0.23 Used for HerculesGlades and Upper 

Buffalo 

BKN03 Background nitrate 

extinction coeff­ east 

0.10 0.10 Used for HerculesGlades and Upper 

I Buffalo 

BKPMC Background sulfate 

extinction coeff­ east 

3.00 3.00 Used for HerculesGlades and Upper 

Buffalo 

BKSOC Background organic 

carbon extinction 

coeff- east 

lAO 1.40 Used for HerculesGlades and Upper 

Buffalo 

BKSSOIL Background soil 

extinction coeff­ east 

0.50 0.50 Used for HerculesGlades and Upper 

Buffalo 

BKSEC Background elemental 

carbon extinction 

coeff- east 

0.02 0.Q2 Used for HerculesGlades and Upper 

Buffalo 

'.•


•

•

•

•

•

•

•

• 

5-5 • Lafarge Building Materials Trinity Consultants

Refined BART Modeling Protocol 063701.0044 

• 
•

•

•

•

•

•

•

•
•

•

•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

• 


• 

CALPOST 

Variable Description 

BKS04 Background sulfate 

extinction coeff- west 

BKN03 Background nitrate 

extinction coeff- west 

BKPMC Background sulfate 

extinction coeff- west 

BKSOC Background organic 

carbon extinction 

coeff- west 

BKSSOIL Background soil 

extinction coeff- west 

BKSEC Background elemental 

carbon extinction 

coeff-west 

BEXTRAY Extinction due to 

Rayleigh scattering 

ClIMm) 

LDOC Print documentation 

image? 

IPTRU Print output units for 

concentrations and for 

deposition 

LlHR Report 1 hr averaging 

times 

L3HR Report 3 hr averaging 

times 

L24HR Report 24 hr averaging 

times 

LRUNL Report run-length 

averaging times 

LT50 Top 50 table 

LTOPN Top N table 

NTOP Number of Top-N 

values at each receptor 

ITOP Ranks of Top-N 

values at each receptor 

LEXCD Threshold 

exceedances counts 

Value Included in 

CENRAP Protocol 

0.12 

0.10 

3.00 

0.47 

0.50 

0.02 

10.0 

F 

3 

F 

F 

T 

F 

F 

F 

4 

1,2,3,4 

F 

Value Trinity will 

use 

0.12 

0.10 

3.00 

0.47 

0.50 

0.02 

10.0 

F 

1 

F 

F 

T 

F 

F 

F 

4 

1,2,3,4 

F 

Notes 

Used for Wichita 

Mountains and 

Caney Creek 

Used for Wichita 

Mountains and 

Caney Creek 

Used for Wichita 

Mountains and 

Caney Creek 

Used for Wichita 

Mountains and 

Caney Creek 

Used for Wichita 

Mountains and 

Caney Creek 

Used for Wichita 

Mountains and 

Caney Creek 

Units preference 

•

•

• 

5-6 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•• 
-


CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity will 

use Notes 

THRESH I Averaging time 

threshold for I hr 

averages 

-I -I 

THRESH3 Averaging time 

threshold for 3 hr 

averages 

-I -I 

THRESH24 Averaging time 

threshold for 24 hr 

averages 

-1 -I 

THRESHN Averaging time 

threshold for NAVOhr averages 

-I -I 

NDAY Accumulation period, 

days 

0 0 

NCOUNT Number of 

exceedances allowed 

I I 

LECHO Echo option F F 

LTIME Time series option F F 

LPLT Plot file option F F 

LORD Use grid format 

instead of DATA 

format 

F F 

LDEBUO Output information for 

debugging? 

F F 

5-7 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 


.~. ". 
APPENDIX A - EMISSION DATA 


Lafarge Building Materials Trinity Consultants 


Refmed BART Modeling Protocol 063701.0044 

•••••••••••••••••••••••••••••••••••••••••••• 

Printed /013012006 

EUG 

Emission 

Point 

(Unit ID) 

Point 

ID 

Unit Description 

(Include MMBtu/hr 

if applicable) 

Which of the 

26 Categories 

does this Unit 

fall into?A 

Year 

Installed 

Was this Unit in 

Existence on 

August 7, 1977? 

(Y or N) 

Was this Unit in 

Operation after 

August 7, 1962? 

(Y or N) 

Potential to Emit (tons/year) 

PM 10 NO. S02 VOC NH3 

2A 

221 EP04 #2 discharJ!e chute, seperator, bucket elev. 3 1964 Y Y 28.90 2.70 0.02 0.20 

28 

222 

223 

EP05 

EP05 

Raw mill M pump hopper 

Raw mill #2 weiJ!h meters and feed belts 

3 

3 

1964 

1964 

Y 

Y 

Y 

Y 

3.20 

4C 

404 EP\3 Kiln #2feed system 3 1964 Y Y 4.20 

58 

502 EPI5 Kiln, screw conveyors, & bucket elevator 3 1964 Y Y 128.40 1261.00 2102.00 5.90 

78 

706 

705 

EP24 

EP26 

#2 bucket elevator and separator 

Finish mill #2 

3 

3 

1964 

1964 

Y 

Y 

Y 

Y 

25.50 

3.20 

7D 

704 EP25 Finish mill #2 feed belt 3 1964 Y Y 3.20 

6A 

602 EPI7 Clinker cooler #2 3 196411997 Y Y 34.50 

II 

1101 

1102 

1103 

EP27 

EP28 

EP28 

Coalkoak unloading 

Coalkoak drag conveyor, bucket elevators 

Coallcoak crusher 

3 

3 

3 

1974 

1974 

1974 

Y 

Y 

Y 

Y 

Y 

Y 

7.00 

2.40 ; 

Tolal 240.50 1263.70 2 \02.02 6.10 0 

A. See attached list of categones 

Lafarge Building Materials Using Title V Trinity Consultants 


BART Modeling A-1 063701.0044 

•••••••••••••••••••••••••••••••••••••••••••• 

Printed 1013012006 

Source Emissions Emission Factors 

Location Point III Unit Description Constr. Datel Heat Duty Process Rate Hours of Operation PM10 

(MMBtu/hr) TPH (hrs/day) (2c/acf) 

Tulsa EP05 Raw mill pump hopper 1964 -­ 100 24 0.025 

Tulsa EPI3 Kiln #2 feed system 1964 -­ 120 24 0.025 

Tulsa EP24 #2 bucket elevator andseparator 1964 -­ 39 24 0.025 

Tulsa EP26 Finish mil/1I2 1964 -­ 39 24 0.025 

Tulsa EP25 Finish mill #2feed belt 1964 -­ 39 24 0.025 

Tulsa EPI7 Clinker cooler 112 196411997 -­ 96 24 0.025 

Tulsa EP27 Coalkoak unloading 1974 -­ 150 24 0.025 

Tulsa EP28 Coalkoak drag conveyor, bucket elevators, crusher 1974 -­ 150 24 0.025 

Construction dates and source parameters as listed in Permit 98-174-TV. Construction ofEP14 (olher kiln) commenced on 1961 

which is outside the 7-Aug-62 to 7-Aug-77 BART eligibility window. 

Source Emissions Emission Data Input 

Location Point III Unit Description NO: PMlO PM (coarsei PM (finel SOl 

Ibslhr Ibs/hr Ibslhr Ibs/hr Ibslhr 

Tulsa EP05 Raw mill pump hopper 0.429 0.300 0.129 -­ Tulsa EP13 Kiln #2 feed system 0.557 0.390 0.167 -­ Tulsa EP24 #2 oucket elevator andseparator 3.402 2.381 1.021 -­ Tulsa EP26 Finivh mill #2 -­ 0.429 0.300 0.129 

Tulsa EP25 Finish mill #2feed belt 0.429 0.300 0.129 

Tulsa EPI7 Clinker cooler #2 -­ 4.607 3.225 1.382 

Tulsa EP27 Coalkoak unloading -­ 0.375 0.263 0.113 -­ Tulsa EP28 Coallcoak drag conveyor, bucket elevators, crusher -­ 0.321 0.225 0.096 ! 

2 PM Speciation based on AP 42 Table 8.2-2 category 4 (Dryers, grinders, etc.), PM2.5 is 30% ofTSP 

PMIO­ 2.

5 (filterable) 70% PM (coarse) 

PM2.5 (condenSible) 30% PM (fine) 

Lafarge Building Materials Fabric Filters Consultants 


BART Modeling A-2 063701.0044 

" 

Printed 1013012006 

Emission Factors for NO. & S~ based on AP-42 (1/95) Table 11.6-8, PM 10 based on baghouse manufacturer exhaust loading ratc. 

2 Construction dates and source parameters as listed in Pennit 98-1 74-TV. Unit EPI4 constructed outside of the 7-Aug-62 to 7-Aug-77 BART eligibility window 

3 Where necessary, Ib/MMSCF were converted to IblMMBtu using a natura! gas heating value (i 1020 BTU/SCF 

• Condensible/filterable PM speciation were based on the following profIles: 

National Park Service guidance for PC Dry Bottom FF. http://www2.nature.nps.gov/air/Pennits/ect/ectCemcntKiln.cfin 

Coarse 6.70% 

Fine Soil 7.70% Fine soil included in model as PM (fine) 

Fine EC 0.00% 

CPM lOR 75.30% Inorganic condensible PM considered sulfates pcr NPS guidance ".. CPM OR 10.30% Organic condensible PM considered seconpary organic aerosols per NPS 

'••~.~.:::I 

~ ., 

~: Lafarge Building Materials Kiln Trinity Consultants 

BART Modeling A-3 063701.0044 

~. 

~ 

~' 

l~ 

Source Emissions' Maximuin Hourh Emission Factors 

Location Point ID Unit Description Constr.Datel Process Rate 

(ton cllnker/brl 

Exbaust Flowrate 

(acfm) 

Max Hours of Opel1ltiol 

(hours/dav) 

NOx 

(Ib/ton d.) 

PM10 

(gr/acfm) 

S01 

(lb/ton d.) 

Tulsa EP15 Kiln, sc:ew conwyors, & bucket elevator 1964 96.0 160,000 24 6 0.025 10 

Source Emissions 

Location Point ID Unit Description S01 SO. NO. 

Emissions' 

PM10 

(Ib/hr) 

PM (coarse)' PM (fine)' EC SOA 

Tulsa EPI5 Kiln, screw conveyors, & bucket elevator 960.00 12.91 576.00 17.143 1.149 1.320 0.000 1.766 
• >;;.': 

Printed 1013012006 

Source Emissions Emission Factorl,l 

Location Point ID Unit Description IConstr. Datel Heat Duty SOl 

MMBtu/hr Ib/MMBtu 

discharge chute, separator, bucket elev. I 1964 1.2 0.001 

NO. & S02 emission factors based on AP-42 (7/98) Table 1.4-1 and 1.4-2,PMIO based on baghouse manufacturer exhaust loading rate. 

2 lb/MMSCF converted to IblMMBtu using a natural gas heating value of 1020 Btu/SCF 

Source Emissiolls Emission Data Inpue 

Location Point ID Unit Description SOl S04 NO. PM10 PMco...... PMfine EC SOA 

(lb/hr) 

Tulsa EP04 #2 discharge chute, separator, bucket elev. 0.001 -­ 0.0&9 3.&57 0.964 2.&93 -­ -­

3 Condensible/filterable PM speciation were based on the following profiles: 

For Natural Gas, National Park Service guidance for NG-fircd combustion turoincs. http://www2.nature.nps.gov/air/Permitslect/cctGasFircdCT.cfm 

PM IO•2.

5(filterable) 25% PM (coarse) 

PM2.5 (condensible) 75% PM (fine) 

Lafarge Building Materials Combustion Sources Trinity Consultants 


BART Modeling A-4 063701.0044 


.. 


.'"•.••,.;••,;-. "'-":il!!IIII ••••ti.:<.~. 
•••••••• \~ 

", ,""" ", ''''~' > 

:~""r.'~" ..'r,'. ';"" ~'.'.

' .': 

Printed 1013012006 

Stack parameters (English) 1 

Location Point ID Stack Flow Exit Velocity Stack Height Stack Dia. Stack Temp. Base Elev.1 

ftlsec feet feet Fahrenheit (feet) 

Tulsa EP04 36,000 11.94 100 4 180 875 

Tulsa EP05 4,000 8.29 55 1.6 70 820 

Tulsa EPI3 5,200 9.55 97 1.7 200 851 

Tulsa EPI5 160,000 17.83 60 6.9 500 768 

Tulsa EP24 31,750 10.53 100 4 200 890 

Tulsa EP26 4,000 33.16 64 0.8 70 854 

Tulsa EP25 4,000 33.16 55 0.8 200 845 

Tulsa EPI7 43,000 3.56 53 8.0 200 788 

Tulsa EP27 3,500 4.64 7 2.0 70 754 

Tulsa EP28 3,000 4.91 135 1.8 70 805 

Slack parameters (Metric) 

Location PointlD Exit Velocity Stack Height Stack Dia. Stack Temp Base Elev. 

(m/s) (m) (m) Kelvin (meters) 

Tulsa EP04 3.64 30.49 1.22 355 267 

Tulsa EP05 2.53 16.77 0.49 294 250 

Tulsa EPI3 2.91 29.57 0.52 366 259 

Tulsa EPI5 5.44 18.29 2.10 533 234 

Tulsa EP24 3.21 30.49 1.22 366 271 

Tulsa EP26 10.1l 19.51 0.24 294 260 

Tulsa EP25 10.11 16.77 0.24 366 258 

Tulsa EPI7 1.09 16.16 2.44 366 240 

Tulsa EP27 1.42 2.13 0.61 294 230 

Tulsa EP28 1.50 41.16 0.55 294 245 

Source parameters obtained from Pennit 98-1 74-lV 


2 Base elevations obtained from visual inspection ofUSGS topographic map! 


Lafarge Building Materials Stack Parameters Trinity Consultants 


BART Modeling A-5 063701.0044 

APPENDIX B - DATA FILE LIST 


LIST OF SURFACE STATIONS 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

Station ID Location 

Station No. 

official ID) 

Longitude 

Winfield/Ark. City. KS ___7:..:2:..:.45"'°"'2'--_..........._.I!.,!.§§......... .......~.9.?.QJ.~....... 

Bartlesville. OK ___7:..:2::.35"'6"'5__........ _.J.§.,1§.?.......... _._~017..__ 

-"C"'h""'ildr=es"'s.:...TX!..!.>..____________7'-'2"'3""60""4'--_.............J.~,.~.JJ.................__:lQQ.281.._ 

Coffevville MunicinaL KS 724519 

...... ChanuteiMartin John. KS 13986 

Clinton-Shennan. OK 723526 

E'!QQ 

J.1.§'§'?' .. 

........3..?•.J.J.3.. 

....3??.97... _ 

........... ~~?,4~~.......... 

..........:'!.9•.WQ .... 

Denton. TX 722589 ....3.J.,.7..QQ....... __~. -97.183 

-"E:o..lD=orad=0"-/G"'oo=d:..::Wl'""·n"'"....A~R'________-<-93"-'9"'9=-2__...............3..1,.fJ7...........•..__:92.817 

-'V..::an""c""eC!.AF=.:oB"' lEoe!D1""·d."'-"O""K"-_____.___7:.:.23"'·""53"'5'--_..•..........••J§.J.n.•...............~_-'!.1.9J)_~ 

ShennanlDennison. TX __--'7:..:2""26""'7c.::3__..............._.I~.,.1J7...................... ..:9§,§§7 ... 

......... F..,re"'d"'en""·c",k.:...:O""K"'-__________""22"'2:::2""2__......._... _J.~,JJ.. ........._.... __~~~.983 

.... F-"-lin"'n"'in."'-"AR~_______,____=_55""'5"'5=_5___•...........}§..JQ9........... _--.:9f,467 

Fort Smith. AR 13964 ...1?.J}.~................... _.........~4.}.§7...._ 

FavettevillelDrake. AR 93993}§.,QQ..Q........._...._..~~4,)§I._ 

,""G"'a""l!e"-.""'O""K'_____________~13"'9'_'_7"_5__,_ ...._......... 3.§,.J9.9..........•.•.....___-~9.767_ 

Hobart OK 93986 J.?,..Q99. ....•.......... _----"-9'-"9.""05"'0'---_-1 

Hot Snrinl!s. AR .__--'7C-"2=.::34"'1-=-5__._............J.~d§.7.....•........•.._. __~.91Q~~_.. 

~H"'arn'-'-"'·so"'nIB"-"!:o"'o"'ne:.e-"'C""o.::...'-'A"'R'--______--"13""9:..<.7.... 1__...............J.§.,.m......................__ ..•.n..LE............ 

.!!M...,c",C",onn"""e""Il..!.AFB",,-,,~.K=S________--"3""92"'3'---__.................3..1:.!i.!.?.•_ ........_ .•....,_~.9.?,.~7..,"_. 

-'-W'-'i"'chi-"-·ta.""'-'K~S~__________-"3'-"92"'8"--__...... _..........E.§..?Q............. _ •._:91A3~._ 

Jonlin Re!!ionaL MO 723495 

Lawton. OK 99999 

LiberaL KS 11111 

.3..1,)'?'Q 

..3.A,'?'!.i.?. 

.....J.7Jl.11. 

....:94.5iJQ._ 

............. :9~Al7... 

.... :100,959 

""L""itt"'le"-'R"'0"'c",ki""A"dam""""-'s...!:,AR=-_______--"13""9""6""3___..........•J.4,.n..3....._..........___-9.b.~d.L........ 

-"L"'itto!!le"-R-"'0"'c"'k'-"AFB~~.-"'AR'-"'-________7.!..:.2"'3~40~5'--_ ..............__~.4,.9J7..._•..........._-..::.-?_'.92"'_.!,,15~0:......._1 

Muskol!eelDavis Field. OK .____..!.1-"1l'-!1""2__.._......_.J...?,§..~.Q .............~_-95.3.§7~ 

McAlester Re!!ional. OK 93950 ...............H,2,Q9.......... .........:9.~...m....... 

Mena International Mun.• AR ___7,-,2",34",3<=5,--_.........._.....3..~•. ?_~9............_...__-9.1:217 

.. -"'O"'k1""ab-"-0"'ma~C"'i""tv~.-"'O"'K'--_________"13"'9""6L.7___..........J..?.,1QQ........_.. __..:27..60Q.. _ 

.!-'N"'onnan'-"-""'-'.:-W"-e""s""th""e"!!irn"'e"'r.'-'O~K"'-_______7~2"'35""'7'-"0:......._ ......_.........,}J.,fJJ........_ .... _ ......:'t.7..42..Q._. 

.::.Mo=e=dl"'·c~ino=.e-=L""od=.:e:"'e"-'.K""S"--________7:..:2:..!.45""2""0'--_...... _ .....J_7.,.m_......_._................. ;~~.?~Q .... 

.!..P"'on""c"'-a~C:!!itvc:.:.'-'O"-'K~_________---"13""'9'-"6"-9__....._.........J.~•.m..................~,lQ.O__ 

.!..P",ars~o""n""sIT-'!.!ri,""C""itv'-!"'-"KS~_________7'-"4c.::54"'3~1___..........._.3..1,.m............_-...:91"lli........... 

PariS/Cox Field. TX __---"'33"-=3"'3"'-3__..............J.~..§.~.~...__........._._...._:95.42Q__. 

...Rol!ers.AR __-'6"'6"'66"'6___.........._J.§.,},§..?...._........_.:94.J..Q() .. 

Snrinl!field. MO 13995 _..........._.E.7.JJ....._.........._ -93.3.8L_. 

Lafarge Building Materials 

Refined BART Modeling Protocol 

Trinity Consultants 

06370l.0044 
• • • • • • • • • • • • • • • • • • • • • • • • • • 

• • • • • • 

--

------------------_......._...................__........-....-----.---..­

-"S"il~o",am~S!!!on",,·n~e:,,-,.AR~_________....!7.!...77~7'-'7__.._ ...J.(!,m_ ..._....._. ......:9.1.A~.~ 

-,-W,-,i",ch""ita",....F.."al""ls,,-.TX~________.__J}2§..§...___.. _.........;3.}.,.9~........................:9.l!,?J.7...... 

-"S"ti",ll",w""ate",r"-.."'O"'K'--_________---'I..!.I'-!1l"'3'___......_.......3.JiJJ9_........__........ .:91J!~~ .. 

.!..F"'ort"""Le"'o"'n"'ar""d....!W-'-'o"'od.""-'M~O'_________'7c:2""44"'5'-'-7__._....._.3.lJ.~.Q.__.... _.........:-""92=-....,15"'0'--_1 

-,T....ink=eC!...rAFB=.;"".'-'O"-'K~_________-'-'72=3""54..:.;0'--_._....._J~,.~n?._ ......_-----.::2.:UllL_. 

-"M:..:c::>.Kinn,.·=ev.:..:.-'TX~__________7,..22""'-54_'_'1'___.............._•.~}J83 .. ...:.96.583 

....T"'u"'ls>:.a... In"'t-'-Airn ..·'-""'ort~.'-'O~K~_ _,----_---"13"'9"'6:.><8--...- ........~(!,7.QQ_... _ ... ...:.95... ~.~3

....T.!::ex"'ar"'kan"""a"..R.,e""I!1"'·o"'n"'a1.!:....!.TX"-'-_______..!.\3"'9'-'7.!...7__._......J},1.~.Q.. _....... ~,000__ 

Bentonville. AR _____--"'88"'8,.,,8""8__...........J.(!,J.?Q...........__.... ..:91.,71.7... 

........ EnidIWoodrine:. OK 44444 .3..(!,J.~3. ... ~n7~} 

Ardmore Municinal. OK 723555 34.300 -97.017 

LIST OF PRECIPITATION STATIONS 

Cooperative Latitude Longitude 

station ID Location/name State (deg.) (deg.) 

30130 ALUM FORK AR 34.480 ·92.510 

30178 ANTOINE AR ....34.020 _._._..._._ -93.250 ..­ ---'~'~ 

30220 ARKADELPHIA 2 N AR 34.090 ·93.030 

30458 BATESVILLE LIVESTOCK AR 35.500 -91.460 

30530 BEEBE AR 35.040 -91.540 

30616 BERRYVILLE 5 NW AR 36.260 -93.370 

30764 '._._-_... 	 BLAKELY MOUNTAIN DAM ..._ 34.340 AR ..H.__..._.H.._ .....__-93.120 

..H.___._.

__ 

30798 ••••__ __ ..._.... 	 BLUE MOUNTAIN DAM AR ......._-_35;070 ..._..... ... ......- ........__.....__.
 •••H. -93.390

•

• 

...­

....

30832 __._. BOONEVILLE 3 SSE AR 35.060 -93.550 

30842 BOTKINBURG 3 NE AR 35.430 -92.280 

_ 

30900 .. BRIGGSVILLE AR 34.560 -93.300 

31020 	 BULL SHOALS DAM AR 36.220 -92.340 --_..."._......__....._.,----_...._-_. ...•_.........._.........._-_..­

1--._._- 31140 	 CALIONL&D AR 33.180 -92.290 ... 

31152 _._-,-_.-	 CAMDEN I AR 33.3~._ -92;~ 

31457 CLARKSVILLE 6 NE AR 35.320 -93.240 

31582 COMPTON AR 36.050 -93.180 

31952 DE QUEEN DAM AR 34.060 ·94.220 '._--_.. ._.._-_.._._-_._-----_._. 

32020 	 _....__._..... DIERKS DAM AR ....._ 34.090 ...__..._ ..._ -94.050 ...__.__._-_. 

DUMAS AR 33.530 -91.290 ..

•

•

• 

32300 EL DORADO GOODWIN FD AR 33.130 -92.490 

32356 EUREKA SPRINGS 3 WNW AR 36.250 -93.470 

32444 FAYETTEVILLE EXP STN AR 36.0~__. -94..!.D0 __ ..- ......---.... 

32489 	 .._.__.._..... FERNDALE 6E AR 34.460 -92.270 --'-"'-'---'-''''--'­

~~-. FOREMAN AR 33.430 ._.-94:.~ •

32574 FORT SMITH MU, OK AR 35.200 ·94.220 •

•

•

• 

Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•

• 
' 

• • • • • • • • • • • • • • • • • • •• 

r. 


•

'.:.'.'.'. 

h •


~~.' ~~l'l' ' .. .. ~ ,:'

I:


.........• 
 :,:'

t~. 

~;~.

i,·

~.". 'l 

•

•

•

•

•

•
•

• 


-.• 

•

•

• 


• 

32794 GILBERT AR 35.590. ·92,430. 

32810 GILLHAM AR 34.120. -94.150. 

AR .... 35310 _.__...._--<_.. :2.~:goO ...._ 

AR 36.170. -91.280 .._ 

33165 HARRISON FAA AP AR 36.160. -93.0.90. .....­

~5_44__..___...::HUN=-'-T:..:S=-V:...:I=LL=E=-:cl-=S.::.S-'-W _______--'.:=- AR __::.::.:.=-:..;: 36.0.40.____-93,450. ~=-= .._ 

34185 LEWISVILLE AR 33.220. -93340. 

34248 ..._ .....__....... __ N_O_R_TH LITTLE ROCKIMUNIC:c;:I;;.:PAc::L:::.._--:-=c::_ AR 34.50.0. -92.160. ..._........._._.__•......_...._. 

34756 AR 34.340. -94.160. ..... _.........._._.. .. ~ .........---.....~.--.


~~83~._...____.;;..M-"I=;;Lc;;:L_W;...;O=-O:..:D'--"-D"-'AM=:..._______.:c='---­ AR 33,410. -93.590. -~------..-.----. 

34988 MOUNT IDA 3 SE AR 34.320. -93.350. 

AR 34.0.90. -93,430. .._._-_......._._...._- ........--......... ~ ....
~ 

AR -_....._._33.560. ......__._..__. :2}.,?!.Q..._ 


AR ... }..±:~?L .._:2.3..:.1..gg__ 


35~_28_...__.__._..::N~O~R~F~O~RK~D~A~M~______._~AR __~~~ 36.150. __-92.150.
 ~~_-1 

35320. NORTH LITTLE ROCKIMUNICIPAL AR 34.50.0. -92.160. 


3560.2 PARTHENON AR 35.570. -93.150. 


35754 PINE BLUFF AR 34.140. -92.0.10. 


_3590.8
 ....__....,._-_•.... 


~~_39_3._•.___-=S...oT..::P.:..A:.:U:.::L=---__________-:..;A::.;R'--_ _'3::..:5... 490_.,.........:23.46.L 


~._____"S~T:..:U:..:TT:..::..:G~A~R~T=-9=-=E~SE~______..::A~R~_~3~4.=28~0'~__-.:..9~1.~25~0.~ 

37488~__.___W~A~L=D~R:c;:O~N~___________~A~R~_~3~4.~54~0'~.__.~94~.~0..6~0.,--~ 

141427 CHANUTE M JOHNSON AP KS 37,40.0. -95.290. 


141699 COLBY 1 SW KS 39.240. -10.1.0.40.
 .--.._--_._--_.__.._--_._-----"':::;.._..........._...........­

142164 DODGE CITY MU, KS KS 37.460. -99.580. _'___

142432 ELKHART KS 37.0.0.0. ·10.1.530. 

142980. GARDEN CITY EXP STN KS 38.0.0.0. -10.0..490. ..._-_._-_.__...__.__._-_._--_....................- ........... 


143954 ...__IN_D_E_P~_N_D_E_N_CE___.______--'K:.:::S:::,-_ ........}1:1iQ._' ......._:2.S.:i~9.._ .. 


-~....--..-..­

_ou..._ •••____ 

!~398:!._...•_ .•••_..:cIO~L::.:A_=_.:..I...:.W'---_._______.__-=K.:::S"-_..J!c~Q._ .._.:2?.J..~.L... 

!.451 ~5_..... MEDICINE LODGE KS 37.170. -98.330. 

147922 SUBLETTE KS 37.290. -10.0..50.0. 

148830. WICHITA KS. KS 37390. ·97.260. 

230.789 BOLIVAR 

~ll?..___JOPLIN MUNICIPAL AP MO 37.0.90. -94.30.0. 

~_45_44_______..._..::K=I~RK~S~V=IL~L~E~_______________~M~O~__~~~~ 40..120. ____~~~ -92340...._ 

234825 LEBANON 2 W MO 37.410. -92.420. 

234919 ._-=L:.:;IC::::K:.:;IN:.:..:.;::G:...4:..;N"'--_______---'M=O __.::.:.:=~ 37.330. ___-9!~ :;.,: ~.3...5~?___._._~O_C_KW_O_O_D____.__._..._____'M=O_ 37.230. -93.570. ••••••• n __••••••• _ ••••_ •••••_ •••_. 

23530.7 MARSHFIELD MO 37.20.0. -92.550. .._--_. ~355~~ ___ .__-=-M-=I=LL=E=R~I_=E'--____________'M:..=O ___=..:.~:........ 37.130. __-'_=:.:::..::.:........_1 -93.50.0.

236460. OZARK BEACH MO 36.40.0. -93~~ 

Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • • • •• 

• • • • 

• • • 

•

•
•

•

•
•

• 


•

•

•

•
•

•

•

• 


237435 ST JOSEPH ROSECRANS MO 3M70 -94.550 

MO 37380 -91.320 

MO 36.590 -94.320 

MO 37.140 -93.230 

MO 37.420 -93.460 

MO 36.360 -93.180 

OK 36.450 -96.000 

I::c.:.=.=-----.='-'.::::'-=----------Q!S----. 34:.?~. -95.050 

LAKE OVERHOLSER 

LEHIGH 

MARSHALL 

MCALESTER MUNI AP 

346130 MUSKOGEE 

~:4.648:5.__.. NOWATA 

346616 

~~~~"'M'~_M 

_ ..._ OILT0!i2.~E 

346638 OKEMAH ""M_~"__......_ .._. 

346661 .__» OKLAHOMA CITY, OK 

346944 PAWNEE5N 

347201 PONCA CITY MUNI AP 

347~___.. RANGE 

RIVERSIDE 4 W 

348501 STILLWATER 2 W 

348992 TULSNINTL OK. 

349404 WAYNOKA 

349629 WICHITA MTN WL REF M""""_""_'"'' 

OK 

OK 34.290 -97.580 

OK ...••... 36.180 "......_...... ,,, 

OK 34350 ·99.020 

OK 35.290 -97.400 

OK 34.280 -96.130 

OK 36.040 -96.340 ........._ •••_ •••••• "_•••••m~•••' __ ""'-.-...-~.~--."".--...­. 

OK ....2.?:.~Q ____:2..(i·18~_~ 

OK 35.230 -97.360 

OK 36.240 -96.490 

OK 36.440 -97.060 

OK 36.330 -101.050 ...--.-~-

OK 36.470 -100.250 -..-.~--.--~.-".. 

-99.530 

OK 36.070 -97.060 

OK 36.120 -95.530 

OK 36.320 -98.530 

OK 34.440 -98.430 

..M. ___.........__...._ ..... 

TX 34.260 -100.170 

•

•

•

•

•

•

•

•

•

•

•

• Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 
• • •• 

CLARKSVILLE 1 W TX 33370 


COMMERCE 


DALLAS-FORT WORTWFORT. TX. TX 32.540 -97.010 


DALLAS LOVE FIELD TX 32.510 -96.510 
 '-'-' 

DENTON2SE 	 TX 33.120 -97.060 


TX 32.490 -97.210 


TX 


TX 33.090 -96.490 
 .."'-,..,,-""~.,,-,...~-. 

TX 32.570 -97.030 


TX 34.490 -102.240 


TX 33.390 -101.490 

MIDLOTHIAN 2 TX 32.290 -97.000 

WELLS MUN AP TX 32.470 -98.040 

____.~...MQUNT PLEASANT TX 33.100 ..,.._._.m..._._ -95.000...

NACOGDOCHES 	 TX 31.370 -94.390 ._._--_....._­

TX 35.340 

TX 32.150 -98.120 

TX 33.090 -95.380 

:~ ., 

TX 33.590 -98.300 

•

• 

• Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Protocol 063701.0044 •.,' 

•

•

• 

.~~ '~hr.." ,. #"". 
---------

-;...........".. 
 ••iifi,.·••••• ~, "it -j


" 

List of USGSTerrain-Files (1-<1eg format) 

37 

E,' 361 

35 

34 

33 

Dodge City-E Pratt-W Pratt-E Wiehila-W Wichita-E Joplin-W Joplin-E Springfield-W Springfield-E Rolla-W 

Perryton-E Woodward-W Woodward-E Enid-W Enid-E Tulsa-W Tulsa-E Hanisan-W Harrison-E Poplar Bluff-W 

PI.inview-E Clinton-W Clinton-E Oklahoma City-W Oklahoma City-E Fort Smith-W Fort Smith-E Russellville-W Russellville-E Memphis-W 

Amarillo-E L.wlon-W Lawton-E Ardmore-W Ardmore-E Me Alesler-W Me Alesler-E Little Rock-W Little Roek-E Helen.-W 

Lubbock-E Wichita Falls-W Wichita Falls-E Shennan-W Shennan-E Texarkana-W Texarkana-E EI Dorado-W EI Dorado-E Greenwood-W 

]0] ]00 99 98 97 96 95 94 93 92 

t, 

t

r; 
,,,_t_t_~' ••••

 , , '; '~!,:'" . ' 

List of USGS Land Use Land Co\ier Le) Files (1 :250~OOO;krh Format) 

38.5 

37.5 

36.5 

35.5 

34.5 

Dodge City Pratt Wichita Joplin Springfield Rolla 

Perryton Woodward Enid Tulsa Harrison Poplar Bluff 

Plainview Clinton Oklahoma City Fort Smith Russellville Memphis 

Amarillo Lawton Ardmore McAlester Little Rock Helena 

Lubbock 

- Wichita Falls Sherman Texarkana El Dorado Greenwood 

102 100 98 96 94 92 

:1 

~ 

t: 

r>: