# •

• 

# • • 

# • • • 

# • • • 

# • • • • • • • 

# • 

# • • • • • 

# • • 

# • 

# • • 

# • • 

# • • • • 

# • • 

# • • 

# • • 

# • • • • CALPUFF 

## BART MODELING PROTOCOL • REFINED ANALYSIS LAFARGE BUILDING MATERIALS INC.

TULSA CEMENT PLANT 

### Prepared by: Vern Choquette • Principal Consultant Don Whitney, PE • Managing Consultant Eugene Chen, PE • Senior Consultant Jeremy Townley • Consultant TRINITY CONSULTANTS 

### 120 East Sheridan Suite 205 Oklahoma City, OK 73104 

# REeEl \J

# E0 NOV 

# n 2 701'1-\ .I 

# A\R QlJJ\LlTY -­ (405) 228-3292
November 2, 2006 

### Project 063701.0044



# n.

# TritlrityA '-LJnsu tants

# • • • 

# • TABLE OF CONTENTS 

# • • 

# • • 

# • • • • • 

# • • 

# • • 

# • • • • 

# • • 

# • • 

# • • 

# • • • • • 

# • • • 

# • 

# • • • 

# • 

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

# • 

# • 1. EXECUTIVE SUMMARY•••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• 

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

### 1-1 1.2 OBJECTIVE
............................................................................................................ 1-1 1.3 LOCATION OF SOURCE AND RELEVANT CLASS 

### I AREAS
....................................... 

### 1-2 

## 2. CALPUFF MODELING SYSTEM.•••••••••..••.•....••.•...•.•••..••.•...•••.•.•.••..••.••••.•••••.•..••. 

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

### 2-1 2.2 MODELING DOMAIN
............................................................................................. 

### 2-2 

## 3. CALMET..••••••••••••••••••••••••••••••.••.•••...••••••.••••.•••••••••••••.••••••••••••••••..•••••••••••••.•••••••• 

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

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

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

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

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

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

### 3-1 3.2.3 UPPER AIR DATA
................................................................................................ 

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

### 3-2 3.3 CALMET PROCESSING CONTROL
........................................................................ 

### 3-3 

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

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

### 4-1 4.2.1 PM IO EMISSIONS
................................................................................................. 

### 4-2 4.2.2 STACK PARAMETERS
.......................................................................................... 

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

### 4-3 4.4 BACKGROUND OZONE
.......................................................................................... 

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

### 4-6 4.6 CALPUFF MODEL CONTROL PARAMETERS
........................................................ 

### 4-7 

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

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

### 5-1 5.3 NATURAL BACKGROUND
...................................................................................... 

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

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

### 5-2 APPENDIX A - EMISSION DATA 

### APPENDIX B - DATA 

### FILE LIST

# •



# • • • • •

LIST OF TABLES
TABLE 2-1. CALPUFF MODELING SYSTEM VERSION NUMBERS
........................................................ 2-1 TABLE 2-2. VERTICAL GRID STRUCTURE
.......................................................................................... 2-3 

# • • • • 

# • • • 

# • TABLE 3-1. SUMMARY oFCALMET INPUTS
.................................................................................... 3-4 TABLE 4-1. CONSTRUCTION DATE SUMMARy
................................................................................... 4-1 

# • 

# • TABLE 4-2. ESTIMATED 24-HOUR AVERAGE EMISSION RATES
......................................................... 4-2 TABLE 4-3. SPECIATED PMto EMISSIONS
.......................................................................................... 4-3 TABLE 4-4. POINT SOURCE STACK PARAMETERS
.............................................................................. 4-3 

# • • • 

# • 

# • • TABLE 4-4. SUMMARY OF CALPUFF INPUTS
................................................................................... 4-7 TABLE 5-1. MONTHLYHUMIDITY FACTORS " ..•... ;.;.; ........ ;.; .... ;..... ;.., .... :, .......**~;**.. :.; .... ;........ ;;;;; .... :,:: .... 5-1 

# • • 

# • TABLE 5-2. DEFAULT AVERAGE ANNUAL NATURAL BACKGROUND LEVELS
.................................. 5-2 TABLE 5-3. SUMMARY oFCALPOST INPUTS
................................................................................... 5-2 

# • 

# • •
• 

# • • • 

# • • • 

# • • • • 

# • Lafarge Building Materials 11 

### Trinity Consultants Refmed BART Modeling Protocol 063701.0044 

# • 

# • • •

# •

• • 

# • • LIST OF FIGURES 

# • FIGURE 

### 1-1. LOCATION OF TuLSA PLANT AND NEARBY CLASS 

### I AREAS
......................................... 

### 1-3 

# • 

# • 

# •

FIGURE 

### 2-1. PROPOSED MODELING DOMAIN
.................................................................................... 

### 2-2 

# • 

# • 

# • FIGURE 

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

### 3-1 FIGURE 

### 3-2. PLOT OF SURFACE STATIONS
......................................................................................... 

### 3-2 

# • • • FIGURE 

### 3-3. PLOT OF PRECIPITATION STATIONS
............................................................................... 

### 3-3 

# • FIGURE 

### 4-1. CANEY CREEK WILDERNESS RECEPTOR LOCAnONS
.................................................... 

### 4-4 FIGURE 

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

### 4-4 

# • 

# • • • FIGURE 

### 4-3. UPPER BUFFALO WILDERNESS RECEPTOR LOCAnONS
................................................ 

### 4-5 FIGURE 

### 4-4. WICHITA MOUNTAINS WILDERNESS RECEPTOR LOCAnONS
........................................ 

### 4-6 

# • • 

# • • • • 

# • • • • • • • • • 

# • • • • • 

# • • 

# • Lafarge Building Materials ill 

### Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# •

# •

• 

# • • • • •

• • 

# • • • • 

# • 

# • • 1. EXECUTIVE SUMMARY 

# • 

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

# • 

# • 

# • Rule (RHR). The objective ofthe 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 BARTrule defmes BART-eligible sources as sources that meet the following criteria: 

# • • • • • • (1) Have potential emissions ofat 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 ofthe 26 listed source categories in the guidance. 

# • 

# • • A BART-eligible source is not automatically subject to BART. Rather, BART-eligible sources are subject-to-BART ifthe 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 ifthe visibility impacts from a source are greater than 0.5 deciviews (dv) when compared against a natural background. 

# • • 

# • Air quaJity modeling is the tool that is used to determine a source's visibility impacts. States have the authority to exempt certain BART-eligible sources from installing BART controls ifthe results 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 ofBART modeling protocols and to 

# • 

# • harmonize the approaches between adjacent RPOs, the Central States Regional Air Planning (CENRAP) organization developed 

### BARTModeling Guidelines (December 15, 2005). The intent of the guidelines is to assist CENRAP states and source operators in the development ofstatewide and source-specific modeling protocols. 

## ·1.2 OBJECTIVE 

# • 

# • The objective ofthis document is to provide a report summarizing the modeling methods and procedures that will be followed byTrinity when conducting refined modeling for Lafarge Building Materials' Inc. (Lafarge) Portland Cement Plant in Tulsa, OK.
1-1 Lafarge Building Materials Trinity Consultants, Refmed BART Modeling Protocol 063701.0044

# •

• • 1.3 LOCATION OF SOURCE AND RELEVANT CLASS 

## I AREAS 

# • • • • Lafarge's Tulsa Plant is located due east ofTulsa, OK and operates as a portland cement manufacturing plant. It 

### operates two dry-process kilns that primarily consume a combination ofcoal and coke. Lafarge has identified the sources listed in Appendix A at the Tulsa Plant as sources that meet the three criteria for being BART-eligible sources. These sources include:

.. Unit 221, Raw Material Grinding #2

.. Unit 222, Auxiliary Raw Material Grinding #2

.. Unit 223, Auxiliary Raw Material Grinding #2

.. Unit 404, Kiln #2 Feed System

.. Unit 502, 

### Kiln #2

.. Unit 706, Finish Grinding #2

.. Unit 705, Finish Grinding #2 System

.. Unit 704, Auxiliary Finish Grinding #2

.. Unit 602, Clinker Cooler#2

.. Unit 1101, Coal Unloading from Railcar

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

.. Unit 1103, Coal Silo (crusher) A detailed summary ofBART-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 ofthe plant will be included in this refilled analysis. The following Class I areas fell within this 300 

### km radius: 

# • • • • • • • • • • • • • • • • • 

# • Class I Area Name 

### Distance from source (km) 

# • 

# • • Caney Creek Wilderness 247 Hercules-Glades Wilderness 258 Upper Buffalo Wilderness 213 Wichita Mountains National Wildlife Refuge 300 

# • 

# • 

# • • 

# • • • • • • 

# • • • • • • 1-2 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# ••••••••••••••••••••••••••••••••••••••••••••

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



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


### 1-3 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •

• • • 

# • • • • • • • 2. CALPUFF MODELING SYSTEM 

# • • • • • • • • • • • • • The main components 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 ofair pollutants emitted from a source in "puffs". CALPUFF calculates hourly concentrations ofvisibility 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 ofthe CALPUFF modeling system and all related programs. The versions ofthe CALPUFF modeling system programs that Trinity will use for Lafarge's modeling are listed in Table 2-1. Table 2-1 also compares the program versions that Trinity used with the program versions recommended by CENRAP. Note that some ofthe program versions are not the same as the program versions recommended by CENRAP. Certain programs will be utilized in order to avoid data compatibility issues with previous versions ofthe CALPUFF system. Trinity proposes the use ofthe following software version numbers listed below. Model selection was primarily driven by a desire to utilize the most recent release ofCALMET. CALPlJFF and CALPOST executable versions were then selected as appropriate for compatibility. 

# • 

# • TABLE 

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



# • • • 

# • Program CENRAP Suggested Lafarge Ref"med Analysis Reason for Difference Version Level Version Level TERREL 3.311 030709 3.311 030709 CTGCOMP 2.42 030709 2.22 030528 Version recommended is not available 

### CTGPROC 2.42 030709 2.42 030709 MAKEGEO 2.22 030709 2.22 030709 CALMM5 2.4 050413 2.4 050413 CALMET 5.53a 040716 6.211 060414 Most recent release of CALMET 

### CALPUFF 5.711a 040716 6.112 060412 Use version compatible with CALMET v6.211 

### POSTUTIL 1.4 040818 1.4 040818 CALPOST 5.51 030709 6.131 060410

- Use version compatible with CALPUFF 6.112 

# • • • 

# • • • • • • 2-1 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •



# • • • 

# • • • 

# • • • • • 

# • 2.2 MODELING DOMAIN Trinity will utilize a modeling domain that extends at least 50 km in all directions beyond the Tulsa Plant and each ofthe four surrounding Class I areas. The map projection for the modeling domain is Lambert Conformal Conic (LCC) and the coordinate system will be North American Datum 1927 (NAD27), which is an LCC projection. The meteorological grid spacing is 3 km. The southwest comer ofthe modeling domain will be Latitude 34.00'N, Longitude 

### 100.5°W which is assigned as the 0, 0 reference point for the domain. The northeast comer ofthe 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 'c u 

### "**~** 250 § 

### " ;; 200 Cll 

# • • • • 

# • • • • • • 

# • • 

# • • • • • • • • • 

# • • • • • • • • 

# • • .c E .::J 150

100
50


### o FIGURE 

### 2-1. PROPOSED MODELING DOMAIN

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

Hercules-GIlldes 

# •**Okiahom~City_OK**

### NDrm~tOK WIChita_Mountaln.l' ED Lawton_OK**~-----------------------**
-50 0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 Lambert Conformal Conic (X) (km 1


### Forthe purposes ofthis analysis, the CALMET and CALPUFF doma41s will be and are singularly referred to herein as the "domain."

2-2 Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044

# •

• • 

# • • • • • • • • 

# • •

The horizontal domain is comprised ofgrid 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 ofeach cell indicates its vertical extent. The vertical domain is also composed ofterrain-following grid cells, the number and size ofwhich are chosen so as to constrain the boundary layer in which dispersion and chemical transformations take place. TABLE 

### 2-2. VERTICAL GRID STRUCTURE 

### Cell Face Height Vertical Grid Cell (meters) 1 

### 20 2 

### 40 3 

### 80 4 

### 160 

### 5 320 6 

### 640 7 

### 1,200 8 

### 2,000 9 

### 3,000 10 

### 4,000 

# • • • • • • • 

# • • • • • • • • • • 

# • • 

# • • 

# • • • • • • 2-3 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 

### 063701.0044 

# • 

# • • •

## 3. CALMET



# • • • • • • • 

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

# • • 

# • • observation meteorological data to generate three years ofCALMET outputs to be used as inputs to CALPUFF. The CALMET model requires the input ofgeophysical data, meteorological data, and model parameter settings. }he CALMET modeling procedures that Trinity will use generally follows the recommendations in CENRAP's protocol. It 

### should be noted that Trinity will use surface and precipitation data but not upper air observations in this CALMET analysis. As a result, some of CENRAP's recommendations may not apply.


## 3.1 GEOPHYSICAL DATA CALMET requires geophysical data to characterize the terrain and land use parameters that 

# • • • • • 

# • potentially affect dispersion. Terrain features affect flows and create turbulence in the atmosphere and are potentially subjected to higher concentrations ofelevated 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 willuse terrain data obtained from the United States Geological Survey (USGS) in I-degree (I :250,000 scale or approximately 90 meter resolution) digital format. A list of 

# • 

# • the USGS terrain files 

### to be used in the analyses is provided in Appendix B. A plot 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 Refmed BART Modeling Protocol 063701.0044 

# • • 

# •

# ••••••••••••••••••••••••••••••••••••••••••••

FIGURE 3-1. PLOT OF MODELING DOMAIN 

### USGS TERRAIN DATA



## tl**~ ~**  35 II "d 

# 'i .....l 34.5

34

315

Longitude (degrees W) 

### BOO 700 600 500 400 300 200 100 

### o Meters (Above MSL)



### 3-1 Lafarge Building Materials Trinity Consultants Refined 

### BART Modeling Protocol 

### 06370l.0044

# •

• 

# • •

• • 

# • • • • • • • • • • 

# • Trinity will use the TERREL program to generate grid-cell elevation averages from the USGS terrain data across the modeling domain. 

# • 

# • 3.1.2 LAND USE DATA 

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

# • • (LULC) data files at 1:250,000 resolutions. Where 1:250,000 land use data is not available, Trinity will use USGS data at 1: 100,000 resolutions. A list ofthe USGS land use files to be used is provided in Appendix 

### B. 

# • • • 

# • Trinity will use the CTGPROC program to generate land use for each grid cell from the LULC data across the modeling domain. Since the USGS CTG format LULC data files must be compressed prior to use in the CTGPROC utility processor, Trinity will compress the files using the program CTGCOMP. 

## 3.2 METEOROLOGICAL DATA 

# • CALMET will be used to assimilate data for 2001, 2002, and 2003 using mesoscale model output and National Weather Service (NWS) surface station observations, upper air station observations, and precipitation station observations to develop the meteorological field. 

### 3.2.1 MESOSCALE MODEL DATA 

# • • • 

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

# • precipitation observation data. The mesoscale data will be used to defme the initial guess field for the CALMET simulations. Trinity will use the following 5th 

### generation mesoscale 

# • model meteorological data (orMM5 data): 

# • • 2001 MM5 data (12 km resolution), • 2002 MM5 data (36 

### km resolution), and • 2003 MM5 data (36 

### km resolution) 

# • • 3.2.2 SURFACE DATA 

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


# • • provides a substantial enhancement over the steady-state treatment ofobservations from a single meteorological station. Parameters affecting turbulent dispersion that are observed hourly at surface stations include wind speed and direction, temperature, cloud cover and. ceiling, relative humidity, and precipitation type. Surface data will be extracted from among the surface stations listed in Appendix B and are illustrated in Figure 3-2. These sources were selected from the available data inventory to optimize spatial coverage and representation ofthe domain. Raw observations from these stations were obtained from the National Climatic Data Center (NCDC), quality assured and were merged using the 

# • • 

# • SMERGE pre-processor to create a single assimilated data file ofsurface observations for each year analyzed. 

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

# • • 

# •
FIGURE 3-2. PLOT OF SURFACE STATIONS

i<SLG



### • KFW 

# . .. KCNU 

# • • 

# • •
• 500 I<I'PF KSGF r- - - - -KP2ii" -110

- - -KJUi - &- 1**~** I<t:FV 

### • i<l..BL**~** I<BVO I<PNC 

### 400 

# • • Hercules-G lades 

## E KFI.P**~G** KHRO 

# • 

# • • • 

# • • 2:S i<SWO**~** u 

### '(5 300 

### u ro KOKIfI1K I<CSM 

## .. KOUN 0 

### •**~** 

# .. KH8R 0 

### 200 

### • u I -e I<M1I1IT KHOT (]) Wichita_Mountainal I<LAW 

## .. 

## .. .0 

# .. 

# • • 

# • E ('iJ KFDR KADM 

# • 

# • • • c 

## E



## .. .. -l 

### 100 I<5PS 

# .. KF3I I<I'RX. 

## .. .. KTXK 

## .. IELD


## .. 

# L
___________ _______________ l**~_~**

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

## a 100 

### 200 300 400 

### 500 600 700 

### 800 

# • • • o

Lambert Conformal Conic 

### (Xl (km) 

# • • • 3.2.3 UPPER AIR DATA 

### Upper air observations

### will not be used in this CALMET analysis (i.e., NOOBS 1). 

# • • • 

# • 3.2.4 PRECIPITATION DATA 

### The effects ofchemical transformations and deposition processes on ambient pollutant 

# • 

# • concentrations will be considered in this analysis. Therefore, it will be necessary to include observations ofprecipitation in the CALMET analysis. The precipitation stations from which Trinity will extract data are listed in Appendix 

### B. The locations ofthe precipitation stations with respect to the modeling domain are shown in Figure 3-3. These 

# • • • stations were selected from the available data inventory to optimize spatial coverage and representation ofthe domain. Data from the stations were processed for use in CALMET 

# • • using Version 5.31, Level 030528 ofEPA's PMERGE program. 

# • 

# • • • 

# • 3-2 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • •

# •

• 

# • • • • 

# • • 

# • • • 

# • • FIGURE 

### 3-3. PLOT OF PRECIPITATION STATIONS

latA 1

### STO~\lA**LE~O** 500 

# • • • 

# • 

# • • 400

........ E 

# £ 

# • • • • •**.~** s:: 300' 0 (.) Ci E .E '­ 

# • 

# • s:: 0 200 

# • (.) 1:: II) .c E 

# • 

# • m ....J 100



## o -e!ikUT- LOCt:;..,.., 

### t - MNfSHF r,1'il.m'i ING**IND~ JO~UN** 

• GRhT WlYNOK**G~Y**

### ~O t:J:CI**CHI~KA** I VVichita Mountainsl I 

# • DUNCAN**~T**

Ct..4KS •

## ~O FOREMII. MILLIJIO 

# • • TEXARK
. LE\M S\I I CHltoR 

### I I 

### L ARDttOFl


### cot4JdiHbtoONT­ J,ocf<ss OelrON**FRI~CO** 

## o 100 200 300 400 500 600 700 800 Lambert Conformal Conic (X) (km)


## 3.3 CALMET PROCESSING CONTROL 

# • • • • • • Table 3-1 provides a listing ofthe CALMET parameters that Trinity will use in Lafarge's modeling analysis. In 

### addition to the parameters that will be used for Lafarge's modeling, the table also lists CENRAP's recommended parameters for comparison. In cases where a parameter to be used is different than what CENRAP recommended, a short explanation as to the difference is proved. Note that most ofthe differences from CENRAP's recommended parameters are due to the inclusion of observation data into the modeling analysis, since CENRAP's parameters are based on a noobservation analysis. 

# • • 

# • • • • 

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

# • •

# •

• • 

# • • • • 

# • • • 

# • 

# • • • 

# • • 

# • TABLE 3-1. SUMMARY 

### OF CALMET INPUTS 

# • 

# • 

# • • • • • 

# • 

# • • • 

# • • • 

# • 

# • 

# • • • •



# • • • • • • 

# • CALMET Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes NUSTA Number ofupper air data sites 0 0 Only surface and precipitation observations will be used NOWSTA Number ofoverwater data sites 0 0 IBYR Starting year 2001**l~ate** Years 2001, 2002, 2003 ar IBMO Starting month 1 Appropriate month Due to file size, analysis will be completed by month IBHR Starting hour 1 

### I IBTZ Base time zone 6 6 IRLG Length ofrun 6 Varies with month Due to file size, analysis will be completed by month IRTYPE Run type (I for CALPUFF) 

### I 1 LCALGRD Compute CALFRID data fields (T run CALGRID) F F

ITEST Stop run after SETUP to do input QA (2 run) 2 2

PMAP Map proiection LCC LCC RLATO Latitude (decimal degrees) ofprojection origin 40N 33.25N Appropriate for domain RLONO Longitude (decimal degrees) ofprojection origin 97W 100.5W Appropriate for domain XLATl Latitude of 1st standard parallel 33N 33N Appropriate for domain XLATI Latitude of2nd 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-4 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •



# • • • • 

# •

• • 

# • • 

# • 

# • • • 

# • • 

# • • • • 

# .. • • 

# • • 

# • • • 

# •

• • 

# • • • • • 

# •

• 

# • CALMET Variable Desc Value Included:ln CENRAP Protocol Value Trinity Will Use Notes NY Number ofY grid cells in meteorological grid 192 155 Appropriate for domain DGRIDKM Grid spacing (lan) 6.0 3.0 Refined grid spacing XORIGKM Ref. coordinate of SW comer ofgrid cell -1008 

# ° Appropriate for domain YORIGKM Ref. coordinate of SW comer ofgrid cell 0.0 

# ° Appropriate for domain NZ Number ofvertical layers 10 10 ZFACE Vertical cell face heights (NZ 

### + 1 values) 0,20,40,80, 160, 320,640,1200,2000, 3000,4000 0, 20,40, 80, 160,320,640,

1200,2000, 3000,4000 LSAVE Save met. data fields in an unformatted file? T T

IFORMO Type ofunformatted output file (I for CALPUFF) 1 1

LPRINT Print met. fields F F IPRINF Print intervals 1 I IUVOUT(NZ) Specify layers ofu,v wind components to print NZ*O NZ*O
IWOUT(NZ) Specify layers ofw wind component to print NZ*O NZ*O
ITOUT(NZ) Specify layers of3D temperature field to print NZ*O NZ*O

LDB Print met data and variables F F NNl First time step for debug data to be printed 1 1

NN2 Last time step for debug data to be printed 1 2 Will generate debug data for a total of2 time steps IOUTD Control variable for 0 

# ° 

# • • 

# • • 3-5 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# •

# •



# • • • • 

# • • 

# • • 

# • 

# • • 

# • • 

# • • 

# • • 

# • • • • • 

# • • • • 

# • • • 

# • 

# • • • 

# • 

# • • • 

# • Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044 

# • 

# • •



# • • CALMET Variable Description writing test/debug wind fields NZPRN2 Number oflevels starting at surface to print IPRO Print interpol wind compo IPRI Print terrain adjusted surface wind components IPR2 Print initial divergence fields IPR3 Print fmal w'**~ ~**  and direction IPR4 Print fmal di9 fields IPR5 Print winds after 

## I kinematic effects IPR6 Print winds after Froude number adjustment IPR7 Print winds after slope flows are added IPR8 Print fmal wind field components NOOBS No observation mode (2 = No surface, overwater, orupper air observations; use MM5 for surface, overwater, and upper air data) NSSTA Number of meteorological stations in SURF.DAT file NPSTA Number of precipitation stations in PRECIP.DAT file ICLOUD Gridded cloud fields (0 = no, 3 Gridded cloud cover from prognostic relative Value Included iIi CENRAP Protocol

0

0 0

0 0 0 0 0

0 0 2

0

0

3



### 3-6 Value Trinity Will Use

I

0 0

0 0 0 0 0

0 0 I

47

l35

3 Notes

Default

Will use surface and precipitation observations

Number ofstations

Number ofstations •

# •



# • • 

# • 

# • • • • • • 

# • 

# • 

# • • 

# • • 

# • • • 

# • • • 

# • • 

# • • 

# • • 

# • • 

# • • • 

# • 

# • • • CAL 

# •

scription Value Included in CENRAP Protocol Value Trinity Will Use Notes humidity) IFORMS Format ofsurface data (2 :; formatted) 2 2 IFORMP Format of precipitation data (2

formatted) 2 2

IFORMC Format ofcloud data (2 formatted) 2 

### I Nt A - No cloud data used in model IWFCOD Generate winds by diagnostic wind module? (1 yes) 1 1
IFRAD]
.Adjust winds using Froude number effects? (1 

### == yes) I 1

IKINE Adjust winds using kinematic effects? (0

= no) 0 1 Will compute kinematic effects in this analysis IOBR Use O'Brien procedure for vertical winds? (0 = no) 0 0

rSLOPE Compute slope flows? (1 = yes) 1 I IEXTRP Extrapolate surface winds to upper layers (-1 = no extrapolation and ignore layer 1 of upper air station data) -1 -1

rCALM Extrapolate surface winds even 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 RMIN2 Minimum vertical extrapolation distance Distance (lan) around an upper air site where vertical extrapolation is excluded (set to -1 if IEXTRP=±4) -1 -1

IPROG 

# • Using prognostic or 14 14 

### 3-7 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • 

# • • • e:

# • 

# • 

# • • 

# • • • 

# • • • • • 

# • • 

# • • • 

# • • • • • • • 

# • 

# • • • 

# • • 

# • • • • 

# • • 

# • • • • • • • • CALMET Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes MM-FDDA data? (14 

## =Use winds from MM5.DAT as initial guess wfud field) ISTEPPG Timestep (hours) of the MM5 data I 1 LVARY Use varying radius of influence to develop surface winds? T T Use FLAG default

RMAXl Maximum radius of influence over land in surface layer (km) 30 100 As included in IWAQM Phase 2 case study RMAX2 Maximum radius of influence over land aloft (km) 30 500 As included in IWAQM Phase 2 case study RMAX3 Maximum radius of influence over water (km) 50 36 Saine as MM5 data spacing RMIN Minimum radius of influence used anywhere (km) 0.1 0.1

TERRAD Radius ofinfluence ofterrain features (km) 12 12

RI Weighting offirst guess surface field (km) 1 50 As included in IWAQM Phase 2 case study R2 Weighting offirst guess aloft field (km) 1 100 As included in IW AQM phase 2 case study RPROG MM5 windfield weighting parameter (k:tD.) 0 0

DlVLIM Maximum acceptable divergence 5.E-6 5.E-6 NITER Max number of passes in divergence minimization 50 50

NSMTH Number 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-8 Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044

# •

• • 

# • • 

# • • 

# • • • • • 

# • • • • 

# • 

# • • • • 

# • • 

# • • • 

# • • • • 

# • • • 

# • • • 

# • • • • 

# • • • • CALMET Variable Description Value Included 

### iIi CENRAP Protocol Value Trinity Will Use Notes NlTR2 Max number of stations used in each layer for the interpolation ofdata to a grid point (NZ values) 

### 5,5,5,5,5,5,5,5,5, 5 5,5,5, 5,5,5, 5,5,5,5

CRlTFM Critical Froude number 1.0 1.0 ALPHA Kinematic effects parameter 0.1 0.1 FEXTR2 Scaling factor for extrapolating surface winds aloft NZ*O.O NZ*O.O

NBAR Number ofterrain barriers 0 0 IDIOTPI Compute temperature from observations (0

= true) 0 0

ISURFT Surface station to use for surface temperature (between 1 and NSSTA) 4 4

IDIOPT2 Domain-averaged wind component switch 0 0

IUPT Station for lapse rates (between 1 and NUSTA) 2 2

ZUPT Depth through which lapse rate is calculated 200 200

IDIOPT3 Domain averaged wind component switch 0 0

IUPWND Upper air station for domain winds -1 -1 ZUPWND Bottom and top of layer through which the domain scale winds are computed 1., 1000. 1., 1000.

IDIOPT4 Observed surface wind component switch 0 

### 0



### 3-9 Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044

# •

• • 

# • • 

# • • 

# • • 

# • • •

• 

# • 

# • • • • • • • 

# • • 

# • • • • 

# • 

# • 

# • • • 

# • • • • 

# • • • Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044 

# •

• • 

# • • CALMET Variable Description IDIOPT5 Observed aloft wind component switch LLBREZE Use lake breen > module? NBOX Number oflake breeze regions NLB Number ofstations in the region METBXID(NLB) Station IDs in the region CONSTB Neutral stability mixing height coefficient CONSTE Convective stability mixing height coefficient CONSTN Stable stability mixing height coefficient CONSTW Overwater mixing 

# t= coefficient


# - FCORIOL lute value of olis parameter lAVEZI Conduct spatial averaging? (1 = yes) MNMDAV Max search radius in averaging process (number ofgrid cells) HAFANG Half-angle ofupwind looking cone for averaging (degrees) ILEVZI Layers ofwind use in upwind averaging (between 1 and NZ) DPTMIN Minimum potential temperature lapse rate in the stable layer above the current convective mixing height DZZI Depth oflayer above current convective mixing height Value Included in CENRAP Protocol 0 F 0 0 0 1.41

0.15

2400

0.16 I.E-4

1 10

30

1

0.001

200



### 3-10 Value Trinity Will Use 0 F 0 0 0 1.41

0.15

2400

0.16 1.E-4

1 10

30

1

0.001

200 Notes

# •

• • • • • • • • • 

# • • • 

# • 

# • 

# • • • • 

# • • 

# • • • 

# • 

# • • 

# • • 

# • 

# • • • • 

# •



# • • 

# • • • CALMET Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes through which lapse rate is computed (m) ZIMIN Minimum overland mixing height (m) 50 50 ZIMAX Maximum overland mixing height (m) 3000 3000 ZIMINW Minimum overwater mixing height (m) 50 50 ZIMAXW Maximum overwater missing height (m) 3000 3000 ITPROG 3D temperature from observations or from MM5? 2 I Will use surface observations and precipitation data IRAD Type ofinterpolation (1 = lIr) I I TRADKM Temperature interpolation radius ofinfluence (lan) 36 36

NUMTS Max number of stations for temperature interpolations 5 5

IAVET Spatially average temperature? (1 = yes) I 1

TGDEFB Temperature gradient below the mixing height over water (KIm) -.0098 -0.0098

TGDEFA Temperature gradient above the mixing height over water (KIm) -.0045 -0.0045

JWATl Beginning land use categories over water 55 55

JWAT2 Ending land use categories for water 55 55

NFLAGP Precipitation interpolation flag (2 = 1I?) 2 2

SIGMAP 

# • • • Radius ofinfluence 50 50 

### 3-11 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 

### 063701.0044 

# •

# •

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

# • • • • 

# • • • • • 

# • • • • • 

# • • • 

# • • • • • 

# • • • • • • CALMET Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes for precipitation interpolation (km) CUTP Minimum precipitation rate cut off(mmIhr) 0.01 0.01



### 3-12 Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 

### 063701.0044

## 4. CALPUFF
Trinity will conduct a three-year CALPUFF analysis. TheCALPUFF model requires the input of meteorological data output by CALMET, source emissions data, receptor data, ozone and ammonia data, and model parameter settings. 

## 4.1 BART-ELIGmLE SOURCES As a Portland Cement plant, the Tulsa Plant is considered BART-eligible. In order to determine what specific units constitute the "BART-eligible" source, the construction dates ofthe each 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. PointID Description Construction Date 221 EP04 #2 discharge chute, separator, bucket elev. 1964 222 EP05 Raw mill pump hopper 1964 223 EP05 Raw mill pumphopper 1964 404 EP13 Kiln #2 feed system 1964 502 EPt5 Kiln No.2, screw conveyors, & 

### bucket elev. 1964 706 EP24 #2 bucket elev. and separator 1964 705 EP26 Finish mill #2 1964 704 EP25 Finish mill #2 feed belt 1964 602 EP17 Clinker cooler #2 1964 1101 

# R! EP27 CoaVcoak unloading EP28 CoaVcoak drag conveyor, bucket elev., crusher 1974 1974 EP28 CoaVcoak drag conveyor, bucket elev., crusher 1974

# ~:j\'\=l-l 

# r-fJh oJ'?SOI 4.2 SOURCE EMISSIONS DATA Trinity will include emissions ofS02, NO", and PMlOin the model. For source-specific analyses, CENRAP has suggested several sources ofemission data in the following priority:i

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

.... CEM data for 2002

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

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

.... Permit allowable emissions, or

....

Refined BART Modeling Protocol 063701.0044 Potential to emit

I CENRAP BARTModeling Guidelines. p. 7-/. 

### 4-1 Lafarge Building Materials Trinity Consultants

# i •• 

# • • 

# • • 

# • • • • • • • 

# • • • 

# • • 

# • • 

# • 4.2.1 PM10 EMISSIONS 

### PM lO 

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

### emissions from the kiln are speciated according to the Federal Land Managers' guidance to include the fol1owing: 2,3


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

### • Fine particulate matter (PMF, 

### considered PM<2.5) • Sulfates (S04) • Secondary organic aerosols (SOA) • Elemental carbon (EC) Based upon FLM guidance, PM IO 

### emissions from natural gas combustion consists ofprimarily coarse 

# • 

# • • 

# • • • 

# • • • • • • (25%) and fine (75%) particulate matter and emissions from fabric filters consists ofprimarily coarse (70%) and fme (30%). A summary 

### ofspeciated PM lO 

### e!Uissions is presented in Table 4-3. Full details ofPMIO 

### speciation can be found in Appendix A. 

# • 

# • • A detailed summary ofemission rates and source data is included in Appendix A. Table 4-2 provides a briefsummary ofthese emission rates. Emissions were calculated using emissions factors that were based on either AP-42 emission factors or baghouse manufacturer emission data. Based upon reported annual hours ofoperation 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 ofoperation. In addition, emissions were based on the maximum achievable production rate (or heat duty in the case ofcombustion units) over a one hour period. TABLE 4-2. ESTIMATED 24-HOUR AVERAGE EMISSION RATES Unit ID NOx 

### PM10 

### S02 Ibs/hr Ibs/hr Ibs/hr EP04 0.089 3.857 0.001 EP05 

## -­ 0.429 

## -­ EP13 

# -­ 0.557 

## -­ EP15 576.00 17.143 960.00 EP24 

## -­ 3.402 

# -­ ···EP26 

# -­ 0.429 

# -­

EP25 

### -­ 0.429 

# -­ EPl7


## -­ 4.607 

# -­ EP27 

## -­ 0.375 

## -­ EP28 

## -­ 0.321 

## -­

2 http://www2.nature.nps.gov/airlPermits/ect/ectCementKiln.cfin 3 As recommended by Visibility Improvement State and Tribal Association ofthe Southeast (VISTAS). 

# • 4-2 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • 

# • • • •TABLE 

### 4-3. SPECIATED PM 10 EMISSIONS 

# • 

# • • • • • • • 

# • • • • • • • • 

# • • • 

# • • • 

# • • 

# • • 

# • • • • 

# • • The base elevation ofeach the units is based upon visual inspection ofUSGS topographic maps. Source locations were based upon site inspections performed by Lafarge personnel. 

# • 

# • 

# • 

# • 

# • 4.3 RECEPTOR LOCATIONS The National Park Service (NPS) has electronic files for each Class I area available on their website 

# • 

# • 

# • containing the locations and elevations ofdiscrete Class I area receptors. The receptor files for Caney Creek Wilderness, Hercules-Glades Wilderness, Upper Buffalo Wilderness, and Wichita Mountains Wildlife Refuge will be downloaded from the NPS website, converted into the LCC NAD27 

# • • • • 804 80A 

### PM (coarse) PM (f"me) EC 

### Ibs/hr Ibs/hr Ibs/hr Ibs/hr Ibs/hr Unit ID

EP04 0.964 2.893 

## -- 

# -- -­ EP05 0.300 0.129 

# -- -­ EP13 0.390 0.167 

## -- -- 

# -­ EP15 1.149 1.320 .. 12.91 1.766 

## -­ EP24 2.381 1.021 

# -- -- -­ EP26 0.300 0.129 

## -- 

# -- -­ EP25 0.300 0.129 

# -- 

# -- -­ EP17 3.225 1.382 

## -- 

## -- -­ EP27 0.263 0.113


## -- -- 

# -­ EP28 0.225 0.096 

## -- 

# -- -­ 4.2.2 STACKPARAMETERS 

### Table 4-4 provides a summary ofsource exhaust parameters. TABLE 

### 4-4. POINT SOURCE STACK PARAMETERS 

### Base Emission Latitnde Longitnde Exit Stack Stack Exit (degrees) (degrees) unit Elevation Velocity Height Diameter Temperatnre (m) (ft/s) (feet) (feet) (Fahrenheit) EP04 36.194777 95.811055 216 11.94 100 4 180 EP05 36.19475 95.810722 216 8.29 55 1.6 70 EP13 36.195472 95.811916 216 9.55 97 1.7 200 EP15 36.195833 95.811777 216 17.83 60 6.9 500 EP24 36.19425 95.811222 216 10.53 100 4 200 36.194166 95.811166 EP26 216 33.16 64 0.8 70 EP25 36.194166 95.811166 216 33.16 55 0.8 200 EP17 36.19375 95.811805 216 3.56 53 8.0 200 36.194027 95.811027 216 4.64 7 2.0 70 EP27 EP28 36.194138 95.811333 216 4.91 135 1.8 70



### 4-3 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •

• • • 

# • • • • 

## • • • • • 

## • • • 

## • • • • 

## • •



# • 

# •
• 

# • 

# •
• • • 

# • •
• • 

# • • • • • •
• 

# • • • • 

# • •
• 

# •
• 

# • •
• 

# •
• 

# •
• 

# •
• 

# • 

# • •**~~~.ili:~6rup** AugU51 22,2003



## • 

# ..

.. • •



# • • • • • • 

# • 

# • 

# • • 

# • • 

# • • 

# • • 

# • • • • • 

# • • • • 

# • • • •

• 

# • • 

# • • • • • • •

• 

# • • 4-4 

### Lafarge Building Materials Refined BART Modeling Protocol 

# • • projection, and incorporated into the CALPUFF model. Receptor locations for these Class I areas are shown in Figures 4-1 to 4-4. FIGURE

### 4-1. CANEY CREEK WILDERNESS RECEPTORLOCATIONS Caney Creek**Wilderne~s** 80 Rt:ceptaTs 1.J SFS Boundary Source:

FS National Co.erag& tile:



# .. .. 

# .. 

# .. • • • • 

## • • • • • 

# .. • • • 

# .. • 

# .. • 

## • • 

# .. • • • •



# .. • 

# .. • • • •

# ~ • • • • 

# .. • 

## •

FIGURE 

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



# • 

# • •

Hercules-Glades Wilderness 80 RSCIlplors USFS Boundary Source:

FS National Coverage file: NRIS - ,ALp group Corvallis, OR Aug usl 22.2003 

### Trinity Consultants 063701.0044

# •



# • • • • 

# • • • 

# • 

# • • • 

# • • 

# • • • 

# • • • • • • 

# • • • • • • 

# • • 

# • • • • • • • • • • 

# • • • 

# • • • • 

# • • • • • • • 

# • • • • 

# • • • • 

# • • • • • 

# • 

# • 

# •
FIGURE 

### 4-3. UPPER BUFFALO WILDERNESS RECEPTOR 

### LOCATIONS 

# • • • 

# •• • 

# • • Uppe r Buffalo Wilderne ss 72 Receptors USFS Boundary Source:

FS National Coverage file: NRIS· ALP group Corvallis, OR August 22,2003 

# • 

# • • • • • • • • 

# • • 

# • • •

•



# • 

# • 

# • • 

# • 

# • • 

# • 4-5 Lafarge Building Materials Trinity Consultants Refined BARTModeling Protocol 063701.0044" 

### '",



# • •

FIGURE 

### 4-4. WICHITAMOUNTAINS WILDERNESS RECEPTOR LOCATIONS 

# • 

# • • • 

# • Wichita Mountains Wilderness ""thin Wichita Mountains NWR 59 Receptors FWS Wilderness Boundary Source: 

# • Region2GIS 500 Gold A.B., SW PO 80.1306 Albuquerque. Nt/! 87103 FWS NWR Boundary Source: 

# • 

# • • • • • 

# • • 

# • • • • • 

# • 4.4 BACKGROUND OZONE 

# • • • • •• • • 

## • hltp :lIwww.fws.govldatalr2gislboundary.hlml SeplBm ber 1

### e, 2003



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

# • • 

# • S02 and NOli. to sulfates (S04) and nitrates (N03)' CALPUFF can use either a single background value representative ofan area orhourly ozone data from one or more ozone monitoring stations. 

# • • • • CENRAP recommends either developing background ozone estimates from ambient monitors located within the particular domain being modeled or developing background ozone estimates from CENRAP's most recent CMAQ or CAMx simulation for the 2002 base year. Trinity proposes the use ofa 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 ofammonium sulfates and 

# • 

# • ammonium nitrates. CENRAP recommends developing background ammonia estimates from CENRAP's most recent CMAQ or CAMx simulation for the 2002 base year. Since 

### CMAQ/CAMx modeled and observed monthly averaged ammonia concentrations exhibit wide spatial variability, 

# • 

# • • CENRAP recommends obtaining separate monthly-averaged ammonia concentrations from CMAQ or CAMx for the CENRAP north, central and south modeling domains, respectively. These would then be used as input to CALPUFF. Since the data from CENRAP's CMAQ and CAMx simulations is not readily available, Trinity proposes the use ofa conservative monthly background concentration 

# • 4-6 Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044 

# • • 

# • 

# •

# ..,



# • 

# • • 

# • 

# • • • 

# • • • 

# • • 

# • • 

# • • 

# • • 

# • • • 

# • 

# • • • • • 

# • • • 

# • 

# • • • • 

# • 

# • • 

# • • • 

# • • of3 ppb. This background concentration is the value included in CENRAP's protocol as a default background value for the CENRAP region. 

## 4.6 CALPUFF MODEL CONTROL PARAMETERS Table 4-4 provides a listing ofthe CALPUFF parameters that Trinity will use in Lafarge's modeling analysis. In addition to the parameters used for Lafarge's modeling, the table also lists CENRAP's recommended parameters for comparison. Incases where a parameter is different than what CENRAP recommended, a short explanation as to the difference is provided. TABLE 

### 4-4. SUMMARY OF 

### CALPUFF INpUTS CALPUFF

Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes METRUN All model periods in met files will be run 0 0

IBYR Starting year 2001 Appropriate met year Years 2001, 2002, 2003 IBMO Starting month 1 

# =:i IBDY Starting day 

### I 

### 1 IBHR Starting hour 1 1 XBTZ Base time zone (6
= CST) 6 6 IRLG Length ofrun 8760 8760 NSPEC Number of MESOPUFF II chemical species 10 10

NSE Number of chemical species to be emitted 8 6 S02, S04, NOx, PMC, PMF, SOA lTEST Program is executed after SETUP phase 2 . 2

MRESTART Do not read or write a restart file during run 0 0

NRESPD File written only at last period 0 0 METFM CALMETbinary file (CALMET.MET) I 1
Averaging time in minutes 60 60 PGTIME PG Averaging 60 60 

### 4-7 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •

• 

# • • 

# • • • 

# • • • • 

# • 

# • • • 

# • • 

# • • 

# • • • • • • • • • 

# • • • • • 

# • • 

# • • • 

# • CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes. time in minutes MGAUSS Gaussian distribution used in near field I I

MCTADJ Partial plUme path terrain adjustment 3 3 MCTSG Sub-grid-scale complex terrain not modeled 0 0

MSLUG Near-field puffs not modeled as elongated 0 0

MTRANS Transitional plume rise modeled 1 1 MTIP tip ash used 1 1 MSHEAR (0, 1) Vertical wind shear (not modeled, modeled) 0 0

MSPLIT Puffs are not split 0 1 Allowed for puff splitting due to significant sourcereceptor distances MCHEM MESOPUFFII chemical parameterization scheme 1 1

MAQCHEM Aqueous phase transfonnation not modeled 0 0

MWET Wet removal modeled 

### I 1 MDRY Dry deposition modeled I 1 MDISP PG dispersion coefficients 3 3 MTURBVW Use both cry and cr

w from PROFILE.DAT to compute cry and cr z (n/a) 3 

### 3

MDISP2 PG dispersion 

### 3 3 

# • • 

# • 

# • 4-8 Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044 

# •

# •

• • 

# • • 

# • • 

# • 

# • 

# •

• 

# • • • 

# • • 

# •

• • • 

# • • 

# • • • • •

• • • • 

# • 

# • • • • 

# • • 

# • 

# • Lafarge Building Materials Trinity Consultants 

# • 

# • • • CALPUFF Value Included in Variable Description CENRAP Protocol coefficients MROUGH PG O'y and O'z not 

### 0 adjusted for roughness MPARTL No partial plume 

### I penetration of elevated inversion MTINV Strength of 

### 0 temperature inversion computed from default gradients MPDF PDF not used for 

### 0 dispersion under convective conditions MSGTIBL Sub-grid TIBL 

### 0 module not used for shoreline MBCON Boundary 

### 0 concentration conditions not modeled MFOG Do not configure 

### 0 for FOG model output MREG Technical options 

### 1 must conform to USEPALong Range Transport (LRT) guidance CSPEC 

### CENRAP Output Group Species Modeled Emitted S02 1 

### 1 S04 1 

### 0 

### NOX 1 1 HN03 1 

### 0 N03 1 

### 0 NH3 

### 1 1 PMC 

### 1 1 PMF 

### 1 1 Dry Deposition 1 

### 2 1 1 2 1 

### 2 2 

### 4-9 Value Trinity Will Use Notes 

### 0

I



### 0

0

0

0

0

1



### Lafarge Output Group Species Modeled S02 1 S04 1 

### NOX 1 HN03 1 N03 

### 1

- NH3 1 PMC 1 PMF 1 

### Dry Emitted Deposition 1 

### 1 1 

### 0 0 0 1 

### 1 1 2 1 

### I 

### 2 1 2 

# 2"

Refined BART Modeling Protocol 063701.0044

# •



# • • • 

# • • • • 

# • • 

# • • • 

# • • 

# • • 

# • • 

# • 

# • • • • 

# • • 

# • • 

# • • 

# • • • • • • • 

# • • 

# • • • • • CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes EC I I 2 EC 1 0 2 SOA I 1 2 SOA 1 1 2 PMAP Map projection Number ofX grid cells in meteorological grid UTM LCC 

### NX 66 256 Based on refmed grid parameters as specified in CALMETruns NY Number ofY grid cells in meteorological grid 66 155 Based on refmed grid parameters as specified in CALMETruns NZ Number ofvertical layers in meteorological grid 10 10

DGRIDKM Grid spacing (km) 6 3.0 0,20,40,80, 160, 320,640, 1200, 4000 Refmed grid size ZFACE Cell face heights in meteorological grid(m) 0, 20, 40, 80, 16O, 320,640,1200,2000, 3000,4000 XORIGKM Reference X coordinate for SW comer ofgrid cell ofmeteorological grid (km) 5 

# ° Appropriate for domain (refmed grid)

YORIGKM Reference Y coordinate for SW comer ofgrid cell ofmeteorological grid (km) 3327 

# ° Appropriate for domain (refmed grid)

IUTMZN UTM zone of coordinates (NAD83) (use only if PMAP=UIM) 12 14 N/A; LCC coordinates used

IBCOMP X index oflower left comer ofthe computational grid I 

### I

JBCOMP Y index oflower left comer ofthe computational grids I 

### I

IECOMP X index ofupper right comer ofthe 66 256 Appropriate for domain (refined grid) 

### 4-10 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# •

• 

# • •

• 

# • 

# • • • • • • 

# • • • 

# • • • 

# • • • 

# • • • 

# • • 

# • 

# • • 

# • 

# •
• • 

# • • • 

# • • • •



# • 

# • 

# • CALPUFF Variable 

# •

Description Value Included in CENRAP Protocol Value Trinity Will Use Notes computational grid JECOMP Y index ofupper right comer ofthe computational grid 

### 66 155 Appropriate for domain (refmed grid) LSAMP Sampling grid is not used F F IBSAMP X index oflower left comer of sampling grid 

### 1 1

JBSAMP Y index oflower left comer of sampling grid 

### 1 1

IESAMP X index ofupper right comer of samplingDid 

### 66 256 Appropriate for domain (refined grid) JESAMP Y index ofupper right comer of sampling grid 

### 66 155 ':Appropriate for domain (refmed grid) MESHDN Nesting factor of sampling grid 

### 1 

### 1 ICON Output file CONC.DAT containing concentrations is created 1 

### 1

IDRY Output file DFLX.DAT containing dry fluxes is created 

### 1 1

IWET



### . Output file WFLX.DAT containing wet fluxes is created 

### 1 1

!VIS Output file containing relative humidity data is created 

### 1 1

LCOMPRS Perform data compression in output file T T

IMFLX Do not calculate mass fluzes across specific 

### 0 0



### 4-11 Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044

# •



# • • • • e e. 

# • • 

# • • 

# • • 

# • e 

# • • 

# • e 

# • • e 

# • • 

# • • • • 

# • • 

# • • • 

# • 

# • • • 

# • • CALPUFF Value Included in Variable Description CENRAP Protocol boundaries 

### IMBAL Mass balances for 0 each species not reported hourly ICPRT Print concentration 1

, fields to output list file IDPRT Do not print dry 0 flux fields to output list file IWPRT Do not print wet 0 flux fields to output list file ICFRQ Concentration 1 fields are printed to output list file every hour IDFRQ Dry flux fields are 1 printed to output list file every 1 hour IWFRQ Wet flux fields are 1 printed to output list file every 1 hour IPTRU Units for line 

### 3 printer output are 

### in glm3 for concentration and g/m2/s for deposition IMESG Messages tracking 

### 2 the progress ofrun written to screen Logical value for F debug output IP First uffto track 

### 1 NP puffs 

### 1 NNI Meteorological 1 period to start output NN2 Meteorological 10 Value Trinity Will Use 0

1

0

0


### I

I

1



### 3

2

F 1 1
1

10 Notes



# • • • 

# • 4-12 

### Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044 

# •

# •

• 

# • • • • 

# • • 

# • • • • 

# • • 

# • • • • 

# • • 

# • • • • 

# • • • 

# • • • • • 

# • • 

# • • • 

# • CALPUFF Variable 

# •

Description Value Included in CENRAP Protocol Value Trinity Will Use Notes period to end output NHILL Number ofterrain features 0 0 NCTREC Number ofspecial 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 (Ian) 0 0

YCTDMKM Y origin ofCTDM system relative to CALPUFF coordinate system (Ian) 0 0



### S02 Diffusivity 0.1509 0.1509 Alpha star 1000 1000 Reactivity 8 8 Mesophyll resistance 0 0 Henry's Law coefficient 0.04 0.04 

### NOX Diffusivity 0.1656 0.1656 Alpha star 1 1 Reactivity 8 8 

# • 

# • • Mesophyll resistance 5 5 

### 4-13 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • •

# •



# • 

# • • 

# • • • • • • • • • 

# • 

# • • 

# • • • 

# • • 

# • • • • 

# • • • 

# • • 

# • • 

# • • • 

# • • 

# • • 

# • 

# • 

# .' CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes Henry's Law coefficient 3.5 3.5 HN03 Diffusivity 0.1628 0.1628 Alpha star I I Reactivity 18 18 Mesophyll resistance 0 0 Henry's Law coefficient 8.e-8 8.e-8 S04-2 Geomatric mass mean diameter of S04-2 (!-lm) 0.48 0.48

N03­ Geometric mass mean diameter of N03- (].lm) 0.48 0.48

PMC Geometric mass mean diameter of PMC (!-lm) 6 6

PMF Geometric mass mean diameter of PMF (!-lm) 0.48 0.48

EC Geometric mass mean diameter of EC (!-lID) 0.48 0.48

SOA Geometric mass mean diameter of SOA (!-lm) 0.48 0.48

RCUTR Reference cuticle resistance (s/cm) 30 30 RGR Reference ground resistance (s/cm) 10 10 REACTR Reference pollutant reactivity 8 8 NINT Number ofparticle size intervals for effective particle deposition velocity 9 9

IVEG Vegetation in nonirregated areas is active and unstressed I I

S02 

# • 

# • Scavenging coefficient for 3.2IE-05 3.2IE-05


### 4-14 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes liquid precipitation (S·I) Scavenging coefficient for frozen precipitation (S·I) 0 0

S04-2 Scavenging coefficient for liquid precipitation (S·I) I.OE-04 I.OE-04

Scavenging coefficient for frozen precipitation (S'I) 3.0E-05 3.0E-05

RN03 Scavenging coefficient for liquid precipitation (S·I) 6.0E-05 6.0E-05

Scavenging coefficient for frozen precipitation (S·I) 0 0

N03­ Scavenging coefficient for liquid precipitation

.(S·I) 1.0E-04 I.OE-04

Scavenging coefficient for frozen precipitation (S·I) 3.0E-05 3.0E-05

NH3 Scavenging coefficient for liquid precipitation (S·I) 8.0E-05 NA

Scavenging coefficient for frozen precipitation (S·I) 0 NA

PMC Scavenging coefficient for liquid precipitation (S·I) 1.0E-4 1.0E-4

Scavenging 3.0E-05 3.0E-05 

### 4-15 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# :.'

• 

# ,.

•

Ie



# .'**~.** 

# • 

# • o


# . e

'- ' 

# • 

# • 4-16 

# •• 

# • • '



### Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044

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

Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes coefficient for frozen precipitation (S-I) PMF Scavenging coefficient for liquid precipitation (S-I) 1.OE-05 l.OE-05

Scavenging coefficient for frozen precipitation (S-I) 3.0E-05 3.0E-05

EC Scavenging coefficient for liquid precipitation (S-I) 1.OE-04 l.OE-04

,Scavenging coefficient for frozen precipitation (s-I) 3.0E-05 3.0E-05

OC Scavenging coefficient for liquid precipitation (S-I) l.OE-04 l.OE-04

" Scavenging coefficient for frozen precipitation (S-I) 3.0E-05 3.0E-05

MOZ Read ozone background concentrations from ozone.dat file (measured values) 1 1

BCK03 Background ozone concentration (ppb) 12*40 12*40

BCKNH3 Background ammonia concentration (ppb) 12*3 12*3

RNITE1 Nighttime S02 loss rate is %/hour 0.2 0.2 RNITE2 Nighttime NOX loss rate is %/hour 2 2

# •• CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes RNITE3 Nighttime HN03 loss rate is %lhour 2 2 MH202 Background H202 concentrations 1 

# ° Need to choose

# ° in order to use monthly background value BCKH202 Background monthly H202 concentrations 1 12*1

BCKPMF Fine particulate concentration for SOA option (llg 

### /m3) 1 1

OFRAC Organic fraction of fme particulate for SOA option .2 0.15,0.15,0.2,0.2,0.2, 0.2,0.2,0.2,0.2,0.2, 0.2,0.15 Irrelevant, since MCHEM not equal to 

### 4 VCNX VOC/NOXratio

! for SOA option 50 50 SYDEP Horizontal size of a puffin meters beyond which the time dependant dispersion equation ofHeffter is used 550 550

MHFTSZ Do not use Heffier formulas for sigma z 

# ° 0

JSUP Stability class used to determine dispersion rates for puffs above boundary layer 5 5

CONK] Vertical dispersion constant for stable conditions 0.01 0.01

CONK2 Vertical dispersion constant for neutral/stable conditions 0.1 0.1

TBD Use ISC transition point for determining the transition point 0.5 0.5 

# • • 

# • • 

# • 

# • 

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

# • 

# •

# •

• • 

# • • 

# • • • • 

# • 

# • • • 

# • • • • 

# • • • • 

# • • 

# • • • 

# • • 

# • • • • • 

# •

• • • • 

# • CALPUFF Variable 

# •

Description Value Included in CENRAP Protocol Value Trinity Will Use Notes between the Schulman-Scire to Huber-Snyder Building Downwash scheme IURBI Lower range of land use categories for which urban dispersion is assumed 10 10

IURB2 Upper range of land use categories for which urban dispersion is assumed, 19 19

ILANDUlN Land use category for modeling domain 

## * *

XLAIIN Leafarea index for modeling domain 

## * * ZOIN Roughness length in meters for modeling domain 

## * *

ELEVIN Elevation above sea level 

## * * XLATIN North latitude of station in degrees 

## - XLONIN South latitude of station in degrees 

### - 

# - meter 

### in meters I 10 10 ISIGMAV Sigma-v is read for lateral turbulence data 1 1

IMIXCTDM Predicted mixing heights are used 0 0 XMXLEN Maximum length ofemitted slug in meteorological grid units 1 1

XSAMLEN 

# • • Maximum travel distance ofslug or 10 10


### 4-18 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • 

# •

# •

• • •



# • • 

# • • • 

# • • 

# • • • • • 

# • • •
• 

# • • • 

# • • • • • • • 

# • 

# • • • 

# • • 

# • • • • 

# • 

# • • CALPUFF Variable 

# •

Description Value Included in CENRAP Protocol Value Trinity Will Use Notes puffin meteorological grid units during one sampling unit MXNEW Maximum number ofpuffs or slugs released from one source during one time step 60 60

MXSAM Maximum number ofsampling steps during one time step for a puffor slug 60 60

NCOUNT Number of iterations used when computing the transport wind for a sampling step that includes ,transitional plume rise 2 2

SYMIN Minimum sigma y in meters for a new puffor slug I I

SZMIN Minimum sigma z in meters for a new puffor slug I I

SVMIN Minimum lateral turbulence velocities (mls) 0.5 0.5

SWMIN Minimum vertical turbulence velocities (mls) 0.20,0.12, 0.08, 0.06, 0.03,0.016 0.20,0.12, 0.08, 0.06,0.03,0.016 WSCALM Minimum noncalm wind speeds (mls) 0.5 0.5

XMAXZI Maximum mixing height (m) 3000 3000 XMINZI Minimum mixing height (m) 20 20 SL2PF Maxmium cry/puff length 10 

### 10


### 4-19 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044

# ~­



# • 

# • • • 

# • • 

# ., • • • 

# • 

# • Lafarge Building Materials Trinity Consultants Refilled BART Modeling Protocol 063701.0044

CALPUFF

Variable Description PLXO Wind speed power-law exponents WSCAT Upper bounds of 1st 5 wind speed classes PGGO Potential temp gradients PG E &

N CDIV Divergence criterion for dw/dz (lis) PPC Plume path coefficients (only ifMCTADJ = 3) NSPLIT Number ofpuffs when split IRESPLIT Hours when puff is eligible to split ZISPLIT Previous hours minimum mixing height,m ROLDMAX Previous max Olixing. height/current height ratio, must be less than this value to allow puff to split NSPLITH Number ofpuffs resulting from a split SYSPLITH Minimum sigma-y ofpuffbefore it may split SHSPLITH Minimum puff elongation rate from wind shear before puffmay split CNSPLITH Minimum species concentration before a pluffmay Value Included in CENRAP Protocol 0.07,0.07, 0.1 0,0.15, 0.35,0.55 1.54,3.09,5.14,8.23, 10.80 0.020, 0.035

0.01
0.5,0.5,0.5,0.5, 0.35,0.35 3 1900 100

0.25

5

1.0

2.0

I.OE-07



### 4-20 Value Trinity Will Use 0.07,0.07,0.10, 0.15,0.35, 0.55 1.54,3.09,5.14, 8.23,10.80 0.020, 0.035

0.01
0.5, 0.5, 0.5, 0.5, 0.35,0.35 3 Hour 17 100

0.25

5

1.0

2.0

I.OE-07 Notes

Should be by hour of day-1900 is hour 17

,

. , .**.~**CALPUFF Variable Description Value Included in CENRAP Protocol Value Trinity Will Use Notes split EPSSLUG Criterion for SLUG sampling 1.0E-04 1.0E-04 EPSAREA Criterion for area Trajectory step length for numerical site algorithm 1.0E-06 1.0E-06 DSRlSE 1.0 1.0

NPTl Number ofpoint sources with constant stack parameters or variable emission rate scale factors Varies by scenario 10 Appropriate for Lafarge

IPTU Units for point source emission rates are gls 1 3 Used different units (3 = lblhr) NSPTl Number ofsourcespecies combinations with variable emissions scaling factors 

# - -

,. NPT2 Number ofpoint sources with variable emission parameters provided in extemal file 

# - -

MISC Other point source inputs include stack height, stack diameter, exit temperature, exit velocity, downwash flag and emissions by species 

# - -

NARI Number of polygon area sources Varies by scenario 0 None modeled

IARU Units for area source emission 1 1


### 4-21 

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

# •



# • • • • • 

# • • • • 

# • • • 

# • • • 

# • • • • • • • • • • • • • • • • • • • • 

# • • 

# • • 

# • • • • CALPUFF Variable**~** rates are g/m2/s NSARI Number ofsource 

# - species combinations with variable emissions scaling factors NAR2 Number of 

# - buoyant polygon area sources with variable location and emission parameters NLN2 Number of 

# - buoyant line sources with variable location and emission parameters NLINES Number of 

# - buoyant line sources ILNU Units for line source emission rates in 

### 't! s NSLNI Number ofsource­ 

# - species combinations with variable emissions scaling factors MXNSEG Maximum number 

# - ofsegments used to model each line NLRISE Number of 

## - distance at which transitional rise is computed XL Average line 

# - source length (m) HBL Average height of 

## - line source height (m) WBL Average building 

# - width(m) 

## -

-

0



## -



# -

-



# -

-

-

-

- Value Trinity Will Use Notes

None modeled



### 4-22 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044CALPUFF Variable Description WML Average line source width (m) DXL Average separation between buildings (m) FPRlMEL Average buoyancy parameter (m4/s3) NVLl Number ofvolume sources IVLU Units for volume source emission rates in grams/second NSVLl Number of sourcespecies combinations with variable emissions scaling factors IGRDVL Gridded volume source data is not used VEFFHT Effective height of emissions (m) VSIGYI Initial sigma-y value VSIGZI Initial sigma-z value NREC Number ofnongridded receptors Value Included in CENRAP Protocol



## -

-

-

-

-

-

-

-

-

- 5630 Value Trinity Will Use



# -

-

- 0



## -

-

-

-

-

- 271 Notes

None modeled



## -

Receptor data provided in Appendix 

# ..• •• /.

# ~.

4-23 Lafarge Building Materials Trinity Consultants 

# '. • 

# • • •

Refined BART Modeling Protocol 063701.0044

# •



# • • •
5. CALPOST 

# • 

# • • • • • Trinity will conduct a three-year CALPOST analysis to determine the change in light extinction. caused by Lafarge's BART-eligible sources when compared to a natural background. The CALPOST model requires the input ofconcentration data output by CALPUFF. 

## 5.1 LIGHT EXTINCTION ALGORITHM Trinity will utilize EPA's currently approved algorithm for reconstructing light extinction (as opposed to the new equation for reconstructing light extinction recommended by the IMPROVE Steering Committee). The light extinction equation is provided below. 

### b ext 

## = 3*f(RH)*[(NH4)2 S04] 

## + 3*f(RH)*[NH4N03] 

## + 4*[OC] 

## + l*[PMrJ 

## + 0.6 * [PMc ] 

## + lO*[EC] 

## + bRay 

# • • • • • 

# • • • • • The algorithm will be used to calculate the daily light extinction attributable to Lafarge's BART­ eligible sources and light extinction attributable to a natural background. The change in deciviews based on the source and background light extinctions was evaluated using the equation below. 

### /), dv 

## = 10* In b ext, background 

### + bext, source 

# 1 [ bext, background 

# • 

# • • • 

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

# • 

# • and background using monthly average relative humidity adjustment factors. Trinity will use monthly Class I area-specific relative humidity adjustment factors based on the centroid ofthe Class I area as included in Table A-3 of EPA's 

### Guidance for Estimating Natural Visibility Conditions Under the Regional Haze Program. The factors for the relevant Class I areas in the analysis are provided in Table 5-1. TABLE 

### 5-1. MONTHLY HUMIDITY FACTORS

Class I Area Jan Feb Mar Apr May Jun Jut Aug Sep Oct Nov Dec Caney Creek 3.4 3.1 2.9 3.0 3.6 3.6 3.4 3.4 3.6 3.5 3.4 3.5 Hercules-Glades 3.2 2.9 2.7 2.7 3.3 3.3 3.3 3.3 3.4 3.1 3.1 3.3

Upper Buffalo 3.3 3.0 2.7 2.8 3.4 3.4 3.4 3.4 3.6 3.3 3.2 3.3 Wichita Mountains 2.7 2.6 2.4 2.4 3.0 2.7 2.3 2.5 2.9 2.6 2.7 2.8 

# • • • • • 

# • • 

# • • 

# • • • • 

# • • 5-1 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 

### 063701.0044 

# • 

# •

# •

• 

# • • • 

# • • • • 5.3 NATURALBACKGROUND For Caney Creek and the Wichita Mountains, Trinity will use default west natural background values. Hercules-Glades and Upper Buffalo will be processed using default east natural background values. 4 The annual average concentrations are included in Table 2-1 ofEPA's 

### GuidanceforEstimating Natural Visibility Conditions Under the RegionalHaze Program and are provided in Table 5-2. TABLE 5-2. DEFAULT AVERAGE ANNUAL NATURAL BACKGROUND LEVELS

Component Average Annual Natural Background West (llg/m3) Average Annual Natural Background East (llg/m3) Ammonium Sulfate 0.12 0.23 Ammonium Nitrate 0.10 0.10 Organic Carbon Mass 

### OA7 lAO

Elemental Carbon 0.02 0.02 Soil 0.50 0.50 Coarse Mass 3.0 3.0 

### 5.4 EVALUATING BART-EXEMPTION Trinity will compare the 2001 through 2003 daily**~dv** values output by CALPOST to a contribution threshold of0.5**~dv.** Ifthe 98 th 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 98 th percentile value 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 ofthe CALPOST parameters that Trinity will use in Lafarge's modeling analysis. In addition to the parameters used for Lafarge's modeling, the table also lists CENRAP's recommended parameters for comparison. In cases where a parameter to be used is different than what CENRAP recommended, a short explanation as to the difference is provided. 

# • • TABLE 5-3. SUMMARY OF CALPOST INPUTS CALPOST

Variable Description Value Included in CENRAP Protocol Value Trinity will use Notes ISYR Starting year 2001 Appropriate met year Years 200 I, 2002, 2003 ISMO Starting month I 1 ISDY Starting day I 1

4 Per comments received from EPA Region 7 to Trinity Consultants' Kansas City office, May 22, 2006. 5-2 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • • • • • • • • • • 

# • •

# •



# • • 

# • • 

# • • • 

# • • • 5-3 Lafarge Building Materials Trinity Consultants CALPOST

Variable Description Value Included in CENRAP Protocol Value Trinity will use Notes ISHR Starting hour 0 0 All CALPUFF periods will be included NPER Number ofperiods to process 8760 8760 NREP Process every hour of data? 1 

### = yes 1 1 ASPEC Process species for visibility VISIB VISIB ILAYER Layer/deposition code; 1 for CALPUFF concentrations 1 1

A Scaling factor, slope 0 0 B Scaling factor, intercept 0 0 LBACK Add hourly background concentrations of fluxes? F F

LG Process gridded receptors? F F LD Process discrete receptors? T T LCT Process complex terrain receptors? F F LDRING Report receptor ring results? F F NDRECP Select all discrete receptors -1 Varies As appropriate for Class I area being analyzed IBGRID X index ofLL comer ofreceptor grid 

### -1 -1 JBGRID Y index ofLL comer ofreceptor grid 

### -1 -1 IEGRID X index ofUR comer ofreceptor grid 

### -1 -1 JEGRID Y index ofURcorner ofreceptor grid -1 -1 NGONOFF Number ofgridded receptor rows 0 0 NGXRECP Exclude specific gridded receptors, Yes =0 0 0



### Refmed BART Modeling Protocol 063701.0044 

# • 

# • 

# •

# •

• 

# • • • 

# •

• • • • 

# • • • 

# • • • • • • • 

# •

• 

# • • 

# • • • • • • • 

# • • 

# • 

# • • 

# • • •

• • 

# • • 

# • CALPOST Variable Description Value Included in CENRAP Protocol Value Trinity will use Notes RHMAX Maximum RH% used in particle growth curve 95 • 95

LVS04 Compute light extinction for sulfate? T T LYN03 Compute light extinction for nitrate? T T LVOC Compute light extinction for organic carbon? T T

LVMPC Compute light extinction for coarse particles? T T

LVMPF Compute light extinction for fine particles? T T

LVEC Compute light extinction for elemental carbon? T T

LVBK Include background in extinction calculation? T T SPECPMC Coarse particulate species PMC PMC SPECPMF Fine particulate species Extinction efficiency for coarse particulates PMIO PMF Notation difference EEPMC 0.6 0.6 EEPMF Extinction efficiency for fine particles? 1.0 1.0 EEPMCBCK Extinction efficiency for coarse part. Background 0.6 0.6

EES04 Extinction efficiency for ammonium sulfate 3.0 3.0 EEN03 Extinction efficiency for ammonium nitrate 3.0 3.0 EEOC Extinction efficiency for organic carbon 4.0 4.0 EESOIL Extinction efficiency for soil 1.0 1.0

, EEEC Extinction efficiency for elemental carbon 10.0 I 10.0

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

# -

CALPOST Variable Description Value Included in CENRAP Protocol Value Trinity will use Notes MVISBK Method 6 for background light extinction 6 6

BEXTBTBK Background extinction for MVISBK=1 12 Not Used Not necessary since MVISBK=6 RHFRAC % ofparticles affected byRH 10 Not Used Not necessary since MVISBK=6 RHFAC Extinction coefficients for modeled and background hygroscopic species computed using EPA (2003) monthly RH adjustment factors Depends on Class I Area For Caney Creek: 3.4,3.1,2.9,3.0, 3.6,3.6,3.4,3.4, 3.6,3.5,3.4,3.5 Hercules-Glades: 3.2,2.9,2.7,2.7, 3.3,3.3,3.3,3.3, 3.4,3.1, 3.1, 3.3 Upper Buffalo: 3.3,3.0,2.7,2.8, 3.4,3.4,3.4,3.4, 3.6,3.3,3.2,3.3 Wichita Mountains: 2.7, 2.6, 2.4, 2.4, 3.0,2.7,2.3,2.5, 2.9,2.6,2.7,2.8 As appropriate for Class I area

BKS04 Background sulfate extinction coeff

­ east 0.23 0.23 Used for HerculesGlades and Upper Buffalo BKN03 Background nitrate extinction coeff

­ east 0.10 0.10 Used for HerculesGlades and Upper I Buffalo BKPMC Background sulfate extinction coeff

­ east 3.00 3.00 Used for HerculesGlades and Upper Buffalo BKSOC Background organic carbon extinction coeff- east lAO 1.40 Used for HerculesGlades and Upper Buffalo BKSSOIL Background soil extinction coeff

­ east 0.50 0.50 Used for HerculesGlades and Upper Buffalo BKSEC Background elemental carbon extinction coeff

- east 0.02 0.Q2 Used for HerculesGlades and Upper Buffalo 

# '. • • • 

# • • • 

# • 

# • • 5-5 Lafarge Building Materials Trinity Consultants 

# • Refined BART Modeling Protocol 

### 063701.0044 

# •

# •

• • • • • • 

# • • 

# • • • • • 

# • • 

# • • • • 

# • 

# • • 

# • • 

# • • • • • • • • • 

# • • • 

# • •

CALPOST Variable Description BKS04 Background sulfate extinction coeff

- west BKN03 Background nitrate extinction coeff

- west BKPMC Background sulfate extinction coeff

- west BKSOC Background organic carbon extinction coeff

- west BKSSOIL Background soil extinction coeff

- west BKSEC Background elemental carbon extinction coeff-west BEXTRAY Extinction due to Rayleigh scattering ClIMm) LDOC Print documentation image? IPTRU Print output units for concentrations and for deposition LlHR Report 1 hr averaging times L3HR Report 3 hr averaging times L24HR Report 24 hr averaging times LRUNL Report run-length averaging times LT50 Top 50 table LTOPN Top N table NTOP Number ofTop-N values at each receptor ITOP Ranks ofTop-N values at each receptor LEXCD Threshold exceedances counts Value Included in CENRAP Protocol 0.12

0.10

3.00

0.47

0.50

0.02

10.0

F 3

F F T F F F 4 1,2,3,4 F Value Trinity will use 0.12

0.10

3.00

0.47

0.50

0.02

10.0

F 1

F F T F F F 4 1,2,3,4 F Notes Used for Wichita Mountains and Caney Creek Used for Wichita Mountains and Caney Creek Used for Wichita Mountains and Caney Creek Used for Wichita Mountains and Caney Creek Used for Wichita Mountains and Caney Creek Used for Wichita Mountains and Caney Creek

Units preference



# • • 

# • 5-6 Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044 

# • •

# -

CALPOST Variable Description Value Included in CENRAP Protocol Value Trinity will use Notes THRESHI Averaging time threshold for I hr averages -I -I

THRESH3 Averaging time threshold for 3 hr averages -I -I

THRESH24 Averaging time threshold for 24 hr averages -1 -I

THRESHN Averaging time threshold for NAVOhr averages -I -I

NDAY Accumulation period, days 0 0 NCOUNT Number of exceedances allowed I I LECHO Echo option F F LTIME Time series option F F LPLT Plot file option F F LORD Use grid format instead ofDATA format F F

LDEBUO Output information for debugging? F F

5-7 

### Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 063701.0044**.~.** ".

### APPENDIX 

### A - EMISSION DATA

Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044

# ••••••••••••••••••••••••••••••••••••••••••••

Printed /013012006

EUG Emission Point (Unit ID) Point ID Unit Description (Include MMBtu/hr if applicable) Which ofthe 26 Categories does this Unit

fall into?A Year Installed Was this Unit in Existence on August 7, 1977? (Y or N) Was this Unit in Operation after August 7, 1962? (Y or N) Potential to Emit (tons/year)
PM 10 NO. S02 VOC NH3 2A 221 EP04 #2 discharJ!e chute, seperator, bucket elev. 3 1964 Y Y 28.90 2.70 0.02 0.20 28 222 223 EP05 EP05 Raw mill Mpump hopper Raw mill #2 weiJ!h meters andfeed belts 3 3 1964 1964 Y Y Y Y 3.20 4C 404 EP\3 Kiln #2feedsystem 3 1964 Y Y 4.20 58 502 EPI5 Kiln, screw conveyors, & bucket elevator 3 1964 Y Y 128.40 1261.00 2102.00 5.90 78 706 705 EP24 EP26 #2 bucket elevator andseparator Finish mill #2 3 3 1964 1964 Y Y Y Y 25.50 3.20 7D 704 EP25 Finish mill #2 feed belt 3 1964 Y Y 3.20 6A 602 EPI7 Clinker cooler #2 3 196411997 Y Y 34.50 II 1101 1102 1103 EP27 EP28 EP28 Coalkoak unloading Coalkoak drag conveyor, bucket elevators Coallcoak crusher 3 3 3 1974 1974 1974 Y Y Y Y Y Y 7.00 2.40 ; Tolal 240.50 1263.70 2\02.02 6.10 0 A. See attached list ofcategones

Lafarge Building Materials Using Title V Trinity Consultants BART Modeling A-1 063701.0044

# ••••••••••••••••••••••••••••••••••••••••••••

Printed 1013012006

Source Emissions Emission Factors Location Point III Unit Description Constr. Datel Heat Duty Process Rate Hours ofOperation PM 10 (MMBtu/hr) TPH (hrs/day) (2c/acf) Tulsa EP05 Rawmillpump hopper 1964 -­ 100 24 0.025 Tulsa EPI3 Kiln #2feedsystem 1964 -­ 120 24 0.025 Tulsa EP24 #2 bucket elevator andseparator 1964 

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

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

## -­ 39 24 0.025 Tulsa EPI7 Clinker cooler 112 196411997 -­ 96 24 0.025 Tulsa EP27 Coalkoak unloading 1974 -­ 150 24 0.025 Tulsa EP28 Coalkoak drag conveyor, bucket elevators, crusher 1974 

## -­ 150 24 0.025

Construction dates and source parameters as listed in Permit 98-174-TV. Construction ofEP 14 (olher kiln) commenced on 1961 which is outside the 7-Aug-62 to 7-Aug-77 BART eligibility window. Source Emissions Emission Data Input Location Point III Unit Description NO: PMlO PM (coarsei PM (finel SOl Ibslhr Ibs/hr Ibslhr Ibs/hr Ibslhr Tulsa EP05 Raw millpump hopper 0.429 0.300 0.129 

## -­ Tulsa EP13 Kiln #2 feed system 0.557 0.390 0.167 

## -­ Tulsa EP24 #2 oucket elevatorandseparator 3.402 2.381 1.021 

## -­ Tulsa EP26 Finivh mill #2 

## -­ 0.429 0.300 0.129 Tulsa EP25 Finish mill #2 feed belt 0.429 0.300 0.129 Tulsa EPI7 Clinker cooler #2 

## -­ 4.607 3.225 1.382 Tulsa EP27 Coalkoak unloading 

## -­ 0.375 0.263 0.113 

# -­ Tulsa EP28 Coallcoak drag conveyor, bucket elevators, crusher 

## -­ 0.321 0.225 0.096 ! 2 PM Speciation based on AP 42 Table 8.2-2 category 4 (Dryers, grinders, etc.), PM2.5 is 30% ofTSP PM IO­2.5 (filterable) 70% PM (coarse) PM2.5 (condenSible) 30% PM (fine)

Lafarge Building Materials Fabric Filters Consultants BART Modeling A-2 063701.0044" Printed 1013012006

Emission Factors for NO. &**S~** based on AP-42 (1/95) Table 11.6-8, PM 10 based on baghouse manufacturer exhaust loading ratc. 2 Construction dates and source parameters as listed in Pennit 98-1 74-TV. Unit EPI4 constructed outside ofthe 7-Aug-62 to 7-Aug-77 BART eligibility window 3 Where necessary, Ib/MMSCF were converted to IblMMBtu using a natura! gas heating value (i 1020 BTU/SCF

• Condensible/filterable PM speciation were based on the following profIles: National Park Service guidance for PC Dry Bottom FF. http://www2.nature.nps.gov/air/Pennits/ect/ectCemcntKiln.cfin Coarse 6.70% Fine Soil 7.70% Fine soil included in model as PM (fine) Fine EC 0.00% CPM lOR 75.30% Inorganic condensible PM considered sulfates pcr NPS guidance ".. CPM OR 10.30% Organic condensible PM considered seconpary organic aerosols per NPS 

# ' ••**~.~.:::I**

# ~ .,**~:** Lafarge Building Materials Kiln Trinity Consultants BART Modeling A-3 063701.0044**~. ~** 

### ~'**l~** Source Emissions' Maximuin Hourh Emission Factors Location Point ID Unit Description Constr.Datel Process Rate (ton cllnker/brl Exbaust Flowrate (acfm) Max Hours ofOpel1ltiol (hours/dav) NOx (Ib/ton d.) PM 10 (gr/acfm) S01 (lb/ton d.) Tulsa EP15 Kiln, sc:ew conwyors, & bucket elevator 1964 96.0 160,000 24 6 0.025 10

Source Emissions Location Point ID Unit Description S01 SO. NO. Emissions' PM 10 (Ib/hr) PM (coarse)' PM (fine)' EC SOA Tulsa EPI5 Kiln, screw conveyors, & bucket elevator 960.00 12.91 576.00 17.143 1.149 1.320 0.000 1.766• >;;.':

Printed 1013012006

SourceEmissions Emission Factorl,l Location Point ID Unit Description 

### I Constr. Datel Heat Duty SOl MMBtu/hr Ib/MMBtu discharge chute, separator, bucket elev. I 1964 1.2 0.001 NO. & S02 emission factors based on AP-42 (7/98) Table 1.4-1 and 1.4-2,PMIO based on baghouse manufacturer exhaust loading rate. 2 lb/MMSCF converted to IblMMBtu usinga natural gas heating value of 1020 Btu/SCF SourceEmissiolls Emission Data Inpue Location Point ID Unit Description SOl S04 NO. PM 10 PMco...... PMfine EC SOA (lb/hr) Tulsa EP04 #2 discharge chute, separator, bucket elev. 0.001 -­ 0.0&9 3.&57 0.964 2.&93 -­ -­ 3 Condensible/filterable PM speciation were based on the following profiles: For Natural Gas, National Park Service guidance for NG-fircd combustion turoincs. http://www2.nature.nps.gov/air/Permitslect/cctGasFircdCT.cfm PM IO•2.5 (filterable) 25% PM (coarse) PM2.5 (condensible) 75% PM (fine)

Lafarge Building Materials Combustion Sources Trinity Consultants BART Modeling A-4 063701.0044



# .. .'" •.••,.;••,;-. "'-":il!!IIII 

# ••••**ti.:<.~.**

# ••••••••**\~** ", ,""" 

## ",**''''~'**
> 'r,'. ';"" ' .'

# :~

# ""r.'~" 

# ..**~'.'**

# .'

:

Printed 1013012006

Stackparameters (English) 1 Location Point ID Stack Flow Exit Velocity Stack Height Stack Dia. Stack Temp. Base Elev. 1 ftlsec feet feet Fahrenheit (feet) Tulsa EP04 36,000 11.94 100 4 180 875 Tulsa EP05 4,000 8.29 55 1.6 70 820 Tulsa EPI3 5,200 9.55 97 1.7 200 851 Tulsa EPI5 160,000 17.83 60 6.9 500 768 Tulsa EP24 31,750 10.53 100 4 200 890 Tulsa EP26 4,000 33.16 64 0.8 70 854 Tulsa EP25 4,000 33.16 55 0.8 200 845 Tulsa EPI7 43,000 3.56 53 8.0 200 788 Tulsa EP27 3,500 4.64 7 2.0 70 754 Tulsa EP28 3,000 4.91 135 1.8 70 805 Slackparameters (Metric) Location PointlD Exit Velocity Stack Height Stack Dia. Stack Temp Base Elev. (m/s) (m) (m) Kelvin (meters) Tulsa EP04 3.64 30.49 1.22 355 267 Tulsa EP05 2.53 16.77 0.49 294 250 Tulsa EPI3 2.91 29.57 0.52 366 259 Tulsa EPI5 5.44 18.29 2.10 533 234 Tulsa EP24 3.21 30.49 1.22 366 271 Tulsa EP26 10.1l 19.51 0.24 294 260 Tulsa EP25 10.11 16.77 0.24 366 258 Tulsa EPI7 1.09 16.16 2.44 366 240 Tulsa EP27 1.42 2.13 0.61 294 230 Tulsa EP28 1.50 41.16 0.55 294 245

Source parameters obtained from Pennit 98-1 74-lV 2 Base elevations obtained from visual inspection ofUSGS topographic map!

Lafarge Building Materials Stack Parameters Trinity Consultants BART Modeling A-5 063701.0044

### APPENDIX 

### B - DATA FILE LIST

LIST OF SURFACE STATIONS 

# • 

# • • • • 

# • • • • 

# • • • • 

# • • • • • • • 

# • • 

# • • • 

# • • • 

# • • 

# • • • • • 

# • • • 

# • • • 

# • 

# • • Station ID Location Station No. official ID) Longitude
Winfield/Ark. City. KS

___ 7:..:2:..:.45"'°"'2'--_
..........._.I!.,!.§§
......... .......**~.9.?.QJ.~** ....... Bartlesville. OK

___ 7:..:2::.35"'6"'5

__........_.J.§.,1§.?
..........**_._~017..** __ -"C"'h""'ildr=es"'s.:...TX!..!.>..

____________ 7'-'2"'3""60""4'--_
.............**J.~,.~.JJ..** ...............__:lQQ.281
.._ Coffevville MunicinaL KS 724519

...... ChanuteiMartin John. KS 13986

Clinton-Shennan. OK 723526 E'!QQ J.1.§'§'?'
..

........3..?• .J.J.3.. ....3??.97.
.. _

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

..........:'!.9•.WQ
.... Denton. TX 722589
....3.J.,.7..QQ.....
..__**~.** -97.183 -"E:o..lD=orad=0"-/G"'oo=d:..::Wl'""·n"'"
....**A~R'__** ______-<-93"-'9"'9=-2

__...............3..1,.fJ7...........•..

__ :92.817 -'V..::an""c""eC!. AF =.:oB"' lE oe!D1""·d."'-"O""K"-

_____. ___ 7:.:.23"'·""53"'5'--_..•..........••J§.J.n.•...............**~_-'!.1.9J)_~** ShennanlDennison. TX

__--'7:..:2""26""'7c.::3

__...............**_.I~.,.1J7......................** ..:9§,§§7
...

.........

F..,re"'d"'en""·c",k.:...:O""K"'-

__________""22"'2:::2""2

__......._...**_J.~,JJ** .. ........._ ....__**~~~.983**

....**F-"-lin"'n"'in."'-"AR~** _______, ____=_55""'5"'5=_5

___•...........}§..JQ9.
.........._--.:9f,467 Fort Smith. AR 13964
...**1?.J}.~** ..................._.........**~4.}.§7** ...._ FavettevillelDrake. AR 93993}§.,QQ..Q
........._....**_..~~4,)§I._** ,""G"'a""l!e"-.""'O""K'__

___________**~13"'9'_'_7"_5** __,_...._ .........3.§,.J9.9..........•.•.....

___**-~9.767_** Hobart OK 93986 J.?,..Q99.....•.......... _----"-9'-"9.""05"'0'---_-1 Hot Snrinl!s. AR
. __--'7C-"2=.::34"'1-=-5

__._............**J.~d§.7.** ....•........•.._.__**~.91Q~~_** ..**~H"'arn'-'-"'·so"'nIB"-"!:o"'o"'ne:.e-"'C""o.::...** '-'A"'R'--

______--"13""9:..<.7
.... 1

__ ...............J.§.,.m
......................__..•.

## n..LE
............**.!!M...,c",C",onn"""e""Il..!.AFB",,-,,~.** K=S

________--"3""92"'3'---

__.................3..1:.!i.!.?.•_ ........_ .•....**,_~.9.?,.~7..,"_. -'-W'-'i"'chi-"-·ta.""'-'K~S~**  __________-"3'-"92"'8"--

__...... _ ..........E.§..?Q
............. _ •.**_:91A3~._** Jonlin Re!!ionaL MO 723495 Lawton. OK 99999 LiberaL KS 11111 .3..1,)'?'Q

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

.....J.7Jl.11. ....:94.5iJQ._

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

.... :100,959 ""L""itt"'le"-'R"'0"'c",ki""A"dam""""-'s
...!:,AR=-

_______--"13""9""6""3

___..........•J.4,.n
..3.
...._..........___**-9.b.~d.L........ -"L"'itto!!le"-R-"'0"'c"'k'-"AFB~~.-"'AR'-"'-**  ________**7.!..:.2"'3~40~5'--_** ..............__**~.4,.9J7.** .._ •...........**_-..::.-?_'.92"'_.!,,15~0:......._1** Muskol!eelDavis Field. OK
. ____..!.1-"1l'-!1""2

__.._......_.J.
.**.?,§..~.Q**.............**~_-95.3.§7~** McAlester Re!!ional. OK 93950

...............H,2,Q9.
......... .........**:9.~** ...

## m....... Mena International Mun.• AR

___ 7,-,2",34",3<=5,--_
.........._ .....**3..~** •.**?_~9.** ..........._ ...__ -9.1:217

..**-"'O"'k1""ab-"-0"'ma~C"'i""tv~.-"'O"'K'--** _________"13"'9""6L.7

___ ..........J..?.,1QQ
........_..__..:27..60Q
.._**.!-'N"'onnan'-"-""'-'.:-W"-e""s""th""e"!!irn"'e"'r.'-'O~K"'-** _______**7~2"'35""'7'-"0:......._** ......_ .........,}J.,fJJ.
......._ .... _ ......:'t.7..42..Q._
.**.::.Mo=e=dl"'·c~ino=.e-=L""od=.:e:"'e"-'.**K""S"--

________ 7:..:2:..!.45""2""0'--_
......_ .....

### J_7.,.m_
......_._.................**;~~.?~Q** ....**.!..P"'on""c"'-a~C:!!itvc:.:.'-'O"-'K~** _________---"13""'9'-"6"-9

__....._ .........**J.~** •.

## m..................**~,lQ.O** __**.!..P",ars~o""n""sIT-'!.!ri,""C""itv'-!"'-" KS~**  _________**7'-"4c.::54"'3~1**___..........._.3..1,.m
............_-...:91"lli........... PariS/Cox Field. TX

__---"'33"-=3"'3"'-3

__..............**J.~..§.~.~** ...__........._._...._:95.42Q

__.

...Rol!ers.AR

__-'6"'6"'66"'6

___.........._J.§.,},§..?
...._........_.:94.J..Q()
.. Snrinl!field. MO 13995

_ ..........._.E.7.JJ.
...._ .........._ -93.3.8L_. 

### Lafarge Building Materials Refined BART Modeling Protocol Trinity Consultants 06370l.0044

# •

• • • • 

# • • • 

# • • 

# • 

# • • • • • • 

# • • 

# • • • 

# • • • 

# • • • • 

# • • • -- ------------------_.......

# _ ...................

## __........- ....-----.---..­**-"S"il~o",am~S!!!on",,·n~e:,,-,. AR~**  _________**....!7.!...77~7'-'7**__ .._ ...J.(!,m_
..._ ....._. ......**:9.1.A~.~** -,-W,-,i",ch""ita",....F.."al""ls,,-.**TX~** ________ . __J}2§..§...

___.._.........**;3.}.,.9~** ........................:9.l!,?J.7.
..... -"S"ti",ll",w""ate",r"-.."'O"'K'--

_________---'I..!.I'-!1l"'3'___
......_.......3.JiJJ9_
........__........**.:91J!~~** ..**.!..F"'ort"""Le"'o"'n"'ar""d....!W-'-'o"'od.""-'M~O'__** _______'7c:2""44"'5'-'-7

__._.....**_.3.lJ.~.Q.** __...._.........:-""92=-
....,15"'0'--_1 -,T

.... ink=eC!...r**AFB=.;"".'-'O"-'K~** _________-'-'72=3""54..:.;0'--_._
.....**_J~,.~n?._** 

### ......_-----.::2.:UllL_.**-"M:..:c::>.Kinn,.·=ev.:..:.-'TX~** __________ 7,..22""'-54_'_'1'___
.............._ •.**~}J83** .. ...:.96.583

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

.... T.!::ex"'ar"'kan"""a"..R.,e""I!1"'· o"'n"'a1.!:....!.TX"-'-

_______..!.\3"'9'-'7.!...7

__._......**J},1.~.Q** .. _ .......**~,000** __ Bentonville. AR

_____--"'88"'8,.,,8""8

__...........J.(!,J.?Q
...........__.... ..:91.,71.7.
..

........EnidIWoodrine:. OK 44444**.3..(!,J.~3.** . ..**~n7~}** Ardmore Municinal. OK 723555 34.300 -97.017

LIST OF PRECIPITATION STATIONS Cooperative Latitude Longitude station ID Location/name State (deg.) (deg.) 30130 ALUM FORK AR 34.480 ·92.510 30178 ANTOINE AR 34.020 -93.250

....

### _._._...

### _._ ..­**---'~'~** 30220 ARKADELPHIA 2 N AR 34.090 ·93.030 30458 BATESVILLE LIVESTOCK AR 35.500 -91.460 30530 BEEBE AR 35.040 -91.540 30616 BERRYVILLE 5 NW AR 36.260 -93.370 30764 

### '._._-_

... BLAKELY MOUNTAIN DAM

..._ 34.340

.....__..H.

___._. AR

..H.

__..._.H
.._ -93.120

__ 30798

••••__ __ ...

# _ .... BLUE MOUNTAIN DAM AR

.......

# _-_ 35;070

...

# _

..... ... ......

## - ........

# __.....

# __. -93.390

••• H. 

# • 

# • ...­

.... 30832



# __._. BOONEVILLE 3 SSE AR 35.060 -93.550 30842 BOTKINBURG 3 NE AR 35.430 -92.280 30900 

## _ .. BRIGGSVILLE AR 34.560 -93.300 -92.340 31020 BULL SHOALS DAM AR 36.220 ...•

# _ ..........

# _ ..........

## _-_..

### ­ --_..."._
......

# __.....

### _.,----_....

## _-_. 31140 CALIONL&D AR 33.180 -92.290 1--._._-
... 31152 CAMDEN I AR _._-,-_.-**33.3~._ -92;~**  31457 CLARKSVILLE 6 NE AR 35.320 -93.240 31582 COMPTON AR 36.050 -93.180 31952 DE QUEEN DAM AR 34.060 ·94.220 

## '._--_

.. 

## ._ ..

### _-_.._._-_._-----_._.

_ ....

# __._..... 32020 DIERKS DAM AR

.....

# _ 34.090

...

# __...

### _ ...

# _ -94.050

...

# __. __._-_. DUMAS AR 33.530 -91.290

.. 

# • 

# • • 32300 EL DORADO GOODWIN FD AR 33.130 -92.490 32356 EUREKA SPRINGS 3 WNW AR 36.250 -93.470 32444 FAYETTEVILLE EXP STN AR**36.0~** __. -94
..!.D0
__ ..- ......

## ---.... 32489 FERNDALE 6E AR 34.460 -92.270

.._.

# __..

# _ ..... --'-"'-'---'-''''--'­**~~-.** FOREMAN AR 33.430**._.-94:.~** 32574 FORT SMITH MU, OK AR 35.200 ·94.220 

# • 

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

# • 

# •'

# •

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

# • • •• 

# r.

• •

h 

## '.:.'.'.'.**~~.' ~~l'l'**  '
.. ..**~** ,:' 

# I:

:,:' 

## .........•**t~. ~;~.**  

# i,·**~.".** 'l

# •

• • • 

# • • 

# • 

# • 

# -. • 

# • • • • 32794 GILBERT AR 35.590. ·92,430. 32810 GILLHAM AR 34.120. -94.150. AR 35310
...._.__...._--<_..**:2.~:goO** ...._ AR 36.170. -91.280
.._ 33165 HARRISON FAA AP AR 36.160. -93.0.90.

.....­**~5_44** __..___. ..::HUN=-'-T:..:S=-V:...:I=LL=E=-:cl-=S.::.S-'-W AR 36.0.40.

____**~=-=** _______--'.:=-

__::.::.:.=-:..;: -93,450.
.._ 34185 LEWISVILLE AR 33.220. -93340. 34248
..._ .....__....... N_O_R_TH LITTLE ROCKIMUNIC:c;:I;;.:PAc::L:::.._--:-=c::_ 34.50.0. -92.160.

__ AR

...

# _ .........

### _._.__•......


### _ ...._. 34756 AR 34.340. -94.160.

.....

# _ ..........

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

_______.:c='---­**-~------** ..

### -.----. 34988 MOUNT IDA 3 SE AR 34.320. -93.350. AR 34.0.90. -93,430.

..

# _._-_
.......

### _._....

## _-
........--......... ....**~ ~**  AR 

## -_.....

# _._ 33.560.

......

# __._..

# __. 

### :2}.,?!.Q
..._ AR 

### ...**}..±:~?L** 

### .._:2.3..:.1
..

### gg __ AR 36.150. -92.150.**35~_28_** ...__. __**._..::N~O~R~F~O~RK~D~A~M~** ______**._~** __**~~~** __**~~_-1** 35320. NORTH LITTLE ROCKIMUNICIPAL AR 34.50.0. -92.160. 3560.2 PARTHENON AR 35.570. -93.150. 35754 PINE BLUFF AR 34.140. -92.0.10. 3590.8



## _ ....

# __....,._-_•....**~~_39_3._** •.

___-=S...oT..::P.:..A:.:U:.::L=---

__________-:..;A::.;R'--
_ _'3::..:5
...490_.,
.........:23.46.L**~.** _____**"S~T:..:U:..:TT:..::..:G~A~R~T=-9=-=E~SE~** ______**..::A~R~_~3~4.=28~0'~** __**-.:..9~1.~25~0.~ 37488~**  __. ___**W~A~L=D~R:c;:O~N~** ___________**~A~R~_~3~4.~54~0'~.** __**.~94~.~0..6~0.,--~** 141427 CHANUTE M JOHNSON AP KS 37,40.0. -95.290. 141699 COLBY 1 SW KS 39.240. -10.1.0.40. 

### .--..

## _--_._--_.


# __.._--_._-----"':::;.._

...........

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

___**w.n'~_N_** 142432 ELKHART KS 37.0.0.0. ·10.1.530. 142980. GARDEN CITY EXP STN KS 38.0.0.0. -10.0..490.

...

# _-_._-_.



# __...

# __. 

# __._-_._--_

....................- ........... 143954

...__**IN_D_E_P~_N_D_E_N_CE** ___.______--'K:.:::S:::,-_
........ }1:1iQ._'
.......**_:2.S.:i~9** 

### .._ ..**-~** ....-- ..- ..­

_ou ... _ •••

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

### ... !.451**~5_.....** MEDICINE LODGE KS 37.170. -98.330. 147922 SUBLETTE KS 37.290. -10.0..50.0. 148830. WICHITA KS. KS 37390. ·97.260. 230.789 BOLIVAR**~ll?..** 

### ___JOPLIN MUNICIPAL AP MO 37.0.90. -94.30.0.**~_45_44** _______...**_..::K=I~RK~S~V=IL~L~E~** _______________**~M~O~** __**~~~~** ____**~~~** -92340.
..._ 40..120. 234825 LEBANON 2 W MO 37.410. -92.420. 234919 ._-=L:.:;IC::::K:.:;IN:.:..:.;::G:...4:..;N"'--

_______---'M=O 37.330.**-9!~** __**.::.:.:=~** 

# ___:;.,:**~.3.**

### ..**5~?** ___**._._~O_C_KW_O_O_D** ____. __._..._____'M=O_ 37.230. -93.570.

••••••• n

__ ••••••• _ •••• _ ••••• _ ••• _. 23530.7 MARSHFIELD MO 37.20.0. -92.550.

..

## _--_.**~355~~** ___ .__**-=-M-=I=LL=E=R~I_=E'--** ____________'M:..=O 37.130. -93.50.0.

__**_=..:.~:........** __-'_=:.:::..::.:........_1 236460. OZARK BEACH MO 36.40.0.**-93~~** 

### Lafarge Building Materials Trinity Consultants Refined BART Modeling Protocol 

### 063701.0044

# •

• • 

# • • • • 

# • 

# • ••


# • 

# • • 

# •

• 

# • • 

# • • • 

# • • • •


# • • 

# • • • 

# • 

# • • 237435 ST JOSEPH ROSECRANS MO 3M70 -94.550 MO 37380 -91.320
MO 36.590 -94.320 MO 37.140 -93.230 MO 37.420 -93.460 MO 36.360 -93.180

OK 36.450 -96.000 I::c.:.=.=-----.='-'.::::'-=----------Q!S----.**34:.?~.** -95.050

LAKE OVERHOLSER LEHIGH MARSHALL MCALESTER MUNI AP 346130 MUSKOGEE**~:4.648:5.** __.. NOWATA 346616
_ ..._ OIL**T0!i2.~E ~~~~"'M'~_M**  346638 OKEMAH**""M_~"** __ ......_ .._. 346661
. __ » OKLAHOMA CITY, OK 346944 PAWNEE5N 347201 PONCA CITY MUNI AP**347~** ___.. RANGE RIVERSIDE 4 W
348501 STILLWATER 2 W 348992 TULSNINTL OK. 349404 WAYNOKA 349629 WICHITA MTN WL REF M""""_""_'"'' OK

OK 34.290 -97.580
OK 36.180 ...••..." ......

# _...... ,,,

OK 34350 ·99.020 OK 35.290 -97.400 OK 34.280 -96.130

OK 36.040 -96.340

......... _ ••• _ •••••• "_•••••**m~** ••• '
__ ""'-.-
...**-~.~--."".--** ...­ . OK 

### ...**.2.?:.~Q** 

### ____:2
..**(i·18~_~** OK 35.230 -97.360 OK 36.240 -96.490 OK 36.440 -97.060 OK 36.330 -101.050

...**--.-~-** OK 36.470 -100.250

-..**-.~--.--~.-"** .. -99.530 OK 36.070 -97.060 OK 36.120 -95.530 OK 36.320 -98.530 OK 34.440 -98.430**,,,~~_~<ff"** .. M.

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

TX 34.260 -100.170 

# • • • • • • • 

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

# •

• •• CLARKSVILLE 1 W TX 33370 COMMERCE
DALLAS-FORT WORTWFORT. TX. TX 32.540 -97.010 DALLAS LOVE FIELD TX 32.510 -96.510 '-'-' DENTON2SE TX 33.120 -97.060 TX 32.490 -97.210 TX TX 33.090 -96.490

.. "'-,..**,,-""~.,,-,** ...**~-.** TX 32.570 -97.030 TX 34.490 -102.240

TX 33.390 -101.490 MIDLOTHIAN 2 TX 32.290 -97.000 WELLS MUN AP TX 32.470 -98.040

____**.~** ...MQUNT PLEASANT TX 33.100 -95.000

..,.._._.m..._._...**~** NACOGDOCHES TX 31.370 -94.390 ._._--_.....

## _­ TX 35.340
TX 32.150 -98.120 TX 33.090 -95.380**:~** ., TX 33.590 -98.300



# • 

# • • Lafarge Building Materials Trinity Consultants Refmed BART Modeling Protocol 063701.0044 

# • .,' • • •**.~~ '~hr.."**  ,. 

### #"".--------- 

# -;

# ...........
"

# .. ••iifi,.·•••• •**~,** -j "it

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

37
E,' 361 35 34 33 Dodge City-E Pratt-W Pratt-E Wiehila-W Wichita-E Joplin-W Joplin-E Springfield-W Springfield-E Rolla-W Perryton-E Woodward-W Woodward-E Enid-W Enid-E Tulsa-W Tulsa-E Hanisan-W Harrison-E Poplar Bluff-W
PI.inview-E Clinton-W Clinton-E Oklahoma City-W Oklahoma City-E Fort Smith-W Fort Smith-E Russellville-W Russellville-E Memphis-W
Amarillo-E L.wlon-W Lawton-E Ardmore-W Ardmore-E Me Alesler-W Me Alesler-E Little Rock-W Little Roek-E Helen.-W
Lubbock-E Wichita Falls-W Wichita Falls-E Shennan-W Shennan-E Texarkana-W Texarkana-E EI Dorado-W EI Dorado-E Greenwood-W ]0] ]00 99 98 97 96 95 94 93 92



# t, 

# t r;

# '**,,,_t_t_~'** 

# ••••

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

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

### Format)

38.5 37.5 36.5 35.5 34.5 Dodge City Pratt Wichita Joplin Springfield Rolla Perryton Woodward Enid Tulsa Harrison Poplar Bluff
Plainview Clinton Oklahoma City Fort Smith Russellville Memphis
Amarillo Lawton Ardmore McAlester Little Rock Helena

Lubbock


# - Wichita Falls Sherman Texarkana El Dorado Greenwood 102 100 98 

### 96 94 92

:1**~** t: r>: