° ° ~ Ttnt:lllty~ Offices nationwide. trinityconsultants.com

\JJllSU tants 

1120 East Sheridan, Suite 205, Oklahoma City, OK 73104 U.S.A.• (405) 228-3292 • Fax (405) 228-3293 

December 1, 2006 

Mr. Eddie Terrill 

Director, Air Quality Division 

Oklahoma Department of Environmental Quality 

P.O. Box 1677 

Oklahoma City, OK 73101-1677 

Re: 	 Visibility Modeling and BART Waiver Request 

Lafarge Building Materials, Inc. 

Tulsa Cement Plant 

Dear Mr. Terrill: 

RECEIVED 


[ DECIl~ 7006 ] 

AIR QUAbiTY 

On behalf of the Lafarge Building Materials, Inc. (Lafarge), Trinity Consultants is pleased to submit this 

BART Wavier Request for the Tulsa Cement Plant. The Tulsa Plant has been identified as potentially 

subject to BART under the Regional Haze Program. Consequently, Lafarge contracted with Trinity 

Consultants to'determine the visibility impact of the facility on the nearby Class I Areas. The report of 

the modeling study is enclosed. The CALPUFF refined modeling conducted by Trinity confirms that the 

visibility impact is less than the 0.5 deciview threshold at Class I receptors located within 300 kilometers 

from the Plant. . 

, 

In accordance with OAC 252: 100-8-73, Lafarge is requesting a waiver from DEQ that a BART 

determination for the Tulsa Cement Plant is not required. If you have any questions or need additional 

information, please call Ms. Pamela Bennett at (918) 437-3902 ext. 271 or me at (405) 228-3292. 

Sincerely, 

/ 

..r 

Vern Choque e 

Principal Consultant 

Encl: 	 Visibility Modeling Report 

Copy: 	 Pamela Bennett, Lafarge Building Materials, Inc. 
••••

•
••

•
••••

•
•••••••••

•
••••••••

•
•
••

•
 '. 


CALPUFF BART MODELING REpORT • REFINED ANALYSIS 

LAFARGE BUILDING MATERIALS INC. 

GE 

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 

(405) 228-3292 

November 29, 2006 

Project 063701.0044 

REC'El'J~-D 

IOEC(1.1 7006 1 

AIR QUALITY 

••••

•
•
• 
-
•
•

•
• 
 • TABLE OF CONTENTS •

• 
 • 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
 •

• 
 1.4 MODEL RESULTS .................................................................................................. 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. CALPUFF ........................................................................................................... 4-1 


•

•

•

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


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


•

• 


4.2.1 PM10 EMISSIONS.................................................................................................. 4-2 


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


4.3 RECEYfOR LOCATIONS.......................................................................................... 4-3 


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


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


4.6 CALPUFF MODEL CONTROL PARAMETERS ........................................................ 4-7 


•

•

•

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


5.1 LIGHT EXTINCTION ALGORITHM........................................................................... 5-1 
 • 5.2 CALPOST PROCESSING METHOD......................................................................... 5-1
 •

• 


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


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


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

• 6. MODEL RESULTS ...•......••....••....•••••.................••....•.....•...............•......•.....••.••....... 6... 1


•

• 
 •

Lafarge Building Materials Trinity Consultants •

•

•

• 

Refined BART Modeling Report 063701.0044 
• •••••

• ApPENDIX A - EMISSION DATA 

APPENDIX B - DATA FILE LIST 

•

APPENDIX C - MODEL INPUT/OUTPUT FILES

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

• 

ii Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

••••

• 
• • • • • • • • • 

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

• • • • • • • 

•

•

•

•

• LIST OF FIGURES 

TABLE 2-1. 

TABLE 2-2. 

TABLE 3-1. 

TABLE 4-1. 

TABLE 4-2. 

TABLE 4-3. •

• TABLE 4-4. 

TABLE 4-4. 

TABLE5-1. 

TABLE 5-2. 

TABLE 5-3. 

TABLE 6-1. 

TABLE 6-2. 

TABLE 6-3. 

TABLE 6-4. 

•

•

•
•

•

• 


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


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


SUMMARY OF CALMET INPUTS ..................................................................................... 3-4 


CONSTRUCTION DATE SUMMARY ...................................................................................4-1 


ESTIMATED 24-HOURAVERAGE EMISSION RATES ..................................................:......4-2 


SPECIATED PMIO EMISSIONS........................................................................................... 4-3 


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


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


MONTHLY HUMIDITY FACTORS ...................................................................................... 5-1 


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


SUMMARY OF CALPOST INPUTS ................................................................................... 5-2 


VISIBLITY IMPACT MODEL RESULTS FROM CANEY CREEK 2001-03 .............................. 6-1 


VISIBLITY IMPACT MODEL RESULTS FROM HERCULES-GLADES 2001-03 ..................... 6-2 


VISIBLITY IMPACT MODEL RESULTS FROM UPPER BUFFALO 2001-03 ........................... 6-3 


VISIBLITY IMPACT MODEL RESULTS FROM WICHITA MOUNTAINS 2001-03 .................. 6-4 


• 

iii Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
• • • • • • 

• • • 

•

•

•

•

•

• • LIST OF TABLES

•

•

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


• 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 LOCATIONS ....................................................4-4 


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

•

• 


FIGURE 4-3. UPPER BUFFALO WILDERNESS RECEPTOR LOCATIONS .................................................4-5 


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

•

•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

Lafarge Building Materials IV Trinity Consultants 

•

• Refined BART Modeling Report 063701.0044 
• • • • 

• • • • • • 

•
•

•

• 


1. EXECUTIVE SUMMARY 

•

• 

• 

•

•

•

• 1.1 BACKGROUND 

On July 1, 1999, the U.S. Environmental Protection Agency (EPA) published the final 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 defines BART-eligible sources as sources that meet the following 

criteria: •

•

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


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

(3) Are listed as one of the 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 quality 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 if the results ofthe 

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 of BART modeling protocols and to

•

• 

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

(CENRAP) organization developed BART Modeling 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 of this document is to provide a report summarizing the model results and the modeling 

methods and procedures that were 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 

Refined BART Modeling Report 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 Apperidix A at the Tulsa Plant as sources that • meet the three criteria for being BART-eligible sources. These sources include:

.it. Unit 221, Raw Material Grinding #2

•

•

• 

.it. Unit 222, Auxiliary Raw Material Grinding #2 


.it. Unit 223, Auxiliary Raw Material Grinding #2 


.it. Unit 404, Kiln #2 Feed System 


.it. Unit 502, Kiln #2 


•

• 

•

•

• 

.it. Unit 706, .Finish Grinding #2 


.it. Unit 705, Finish Grinding #2 System 


.it. Unit 704, Auxiliary Finish Grinding #2 


.it. Unit 602, Clinker Cooler #2 


.it. Unit 1101, Coal Unloading from Railcar 


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

.it. 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 ofthe Tulsa Plant in relation to its neighboring Class I areas. 

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

included in this refined 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.4 MODEL RESULTS 

Model results from the refined analysis are included in Section 6. The tables in Section 6 indicate

•

•

• 

that the 98th percentile visibility impairment values are below the 0.5 dv threshold at each of the four 

Class I areas within 300 km. As described in CENRAP's BART Modeling Guidelines (December 15, 

2005), these results indicate that the Lafarge Tulsa Plant is eligible for a waiver from BART 

requirements.

•

•

•
•

• 

1-2 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

•

•

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

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

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

1-1 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 
• • • • • • • • • 

•

•

•

•

•

• 2. CALPUFF MODELING SYSTEM 

•

•

• 

• The main components ofthe 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 ofthe CALPUFF modeling system programs that Trinity used 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 ofthe program versions are not the 

•

•

• 

same as the program versions recommended by CENRAP. Certain programs are being utilized in 

order to avoid data compatibility issues with previous versions ofthe CALPUFF sy stern. Trinity used 

the following software version numbers listed below. Model selection was primarily driven by a 

desire to utilize the most recent release ofCALMET. CALPUFF and CALPOST executable versions 

were then selected as appropriate for compatibility. 

TABLE 2-1. CALPUFF MODELING SYSTEM VERSION NUMBERS 

•

'. 

•

•

•

•

•

•

•

•

• 

•

•

• 

Program CENRAP Suggested 

Version Level 

TERREL 3.311 030709 

CTGCOMP 2.42 030709 

CTGPROC 2.42 030709 

MAKEGEO 2.22 030709 

CALMM5 2.4 050413 

CALMET 5.53a 040716 

CALPUFF 5.711a 040716 

POS~18 CALPOST 5.51 030709 

Lafarge Refined 

Analysis 

Version Level 

3.311 030709 

2.22 030528 

2.42 030709 

2.22 030709 

2.4 050413 

6.211 060414 

6.112 060412 

1.4 040818 

6.131 060410 

Reason for Difference 

Version recommended is not 

available 

Most recent release of 

CALMET 

Use version compatible with 

CALMET v6.211 

Use version compatible with 

CALPUFF 6,112 

2-1 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 
• • • • 

• • • • • • • • • 

•
•

•
•

•

•

• 


•

•
•

•

•
•

•
•

•

.'
•

•

•

•

•
•

• 


-.• 

•

•

• 


•

• 


2.2 MODELING DOMAIN 

Trinity utilized a modeling domain that extends at least 50 km in all directions beyond the TulsaPlant 

and each ofthe fo'ur surrounding Class I areas. The map projection for the modeling domain is 

Lambert Conformal Conic (LCe) 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 corner of the modeling domain is Latitude 33.25°N, Longitude 100.5°W which was 

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

approximately Latitude 37.25°N, Longitude 92°W. At a grid spacing of3 km, the number ofX 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 .§ 

u 

§ 250 o 

C 

<3 200 " J:> '" E 

.'3 150 

100 

50 

FIGURE 2-1. PROPOSED MODELING DOMAIN 

r----------------------------------------I 

I m Joplin_MO 

Hercules-Glades • 

Lafarge_itulsa_PIant 

Ell EB Tulsa_OK Fayettevllle_AR • 

Upper_Buffalo 

~---------------------------------- ______I 

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

Lambert Confonnal Conic (X) (km) 

For the purposes of this analysis, the CALMET and CALPUFF domains are identical for each 

analysis and are singularly referred to herein as the "domain." 

2-2
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 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 of3 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 trarlsformations take place. 

TABLE 2-2. VERTICAL GRIn 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 Report 063701.0044 

• 

•

• 
•

•

•

•

• • 3. CALMET • 

Trinity conducted a three-year CALMET analysis incorporating both mesoscale model and 

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

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

•

•

•

•

• 

model parameter settings. The CALMET modeling procedures that Trinity used generally follows 

• 

the recommendations in CENRAP's protocol. It should be noted that Trinity used 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 used terrain data obtained from the United States Geological Survey (USGS) in 1­

degree (1 :250,000 scale or approximately 90 meter resolution) digital format. A list ofthe 

• 

USGS terrain files used in the analyses is provided in Appendix B. A plot ofthe land •

• elevation for the modeling domain based on the referenced files is provided in Figure 3-1. 

•

•
•

•

•

•
•

•

•

•

•

•

• 

3-1 • Lafarge Building Materials Trinity Consultants • Refined BART Modeling Report 063701.0044 •

•

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

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

Longitude (degrees W) 

300 

200 

100 

o 

800 


700 


600 


500 


Meters (Above MSL) 

3-1
 Lafarge Building Materials Trinity Consultants 


Refilled BART Modeling Protocol 063701.0044 

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

•
•

•

•

•

•

•

•

•

•
•

•

•
•

•

•

•

•
•

•
•

•

•

• 


Trinity used the TERREL progralll to generate grid-cell elevation averages from the USGS 

terrain data across the modeling domain. 

3.1.2 LAND USE DATA 

Trinity used USGS Composite Theme Grid (CTG) format Land Use and Land Cover 

(LULC) data files at I:2S0,000 resolutions. Where I :2S0,000 land use data was not 

available, Trinity used USGS data at I: 100,000 resolutions. A list of the USGS land use 

files that were used is provided in Appendix B. 

Trinity used 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 compressed the files 

using the program CTGCOMP. 

3.2 METEOROLOGICAL DATA 

CALMET was 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 was used to supplement the hourly surface, upper air, and 

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

field for the CALMET simulations. Trinity used the following Sth generation mesoscale 

model meteorological data (or MMS data): 

• 2001 MMS data (12 km resolution), 

• 2002 MMS data (36 km resolution), and 

• 2003 MMS data (36 km resolution) 

3.2.2 SURFACE DATA 

The use ofmultiple stations for meteorological observations in CALMETICALPUFF 

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 were 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 of the domain. Raw observations from these stations were obtained from 

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

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

each year analyzed. 

3-3 Lafarge Building Materials Trinity Consultants • Refined BART Modeling Report 063701.0044 •

•
• 
~. 

FIGURE 3-2. PLOT OF SURFACE STATIONS 

KTBN . .~. 

500 I<I()iAS .. I<PPF --~- ....... -- - --1 
 KWLI) . !<I.BL . 400 . ~O KPNC Hercules-Glades 

•

•

•

•

•

•

•

•

•

•

•

•

• 

E ~G KHRO ~.L.P 6 ~G

. 
 KSWO 2:. • u 

c 300 0 

() • (i} KFSM I<IJK1S11K .. . E I KO.,\JN 0 

c I I<MLC 8 200 I . 1::: 

•

•

•

•

•

• (I) 

. ..CI Wichita_Mountainsl ~ 

E 

<ll ~DM 

-l 

100 KBPS . KPRX . KTXK

• !<ELI)

. __________~_~ ______________ 1

•

•

o L

•

•

•

• 

o 100 200 300 400 500 600 700 800 

Lambert Conformal Conic (Xl (km) 

3.2.3 UPPER AIR DATA 

Upper air observations were not used in this CALMET analysis (i.e., NOOBS = I). 

3.2.4 PRECIPITATION DATA 

The effects of chemical transfonnations and deposition processes on ambient pollutant 

concentrations were considered in this analysis. Therefore; it was necessary to include 

observations ofprecipitation in the CALMET analysis. The precipitation stations from 

which Trinity extracted 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 

of the domain. Data from the stations were processed for use in CALMET using Version 

5.31, Level 030528 ofEPA's PMERGE program. 

•

•

•

•

•

•

•

•

•

•

•

• 

•

• 

3-4 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

•

• 

•

• 

i<MKO 

KSL() . . 

I4AWT . 
• • • 

- --

1 

• 


•

• 

•

•

•

•

•

•

• 	 500 

400 

E 

.2 

c: 300 

•

•

•

•

•

• 0 

() 

ii 

E 

.s 

c: •

• 

() 

0 200 


1:: 


Q) 


.0 

E 

IV 

...I 

100 

o 

•

•

•

•

•

•

•

• 

FIGURE 3-3. PLOT OF PRECIPITATION STATIONS 

IO~l ST6~VA LE~O - -IMCHIT- - - - - -€~ LOC1Mn.. - r;w!SHF ­

JO~N M"lllJfliNG 

MEbICI 	 INDEPe 

Gltl!~ I 

EIOng~lbs 

I 

I 

Ge.!Ry 

~O Oi<E~

CHI~KA I 

Wichita Mountaini 

I 

LE~GH I CHltoR C.llNEY 

I 

FOR~MA MILLw) 

~~O

I TEXAAK ~sv

J..cf<ss OB>1rON 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 of the CALMET parameters that Trinity used 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 of the 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-5 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 
..,


•

•

•

•

•

• 


TABLE 3-1. SUMMARY OF CALMET INPUTS 


•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

•
•

•
•

•

•

•

•

•

•

•

•
•

•

•

•

•

•

•

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Used Notes 

NUSTA Number ofupper air 

data sites 

0 0 Only surface and 

precipitation 

observations will be 

used 

Number ofoverwater 

data sites 

0 0 

IBYR Starting year 2001 Appropriate 

met year 

Years 2001, 2002, 

2003 

IBMO Starting month 1 Appropriate 

month 

Due to file size, 

analysis will be 

completed by month 

IBHR Starting hour 1 1 

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) 

1 1 

LCALGRD Compute CALFRlD 

data fields (T run 

CALGRlD) 

F 

" 

F 

!TEST Stop run after SETUP 

to do input QA (2 == 

run) 

2 2 

PMAP Map projection LCC LCC 

RLATO Latitude (decimal 

degrees) ofprojection 

origin 

40N 3325N 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 

XLAT2 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 ofX grid 

cells in 

meteorological grid 

300 256 Appropriate for 

domain 

3-6 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

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

• 

CALMET Variable Description 

NY Number ofY grid 

cells in 

meteorological grid 

DGRlDKi\1 Grid spacing (krn) 

XORIGKM Ref. coordinate of 

SW corner of grid 

cell 

YORIGKM Ref. coordinate of 

SW corner of grid 

cell 

NZ Number of vertical 

layers 

ZFACE Vertical cell face 

heights (NZ + I 

values) 

LSAVE Save met. data fields 

in an unformatted 

file? 

IFORMO Type of unformatted 

output file (I for 

CALPUFF) 

LPRINT Print met. fields 

IPRINF Print intervals 

IUVOUT(NZ) SpecifY layers ofu,v 

wind components to 

print 

IWOUT(NZ) SpecifY layers ofw 

wind component to 

print 

ITOUT(NZ) SpecifY layers of 3D 

temperature field to 

print 

LOB Print met data and 

variables 

NNI First time step for 

debug data to be 

printed 

NN2 Last time step for 

debug data to be 

printed 

IOUTD Control variable for 

Value Included in 

CENRAP Protocol 

192 

6.0 

-1008 

0.0 

10 

0, 20, 40, 80, 160, 

320,640,1200,2000, 

3000,4000 

T 

I 

F 

I 

NZ*O 

NZ*O 

NZ*O 

F 

1 

1 

0 

Value Trinity 

Used 

155 

3.0 

0 

0 

10 

0, 20, 40, 80, 

160, 320, 640, 

1200,2000, 

3000,4000 

T 

I 

F 

I 

NZ*O 

NZ*O 

NZ*O 

F 

1 

2 

0 

Notes 

Appropriate for 

domain 

Refined grid 

spacing 

Appropriate for 

domain 

Appropriate for 

domain 

Will generate debug 

data for a total of2 

time steps 

3-7 

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

• Lafarge Building Materials Trinity Consultants • Refined BART Modeling Report 063701.0044 •

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

•
••••

•
••••

•
••

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

• 

CALMET Variable 

.. 

Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Used Notes 

writing test/debug 

wind fields 

NZPRN2 Number of levels 

starting at surface to 

print 

0 1 Default 

IPRO Print interpolated 

wind components 

0 0 

IPRI Print terrain adjusted 

surface wind 

components 

0 0 

IPR2 Print initial 

divergence fields 

0 0 

IPRJ Print final wind speed 

and direction 

0 0 

IPR4 Print final divergence 

fields 

0 0 

IPR5 Print winds after 

kinematic effects 

0 0 

IPR6 Print winds after 

Froude number 

adjustment 

0 0 

IPR7 Print winds after 

slope flows are added 

0 0 

IPR8 Print final wind field 

components 

0 0 

NOOBS No observation mode 

(2 No surface, 

overwater, or upper 

air observations; use 

MM5 for surface, 

overwater, and upper 

air data) 

2 1 Will use surface and 

precipitation 

observations 

NSSTA Number of 

meteorological 

stations in 

SURFDAT file 

0 47 Number of stations 

NPSTA Number of 

precipitation stations 

in PRECIP.DAT file 

0 135 Number of stations 

ICLOUD Gridded cloud fields 

(0 =: no, 3 Gridded 

cloud cover from 

prognostic relative 

3 3 

••

• 

3-8 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 

•

• :. .. ':" 
• • • • 

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

• • • • • • • 

• •••••

• 

CALMET V,ri,bl, k~Ption Value Included in 

CENRAP Protocol 

Value Trinity 

Used Notes 

IFORMS Fonnat ofsurface 

data (2 = fonnatted) 

2 2 

IFORMP Forrnatof 

precipitation data (2 

= fonnatted) 

2 2 

IFORMC Fonnat of cloud data 

(2 = fonnatted) 

2 1 NIA - No cloud data 

used in model 

IWFCOD Generate winds by 

diagnostic wind 

module? (1 = yes) 

1 1 

IFRADJ Adjust winds using 

Froude number 

effects? (1 

=yes) 

1 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 

ISLOPE Compute slope 

flows? (1 = yes) 

1 1 

IEXTRP Extrapolate surface 

winds to upper layers 

(-I = no extrapolation 

and ignore layer 1 of 

upper air station data) 

-1 -1 

ICALM Extrapolate surface 

winds even ifcalm? 

(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,0 

RMIN2 Minimum vertical 

extrapolation distance 

Distance (km) 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-9 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

• •

• 

e, 
••

•
••••••

•
••

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

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Used Notes 

MM-FDDA data? (14 

= Use winds from 

MM5.DAT as initial 

guess wind field) 

ISTEPPG Timestep (hours) of 

the MM5 data 

1 1 

LVARY Use varying radius of 

influence to develop 

surface winds? 

T T Use FLAG default 

RMAXI Maximum radius of 

influence over land in 

surface layer (lan) 

30 100 As included in 

IWAQM Phase 2 

case study 

RMAX2 Maximum radius of 

influence over land 

aloft (Ian) 

30 500 As included in 

IWAQM Phase 2 

case study 

RMAX3 Maximum radius of 

influence over water 

(Ian) 

50 36 Same as MM5 data 

spacing 

RMIN Minimum radius of 

influence used 

anywhere (lan) 

0.1 0.1 

TERRAD Radius ofinfluence 

of terrain features 

(km) 

12 12 

Rl Weighting offirst 

guess surface field 

(lan) 

1 50 As included in 

IWAQM Phase 2 

case study 

R2 Weighting offirst 

guess aloft field (Ian) 

1 100 As included in 

IWAQM Phase 2 

case study 

RPROG MM5 windfield 

weighting parameter 

(lan) 

0 0 

DIVLIM Maximum acceptable 

divergence 

5.E-6 5.E-6 

NITER Max number of 

passes in divergence 

minimization 

50 50 

NSMTH 

• 

Number ofpasses 

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-10 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 

• 
• • ••

•
••

•
•••

•
••

•
•

•
••

•
•

•
•••••

•
••••

•
•

• 
 •

• 


CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Used Notes 

NITR2 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 

LO 1.0 

ALPHA Kinematic effects 

parameter 

0.1 O.l 

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 

I and NSSTA) 

4 4 

-

IDIOPT2 Domain-averaged 

wind component 

switch 

0 0 

IUPT Station for lapse rates 

(between I 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-11 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 

•

• !'--.'-•. 
••

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

•
••

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

• 

Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 

CALMET Variable Description 

IDIOPT5 Observed aloft wind 

component switch 

LLBREZE Use lake breeze 

module? 

NBOX Number oflake 

breeze regions 

NLB Number of stations 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 

height coefficient 

FCORIOL Absolute value of 

Coriolis parameter 

IAVEZI Conduct spatial 

averaging? (1 = yes) 

MNMDAV Max search radius in 

averaging process 

Jnumber of grid cells) 

HAFANG Half-angle of upwind 

looking cone for 

averaging (degrees) 

ILEVZI Layers of wind 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 of layer above 

current convective 

mixing height 

Value Included in 

CENRAP Protocol 

0 

F 

0 

0 

0 

1.41 

0.15 

2400 

0.l6 

l.E-4 

1 

10 

30 

1 

0.001 

200 

3-12 

Value Trinity 

Used 

0 

F 

0 

0 

0 

1.41 

0.l5 

2400 

0.16 

l.E-4 

1 

10 

30 

1 

0.001 

200 

Notes 

••· ,~.~.L",;...... 
•

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

• 

.'••••••••••••••••

• 

CALMET Variable Description 

Value Included in 

CENRAP Protocol 

Value Trinity 

Used 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 1 Will use surface 

observations and 

precipitation data 

lRAD Type of interpolation 

(l = lIr) 

1 1 

TRADKM Temperature 

interpolation radius 

of influence (km) 

36 36 

NUMTS Max number of 

stations for 

temperature 

interpolations 

5 5 

IAVET Spatially average 

temperature? (1 = 

yes) 

1 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 

IIr) 

2 2 

SIGMAP 

•

• 

Radius of influence 50 50 

•

• 

3-13 • Lafarge Building Materials Trinity Consultants

Refined BART Modeling Report 063701.0044 
• • • 

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

• Value Included in ••

• 

Value Trinity 

CALMET Variable Description CENRAP Protocol Notes 

for precipitation 

interpolation (km) 

Used 

•

• 

CUTP Minimum 0.01 

precipitation rate cut 

off(mmIhr) 

0.01 

••

• 
 ••

•
•

• •••••••

• 

3-14 Lafarge Building Materials Trinity Consultants 

Refmed BART Modeling Report 063701.0044 

•

• 
•

•
•
•

•

•

•
•

•

•
•
•

•
•

•
•
•

•

•

•
•

•

•
•
•
•

•
•
•

•

•
•

•

•
•
•
•


•

­

4. CALPUFF 


Trinity conducted a three-year CALPUFF analysis. The CALPUFF 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-ELIGIBLE 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 ofthe 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. Point ID 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 EP15 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 EP27 Coallcoak unloading 1974 

1102 EP28 Coallcoak drag conveyor, bucket elev., crusher 1974 

1I03 EP28 Coallcoak drag conveyor, bucket elev., crusher 1974 

4.2 SOURCE EMISSIONS DATA 

Trinity included emissions of802, NO", and PM10 in the model. For source-specific analyses, 

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

.. 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 


.. Potential to emit 


1 CENRAP BART Modeling Guidelines, p. 7-1. •

•

• 

4-1 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
• • • • • 

-- -- -- --

-- -- -- -- -- -- -- -- -- -- -- --

•

•

•

• CEMS data over the 2001-03 time period is not available (as the Tulsa Plant does not operate with 

CEMS units), so model emissions are primarily based on emission factors and maximum 24-hour 

production rates. 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 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 ofthese 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 

v. 

NO. PM10 S02 

Ibslhr Ibslhr Ibslhr 

EP04 0.089 3.857 0.001 

EP05 0.429 

EPl3 0.557 

EP15 576.00 17.143 960.00 

EP24 3.402 

EP26 0.429 

EP25 0.429 

EP17 4.607 

EP27 0.375 

EP28 0.321 

4.2.1 PMio EMISSIONS 

•

•

•

•

•

• 

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

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

• Managers' guidance to include the following:2

• 

3 

•

• 


• Coarse particulate matter (PMc, considered PM25• l0) 


• Fine particulate matter (PMF, considered PM<2.s) 

• 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 of primarily coarse 

(70%) and fine (30%). A summary of speciated PM 10 emissions is presented in Table 4-3. Full 

• 
 details of PM10 speciation can be found in Appendix A. 


2 http://www2.nature.nps.gov/airlPennits/ectlectCementKiln.cfin 


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


•

•

•

•

• 

4-2 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
• • • • • • • • • 

• • • 

.•
•

•

•
•

•

•

• 


4.2.2

-

Unit ID 

EP04 

EP05 

EP13 

EP15 

EP24 

EP26 

EP25 

EP17 

EP27 

EP28 

TABLE 4-3. SPECIATED PMIO EMISSIONS 

PM (coarse) PM (fine) S04 SOA EC 

Ibslbr Ibslbr Ibslhr Ibslbr Ibslhr 

0.964 2.893 -- -- -­ 0.300 0.129 -- -- -­ 0.390 0.167 -- -- -­ 1.149 1.320 12.91 1.766 -­ 2.381 1.021 -- -- -­ 0.300 0.129 -- -- -­

0.300 0.129 -- -- -­ 3.225 1.382 -- -- -­ 0.263 0.113 -- -- -­ 0.225 0.096 -- -- -­ • 

STACK PARAMETERS 

•

•

• 

Table 4-4 provides a summary ofsource exhaust parameters. 

TABLE 4-4. POINT SOURCE STACK PARAMETERS 

•

•

•

•

•

• 

Emission Latitude Longitude Base Exit Stack Stack Exit 

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

(m) (ftls) (feet) (feet) (Fahrenheit) 

EP04 36.194777 95.811055 267 11.94 100 4 180 

EP05 36.19475 95.810722 250 8.29 55 1.6 70 

EP13 36.195472 95.811916 259 9.55 97 1.7 200 

EPl5 36.195833 95.811777 234 17.83 60 6.9 500 

EP24 36.19425 95.811222 271 10.53 100 4 200 

EP26 36.194166 95.811166 260 33.16 64 0.8 70 

EP25 36.194166 95.811166 258 33.16 55 0.8 200 

EP17 36.19375 95.811805 240 3.56 53 8.0 200 

EP27 36.194027 95.811027 230 4.64 7 2.0 70 

EP28 36.194138 95.811333 245 4.91 135 1.8 70 

•

•

•

•

• Base elevations are based upon visual inspection of USGS topographic maps, while source locations 

are 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 were downloaded from the NPS we bsite, converted into the LCC NAD27 projection, 

4-3

•

Lafarge Building Materials Trinity Consultants •

• 

Refined BART Modeling Report 063701.0044 

•

• 
• • 

•
•

•

.. 


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 

CanBY C{i!tkWildomess

aD ReeeplCfS 

USFS BQlJnl3arl' SOI.Hte: 

FS NaticM! Co\l~rage tl.:l: 

NR 'S - ALP group Corvalli'9. OR 

August 22,2003 

.

• 

•

•

•

•

•

•

•

•

• .. .. .. .. '" 

.. ,. .. .. .. .. .. 

• 

• ' 

.. .. • .. ~ ..

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

L" .. · .. .. 

FIGURE 4-2. HERCULES-GLADES WILDERNESS RECEPTOR LOCATIONS 

•

•

•

•

•

• ~ 

.. •

•

•

•

• .. 

•

•

• 
 .. .. .. 


.. .. •

• .. .. .. .. .. 

.. •

•

•

.

•

' 

.. Hercules-Glades W~d.rness 

60 Receptors

-

USFS Boundary Source: 

FS National Coverage ule: 

NIlIS - I4..P grQup 

.. COIII.llis. OR 

August 22.2003 

-

• 

4-4 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
• • • 

• • 

• • 

.,

•

•


• 


•

•
•

•

•
•

•
• 


•

•

•
•

•

•
•

•

.'
•

•

•
•

•
•

• 


•

•
•

•
•
•

•
•

• 

FIGURE 4-3. UPPER BUFFALO WILDERNESS RECEPTOR LOCATIONS 

.. 

Upper Buffalo Wilderness 

72 Re ce ptors 

USFS Boundary Source: 

FS National Coverage file: 

NRIS - ALP group 

Corvallis, OR 

Aug us! 22,2003 

4-5 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

• • • • 

•

•
•

•

•

•
•

•
•

• 


• 

••

• 

•

• 

.. 


•

•

•
•

•

•
•

•
•

•

•
•

•
•

•

•
•

•
• 


•

•


•

Ct 

FIGURE 4-4. WICHITA MOUNTAINS WILDERl'JESS RECEPTOR LOCATIONS 

Wichita MountainsWilderness ;;.ithin Wichita MountainsNWR 

59 Receptors 

FWS Wilderness Boundary S.urce: 

Region 2 GIS 

500 Gold Ave., SW 

PO S•• 1306 

Albuquerque, NM 87103 

fWS ;JWR Boundary Source: 

htlp:llwwwJws.govldatalr2gislboundary.hlml 

Septem ber 18, 2003

· . . . . . . · .. .,. . ,. . .. 

·1 

· 

• 

.! 

.. :-:l 

+ 

4.4 BACKGROUND OZONE 

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

S02 and NO" to sulfates (SOJ 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 de~eloping background ozone estimates from 

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

conservative value of40 ppb (as recommended in CENRAP protocol Appendix) in lieu ofthe 

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 toCALPUFF. Since the data from CENRAP's CMAQ and CAMx simulations 

is not readily available, Trinity used a conservative monthly background concentration of 3 ppb. This 

4-6
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

• • • • • • 

• • • 

• • • • • 

• • • • • • • • • • • 

• • • • 

.. 


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 •

• 


• 

•

•

4' 


-


•

• 


•

•

• 


• 

Table 4-4 provides a listing ofthe CALPUFF parameters that Trinity used 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 

METRUN All model periods 

in met files will be 

run 

IBYR Starting year 

IBMO Starting month 

IBDY Starting day 

IBHR Starting hour 

XBTZ Base time zone (6 

=CST) 

IRLG Length of run 

NSPEC Number of 

MESOPUFF II 

chemical species 

NSE Number of 

chemical species 

to be emitted 

ITEST Program is 

executed after 

SETUP phase 

MRESTART Do not read or 

write a restart file 

during run 

NRESPD File written only at 

last period 

METFM CALMET binary 

file 

(CALMET.MET) 

AVET Averaging time in 

minutes 

PGTlME PG Averaging 

Value Included in 

CENRAP Protocol Value Trinity Used 

0 0 

2001 Appropriate met year 

1 I 

I I 

1 1 

6 6 

8760 8760 

10 10 

8 6 

2 2 

0 0 

0 0 

1 1 

60 60 

60 60 

Notes 

Years 2001, 2002, 

2003 

S02, S04, NOx, 

PMC, PMF, SOA 

4-7 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

• • 

• • • • • • • • • • .

' .

•
' •••it••

• 
 ••.'•
•••••

•
••••

• 
 •

• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

time in minutes 

MGAUSS Gaussian 

distribution used 

in near field 

1 1 

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 Stack tip 

downwash 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 MESOPUFF II 

chemical 

parameterization 

scheme 

1 1 

MAQCHEM Aqueous phase 

transformation not 

modeled 

0 0 

MWET Wet removal 

modeled 

1 1 

MDRY Dry deposition 

modeled 

1 1 

MDlSP PG dispersion 

coefficients 

3 ,., 

.J 

MTURBVW Use both a v and 

a w from 

PROFILE.DAT to 

compute a 

yand a, 

(n/a) 

3 ,., 

.J 

MDlSP2 

•

• 

PG dispersion 3 

3

••

• 

4-8 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 
• • • • • 

• • • 

.' 

•

•

•

•

•

• 

•

•

•

•(e'

•

• 

'..'•

• 

•

• 

•

•

• 

'.

• 

•

•

•

• 

•

•

•

• 

4-9 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

coefficients 

MROUGH PG cry and crz not 

adjusted for 

roughness 

0 0 

MPARTL No partial plume 

penetration of 

elevated inversion 

1 1 

MTINV Strength of 

temperature 

inversion 

computed from 

default gradients 

0 0 

MPDF PDF not used for 

dispersion under 

convective 

conditions 

0 0 

MSGTIBL Sub-grid TIBL 

module not used 

for shoreline 

0 0 

MBCON Boundary 

concentration 

conditions not 

modeled 

0 0 

MFOG Do not configure 

for FOG model 

output 

0 0 

MREG Technical options 

must conform to 

USEPA Long 

Range Transport 

(LRT) guidance 

1 1 

CSPEC CENRAP Lafarge 

Output 

Group 

Species Modeled Emitted 

Dry 

Deposition 

Output 

Group 

Species Modeled Emitted 

Dry 

Deposition 

S02 1 1 1 S02 1 1 1 

S04 1 0 2 S04 1 1 2 

NOX 1 1 1 NOX 1 1 1 

HN03 1 0 1 HN03 1 0 1 

N03 1 0 2 N03 1 0 2 

NH3 1 1 1 NH3 1 0 1 

PMC 1 1 2 PMC 1 1 2 

PMF 1 1 2 PMF 1 1 2 

Lafarge Building Materials Trinity Consultants 

•

•

• 

Refined BART Modeling Report 063701.0044 
• • • 

•

•

•

•


•

.~ 

• 


•

•

• 

•

• 

• 

«I 


•
• 


•

•
•

•

•
•

•

•
•

•

•

•

•

•
•

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

~ Map projection 

12 EC I 1 0 12 

-+ 12 SOA 1 1 I 12 

PMAP UTM LCC 

NX Number of X grid 

cells in 

meteorological 

grid 

66 256 Based on refined grid 

parameters as 

specified in 

CALMETruns 

NY Number of Y grid 

cells in 

meteorological 

grid 

66 155 Based on refined grid 

parameters as 

specified in 

CALMETruns 

NZ Number of vertical 

layers in 

meteorological 

grid 

10 10 

DGRIDKM Grid spacing (km) 6 3.0 Refined grid size 

ZFACE Cell face heights 

in meteorological 

grid (m) 

0,20, 40, 80, 160, 

320,640,1200,2000, 

I 3000,4000 

0,20,40,80,160, 

320,640, 1200, 

2000,3000,4000 

XORIGKM Reference X 

coordinate for SW 

comer ofgrid cell 

of meteorological 

grid (km) 

5 0 Appropriate for 

domain (refined grid) 

YORIGKM Reference Y 

coordinate for SW 

comer of grid cell 

of meteorological 

grid (km) 

3327 0 Appropriate for 

domain (refined grid) 

IUTMZN UTM zone of 

coordinates 

(NAD83)(use 

only if 

PMAP=UTM) 

12 14 N/A;LCC 

coordinates used 

IBCOMP X index of lower 

left comer ofthe 

computational grid 

1 1 

JBCOMP Y index oflower 

left comer ofthe 

computational 

grids 

1 1 

IECOMP X index ofupper 

right comer ofthe 

66 256 Appropriate for 

domain (refined grid) 

Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

•

•

•

•

•

•

•

• 
• • • 

• • • • 

•

•

•

•

.' 


•

•

•

• 


•

•

•

•

•

•

•
•

•

•

•

•

•
•

•

•

•

•

•

•

•

•

• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

computational grid 

JECOMP Y index ofupper 

right comer ofthe 

computational grid 

66 155 Appropriate for 

domain (refined grid) 

LSAMP Sampling grid is 

not used 

F IF 

IBSAMP X index of lower 

left comer of 

sampling grid 

1 1 

JBSAMP Y index oflower 

left comer of 

sampling !!rid 

1 I 

IESAMP X index of upper 

right comer of 

sampling grid 

66 256 Appropriate for 

domain (refined grid) 

JESAMP Y index ofupper 

right comer of 

sampling grid 

66 155 Appropriate for 

domain (refined grid) 

MESHDN Nesting factor of 

sampling grid 

I 1 

ICON Output file 

CONC.DAT 

containing 

concentrations is 

created 

I 1 

IDRY Output file 

DFLX.DAT 

containing dry 

fluxes is created 

I 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 boundaries 

0 0 

•

•

• 

4-11 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 
•••••••

•
•••••••

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

•
••••••

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

IMBAL Mass balances for 

each species not 

reported hourly 

0 0 

ICPRT Print concentration 

fields to output list 

file 

1 1 

IDPRT Do not print dry 

flux fields to 

output list file 

0 0 

IWPRT Do not print wet 

flux fields to 

output list fil e 

0 0 

ICFRQ Concentration 

fields are printed 

to output list file 

every hour 

1 1 

IDFRQ Dry flux fields are 

printed to output 

list file every 1 

hour 

1 1 

IWFRQ Wet flux fields are 

printed to output 

list file every 1 

hour 

1 1 

IPRTU Units for line 

printer output are in 

g/mJ for 

concentration and 

g/m2/s for 

deposition 

3 3 

IMESG Messages tracking 

the progress ofrun 

written to screen 

2 2 

LDEBUG Logical value for 

debug output 

F F 

IPFDEB First puffto track 1 1 y 

NPFDEB Number of puffs 1 1 

NNI Meteorological 

period to start output 

1 1 

NN2 Meteorological 

period to end 

output 

10 10 

4-12 

••

• Lafarge Building Materials Trinity Consultants • Refined BART Modeling Report 063701.0044 •

•
• 
•••••••••

•
••

•
•••

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

• 

CALPUFF 

Variable 

• 

Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

NHILL Number ofterrain 

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 ofCTDM 

system relative to 

CALPUFF 

coordinate system 

(km) 

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 

•••

• 

Henry's Law 

coefficient 

3.5 3.5 

4-13 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

• 
CALPUFF 

Variable Description 

HN03 Diffusivity 

Alpha star 

Reactivity 

Mesophyll 

resistance 

Henry's Law 

coefficient 

S04-2 Geomatric mass 

mean diameter of 

S04-2 (flm) 

N03­ Geometric mass 

mean diameter of 

N03- ().tm) 

PMC Geometric mass 

mean diameter of 

PMC (flm) 

PMF Geometric mass 

mean diameter of 

PMF (flm) 

EC Geometric mass 

mean diameter of 

EC (flm) 

SOA Geometric mass 

mean diameter of 

SOA (flm) 

RCUTR Reference cuticle 

resistance (stcm) 

RGR Reference ground 

resistance (stcm) 

REACTR Reference 

pollutant reactivity 

NINT Number of particle 

size intervals for 

effective particle 

deposition velocity 

IVEG Vegetation in nonirregated areas is 

active and 

unstressed 

S02 Scavenging 

coefficient for 

liquid precipitation 

(S'I) 

Value Included in 

CENRAP Protocol 

0.1628 

1 

18 

0 

8.e-8 

0.48 

0.48 

6 

0.48 

0.48 

0.48 

30 

10 

8 

9 

1 

3.21E-05 

Value Trinity Used 

0.1628 

1 

18 

0 

8.0e-8 

0.48 

0.48 

6 

0.48 

0.48 

0.48 

30 

10 

8 

9 

1 

3.21E-05 

Notes 

4-14
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 


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

• 
• • • • 

•
••

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

•
•••••

•
•••••

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

Scavenging 

coefficient for 

frozen 

precipitation (5.

1

) 

0 0 

S04-2 Scavenging 

coefficient for 

liquid precipitation 

(5.

1) 

l.OE-04 1.0E-04 

Scavenging 

coefficient for 

frozen 

precipitation (5. 

1

) 

3.0E-05 3.0E-05 

HN03 Scavenging 

coefficient for 

liquid precipitation 

(5.

1

) 

6.0E-05 6.0E-05 

Scavenging 

coefficient for 

frozen 

precipitation (5.

1

) 

0 0 

N03­ Scavenging 

coefficient for 

liquid precipitation 

(5.

1

) 

l.OE-04 l.OE-04 

Scavenging 

coefficient for 

frozen 

precipitation (S'l) 

3.0E-05 3.0E-05 

NH3 Scavenging 

coefficient for 

liquid precipitation 

(5. 

1

) 

8.0E-05 NA 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

0 NA 

PMC Scavenging 

coefficient for 

liquid precipitation 

(5.

1) 

1.0E-4 1.0E-4 

••

• 

Scavenging 

coefficient for 

frozen 

3.0E-05 3.0E-05 

•

• 


4-15
 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 
• • • • • • • • •••••••

•
••

•
•

•
•••••

•
•
•

•
••

•
•

•
•

•
•••

•
•
•

• 


CALPUFF 

Variable Description 

precipitation (s") 

PMF Scavenging 

coefficient for 

liquid precipitation 

(s") 

Scavenging 

coefficient for 

frozen 

precipitation (s") 

EC Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

Scavenging 

coefficient for 

frozen 

precipitation (S·I) 

OC Scavenging 

coefficient for 

liquid precipitation 

(S·I) 

Scavenging 

coefficient for 

frozen 

precipitation (s") 

MOZ Read ozone 

background 

concentrations 

from ozone.dat file 

(measured values) 

BCK03 Background ozone 

concentration 

(ppb) 

BCKNH3 Background 

ammonia 

concentration 

(ppb) 

RNITEI Nighttime S02 

loss rate is o/olhour 

Nightti"'" NOX 

loss rate is %/hour 

Nighttime HN03 

loss rate is o/olhour 

Value Included in 

CENRAP Protocol 

I.OE-05 

3.0E-05 

1.0E-04 

3.0E-05 

I.OE-04 

3.0E-05 

1 

12*40 

12*3 

0.2 

2 

2 

Value Trinity Used 

I.OE-04 

3.0E-05 

I.OE-04 

3.0E-05 

I.OE-04 

3.0E-05 

1 

12*40 

12*3 

0.2 

2 

2 

Notes 

4-16 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

• • • 

• • • 

•
•
• 


•

•
•

•
•

•

•

•

•

•

•

•

•
•

•

•

•
•

•
•
•

•

•

•

•
•

•

•

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

MH202 Background H202 

concentrations 

I 0 Need to choose 0 in 

order to use monthly 

background value 

BCKH202 Background 

monthly H202 

concentrations 

1 12*1 

BCKPMF Fine particulate 

concentration for 

SOAoption 

(llg/m3) 

1 1 

OFRAC Organic fraction of 

fine 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 

SYTDEP Horizontal size of 

a puff in meters 

beyond which the 

time dependant 

dispersion 

equation ofHeffter 

is used 

550 550 

MHFTSZ Do not use Heffter 

formulas for sigma 

z 

0 0 

JSUP Stability class used 

to determine 

dispersion rates for 

puffs above 

boundary layer 

5 5 

CONK I Vertical dispersion 

constant for stable 

conditions 

0.01 om 

CONK2 Vertical dispersion 

constant for 

neutral/stable 

conditions 

0.1 0.1 

TBD Use ISC transition 

point for 

determining the 

transition point 

between the 

Schulman-Scire to 

0.5 0.5 

•

•

•

• 

4-17 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
• • • • • • • • 

.,
 •••••••••••

•
•••••••••

• 
 ••

•
•

•
•••••••

• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

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 

ILANDUIN 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 

. . 

ANEMHT Anemometer 

height in meters 

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 of slug or 

puffin 

meteorological 

10 10 

4-18 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

••••••••

•
•••

•
•

•
•

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

•
•••••••

• 


CALPUFF 

Variable 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

grid units 

one sampi

[)e""~ MXNEW Maximum number 

of puffs or slugs 

released from one 

source during one 

time step 

60 60 

MXSAM Maximum number 

of sampling steps 

during one time 

step for a puff or 

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 

1 1 

SZMIN Minimum sigma z 

in meters for a 

new puffor slug 

1 1 

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 Maximummi 

height (m) 

3000 3000 

XMlNZI Minimummi' 

height(m) 

20 20 

SL2PF Maxmium cry/puff 

length 

10 10 

PLXO Wind speed 

power-law 

0.07,0.07,0.10,0.15, 

0.35; 0.55 

0.07, 0.Q7, 0.10, 

0.15,0.35,0.55 

4-19 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•••••

• 
• • • • • 

• • • • 

•

•

•

•

•

•

•

•

•

•

•

• 

•

•

•

•

•

•

• 

•

•

•

•

•

•

•
•

•

• 


CALPUFF 

Variable 

• 

" 

Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

exponents 

WSCAT Upper bounds of 

1

st 5 wind speed 

classes 

1.54,3.09,5.14,8.23, 

10.80 

1.54,3.09,5.14, 

8.23,10.80 

PTGO Potential temp 

gradients PG E & 

F (deg/km) 

0.020, 0.035 0.020,0.035 

CDIV Divergence 

criterion for dw/dz 

(lIs) 

0.01 0.01 

PPC Plume path 

coefficients (only 

if MCTADJ = 3) 

. 0.5,0.5,0.5,0.5, 

0.35,0.35 

0.5,0.5,0.5,0.5, 

0.35,0.35 

NSPLIT Number of puffs 

when split 

3 3 

IRESPLIT Hours when puff 

is eligible to split 

1900 Hour 17 Should be by hour of 

day - 1900 is hour 17 

ZISPLIT Previous hours 

minimum mixing 

height, m 

100 100 

ROLDMAX Previous max 

mixing 

height/current 

height ratio, must 

be less than this 

value to allow puff 

to split 

0.25 0.25 

NSPLITH Number of puffs 

resulting from a 

split 

5 5 

SYSPLITH Minimum sigma-y 

of puff before it 

may split 

1.0 1.0 

SHSPLITH Minimum puff 

elongation rate 

from wind shear 

before puff may 

split 

2.0 2.0 

CNSPLITH Minimum species 

concentration 

before a pluff may 

split 

1.0E-07 1.0E-07 

EPSSLUG Criterion for 1.0E-04 1.0E-04 

• 

4-20 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

•

•

• 
• • • • 

•
••

•
••••

• ••

•
••

•
•••••••

•
•••

•
•••

•
•••

• 

CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

SLUG sampling 

EPSAREA Criterion for area 

source integration 

1.0E-06 I.OE-06 

DSRISE Trajectory step 

length for 

numerical site 

algorithm 

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 = Ib/hr) 

NSPTl Number of sourcespecies 

combinations with 

variable emissions 

scaling factors 

--

NPT2 Number of point 

sources with 

variable emission 

parameters 

provided in 

external 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 

rates are glm2/s 

1 1 

NSARI Number of source 

--

•

• 

4-21 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

• •

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

•
••••

• 
 ••••

•
•

•
••••••••

•
•

•
•

•
•

•
••••

• 

CALPUFF 

Variable Description 

species 

combinations with 

variable emissions 

scaling factors 

NAR2 Number of 

buoyant polygon 

area sources with 

variable location 

and emission 

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 gls 

NSLNl Number ofsourcespecies 

combinations with 

variable emissions 

scaling factors 

MXNSEG Maximum number 

of segments used 

to model each line 

NLRlSE 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) 

WML Average line 

source width (m) 

Value Included in 

CENRAP Protocol 

-----------

--

. ---

. ----

Value Trinity Used Notes 

4-22
 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

••••

•
••

•
••

•
•••

•
••••••

•
••••••

•
••••••

• 


CALPUFF 

Variable Description 

Value Included in 

CENRAP Protocol V Notes 

DXL Average 

separation 

between buildings 

(m) 

· · 

FPRlMEL Average buoyancy 

er (m4/s3) 

· · 

NYLl of volume · 0 None modeled 

IVLU Units for volume 

source emission 

rates in 

grams/second 

· 3 . 

NSVLl Number ofsource· 

species 

combinations with 

variable emissions 

scaling factors 

· · 

rGRDVL Gridded volume 

source data is not 

used 

--

VEFFHT Effective height of 

emissions (m) 

-

-

VSIGYI Initial sigma-y 

value 

-

-

VSIGZI Initial signia-z 

value 

-

-

NREC 

•••

• 

Number of nongridded receptors 

5630 291 Receptor data 

provided in 

Appendix 

••

• 

4-23 Lafarge Building Materials Trinity Consultants 

RefIned BART Modeling Report 063701.0044 

• 
5. CALPOST 


Trinity conducted 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 utilized 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. 

= 3*f(RH)*[(NH4)2S04] + 3*f(RH)*[NH4N03] + 4*[OC] + l*[PM f b ] ext 

+ 0.6 *[PMe ] + 10*[EC] + bRay 

The algorithm was 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 used CALPOST Method 6, which calculates hourly light extinction impacts for the source and 

background using monthly average relative humidity adjustment factors. Trinity used monthly Class 

I area-specific relative humidity adjustment factors based on the centroid of the Class I area as 

included in Table A-3 of EPA's Guidance/or 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 Jul 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 Report 06370l.0044 

•

• 
•
• 


5.3 NATURAL BACKGROUND 

•

•

•

•

•

• For Caney Creek and the Wichita Mountains, Trinity used default west natural background values. 

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

•

• 

The annual average concentrations are included in Table 2-1 ofEPA's Guidance/or 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 •

•

•

•
•

•

•

•

• • 5.4 EvALUATING BART-EXEMPTION

•

•

• 

Trinity compared the high ofthe 2001 through 2003 daily ~dv values output by CALPOST to a 

contribution threshold of 0.5 ~dv. Ifthe 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 highest ofthe ~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 of the CALPOST parameters that Trinity used 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 

Component 

Average Annual 

Natural Background 

West (j.lglm3 

) 

0.12 

Average Annual 

Natural Background 

East Olglm3 

) 

Ammonium Sulfate 0.23

Ammonium Nitrate 0.10 0.10 

Organic Carbon Mass 0.47 lAO 

Elemental Carbon 0.02 0.02 

Soil 0.50 0.50 

Coarse Mass 3.0 3.0 

CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

ISYR Starting year 2001 Appropriate met 

year 

Years 200 I, 2002, 

2003 

ISMO Starting month I I 

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 Report 063701.0044 •

•

• 
•
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 


CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol .Value Trinity Used Notes 

ISHR Starting hour 0 1 All CALPUFF 

periods will be 

included 

NPER Number of periods to 

process 

8760 8760 

NREP Process every hour of 

data? 1 yes 

1 1 

ASPEC 

ILAYER 

Process species for 

visibility 

VISIB VlSIB 

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? 

Process complex 

terrain receptors? 

T T 

LCT F F 

LORING Report receptor ring 

results? 

F F 

NDRECP Select all discrete 

receptors 

-1 Varies As appropriate for 

Class I area being 

analyzed 

IBGRID 

mGRID 

X index ofLL comer 

of receptor grid 

-1 -1 

Y index ofLL comer 

of receptor grid 

-1 -1 

IEGRID X index of UR corner 

of receptor grid 

-1 -1 

JEGRID Y index of UR corner 

of receptor grid 

-1 -1 

NGONOFF Number of gridded 

receptor rows 

0 0 

NGXRECP Exclude specific 

gridded receptors, Yes 

=0 

0 0 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

5-3 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•
•

•
•

•

•
•
•

•
•

•

•

•
•

•
•

•

•

• 


CALPOST 

Variable Description 

RHMAX Maximum RH% used 

in particle growth 

curve 

LVS04 Compute light 

extinction for sulfate? 

LVN03 Compute light 

extinction for nitrate? 

LVOC Compute light 

extinction for organic 

carbon? 

LVPMC Compute light 

extinction for coarse 

particles? 

LVPMF Compute light 

extinction for fme 

particles? 

LVEC Compute light 

extinction for 

elemental carbon? 

LVBK Include background in 

extinction calculation? 

SPECPMC Coarse particulate 

species 

SPECPMF Fine particulate 

species 

EEPMC Extinction efficiency 

for coarse particulates 

EEPMF Extinction efficiency 

for fine particles? 

EEPMCBCK Extinction efficiency 

for coarse part. 

Background 

EES04 Extinction efficiency 

for ammonium sulfate 

EEN03 Extinction efficiency 

for ammonium nitrate 

EEOC Extinction efficiency 

for organic carbon 

EESOIL Extinction efficiency 

for soil 

EEEC Extinction efficiency 

for elemental carbon 

Value Included in 

CENRAP Protocol 

95 

T 

T 

T 

T 

T 

T 

T 

PMC 

PMIO 

0.6 

1.0 

0.6 

3.0 

3.0 

4.0 

1.0 

10.0 

Value Trinity Used 

95 

T 

T 

T 

T 

T 

T 

T 

PMC 

lPMF 

0.6 

1.0 

0.6 

3.0 

3.0 

4.0 

1.0 

10.0 

Notes 

Notation difference 

•

•

•

•

•

•

• 

5-4 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

• 
•
•

•

•

•

•

•

•

•

•

•

•

•

•
•

•

•

•

•

•

•

•

•

• 


CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used 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 % of particles 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: 

304,3.1,2.9,3.0, 

3.6,3.6,304,304, 

3.6, 3.5, 304, 3.5 

Hercules-Glades: 

3.2,2.9,2.7,2.7, 

3.3,3.3,3.3,3.3, 

304,3.1,3.1,3.3 

Upper Buffalo: 

3.3,3.0,2.7,2.8, 

304,304,3.4,304, 

3.6, 3.3, 3.2, 3.3 

Wichita Mountains: 

2.7,2.6,2.4,204, 

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 

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 lAO 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.02 Used for HerculesGlades and Upper 

Buffalo 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

5-5 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

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

• • • • • • • • • 

•

•

•

•
•

•
•

•
•

•

•
•

•

•

•

•
•

•

•
•
•

•

•

• 


CALPOST 

Variable Description 

BKS04 Background sulfate 

extinction coeff- west 

BKN03 Background nitrate 

extinction coeff- west 

BKPMC Background sulfate 

extinction coeff - west 

BKOC Background organic 

carbon extinction 

coeff-west 

BKSOIL Background soil 

extinction coeff- west 

BKEC Background elemental 

carbon extinction 

coeff- west 

BEXTRAY Extinction due to 

Rayleigh scattering 

(IIMm) 

LDOC Print documentation 

image? 

IPRTU 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 Used 

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 Report 063701.0044 

• • • • 

• • • 

• • • • • • • • 

•
•

•
•

•
•

• 


•

•
•

•

•

• 


•

• 


CALPOST 

Variable Description 

Value Included in 

CENRAP Protocol Value Trinity Used Notes 

THRESHl Averaging time 

threshold for 1 hr 

averages 

-1 -1 

THRESH3 Averaging time 

threshold for 3 hr 

averages 

-1 -1 

THRESH24 Averaging time 

threshold for 24 hr 

averages 

-1 -1 

THRESHN Averaging time 

threshold for NAVOhr averages 

-1 -1 

NDAY Accumulation period, 0

days 

0 

NCOUNT Number of 

exceedances allowed 

1 1 

LECHO Echo option F F 

LTIME Time series option F F 

LPLT Plot file option F F 

LORD Use grid format 

instead ofDATA 

format 

F F 

LDEBUO Output information for 

debugging? 

F F 

•

•

•

•

•

•

•

•

•

•

•

•

• 

5-7 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

-
• • • • 

• • • 

6. MODEL RESULTS 


A summary ofthe top 25 visibility-impairing days at each ofthe relevant Class I areas over the 

•

•

•

•

•

•

•

•

•

•

• 

modeled 2001-03 period is summari zed in Table 6-1 to 6-4 below. Shaded values indicate visibility 

impacts that exceed the 0.5 Lldv threshold for contribution to visibility impairment. As seen in these 

tables, the 98th percentile values indicate that the Tulsa Plant's contribution to visibility impairment is

•

• 

below the 0.5 dv threshold at each ofthe neighboring Class I areas. These results indicate that the 

Tulsa Plant therefore qualifies for exemption from BART. 

TABLE 6-1. VISIBLITY IMPACT MODEL RESULTS FROM CANEY CREEK 2001-03 

Rank Modeled Background Extinction by Species (%) Year Day 

Visibility 

____________~Im~p_a_ct_(~Ll_d_v~)1-~S~O~4L-__~N~O~3__~O~C~~E~C~~PM~C~~PMF 

"-peak 24-hr A v~rage 45.59 _ 53.9_~,----=.:0.:::.:21 0 • 2

•

• 

6.6_.5_8__ ..._.~3.25 O.O?~.. __..:..;.;:.____.....;;. 

68.81_--.:;.;=:::....__ 

4 

8 

•

•

•

•

• 

9 


10 
....._._-+--_..............."'..:...:::::........--_.. ---=~-... -....".".,.~.............-...".:..'-'-.....--....."..-....__.--'=-"_..........::0.:.:..:.._+._...::".".:...._._......"..:::..,'.......1 


11

•

• 

12 L........................~.:...: : ....:.................. +_._36_.3_2_...§~:~~_... Q.,~? ...__.... . o 0.04 0.14 .........--.. ~ ....­ ..,,---.---~~-, 

1-__..._......:;:.13:.....__.__1...__:::.:...:..:....:......__1...........:2=2:.:.:.2::::5_......2I:.~.~_ ~.. 0.04 0.11 


1_........___._....._....._:..'-....__........_.+ __....................".:.:.. :~.:.................. _1_.::::55:::.:.:... 1 - .._ 44:~J2 ..._Q:g§__
 ...!L__~~. 

1___.._.............:.::........................_ ..+··...···..··..·....·"..:...:.,~·"····....·..·..···_·1_.:...73:::.:.2=--_........-"'.,,.:."'--"........._Q.:.!L_...o___0._05_... 


1------___......:;:.::..______+ .._...__.::.:..:.:~___.~1...........:2::.::9::.::.8:.:::9___§9.4§..._._0_.2_9___0_ 0.06 


17 39.96 0.15 0.01 •••. " .._.__....w_~_""'·~."""'_~, ° 18 

20 

•

•

•

•

•

•

•

• 

0.04 

98th Percentile 47.42 o 0.04 0.2 

..=2.::..3____-I__-------'0"-'".1.:.:;2:.....__ .........:::66;::.:.::..54-'--._32_.7_4 __0_.3_3_....... ° 0.04 ....Q,,~L 2003 110 r

._11_._ 0.115 1--3=2:.:..1:..:1_....~~!.:<!.t.. _Q.._38_.__ !L___ o.!..!_._ 0.35 2001_..~§_

•

• 

25 0.104 64.25 35.5 0.12 0.01 0.12 2001 89 ° 

•

• 

6-1 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 

•

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

•

• 


TABLE 6-2. VISIBLITY IMPACT MODEL RESULTS FROM HERCULES-GLADES 2001-03 

•

•

•

•

•

•
•

•

•
•

•

•

•

•

•

•

•
•

•

•
•

•

•
•

•

•

•

• 


Rank Modeled Background Extinction by Species (%) Year Day 

Visibility 

Impact (Adv) 504 N03 OC EC PMC PMF 

Peak 24·hr Avera~I-__0:::..3;:c8:.-01_ 30.77 68.58 0.32 0 0.05 0.28 2001 362 

2 30.45 ....._---..__.._-_...... + .....--..--==:::=-...--+---=.::::...:..=-­-,....__69 ....._._0.26 .. 0 0·.. ___9·2f'L ....,_ 2003 351 .-.~"""-~,--.~-.-.-

3 0.232 85.49 14.2 0.14 0 0.02 0.15 2003 256 

0.191 32.92 2002 12 --""-'-' 

39.7 nn___60.01 """" __h 0.16 0 0.01 0.13 2001 22 "-'~-' -.~.---.".--....... 

64.13 ~~L__O.13 0 0.03 0.15 ?~_.. ?L. 

_ ...._.+.....:6:::.::6"-.1:..:::5__.._~~:!?1 __.__QJ_L____O___.. ...Q,QL_QJ~.. 2.9.92.... _.......___~L 

44.92 

1­__.__,._._...... _____ +_.___....~..:..::....:....__...._1__7:...:3"".2::..::6'-­ ......?f3:1?......._o'Jl.._ .. -.2_...Q,.9..L_ ..._9.J.;L_ ....2..oq? ___.?~?_. 

f____7..,-'-"--:::._......_+--=2:..:.7.:..:.0c::3_.......1?:3~......_O'.27._.__0 _ .. O.O§._ ...JL?L._......fQ9.L_.._~_ .. 

0.128 27.6 71.62 0.38 0 0.05 0.35 2002 47 - .. ....-.. ----~"--

1-.....................___.....:.............__...___ +...................c.:...:......::.._................... f...._...~Q:.!5.9.._......... 48. 86 ......_.....9..:.?.~._......_ ..._Q....... ....Q,g.f3..._.9...27.... _... _ 2003 ...............!l.?....._ 

l3 

14 

16 

17 

18 

19 

1--_---=2::..::0..___ 

21 

98th Percentile 

23 

24 

25 

41.77 57.85 

61.73 37.32 ."..~~......-~ 

54.45 45.18 

41.95 57.61 

52.86 46.71 - ..,,~,~,,-..---. 

38.76 

52.28 

:.:=-:_-----'3::..::9.94 

5.27 

0.16 0 0.02 

....__0.54 ...._-_0 .•._",_ 0.06 . 

0.2 0 0.01 '''''-­

0.21 0 0.02 

0.15 0 0.03 

0 

0.22 o 

0.13 o 

0.46 o 

0.23 o 0.08 0.26 

6-2 Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Report 063701.0044 

• • • • • • • • • 

• • • 

• • • • 

• • • • • • • • • 

TABLE 6-3. VISIBLITY IMPACT MODEL RESULTS FROM UPPER BUFFALO 20Ot-03 

Rank Modeled Background Extinction by Species (%) Year Day 

•

•

•

• 


•

•
•

•

• 

•

• 98th Percentile 

0.114 40.74 58.94 0.16 0 0.01 0.14 

24 0.113 58.63 40.9 0.21 0 0.02 0.24 2001 .............. 


•

25 0.113 73.93 •

•

•

•

• 

6-3 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Report 063701.0044 •

• 
•• 

•• 

•

•

•

•

•

•

• TABLE 6-4. VISIBLITY IM:PACT MODEL RESULTS FROM WICHITA MOUNTAINS 2001-03 


Background Extinction by Species (%) 
 Year Day 

Visibility Impact 1----------------------------1 

(L\dv) 

Rank Modeled 

504 N03 OC EC PMC PMF 

28.14 0.15 0 0.03 0.15 

•

•

•

•

• 

19 0.34 0 0.1 0.34 

66.34 0.27 0 0.07 0.26 2002 306 
 -,.~~-.. 

,42.27 57.04 0.31 0 (l.07 2003 43 


5 0.24 78.29 .. 21.38 .~_. 0 0.04 2003 m"".....""_•• ~.,.,".._''',',' ...J§L

r-_____6~.. _ ............_ ......_.~0.2=:2=8_..__....+-~~~.~3~1___4~5~.5~5__~0~.0~7___~0~__~0.~0~1 0.07 2003 59
 M_·'W~..___.,.. 

1-- ..__...7.:.__..__..._.1 ....................0:"..=2:""'26-:.........._ .....--l... ~0=.0:.::2'--_4.:..:9:.:::.5:..:..7_-=0.=.2.3_.()._-.!l:OL_0J..i......_.2001..~308 ..... 

1-___..._.8.:::._____+.._ ....._....:::.:0..:::2.:17.:._..__+......:-:3:-::::8.29"----"6:..:.1:.::..2:::..9_...:0:.:.:.1:.::9.__.::0..__.J1.:.0~. 0.18 ~.__.~_ .. 

9 0.213 52.47 47.07 0.23 0 0.05 0.19 2003 145 


•

•

•

•

•

• 

1--__....c1;..;.4 __----i___... __ QJ 38......_ .._ .....1_9;:..;1:..:,:.7..:..1_--'-7.:.::.94-'.- --=0.:..;..1..:..4_--=0__.....:0:.:.;.0:....:.4.__.....::.:0 ...:..;17_:--=-=_ 2002 149 


1--__.....::. 1-=---.5_._-'._..._ ..._.__ 0-::..:.1..:.;12=6 .. ...._.__..t_._6_1"Q4 ___..::3c:::8.c:::.,61'--_-=0 ..:.:15'---.-._ ....Q __0,QIL_..;;.;0.:.:.;15,,-­

16 0.126 85.29 14.22 0.21 0 0.06 0.22 

17 0.12 4§l.:.t? .......__~.9_.B_5_0_.2§_ 0 0.06 0.24 
 .............. --,.,-.'''~~'''''.. ~". 

18 0.12 BO.19 19.39 0.19 0 0.04 0.21 I----~--·--I······················..··:.,··~=·-·-··-·i--·~~-~~ ......­

19 0.119 42.63 56.9 0.22 0 0.02 0.22 

20 0.11 2.B2 0.34 0.07 0.37 ._-'-,,--" ° 21 0.102 95.05 4.09 0.38 0 0.08 0.4 

98th Percentile 0.102 50.45 49.12 0.2 0 0.02 0.21 

23 0.102 26.84 73.06 0.05 0 0.Q1 0.04 

__._..._~4,_""""_"_,___""",,,,,,_, 37.24 0.24 

25 59.38 0.45 o 0.19 0.44 

•

•

•

•

•

•

•

•

•

•

•

•

•

• 

,.•

•

•

•

• 

6-4
 • Lafarge Building Materials

2002 156 


2003 279 


2002 357 


2003 278


N_~D __ 

2001 16 


2003 223 


2001 223
 ,"-­

2001 78 


2002 123 
 '~_ ••M •• •• _ ••• _. 

2001 268 


•

• 

Trinity, Consultants 

Refined BART Modeling Report 063701.0044 

-

• 
•• 

ApPENDIX A - EMISSION DATA 

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

•
••••••

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

• Lafarge Building Materials Trinity Consultants • Refined BART Modeling Protocol 063701.0044 •

•
-
• ••• ••••

,••:1 ..~ lit•••••••••·•••• " .\ 

Printed 1013012006 

Year 


Installed 


1964 

1964 

1964 

1964 

1964 

1964 

1964 

1964 

1964/1997 

1974 

1974 

1974 

Was this Unit in 

Existence on 

August 7, 1977? 

(Y or N) 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Was this Unit in I Potential to Emit (tons/year) 

Operation after 

August 7, 1962? 

(Y or N) I PM IO I NO. I SOl I VOC I NH) 

Y 28.90 I 2.70 I 0.02 I 0.20 

I I I Y 3.20 Y 

Y 4.20 

Y 128.4011261.0012102.001 5.90 

I I Y 25.50 I

Y 3.20 

Y 3.20 I I I I I 

I I I I Y 34.50 I

Y I 7.00 

Y 2.40 

Y 

240.50 [1263.7012102.021 6.10 I 0 

Emission P . t 

L I P' OlD EUG Oint ID I(UnitID) ~ 

2A I I I

11 I I 7 Coallcoak unloading 

Unit Description 

(Include MMBtulhr 

if applicable) 

Which of the 

26 Categories 

does this Unit 

fall into?A 

I 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

EP28 Coallcoak drag cOlweyor, bucket elevators 3 

EP28 Coallcoak crusher 3 

Total 1 1 

A. See attachc<nist of categories 

Lafarge Building Materials Using Title V Trinity Consultants 


BART Modeling A-1 063701.0044 


,'•.•I.;<i 
•••••••••••••••••••••••••••••••••••••••••••• 

Printed 1013012006 

Source Emissions 

Location Point III Unit Description , Constr.Datel Heat Duty 

(MMBtulhr) 

Process Rate 

TPH 

Hours of Operation 

(hrs/day) 

Emission Factors 

PM10 

(Ilr/ad) 

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

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

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

Tulsa EP26 Finish mill #2 1964 -­ 39 24 0.025 

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

Tulsa EPI7 Clinker cooler #2 1964/1997 -­ 96 24 0.025 

Tulsa EP27 Coalkook unloading 1974 -­ ISO 24 0.025 

Tulsa EP28 Coalkook drag conveyor, bucket elevators, crusher 1974 -­ ISO 24 0.025 

Construction datcs and sourcc parameters as listed in Permit 98-174-TV. Construction of EP14 (odlcr kiln) commcnced on 1961 

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

Source Emissions 

Location Point III 

Tulsa EPOS 

Unit Description 

Raw mill pump hopper 

N02 

Ibs/hr 

-­

PM10 

Ibslbr 

0.429 

Emission Data Input 

PM (coarse)2 PM (fine)l 

Ibslhr Ibslbr 

0.300 0.129 

S02 

Ibs/hr 

-­ Tulsa EPI3 Kiln #2 feed system 0.557 0.390 0.167 -­ Tulsa EP24 #2 bucket elevator and separator 3.402 2.381 1.021 

Tulsa EP26 Finish mill #2 0.429 0.300 0.129 -­ Tulsa EP2S Finish mill #2 feed belt 0.429 0.300 0.129 -­ Tulsa EPI7 Clinker cooler #2 -­ 4.607 3.225 1.382 -­ Tulsa EP27 Coo/kook unloading -­ 0.375 0.263 0.113 -­

Tulsa EP28 Coo/kook 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 

PM IO-Z.

5 (filterable) 70% PM (coarse) 

PMz.

5 (condensible) 30% PM (fine) 

Lafarge Building Materials Fabric Filters Trinity Consultants 


BART MorlAlino A-2 063701.0044 

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

Printed /013012006 

Source Emissions 

Location Point ID Unit Description 

Tulsa EP 15 IKiln, screw conveyors, & bucket elevator 1964 

Emission Factors for NOx & S02 based on AP-42 (1/95) Table 11.6-8, PM 10 based on baghouse manufacturer exhaust loading rate. 

2 Construction dates and source parameters as listed in Permit 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 natural gas heating value 0: 1020 BTUlSCF 

Source Emissions 

Location Point ID Unit Description SOl SO. NG. 

Emissions· 

PMIO 

(lb/hr) 

PM (coarse)4 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 

4 Condensible/filterable PM speciation were based on the following proflles: 


National Park Service guidance for PC Dry Bollom FF. http://www2.naturc.nps.gov/air/Pcrmits/cct/cctCcmcntKiln.cfm 


Coarse 6.70% 

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

FineEC 0.00% 

CPMlOR 75.30% Inorganic condensible PM considered sulfates per NPS guidance 

CPMOR 10.30% Organic condensible PM considered secondary organic aerosols per NPS 

Lafarge Building Materials Kiln Trinity Consultants 


BART Modeling A-3 063701.0044 

• ••• "~ - ....... '. ," 


'''~,•.•' •••••. 

Printed 1013012006 

Source Emissions Emission Factors :,L 

Location Point III Unit Description Constr. Date Heat Duty Exhaust FI,owrate Hours of Operation NOx PMIO SOl 

i(MMBtu/hr (aefm) (hrs/day) i(lh/MMBtu)i (et/acC) (Ib/MMBtul ----------------­

Tulsa EP04 #2 discharge chute, separator, bucket elev. 1964 1.2 ' 36,000 24 0,075 0.025 0.001 

I NO. & SDz emission factors based on AP-42 (7198) Table 1.4-1 and 1.4-2,PMIO based on baghouse manufacturer exhaust loading rate. 

2 Ib/MMSCF eonvll1ed to Ib/MMBtu using a natural gas heating value of 1020 BtulSCF 

Source Emissions , Emission Data Inpue 

Location Point ID U nit Description SOl S04 NO. PMIO PM...,.. PMn.. EC SOA 

(Ib/hrl 

Tulsa EP04 #2 discharge chide. separator, bucket elev. 0.001 .­ 0.089 3.857 0.964 2.893 _. 

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

For Natural Gas, National Park Service guidance for NG·fired combustion turbines. http://www2.nature.nps.gov/air/Pennits/cct/cctGasFiredCT.c!in 

PMIO•2.5 (filtcrablc) 25% PM (coarse) 

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

Lafarge Building Materials Combustion Sources Trinity Consultants 


BART Modeling A-4 063701.0044 


j 
It•••':.•;.·•••'••••;. ~,r: . ':,:••·.,·.,1 

Printed 1013012006 

Stack parameters (Eng/ish) 1 

Location Point ID Stack Flow 

Tulsa EP04 36,000 

Tulsa EP05 4,000 

Tulsa EPI3 5,200 

Tulsa EPI5 160,000 

Tulsa EP24 31,750 

Tulsa EP26 4,000 

Tulsa EP25 4,000 

Tulsa EPI7 43,000 

Tulsa EP27 3,500 

Tulsa EP28 3,000 

Stack parameters (Metric) 

Location PointID 

Tulsa EP04 

Tulsa EP05 

Tulsa EPI3 

Tulsa EPI5 

Tulsa EP24 

Tulsa EP26 

Tulsa EP25 

Tulsa EPI7 

Tulsa EP27 

Tulsa EP28 

Exit Velocity 

ft/sec 

11.94 

8.29 

9.55 

17.83 

10.53 

33.16 

33.16 

3.56 

4.64 

4.91 

Exit Velocity 

(m/s) 

3.64 

2.53 

2.91 

5.44 

3.21 

10.11 

10.11 

1.09 

1.42 

1.50 

Stack Height 

feet 

100 

55 

97 

60 

100 

64 

55 

53 

7 

135 

Stack Height 

(m) 

30.49 

16.77 

29.57 

18.29 

30.49 

19.51 

16.77 

16.16 

2.13 

41.16 

Stack Dia. 

feet 

4 

1.6 

1.7 

6.9 

4 

0.8 

0.8 

8.0 

2.0 

1.8 

Stack Dia. 

(m) 

1.22 

0.49 

0.52 

2.10 

1.22 

0.24 

0.24 

2.44 

0.61 

0.55 

Stack Temp. 

Fahrenheit 

180 

70 

200 

500 

200 

70 

200 

200 

70 

70 

Stack Temp 

Kelvin 

355 

294 

366 

533 

366 

294 

366 

366 

294 

294 

Base Elev.2 

(feet) 

875 

820 

851 

768 

890 

854 

845 

788 

754 

8j}5 

Base Ek-v. 

(meters) 

267 

250 

259 

234 

271 

260 

258 

240 

230 

245 

I Source parameters obtained from Permit 98-1 74-TV 


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 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

_.

•

•••

••

•

•

•

••

•

•

•

•

•

•

•

•

• 

Station ID 

LIST OF SURFACE STATIONS 

Location 

Winfield/Ark. Ci""tv'-'..~KS",,-_______ 

Station No. 

(official ID Latitude (d .) 

7..7.4502 .___ .._.Jll(i§._ 

Longitude 

-97.033 

Bartlesvill""e.-"O~K-,--__________-,-7~23"-,5",6,,,-5__......_ ...3.!iL(i7...._ .... _._.._:2!i'Q.l1...._. 

Childress. TX 723604 ._....__.._..3.~Hn.._•._... -100.283 

CoffeYVille MunicioaL KS 724519.17JQQ....._.__.2~ 

.ChanutelM~n John. KS 13986.n,§§1. .. _.._ ......::2?,4.8.J ...._.. 

.... Clil).ton-SherIlllll!,_OK _____. 723526.3.5..Jn._._... __ .:22,~00 

..... _DentO!LTX _______ 722589 ...........:27.J.~.3...._ 

El Dorado/Goodwin_ AR.-'--_______---"9""3"'99"-=2'--_ 11l.!L... _........----=2.b..817 

.......................Yanfll AF!JlEnid. OK 723535 __ ......•. }!i13.3._......._._ -97.91_7_ 

...........She!:!!]lllllI?.ennison. TX _____.___---'-7-=22""6"-'7-"'3__........ _.}11J7._._...... _.~2§J:i§7._ 

......_..Fred~Jick._~_____.___. 22222 ..._3.413.3. __ .... __ ~2~,2~L 

Fliooin. AR 55555 ...3.!iJQQ...___... -92.467 

...... Fort~mith·_M !3964}_~13.3.. -94.367 

.FavettevilleLPrak~ AR __.__.. 93993 .. }_!i99Q. ......:24J!i?......._... 

Gae:e. OK 13975 .._ ..3.!iJQQ..._ ...___ :22]§.Z ..._ 

Hobart OK 93986. __...3.2,.9.QQ....._.___..:-2.2,.9..29.__ 

....,H"'o-"tS",o,..nn",,·",e:s",­....,AR=--__________7'-'2"'3"'41'-=5'___.....___.?_4,4§7.._._._~·9'_=3"'.0'"8",-3_ _1 

...,Harri=-,-·",·so",nIB~o",o",n",-e-",C",-o.,,-.AR'-='--_______--"13"'9'-'7.!.1__..___..3.§.,ln___...__-~93"'._"11.'_'7_ _I 

McConnell AFB. KS 3923 .___.3.7.61L__...__~2n§L.. 

Wichita KS""'--___________ 39~8 __._._._Jl&?Q. __...__ .. ~2H?3 

..._.l::!........_........_}oolin Rel!ionaL MO ._______--'7-=2~34"'9"'5___.._._._..3.7-15Q.._ ..__.__~94..299.._.._ 

Lawto"'n.-"O"-'K"'-__________~99"'9"'99~__.....__._...}4,.2!i7.... ___._._~-9'-'8"-.4'-'-1_'_7__1 

Liberal. KS 11111 ......_._ ..3.l,m~......_._.... _ -100.950 

Little Rock/Adams. AR________--'1-"'39~6""3'___ ......__..}4.21~......__.... _ ......:27,7..3.? .._. 

Little Rock AFB. AR 723405 ...._ .. .3.4.JE..._._....._-=27.J19 ..._.. 

..._._Q.. __.... Mu~koe:ee/Q!lvis Field. OK 11112. __ .3.~§..?Q. __._...____:2.2d§7 ___ 

McAlester Ree:ionaL OK"-______----"9-=-39"-'5"'0'___ _.._.~4.JQQ..__ ..._~-9'_=5"--'.7'_"8"'_3_ _I 

Mena International Mun.. AR 723435 ........._.._J!JJ.Q.._ .........._--'-'-9"'4"'.2C'-17'--. 

.......'OklaJ:\.oma City. OK ____-"'13<-"9"'6.!...7__.........._..~5..~4QQ. __....._.._.~2I§QQ.. 

_J-l·orman-vtesthe!m!l!:. O~__. ___7'-'2"'3"'57'_"0'--__..... _._ ..3..5...7.1~.._................._.:2H?Q.... .. 

Medicine Lode:e. K""S_______--'7""2""45""2"'0'--_...._._ .... E.7.J!3...._............._---'-9-=8=.5",50,,--_ 

.Poll~.£)tv·_91L.__._._..__....___-'1-=-39"-'6""9__ 

'par~on~ri <;:itv.JS.L_ 745431 ..J7}~? 

........ ParislCo?' Fi~ld. IX__________-=3"'33"'3""3'___ ......_........It§3.3... 

.... Ro~ers. AR _. 66666' ......._......J!iJ.Q!.. 

_.Sori!Jgfielq. MO _____ 13995 ........._.JV~.3..... __ ....... 

...........~95.517 

..........~95.450 

_:24J99..._ 

-93.383 

Lafarge Building Materials Trinity Consultants 


. Refined BART Modeling Protocol 063701.0044 

.... Siloam SoriJ:!l!. ~'--___ ___..!...77"-,7,-,-7..!...7__.......__._ ...l(Lm..._.___ ..._._:2~A~J. __.... 

Wichita Falls .. TX ________--'1"""3"'966""'-__....__.... IL9~~..._._......__-_"_9"'8.=_51"_'7'___I 

Stillwater. O""K"'-___ _____~1l~1~13'---__...... _~§..J..?:Q.._ .. _.._-97.083 

.....__.Fort Leonl!!d Wood~.. ___---'-'72"'44~57'--_...... _..._._J?7~Q_ ......___ ~~!?L 

........3ink!!r AFB".....-"'O"'K'--________ ._.. 72~?~0 ...__ ....._J..?41L_......._......._...._~.97.cJ.~L... 

____._____-'7-=2=.o25'-'4..... 1___....._._ ...J?cHJ_.__ .... _... __ ~9§..583 

Tulsa lnt. Airnort. OK 13968 ...... _._._..J§,;WO...___...._~-9'-=5"'.8""8"'-3_ _I 

Texarkana ReeiolJll"'I.""TX'-'-"-_______.....1"""3""97""7__.......__._.....lH..?JL...._ .... _~-9'-'4"-'.0'_"0"'_0_ _I 

Bentonville. AR _~_.....;8"-'8"'8"'-88"____......_ ......J§,}..?...Q...._..__.___.... -94.~E __ 

M Ardmore Municioal. OK 723555 34.300 -97.017 

LIST OF PRECIPITATION STATIONS 

Cooperative 

station ID Location/name State 

Latitude 

(deg.) 

Longitude 

(deg.) 

30130 ALUM FORK AR 34.480 -92.510 

~O 17~__.._ .._... ANTOINE ._____.___.. _...;..A.c:..R'--__.}±.2.~0_..__-_9_3._25_0_..._ 

30220 ARKADELPHIA 2 N AR 34.090 -93.030 

30458 BATESVILLE LIVESTOCK AR 35.500 -91.460 

30530 BEEBE _______........:..AR=-=-_ .......J...s...:Q:I.<L.._____...~2!_.5_40__.._ 

~.Q.~.~....__._. ____..!'!~~YVI:....:..:=LL""E=..::.5..:..:NW:..:..:....._ ...__.__..__..:..:A:.::R'-­..._~~.,~Q...._ ......_.:~3.3:?()_._ 

3.Q_764_.. _._..._,B=.:L=.:Ac.::K:.:;::EL=.Y:.....:.:M.:..:O:..:UN::.:..:.T.:..:Ac.:::I:.:..N:...:D::.:Ac..:.M:.:..:..-___._:...:A::.:R=--_.....;3=...4;.:,:.3:..:4c::O____-9_3._12_0_ 

~Q79~_ .._____B=-L=.UE=_"M.:..:O:..:UN::.:..:.T.:..:A...::IN:.:..:...:D::.:Ac.::M=-_____---=-A.::.:Rc:......._ _=3.::..5.:.:.0.:..;70'---__-...;..9::.:3..::..39:...:0:........._1 

~.083?._._.. 

30842 

..._._B_O_O~_EV_I_L_L~_3_S_SE__.______._A~R___ }2:.Q~.Q...._...._ ..._...~~:.s2Q_... _ 

BOTKINBURG 3 NE AR 35.430 -92.280 ............................................_-_.._---­ ................................................_....... 

30900 BRIGGSVILLE AR 34.560 -93.300 ..... _----_...._.._....._.........._....... - ..........­ .........­

31020 BULL SHOA."'L:..;S....;;D,~A...;..M _______---'A"'R:..:....._ ...:3....:c6':"c2:,;;;2,;:...0___-9-=.2.;.c.3...;..40"----I 

~3-ll-4-0-_.--~C..:..:A~L.:..:10::.:N~L~&:..:D~--------:...:A~R~--=33::.:..:..:18~0~_~-9~2~.2~90~_1 

3JJ.~....._ ..___.__ CI\MD§~_ ..__.. __._____.__........:...:A=-=R_..._2l.}~Q_ ... ___ ._~f.,:I.?Q_ 

31457 ..........._._......_ ... -_....__CLARKSVILLE 

........_ ...__._.__.........6 _-_. NE __...._._.._._._._:..=.:_ AR .....__35.320 ....- .........._._......_ -93.240 .... __......_....... 

~!582 COMPTON AR 36.050 -93.180 

~!952..___.....__..__p=..:E::...::sQ.::U=E=EN:..:...::D.:..:A:.:cM-"-_______--=..:A::.:R'--_....:3:..:4-'.:.0.:::,6°"--__-::..94-".2:::2:.::0__ 

32020 DIERKS DAM AR 34.0?Q....___._ -94.95o__.._ 

~.?~.~..__ I?Y._~~.........__.___________ ......A~....__._.....33 .530=~:f~<L..... 

~230Q_.. __......_ ..E_L_D_O_RA__D_OgOo.'2..."Y..IN_FD___ 

P2356 

32444 

EUREKA SPRINGS 3 WNW 

FAYETTEVILLE EXP STN 

AR ...~J~g__.__ ..-92.:I.~__ 

AR 36.250 -93.470 

AR 36.060 -94.100 

~j.~_.___._cF:...:E:.:RN.=..::D=.:A...::L=E=...6=-=E________....:AR-=.:..._ ..._.~_::I.~Q_____._:?_2.:_27_0__ 

~.s._"±...._.... __._ FOREl'1AN __________...__---=-A.::R-=--_ _~:~Q.._._.........::?4_.2_3Q......._ 

FORT SMITH MU, OKc:....._____..__.:..:A=R'---_.._35.:?QQ.._............::?~:229.........._. 

1 

Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • 

• • • • • • 

• • •• 

• •• 

• • • 

• 


•

•

•
•

•

• 


•

•
•

•

•

•

• 


•

•
•

• 


•

•

•

•

• 

~~.22i .. __ .__-,G::.::I=LB::.:E=Rc:.:T,---____._____---,-A""R_,_____3_5._59_0_ ___-92_A_3_°_ 

3~81 O..._____~G",I:::::L::::LHA:o:.:~M'--___ ._____.__.:.:A::.:R'--_..:3:c.:4~.1=20=--____-..::.9.-"4..:;.:15:.;0'-­

32978 ..._ .._g_RE_ERS .~E~Y I?~_M__.__..____-,-,A::.:R,--_ _ ~}l.L_. ____::2~-,QQQ_ ... 

~}1}~. HARDY AR 36.170 -91.280 •••••••••••_ ••••H ••••••••_ •••••••••••••••••••• _ ••••••••_ •••••••••••••••••••••••••••••• _ •••••• 

33165 HARRISON FAA AP 

~¥__...___.._.::.:HUN:.::.:..:..::T=-SVI-'-=L:.::L:::E--"I...::.Sc::S-'-W'--_._____-'A-=R"--__}_6_.0_4_0___•-93A50 

34185 LEWISVILLE AR 33.220 -93.340 

34248 NORTH LITTLE ROCKfMUNICIPAL AR 34.500 -92.160 

~.~?56.__..._....... _.MENA._._..._ ...___..__._._.. __ ...._._A:...:R~_ ..._~~:~Q..._ ....__..:.~±l~.9.._...... 

3~839 __ .__._._:.:..M:.::IL=L~..:Wc.:..O::::.O:;;;.::::.D...::D:.:.AM:::.:..__________=_A.::.:R~_ __:::.:33::.:.A'_'I..:::.0___----'-~_3._59_0_ 

3~988._._._"'M:-=-O::::.UN.::.:..:.T.:....:.:ID::.:A..:..::..3..:::S::::E_______._:..:A.:.:R'____ _=3:..:.4:.=..3=2.0:=--__-9::..:3:.:.;.3:.:5..:::.0_ 

35110 ..._ ..._N_A_~_OW_S_D_A_M__....___._.._._.:..A:.::R=---__ .}..±:Q~Q.. .._._::2}~30 .._ 

35112 NASHVILLE AR 33560 -93510 ......._-_._---_......... -... 

35200 _... .N_I_MROp. D~~___.____.. ___ AR.}..±??Q.. ..:.~~JQ9_.__ 

35228 NORFORK DAM AR 36.150 -92.150 

~2_0_.._.___.::..N:.::O:.::..R:.::T=..:H:..::L:;:.ITT.:...:.;:L::.::E:...;R;..:.O::...C;;;.:KfMUNI=:.;;;;.::..;.:.C.;;;.;I:.:.P:...;A:::L._....:A..::R:..:.-_--"--34.cc.:...;50;..;;0__---'--9..;;;2.;..;;.1.::..60'----I 

35602 PARTHENON AR 35570 -93.150 

~??~'! .......__...._ .....__Xr..NE .!3.!::.IB:.~___....._ ....__...___._..::.A""R_'___ .......l.~!..±9.....__ ... __:.?~1..9_.._. 

3?2.Q§. ...____._cP:...;RE..;;;::.::S-"C-"O..::.T..::.T__________---'-A:.:;R=---_....:3::.::3c:..A"'8::::.0___-..::.93::.:."'23::.::0'-----1 

3?393....•_.__.._.:S_T-"P.c..A-"UL-"-'..___________A'-'R'-'--__3~5_._49_0____-~346.9._ 

36920 

...__....._....__.....__._.. STUTTGART 9 ESE ..._._....__...._--'----_ AR ......_._34.280 .........._.._......__....-91.250 _........._......_­

37488 WALDRON AR 34.540 -94.060 .........................._...... 

141427 ........_ ......_.C_HA_.Nl[T_E_l'v!:!Q.Jil'I~QN ~.!'_...._ ..._..:.K:.::S'--_......}2:~.9..L ..._ ......::2.?.~2.9.___ 

._14_16_99__..._.__....:C:..:O:.:::;L=B..:;.Y....:1...:S:...;W-'--________---CK:::S=--_......;;:..39:...;.::..24.:..;:0___-..::.1.::..01;::..0:..;4:..:;.0_-1 

~.~4----._-=D~O-=D...:G:.::::E....:C:::IT.::..Y.::...::..MU~,.:.:K=S--------'K:.::::S=---~37:..:.A~6:.::..0__~-9.::..9:.=.5.::..80~_ 

142432 ......_.. __ ... EL_KHA_.RT __.______ _",",K=.S_ .... _}?:.929.._.:.!9.1.:.?}.9.._ 

142980 ..__QARJ?JO:1i..£!.E EXR STI.'I___._.__.......;:;K=S_ .... _~.?QQ9..._ ..... _~}9..9..:.~~L 

1~}95~____..._cIND,"-'-":.=E"-P::::ENDo..:.o::EN"'-'-C=E"-_______---"K"'S'---_..:3:.:.7.:.:.1:..:. 40~___-:.;95:::.A=20, 

143984 lOLA 1 W KS 37.550 -95.260 r-...._ ....--...- ..._­

145175_MEDICINE LODGE________K..-S"___.}?:!.?L.._.:2.?}30_ 

147922 SUBLETTE KS 37.290 -100.500 

148830 WICHITA KS. KS 37.390 -97.260 ..- ..._.._......_..._ ...._._..._..__.....__._...._ ...•----'--_..•....._........._.•....--.....•.-..........- ......•-­

~_:3_078_9_ .._____"'B::..:O::.=L::.:I..:..VA:...:.:.:R__________....:M=O__...::3.:..7::::.3..:.:70:.... ___....:-9...::3...::.2::.;30'-_--1 

~32302 DORA MO 36A70 -92.140 

~~315._____. ___ JOPLll-!.. ICIPAL..Af....___. MO 37.090 -94.300 

@}~~~.~.. _ ...._ .. _.!!RK.~Yl~!-E__..._..._ ... ___ .._._. MO 40.120 -92.340 

~.~~~___.__L:::E:::B.:.:A::..:Nc:::O.:..N:...:2=-W'-'-________....:M:.::.O::::....._ 37.410 •.__:?24~~ 

~~191 ~ __ .._._.. LICKING 4 N MO 37.330 -91.530 

@350~2......_..!::OCKWQOD _.______. MO}7.23.L._._.:2},5l9... __ 

@~?}9.?. MARSHFIELD MO 37.200 -92.550 

@}??~'!... _......MILLER _-_....__..__1 E ..._..._..._._.__._..._.......- ......_.:.:..:..::'--_ MO .......37.130 

__..._ ......_........._.._...•._........
 -93.500 

~36460 OZARK BEACH MO 36AOO -93.080 

Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044 

• • • • • • • • .._--.__..--.-­ ..............-~ 

SPRINGFlELDIMUN. MO. MO 37.140 .93.230 ..-.~--...­

STOCKTON DAM MO 37.420 -93.460 

TABLE ROCK DAM MO 36.360 -93.180 

MO 36.450 

OK 36.450 -96.000 

-~--. 

BENGAL ...~-......­ ..--~----------~-.....--. _c::0.:.::K____..._34_5!_o___ -95.050 

.102.290 

ClllCKASHA EXP STN OK 35.030 -97.550 

DUNCAN AIRPORT OK 34.290 -97.580 

-'---__________O:::;K"'-_ 35.230 

OK 36.450 -98.080 

OK 36.090 -97.370 

OK 34.530 -95.470 

OK 35.460 -95.200 

, ....._ ••••_ •••m .._._••••_____ NOWATA • ___••_~~__""~__..__._. OK 36.420 -95.380

-"'O:::IL:o.;T:.oO;:.:N.:c2=.S""E:::.......________---'O::.:K-'---_ 36040 -96.340 

OKEMAH OK 35.260 -96.180 

OKLAHOMA CITY, OK OK 35.230 -97.360 

______.!'_A_WNE_._E_5_N____ ....______9JS.._..._---.?.~~40 ............... ~?~~?(L.. 

'--___---':.==:...=.::.:..:...::.:..o::.:..:.:..:..:.:'----_____...::::O.:.::K'-----','--1§.:.±:IQ....~g60 .._ 

OK 36.330 -101.050 

OK 36.470 -100.250 

OK 36.170 -99.530 ._--_.......- ...... ~~,.. ,-,~... ,-,-''''''-,,~,~. 

_____.----'O""K"'--_. 36.070 

OK 

OK 

36.320 

34.440 

32.250 

-95.530 

-98.530 

-98.430 

-99.410 

_ ....____....___A_MA__RI_L_L_O_m:!!!::J.:..;:..X_____...::T.:.:X=--__}_5_.I~L...___._...~19._~.:4_20__ 

ClllLDRESS MUNI AP TX 34.260 -100.170 

•

•

•

•

•

•

•

•

•

•

•

•

•

•

•

• 


-. 

•

•

• 


-e

•

• 

•

• 


• 


". 

•

• 


_. 

•

•

• 
 Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044

_.••

• 
• • • • 

• 


•

c. 
 r-...._ .....__......-.-.~.---------------------

~! 1773 CLARKSVILLE I W 

~U~L____.__-=.C.=O=M:.::M..::E=R::.:C:.=E,--_________:.:..:..__=.:..=__--'...:c==-_I 

~.'-!.~.?L ...__<;:ON~OE__._....__.____.___ __ 

~l?~.~L DALLAS-FORT WORTWFORT. TX. ........- ............_.._........._..... .................._..............__..... 

~-'-~.~ ______._ DA_l:LA~.LOYE ~!~LQ___.______--"",,-_ 

~12404 DENTON 2 SE 

~13284 FORT WORTH MEACHAM F 

~.~285____.... FORT WORTH WSFO 

~}370 ___._£.~::.::S:..::Cc:::O____.______.____-==--_ 

~.. GRAPEVINE DAM 

~1409~_.._. HEREFORD 

~143??..._ .. ___..!:fUN_.!~~_______...__._.::..:..:..__ 

~l~?IL__ JACKSBORO TX 33.140 -98.090 

~IJ..'.~~.___._'=§.?9_N_GTO~....______.._..__._...__..:..TX""-_..... ~,??L __...-97:.Q!Q_ 

~15411 LUBBOCKfLUBBOCK INTL TX. TX 33.390 -101.490 

~15897 MlDLOTIDAN 2 TX 32.290 -97.000 

~1958_._._ MINERAL WELLS MUN AP ___..:..TX",,-__ ..1~j70 .._ ..__..J80'!,9..._._. 

~_?.108__....._._ MOUNT PLEA~ANT__.__._____---'T;.:.X'--_ _2JQQ.. __--=2.5_0_09_. 

~!6177 NACOGDOCHES TX 31.370 -94.390 

~!_6776 _ ...;:.P_A_M...;:P_A...;:2'--________._..;:..TX"'-_--'.3c:...5..;:..34..:..:0'--_--=-I..:c00"".c:...58",,0_1 

~170~ PLAINS ________....:.T.:,:X'--__..:3:.::.;.:31_10___:102500__ 

4.!_~.~L STEPHENVILLE I N TX 32.150 -98.120 

~!!?:!.~ .. __..!Lu.r:-_PHl!.~Jf~(J~._.___.._ .....___.. _....:.T.:,:X'--_.. _}l:.Q~L___~?}~_ 

418942 TEXARKANA TX 33.250 -94.050 

419532 WEATHERFORD TX 32.450 -97.460 

419729 WlCIDTA FALLS/SHEPS AFB TX TX 33590 -98.300 

Lafarge Building Materials Trinity Consultants 

Refined BART Modeling Protocol 063701.0044 

•• ....ot ,.. ~ 
l 

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

••;•.•.".~

'''>'. 

37 


36 


35 


34 


33 


.; .';<" ., Dodge Cily-E Prall-W Prall-E Wichita-W Wichita-E Joplin-W Joplin-E Springfield-W Springficld-E RolI.-W 

'. 

Pcrry1on-E Woodward-W Woodward-E. Enid-W Enid-E Tulsn-W Tulsa~E Harrison-W llarrison-E Poplar Bluff-W 

Plainvlcw-E Chnlon-W Clinton-E Oklahoma Cily-W Oklahoma Cily-E FortSmilh-W FortSmilh-E Russcllvillc-W Russcllvillc-E Mcmphis-W 

AmariUo-E Lawton-W Lawton-E Ardmorc-W Ardmorc-E Me Alcslcr-W Me Alcslcr-E LIllie Rock-W Lillie Rock-E Helen.-W 

Lubbock-E Wichila Fall,·W Wichita Falls-E Shcrman-W Shennan-E Texarkana-W Texarkana-E EIDorado-W EI Dorado-E Greenwood-W 

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


! 


~ 

~ ~ 

,

~. 
. . ~ " " .,~\~~:~~~ 

' .......... . .:'. ,." ••'., ·'.J'.~I -.. ;:' >':':'(~~~~~

List of US Uind Cover (LULC) File,s (1 :250,000 km 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 

Amarillo 

Clinton 

Lawlon 

Oklahoma City 

Ardmore 

Fort Smith 

Me Alester, Antlers, 

DeQueen, Mena 

Russellville, Conway, 

Flygap Mountain, 

Mountain View 

Little Rock 

Memphis 

Helena 

Lubbock Wichita Falls Sherman Texarkana EI Dorado Greenwood 

102 100 98, 96 94 92 
ApPENDIX C - MODEL INpUT/OUTPUT FILES 

Lafarge Building Materials Trinity Consultants 


Refined BART Modeling Protocol 063701.0044