### RECElVEU

-- 

## &;'EL 

## AIR PLANNING SEC.



### 08SEP23 

## PW l: 26

STATE OF ARKANSAS 

### MIKE 

### BEEHI: (;OVFRNOR



### July 29,2008

Richard Greene, Regional Administrator United States Environmental Protection Agency- Region 6 1445 Ross Avenue, Suite 1200 Mail Code 6RA Dallas, TX 75202-2733

Dear Mr. Greene:

The State of Arkansas is hereby submitting a State Implementation Plan (SIR$ to? address requirements for regional haze visibility protection in mandatory federal Clm I areas. This SIP contains the necessary measures to ensure that the tate of ~rkanss is working toward the national goal for visibility contained in 42 U. 

# C .C. 7491-Visibility protection for federal Class I areas. It fulfills requirements under 40 C.F.R. Part 51, Requirements for Preparation, Adoption, and Submittal of Implementation Plans, Subpart
P- Protection of Visibility.

The Clean Air Act (CAA) and the United States Environmental Protection Agency (EPA) regulations require states to submit SIPs to make "reasonable progress" in reducing visibility impairment in federal Class I areas resulting from anthropogenic pollution. The CAA, 169A(a)(l), "declares as a national goal the prevention of any future, and the remedying of any existing, impairment of visibility in mandatory federal Class I areas, where impairment results from man-made air pollution." Class I areas are national parks, containing more than 6,000 acres, and wilderness areas, containing more than 5,000 acres, which Congress has recognized as significant sites. Arkansas has two such areas: the Upper Buffalo Wilderness Area in the Ozark National Forest and the Caney Creek Wilderness Area in the Ouachita National Forest.

The Regional Haze Rule focuses on older emission sources that have not been regulated under other provisions of the Clean Air Act. Such older sources that could contribute to visibility impairment in Class I areas may be required to install Best Available Retrofit Technology (BART). Regional Haze SIPs must contain such emission limits, schedules of compliance, and other measures as may be necessary to make reasonable progress toward meeting the national goal, including requiring installation, operation, and maintenance of BART, as determined by the State on certain existing stationary sources.

### The EPA Regional Haze Rule strongly encourages states to work together in regional partnerships to reduce haze. As a member of the Central Regional Air Planning Association (CENRAP), Arkansas has worked closely with Kansas, Louisiana, Iowa, Missouri, Nebraska, Minnesota, Oklahoma, and Texas, as well as with federal agencies, tribal areas, and other interested parties to improve air quality in our region and the country.

The State of Arkansas hereby submits to the EPA that it has met the requirements for regional haze visibility protection in mandatory federal Class I areas under 40 C.F.R. 551.308. We are enclosing two hard copies and one electronic copy of the State Implementation Plan for your review.

We look forward to your prompt review of this SIP submission. If you have any questions, please contact Mike Bates, Air Division Chief of the Arkansas Department of Environmental Quality, at 5301 Northshore Drive, North Little Rock, AR 7221 8-53 17,

or by phone at (501) 682-0750. Thank you for your attention to this matter.

Sincerely,

ike Beebe

cc: Mike Bates, Chief

- 

### Air Division, Arkansas Department of Environmental Quality

Enclosure

## Arkansas Regional Haze SIP Appendices - Volume 4 of 6 September 09, 2008



## Arkansas Department of Environmental Quality

Appendices for the

State of Arkansas Regional Haze Rule State Implementation Plan



# Appendix Volume 4 Appendices 9.1A -9.2C

# Appendix 9.1A Detailed Description of BART-eligible Emission UnitsSubject-to-Bart Unit

American Electric Power SWEPCO Flint Creek Power Plant, AFIN: 04-00107

No. 1 Power Boiler (SN-01) The Flint Creek Power Plant produces power using a 6324 million BTU per hour, dry bottom, wall fired boiler (SN-01) to produce sufficient steam to operate the turbine generator at the 558 megawatt gr oss electrical out put capability of the unit. The boiler burns primarily low sulfur we stern coal, but can also com bust fuel oil and tire derived fuels. Fuel oil firing is only allowed during startup and shutdown of the boiler, startup and shutdown of the pulverizer mills, for flame stabilization when the coal is frozen, fuel oil tank maintenance, to prevent boiler tube fa ilure in extreme cold weather, and when the unit is offline for maintenance. Fly ash re sulting from the coal combustion process is collected by two hot side elec trostatic precipitators.

SWEPCO maintains Continuous Emissions Monitoring Systems (CEMS) for NO x , SO 2 , opacity, and carbon dioxide (CO 2 ) on the boiler.

Stack (Emission Point) Parameters

SN-01

Base Elevation (1) , ft 1148 Stack Height, ft 540 Stack Diameter, ft 20.0 Stack Velocity, ft/sec 112 Temperature, F° 275 UTM Northing, km 4013.449 UTM Easting, km 363.123 (1) Above Sea Level

9.1A-1Subject-to-Bart Unit

Arkansas Electric Cooperative Corporation Carl E. Bailey Generating Station, AFIN: 74-00024

Boiler (SN-01) The Carl E. Bailey Generating Station produces power using a 1350 million BTU per hour Riley Stoker boiler (SN-01) to drive a 122 megawatt generator. The primary fuel is natural gas but the facility is also permitted to combust fuel oil. Before BART determination the facility was permitted to us e any grade fuel oil with a sulfur content equal to or below 2.3 per cent. For BART compliance the facility will be limited to combusting fuel oil containing a maximum of one percent sulfur by weight. There is no emissions control equipment connected to the boiler.

AECC maintains Continuous Emissions Monitoring Systems (CEMS) for NO x , SO 2 , opacity, and carbon monoxide (CO) on the boiler.

Stack (Emission Point) Parameters

SN-01

Base Elevation (1) , ft 201 Stack Height, ft 167 Stack Diameter, ft 9.83 Stack Velocity, ft/sec 92.0 Temperature, F° 340 UTM Northing, km 3903.072 UTM Easting, km 648.831 (1) Above Sea Level

9.1A-2Subject-to-Bart Unit

Arkansas Electric Cooperative Corporation John L. McClellan Generating Station, AFIN: 52-00055

Boiler (SN-01) The John L. McClellan Generating Station pr oduces power using a 1436 million BTU per hour Riley Stoker boiler (SN-01) to drive a 134 megawatt generator. The primary fuel is natural gas but the facility is also permitted to combust fuel oil. Before BART determination the facility was permitted to us e any grade fuel oil with a sulfur content equal to or below 2.8 per cent. For BART compliance the facility will be limited to combusting fuel oil containing a maximum of one percent sulfur. There is no emissions control equipment conn ected to the boiler.

AECC maintains Continuous Emissions Monitoring Systems (CEMS) for NO x , SO 2 , opacity, and carbon monoxide (CO) on the boiler.

Stack (Emission Point) Parameters

SN-01

Base Elevation (1) , ft 110 Stack Height, ft 160 Stack Diameter, ft 10.83 Stack Velocity, ft/sec 92.0 Temperature, F° 340 UTM Northing, km 3713.898 UTM Easting, km 519.310 (1) Above Sea Level

9.1A-3Subject-to-Bart Unit

Entergy Arkansas, Inc. Lake Catherine Plant, AFIN: 30-00011

Unit 4 Boiler (SN-03) Unit 4 is a Combustion Engineering tilting tangential fired 5,850 million BTU per hour boiler powering a 552 MW generator. The primar y fuel is natural gas with No.6 fuel oil as the secondary fuel. This boiler came on line in 1970. There is no emissions control equipment connected to the boiler.

Stack (Emission Point) Parameters

SN-01

Base Elevation (1) , ft 328 Stack Height, ft 195 Stack Diameter, ft 17.0 Stack Velocity, ft/sec 10.1 Temperature, F° 254 UTM Northing, km 3810.324 UTM Easting, km 508.781 (1) Above Sea Level

9.1A-4Bart Eligible Unit

Entergy Arkansas, Inc. Ritchie Plant, AFIN: 54-00017

Unit 2 Boiler (SN-02) Unit 2 Boiler is a 6,000 million BTU per hour bo iler powering a 544 MW generator. The primary fuel is natural gas with No.6 fuel oil as the secondary fuel. The boiler is permitted to operate combusting natural gas only or a combination of natural gas and No. 6 fuel oil. This boiler came on line in 1968.

There is no emissions control equi pment connected to the boiler.

Stack (Emission Point) Parameters

SN-02

Base Elevation (1) , ft 180 Stack Height, ft 236 Stack Diameter, ft 12.0 Stack Velocity, ft/sec 93.9 Temperature, F° 242 UTM Northing, km 3816.084 UTM Easting, km 720.501 (1) Above Sea Level

9.1A-5Subject-to-Bart Unit

Entergy Arkansas, Inc. White Bluff Plant, AFIN: 35-00110

Units No.1 and 2 (SN-01 and SN-02) Units Nos. 1 and 2 are identical Combus tion Engineering tilti ng tangential 8700 million BTU per hour coal fired boilers with a ma ximum power rating of 850 MW each. The boilers use sub-bituminous or bituminous coal as the primary fuel and No.2 fuel oil as a start-up fuel. Particulate matter is controll ed by an electrostatic precipitator on each boiler.

Auxiliary Boiler (SN-05) The Auxiliary Boiler is a 183 million BTU per hour boiler burning No. 2 fuel oil as its only fuel type. The original purpose of the Auxiliary Boiler was to provide steam for the start-up of the two primary boilers, SN-01 a nd SN-02. Entergy Arkansas is proposing to use SN-01 and SN-02 as peaking units. The A uxiliary Boiler will be used to keep the two primary boilers hot during the down-time periods to reduce the start-up time to a minimum. There is no emissions control e quipment connected to the Auxiliary Boiler.

Entergy maintains Continuous Emissions Monitoring Systems (CEMS) for NO x , SO 2 , and opacity on both primary boilers. Ther e are no CEMS on the Auxiliary Boiler.

Stack (Emission Point) Parameters

SN-01 SN-02 SN-05

Base Elevation (1) , ft 309 309 309 Stack Height, ft 1000 1000 15 Stack Diameter, ft 25.70 25.70 3.0 Stack Velocity, ft/sec 90.0 90.0 65.0 Temperature, F° 321 321 475 UTM Northing, km 3809.437 3809.447 3809.357 UTM Easting, km 579.101 579.101 579.113 (1) Above Sea Level

9.1A-6Subject-to-Bart Unit

Domtar Industries, AFIN: 41-00002

No. 1 Power Boiler (SN-03) The No.1 Power Boiler was installed in 1968 as part of the original construction of the Ashdown Mill. It has a heat input rating of 580 million BTU per hour and an average steam generating rate of 120,000 pounds of steam per hour (lb/hr) at 850 psig. It combust primarily bark, but it is also permitted to bur n wood chips, wood waste, recycled sanitary products composed of cellulo se and polypropylene, pelletiz ed paper fuel (PPF), tirederived fuel (TDF), municipal yard waste, No. 6 fuel oil, reprocessed fuel oil, used oil generated on site, and natural gas. Natura l gas is only used to supplement other fuels during high steam demand periods.

The No.1 Power Boiler is equipped with a trav eling grate and a combustion air system. To meet applicable Boiler MACT PM em issions standard of 0.07 lb/MMBTU Domtar Industries installed a wet elec trostatic precipitator (WESP) during the spring of 2007.

No. 2 Power Boiler (SN-05) The No. 2 Power Boiler started operation in Fe bruary 1976. It has a heat input rating of 820 million BTU per hour and an average steam generating rate of approximately 600,000 lb/hr. It combust primarily bituminous co al (over 80 percent of the heat input is supplied by coal), but it is also permitted to burn bark, bark and wood chips used to absorb oil spills, wood waste, petroleum c oke, recycled sanitary products based on cellulose and polypropylene, PPF, TDF, municipal waste, No. 6 fuel oil, reprocessed fuel oil, used oil generated on site, natural ga s, and non-condensable gases (NCGs). The NCGs are produced in the pulp and evaporator ar eas. It consist of n itrogen, total reduced sulfur (TRS) compounds, methanol, SO 2 , and minor quantities of other compounds such as methyl ethyl ketone (MEK). Under nor mal conditions, natural gas is not combusted.

The No.2 Power Boiler is equipped with a traveling grate, combustion air system including overfire air, mu lticlones, and two parallel venture scrubbers. The SO 2 loading to the boiler is significant since the boile r burns coal and NCGs. Therefore, the scrubbing fluid includes water and a source of alkali, such as Sodium hydroxide (NaOH) and/or pulp mill extraction stage filtrate.

The No. 2 Power Boiler is equipped with Continuous Emissions Monitoring Systems (CEMS) for NO X , SO 2 , and carbon monoxide (CO).

Stack (Emission Point) Parameters

SN-03 SN-05

Base Elevation (1) , ft 320 320 Stack Height, ft 217 235 Stack Diameter, ft 6.2 12.0 Stack Velocity, ft/sec 87.8 39.1 Temperature, F° 480 125 UTM Northing, km 3722.716 3722.640 UTM Easting, km 396.754 396.759 (1) Above Sea Level

9.1A-7Bart Eligible Unit

Delta Natural Kraft and Mid-America Packaging, LLC Pine Bluff Plant, AFIN: 35-00017

Recovery Boiler (SN-02) The Recovery Boiler is a 100 million BTU pe r hour boiler used to produce steam for various heat requiring processe s throughout the plant. The primary fuel is black liquor, a by-product of the pulping process. Natu ral gas is used as a secondary fuel.

The boiler is equipped with a scrubber to co ntrol emissions of sulfur compounds and an electrostatic precipitator to contro l particulate matter (PM) emissions.

Delta Natural Kraft maintains a Continuous Emissions Monitoring System (CEMS) for Total Reduced Sulfur (TRS) on this recovery boiler.

Stack (Emission Point) Parameters

SN-02

Base Elevation (1) , ft 218 Stack Height, ft 165 Stack Diameter, ft 7.0 Stack Velocity, ft/sec 43.6 Temperature, F° 166 UTM Northing, km 3791.866 UTM Easting, km 589.735 (1) Above Sea Level

9.1A-8Bart Eligible Unit

Evergreen Packaging Inc. Pine Bluff Plant, AFIN: 35-00016

#4 Recovery Boiler (SN-04) The #4 Recovery Boiler is a 1,085 million BTU pe r hour boiler used to produce steam for various heat requiring processes throughout th e plant. The #4 Recovery Boiler burns black liquor solids, a by-product of the pulping process, as its primary fuel. The boiler also has the capacity to burn natural gas a nd No. 6, 5, 4, and 2 fuel oils with a sulfur content not to exceed 1%. Black liquor solids and fuel oil are not permitted to be burned simultaneously. The #4 Recovery Boiler was built in 1968.

The boiler is equipped with an electrostatic pr ecipitator to control particulate matter (PM) emissions.

Evergreen Packaging maintains a Continuous Emissions Monitoring System (CEMS) for Total Reduced Sulfur (TRS) on this recovery boiler.

Stack (Emission Point) Parameters

SN-04

Base Elevation (1) , ft 207 Stack Height, ft 275 Stack Diameter, ft 8.5 Stack Velocity, ft/sec 58.5 Temperature, F° 305 UTM Northing, km 3787.135 UTM Easting, km 600.660 (1) Above Sea Level

9.1A-9Bart Eligible Unit

Georgia-Pacific LLC – Crossett Paper Operations AFIN: 02-00013

6A Boiler (SN-19) The 6A Boiler (SN-19) is a 357 million BTU pe r hour boiler. The boiler burns natural gas and specification grade oil. Specification grade oil consists of new oil, used oil, and pitch from the production of tall oi l. The boiler wa s installed in 1962.

There are no emissions controls associated with the boiler.

9A Boiler (SN-22) The 9A Boiler is a 720 million BTU per hour co mbination fuel boiler that is used to generate steam for general use throughout the facility. The boiler may serve as a backup combustion unit when the incinerator (SN-83) is offline. The boiler was installed in 1973.

The combination of fuels permitted for the boiler are tire derived fuel (TDF), agriculture derived fuel (ADF), refuse derived fu el (RDF), non-condensable gases (NCGs), woodwaste, specification grade oil, natural gas, and sludge.

The 9A Boiler is equipped with a wet venturi scrubber to control sulfur compound emissions. The scrubber was installed in 1980.

Georgia-Pacific maintains a Continuous Em issions Monitoring System (CEMS) for carbon monoxide (CO) on Boiler 9A.

Stack (Emission Point) Parameters SN-19 SN-22

Base Elevation (1) , ft 148 151 Stack Height, ft 140 175 Stack Diameter, ft 13.46 12.0 Stack Velocity, ft/sec 39.7 34.3 Temperature, F° 832 155 UTM Northing, km 3667.358 3667.526 UTM Easting, km 596.129 596.143 (1) Above Sea Level

9.1A-10Bart Eligible Unit

Green Bay Packaging Inc. – Arkansas Kraft Division AFIN: 15-00001

Recovery Boiler (SN-05A) The Recovery Boiler is a 730 million BTU pe r hour boiler used to produce steam for various heat requiring processe s throughout the plant. The primary fuel is black liquor, a by-product of the pulping process. Natural gas is used as a secondary fuel. The boiler was installed in 1975.

The boiler is equipped with an electrostati c precipitator (ESP) to control particulate matter (PM) emissions.

Green Bay Packaging maintains Continuous Emissions Monitoring Systems (CEMS) for Total Reduced Sulfur (TRS) and oxygen (O 2 ) concentration on the recovery boiler.

Stack (Emission Point) Parameters

SN-05A

Base Elevation (1) , ft 323 Stack Height, ft 101 Stack Diameter, ft 5.9 Stack Velocity, ft/sec 84.0 Temperature, F° 362 UTM Northing, km 3883.622 UTM Easting, km 523.874 (2) Above Sea Level

9.1A-11Bart Eligible Unit

Potlatch forest Products Cor poration – Cypress Bend Mill AFIN: 21-00036

Power Boiler (SN-04) The Power Boiler is a 479 million BTU per hour boiler using natural gas as the primary fuel and fuel oil as a secondary fuel. The steam from this unit is used to co-generate electricity in the turbine generator and to provide steam for various mill processes. The boiler was installed in 1977.

There are no emissions controls associated with the boiler.

Potlatch maintains a CEMS to monitor and record the NO X and CO 2 emissions when the boiler is operating with either fuel. A C ontinuous Opacity Monitoring System (COMS) and CEMS for SO 2 are required when the boi ler is operating on fuel oil. The CEMS eliminate the requirement to limit the fuel oil type and sulfur content.

Stack (Emission Point) Parameters
SN-04

Base Elevation (1) , ft 144 Stack Height, ft 2994 Stack Diameter, ft 9.0 Stack Velocity, ft/sec 48.5 Temperature, F° 340 UTM Northing, km 3730.970 UTM Easting, km 663.260 (1) Above Sea Level

9.1A-12Bart Eligible Unit

Lion Oil Co. AFIN: 70-00016

#7 Catalyst Regenerator (SN-809) In the #7 Fluid Catalytic Cracking Unit (FCCU) crude oil is heated to 675°F. The hot crude oil is then contacted with a hot (approximately 1350°F) fluidized catalyst which causes the gas oil to crack into lighter petrol eum products. The catalyst is separated from the products in the reactor and se nt to the regenerator. In th e regenerator, coke which has deposited on the catalyst is burned off and the catalyst is recycled. The hot flue gas is passed through a series of cycl ones to remove catalyst fines and then the hot gas is used to produce steam in the waste heat boiler. Th e hot gases are cooled in the boiler to less than 500°F before passing through the scrubbe r and exiting through the #7 Catalyst Regenerator Stack (SN-809).

The FCCU Catalyst Regenerator Stack was equipped in 2004 with a wet gas scrubber (WGS) for the control of SO 2 , O 2 , and CO emissions.

Lion Oil maintains CEMS to monitor and record the SO 2, O 2 and CO emissions of the Catalyst Regenerator Stack.

Stack (Emission Point) Parameters (SN-809)

Base Elevation (1) , ft 148 Stack Height, ft 200 Stack Diameter, ft 5.75 Stack Velocity, ft/sec 32.0 Temperature, F° 500 UTM Northing, km 3673.684 UTM Easting, km 530.339 (1) Above Sea Level

9.1A-13Bart Eligible Unit

Big River Industries AFIN: 18-00082

Kiln A (SN-01) Big River Industries, formally General Shale Products, LLC, operates a lightweight aggregate facility in West Memphis, Arkans as. Locally mined clay is stockpiled during the summer months and moved to covered storag e areas before processing. Kiln A (SN01) is a coal-fired ro tary kiln that maintains a temper ature of approximately 2,200°F. The clay is transferred to the elevated end of th e kiln by a system of conveyors and it tumbles to the discharge end where it exits as a ggregate. The kiln was installed in 1975.

The discharge gases are pro cessed through a wet scrubber.

Stack (Emission Point) Parameters (SN-01)

Base Elevation (1) , ft 197 Stack Height, ft 100 Stack Diameter, ft 5.0 Stack Velocity, ft/sec 71.8 Temperature, F° 134 UTM Northing, km 3888.519 UTM Easting, km 746.849 (1) Above Sea Level

9.1A-14Bart Eligible Unit

Albermarle Corporation – South plant AFIN: 14-00028

Tail Gas Incinerator (SR-01) Acid gas from the Gas Sweetening Plant is sent to a sulfur recovery plant where 93% of the sulfur is removed from the acid (sour) gas. This sour gas along with the exhaust from the diethylchlorothiophosphate (DECTP) proces s scrubber vent is burned in the Tail Gas Incinerator (SR-01). The tail gas incine rator is designed for a minimum exhaust temperature of 1200°F.

Boiler #1 (BH-01) and Boiler #2 (BH-02) There are two boilers at the Albemarle South facility. Each boiler has the capacity to produce 200,000 pounds of 225 psig steam per hour. Th is is equivalent to a heat input of 340 million BTU per hour. The boilers burn natura l gas, which has been treated either in the sulfinol or the MDEA plants. They may al so burn pipeline quality natural gas. They are not permitted to burn any other fuel.

Emissions generated by the two boilers are permitted under a single bubble.

Stack (Emission Point) Parameters
SR-01 BH-01 BH-02

Base Elevation (1) , ft 285 290 290 Stack Height, ft 190 21 21
Stack Diameter, ft 2.67 8.5 8.5 Stack Velocity, ft/sec 50.0 30.0 30.0
Temperature, F° 1200 450 450
UTM Northing, km 3670.700 3670.758 3670.757
UTM Easting, km 479.964 479.823 479.806 (1) Above Sea Level

9.1A-15Bart Eligible Unit

FutureFuel Chemical Company AFIN: 32-00036

Coal Boilers, 3ea. (6M01-01) The three coal fired boilers are balanced draf t, coal-fired steam generation boilers that have been fitted with atomizing nozzles to fa cilitate burning of liqui d chemical wastes. Each coal fired boiler system is designed as a 70 million BTU per hour unit. The three coal fired boilers share a common primary fuel conve ying system, a common ash handling system, and a common 200 foot tall stack. The boilers are independently controlled by a Distributed Control System (DCS ). All interactions from the operator to the burners are made through th is computer system. The th ree coal fired boilers were installed in 1975.

The spent solvent from the 2,000 gallon liquid wa ste process tank is r outed to either the coal-fired boiler auxiliary waste chemical burne rs or to the burner of the chemical waste destructor (6M03-05). Emissions from tank venting are collected and routed to the coalfired boilers (6M01-01).

Each boiler is equipped with its own electrostatic precipitator (ESP) to control particulate emissions.

Stack (Emission Point) Parameters (6M01-01)

Base Elevation (1) , ft 270 Stack Height, ft 200 Stack Diameter, ft 9.0 Stack Velocity, ft/sec 31.0 Temperature, F° 300 UTM Northing, km 3953.681 UTM Easting, km 633.332 (1) Above Sea Level

9.1A-16Bart Eligible Unit

El Dorado Chemical Company AFIN: 70-00040

West Nitric Acid Plant (SN-08) a nd East Nitric Ac id Plant (SN-09) The East and West Regular Nitric Acid Plants produce weak nitric ac id at concentrations ranging from 52% to 58%. These nitric acid plants employ the DuPont single (high)

pressure process.

Liquid ammonia (NH 3 ) is received from a pipeline and pl aced in pressurized storage at a pressure of 65 psig, or in an atmospheric storage tank. Ambient air is compressed and preheated to approximately 125 psig and 475 o F. A mixture of appr oximately ten percent ammonia and the hot air are reacted over a pl atinum gauze catalyst where the ammonia is oxidized to nitrogen oxide(s) and water vapo r. The nitrogen oxides are then absorbed into water in a cooled absorption process forming nitric acid (HNO 3 ).

Nitrogen oxide emissions from these plants are then passed through Selective Catalytic Reduction (SCR) Units. These SCR Units remove most of the remaining nitrogen oxides by reacting them with ammonia in the presence of a catalyst to form elemental nitrogen and water. The tail gases are then vented to the atmosphere.

The East and West Nitric Acid Plants we re built in 1962. The SCR control equipment was installed in 1996.

Nitric Acid Concentrator (SN-10) The Nitric Acid Concentrator (SN-10) uses concentrated sulfuric acid as a dehydrating agent to remove water from the weak acid nitr ic acid. The use of a dehydrating agent is necessary since nitric acid forms an azeotrope with water at approximately 70% concentration so concentrated nitric acid cannot be produced by using normal distillation. Sulfuric acid at 93% concentration and nitr ic acid at 62% concentration are fed to a distilling tower. Nitric acid at 98%+ concentration distills out the top and sulfuric acid at 70% concentration exits the bottom of the distillation tower.

Emissions from the process are controlled by a hydrogen peroxide scrubber and a venture and packed tower scrubber. These two scrubbers are part of the Nitric Acid Vent Collection System (SN-10).

Stack (Emission Point) Parameters

SN-08 SN-09 SN-10

Base Elevation (1) , ft 207 207 204 Stack Height, ft 75 75 78
Stack Diameter, ft 4.0 4.0 0.5 Stack Velocity, ft/sec 110 105 78.0
Temperature, F° 450 440 103
UTM Northing, km 3680.940 3680.940 3680.914
UTM Easting, km 528.761 528.777 528.838 (1) Above Sea Level

9.1A-17

# Appendix 9.1B Method of Identifying BART-Eligible Sources in the State of ArkansasMethod of Identifying BART-Eligible Sources in the State of Arkansas

The Air Division of ADEQ began the deve lopment of BART-eligible sources by searching the Air Divisions data base for facili ties that fit into one of the 26 BART source categories (Table 9.1B-1). The original list of possible-BART facilities was compiled of those that had the potential to emit 250 or more tons per year (tpy) of the following visibility impairing pollutants, sulfur dioxide (SO 2 )

, nitrogen oxides (NO X ), particulate matter equal to or smaller than ten microns (PM 10 ), volatile organic compounds (VOC) or ammonia (NH 3 ). These facilities were listed in an Excel spreadsheet along with the name and source number of their possible BART -eligible units. The potential to emit emissions of the five above listed pollutant s were also listed. (See Table 9.1B-2)

The units listed as “in operation” before A ugust 7, 1962 were classifi ed as pre-BART and those that were not “in existence” on August 7, 1977 were classified as post-BART. The BART-eligible units were those not “i n operation” on August 17, 1962 but were “in existence” on August 7, 1977.

The term “in operation” means engaged in ac tivity related to the primary design function of the source. ADEQ interprets this as the source is physically performing the task for which it was designed.

The term “in existence” means that the ow ner or operator had obtained all necessary preconstruction approvals or permits. The owner or operator had either (1) begun or caused to begun, a continuous program of physical on-site construction of the facility, or (2) entered into binding agreemen ts or contractual obligations.

In 2004 the Environmental Prot ection Administration (EPA) removed ammonia from the list of visibility-impairing pollutants. The reasoning for the removal of ammonia was the complexity and variability of ammo nia’s roll in the formation of PM 2.5 in the atmosphere. Also, ammonia is the reduction agent used in Selective Catalytic Reduction (SCR) type NO x emissions control equipment that converts NO x to nitrogen and water.

ADEQ did not use VOC emissions from facilities for subject-BART analysis. Only specific VOC compounds from secondary orga nic aerosols affect visibility. These compounds are a fraction of th e total VOCs reported in th e emissions inventory, and ADEQ does not have the breakdown of VOC em issions necessary to model those that only impair visibility. Further, the pr escribed screening model (CALPUFF) cannot simulate formation of particles from anth ropogenic VOCs, nor their visibility impact.

9.1B-1Table 9.1B-1: BART Categories

BART Category ID BART Category BART Subcategory

BART-01 Fossil fuel-fired steam electric plants > 250 MMBTU per hour

BART-02 Coal Cleaning Plants (thermal dryers) Thermal Dryers

BART-03 Kraft Pulp Mills

BART-04 Portland Cement Plants

BART-05 Primary Zinc Smelters

BART-06 Iron and Steel Mill Plants

BART-07 Primary Aluminum Ore Reduction Plants

BART-08 Primary Copper Smelters

BART-09 Municipal Incinerators Capable of Charging > 250 tons of Refuse Per Day

BART-10 Hydrofluoric, Sulfuric, and Nitric Acid Plants

BART-11 Petroleum Refineries

BART-12 Lime Plants

BART-13 Phosphate Rock Processing Plants

BART-14 Coke Oven Batteries

BART-15 Sulfur Recovery Plants

BART-16 Carbon Black Plants Furnace Process

BART-17 Primary Lead Smelters

BART-18 Fuel Conversion Plants

BART-19 Sintering Plants

BART-20 Secondary Metal Production Facilities

BART-21 Chemical Process Plants

BART-22 Fossil fuel-fired boilers > 250 MMBTU per hour

BART-23 Petroleum Storage and Transfer Facilities with a Capacity > 300,000 Barrels

BART-24 Taconite Ore Processing Plants

BART-25 Glass Fiber Processing Plants

BART-26 Charcoal Production Facilities

9.1B-2POTENTIAL BART SOURCES LOCATED IN ARKANSASPOLLUTANTS (TONS PER YEAR FOR CY 2001) PER UNITNH3NOXPM10SO2VOCsDATESINBEGANAIRSEXISTENCE OPERATIONEXPECTED BARTSTATIONARYEMISSIONEMISSION UNITCATEGORY8/7/77? <8/7/62?CLASSIFICATIONSTATIONARY SOURCE NAME/LOCATION SOURCE ID2-DIGIT SICUNIT IDDESCRIPTIONCODEPOTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUALAlbermarle-South Plant / Magnolia05-027-00028 28 SR-01 Tail Gas Incin. 152.6 2.6 0.3 0.3 3,184.0 900.5 0.3 0.3 yes no, 1974 BARTAlbermarle-South Plant / Magnolia05-027-00028 28 BH-01 & 02 Boilers 1 & 2 21417.0 26.2 149.0 1.4 24.5 3.2 8.4 0.5 yes no, 1968 BARTAmerican Elect. Power (SWEPCO) / Gentry 05-007-00107 49 SN-01 Boiler116,379.0 5,682.3 2,340.0 125.8 28,078.0 14,531.2 92.0 69.3 yes no, 1978 BARTArk. Elect. Coop - Bailey Plant / Augusta 05-147-00024 49 SN-01 Boiler, 1350mm 13,069.0 324.2 572.0 202.4 7,118.0 1,560.0 41.2 9.1 yes no, 1966 BARTArk. Elect. Coop - Fitzhugh / Ozark05-047-00012 49 SN-06 Turbine1447.0no no, 2003 post-BARTArk. Elect. Coop - McClellan Plant / Camden 05-103-00055 49 SN-01 Boiler13,806.0 464.4 607.0 224.7 9,219.0 2,369.1 43.8 10.0 yes no, 1971 BARTArkansas Lime Co. / Batesville05-063-00014 32 SN-11Q Lime Kiln12399.2 374.6 31.6 10.8 226.6 50.0no no, 2000 post-BARTAsh Grove Cement / Foreman05-081-00001 32 SN-P1 Kiln #142,400.0 912.7 85.4 76.2 1,960.0 290.2 42.1 25.1 yes yes, 1958 pre-BARTAsh Grove Cement / Foreman05-081-00001 32 SN-P2 Kiln #242,453.0 993.0 85.4 73.9 1,686.0 270.6 42.1 24.8 yes yes, 1960 pre-BARTAsh Grove Cement / Foreman05-081-00001 32 SN-P3 Kiln #344,231.0 1,368.7 118.3 105.1 2,089.0 196.8 57.0 34.1 yes yes, 5/1962 pre-BARTColumbian Chemical / E Dorado05-139-00014 28 SN-23 Boiler & Dryers 16178.0 82.0 5.3 3.0 1,672.0 513.0 23.0 3.0 no no, 1989 post-BARTColumbian Chemical / E Dorado05-139-00014 28 SN-37 Thermal Oxidizer 16608.0 377.0 17.9 6.0 5,690.0 1,151.0 78.0 7.0 no no, 2000 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-01 #3 Power Boiler 31,038.1 523.7 86.5 61.3 271.6 37.3 93.3 5.0 no no, 1991 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-02 #3 Lime Kiln3291.7 72.1 37.7 2.1 58.3 12.0 63.9 6.8 no no, 1991 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-03 #1 Power Boiler 31,085.9 275.3 1,504.1 711.0 214.0 128.5 214.6 4.5 yes no, 1968 BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-05 #2 Power Boiler 32,514.1 1,796.3 359.2 188.9 4,305.5 2,433.5 206.5 61.1 yes no, 1975 BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-06 #2 Rec. Boiler 31,354.3 499.8 369.7 106.0 1,252.7 37.0 204.6 43.4 no no, 1989 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-09 #2 Lime Kiln3300.5 31.8 223.4 81.6 73.2 5.9 74.9 7.2 no no, 1979 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-14 #3 Rec. Boiler 31,182.6 580.6 409.5 77.9 1,861.5 33.1 600.1 14.7 no no, 1989 post-BARTDomtar, Inc. / Ashdown05-081-00002 26 SN-22 Brownstock Vent 3259.1 197.0 no no, 1991 post-BARTEastman Chemical / Batesville05-063-00036 28 6M01-01 3 Coal Boilers 21488.2 488.2 205.3 205.3 6,213.8 6,213.8 2.3 2.3 yes no, 1976 BARTEl Dorado Chemical / El Dorado05-139-00040 28 SN-07 Sulfuric Acid Plan 102,520.0 1,057.0yes yes, 1955 pre-BARTEl Dorado Chemical / El Dorado05-139-00040 28 SN-08 W. Nitric Acid Plan10840.0 507.0yes no, 10/1962 BARTEl Dorado Chemical / El Dorado05-139-00040 28 SN-09 E. Nitric Acid Plan10840.0 372.0yes no, 10/1962 BARTEl Dorado Chemical / El Dorado05-139-00040 28 SN-10 Nitric Acid Consen10314.5 8.5yes no, 1964 BARTEl Dorado Chemical / El Dorado05-139-00040 28 SN-16A #2 Boiler10349.5 212.0 4.0 2.0 0.4 0.3 1.1 0.7 yes yes, 1944 Pre-BARTEntergy - Blytheville05-093-00113 49 SN-01 Turbine 11265.7yes no, 1974 Permit VoidedEntergy - Blytheville05-093-00113 49 SN-02 Turbine 21265.7yes no, 1974 Permit VoidedEntergy - Blytheville05-093-00113 49 SN-03 Turbine 31265.7yes no, 1974 Permit VoidedEntergy - Couch Plant / Stamps05-073-00004 49 SN-01 Unit 11272.8 38.0 2.7 0.7 19.3 0.4 1.1yes yes, 1943 pre-BARTEntergy - Couch Plant / Stamps05-073-00004 49 SN-02 Unit 211,780.0 96.0 17.5 1.7 69.0 0.4 7.3 3.0 yes yes, 1954 pre-BARTEntergy - Independence / Newark05-063-00042 49 SN-01 Coal Boiler #1 126,674.2 7,792.0 284.7 19.0 35,438.6 11,186.0 153.3 117.0 no no, 1978 post-BARTEntergy - Independence / Newark05-063-00042 49 SN-02 Coal Boiler #2 126,674.2 10,139.0 284.7 20.0 35,438.6 11,830.0 153.3 119.0 no no, 1978 post-BARTEntergy - Independence / Newark05-063-00042 49 SN-05 Aux. Boiler1464.0 5.2 2.6 1.0 0.5 16.0 1.2 0.0 no no, 1978 post-BARTEntergy - Independence / Newark05-063-00042 49 SN-06 Feed Conveyor 1706.5 579.0no no, 1978 post-BARTEntergy - Lake Catherine / Jones Mill 05-059-00011 49 SN-01 Unit 1Boiler11,669.0 107.0 9.1 0.9 2.0 0.2 4.3 2.0 yes yes, 1950 pre-BARTEntergy - Lake Catherine / Jones Mill 05-059-00011 49 SN-01 Unit 2 Boiler11,669.0 123.0 9.1 0.9 2.0 0.2 4.3 2.0 yes yes, 1950 pre-BARTEntergy - Lake Catherine / Jones Mill 05-059-00011 49 SN-02 Unit 3 Boiler12,764.0 92.0 15.4 1.7 4.4 0.3 7.1 2.0 yes yes, 1953 pre-BARTEntergy - Lake Catherine / Jones Mill 05-059-00011 49 SN-03 Unit 4 Boiler114,016.0 1,414.0 72.0 22.0 14.4 4.3 37.9 40.0 yes no, 1970 BARTEntergy - Moses Plant / Forrest City 05-123-00010 49 SN-01 Unit 11993.2 36.0 78.2 0.6 478.4 0.1 3.4 1.0 yes yes, 1951 pre-BARTEntergy - Moses Plant / Forrest City 05-123-00010 49 SN-02 Unit 21993.2 27.0 78.2 0.5 478.4 0.1 3.4 1.0 yes yes, 1951 pre-BARTEntergy - Ritchie Plant / Helena05-107-00017 49 SN-01 Unit 1113,140.0 182.0 86.6 2.5 819.6 0.5 20.0 5.0 yes yes, 1961 pre-BARTEntergy - Ritchie Plant / Helena05-107-00017 49 SN-02 Unit 217,884.0 3.0 132.7 0.2 1,173.1 0.1 35.0 0.0 yes no, 1967 BARTEntergy - White Bluff / Redfield05-069-00110 49 SN-01 Unit 1126,674.2 7,452.0 284.7 16.0 45,727.2 16,182.0 153.3 97.0 yes no, 1974 BARTEntergy - White Bluff / Redfield05-069-00110 49 SN-02 Unit 2126,674.2 10,725.0 284.7 19.0 45,727.2 20,516.0 153.3 113.0 yes no, 1974 BARTEntergy - White Bluff / Redfield05-069-00110 49 SN-05 Aux Boiler1140.5 3.2 5.9 1.0 420.1 9.4 1.2 0.0 yes no, 1974 BARTGaylord Container / Pine Bluff05-069-00017 26 SN-02 Recovery Boiler 367.9 55.9 122.0 13.9 92.7 2.5 1,638.9 489.0 yes no, 1976 BARTGeneral Shale Products / W. Memphis 05-035-00082 19 SN-01 Kiln A19294.8 197.3 131.4 78.6 237.4 151.5 10.6 2.3 yes no, 1975 BARTGeorgia-Pacific Paper / Crossett05-003-00013 26 SN-03 10A Boiler32,192.2 684.0 438.4 246.0 92.0 47.9 661.4 344.0 yes no, 1984 post-BARTGeorgia-Pacific Paper / Crossett05-003-00013 26 SN-18 5A Boiler3403.0 92.8 125.6 0.3 2,413.4 0.1 14.2 0.3 yes yes, 1953 pre-BARTGeorgia-Pacific Paper / Crossett05-003-00013 26 SN-19 6A Boiler3722.7 178.0 262.8 3.8 4,345.0 18.0 12.0 0.6 yes yes, 1961 pre-BARTGeorgia-Pacific Paper / Crossett05-003-00013 26 SN-22 9A Boiler31,511.1 987.0 394.2 112.0 2,601.7 1,700.0 217.7 142.0 yes no, 1975 BARTGeorgia-Pacific Paper / Crossett05-003-00013 26 SN-26 8R Boiler3832.2 724.0 407.3 45.0 1,314.0 1,230.0 172.5 161.0 no no, 1981 post-BARTGreat Lakes Chemical / El Dorado05-139-00012 28 SN-301A Boiler, 150mm 21361.4 13.4 3.5 1.8 525.6 0.1 3.5 0.7 yes no, 1966 BARTGreat Lakes Chemical / El Dorado05-139-00012 28 SN-302A Boiler, 113mm 21272.5 13.1 2.7 1.3 25.0 0.2 0.9 0.7 yes no, 1970 BARTGreat Lakes Chemical-West Plant / El Dorado 05-139-00101 28 SN-102&20 #2 Boiler&Flare 2129.0 28.5 2.3 1.8 312.0 0.0 3.2 2.7 no no, 1983 post-BARTGreen Bay Packaging / Morrilton05-029-00001 26 SN-04 #3 Wood Boiler 3521.0 436.4 145.0 145.5 42.0 35.2 98.7 82.7 no no, 1997 post-BARTGreen Bay Packaging / Morrilton05-029-00001 26 SN-05 & 06 Recovery Boiler 3206.0 537.9 910.0 637.9 494.0 1,786.1 61.0 63.2 yes no, 1975 BARTGreen Bay Packaging / Morrilton05-029-00001 26 SN-25A Paper Machine 3339.8 261.8 yes no, 1965 BARTGreen Bay Packaging / Morrilton05-029-00001 26 SN-25B Paper Machine 3270.1 208.0 no no, 1998 post-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-01 Bark Boiler3894.5 191.6 526.0 57.7 219.0 6.3 447.0 24.8 no no, 1989 post-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-02 #2 Rec. Boiler 3458.6 180.0 202.0 122.3 876.0 531.9 153.3 6.7 no no, 1990 post-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-03 #3 Rec. Boiler 3545.0 187.5 240.0 141.4 876.0 516.1 153.3 1.1 no no, 1990 post-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-04 #4 Rec. Boiler 3795.0 217.6 350.8 274.8 2,376.2 1,861.9 389.8 23.6 yes no, 1968 BARTInternational Paper / Pine Bluff05-069-00016 26 SN-13 #1 Power Boiler 31,489.0 619.5 56.9 24.7 746.0 600.7 6.4 2.8 yes yes, 1958 pre-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-15 #2 Power Boiler 31,489.0 695.9 56.9 26.0 746.0 628.0 6.4 3.0 yes yes, 1958 pre-BARTInternational Paper / Pine Bluff05-069-00016 26 SN-21 Black Liquor Ox.T3873.2 486.3 no no, 1993 post-BARTInternational Paper / Pine Bluff 05-069-00016 26 SN-22,23,24Pulping Area 3349.5 251.8 no no, 1986 post-BARTInternational Paper / Pine Bluff 05-069-00016 26 SN-41, 42 Off Mach. Coater 3 27.7 7.0 3.1 0.8 0.2 0.1 361.3 251.8 no no, 1986 post-BARTInternational Paper / Pine Bluff 05-069-00016 26 SN-54 #2 Paper Mach. 3278.5 223.7 no no, 1990 post-BARTLion Oil / El Dorado 05-139-00016 29 SN-809 #7 Calalyst Regen112.3 260.1 178.7 329.4 230.1 1,945.3 691.8 805.1 1.9 yes no, 1973 BARTLongview Gas / Stamps05-073-00005 13 SN-04 Sulfur Rec. Unit 151.7 0.5 0.2 0.1 938.6 613.0 0.1 0.1 no no, 1982 post-BARTPOLLUTANTS (TONS PER YEAR FOR CY 2001) PER UNITNH3NOXPM10SO2VOCsDATESINBEGANAIRSEXISTENCE OPERATIONEXPECTED BARTSTATIONARYEMISSIONEMISSION UNITCATEGORY8/7/77? <8/7/62?CLASSIFICATIONSTATIONARY SOURCE NAME/LOCATION SOURCE ID2-DIGIT SICUNIT IDDESCRIPTIONCODEPOTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUAL POTENTIAL ACTUALNucor Steel / Hickman05-093-00233 33 SN-01 EAF Baghouse 61,061.3 233.8 162.1 9.2 416.1 190.5 183.1 32.6 no no, 1991 post-BARTNucor-Yamato Steel / Blytheville05-093-00202 33 SN-01 EAF Baghouse 6749.0 355.0 108.2 3.9 295.7 102.0 256.2 111.0 no no, 1991 post-BARTNucor-Yamato Steel / Blytheville05-093-00202 33 SN-03 LMF Baghouse 6153.3 16.7 40.2 13.1 394.2 295.0 21.9 15.6 no no, 1991 post-BARTNucor-Yamato Steel / Blytheville05-093-00202 33 SN-04 Beam Reheat 6297.8 121.0 29.4 18.5 1.1 0.8 9.6 6.6 no no, 1993 post-BARTPotlatch Corp. / McGehee05-041-00036 26 SN-02 Rec. Boiler3651.6 446.3 208.5 129.4 164.8 17.5 38.1 6.3 no no, 1992 post-BARTPotlatch Corp. / McGehee05-041-00036 26 SN-04 Power Boiler3419.6 297.0 209.8 6.7 28.7 0.8 3.1 2.3 yes no, 1976 BARTPotlatch Corp. / McGehee05-041-00036 26 SN-10 Brown Stack Was3658.7 50.1 yes no, 1976 BARTPotlatch Corp. / McGehee 05-041-00036 26 SN-15 Digesters 3288.6 58.7 no no, 1991 post-BARTQuanex Corp (MacSteel) / Fort Smith 05-131-00274 33 SN-01 EAF 6 161.0 72.6 100.0 4.5 332.0 94.5 41.0 20.1 no no, 1984 post-BARTTerra Nitrogen / Blytheville 05-093-00012 28 SN-01 #1 Packg. Boiler 21 282.1 13.6 2.6 0.1 1.5 0.2 0.9 0.1 yes no, 1965 Permit VoidedTerra Nitrogen / Blytheville 05-093-00012 28 SN-02 #2 Packg. Boiler 21 375.8 134.8 3.5 0.8 4.4 1.8 1.3 1.0 yes no, 1965 Permit VoidedTerra Nitrogen / Blytheville 05-093-00012 28 SN-03 Reformer Furnace211,806.7 485.2 31.5 2.8 8.6 6.6 6.2 3.6 yes no, 1965 Permit VoidedTerra Nitrogen / Blytheville05-093-00012 28 SN-07 CO2 Vent21624.2 29.360.1 4.7 yes no, 1965 Permit VoidedTerra Nitrogen / Blytheville05-093-00012 28 SN-12A,B,CUrea Granulation 21876.0 103.4166.4 23.0yes no, 1974 Permit VoidedTerra Nitrogen / Blytheville05-093-00012 28 SN-13 Urea Boiler21259.3 47.5 17.5 1.4 13.1 3.1 4.8 1.7 no no, 1995 Permit VoidedTerra Nitrogen / Blytheville05-093-00012 28 SN-14 Urea Production 211,292.1 85.9yes no, 1974 Permit VoidedU. S. Vanadium / Hot Springs05-051-00002 33 SN-12 MVO Product Ven21832.0 574.0 2.8 0.7 17.5 11.2 1.0 0.0yes no, 1967 BARTWilliams Refining / W. Memphis05-035-00120 51 SN-22 Barge Loading 23480.3 26.7 yes yes, 1954 pre-BARTEastman Chemical name changed to FutureFule ChemicalGeneral Shale Products name changed to Big River IndustriesEntergy - Independence, Previously and erroneously listed as BARTInternational Paper/Pine Bluff name changed to Evergreen PackagingTerra Nitrogen out of Business, Permit Voided 11/19/2004Entergy-Blytheville, Permit Voided 12/02/2000

# Appendix 9.2A Modeling Protocol Used to Determine Subject-to-BART SourcesSTATE OF ARKANSAS DEPARTMENT OF ENVIRONMENTAL QUALITY 8001 NATIONAL DRIVE P.O. BOX 8913 LITTLE ROCK, ARKANSAS 72219-8913 PHONE: (501) 682-0744

FAX: (501) 682-0753

Draft BEST AVAILABLE RETROFIT TECHNOLOGY (BART) MODELING PROTOCOL

June 7, 2006

Prepared by: Mary Pettyjohn, Senior Epidemiologist Planning and Air Quality Analysis Branch Air Divisionii 

## Table of Contents

List of Figures..................................................................................................... iii List of Tables ...................................................................................................... iv

I. Introduction...........................................................................................................1 II. Background...........................................................................................................2 III. BART-Eligible Sources ........................................................................................4 IV. CAIR and Arkansas ..............................................................................................8 V. BART Air Quality Modeling Approach ...............................................................8 VI. BART-Eligible Units Physical Parameters.........................................................10 A. Stack parameters ............................................................................................10 B. Emission rates ................................................................................................10 VII. Air Quality Model and Inputs.............................................................................11 A. CALPUFF No-obs and refined screening modeling......................................12 1. Modeling domain......................................................................................12 2. CALPUFF system implementation...........................................................13 3. Meteorological data modeling (CALMET) ..............................................13 4. Dispersion modeling (CALPUFF)............................................................14 5. Post-processing (POSTUTIL/CALPOST)................................................16 6. Measuring visibility impacts.....................................................................17 VIII. Change in Visibility Due to BART Controls......................................................19 IX. References...........................................................................................................20

Appendix A. Proposed Time-line.............................................................................. A-1 Appendix B. Map of Receptors..................................................................................B-1 Appendix C. BART-Eligible Sources’ Stack Parameters, Base Elevation at Ground Level and Stack Coordinates ...............................................C-1 Appendix D. BART-Eligible Emission Rates used for the No-obs Modeling Run.......................................................................... D-1 Appendix E.Chapter 5 of the CENRAP BART Modeling Guidelines ......................E-1 Appendix F.CALMET Input Control Files................................................................F-1 Appendix G.CALPUFF Input Control Files............................................................. G-1 Appendix H.POSTUTIL Input Control Parameters...................................................F-1 Appendix I. CALPOST Input Control Parameters .....................................................I-1 Appendix J. Chapter 3 of the CENRAP BART Modeling Guidelines.......................J-1iii 

## List of Figures

Figure 1 Mandatory Class I federal areas in the U.S.....................................................3

Figure 2 Arkansas’ Class I areas ...................................................................................3

Figure 3 Map showing the 300 km radius buffer zones around five separate receptors (north, south, east, west , and center) located in the following Class I areas: Upper Buffalo, Cane y Creek, Hercules Glade, Mingo, and Sipsey. This map was developed to determine which Class I areas will Bee assessed during the BA RT determination modeling ...............................4

Figure 4 Map indicating the locations of Upper Buffalo, Caney Creek, and the eighteen BART-eligible fac ilities located in Arkansas...................................7

Figure 5 Map indicating the locations of Upper Buffalo, Caney Creek, Hercules Glade, Mingo, Sipsey, and the ei ghteen BART-eligible facilities located in Arkansas.........................................................................................8

Figure 6 Map of the 6 km CENR AP Central CALPUFF domain...............................12iv 

## List of Tables

Table 1 Class I areas and their locations as well as the supervising agencies ADEQ will be assessing..................................................................................4

Table 2 Fossil fuel-fired steam elect ric plants > 250 MMBtu/hr and Kraft pulp mills facilities with BART-eligible emissions units...............................5

Table 3 Petroleum refineries, sinteri ng plants and chemical processing plants facilities with BART-eligible emission units .......................................6

Table 4 CALPUFF modeling components.................................................................13

Table 5 Species modeled in BART screening analyses.............................................15

Table 6 EPA recomme nded monthly average f ( RH ) for the five Class I areas ADEQ are assessing......................................................................................17

Table 7 Average annual natural levels of aerosol components ( μ g/m 3 ) ....................181 

## I. Introduction

On 6 July 2005, the U.S. Environmental Protection Agency (EPA) published final amendments to its 1999 Regional Haze Rule in the Federal Register , including Appendix Y, the final guidance for Best Available Re trofit Technology (BART) determinations (70 FR 39104-39172). The BART rule requires the in stallation of BART on emission sources that fit specific criteria and “may reasonably be anticipated to cause or co ntribute” to visibility impairment in any Class I area. Air quality modeling is the preferred method for establishing which emission sources cause or contribute to visi bility impairment. Arkansas’ BART modeling prot ocol is provided herein.

According to the Regional Haze Regulations a nd Guidelines for Best Available Retrofit Technology (BART) Determination; Final Ru le (40 CFR Part 51, p 39125), each state is required to develop a BART Modeling Protocol that describes the required methodology to assess the levels of controls needed on sources subject to BART. The aforementioned regulation also requires states to work in partnership with all st akeholders including Tribes, EPA, Federal Land Managers (FLMs), Regional Planning Organizations (RPOs) and the various source operators. Although stat es are required to wo rk in concert with the previously mentioned stakeholders, EPA has the ultimate authority to approve or disapprove a state’s State Im plementation Plan (SIP).

The main objective of this protocol is co mpliance with the RHR visibility improvement goals. To accomplish this goal, the Arkans as Department of Environmental Quality (ADEQ) has set forth three functions of this protocol. First, ADEQ will use the protocol to determine which BART-eligible units are subject-to-BART and must perform a BART-analysis. Second, facilities that ADEQ notifies that are subject-to-BART will use this protocol to conduct post-control modeli ng required for their BA RT-analysis. Third, the results from this protocol will be used to conduct cumulative modeling to show the change in visibility impact on Class I ar eas based on ADEQ’s BART determination and the BART emission limits for facilities base d on their BART-analysi s. The subject-toBART and final modeling will be submitted to the EPA as part of the BART section of the Arkansas State Implementation Plan (SIP) for Regional Haze.

The AR RH SIP submittal deadline to EPA as set forth in the Regional Haze Regulations and Guidelines for Best Available Retrof it Technology (BART) De termination; Final Rule (40 CFR Part 51, p 39156) is December 17, 2007. To meet this deadline, ADEQ has developed a schedule for completing BART determinations and implementing the BART strategy in order to meet the mandato ry SIP submittal deadline (Appendix A). As shown in Appendix A, the modeling results must be completed no later than March 1, 2007.

The Central States Regional Planning Asso ciation (CENRAP) cont racted with Alpine Geophysics, LLC to develop a modeling protocol for the states within CENRAP’s region of which the state of Arkansas is a memb er. On December 22, 2005, Alpine Geophysics, LLC delivered the final version of the CENR AP BART Modeling Guidelines (Tesche, et2 al, 2005). However, comments from EPA Re gions VI and VII and the Federal Land Managers (FLMs) were not incorporated in to the guidelines; t hus, Alpine Geophysics, LLC rewrote the guidelines to reflect the comm ents from Regions VI and VII and FLMs. These guidelines were re-issued February 3, 2006. Hence, CENRAP’s BART Modeling Guidelines (Tesche, et al, 2005) have been approved by Regions VI and VII and the FLMs. Therefore, the Planning and Air Qual ity Analysis Branch, Air Division, Arkansas Department of Environmental Quality ha s chosen to adopt the CENRAP BART Modeling Guidelines as ADEQ’s BART Modeli ng Protocol. Additionally, in preparing this draft protocol, ADEQ also consulted th e following draft BART modeling protocols:

1. Best Available Retrofit Technology (BART) Modeling Protocol to Determine Sources Subject to BART in the State of Kansas draft version February 24, 2006 2. Best Available Retrofit Technology (BART) Modeling Protocol to Determine Sources Subject to BART in the State of Minnesota draft version February 24, 2006

This draft protocol is most similar to th e CENRAP BART Modeling Guidelines. These guidelines were developed to ensure “consiste ncy between states in the development of BART modeling protocols and to harmonize the approaches between adjacent RPOs” (Tesche, et al, 2005).

Soon after the finalization of this modeli ng protocol, ADEQ will notify sources subjectto-BART. For those facilities subject-to -BART, ADEQ will provide guidance for conducting their BART-analyses.



## II. Background

The Clean Air Act Amendments (CAAA) of 1977 established 156 Class I areas where visibility was determined to be an important value (Figure 1). Areas designated as Class I areas are those national parks exceeding 6000 acres, wi lderness areas and national memorial parks exceeding 5000 ar eas, and all international park s that were in existence on August 7, 1977. While Rainbow Lake W ilderness Area, Wisconsin has been designated as a Class I area, the FLMs have in dicated that visibility is not a valuable characteristic and therefore, is not included in BART or other RH analyses.

The state of Arkansas has w ithin her boundary two mandatory Class I federal areas (Class I area), Upper Buffalo Wilderness Area a nd Caney Creek Wilderness Area which are managed by the United States Forest Service (F igure 2). However, there are two Class I areas in southern Missouri that are located downwind of facilities operating in Arkansas. The Missouri Class I areas are Hercules-Glade Wilderness Area (US Forest Service) and Mingo National Wildlife Refuge (US Fish and Wildlife). While EPA has not listed the maximum distance from a Class I area to model, this criteria has been set by CENRAP as 300 km. As shown in Figure 3, the eastern po rtion of Arkansas is within the 300 km radius of Sipsey Wilderness Area (US Forest Service), Alabama. Therefore, there are3 five Class I areas Arkansas will be perfor ming BART determination/exemption modeling (Table 1).

Figure 1 Mandatory Class I federal areas in the United States of America

Figure 2 Arkansas’s Class I areas4 Figure 3 Map showing the 300 km radius buffe r zones around five separate receptors (north, south, east, west, and center) locat ed in the following Class I areas: Upper Buffalo, Caney Creek, Hercules Glade, Mingo, and Sipsey. This map was developed to determine which Class I areas will be asse ssed during the BART determination modeling

Table 1 Class I areas and the State they are locat ed in as well as the supervising agencies ADEQ will evaluate during the BART determination/exemption modeling Class I Area State Supervising Agency Upper Buffalo Wilderness Area AR U.S. Forest Service Caney Creek Wilderness AR U.S. Forest Service
Hercules Glade Wilderness Area MO U.S. Forest Service Mingo NWS MO U.S. Fish and Wildlife Sipsey Wilderness Area AL U.S. Forest Service


## III. BART-Eligible Sources

The BART requirements in the RHR are inte nded to reduce emissions specifically from large emission units that, due to age, were exempted from other control requirements of the CAAA. For an emissions unit to be cons idered eligible for BART, it must fall into one of 26 specified categories, must have the potential to emit at l east 250 tons per year of certain haze-forming pollutants, and must have been in existence on August 7, 1977, but not in operation before August 7, 1962.

ADEQ staff determined Arkansas’ BART-eligib le sources by first identifying which of Arkansas’ stationary sources fit the first crit eria of being listed in the BART 26 specific categories. After identifying the sources which fit the first criteria, a database search of5 these facilities was performed to determine whether or not these emitting units’ potential to emit were at least 250 tons pe r year of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ), and particulate matter (PM). The next stage of determining BART-eligibility was to research the permit applications for the year the point source was placed into operation. The final step in the process was to contact facilities for the exact date of operation especially for sources that were placed into operation in the years 1977 and 1962. Tables 2 and 3 contain the list of BART-e ligible facilities ( 18) by BART source category and the number of BART-eligible emitting units (27) within each facility. Figure 4 is a map of Arkansas which shows the location of the 18 BART-eligible facilities located in Arkansas. Figure 5 depicts the five Class I areas Arkansas will be assessing and the BART-eligible sources in Arkansas. Appendix B contains maps showing the receptors at each Class I area ADEQ will be assessing.

Table 2 Fossil fuel-fired steam electric plants > 250 MMBtu/hr and Kraft pulp mills facilities with BART-eligible emission units BART Source Category Number and Name Facility Name/Location Facility
ID AFIN Unit
ID Unit Description

American Electric Power (SWEPCO)/Gentry 05-007- 00107 04- 0017 SN-01 Boiler

AR Electric Cooperative/Augusta 05-147- 00024 74- 00024 SN-01 Boiler 1350mm

AR Electric Cooperative/Camden 05-103- 00055 52- 00055 SN-01 Boiler

Entergy – Lake Catherine/Jones Mill 05-059- 00011 30- 00011 SN-03 Unit 4 Boiler

Entergy – Ritchie Plant/Helena 05-107- 00017 54- 00017 SN-02 Unit 2

Entergy – White Bluff/ Redfield 05-069- 00110 35- 00110 SN-01 Unit 1

Entergy – White Bluff/Redfield 05-069- 00110 35- 00110 SN-02 Unit 2 1. Fossil fuel-fired Electric Plants > 250 MMbtu/hour – Electric Generating Units (EGUs)

Entergy – White Bluff/Redfield 05-069- 00110 35- 00110 SN-05 Auxiliary Boiler

Domtar, Inc./Ashdown 05-081- 00002 41- 00002 SN-03 #1 Power Boiler

Domtar, Inc./Ashdown 05-081- 00002 41- 00002 SN-05 #2 Power Boiler

Delta Natural Kraft/Pine Bluff 05-069- 00017 35- 00017 SN-02 Recovery Boiler

Evergreen/Pine Bluff 05-69- 00016 35- 00016 SN-04 #4 Recovery Boiler Georgia – Pacific Paper/Crossett 05-003- 00013 02- 00013 SN-22 9A Boiler

Green Bay Packing/ Morrilton 05-029- 00001 15- 00001 SN05A Recover Boiler 3. Kraft Pulp Mills

Potlatch/McGehee 05-041- 00036 21- 00036 SN-04 Power Boiler6 Table 3 Petroleum refineries, sintering plants and chemical proce ssing plant facilities with BART-eligible emissions units
BART Source Category Number and Name Facility Name/Location Facility
ID AFIN Unit
ID Unit Description

11. Petroleum Refineries Lion Oil/El Dorado 05-139- 00016 70- 00016 SN809 #7 Catalyst Regenerator 19. Sintering Plants Big River Industries /West Memphis 05-035- 00082 35- 00082 SN-01 Kiln A

Albermarle – South Plant/Magnolia 05-027- 00028 14- 00028 SR-01 Tail Gas Incinerator
Albermarle – South Plant/Magnolia 05-027- 00028 14- 00028 BH-01 Boiler #1

Albermarle – South Plant/Magnolia 05-027- 00028 14- 00028 BH-02 Boiler #2

FutureFuel Batesville 05-063- 00036 32- 00036 6M01- 01 3 Coal Boilers

El Dorado Chemical/El Dorado 05-139- 00040 70- 00040 SN-08 West Nitric Acid Plant
El Dorado Chemical/El Dorado 05-139- 00040 70- 00040 SN-09 East Nitric Acid Plant 21. Chemical Processing Plants

El Dorado Chemical/El Dorado 05-139- 00040 70- 00040 SN-10 Nitric Acid Concentrator7 Figure 4 Map indicating the locations of U pper Buffalo, Caney Creek and the eighteen BART-eligible facilities located in Arkansas8 Figure 5 Map indicating the locations of U pper Buffalo, Caney Creek, Hercules Glade, Mingo and the eighteen BART-eligible facilities located in Arkansas



## IV. CAIR and Arkansas

The Clean Air Interstate Rule was finalized in May 2005 by EPA and applies to states in the eastern U.S. Reconsiderations were fi nalized March 2006. This rule address air pollution transport across state borders. EPA determined which states must reduce which pollutants based on modeling which showed how the travel of pollution affects nonattainment in other states. CAIR requires states to reduce NO x and/or SO 2 emissions. Of the three programs in CAIR, Arkansas is re quired to participate in only the Ozone-Season NO x reductions program. Although EPA’s BA RT Modeling Guidance allows CAIR states to participate in the CAIR cap and trade program, the state of Arkansas is not eligible for the aforementioned trading progr am because Arkansas is in CAIR only for NO x during the ozone season. Therefore, in Arkansas CAIR is not better than BART. Thus BART-eligible EGUs will be modele d for BART determination/exemption by ADEQ.9 

## V. BART Air Quality Modeling Approach

According to EPA’s BART Modeling Guid ance, “CALPUFF is the best regulatory modeling application currentl y available ... and is curre ntly the only EPA-approved model...” (p 45); therefore, ADEQ and CENRAP have chosen to use CALPUFF in the BART determination process as well as in the post-control analys is. One of the air quality modeling approaches suggested by EP A in the BART guidan ce is an individual source attribution approach. This is the appr oach ADEQ proposes to take. Specifically, this entails modeling source-sp ecific units and comparing mode led impacts to a particular deciview threshold (described below). ADEQ has decided to conduct the subject-toBART modeling, rather than have each BART -eligible facility either conduct the modeling or hire a contractor. This plan will eliminate the need for ADEQ to quickly review many air quality mode ling analyses conducted usi ng varying approaches. This plan will also satisfy the need to use a c onsistent approach among the modeling analyses. Once the subject-to-BART modeling is co mplete, all the modeling inputs will be available to facilities subject to BART for them or their consultants to conduct modeling for making BART analyses.

ADEQ will follow EPA’s BART Modeling Guidance (p 42) in sitting a threshold limit in determining whether a BART-eligible source is either subject-to-BART or exempt. According to the aforementioned modeling guidance, an individual source will be considered to “cause visibility impairment” if the emissions results in a change (delta Δ ) in deciviews (dv) that is greater than or equa l to 1.0 deciview on the visibility in a Class I area. Additionally, if the emissions from a sour ce results in a change in visibility that is greater than or equal to 0 .5 dv in a Class I area the sour ce will be considered to “contribute to visibility impairment” (B ART Final Rule, 40 CFR 51 p 39113). Thus, ADEQ has set the threshold limit at 0.5 dv

.
The modeling approach discussed here is specifically designed for conducting the subject-to-BART screening analyses. There may be differences between modeling for conducting BART analyses and that for conducti ng a visibility analysis for a New Source Review permit, which may involve similar em ission sources and the same air dispersion model used here.

To ensure that no sources pass the screening test when they should fail, the simple approach, by its nature, must be the most c onservative of all the c onditions likely to be examined for the source in question. Fo r example, many factors influence the contribution of a source to the Class I area other than distance. The frequency of winds transporting the pollutants toward the Class I ar ea may often be important to include for a reliable screening analysis . Also, a more distant Class I area downwind in the predominant wind direction from a source may receive a higher visibility impact than a closer Class I area that is infrequently downwind of the source. Another example of conservatism in the screening process is th e use of the latest beta version of the CALMET/CALPUFF modeling system using the no-observation (no-obs) mode (the prognostic meteorological model MM5). Thus, the maximum impact instead of the 98 th percentile will be used to determine if a sour ce has an impact on visibility in a Class I.10 Additionally, the BART analysis process includes several othe r steps in addition to the modeling described in this protocol (EPA, 2005). These steps, none of which are addressed in this document, include detailed analysis of:

 Costs of compliance among the various retrofit control options

 Energy and non-air quality impacts

 Existing pollution control technologies in use at the BART-eligible unit particularly with respect to their a ffecting the choice of retrofit options

 Remaining useful life of the units and/or facility

 Improvements in visibility expected from the use of BART controls.



## VI. BART-Eligible Units Physical Parameters

The physical characteristics of the BART-elig ible point sources to be used for the screening stage one analysis will be pr ovided by ADEQ staff. For the stage two screening analysis, ADEQ staff will work w ith the BART-eligible facilities in the development of actual emissions.



### A. Stack Parameters

Stack parameters required for modeling BART -eligible units were extracted from the permit applications. Stack parameters include height of the stack opening from ground in meters, inside diameter in meters, exit veloci ty in meters per second, exit gas temperature in Kelvin, ground elevation of the stack base in meters, and location coordinates of the stack in Lambert Conformal Conical (LCC). The stack coordinates were taken (in Universal Transverse Mertcator, UTM, and then converted to LCC) by ADEQ staff and then verified using ArcMap. Because the BART modeling focuses on mesoscale transport to Class I areas, other source term parameters (needed to calculate localized impacts) such as building heights and widths for calculating downwash will not be used. Appendix C contains tables indicating the stack parameters and coordinates for each BART-eligible emitting unit.



### B. Emission rates

ADEQ notified by email the BART-eligible facilities to provide the 24-hour average actual emission rate with normal operations from the highest emitting day of the year. Excluded from consideration are days where st art-up, shutdown or ma lfunctions occurred unless these activities are re gular, frequently occurring components of the source’s operation cycle.

ADEQ does not intend to use emissions of VOCs and ammonia from facilities for subject-to-BART analysis. Only specific VOC compounds form secondary organic aerosols that affect visibility. These compounds are a fraction of th e total VOCs reported in the emissions inventory, and ADEQ does not have the breakdown of VOC emissions11 necessary to model those that only impair visibility. Furthe r, the prescribed screening model (CALPUFF) cannot simulate formation of particles from anthropogenic VOCs, nor their visibility impacts. Ammonia from sp ecific sources will not be evaluated in this process, although ammonia is included in the modeling as a background concentration— this will be discussed later in this modeling protocol. The appropriate VOCs and ammonia emission data can, and will be, incl uded in regional scale modeling used for the Regional Haze SIP.


## VII. Air Quality Model and Inputs

As stated in the previous section, CALPUFF is the prefer red regulatory air dispersion model for long distance and th erefore is the model ADEQ will be using in the BART determination process. ADEQ recognizes that CALPUFF has limited ability to simulate the complex atmospheric chemistry involved in the estimation of secondary particulate formation. However, for purposes of the subj ect-to-BART analysis, ADEQ intends to use CALPUFF for the following reasons:

1. The increased level of effort required fo r conducting particulate apportionment in the regional scale, full-chemistry Eulera in model (CAMx) to acquire individual source contributions to Cla ss I areas, relative to the simplicity of the CALPUFF model 2. The lack of a plume-in-grid feature with the particulate apportionment technique currently available in CAMx 3. The desire to be consistent with othe r CENRAP states, which all (except Texas and Iowa) appear to be using CALPUFF 4. The limited scope of what this modeling is to determine 5. The additional modeling of BART controls that will be conducted as part of the Regional Haze SIP with the CAMx or CM AQ model(s). EPA’s BART guidance states that States should follow the EPA’s Interagency Workgroup on Air Quality Modeling (IWAQM) guidance, Phase 2 recommendations for long-range transport. The IWAQM guidance was devel oped to address air quality impacts as assessed through the Prevention of Signifi cant Deterioration (PSD) program at Class I areas, where the source generally is located beyond 50 km of the Class I area. The IWAQM guidance doe s not specifically addre ss the type of assessment that will occur with the BART analysis.

EPA recommends in their BART modeling gu idelines (2005) that States follow theInteragency Workgroup on Air Quality Modeling (IWAQM) Phase II (1998) for longrange transport. The IWAQM guidance was developed to address air quality impact – as assessed through the Prevention of Significant De terioration (PSD) program – at Class I areas, where the source generally is locat ed beyond 50 km of the Class I area. The IWAQM guidance does not specifically address the type of assessment that will occur with the BART modeling.12 

### A. CALPUFF Screening Modeling:

CALPUFF modeling will be pe rformed on all Arkansas BART-eligible sources. ADEQ intends to closely follow the CENRAP BART mo deling protocol for most of the settings and inputs. Kansas attempted puff splitting and found this method to be computationally prohibitive on the current domain (State of Ka nsas, 2006). Also, according to Tesche, et al (2005),

“There is no quantitative evidence that th e horizontal and vertical puff-splitting algorithms in CALPUFF yield improved accuracy and precision in model estimates of inert or linearly react ive pollutants although conceptually the methods have appeal in that they attemp t to mimic lateral and vertical wind speed and direction shears.” (p 6-6)

Therefore, ADEQ will not invoke puff splitting in the no-obs screening analysis nor in the refined screening analysis. However, if a potentially subject-to-BART facility wishes to invoke the puff splitting mode, they will be required to notify ADEQ in writing of their intent and provide a protocol for appr oval prior to performing the analysis.



### 1. Modeling domain The CALPUFF modeling will be conducted on the CENRAP central 6 km grid. The extent of the proposed CALPUFF domain is shown in Figure 2.

Figure 6 6 km CENRAP Central CALPU FF domain (Tesche, et al, 2005)13 CALPUFF will be applied to each source for three annual simulations spanning the years 2001 through 2003. The IWAQM Phase II guidance allows the use of fewer than 5 years of meteorological data if a meteorological model using four-dimensional data assimilation is used to supply data. This is the case in this modeling analysis. See the section on meteorology for more information.

2. CALPUFF system implementation There are three main compone nts to the CALPUFF model:

1. Meteorological Data Modeling (CALMET); 2. Dispersion Modeling (CALPUFF); and 3. Post-proce ssing (CALPOST)

Versions of the modeling components to use in this BART analysis are shown in Table 4. Table 4 CALPUFF Modeling Components

Processor Version Level
TERREL 3.311 030709
CTGCOMP 2.42 030709

CTGPROC 2.42 030709
MAKEGEO 2.22 030709
CALMM5 2.4 050413
CALMET 5.53a 040716
CALPUFF 5.753 051130
POSTUTIL 1.4 040818
CALPOST 5.6392 051130 The specific use of each of these components in the BART analysis is described in more detail below.

For screening applications, ADEQ will use the VI STAS version which is the latest ‘beta’ versions of the CALMET/CALPUFF mode ling system. Note that these are not the EPA guideline codes but rather an updated version containing recent (as of this writing) science improvements and bug fixes. The curre nt guideline CALPUF F code is version 5.7, level 030402. This substitution results fr om EPA phasing out the use of the legacy Pasquill-Gifford (P-G) dispersion paramete rs with the introduction of AERMOD as a new guideline model. CALPUFF employs the AERMOD turbulence-based dispersion coefficients and probability density function (pdf) dispersion methods scheme instead of P-G.

The appropriate model codes may be downloaded from www.src.com or purchased with the latest graphical user in terface (GUI) from the model de veloper. The sequence of model processors listed in Table 4 corres ponds to the order in which the programs are typically run.

3. Meteorological data modeling (CALMET) ADEQ will use the 2001-2003 CENRAP deve loped no-obs CALMET dataset for the screening analysis.* This decision was ba sed on EPA Regions VI and VII written comments on the CENRAP BART Modeling G u idelines (Tesche, et al, 2005) which
state,14 “Normally, in accordance with Section 8.3.1.2 (d) of the Guideline on Air Quality Models, the EPA would require that observations be in corporated in conjunction with prognostic meteorological data. Wh ile the idea of use of prognostic data alone holds promise, it is our opinion th at this option require s further evaluation to insure that this approach does not bias CALPUFF towards underestimation (Guideline on Air Quality Models, S ection 3.2.2 (d)(iv)). While we have significant concern regarding the use of the CALMET fields as they have been developed under the procedures documented in this protocol, we would consider the use of the CALMET meteorological fields provided the screening methodology described in Section 6.1 of the protocol is strictly adhered to. In this case, we feel that the use of the maximu m visibility impact rather than the 98 th percentile value is conservative in its application, a nd would overcome concerns of a potential bias towards underprediction [sic] of the “no-observation” mode. Under these circumstances, we would c onsider the use of the CALMET fields acceptable for the CALPUFF screening procedure.” (EPA, 2005)

As stated in Section V. BART Air Qu ality Modeling Approach, ADEQ will use the maximum impact instead of the 98 th percentile to determine if a source has an impact on visibility in a Class I.

However, subject-to-BART facilities have the option of using the CENRAP CALMET processed data or incorporating obser vational meteorological data into the aforementioned CALMET data. If a subject- to-BART facility opts to use the CENRAP CALMET processed data, then the facility will be required to use the maximum impact instead of the 98 th percentile (8 th day). If a subject-to-BART facility decides they would rather use the 98 th percentile, then the facility will be required to incorporate observational data and provide a protocol as well as a performance evaluation which will need to be approved by ADEQ, EPA, and the FLMs.

Appendix F contains the detailed informati on on all CALMET setting that was used to develop the post-processe d no-obs data fields.

4. Dispersion modeling (CALPUFF) The CALMET output is used as input to the CALPUFF model, which simulates th e effects of the meteorological conditions on the transport and dispersion of pollutants from an individual source. In gene ral, ADEQ proposes to use the recommended default options in the CALPUFF model. There are some devia tions, which are discussed below. Table 5 indicates the species that will be modeled an d/or emitted in the no-obs and refined BART analyses.15 Table 5 Species modeled in BART screening analyses Species Modeled Emitted Dry Deposited

SO 2 Yes Yes Computed-gas
SO 4 -2 Yes No Computed-particle
NO
x Yes Yes Computed-gas
HNO 3 Yes No Computed-gas
NO 3 - Yes No Computed-particle PM-fine* Yes Yes Computed-particle
PM 10 * Yes Yes Computed-particle *Please refer to Section VI subsection B for a detailed discussion on PM-fine and PM 10

.
Emissions Speciation: ADEQ does not intend to model sulfate (SO 4 -2 ), nitrate (NO 3 - ), elemental carbon (EC), and secondary orga nic aerosols (SOA) during the screening analyses. However, ADEQ recognizes the impa ct EC and SOA have on visibility. For instance, the light extinction ( β ext ) coefficient for EC is 10 and for SOA it is 4. Currently, data are quite limited on appropriate speciation of organic/inorganic and filterable/condensable emissions by source category. Although there are speciation profiles available for gas- and oil-fired combustion turbines and coal combustion processes, currently there are no detailed profiles for the full range of BART-eligible sources. Thus, in the case of a subject- to-BART source where the PM profile for SO 4 2-

, EC, and SOA are known, ADEQ recommends the aforementioned species be modeled as separate species in CALPUFF in the post-control modeling analysis.

Condensable Emissions : According to Tesche, et al (p 6-5 2005), “condensable emissions are considered primary fine partic ulate.” ADEQ is aware of the inability to
measure PM 2.5 emissions. Thus, BART-eligible fac ilities will be most likely use AP-42 emission factors to develop the “actual” highest average 24-hr emission rate for this pollutant. In the development of this emissi on rate, ADEQ will require these facilities to use the AP-42 emission factors for condensable PM 2.5 . For sources where AP-42 factors are not available, assumptions for partitio ning need to be resolved with ADEQ.

Size Classification of Primary PM Emissions : Particle size parameters are entered in the CALPUFF input file for dry deposition of partic les. There are default values for “aerosol” species (i.e., SO 4 -2 , NO 3 - , and PM 2.5 ). The default value for ea ch of these species is 0.48 μm geometric mass mean diameter and 2.0 μm geometric standard deviation. The main sources of these particles are fuel com bustion. A way to account for this, without including EC and SOA in the modeling, is to use particle speciation in the postprocessing step. This is discussed below in the CALPOST section.

As stated in a prev ious section, all PM 10 emissions will be modeled as PM 2.5 for the noobs model simulations (Tesche, et al, 2005).

Background Ozone concentrations: Ozone (O 3 ) can be input to CA LPUFF as hourly or monthly background values. Hourly values of ozone concentrations were obtained from two rural monitoring sites in Arkansas: Deer , Newton County monitoring site and Eagle Mountain, Montgomery County monitoring site. The hourly ozone concentrations were adjusted for the time differences between the post-processed pr ognostic meteorological16 file (0 GMT) and the collection time of the ozone (LST). Also, the concentrations were adjusted from parts per million (ppm) to parts per billion (ppb). These hourly ozone values will be used in this modeling.

Background Ammonia concentrations: Background ammonia concentration is assumed to be temporally and spatially invariant a nd will be fixed at 3 ppb across the entire domain for all months. It may be possible to derive NH 3 concentrations from regional modeling outputs that CENRAP is curre ntly developing. At this time these NH 3 values are not available in a model ready form.

Receptors: Receptors are locations where model results are calculated and provided in the CALPUFF output files. Receptor locations were derived from the National Park Service’s Class I area receptor database at http://www2.nature.nps.gov/ai r/maps/receptors/index.cfm . Only these discrete NPS receptors will be modeled in CALPUFF. Th e discrete receptors are necessary for calculating visibility impacts in the nine sel ected Class I areas that will be evaluated by ADEQ. All the discrete receptors will be pl aced with enough dens ity that the highest visibility impacts should be evident. The NPS provides recepto rs in all the Class I areas on a 1 km basis. These receptors will be kept at the 1 km spacing for the BART modeling, and all receptors will be retained. NPS also provides a conversion program to convert the coordinates of the receptors fr om latitude/longitude (lat/long) to Lambert Conformal Conical (LCC). ADEQ used this conversion program to convert the receptors located in the five Class I areas it is assessing from lat/long to LCC.

Outputs: The CALPUFF modeling resu lts will be displayed in units of micrograms per cubic meter (μg/m 3 ). In order to determine visibility impacts, the CALPUFF outputs must be post-processed.

Detailed information on all CALPUFF setting to be used in this sc reening analysis is located in Appendix G.

5. Post-processing (POSTUTIL/CALPOST) Hourly concentration out puts from CALPUFF are pro cessed through POSTUTIL and CALPOST to determine visibility cond itions. Specifically, PO STUTIL takes the concentration file output from CALPUFF and recalculates the nitric acid and nitrate partition based on total available sulfate and ammonia. The ammonia-limiting method (ALM) in CALPUFF repartitions nitric acid and nitrate on a receptor-by-receptor and hour-by-hour basis to account for the mode ls systematic over-prediction due to overlapping puffs. For both screening applic ations, the parameter MNIRATE=1 is set in POSTUTIL to implement this approximat e correction in its simplest form. The background ammonia concentration that was obtained from CENRAP’s regional modeling effort will be used to maintain re gional consistency in the CENRAP region. CALPOST uses the concentration file proc essed through POSTUTIL, along with relative humidity (RH) data, to perform visibility ca lculations. For the BART analysis, the only modeling results out of the CA LPUFF modeling system of in terest are the visibility impacts.17 Please see Appendix H and I for detailed settings for POSTUTIL and CALPOST.

Light extinction: Light extinction must be computed in order to calculate visibility. CALPOST has seven methods fo r computing light extinction. This BART screening analysis will use Method 6, which computes extinction from speciated particulate matter with monthly Class I area-specific relative hum idity adjustment factor s, and is implied by the BART guidance. Relative humidity (RH) is an important factor in determining light extinction (and therefore visibility) because SO 4 -2 and NO 3 - aerosols, which absorb moisture from the air, have greater extinc tion efficiencies with greater RH. All BART analyses will apply relative humidity correction factors ( f

(RH) s) to SO 4 -2 and NO 3 - concentrations outputs from CALPUFF, which were obtained from EPA’s “Guidance for Estimating Natural Visibility Conditions under the Regional Haze Rule (EPA, 2003). The f

(RH) values for the Class I areas that will be assessed are provided in Table 5.

Table 6 EPA recommended monthly averaged f ( RH ) for the five Class I areas ADEQ is assessing (EPA, 2003) Class I Area Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Caney Creek 3.4 3.1 2.9 3.0 3.6 3.6 3.4 3.4 3.6 3.5 3.4 3.5
HerculesGlades 3.2 2.9 2.7 2.7 3.3 3.3 3.3 3.3 3.4 3.1 3.1 3.3 Mingo 3.3 3.0 2.8 2.6 3.0 3.2 3.3 3.5 3.5 3.1 3.1 3.3 Sipsey 3.3 3.0 2.8 2.7 3.1 3.4 3.5 3.5 3.5 3.3 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

The PM 2.5 concentrations are considered part of the dry light extinction equation and do not have a humidity adjustment factor. The lig ht extinction equation is the sum of the wet
SO 4 -2 and NO 3 - and dry components PM 2.5 plus Rayleigh scattering ( β Ray ), which is 10 inverse megameters (Mm -1 ).

To account for sources modeled with a known PM speciation profile for EC, SOA, and SO 4 , an adjustment to the extinction coefficient for the PM components will be made in CALPOST. ADEQ intends to follow the method outlined in the FLM CALPUFF Reviewer’s Guide (Gebhart, 2005) wh ich is located in Appendix K.

6. Measuring visibility impacts The recommended procedure for quantifying visib ility impacts can be found in Chapter 3 of the CENRAP BART Modeling Guidelines (Tesche, et al, 2005) which is located in Appendix J. The key point is that the light extinction coefficient ( β ext ) can be calculated from the IMPROVE equation as:

β ext = 3 f ( RH ) [(NH 4 ) 2 SO 4 ] + 3 f ( RH ) [NH 4 NO 3 ] + 4[OC] + 1[Soil] +

+ 0.6[Coarse Mass] + 10[EC] + β Ray18 The monthly site-specific f ( RH ) values were obtained for th e five Class I Area ADEQ is assessing from Table A-3 in the EPA (2003) guidance document. Then, the haze index (HI), in dv, is calculate d in terms of the extin ction coefficient via:

HI = 10 ln ( β ext /10)

The change in visibility (measured in terms of Δ dv) is then compared against background conditions. The Δ dv value is calculated from the source’s contribution to extinction, β source , and background extinction, β background , as follows:

Δ dv = 10 ln ({ β background + β

source }/ β background )

If the Δ dv value is greater than or equal to 0.5 dv , the source is said to contribute to visibility impairment and is thus subject-t o-BART controls. If not, it is BART-exempt.

The annual average natural levels of aeros ol components at each Class I area being evaluated by ADEQ are shown in Table 7. Natural conditions by component in Table 6 are based on whether the Class I area is in the eastern or the western part of the United States. In this BART analysis, all Class I areas are located in the East. The source of this data is from EPA’s Guidance for Estimati ng Natural Visibility Conditions Under the Regional Haze Rule (EPA, 2003).

Table 7 Average annual natural le vels of aerosol components ( μ g/m 3 ) (EPA, 2003) Class I Area Region SO 4 NO 3 OC EC Soil Coarse Mass Caney Creek East 0.23 0.10 1.40 0.02 0.50 3.00 Hercules-Glades East 0.23 0.10 1.40 0.02 0.50 3.00 Mingo East 0.23 0.10 1.40 0.02 0.50 3.00 Sipsey East 0.12 0.10 0.47 0.02 0.50 3.00 Upper Buffalo East 0.23 0.10 1.40 0.02 0.50 3.00

As stated in section V, in a cooperative agreement with EPA Regions VI and VII and FLMs, CENRAP guidance deviates from use of the 98 th percentile impact. The CALMET datasets as described in this protoc ol were processed w ith the no-obs options (i.e., surface observations were not used in the CALMET wind field interpolation). Aware that exercising CALMET with no-obs may lead in some applications to potentially less conservatism in the CALPUFF vi sibility results compared with the use of CALMET with observations, CENRAP has agre ed to EPA’s recommendation that the maximum visibility impact, rather than the 98 th percentile value, should be used for the no-obs screening analysis using the CENRAP-developed CALMET datasets.

If the no-obs screening analys is results indicate a BART-e ligible facility’s maximum Δ dv on a Class I area is less that 0.5 dv, then they will be considered exempt from BART and will be notified by ADEQ of their status. However, if the maximum Δ dv is equal to or greater than 0.5 dv, the source will be c onsidered to be subject-to-BART. ADEQ will notify these subject-to-BART facilities.19 

## VIII. Change in Visibility Due to BART Controls

Once a facility is determined to be subj ect-to-BART, this facility must perform an engineering analysis and a post-control m odeling analysis using CALPUFF. This modeling analysis must be compared to the pr e-control modeling result s. Please note that this will be a source specific (i.e. emitting unit specific) and pollutant specific modeling analysis using CALPUFF. If a subjec t-to-BART facility opts to use the 98 th percentile rather than the maximum impact, the subjec t-to-BART facility will be required to be incorporate observational data with the post processed CALMET prognostic meteorological data. Also these facilities wi ll be required to submit their meteorological modeling protocol, model performance eval uation, and CALPUFF modeling protocol to ADEQ, EPA Region VI, and FLMs for approva l. However, if the subject-to-BART facility opts to use the maximum impact rather than the 98 th percentile, these facilities may use the post-processed CALMET MM5 data.

Additionally, one control measure that a source may opt to use is to revise their Title V permit to provide for synthetic minor limits so that it falls under the BART emission cap. That permit modification must be done prior to the State going to public hearing on its RH SIP. The limits must be in place for as long as the RH SIP is applicable or for as long as the source is operational. However, the source will still need to do a post-control CALPUFF modeling analysis to determine the amount of emissi ons it needs to reduce for visibility improvement. (Note: ADEQ strongly recommends that all subject-to-BART facilities work closely with ADEQ in their engi neering analyses.) Also, after all of the post-control results are submitted to and approved by ADEQ, these results will then be inputted into either CAMx or CMAQ for a cumulative model run. If the control measures proposed by the BART facilities still impact a Class I area, the BART facilities will need to implement additio nal control. Please note that all post-modeling results are due to ADEQ no later than October 23, 2006.20 

## IX. References

Gebhart, Howard (2005). Draft CALPUFF REVIEWER’S GUIDE

. Air Resource Specialists, Inc., Fort Collins, Colorado.

Kansas Department of Health and Environm ent, Bureau of Air and Radiation (March, 2006). Best Available Retrofit Technology (B ART) Modeling Protocol to Determine Sources Subject to BART in the State of Kansas

Minnesota Pollution Control Agency (2005). Best Available Retrofit Technology (BART) Modeling Protocol to Determine S ources Subject to BART in the State of Minnesota.

Scire, J.S., D.G. Strimaitis, and R.J. Yamartino. (2000). A User’s Guide for the CALPUFF Dispersion Model (Version 5). Earth Tech, Inc., Concord, Massachusetts.

Scire, J.S., D.G. Strimaitis, and R.J. Yamartino. (2002). A User’s Guide for the CALMET Dispersion Model (Version 5). Earth Tech, Inc., Concord, Massachusetts.

Tesche, T.W, et al (2005) CENRAP BART Modeling Guidelines. Alpine Geophysics, LLC. Ft. Wright, KY

U.S. EPA. (1998). Interagency Workgroup on Air Quality Modeling (IWAQM) Phase2—Summary Report and Recommendations for Modeling Long Range Transport Impacts

. EPA-454/R98-019.

U.S. EPA (2003). Regional Haze: Estimating Natural Visibility Conditions Under the Regional Haze Rule. EPA-454/B-03-005.

U.S. EPA. (2005). Regional Haze Regulations and Guidelines for Best Available Retrofit Technology (BART) Determinations. Federal Register Vol. 70, No. 128.

U.S. EPA, Regions VI & VII ( 2005). BART Protocol Comments.

U.S. EPA, Region VI (2006). Personal communication.A-1 

## Appendix A. Proposed Time-line

Figure A-1 ADEQ’s proposed time-line to m eet the RH SIP deadline of December 17, 2007 as set forth by EPA in its Regional Haze Regulations and Guidelines for Best Available Retrofit Technology (BART) De termination (40 CFR Part 51, p 39156)B-1 

## Appendix B. Map of receptors

Figure B-1 Receptors located in Cane y Creek Wilderness Area, Arkansas

Figure B-2 Receptors located in Herc ules-Glade Wilderness Area, MissouriB-2 Figure B-3 Receptors located in Mingo Wilderness Area, Missouri

Figure B-4 Receptors located in Sipsey Wilderness, AlabamaB-3 Figure B-4 Receptors located in Uppe r Buffalo Wilderness Area, ArkansasC-1 

## Appendix C. BART-Eligible Sour ces’ Stack Parameters, Base Elevation at Ground Level, and Stack Coordinates

Table C-1 BART-eligible sources’ stack parameters

STATIONARY SOURCE NAME/LOCATION Emission Unit ID Stack Height Meters Stack Diameter Meters Exit Velocity m/sec Temperature

o K Albermarle-South Plant / Magno lia SR-01 57.9 0.814 15.24 922 Albermarle-South Plant / Magno lia BH-01 6.4 2.591 9.14 505 Albermarle-South Plant / Magno lia BH-02 6.4 2.591 9.14 505 American Elect. Power (SWEPCO) / Gentry SN-01 164.6 6.096 34.14 408 AR Elect. Coop - Bailey Plant / Augusta SN-01 50.9 3.000 28.04 444 AR Elect. Coop - McClellan Plant / Camden SN-01 48.8 3.301 28.04 444 Big River Industries / W. Memphis SN-01 30.5 1.524 21.88 330
Delta Natural Kraft / Pine Bluff SN-02 50.3 2.134 13.29 348 Domtar, Inc. / Ashdown SN-03 66.1 1.890 26.76 522 Domtar, Inc. / Ashdown SN-05 71.6 3.659 11.92 325
FutureFuel / Batesville 6M01-01 61.0 2.743 9.45 422
El Dorado Chemical / El Dorado SN-08 22.9 1.219 33.53 505 El Dorado Chemical / El Dorado SN-09 22.9 1.219 32.00 500
El Dorado Chemical / El Dorado SN-10 23.8 0.152 23.77 313 Entergy - Lake Catherine / Jones Mill SN-03 59.4 5.182 3.08 396 Entergy - Ritchie Plant / Helena SN-02 71.9 3.658 28.62 390 Entergy - White Bluff / Redfield SN-01 304.8 7.833 27.43 434 Entergy - White Bluff / Redfield SN-02 304.8 7.833 27.43 434 Entergy - White Bluff / Redfield SN-05 4.6 0.914 19.81 519 Evergreen/Pine Bluff SN-04 48.8 2.438 14.33 474 Georgia-Pacific Paper / Crossett SN-22 53.3 3.658 10.45 341
Great Lakes Chemical / El Dorado SN-302A 9.1 0.762 40.54 555 Green Bay Packaging / Morrilton SN-05A 30.8 1.798 25.60 456
Lion Oil / El Dorado SN-809 61.0 1.753 9.75 533 Potlatch Corp. / McGehee SN-04 89.6 2.743 14.78 444C-2 Table C-2 BART-Eligible Emission Units’ Base Elevation and Lambert Conformal Conical (LCC) Coordinates

STATIONARY SOURCE NAME/LOCATION (BART File Name) Emission UNIT ID Base Elevation,
meters (m) X Easting LCC x Y Northing LCC y Albermarle-South Plant / Magnolia SR-01 86.9 352.81836 -747.03381 Albermarle-South Plant / Magnolia BH-01 88.4 352.67618 -746.98114 Albermarle-South Plant / Magnolia BH-02 88.4 352.65801 -746.98190 American Elect. Power (SWEPCO) / Gentry SN-01 349.9 221.58128 -410.39077 Ark. Elect. Coop - Bailey Plant / Augusta SN-01 61.3 510.86643 -507.71488 Ark. Elect. Coop - McClellan Plant / Camden SN-01 33.5 390.21870 -702.15534 Big River Industries (General Shale)/ W. Memphis SN-01 60.0 609.12652 -517.70639
Delta Natural Kraft / Pine Bl uff SN-02 66.4 457.00824 -621.20692 Domtar, Inc. / Ashdown SN-03 97.5 267.47491 -698.66686 Domtar, Inc. / Ashdown SN-05 97.5 267.48245 -698.74355
FutureFuel/ Batesville 6M01- 01 82.3 493.14724 -458.02938
El Dorado Chemical / El Dorado SN-08 63.1 401.11728 -734.65321
El Dorado Chemical / El Dorado SN-09 63.1 401.13533 -734.65236
El Dorado Chemical / El Dorado SN-10 62.2 401.19594 -734.67412 Entergy - Lake Catherine / Jon es Mill SN-03 100.0 375.45658 -606.40861 Entergy - Ritchie Plant / Helena SN-02 54.9 586.25363 -591.07129 Entergy - White Bluff / Redfield SN-01 94.2 446.73457 -625.11197 Entergy - White Bluff / Redfield SN-02 94.2 445.61252 -604.15523 Entergy - White Bluff / Redfield SN-05 94.2 445.61539 -604.25671 Evergreen/Pine Bluff SN-04 63.1 468.11868 -625.3954 Georgia-Pacific Paper / Crossett SN-22 46.0 469.03486 -745.02133 Green Bay Packaging / Morrilton SN-05A 98.5 387.29077 -532.44265
Lion Oil / El Dorado SN-809 75.6 403.01817 -741.82948 Potlatch Corp. / McGehee SN-04 43.9 533.13136 -678.59798D-1 

## Appendix D. BART-Eligible Emis sion Rates used for the NoObs Modeling Run

Table D-1 BART-eligible units’ highes t 24-hour actual emi ssion rates for SO 2, NO x , PM 10 * and PM 2.5 in grams per second (g/sec)

Highest 24-Hour Actual Emission Rates (g/sec)

BART-Eligible Facilities/ Locations Emission Unit ID Number SO 2 NO
x PM 10 PM 2.5 Albermarle-South Plant / Magnolia SR-01 48.126 0.076 0.000 0.009 Albermarle-South Plant / Magnolia BH-01 0.353 2.075 0.000 0.136 Albermarle-South Plant / Magnolia BH-02 0.535 2.578 0.000 0.128

American Elect. Po wer (SWEPCO) / Gentry SN-01 595.781 245.066 21.725 5.531 AR Elect. Coop - Bailey Plant / Augusta SN-01 299.344 36.933 21.729 21.729 AR Elect. Coop - McClellan Plant / Camden SN-01 346.189 47.124 28.764 28.764 Big River Industries/ W. Memphis SN-01 0.000 8.589 0.000 7.076

Delta Natural Kraft / Pine Bluff SN-02 0.239 1.701 1.058 0.529
Domtar, Inc. / Ashdown SN-03 0.774 22.632 0.000 21.354
Domtar, Inc. / Ashdown SN-05 70.175 52.008 0.000 7.881

FutureFuell / Batesville 6M01-01 54.046 11.045 0.290 0.217

El Dorado Chemical / El Dorado SN-08 0.000 20.060 0.000 0.000

El Dorado Chemical / El Dorado SN-09 0.000 15.645 0.000 0.000

El Dorado Chemical / El Dorado SN-10 0.000 0.415 0.000 0.000 Entergy - Lake Catherine / Jones Mill SN-03 0.420 309.535 0.365 0.246 Entergy - Ritchie Plant / Helena SN-02 0.105 17.640 0.997 0.997 Entergy - White Bluff / Redfield SN-01 978.164 550.821 15.592 11.802 Entergy - White Bluff / Redfield SN-02 985.933 596.075 16.653 12.915 Entergy - White Bluff / Redfield SN-05 4.095 3.811 0.365 0.246 Evergreen/Pine Bluff SN-04 21.42 28.35 1.008 0.000 Georgia-Pacific Paper / Crossett SN-22 77.275 182.677 0.000 9.310 Green Bay Packaging / Morrilton SN-05A 4.934 8.771 0.000 1.165

Lion Oil / El Dorado SN-809 23.142 5.980 0.000 7.696 Potlatch Corp. / McGehee SN-04 6.942 10.533` 2.752 2.752E-1 

## Appendix E. Chapter 5 of the CENRAP BART Modeling Guidelines (Tesche, et al, 2005)

5.0 DATA BASES FOR CA LPUFF MODELING

To support BART modeling by the states and source operators, both meteorological and aerometric data sets are re quired. Regional meteorological data sets generated by the CALMET model suitable fo r direct input to the CALPUFF modeling system have been developed and archived. These data sets cover calendar years 2001, 2002, and 2003 for three sub-regional grid dom ains shown in Figures 5-1 through 5-4. The procedures used in developing the CA LMET data sets generally follow the IWAQM recommendations (EPA, 1998), except for a few notable refinements. The processed CALMET files, in CALPUFF-ready input fo rmat, are available from CENRAP on hard disk drives to interested states and stakeholders.

This chapter describes how these meteor ological modeling sets were developed and evaluated. The basic CALMET model confi guration used to generate the three years of CALPUFF-ready meteorology is described in detail so that users of this information have a clear understanding of the data sets and their applicability. In addition, for those states or source operators who elect to conduct more source-specific CALMET/CALPUFF modeling, the information in this chapter may be helpful in guiding specification of revised CALM ET model inputs and generati on of revised CALMET data

sets.

Also included in Section 5 .2 is a disc ussion of routinely available air quality monitoring data sets available to the states and source operators in support of screening and source-specific BART modeling exercises.

5.1 Development of CALMET Meteorological Files

5.1.1 MM5 Data Sets

Alpine Geophysics developed a c onsistent set of CALMET regional meteorological modeling data sets for use by the CENRAP States, BART eligible sources within the region and others. These meteorol ogical modeling data sets were constructed through the joint use of the CALMET processo r and results from existing annual threedimensional MM5 meteorological simulati ons. The specific annual prognostic model simulations available for CENR AP BART modeling included:

> 2001 MM5 data set at 36/12 km re solution developed for EPA by Alpine Geophysics (McNally and Te sche, 2002; McNally 2003);

> 2002 MM5 data set at 36 km resolution developed for CENRAP by Iowa DNR (Johnson, 2003a,b),E-2 > 2003 MM5 data set at 36 km resolution developed for the Midwest RPO (Baker, 2005; Baker et al., 2004; Kembell-Cook et al., 2005)

Each of these studies included a performance evaluation of the MM5 generated data sets against surface meteorological observations and the results of these evaluations are contained in the reports or pr esentations cited above. While there exists a set of annual 12 km MM5 meteorology for 2002, this data set was developed by four independent CENRAP modeling centers and th ese data sets have not been concatenated into one master data base. More im portantly, there has been no systematic, rigorous model performance evaluation performed on the CENRAP 2002 12 km MM5 data yet. Accordingly, until such time as the 2002 12 km data set has been evaluated and shown to be of comparable reliability as the aforementioned MM5 data sets, it’s use is contraindicated.

5.1.2 CALMET Model Configuration

The CALMET modeling procedures used to construct meteorological inputs to CALPUFF for visibility scr eening of BART eligible sources generally follows the IWAQM recommendations (EPA, 1998), except as noted below.

CALMET Model Options

. The CALMET model has a numbe r of user-selected options, parameter settings, and ‘switches’ that must be defined prior to exercising the processing system. These options and settings are we ll-described in the CALMET User’s Guide (Scire et al., 2000a) and in the CALMET input file to the executable code. Appendix A of this protocol summarize s the CALMET configurations used in developing the processed 6 km meteorological fields over the three CENRAP BART modeling domains. Also included in the tables in Appe ndix A are the default CALMET options and parameter settings recommended in the IWAQM Phase 2 Report (EPA, 1998).

CALMET Domain. Three slightly overlapping m odeling domains were defined by CENRAP to support BART modeling. These do mains are shown in Figures 5-1 through 5-4 and Table 5-1. The processors used to generate the domain, land use, and elevation data for the CALMET/CALPUFF system include TERREL, CTGPROC, and MAKEGEO, as described below.E-3 > TERREL is the terrain pre-processor that averages terrain features to the modeling grid resolution; TERREL constr ucts the basic pr o perties of the gridded domain and defines the c oordinates upon which meteorological data are stored. Key parameters in clude specification of grid type, location, resolution and terrain elevation.

> CTGPROC computes the fractional land use for the modeling grid resolution. Land use characteristics for each grid cell are assigned usin g CTGPROC. The primary variable adju stment associated with CTGRPOC is selection of an appr opriate land use database. Version 2.0 of the North American Land Cover Characteristics database is used.

> MAKEGEO is the final pre-processor th at comb ines the terrain and land use data for input to CALMET. Generating the appropriate MAKEGEO.INP control file require s only minimal alteration of the default assignments. Key modificat ions include specifying domain attributes and ensuring input f iles are correctly referenced.

Terrain . CALMET requires both te rrain height and land us e/land cover for the application region. These are generated using the CALMET CTGPROC, TERREL and MAKEGEO processors. The terrain data were created using the TERREL (version 3.311, level 030709) processor and the Shuttle Radar Topography Mission (SRTM)-GTOPO 30 second (~1 km) reso lution dataset.

Land Use . The landuse data set was created using the Composite Theme Grid CTGROC processor (version 2.42, level 030709) and the United States Geological Survey (USGS) Global Land Cover Characterization (GLCC) version 2.0 database. The GLCC database is available at 30 second (~1km) resolution. References for these and other modeling datasets can be found at www.src.com

.
Vertical Layer Structure. The vertical layer stru cture for the CALMET/CALPUFF screening applications is more refined than the ge neral suggestions of IWAQM. The CENRAP vertical structure wa s designed to reduce the need for vertical interpolation while simultaneously improving vertical re solution within the pl anetary boundary layer (PBL). Table 5-2 identifies the 11 layer interface s required to define the 10 layer vertical CALMET grid structure. The top interface in the CALMET simulation is 4000 meters.

Use of Observations

. Based on considerable disc ussions with State and Federal managers and agency personnel, CENRAP has elected to use the No-Obs mode in CALMET for constructing the 6 km meteor ological fields for CALPUFF screening exercises. The three annual MM5 simula tions (2001, 2002, and 2003) will be used as the sole source for meteorological data within CALMET. Blending observational data with the MM5 data within CALMET (i.e., use of th e “OBS” option is esse ntially a redundant use of the same data. Substantial improvement in the MM5 initialization data and in the use of four dimensional data assimilation (FDDA) has been achieved in recent years using observational data. Th e ETA analysis data used in initial and boundary conditions estimates as well as within the FDDA fiel ds derive from 3-hourly, 40 km objectiveE-4 analysis fields computed using an extens ive supply of observational data (National Weather Service surface and upper air data, GOES satellite precip itable water; VAD wind profiles from NEXRAD; ACARS aircraft temperature data; SSM/I oceanic surface winds; daily NESDIS snow cover and sea-ic e analysis data; RAOB balloon drift; GOES and TOVS-1B radiance data; 2D-VAR SST fr om NCEP Ocean Modeling Branch; radar estimated rainfall; and surface rainfall). The complexity, resolution, and accuracy of the ETA data that is used to in itialize and ‘nudge’ the MM5 for ecasts is extensive indeed. Particularly at the 12-36 km horizontal grid scales over the flat to modestly rolling topography of the CENRAP domain, there is no need to introduce local meteorological observations in order to retrieve local terrain effects, for example. Thus, mesoscale wind patterns are likely to be adequately characterized by the MM5 simulations.

Many observations, especially surface obser vations, reflect local conditions on a scale smaller than the 6 km CENRAP CALM ET fields. The introduction of the local observations into the regional modeling domain may extend the influence of the observational data beyond its true representati veness and result in internally inconsistent flow features. In particular the time interpolation of the 12-hourly upper air sounding data may wash out structure in the MM5 fields that are appropr iate to retain. Given that the CENRAP domain as a whole includes area s of moderately rolling terrain, coastal regions and relatively flat terrain, a single set of representative weights 1 that allows significant influence of the observations where appropriate, will involve a considerable effort and substantial testing. The internally consistent MM5 fields are considered likely to be appropriate for the regi onal simulations, and the increm ental benefit of adding the observational data into the regi onal CALMET simulations is not considered worthwhile.

However, on the smaller domains likely to be considered in source-specific modeling (e.g., 1-4 km in scale) with the hi gher CALMET grid resolu tion and the smaller domain size, more control over the region of influence of the mete orological observations can be achieved. It is easier for the diagnostic model to a llow the local flow observations to have appropriate influence in the vicinity of the observation, but allow terrain-adjusted flow to dominate away from the observations. Given that the fine scale source-specific domains will be used especially in irregular and/or meteorologically complex settings, the relatively coarser-scale MM5 simulations are less likely to be fully adequate, and the introduction observational data into CALMET is more likely to achieve improvements in the resulting meteorological fields.

Diagnostic Model Settings

A number of diagnostic m odel settings must be selected for CALMET to properly process representative diagnostic meteorologi cal data sets. Thes e are summarized in Appendix A, compared to the default CALMET settings, and discussed in the following:

1 Weights are assigned in CALMET to control the ‘blending’ of observations and MM5 predictions

.E-5 > CALMET options dealing with radius of influence parameters (R1, R2, RMAX1, RMAX2, RMAX3), BIAS, ICALM parameters are not used in No-Observations mode;

> Gridded cloud data were inferred from the MM5 relative humidity fields (ICLOUD=3);

> Given that all state variables are MM5-derived (IPROG=14; ITPROG=2), surface layer winds were not extrapolat ed to the upper layers (IEXTRP = - 1);

> The IWAQM recommendation for disa bling the computation of kinematic effects in the wind field options and parameters was selected. This was selected in light of the very mode st elevated terrain in the CENRAP domain, relative to the mountainous re gions in the U.S. and Alps where the kinematic parameterizations were originally developed. Thus, the option for computing kinematic effects was disabled (IKINE = 0).

> The BIAS array was set to 0. in the CALMET control file because surface and upper air data were not used (NOOBS = 2);

> Because the MM5 wind fields supply CALMET with the initial guess fields to the diagnostic wind mo del (IWFCOD =1, IPROG = 14) and observational data are not reintroduced, the following variables were set to nominal values:

 The minimum distance for which extrapolation of surface winds should occur was set to -1 (RMIN2 = -1.).

 RMIN was left at the IWAQM recommendation of 0.1 km.

 RMAX1 and RMAX2 were each assigned a value of 30 km. RMAX3 was assigned a value of 50 km.

 R1 and R2 were each assigned the value of 1.0.

 ISURFT and IUPT were assigned placeholder values of 4 and 2, respectively.

> The radius of influence regarding terrain features is comparable to the resolution of the processed terrain data: 12 km.

> The radius of influence for temp erature interpolation is set to 36 km (TRADKM), a value considered appr opriate given the 6 km CALMET domain and 36/12 km MM5 domain.

> The beginning/ending land use cate gories for temperature interpolation over water are assigned category 55: (JWAT1 = JWAT2 = 55).E-6 > SIGMAP was set to 50 km, while the IWAQM recommendation is 100 km, but with no supporting documentat ion. Because precipitation rates are explicitly incorporated from the MM5 data, a lower radius of influence was deemed appropriate.

> Diagnostic options: IWAQM default values were used (see Appendix A);

> TERRAD (terrain scale) is required for runs with diagnostic terrain adjustments (i.e., the 2003 simulations ). Values of ~10-20 km were tested, and an appropriate value determined.

> Land use defining water: JWAT1 = 55, JWAT2 = 55 (large bodies of water). This feature allows the te mperature field over large bodies of
water such as the Gulf of Mexico and the Great lakes to be properly characterized by buoy observations.

> Mixing height averaging parame ter (MNMDAV) were determined sensitivity tests. The purpose of the testing is to optimize the variable to allow spatial variability in the mixing height field, but without excessive
noise.

Obviously, there are some instances where mo re advanced and/or recently developed procedures for constructing the CALMET fiel ds have been used compared with the IWAQM (1998) guidance. For example, one ag ency expressed concer n about the choice to employ prognostic model-derived gridded cloud cover data in CALMET (ICLOUD = 3). While this is admittedly a ‘non-guideli ne’ option, in our view it represents the best science option currently available. In par ticular, the EPA CAIR and CAMR rulemaking modeling and the CAMx/CMAQ modeling being performed by the RPOs for regional haze all utilize the gridded moisture fields in the MM5 model as a basis for estimating cloud. Presumably, if the method is suitable fo r such advanced visib ility modeling, it is adequate for CALPUFF modeling. Of course , in the protocol negotiation, the States, source operators, and regulat ory agencies have an oppor tunity to re-examine the CALMET diagnostic model settings used in creating the CENRAP gridded fields and modify them if warranted.

In summary, the development of the regional CALMET meteorological fields from MM5 data was conducted in No-Observations (“No-Obs”) mode. CALMET’s boundary layer modules were used to comput e mixing heights, turbulence coefficients and other meteorological parameters required as input to CALPUFF.

5.1.3 MM5/CALMET Processing

Construction of the CALPUFF-ready mete orological fields en tails a two-step process. First, the MM5 prognostic model out put fields are extrac ted and processed for input to CALMET. This step entails runn ing various extraction software routinesE-7 followed by the CALMM5 code. Then, CALM ET is exercised for the full three year period over each sub-regional CENRAP domain.

CALMM5. Previous applications of the prognostic Mesos cale Meteorological model version 5 (MM5) served as th e source of the gridded meteor ological fields for calendar years 2001, 2002, and 2003. The actual CAL MM5 configuration entailed modification of a few user-specified variables. Howeve r, two setting are of primary importance:

> All vertical layers from MM5 were extracted, providing CALMET configuration flexibility, and

> Vertical velocity, relative humidity, cloud/rain fields, an d ice/snow fields were extracted. (Graupel was extrac ted for 2001, the only year where the data were available in the MM5 datasets.)

CALMET. CALMET (v5.53a, lev 040716) was appl ied consistent with CENRAP’s recommendation that the 6 km be generate d using the ‘No-Obs’ option. The specific options used have been discussed abov e and are summarized in Appendix A.

5.1.4 Evaluation of the CALPUFF-Re ady Meteorological Data Sets

In typical applications th e adequacy of the CALMET fields is seldom evaluated using independent measurements. Often, only cursory visual examin ation of wind vector plots or time series is considered. This evaluation is important because the CALMET performance analysis gives dir ect insight into the adequacy of the model-processed fields on a subregional basis. It also serves as an independent quality assurance tool. Alpine’s MAPS evaluation software to perform an independent evaluati on of the processed CALMET data bases. MAPS was used in conjunction with the NCAR DS472 TDL data sets to evaluate the surface winds and temperatures for 2001-2003 across all three domains. Since only a sma ll portion of the meteorological content of these data were ingested in the MM5 data as similation routines (see Johnson, 2003a), these data sets are essentially an independent, quantitative means for evaluating the adequacy of the meteorological fields input to CALPUFF.

CALMET Evaluation Methodology

Several statistical measures were calculated as part of the CALMET meteorological evaluation using established pr ocedures (e.g., Tesche et al., 1990; Emery et al., 2001). Additional plots and graphs ar e used to present these statistics on both hourly and daily time frames over the full annual cycle. For th is study, evaluation measures were calculated for wind, temperature, and relative humidity because these parameters are the principal meteorological inputs to CALPUFF. The full set of CALMET evaluation statistics and graphica l displays generated with the AG-MAPS software (McNally and Tesche,1994) are cont ained on a DVD available from CENRAP.E-8 The statistics used to evaluate th e meteorological fields for 2001-2003 are generated in both absolute terms (e.g., wind speed error in m/s), and relative terms (percent error) as is commonly done for ai r quality assessments. Obviously, a very different significance is associated with a give n relative error for di fferent meteorological parameters. For example, a 10% error for wind speed measured at 10 m/s is an absolute error of 1 m/s, a minor error. Yet a 10% e rror for temperature at 300 K is an absolute error of 30 K, a ridiculously large error. On the other hand, pollutant concentration errors of 10% at 1 ppb or 10 ppm carry pr actically the same significance.

Three key meteorological metrics include th e bias, error, and index of agreement (IOA) for wind speed, temperature and relative humidity. These measures are defined as follows:

Bias (B) : Calculated as the mean difference in prediction-observation pairings with valid data within a given analysis region and fo r a given time period (hourly or daily):

Error (E) : Calculated as the mean absolute difference in predictio n-observation pairings with valid data within a given analysis region and for a given time period (hourly or

daily).

Note that the bias and gross error for wi nds are calculated from the predicted-observed residuals in speed and directio n (not from vector components u and v ). The direction error for a given prediction-observation pa iring is limited to range from 0 to ± 180 °

.
Index of Agreement (IOA) : calculated following the approa ch of Willmont (1981). This metric condenses all the differences between model estimates and observations within a given analysis region and for a given time pe riod (hourly and daily) into one statistical quantity. It is the ratio of the total RMSE to the sum of two differences – between each prediction and the observed mean, and each observation and the observed mean:

Viewed from another perspective, the inde x of agreement is a measure of the match between the departure of each prediction from the observed mean and the departure of each observation from the observed mean. Thus, the correspondence between predicted 

# ( ) ∑∑

= = − = J

j I

i i j i j OP IJ B

1 1 1



## ∑∑

= = − = J

j I

i i j i j OP IJ E

1 1 1

⎥ ⎥ ⎥ ⎥

⎦ ⎤

⎢ ⎢ ⎢ ⎢

⎣ ⎡

−+− ⋅

−= 

## ∑∑

= = J

j I

i o i jo i j MOMP

RMSEIJ IOA

1 1 2 1E-9 and observed values across the domain at a given time may be quantified in a single metric and displayed as a time series. The i ndex of agreement has a theoretical range of 0 to 1, the latter score suggesting perfect agreement.

CALMET Evaluation Results

Table 5-5 summarizes the statistical measures, averaged over the month, for temperature, wind speed, and relative humidity for all three years. The CALMET evaluation DVD contains a full compilation of the statistic al and graphical results. Figures 5-7 through 5-31 present a variety of graphical displays of processed and observed surface temperature, relative humidity, and wind across the three CENRAP subdomains for the three-year period 2001- 2003. Figures 5-28 through 5-31 provide convenient summaries of the bi as and error in the relative humidity, temperature, and wind speed fields across the conti nuous 36 month period by subdomain.

Thorough discussion of the performance fi ndings is beyond the scope of these guidelines. However, a few ke y findings of the evaluation ar e worth noting here. From Table 5-5, the wind speed index of agreement, a general measure of correlation between measured and observed winds, is systematically greater than a value of 0.8 for virtually every month. These values are typically bette r than those generally achieved in urbanand regional-scale model applications for ozone SIPs. For example, the statistical benchmark for IOA suggested by Emery et al., (2001) is IOA > 0.6. Thus, the wind speed agreement for all three domains and al l three years appears quite good relative to other MM5/RAMS model applic ations. From Figure 5-11, the wind speed root mean square error for the Central domain for 2002 is generally below 2.0 m/s, the performance goal for this parameter. From Figure 5-29 (as well as in Table 5-5), the temperature bias results for the 36 month are ge nerally quite close to the + 0.5 deg C performance goal. As shown in Figure 5-30 the temperature error re sults are slightly poor er than the 2 deg C performance goal for 2001 and 2003, but are below the 2.0 deg C threshold for 2002. Note that the benchmarks were developed not to provide a pass/fail standard to which all modeling results should be held, but rather to put the results in to an historical context.

In summary, we find that:

 Relative Humidity

 Bias over three-year peri od near zero all domains

 For some months over- and under-predi ction (up to 10% or more) is evident – no discernable trend

 Errors typically diminish from 2001 through 2003, and are generally < 12% after 1st quarter of 2001.

 Surface Temperatures

 Monthly averaged temperatures are sy stematically biased low (cooler) by 0.25 to 1.25 deg C.

 The errors in monthly averaged temperatures typically range between 1.8 and 2.6 deg CE-10  Average error over all mont hs is about 2.2 deg C.

 Surface Wind Speeds

 IOA typically between 0.8 0-0.9

 Seasonally variable

 Central subdomain gives best correlation

 Results from MM5/CALMET evaluation prov ide potentially useful information for diagnosing BART visibility modeling analyses

 MM5/CALMET fields exhibit good statisti cal agreement with observations, in part because observations figure prominently in the construction of the interpolated CALMET fields.

 MM5/CALMET fields for the three CENRAP subdomains are quite sufficient for use in CALPUFF modeling.

5.1.5 Meteorological Data Archive and Distribution

All models, scripts and CALMET data (e xcepting MM5 outputs) are available from CENRAP on appropriate external combination Firewire/USB drives.

5.2 Aerometric Monitoring Networks

Data from ambient monitoring networks for both gas-phase and aerosol species are available for use in CENRAP BART modeling analyses. Table 5-4 summarizes ambient monitoring networks. Data for 2002 have been compiled for all networks covering the CENRAP domain with the exce ption of the PAMS and PM Supersites. These data sets may be obtained from CE NRAP. Figures 5-5 and 5-6 display the locations of monitori ng sites in and near the CENRAP States.E-11 Table 5-1. CENRAP Lambert Conic Conf ormal Modeling Domain Specifications (40.97 degree projection origin; 33 and 45 degree matching parallels).

Domain Southwest Coordinate (km) Number of X grid cells Number of Y grid cells Horizontal Resolution

CALMET South -1008, -1620 306 246 6 km Central -1008, -864 388 234 6 km North -1008, 0 300 193 6 km

Table 5-2. Vertical Layer Structure in CALMET Fields. (Heights are in meters.)

L AYER

NUMBER L AYER H EIGHT L AYER

NUMBER L AYER H EIGHT 0 0. 6 640. 1 20. 7 1200. 2 40. 8 2000. 3 80. 9 3000. 4 160. 10 4000. 5 320.

Table 5-3. Meteorological Model File Sizes for CENRAP BART Modeling.

Domain Monthly A nnual 3 Years Domain Grid 3 years North 4.6 55.2 165.6 2001 12 km 1370 Central 6.6 79.2 237.6 2002 36 km 430 South 6.0 72.0 216.0 2003 36 km 430 total 17.2 206.4 619.2 total 2230 CALMET 6 km File Sizes, (Gbytes) MM5 File Sizes, (Gbytes)E-12 Table 5-4. Statistical Evaluation of th e CALMET Meteorological Fields for 2001- 2003.

Jan Feb Ma r Ap r May Jun Jul A ug Sep Oct Nov Dec Mean RH Bias (%)
North 4.54 3.19 0.17 -14.55 -12.09 -4.35 -0.62 1.17 -2.07 -7.98 -6.62 -4.22 -3.62 Central -2.60 -7.28 -11.38 -10.69 -8.62 -2.90 0.66 1.07 -1.44 -5.46 -6.16 -7.78 -5.21 South -10.23 -11.53 -13.78 -4.24 -2.08 0.99 4.12 3.16 -0.12 -2.12 -3.44 -9.76 -4.09 RH Error (%) North 10.06 10.31 14.03 18.77 16.28 12.39 11.82 11.76 13.26 15.54 13.53 12.89 13.39 Central 13.32 15.86 17.45 17.05 14.50 11.67 11.52 11.32 12.26 15.52 14.79 14.95 14.18 South 16.22 18.37 18.17 13.26 12.15 11.51 12.09 12.40 11.82 14.85 14.73 16.19 14.31 Temp Bias ( 0 C)
North -1.63 -1.23 -1.23 -0.24 0.08 -0.29 -0.23 -0.54 -0.55 -0.09 -0.40 -1.27 -0.64 Central -0.99 -0.65 -0.54 -0.16 0.13 -0.23 -0.43 -0.54 -0.36 -0.34 -0.30 -0.74 -0.43 South -0.47 -0.42 0.03 -0.31 -0.33 -0.63 -0.99 -0.85 -0.52 -0.36 -0.19 -0.21 -0.44 Temp Error ( 0 C)
North 3.10 2.88 2.54 2.49 2.44 2.43 2.42 2.49 2.58 2.48 2.89 2.55 2.61 Central 2.38 2.25 1.99 2.18 1.99 2.01 2.07 2.11 2.21 2.52 2.61 2.42 2.23 South 2.31 2.28 1.92 2.13 2.01 2.17 2.19 2.21 2.19 2.70 2.49 2.50 2.26 Wind Speed IO A North 0.79 0.83 0.83 0.87 0.86 0.85 0.81 0.84 0.84 0.82 0.81 0.79 0.83 Central 0.85 0.87 0.88 0.88 0.89 0.86 0.84 0.86 0.87 0.86 0.85 0.84 0.86 South 0.81 0.80 0.85 0.79 0.83 0.83 0.78 0.80 0.82 0.82 0.80 0.82 0.81

Jan Feb Ma r Ap r May Jun Jul A ug Sep Oct Nov Dec Mean RH Bias (%)
North 8.33 9.52 6.63 0.95 -2.42 1.25 2.43 1.60 0.57 0.47 4.47 7.73 3.46 Central 7.43 5.13 4.60 1.65 -1.02 1.52 2.50 1.88 -0.27 -1.40 -0.01 4.35 2.20 South 3.08 -1.19 2.53 2.32 1.26 1.98 2.51 2.62 -0.80 -2.42 -4.45 -1.03 0.53 RH Error (%) North 11.85 13.18 11.61 11.13 11.90 10.04 9.54 9.08 10.26 10.26 11.55 11.61 11.00 Central 12.21 12.43 11.26 10.58 10.72 9.89 9.55 9.54 10.22 10.25 11.42 11.26 10.78 South 11.24 11.76 10.34 8.95 9.30 9.49 9.46 9.61 9.68 9.33 11.63 10.95 10.14 Temp Bias ( 0 C)
North -0.70 -0.82 -0.96 -0.52 -0.25 -0.36 -0.53 -0.49 -0.44 -0.67 -0.76 -0.69 -0.60 Central -0.57 -0.65 -0.79 -0.62 -0.41 -0.68 -0.81 -0.74 -0.49 -0.54 -0.55 -0.52 -0.61 South -0.23 -0.13 -0.52 -0.61 -0.61 -0.94 -0.94 -1.07 -0.65 -0.47 0.04 -0.13 -0.52 Temp Error ( 0 C)
North 2.15 2.07 2.04 1.89 1.86 1.83 1.86 1.80 1.95 1.78 1.99 2.15 1.95 Central 2.12 2.05 2.14 1.95 1.91 1.93 1.93 1.92 2.02 1.77 2.00 2.00 1.98 South 2.18 2.05 2.17 1.83 1.89 1.91 1.88 2.00 1.92 1.68 2.06 1.93 1.96 Wind Speed IO A North 0.82 0.84 0.86 0.88 0.86 0.85 0.85 0.83 0.85 0.85 0.81 0.78 0.84 Central 0.87 0.88 0.90 0.90 0.88 0.87 0.84 0.84 0.87 0.88 0.85 0.85 0.87 South 0.86 0.86 0.85 0.85 0.84 0.82 0.79 0.80 0.83 0.83 0.83 0.82 0.83

Jan Feb Ma r Ap r May Jun Jul A ug Sep Oct Nov Dec Mean RH Bias (%)
North 10.15 7.40 6.01 0.93 -3.76 -0.38 1.38 2.04 -1.66 -1.99 2.96 7.68 2.56 Central 6.94 4.76 4.15 0.42 -2.18 0.17 2.08 2.13 -2.05 -4.13 0.00 5.47 1.48 South 0.00 0.00 0.47 -1.10 -0.37 0.54 1.77 2.89 -3.31 -6.01 -3.66 -0.33 -0.76 RH Error (%) North 13.30 11.21 12.32 11.70 11.65 10.03 9.70 9.57 11.13 12.68 11.53 11.85 11.39 Central 12.77 10.95 11.61 11.18 10.33 9.91 9.49 9.50 10.70 12.69 12.10 12.43 11.14 South 11.18 10.00 9.85 10.17 9.20 9.54 8.90 9.91 10.21 12.12 12.15 12.39 10.47 Temp Bias ( 0 C)
North -1.24 -0.99 -0.63 -0.29 -0.11 -0.10 -0.22 -0.49 -0.34 0.29 -0.85 -1.34 -0.53 Central -0.84 -0.80 -0.64 -0.47 -0.27 -0.36 -0.60 -0.66 -0.32 0.30 -0.54 -0.89 -0.51 South -0.17 -0.27 -0.36 -0.43 -0.46 -0.62 -0.91 -0.98 -0.28 0.53 0.00 -0.03 -0.33 Temp Error ( 0 C)
North 2.31 2.15 2.14 2.02 1.81 1.77 1.91 1.98 2.25 2.57 2.30 2.67 2.16 Central 2.14 2.03 2.15 2.13 1.80 1.81 1.96 1.99 2.16 2.54 2.31 2.45 2.12 South 2.10 1.90 2.00 2.08 1.84 1.81 1.88 2.06 1.94 2.40 2.28 2.48 2.06 Wind Speed IO A North 0.79 0.81 0.83 0.86 0.87 0.85 0.86 0.87 0.84 0.82 0.80 0.82 0.83 Central 0.85 0.88 0.87 0.89 0.90 0.87 0.87 0.86 0.87 0.87 0.86 0.86 0.87 South 0.83 0.83 0.85 0.83 0.85 0.81 0.83 0.82 0.85 0.84 0.82 0.82 0.83 CALMET Model Evaluation Statistics for 2001.

CALMET Model Evaluation Statistics for 2002.

CALMET Model Evaluation Statistics for 2003.E-13 Table 5-5. Overview of Ambient Data Moni toring Networks Covering the CENRAP Domain.

Monitoring Network Chemical Species Measured Sampling Period Data Availability/Source

The Interagency Monitoring of Protected Visual Environments (IMPROVE) Speciated PM25 and PM10 (see species mappings) 1 in 3 days; 24 hr

average http://vista.cira.colostate.e du/improve/Data/IMPRO VE/improve_data.htm

Clean Air Status and Trends Network (CASTNET) Speciated PM25, Ozone (see species mappings) Approximately 1- week average http://www.epa.gov/castne t/data.html

National Atmospheric Deposition Program (NADP) Wet deposition (hydrogen (acidity as pH), sulfate, nitrate, ammonium, chloride, and base cations (such as calcium, magnesium, potassium and sodium)), Mercury 1-week average http://nadp.sws.uiuc.edu/

Air Quality System (AQS) Aka Aerometric Information Retrieval System (AIRS) CO, NO2, O3, SO2, PM25, PM10, Pb Typically hourly

average http://www.epa.gov/air/dat a/

Speciation Trends Network (STN) Speciated PM 24-hour average http://www.epa.gov/ttn/am tic/amticpm.html

Southeastern Aerosol Research and Characterization (SEARCH) (Southeastern US only) 24-hr PM25 (FRM Mass, OC, BC, SO4, NO3, NH4, Elem.); 24-hr PM coarse (SO4, NO3, NH4, elements); Hourly PM2.5 (Mass, SO4, NO3, NH4, EC, TC); Hourly gases (O3, NO, NO2, NOy, HNO3, SO2, CO) Hourly or 24-hour average, depending on parameter. Electric Power Research Institute (EPRI), Southern Company, and other companies.http://www.atmosphericresearch.com

EPA Particulate Matter Supersites (Includes St. Louis in the CENRAP region) Speciated PM25 http://www.epa.gov/ttn/am tic/supersites.html

Photochemical Assessment Monitoring Stations (PAMS) Varies for each of 4 station types. http://www.epa.gov/ttn/am tic/pamsmain.html

National Park Service Gaseous Pollutant Monitoring Network Acid deposition (Dry; SO4, NO3, HNO3, NH4, SO2), O3, meteorological data Hourly http://www2.nature.nps.go v/ard/gas/netdata1.htmE-14 Figure 5-1. CENRAP North, Central, a nd South 6 km Meteorological Domains.E-15 Figure 5-2. CENRAP South Domain.E-16 Figure 5-3. CENRAP Central Domain.E-17 Figure 5-4. CENRAP North Domain.E-18 Figure 5-5. Locations of IMPROVE, CASTNe t, SEARCH, STN and NADP Monitoring Sites in and Near the CENRAP States.E-19 Figure 5-6. Locations of AQS Monitoring Si tes in and Near the CENRAP States.E-20 Figure 5-7. Spatial Mean Relative Humidity (%) over the Central Domain: July 2002.

Figure 5-8. Spatial Mean Surface Temperature (deg C) over the Central Domain: July 2002.E-21 Figure 5-9. Wind Speed Index of Agreement over the Central Domain: July 2002.

Figure 5-10. Standard Deviation in Wind Sp eed (m/s) over the Central Domain: July 2002.E-22 Figure 5-11. Root Mean Square Error in Wi nd Speed (m/s) over the Central Domain: July 2002

.

Figure 5-12. Scalar Mean Wind Speed (m /s) over the Central Domain: July 2002.E-23 Figure 5-13. Vector Mean Wind Speed (m/s) over the Central Domain: July 2002.

Figure 5-14. Normalized Bias in Relative Hu midity (%) over the Central Domain: July 2002.E-24 Figure 5-15. Normalized Error in Relative Hu midity (%) over the Central Domain: July 2002.

Figure 5-16. Relative Humidity (%) at Kenosha, WI: July 2002.E-25 Figure 5-17. Relative Humidity (%) at Topeka, KS: July 2002.

Figure 5-18. Normalized Bias in Surface Temperature (%) over the Central Domain: July 2002.E-26 Figure 5-19. Normalized Error in Surface Temperature over the Central Domain: July 2002.

Figure 5-20. Surface Temperature (deg C) at Kenosha, WI: July 2002.E-27 Figure 5-21. Surface Temperature at Topeka, KS: July 2002.E-28 Figure 5-22. MM5/CALMET Relative Humidity Bias (%) by Month for Three BART Modeling Years (2001, 2003, and 2003).E-29 Figure 5-23. MM5/CALMET Relative Humidity Error (%) by Month for Three BART Modeling Years (2001, 2003, and 2003).E-30 Figure 5-24. MM5/CALMET Temperature Bias (deg C) by Month for Three BART Modeling Years (2001, 2003, and 2003).E-31 Figure 5-25. MM5/CALMET Temperature Erro r (deg C) by Month for Three BART Modeling Years (2001, 2003, and 2003).E-32 Figure 5-26. MM5/CALMET Wind Speed Inde x of Agreement by Month for Three BART Modeling Years (2001, 2003, and 2003).E-33 Figure 5-27. MM5/CALMET Relative Humidity Bias (%) over Three Years in All CENRAP Domains.E-34 Figure 5-28. MM5/CALMET Relative Humidity Error (%) over Three Years in All CENRAP Domains.E-35 Figure 5-29. MM5/CALMET Surface Temperature Bias (deg C) over Three Years in All CENRAP Domains.E-36 Figure 5-30. MM5/CALMET Surface Temperature Error (deg C) over Three Years in All CENRAP Domains.E-37 Figure 5-31. MM5/CALMET Wind Speed Index of Agreement over Three Years in All CENRAP Domains.F-1 

## Appendix F. CALMET Input Control Parameters

Table F-1 Input Groups in the CALMET Control File. Input Group Description Applicable
0 Input and output file names Yes
1 General run control parameters Yes
2 Map Projection and Grid Control Parameters Yes
3 Output Options Yes
4 Meteorological Data Options Yes
5 Wind field Options and Parameters Yes
6 Mixing Height, Temperature and Precipitation Parameters Yes
7 Surface Meteorological Station Parameters Yes
8 Upper Air Meteorological Station Parameters Yes
9 Precipitation Station Parameters Yes

Table F-2 CALMET Model Input Group 0: Input and Output File Names

Parameter Default Value Comments Input GEO.DAT GEO.DAT Input SURF.DAT SURF.DAT Input CLOUD.DAT CLOUD.DAT Input PRECIP.DAT PRECIP.DAT Input MM4.DAT MM4.DAT Input WT.DAT WT.DAT Output CALMET.LST CALMET.LST Output CALMET.DAT CALMET.DAT Output PACOUT.DAT PACOUT.DAT
NUSTA -- 0 Number of upper air stations
NOWSTA -- 0 Number of over water stations Input UP1.DAT UP1.DAT Input UP2.DAT UP2.DAT Input UP3.DAT UP3.DAT Input SEA1.DAT SEA1.DAT Input DIAG.DAT DIAG.DAT Input PROG.DAT PROG.DAT Output TEST.PRT TEST.PRT Output TEST.OUT TEST.OUT Output TEST.KIN TEST.KIN Output TEST.FRD TEST.FRD Output TEST.SLP TEST.SLPF-2 Table F-3 CALMET Model Input Group 1: General Run Control Parameters

Parameter Default Value Comments IBYR - 2001 Starting year
IBMO - 1 Starting month IBDY - 1 Starting day IBHR - 1 Starting hour
IBTZ - 6 Base time zone IRLG - 8760 Length of run
IRTYPE 1 1 Run type (must = 1 to run CALPUFF)
LCALGRD T F Compute CALGRID data fields
ITEST 2 2 Stop run after SETUP to do input QA

Table F-4 CALMET Model Input Group 2: Map Projection and Grid Control Parameters

Parameter Default Value Comments PMAP UTM LCC Map Projection RLATO -- 40N Latitude (dec. degrees) of projection origin RLONO -- 97W Longitude (dec. degrees) of projection origin XLAT1 -- 33N Matching parallel(s) of latitude for projection XLAT2 -- 45N Matching parallel(s) of latitude for projection
DATUM WGS-G WGS-G NX -- 300 Number of X grid cells in meteorological grid NY -- 192 Number of Y grid cells in meteorological grid DGRIDKM -- 6.0 Grid spacing, km XORIGKM -- -1008. Ref. Coordinate of SW corner of grid cell (1,1) YORIGKM -- 0.0 Ref. Coordinate of SW corner of grid cell (1,1) NZ -- 10 No. of vertical layers ZFACE -- 0, 20 40, 80, 160, 320, 640, 1200, 2000, 3000, 4000 Cell face heights in arbitrary vertical grid, mF-3 Table F-5 CALMET Model Input Group 3: Output Options

Parameter Default Value Comments LSAVE T T Disk output option IFORMO 1 1 Type of unformatted output file
LPRINT F F Print met fields
IPRINF 1 1 Print intervals IUVOUT(NZ) NZ*0 NZ*0 Specify layers of u,v wind components to print IWOUT(NZ) NZ*0 NZ*0 Specify layers of w wind component to print ITOUT(NZ) NZ*0 NZ*0 Specify levels of 3-D temperature field to print
LDB F F Print input met data and variables
NN1 1 1 First time step for debug data to be printed
NN2 1 1 Last time step for debug data to be printed
IOUTD 0 0 Control variable for writing test/debug wind fields NZPRN2 1 0 Number of levels starting at surface to print
IPR0 0 0 Print interpolated wind components
IPR1 0 0 Print terrain adjusted surface wind components IPR2 0 0 Print initial divergence fields
IPR3 0 0 Print final wind speed and direction
IPR4 0 0 Print final divergence fields
IPR5 0 0 Print winds after kinematic effects
IPR6 0 0 Print winds after Froude number adjustment
IPR7 0 0 Print winds after slope flows are added
IPR8 0 0 Print final wind field components

Table F-6 CALMET Model Input Group 4: Meteorological Data Options

Parameter Default Value Comments NOOBS 0 2 2 = No surface, over water, or upper air observations; use MM5 for surface, over water, and upper air data NSSTA -- 0 Number of meteorological surface stations NPSTA -- 0 Number of precipitation stations
ICLOUD -- 3 Gridded cloud fields IFORMS 2 2 Formatted surface meteorological data file IFORMP 2 2 Formatted surface precipitation data file
IFORMC 2 2 Formatted cloud data fileF-4 Table F-7 CALMET Model Input Group 5: Wind field Options and Parameters

Parameter Default CENRAP Comments

IWFCOD 1 1 Model selection variable
IFRADJ 1 1 Compute Froude number adjustment effects?
IKINE 0 0 Compute kinematic effects?
IOBR 0 0 Use O’Brien (1970) vertical velocity adjustment?
ISLSOPE 1 1 Compute slope flow effects?
IEXTRP -4 -1 Extrapolate surface wind obs to upper levels?

ICALM 0 0 Extrapolate surface winds even if calm?

BIAS NZ*0 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 Layer-dependent biases weighting aloft

measurements RMIN2 4. -1.0 Minimum vertical extrapolation distance
IPROG 0 14 14 = Yes, use winds from MM5.DAT file as initial guess field [IWFCOD = 1
ISTEPPG 1 1 MM5 output time step
LVARY F T Use varying radius of influence
RMAX1

-- 30. Maximum radius of influence over land in sfc layer
RMAX2 -- 30. Maximum radius of influence over land aloft
RMAX3 -- 50. Maximum radius of influence over water
RMIN 0.1 0.1 Minimum radius of influence used anywhere
TERRAD -- 12. Terrain features radius of influence
R1 -- 1. Weighting of first guess surface field
R2 -- 1. Weighting of first guess aloft field
RPROG -- 0. MM5 windfield weighting parameter DIVLIM 5.E-6 5.E-6 Minimum divergence criterion
NITER 50 50 Number of divergence minimization iterations NSMMTH 2, 4, 4, 4, 4, 4, 4 2, 4, 4, 4, 4, 4, 4 Number of passes through smoothing filter in each layer of CALMET
NITR2 99. 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 Maximum number of stations used in each layer for the interpolation of data to a grid point
CRITFN 1.0 1.0 Critical Froude number ALPHA 0.1 0.1 Kinematic effects parameter FEXTR2 NZ*0.0 NZ*0.0 Scaling factor for extrapolating sfc winds aloft
NBAR 0 0 Number of terrain barriers
IDIOTP1 0 0 Surface temperature computation switch
ISURFT -- 4 Number of sfc met stations to use for temp calcs
IDIOPT2 0 0 Domain-avera ged lapse rate switch
IUPT 0 2 Upper air stations to use for lapse rate calculation
ZUPT 200. 200. Depth through which lapse rate is calculated
IDIOPT3 0 0 Domain-averaged wind component switchF-5 IUPWND -1 -1 Number of aloft stations to use for wind calc ZUPWND 1., 1000. 1., 1000. Bottom and top of layer through which the domain-scale winds are computed
IDIOPT4 0 0 Observed surface wind component switch
IDIOPT5 0 0 Observed aloft wind component switch
LLBREZE F F Use Lake Breeze Module
NBOX 0 0 Number of lake breeze regions
NLB -- 0 Number of stations in the region METBXID(NLB) -- 0 Station ID’s in the region

Table F-8 CALMET Model Input Group 6: Mixing Height, Temperature and Precipitation

Parameter Default Value Comments CONSTB 1.41 1.41 Neutral stability mixing height coefficient CONSTE 0.15 0.15 Convective stability mixing height coefficient CONSTN 2400. 2400. Stable stability maxing height coefficient
CONSTW 0.16 0.16 Over water mixing height coefficient FCORIOL 1.E-4 1.E-4 Absolute value of Coriolis parameter IAVEZI 1 1 Conduct spatial averaging? Yes = 1 MNMDAV 1 10 Maximum search radius in averaging process HAFANG 30. 30. Half-angle of upwind looking cone for averaging ILEVZI 1 1 Layers of wind use in upwind averaging DPTMIN 0.001 0.001 Minimum potential temperature lapse rate in the stable layer above the current convective mixing ht DZZI 200. 200. Depth of layer a bove current conv. mixing height through which lapse rate is computed ZIMIN 50. 50. Minimum overland mixing height ZIMAX 3000. 3000. Maximum overland mixing height ZIMINW 50. 50. Minimum over water mixing height ZIMAXW 3000. 3000. Maximum over water mixing height ITPROG 0 2 3D temperature from observations or from MM5? IRAD 1 1 Type of interpolation; 1 = 1/R TRADKM 500. 36. Temperature interpolation radius of influence NUMTS 5 5 Max number of stations for temp interpolation IAVET 1 1 Spatially average temperatures? 1 = yes TGDEFB -.0098 -.0098 Temp gradient below mixing height over water TGDEFA -.0045 -.0045 Temp gradient above mixing height over water JWAT1 -- 55 Beginning land use categories over water JWAT2 -- 55 Ending land use categories for water NFLAGP 2 2 Precipitation interpolation flag; 2 = 1/R-squared

SIGMAP 100. 50. Radius of influence for precipitation interpolation CUTP 0.01 0.01 Minimum precipitation rate cutoff (mm/hr)G-1 

## Appendix G. CALPUFF Input Control Parameters

Table G-1 Input Groups in the CALPUFF Control File Input Group Description Applicable 0 Input and output file names Yes 1 General run control parameters Yes 2 Technical options Yes 3 Species list Yes 4 Grid control parameters Yes 5 Output options Yes 6 Sub grid scale complex terrain inputs Yes 7 Dry deposition parameters for gases Yes 8 Dry deposition parameters for particles Yes 9 Miscellaneous dry deposition for parameters Yes 10 Wet deposition parameters Yes 11 Chemistry parameters Yes 12 Diffusion and computational parameters Yes 13 Point source parameters Yes 14 Area source parameters No 15 Line source parameters No 16 Volume source parameters No 17 Discrete receptor information Yes

Table G-2 CALPUFF Model Input Group 0: Input and Output File Names
Parameter Default Value Comments

METDAT CALMET.DAT Not used Input file name PUFLST CALPUFF.LST Varies with facility CALPUFF output file name CONDAT CONC.DAT Varies with facility Concentration output file

name DFDAT DFLX.DAT Varies with facility Dry flux output file name WFDAT WFLX.DAT Varies with facility Wet flux output file name VISDAT VISB.DAT Varies with facility Visibility output file name OZDAT OZONE.DAT Varies with year Ozone input file name
LCFILES - T File names converted to lower

case NMETDAT 1 12 Number of CALMET.DAT files for run
CALMET.DAT - METDAT=/location of CALMET.DAT files 12 entries one for each monthG-2 Table G-3 CALPUFF Model Input Group 1: General Run Control Parameters
Parameter Default Value Comments METRUN 0 0 All model periods in met file(s) will be run
IBYR

- See note 1 below Starting year

IBMO - 1 Starting month IBDY - 1 Starting day IBHR - 1 Starting hour XBTZ - 0 Time zone for met files (0 = GMT)
IRLG

- See note 2 below Length of run

NSPEC 5 10 Number of MESOPUFF II chemical species
NSE 3 See note 3 below Number of chemical species to be emitted

ITEST 2 2 Program is executed after SETUP phase MRESTART 0 0 Do not read or write a restart file during run NRESPD 0 0 File written only at last period METFM 1 1 CALMET binary file (CALMET.MET) AVET 60 60 Averaging time in minutes PGTIME 60 60 PG Averaging time in minutes Note 1: Enter the year being modeled (i.e. 2001, 2002, or 2003) Note 2: Enter 8760 for the years 2001 and 2002 but enter 8748 for the year 2003 Note 3: Enter 6 for the no-obs run and 7 for the refined run

Table G-4 CALPUFF Model Inpu t Group 2: Technical Options
Parameter Default Value Comments

MGAUSS 1 1 Gaussian distribution used in near field MCTADJ 3 3 Partial plume path terrain adjustment MCTSG 0 0 Sub-grid-scale complex terrain not modeled MSLUG 0 0 Near-field puffs not modeled as elongated MTRANS 1 1 Transitional plume rise modeled MTIP 1 1 Stack tip downwash used MSHEAR 0 0 (0, 1) Vertical wind shear (not modeled, modeled) MSPLIT 0 0 Puffs are not split MCHEM 1 1 MESOPUFF II chemical parameterization scheme MAQCHEM 0 0 Aqueous phase transformation not modeled
MWET 1 1 Wet removal modeled MDRY 1 1 Dry deposition modeled MDISP 3 2 AERMOD dispersion coefficients
MTURBVW 3 3 Use both σ v and σ w from PROFILE.DAT to compute σ y and σ z (n/a) MDISP2 3 2 AERMOD dispersion coefficients

MROUGH 0 0 PG σ y and σ z not adjusted for roughness
MPARTL 1 1 No partial plume penetration of elevated inversionG-3 MTINV 0 0 Strength of temperature inversion computed from default gradients
MPDF 0 0 PDF not used for dispersion under convective conditions MSGTIBL 0 0 Sub-grid TIBL module not used for shoreline MBCON 0 0 Boundary concentration conditions not modeled
MFOG 0 0 Do not configure for FOG model output MREG 1 1 Technical options must conform to USEPA Long Range Transport (LRT) guidance

Table G-5 CALPUFF Model Input Group 3: Species List-Chemistry Options

CSPEC Modeled 1 Emitted 2 Dry Deposition 3 Output Group Number

SO 2 1 1 1 0
SO 4 -2 1 0 2 0
NOx 1 1 1 0
HNO 3 1 0 1 0
NO 3 - 1 0 2 0
NH 3 0 0 1 0
PM 10 4 1 1 2 0
PMF 4 1 1 2 0
EC 5 1 1 2 0
SOA 5 1 1 2 0 Note 1: 0 = No, 1 = Yes Note 2: 0 = No, 1 = Yes (Depends on if species is being modeled or not) Note 3: 0 = none, 1 = computed gas, 2 = computed particle, 3 = user specified Note 4: Only PMF will be modeled and emitted in the no-obs run; however, both PM 10 and PMF will be modeled and emitted in the refined analysis Note 5: EC and SOA will not be modeled nor will it be emitted during the no-obs and the refined runsG-4 Table G-6 CALPUFF Model Input Group 4: Map Projection and Grid Control Parameters
Parameter Default Value Comments PMAP UTM LCC Map Projection FEAST 0.0 0.000 False Easting FNORTH 0.0 0.000 False Northing
RLATO None 40N Latitude and L ongitude of projection origin RLONO None 97W Latitude and Long itude of projection origin
XLAT1 None 33N Matching parallel of latitude for map projection
XLAT2 None 45N Matching parallel of latitude for map projection
DATUM WGS-84 WGS-G Datum region for output coordinates
NX None 366 Number of X grid cells in meteorological grid
NY None 234 Number of Y grid cells in meteorological grid
NZ None 10 Number of vertical layers in meteorological grid
DGRIDKM None 6 Grid spacing (km)
ZFACE None 0, 20 40, 80, 160, 320, 640, 1200, 2000, 3000, 4000 Cell face heights in meteorological grid (m)

XORIGKM None -1008 Reference X coordinate for SW corner of grid cell (1,1) of meteorological grid (km)
YORIGKM None -864 Reference Y coordinate for SW corner of grid cell (1,1) of meteorological grid (km)
IBCOMP None 1 X index of lower left corner of the computational grid
JBCOMP None 1 Y index of lower left corner of the computational grids
IECOMP None 366 X index of the upper right corner of the computational grid
JECOMP None 234 Y index of the upper right corner of the computational grid
LSAMP T F Sampling grid is not used
IBSAMP None 1 X index of lower left corner of the sampling grid
JBSAMP None 1 Y index of lower left corner of the sampling grid
IESAMP None 366 X index of upper right corner of the sampling grid
JESAMP None 234 Y index of upper right corner of the sampling grid
MESHDN 1 1 Nesting factor of the sampling gridG-5 Table G-7 CALPUFF Model Input Group 5: Output Options
Parameter Default Value Comments

ICON 1 1 Output file CONC.DAT containing concentrations is created
IDRY 1 1 Output file DFLX.DAT containing dry fluxes is created
IWET 1 1 Output file WFLX.DAT containing wet fluxes is created
IVIS 1 1 Output file containing relati ve humidity data is created LCOMPRS T T Perform data compression in output file IMFLX 0 0 Do not calculate mass fluxes across specific boundaries IMBAL 0 0 Mass balances for each species not reported hourly ICPRT 0 1 Print concentration fields to the output list file IDPRT 0 0 Do not print dry flux fields to the output list file IWPRT 0 0 Do not print wet flux fields to the output list file ICFRQ 1 1 Concentration fields are printed to output list file every hour (hr) IDFRQ 1 1 Dry flux fields are printed to output list file every 1 hour IWFRQ 1 1 Wet flux fields are printed to output list file every 1 hour IPRTU 1 3 Units for line printer output are in g/m3 for concentration and g/m 2 /s for deposition IMESG 2 2 Messages tracking the progress of run written to

screen LDEBUG F F Logical value for debug output IPFDEB 1 1 First puff to track NPFDEB 1 1 Number of puffs to track NN1 1 1 Meteorological period to start output NN2 10 10 Meteorological period to end outputG-6 Table G-8 CALPUFF Model Input Group 6: Sub-Grid Scale Complex Terrain Input
Parameter Default Value Comments

NHILL 0 0 Number of terrain features NCTREC 0 0 Number of special complex terrain receptors MHILL - 2 Input terrain and receptor data for CTSG hills input in CTDM format XHILL2M 1 1 Conversion factor for changing horizontal dimensions to meters
ZHILL2M 1 1 Conversion factor for changing vertical dimensions

to meters XCTDMKM None 0.0 E+00 X origin of CTDM system relative to CALPUFF coordinate system (km)
YCTDMKM None 0.0 E+00 Y origin of CTDM system relative to CALPUFF coordinate system (km)

Table G-9 CALPUFF Model Input Group 7: Dry Deposition Parameters for Gases Species Default Value Comments 0.1509 0.1509 Diffusivity 1000. 1000. Alpha star 8.0 8.0 Reactivity

0.0 0.0 Mesophyll resistance SO 2

0.04 0.04 Henry’s Law coefficient 0.1656 0.1656 Diffusivity

1.0 1.0 Alpha star 8.0 8.0 Reactivity

5.0 5.0 Mesophyll resistance NO
x

3.5 3.5 Henry’s Law coefficient 0.1628 0.1628 Diffusivity

1.0 1.0 Alpha star 18.0 18.0 Reactivity

0.0 0.0 Mesophyll resistance 8.0E-8 8.0E-8 Henry’s Law coefficient HNO 3

0.000359 0.000359 Henry’s Law coefficientG-7 Table G-10 CALPUFF Model I nput Group 8: Dry Deposition Parameters for Particles Species Default Value Comments

SO 4 -2 0.48 0.48 Geometric mass mean diameter of SO 4 -2 (μm)
NO 3 - 2.0 0.48 Geometric mass mean diameter of NO 3 - (μm)
PM 10 2.0 6.0 Geometric mass mean diameter of PMC (μm)
PMF 2.0 0.48 Geometric mass mean diameter of PMF (μm) EC 2.0 0.48 Geometric mass mean diameter of EC (μm)
SOA 0.48 0.48 Geometric mass mean diameter of SOA (μm) (Geometric Standard Deviation for all species assumed to be 2.0 μm).

Table G-11 CALPUFF Model I nput Group 9: Miscellaneous Dry Deposition Parameters
Parameter Default Value Comments RCUTR 30 30 Reference cuticle resistance (s/cm) RGR 10 10 Reference ground resistance (s/cm) REACTR 8 8 Reference pollutant reactivity NINT 9 9 Number of particle size intervals for effective particle deposition velocity IVEG 1 1 Vegetation in non-irrigated areas is active and unstressed

Table G-12 CALPUFF Mode l Input Group 10: Wet Deposition Parameters Species Default Value Comments 3.21E-05 3.21E-05 Scavenging coefficient for liquid precipitation (s -1 ) SO 2 0.0 0.0 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) SO 4 -2 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 ) 6.0E-05 6.0E-05 Scavenging coefficient for liquid precipitation (s -1 ) HNO 3 0.0 0.0 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) NO 3 - 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 ) 8.0E-05 8.0E-05 Scavenging coefficient for liquid precipitation (s -1 ) NH 3 0.0 0.0 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) PM 10 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) PMF 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) EC 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 ) 1.0E-04 1.0E-04 Scavenging coefficient for liquid precipitation (s -1 ) OC 3.0E-05 3.0E-05 Scavenging coefficient for frozen precipitation (s -1 )G-8 Table G-13 CALPUFF Model Input Group 11: Chemistry Parameters
Parameter Default Value Comments

MOZ 1 1 Read ozone background concentrations from ozone.dat file (measured values).
BCKO 3 12*80 12*40 Background ozone concentration (ppb)
BCKNH 3 12*10 12*3 Background ammonia concentration (ppb)
RNITE1 0.2 0.2 Nighttime NO 2 loss rate in percent/hour
RNITE2 2 2 Nighttime NO X loss rate in percent/hour
RNITE3 2 2 Nighttime HNO 3 loss rate in percent/hour
MH 2 0 2 1 1 Background H 2 O 2 concentrations (Aqueous phase transformations not modeled )

BCKH 2 0 2 1 1 Background monthly H 2 O 2 concentrations (Aqueous phase transformations not modeled )G-9 Table G-14 CALPUFF Model Input Group 12: Dispersion/Computational Parameters
Parameter Default Value Comments SYDEP 550 550 Horizontal size of a puff in meters beyond which the time dependant dispersion equation of Heffter (1965) is used
MHFTSZ 0 0 Do not use Heffter formulas for sigma z
JSUP 5 5 Stability class used to determine dispersion rates for puffs above boundary layer CONK1 0.01 0.01 Vertical dispersion constant for stable conditions CONK2 0.1 0.1 Vertical disp ersion constant for neutral/stable conditions TBD 0.5 0.5 Use ISC transition point for determining the transition point between the Schulman-Scire to Huber-Snyder Building Downwash scheme
IURB1 10 10 Lower range of land use categories for which urban dispersion is assumed
IURB2 19 19 Upper range of land use categories for which urban dispersion is assumed
ILANDUIN 20 * Land use category for modeling domain
XLAIIN 3.0 * Leaf area index for modeling domain
ZOIN -0.25 * Roughness length in meters for modeling domain
ELEVIN 0.0 * Elevation above sea level
XLATIN -999 - North latitude of station in degrees
XLONIN -999 - South latitude of station in degrees
ANEMHT 10 10 Anemometer height in meters
ISIGMAV 1 1 Sigma-v is read for lateral turbulence data
IMIXCTDM 0 0 Predicted mixing heights are used
XMXLEN 1 1 Maximum length of emitted slug in meteorological grid units
XSAMLEN 1 10 Maximum travel distance of slug or puff in meteorological grid units during one sampling unit
MXNEW 99 60 Maximum number of puffs or slugs released from one source during one time step
MXSAM 99 60 Maximum number of sampling steps during one time step for a puff or slug
NCOUNT 2 2 Number of iterations used when computing the transport wind for a sampling step that includes transitional plume rise
SYMIN 1 1 Minimum sigma y in meters for a new puff or slug
SZMIN 1 1 Minimum sigma z in meters for a new puff or slug
SVMIN .50 .50 Minimum lateral turbulence velocities (m/s) SWMIN 0.20, 0.12, 0.08, 0.06, 0.20, 0.12, 0.08, 0.06, Minimum vertical turbulence velocities (m/s)G-10 0.03, 0.016 0.03, 0.016
WSCALM 0.5 0.5 Minimum non-calm wind speeds (m/s) XMAXZI 3000. 3000. Maximum mixing height (m) XMINZI 50. 20. Minimum mixing height (m)
SL2PF 10. 10. Maximum Sy/puff length PLXO 0.07, 0.07, 0.10, 0.15, 0.35, 0.55 0.07, 0.07, 0.10, 0.15, 0.35, 0.55 Wind speed power-law exponents

WSCAT 1.54, 3.09, 5.14, 8.23, 10.80 1.54, 3.09, 5.14, 8.23, 10.80 Upper bounds of 1 st 5 wind speed classes

PGGO 0.020, 0.035 0.020, 0. 035 Potential temp gradients PG E & F (deg/km) CDIV 0.01 0.01 Divergence criterion for dw/dz (1/s)
PPC 0.5, 0.5, 0.5, 0.5, 0.35, 0.35 0.5, 0.5, 0.5, 0.5, 0.35, 0.35 Plume path coefficients (only if MCTADJ=3)

NSPLIT 3 3 Number of puffs when puffs split
IRESPLIT - 1900 Hour(s) when puff is eligible to split ZISPLIT 100 100 Previous hour’s minimum mixing height, m
ROLDMAX 0.25 0.25 Previous Max mixing height/current mixing height ratio, must be less than this value to allow puff to split
NSPLITH 5 5 Number of puffs resulting from a split
SYSPLITH 1.0 1.0 Minimum sigma-y of puff before it may split
SHSPLITH 2.0 2.0 Minimum puff elongation rate from wind shear before puff may split CNSPLITH 1.0E-07 1.0E-07 Minimum species concentration before a puff may split EPSSLUG 1.0E-04 1.0E-04 Criterion for SLUG sampling EPSAREA 1.0E-06 1.0E-06 Criterion for area source integration
DSRISE 1.0 1.0 Trajectory step length for numerical rise algorithm Note: Values indicated by an asterisk (*) were allowed to vary spatially across the domain and were obtained from CALMET

Table G-15 CALPUFF Mode l Input Group 13: Point Source Parameters
Parameter Default Value Comments NPT1 None Varies by scenario Number of point sources with stack parameters

IPTU 1 1 Units for point source emission rates are g/s
NSPT1 0 0 Number of source-species combinations with variable emissions scaling factors
NPT2 None 0 Number of point sources with variable emission parameters provided in external file
MISC None Point source parameters and emission data Point source inputs include stack height (H), stack diameter (d), exit temperature (T), exit velocity (v) emissions by species, and coordinate of stack (LCC)G-11 Table G-16 CALPUFF Model Input Group 17: Discrete Receptor Information
Parameter Default Value Comments NREC None 427 Number of discrete receptors

Please note that ADEQ will not be modeli ng area, line and volume sources which are input groups 14, 15, and 16 respectively.H-1 

## Appendix H. POSTUTIL Input Control Parameters

Table H-1 Input Groups in the PO STUTIL Processor Control File
Sub Group Description Applicable 0a Input and output file names Yes 1 NMET – Number of CALMET data files (365) Yes
2 NFILES – Number of CALPUFF data files Yes

Table H-2 POSTUTIL Processor Input Group 1: General Run Control Parameters
Parameter Default Value Comments

ISYR None See note 1 below Starting year

ISMO None 1 Starting month
ISDY None 1 Starting day
ISHR None 0 Starting hour
NPER None See note 2 below Number of periods to process

NSPECINP None See note 3 below Number of CALPUFF species to process

NSPECOUT None See note 3 below Number of species to output

NSPECCMP None 0 Number of species to derive
MDUPLCT None 1 Stop run if duplicate name
NSCALED None 0 Number of CALPUFF files to ‘scale’
MNITRATE None 1 Re-compute the HNO 3 /NO 3 partition for CALPUFF modeled concentrations? 1 = yes for all sources combined
BCKNH3 10. 3. Default NH 3 concentration (ppb) for HNO 3 /NO 3 partitioning Note 1: Enter the modeled year for the CALPUFF run Note 2: Enter 8760 for years 2001 and 2002, but enter 8748 for the year 2003 Note 3: Enter 6 for the no-obs run and 7 for the refined runH-2 Table H-3 POSTUTIL Processor Input Group 2: Species Processing Information
Parameter Default Value Comments

ASPECI None SO 2 , SO 4 , NO
x , HNO 3 , NO 3 , PM 10 , PMF See Note 1 Below Species to post-process

ASPECO None SO 2 , SO 4 , NO
x , HNO 3 , NO 3 , PM 10 , PMF See Note 4 Below Species to output

CSPECCMP None CSPECCMP = N SO 2 = 0.0 SO 4 = 0.291667 NO = 0.466667 NO 2 = 0.304348 HNO 3 = 0.222222 NO 3 = 0.451613
PM 10 = 0.0 Nitrogen species to be computed by scaling and summing one or more of the processed input species using the scaling factors for each of the NSPECINP input species

CSPECCMP None CSPECCMP = S SO 2 = 0.50 SO 4 = 0.333333 NO = 0.0 NO 2 = 0.0 HNO 3 = 0.0 NO 3 = 0.0
PM 10 = 0.0 Sulfur species to be computed by scaling and summing one or more of the processed input species using the scaling factors for each of the NSPECINP input species

MODDAT None A (Default=1.0) SO 2 = 1.1 SO 4 = 1.5 HNO3 = 0.8 NO 3 = 0.1

B (Default=0.0) SO 2 = 0.0 SO 4 = 0.0 HNO 3 = 0.0 NO 3 = 0.0 Each species in NSCALED CALPUFF data files may be scaled before processing (e.g., to change the emission rate for all sources modeled in the run that produced a data file). For each scaled species the scaling factors are A and B where x' = Ax +
B.

Note 4: In the no-obs run just enter PMF, but in the refined run enter PM 10 and PMFI-1 

## Appendix I. CALPOST Input Control Parameters

Table I-1Input Groups in the CALPOST Processor Control File Group Description Applicable 0 Input and output file names Yes
1 General Run Control Parameters Yes 2 Visibility Parameters Yes 3 Output Options Yes

Table I-2 CALPOST Processor Input Gr oup 1: General Run Control Parameters
Parameter Default Value Comments

METRUN 0 1 1 = Run all met pe riods in CALPUFF data file ISYR None 2001, 2002, 2003 Starting year

ISMO None 1 Starting month
ISDY None 1 Starting day
ISHR None 0 Starting hour
NPER None See note 1 below Number of periods to process

NREP 1 1 Process every hour of data? Yes = 1 ASPEC None VISIB Process species for visibility
ILAYER 1 1 Layer/deposition code; 1 for CALPUFF concentrations
A 0.0 0.0 Scaling factor, slope
B 0.0 0.0 Scaling factor, intercept
LBACK F F Add hourly background concentrations or fluxes
MSOURCE 0 0 Process only total reported contribution
LG F F Process gridded receptors
LD F T Process discrete receptors
LCT F F Process complex terrain receptors
LDRING F F Report receptor ring results
NDRECP -1 See note 2 below To select the Class I area’s receptors enter *1 after the number of receptors otherwise enter *0 IBGRID -1 -1 X index of LL corner of receptor grid JBGRID -1 -1 Y index of LL corner of receptor grid IEGRID -1 -1 X index of UR corner of receptor grid JEGRID -1 -1 X index of UR corner of receptor grid
NGONOFF 0 0 Number of gridded receptor rows
NGXRECP 1 0 Exclude specific gridded receptors, Yes = 0 Note 1: Enter 8760 for the years 2001 and 2002, but enter 8748 for the year 2003. Note 2: CALPOST is to be run for each Class I area assessed. 3⁄4 The following are the number of receptors for each Class I area being assessed: 1. Caney Creek = 80 2. Hercules-Glade = 47 3. Mingo Wilderness = 80 4. Sipsey = 148 5. Upper Buffalo = 72I-2 Table I-3 CALPOST Processor Input Gr oup 2: Species Processing Information
Parameter Default Value Comments RHMAX 98 95 Maximum RH (%) used in particle growth curve
LVSO4 T T Compute light extinction for sulfate?
LVNO3 T T Compute light extinction for nitrate?
LVOC T T Compute light extinction for organic carbon?
LVMPC T T Compute light extinction for coarse particles?
LVMPF T T Compute light extinction for fine particles?
LVEC T T Compute light extinction for elemental carbon?
LVBK T T Include background in extinction calculation? SPECPMC PMC PMC Coarse particulate species SPECPMF PMF PMF Fine particulate species EEPMC 0.6 0.6 Extinction efficiency for coarse particulates EEPMF 1.0 1.0 Extinction efficiency for fine particulates EEPMCBK 0.6 0.6 Extinction efficiency for coarse part. background EESO4 3.0 3.0 Extinction efficiency for ammonium sulfate EENO3 3.0 3.0 Extinction efficiency for ammonium nitrate EEOC 4.0 4.0 Extinction efficiency for organic carbon EESOIL 1.0 1.0 Extinction efficiency for soil EEEC 10.0 10.0 Extinction efficiency for elemental carbon
MVISBK 2 6 Method 6 for background light extinction: Compute extinction from speciated PM measurements. FLAG RH adjustment factor applied to observed & modeled sulfate and nitrate
BEXTBTBK -- 12 Background extinction for MVISBK=1 (1/Mm)
RHFRAC -- 10 Percentage of particles affected by RH RHFAC 12*value Depends on Class I Area Extinction coefficients for modeled and background hygroscopic species computed using EPA (2003) monthly RH adjustment factors
BKSEC 0.02 0.02 Eastern background elemental carbon β ext BKSO4 0.23 0.23 Eastern background sulfate β ext BKNO3 0.10 0.10 Eastern background nitrate β ext BKPMC 3.00 3.00 Eastern background PMC β ext BKSOC 1.40 1.40 Easter background organic carbon β ext BKSSOIL 0.50 0.50 Eastern background soil β ext BKSEC 0.02 0.02 Eastern background elem. β ext BEXTRAY 10.0 10.0 Extinction due to Rayleigh scattering (1/Mm)I-3 Table I-4 CALPOST Processor Input Group 3: Output Options
Parameter Default Value Comments

LDOC F F Print documentation image
IPRTU 1 3 Print output units (μg/m 3 ) for concentrations and (μg/m 2 /sec) for deposition
L1HR T F Report 1 hr averaging times
L3HR T F Report 3 hr averaging times
L24HR T T Report 24 hr averaging times
LRUNL T F Report run-length (annual) averaging times
LT50 T F Top 50 table
LTOPN F F Top ‘N’ table
NTOP 4 4 Number of ‘Top-N’ values at each receptor ITOP 1,2,3,4 1,2,3,4 Ranks of ‘Top-N’ values at each receptor
LEXCD F F Threshold exceedances counts THRESH1 -1.0 -1.0 Averaging time threshold for 1 hr averages THRESH3 -1.0 -1.0 Averaging time threshold for 3 hr averages THRESH24 -1.0 -1.0 Averaging time threshold for 24 hr averages THRESHN -1.0 -1.0 Averaging ti me threshold for NAVG-hr averages
NDAY 0 0 Accumulation period, days
NCOUNT 1 1 Number of exceedances allowed LECHO F F Echo option
LTIME F F Time series option
LPLT F F Plot file option
LGRD F F Use grid format instead of DATA format
LDEBUG F F Output information for debugging?J-1 

## Appendix J. Chapter 3 of the CENRAP BART Modeling Guidelines (Tesche, et al, 2005)

3.0 CALPUFF FORMULATION AND IMPLEMENTATION

The RHR relates visibility attenuation to extinction coefficient (b ext ) which is a measure of light scattering and absorption due to atmospheri c constituents. Values for b ext are estimated using an empirically derived equation which relates the extinction coefficient to relative humidity and the follo wing components of particulate matter mass: (a) sulfates (SO 4 ); (b) nitrates (NO 3 ); (c) organic carbon (OC); (d) elemental carbon (EC); (f) particulate matter (IP) (“crustal mate rial”); and (g) coarse mass (CM) (i.e., PM 10 – PM 2.5 ). The BART guidance requires the us e of modeled concentrations of these components, together with a “humidity corre ction factor”, to es timate values for b ext on all days within a three year period. These estimates, when compared with naturally occurring background extinction, are used to determine whether a source is causing or contributing to visibility impairment and also to measure the effectiveness of emissions controls on the source aimed at mitigating such effects. EPA notes that secondary particulate matter constitutes an important fraction of PM 2.5 and that the modeling requirements for secondary and primary particul ate matter differ in their need to consider atmospheric chemistry and in the degree of spatial resolution needed for the modeling (EPA, 2001, pg 22).

This chapter introduces the formulation of the CALPUFF modeling system. We summarize the model capabilities as describe d in the user’s manuals (Scire et al., 2000a,b) and discuss the capabilities and limitations of the model. Equipped with this information, states and source operators can id entify those situations for which screening and/or source-specific applications of CALPUFF are appropriate.

In most cases, we expect that application of the CALPUFF system will be sufficient to meet the BART Rule requirements. For that subset of conditions requiring advanced methods, Chapter 5 provides details on full-science alternative models and available data bases for BART modeling. Such conditions might include a situation where the default modeling shows that a sour ce just barely causes or contributes to visibility degradation or in negotiations over the final BART determination that weighs technical and economic feasibility against expected air quality benefits. In both situations, a more accurate estimate of a s ource’s impacts may be very important to source operators.

3.1 Original Model Development

The CALPUFF modeling system was origin ally developed as a component of a three-part modeling system sponsored by th e California Air Resources Board (ARB) in the mid-1980s. The ARB sought to develop a new puff-based model, a new grid-based model and an improved meteorological proce ssor that would support application of the two. CALGRID was the urban-scale photochemical grid model resulting from theJ-2 project (Yamartino et al., 1992) comparable in science a nd capabilities to the Urban Airshed Model (UAM-IV) (Scheffe and Morris, 1993). The model formulation was aimed at overcoming the deficiencies in EPA’ s steady-state Gaussian plume models that were routinely used in California for inert and linearly reactive ma terials (principally SO 2 ) from elevated point sources. Thus, th e CALGRID model was designed to treat the complexities of urban-scale photochemical pr ocesses while CALPUFF was formulated to treat the non-steady state transport, diffusion, linear reaction, and deposition of primary pollutants from point sources. CALPUFF wa s not designed to address photochemical oxidants or and secondary aerosol formation production processes in a scientifically rigorous manner.

In recent years, CALPUFF and its meteorological pre-processor (CALMET) have been used in a range of regu latory modeling studies to addr ess point source issues that include complexities posed by complex terrain, large source-receptor distances, parameterized chemical transformation and deposition, and issues related to Class I visibility impacts. These applications ar e more complex than the California ARB’s nonsteady-state, linear chemistry fo rmulation of the mid-1980s.

The CALPUFF modeling system has been adopted by the EPA as a guideline model for source-receptor distances greater than 50 km, and for use on a case-by-case basis in complex flow situati ons for shorter distances. It was recommended for Class I impact assessments by the FLM Workgr oup (FLAG, 2000) and the Interagency Workgroup on Air Quality Modeling (IWAQM) (E PA, 1998). As directed in the BART guidance, CALPUFF is the primary modeling system for screening and source-specific BART applications in the CENRAP re gion. Thus, examination of the model’s formulation provides the context for assessing the extent to which it suitable for simulating the various physical processes an d gas-phase, aerosol , and aqueous-phase chemical processes that influences visibility.

3.2 CALPUFF Model Formulation

The CALPUFF user’s guide (Scire et al ., 2000a) depicts the modeling system as shown in Figure 3-1. CALMET is a diag nostic/interpolation model that provides meteorological inputs to CA LPUFF. These fields in clude hourly-averaged threedimensional wind and temperatur e fields and two-dimensiona l fields of mixing heights and other meteorological parameters. CALMET uses routine surface and aloft meteorological observations and/or three-dimensional output from prognostic numerical models such as MM5 (Grell et al., 1995) or RUC (Benjamin et al., 2004) to construct the meteorological inputs. Other inputs to the air quality program include emissions information, receptor locations, ancillary geophysical information, and estimated concentrations of ambient pollutants that ar e entrained by the modeled puffs as each is carried downwind. Tabl es 3-1 and 3-2 summarize the key features of the CALMET/CALPUFF models as desc ribed in the user’s guides.

Two post-processor routines are included to facilitate cumulative source impacts (POSTUTIL) and estimates of light extinc tion and visibility attenuation at Class IJ-3 receptors of interest (CALPOST). In par ticular, CALPOST contains several options for computing change in extinction and decivi ews for visibility assessments while the POSTUTIL postprocessor includes options fo r summing contributions of individual sources or groups of sources to assess cumula tive impacts. POSTUTIL also contains an empirical nitric acid-nitrate chemical e quilibrium module to estimate the cumulative effects of ammonia consumption by b ackground sources once the simulation is completed.

3.2.1 Model Concept and Governing Equations

The starting point for the CALPUFF development was the choice of the fundamental reference system of which ther e are two: Eulerian and Lagrangian. Consistent with the original ARB design crite ria, the Lagrangian (moving puff) reference system was chosen for CALPUFF. In the Eule rian approach, the behavior of pollutants is described relative to a fixed coordinate sy stem. The Lagrangian reference frame, in contrast, relates the behavior of pollutants relative to a coordi nate system that moves with the average wind. These two approaches yield different mathematical relationships for pollutant concentrations that are equally valid . The choice of which approach to adopt depends upon the specific design goals of the modeling system.

The advantages and drawbacks of each a pproach are thoroughly discussed in the literature (Tesche, 1983; Seinfeld and Pandis, 1998; Jacobson, 1999; Russell and Dennis, 2000). One of the criticisms of early Eulerian grid models was their ‘over-dilution’ of point source emissions into the fixed grid cells; but for the past twenty years, this limitation has been overcome through with the development of sub-grid-scale, plume-ingrid algorithms (Seigneur, et al., 1981; Godowitch, 2004; Karamcha ndani et al., 2005; Emery and Yarwood, 2005) and the use of multi- scale nested grids (Russell and Dennis, 2000). While the Lagrangian approach is conceptually simple, flexible, and computationally inexpensive, the governing e quations are not directly applicable to situations involving non- linear chemical reactions (Seinfel d and Pandis, (1998) and it is awkward to handle a large numbe r of sources realistically.

3.2.2 Transport and Dispersion

Adopting the Lagrangian concept, CALPU FF simulates the transport, dispersion, linear chemical transf ormation, and deposition of indivi dual puffs carried downwind by the three-dimensional fields generated by CALMET. The model’s implementation follows puffs from the near source region (a few tens of meters) to hundreds of kilometers downwind. Its puff-based formula tion, in conjunction with three-dimensional hourly meteorological data, allow CALPUFF to simulate the effects of time- and spacevarying meteorological conditions on pollutant s emitted from a variety of source types. The major features and options of the CALPUFF model are summarized below:

Building Downwash: The Huber-Snyder and Schulman-Scire downwash models are both incorporated into CALPUFF. An option is provided to use eith er model for all stacks, or make the choice on a stack-by-stack and wi nd sector-by-wind sector basis. BothJ-4 algorithms have been implemented in such a way as to allow the use of wind direction specific building dimensions. The PRIME bu ilding downwash model (Schulman et al., 2000) is also included in CALPUFF as an option.

Dispersion Coefficients: Turbulent dispersion in CALPUF F is treated with the K-theory (flux-gradient) closure scheme, defined for a Lagrangian frame of reference. Several options are provided in CALPUFF for the computation of dispersion coefficients, including the use of turbulence measurements ( σ v and σ w ), the use of similarity theory to
estimate σ v and σ w from modeled surface heat and mo mentum fluxes, or the use of Pasquill-Gifford (PG) or McElroy-Pooler (M P) dispersion coefficients, or dispersion equations based on the CTDM. Options ar e provided to apply an averaging time correction or surface roughness length adjustments to the PG coefficients. Recently, the EPA AERMOD dispersion parameters have b een included in CALPUFF and are used regularly.

Puff Sampling Functions: Puff sampling routines are included in CALPUFF to address computational difficulties encountered when a pplying a puff model to near-field releases. For near-field applications during rapi dly-varying meteorolog ical conditions, an elongated puff (slug) sampling function may be used. An integrated puff approach may be used during less demanding conditions. Both techniques repr oduce continuous plume results under the appropriat e steady state conditions.

Wind Shear Effects : A key underpinning of the Lagrangi an concept is that the modeled puffs retain their identity over the time- and spatial-scale asso ciated with the effects the model is attempting to predict (i.e., visibi lity impairment at 200 km or beyond) While discrete puffs emitted from a source retain thei r physical integrity for a period of time, at some point the action of horizontal and vertic al variations in wi nd speed and direction (i.e. ‘wind shear’) shred the puff into multiple elements. These new puff parcels, composed of remnants of the old puff, continue to be diffused and dispersed by the wind. The point where significant pu ff shredding occurs is difficult to define sin ce it depends substantially upon the complexity of the me teorological conditions and the underlying terrain. But when shredding occurs, the Lagr angian concept in CA LPUFF breaks down. By ignoring puff shredding (i.e., by keeping puffs intact), the model will systematically over-predict pollutant concentrations.

To deal with this conceptual limita tion, CALPUFF contains an optional puff splitting algorithm to simulate vertical wind shear effects across individual puffs. Differential rates of dispersion and transpor t among the “new” puffs generated from the original, well-mixed puff act to increase the e ffective rate of horiz ontal spread of the material as would be expected in the real atmosphere. Puffs may also be split in the horizontal when the puff size b ecomes large relative to the gr id size to account for wind shear across the puffs. Detailed guidance on when and how the puff-splitting algorithm should be used and actual verification studies demonstrating that the technique operates as intended are not discussed in the mode l documentation or presented in the science literature.J-5 Complex Terrain: Effects of complex terrain on puf f transport are derived from the CALMET winds. In addition, puff-terrain intera ctions at gridded and discrete receptor locations are simulated using one of two algor ithms that modify the puff-height (either that of ISCST3 or a general “plume path coef ficient” adjustment), or an algorithm that simulates enhanced vertical dispersion deri ved from the weakly-stratified flow and dispersion module of the Complex Terrain Disp ersion Model (CTDMPLU S) (Perry et al., 1989). The puff-height adjustme nt algorithms rely on the receptor elevation (relative to the elevation at the source) a nd the height of the puff above the surface. The enhanced dispersion adjustment relies on the slope of the gridded terrain in the direction of transport during the time step.

Subgrid Scale Complex Terrain (CTSG): An optional module, CTSG treats terrain features that are not resolved by the gr idded terrain field, a nd is based on the CTDMPLUS (Perry et al., 1989). Plume impi ngement on subgrid-scale hills is evaluated at the CTSG subgroup of receptors using a dividing streamline height (H d ) to determine which pollutant material is deflected around the sides of a hill (below H d ) and which material is advected over the hill (above H d ). The local flow (near the feature) used to define H d is taken from the gridded CALMET fiel ds. As in CTDMPLUS, each feature is modeled in isolation with its own set of receptors.

Overwater and Coastal Interaction Effects : The CALMET processor contains overwater and overland boundary layer parameterizations allowing certain of the effects of water bodies on plume transport, dispersion, and de position to be estimated. In a sense, CALPUFF operates as a hybrid model, by uti lizing gridded fields of meteorology and dispersion conditions as well as grid-based de scriptions of underlying land use. This includes the abrupt changes that occur at the coastline of a major body of water.

Dry Deposition : A resistance model is used for th e computation of dr y deposition rates of gases and particulate matter as a functi on of geophysical parameters, meteorological conditions, and pollutant species. For partic les, source-specific mass distributions may be provided for use in the resist ance model. Of particular interest for BART analyses is the ability to separately model the de position of fine particulate matter (< 2.5 μ m diameter) from coarse particulate matter (2.5-10 μ m diameter).

Wet Deposition : An empirical scavenging coefficient approach is used to compute the depletion and wet deposition fluxes due to precipitation scavenging. The scavenging coefficients are specified as a function of the pollutant and precipitation type (i.e., frozen vs. liquid precipitation).

3.2.3 Primary Particulates

CALPUFF is designed to simulate PM 10 or PM 2.5 or other user defined size distributions of particles. The smaller the part icles, the more they disperse like an inert gas. In most cases, the dispersion of inert PM 2.5 particles will differ only slightly from that of an inert gas. A key primary PM 2.5 emission from coal-fired electric generating units (EGUs) of relevance to visibility calculations is particulate sulfate. AlthoughJ-6 primary sulfate emissions accoun t for only a small fraction of the total sulfur emissions from such sources, it is appropriate to incl ude their effect if r easonable estimates of primary sulfate emissions from the source are available. Treat ing primary sulfate emissions is likely to be most important at short distances from the stack before significant SO 2 to secondary sulfate conversion has taken place.

3.2.4 Gas-Phase Chemistry

Chemical reactions in the ga s-phase play an important role in secondary aerosol formation by generating radical concentrations (e.g., the hydroxyl radical). These radical species oxidize SO 2 and NOx, providing the precursors to aqueous–phase chemistry (i.e., chemistry in liquid water droplets) that convert SO 2 to sulfate (e.g., H 2 O 2 and O 3 ), and form condensable gases from some volatile organic compounds (VOCs) that can then condense into particulate secondary organic aerosols (SOA). The levels of NOx, VOC, and O 3 concentrations along with the reactivity of the VOCs, sunlight, temperature, and water vapor are all key variables that influe nce the radical cycle and consequent sulfate and nitrate formation rates.

CALPUFF neglects realistic gas-phase pr ocesses entirely. The chemistry in CALPUFF parameterizes chemical transfor mation effects using five species (SO 2 , SO 4 =

, NO x , HNO 3 , and NO 3 - ) via a set of user-specified, diur nally-varying transf ormation rates. The model estimates secondary fine partic ulate matter (sulfate and nitrate) from emissions of gas-phase SO 2 and NO x . Rather than simulating important non-linear gas phase oxidant chemistry, the model employs a user-supplied hourly ozone concentration as a surrogate for the hydroxyl radical and other oxidizing radical species. Ambient ammonia concentrations are also a user input along with temperature and relative humidity.

Although simplifications of photochemistry have been attempted in the past, correct representation of th e gas-phase photochemistry a nd the radical cycles are critically important in order to properly char acterize sulfate and nitrate formation in the real atmosphere. Seigneur et al., (2000) demonstrated this fact in their evaluation of fullscience representations of photochemistry against simplif ied representations (but more advanced than CALPUFF). They conclude d that simplified linearized transformation schemes are inadequate for describing sulfate and nitrate formation processes:

“These results indicate that the accu rate prediction of source-receptor relationships for PM 2.5 requires a comprehens ive treatment of PM 2.5 formation from gaseous precursors for the secondary components of PM 2.5 and a spatially resolved treatment of trans port processes for primary PM 2.5 . Simplified treatments of either atmospheric chemistry or transport are appropriate only when the secondary or primary components of PM 2.5 , respectively, are not significant. Therefore, the development of source-receptor relationships for PM 2.5 should be based on air quality models that pr ovide comprehensive descriptions of atmospheric chemistry and transport.”J-7 Morris et al., (1998) also compared the sulfat e and nitrate particulate estimates from a comprehensive full-science regional model w ith those from a model incorporating a simplified empirical chemical mechanism developed in a manner similar to the mechanism in CALPUFF. Evaluating the fu ll-science and empirical chemistry models against observed concentrations, Mo rris and co-workers concluded:

“Given the importance of the radical cycle for determining secondary PM formation rates, it appears that empirical gas-phase algorithms are inadequate for determining secondary PM formation.”

The uncertainty and potential bi ases introduced into the CA LPUFF visibility estimates due to neglect of gas phase oxidant chemistry remain unknown.

3.2.5 Aerosol Chemistry

Formation of secondary fine particulate matter (e. g., nitrates, sulfates, organic aerosols) in point source plumes is strongly dependent on the rate of mixing with ambient (background) air and the chemical compositi on of this background. The rates of oxidation of sulfur dioxide (SO 2 ) and nitrogen dioxide (NO 2 ) to sulfate and nitric acid can be very different within a power plant or industrial plume compared to that in the background air (Gillani and Godowitch, 1999; Karamchandani et al., 2000). Similarly, the formation of secondary organic aerosol s from emitted VOCs and those from other anthropogenic and biogenic sources, adds yet another pathway in the formation of visibility-impairing aerosols. The presence of atmospheric amm onia introduces further nonlinearities into the gas phase and aerosol reactions. Accordingly, for a model to realistically simulate the production of secondary particulate sulfate, nitrate, and organic aerosols from a potential BART source, th e mixing processes and chemical reactions within and outside of the plume must be treated realistically. If the chemical interactions between these two fundamentally different an d interactive chemical environments are overly-simplified or neglected altogether, the ability of the model to correctly calculate plume concentrations, deposition, or visibility impacts is lost.

Sulfate and Nitrate Formation. Two SO 2 and NOx chemical transformation schemes are available in CALPUFF: the MESOPUFF-II algorithm (Scire et al., 1983; Atkinson et al., 1982) and the RIVAD algorithm (Latimer et al ., 1986). These algorithms calculate sulfate and nitrate formation rates based on the puff concentr ations, background environmental parameters provided by CA LMET, and background ozone and ammonia concentrations provided as input by the user . SOA particulates are not treated by either mechanism. The parameters used are as fo llows (note that each method does not use all of these parameters).

Puff Average Concentrations (from CALPUFF)

• NOx concentration

• SO 2 concentration

Environmental Parameters (from CALMET)J-8 • Temperature

• Surface Relative Humidity (RH)

• Atmospheric Stability

• Solar Radiation

Background Concentrations (User Input)

• Ozone (O 3 )

• Ammonia (NH 3 )

The MESOPUFF-II chemical transformation scheme is EPA’s recommended approach for Class I area impact assessment (IWAQM, 1998). It entails pathways for five active pollutants (SO 2 , SO 4 , NOx, HNO 3 , and NO 3 ) as follows:

k 1 SO 2 Æ SO 4 k 2 NOx Æ HNO 3 (+RNO 3 ) k 3 NOx Æ HNO 3 NH 3 HNO 3 (g) ÅÆ NO 3 (PM)

where,

SO 2 is the puff average sulfur dioxide concentration; NOx is the puff average oxides of nitrogen concentrations; SO 4 is sulfate concentrations formed from the SO 2 ; HNO 3 is the nitric acid formed from the NOx; NO 3 is the particulate nitrate that is in equilibrium with the nitric acid; and
NH 3 is the background ammonia concentration.

Daytime Rates

k 1 = 36 x R 0.55 x [O 3 ] 0.71 x S -1.29 + k 1(aq)

k 1(aq) = 3 x 10 -8 x RH 4 (added to k 1 above during the day)

k 2 = 1206 x [O 3 ] 1.5 x S -1.41 x [NOx] -0.33

k 3 = 1261 x [O 3 ] 1.45 x S -1.34 x [NOx] -0.12

Nighttime Rates

k 1 = 0.20 (%/hr) k 2 = 0.00 (%/hr) k 3 = 2.00 (%/hr)

with,J-9 k 1 is the SO 2 to SO 4 gas-phase transformation rate (%/hr) k 1(aq) is the SO 2 to SO 4 aqueous-phase transformation rate (%/hr) k 2 is the NOx to HNO 3 +RNO 3 transformation rate (%/hr) k 3 is the NOx to HNO 3 (only) transforma tion rate (%/hr) S is the stability index ranging from 2 to 6 (PGT class A&B=2, C=3, D=4, E=5, F=6) R is the total solar radiation intensity (kw/m 2 ) RH is the relative humidity (%) [O 3 ] is the user provided backgr ound ozone concentrations (ppm) [NOx] is the plume average NOx concentration (ppm)
NH 3 is the user provided bac kground ammonia concentrations

Daytime chemical transformations are based on statistically analyzed hourly transformation rates (Scire et al., 1983) obtained from box model simulations using the Atkinson et al., (1982) photochemical mechanis m. In this scheme, gas-phase oxidation of SO 2 and NOx depends on the hydroxyl (OH) radical conc entrations for which background ozone, solar intensity (R), and stab ility index are used as surrogates. At night, OH concentrations are much lower and default SO 2 and NOx oxidation rates of 0.2 %/hr and 2.0 %/hr are assumed. The k 1(aq) sulfate formation rate is added to the k 1 rate during the day as a surrogate for aqueous-phase sulfate formation which begins to assume importance above approximately 50% RH (~0.2 %/ hr sulfate formation rate) and peaks at 100% RH (3%/hr sulfat e formation rate).

The sulfate and nitrate formation rate equations used in the MESOPUFF II scheme were originally generated by deve loping regression equations for a few key variables on the resu lts of 144 box model simulations th at used the 1982 photochemical mechanism of Atkinson et al. These box mode l simulations varied ambient temperature, ozone concentration, sunlight intensity, VOC concentrations , atmospheric stability, and plume NOx concentrations as shown in Tabl e 3-1. The actual environmental conditions used to generate the sulfate and nitrate transformation equations were extremely limited. For example, the transformation rates did not cover temperatures below 10 deg C (50 deg F) or cleaner rural atmospheric conditions with VOC concentrati ons less than 50 ppbC.

The CALPUFF MESOPUFF-II chemistry clea rly neglects several environmental parameters and chemical processes that are important in simulating sulfate and nitrate formation in NO X /SO 2 emissions source plumes. In many cases these deficiencies lead to an overestimation bias of the source’s sulfate and nitrate impacts. Factors that lead such a bias include:

Lack of Temperature Effects : Photochemistry is known to be highly temperature sensitive, as evidenced by the fact that elevated ozone concentrations tend to occur on hot summer days. Lower temperatures produce lower OH and other radical concentrations and consequently lower sulfate and nitrate formation rates. The CALPUFF sulfate and nitrate forma tion rates, however, do not adequately incorporate temperature effects. Th e MESOPUFF-II chemical transformation algorithm was developed under conditions with a minimum temperature of onlyJ-10 10 ° C (50 ° F). Thus, under conditions colder than 10 ° C, CALPUFF will overpredict sulfate and nitrat e formation rates and impacts. CALPUFF typically estimates maximum sulfate and visibility impacts during the late fall/early spring and winter months; these are the same months when the CALPUFF overestimation bias from not considering temperature effects will be greatest. In addition, under colder temperatures, NOx will be converted to peroxyacetyl nitrate (PAN) so that the NOx is no longer available to be conve rted to nitrate. Since the CALPUFF chemistry ignor es the PAN sink for NOx, it will systematically overpredict nitrate impacts.

Effects of NOx Emissions on Sulfate Chemistry : Downwind of a point source with significant NOx/SO 2 emissions, high NOx and SO 2 concentrations co-exist. Under high NOx concentrations, radical concentrations are greatly reduced, resulting in very low ozone, sulfate, and nitr ate formation rates. This is due to the NOx inhibition effect on photochemistry whereby: (1) the titration of NO with ozone eliminates ozone and its source as a radical generator; and (2) the high NO 2 concentrations eliminate the OH radical via the NO 2 + OH reaction thereby effectively shutting down photochemistry. Thus, in a NOx/SO 2 point source plume near the source, there will be very low OH radical and ozone concentrations and consequently very low sulfate and nitrate formation. Since the simple MESOPUFF-II transformation equations cannot account for the NOx effect on the sulfate formation, CALP UFF will tend to over-predict sulfate formation rate in a NOx/SO 2 point source plume near the source, which in turn leads to overstating the sulfate formation rate. Because NOx/SO 2 point sources are typically buoyant, they are frequently be emitted aloft in a stable layer where the high NOx concentrations and inhibi ted sulfate and nitrate formation rates could persist 100 km or more downwind.

Aqueous-Phase Sulfate Formation Algorithm. CALPUFF’s MESOPUFF-II chemistry treats aqueous-phase sulfate formation solely as a function of re lative humidity (RH), which actually has no direct affect on aqueous -phase sulfate formation chemistry. The CALPUFF MESOPUFF-II aqueous-phase sulfate formation rate ranges from values of approximately 0.2 %/hr at 50% RH to 3.0 %/hr at 100% RH. Relative humidity (RH) is a measure of the content of water vapor in th e atmosphere. However, in reality aqueousphase sulfate formation will depend on the amount of atmospheric liquid water content (LWC) in cloud or fog droplets, the pH of the water droplets, and the level of H 2 O 2 , ozone, and SO 2 concentrations. Accordingly, in the atmosphere, aqueous-phase sulfate formation chemistry is not affected by RH. Thus, the CALPUFF aqueous-phase chemistry parameterization is incorrect. Although under conditions of clouds and fog there will be high RH, the occurrence of high RH with very little or no clouds or fog can be quite frequent.

In a liquid water droplet, the reaction of SO 2 with H 2 O 2 to form sulfate is essentially instantaneous and is usually limited by the amount of H 2 O 2 present (i.e., oxidant limited) for a NOx/SO 2 point source. Once the H 2 O 2 is reacted away within the water droplet, sulfate formation via th is pathway slows to the rate of H 2 O 2 formation,J-11 which would be extremely slow to nonexisten t in a large point source plume due to the scavenging of radicals by th e high NOx concentrations. This introduces an inaccurate representation of sulfate formation in CALP UFF that creates uncertainties and bias in modeled visibility impacts. Whether this unc ertainty results in an under- or overestimate of sulfate formation is difficult to determine since the approach is scientifically invalid. Under conditions of high RH and little clouds or little pl ume interaction with clouds, it will clearly overstate sulfat e formation. However, under conditions of cloudy conditions with available photochem ical oxidants (i.e., H 2 O 2 and O 3 ) and a dilute NOx/Sox point source plume, it may understate sulfate formation. Near large NOx/SO 2 point source where the elevated NOx concentrations s cavenge and limit photochemical oxidants, the MESOPUFF-II algorithm will likely overstates sulfate formation.

Thus, the CALPUFF aerosol chemistry fa ils to account for many environmental parameters that are necessary to simulate sulfate and nitrate formation rates, including VOCs and their reactivity, temperature, liquid wa ter content, and NOx concentrations. In their evaluations against full-sc ience PM models and observati ons, Seigneur et al., (2000) and Morris et al., (1998) both independently fo und that the empirical chemistry modules, such as employed by CALPUFF, are inadeq uate for estimating sulfate and nitrate formation

. These findings are supported by EPA’s PM 2.5 and Regional Haze SIP modeling guidance (EPA, 2001) that recommends against using Lagrangian models such as CALPUFF for simulating secondary PM.

From the foregoing, it is clear that the CALPUFF chemical transformation algorithms neglect important chemical processes necessary to accurately estimate the sulfate and nitrate impacts due to SO 2 and NO X emissions. Given that EPA recommends the model for BART determinations, a key question is “What is the influence of the simplified chemistry on modeled estimates of vi sibility impacts from BART sources? In some cases, the inadequacies in the CA LPUFF chemistry algorithms may simply introduce broader uncertainties into the cal culation of estimated sulfate and nitrate impacts. In many cases, however, the simplif ications made in the CALPUFF description of chemical processes result in a systematic bias in the estimated concentrations and visibility impacts due to SO 2 and NO X emissions sources. For large point sources that emit SO 2 and NOx emissions, such as EGUs, petroc hemical process heaters, cement plant kilns, etc., many of the limitations in the CALPUFF MESOPUFF-II SO 2 and NOx transformation algorithms would result in an overestimation bias. While models that are systematically biased high (i.e., over-predict impacts) may be appealing to regulatory decision-makers because they are ‘conserva tive’, the overpredic tion tendency may well lead to unwarranted and excessive control of emissions from some sources. Thus, the tradeoff between simplicity and conservativism on the one hand and te chnical credibility and unbiased answers on the other is a key element in the negotiation of modeling protocols developed by the states or source operators.

3.2.6 Surface Removal

An especially important contributor to pa rticulate concentrations is the rate of deposition to the surface. PM 2.5 particles, which have a ma ss median diameter around 0.5J-12 μm, have an average net deposition velocity of about 1 cm/min (or about 14 m/day) and thus the deposition of fine particles is not usually significant except for ground-level emissions. On the other hand, coarse particles (those PM 10 particles larger than PM 2.5 ) have an average deposition ve locity of more than 1 m/ min (or 1440 m/day), which is significant, even for emissions from elevated stacks.

CALPUFF includes parametric representati ons of particle and gas deposition in terms of atmospheric, deposition layer, and vegetation layer “r esistances” and, for particles, the gravita tional settling speed. Gr avitational settling, which is of particular importance for the coarse fraction of PM 10 , is accounted for in th e calculation of the deposition velocity. Effects of inertial imp action (important for the upper part of the PM 10 distribution) and Brownian motion (important for small, sub-micron particles) and wet scavenging are also addressed. The BART guidance recommends that fine particulate matter (less than 2.5 μm diameter), which has higher light extinc tion efficiency than coarse particulate matter (2.5-10 μm diameter s), should be treated separately in the model. CALPUFF allows for user-specified size categories to be treated as separate species, which includes calculating size-specif ic dry deposition velocities for each size category.

3.3 CALMET Meteorological Preprocessor

The CALMET meteorological model consis ts of a diagnostic wind field module and micrometeorological modules for over- water and overland boundary layers. When modeling a large geographical area such as th e CENRAP domain, the user has the option to use a Lambert Conformal Projection coordinate system to account for Earth’s curvature. The major features and options of the meteorological model are summarized in Table 3-1. The techniques used in the CA LMET model are briefly described below.

3.3.1 Boundary Layer Modules

The CALMET processor contains two boundary layer modules for application to overland and overwater grid cells.

Overland Boundary Layer Module : Over land surfaces, the energy balance method of Holtslag and van Ulden (1983) is used to comp ute hourly gridded fiel ds of the sensible heat flux, surface friction velocity, Moni n-Obukhov length, and convective velocity scale. Mixing heights are determined from the computed hourly surface heat fluxes and observed temperature soundings using a modi fied Carson (1973) method based on Maul (1980). The module also determines gridded fi elds of PGT stability class and hourly precipitation rates.

Overwater Boundary Layer Module: The aerodynamic and thermal properties of water surfaces suggest that a different method is needed for estimating boundary layer parameters in the marine environment. A profile technique, using air-sea temperature differences, is used in CALMET to compute the micro-meteorological parameters in the marine boundary layer. An upwind-looking spatial averaging scheme is optionallyJ-13 applied to the mixing heights and three-dimensional temperature fields in order to account for important advective effects.

3.3.2 CALMET Diagnostic Wind Field Module

The CALMET wind model was constructe d from two other meteorological models used in California in the late 1970s. One was the Calif ornia Institute of Technology (CIT) mass consistent interp olation model described by Goodin et al., (1980). The other was the Complex Terrain Wind Model (CTWM) developed at Systems Applications, Inc. (Tesche and Yocke, 1978; Yocke and Liu, 1978). The CTWM terrain adjustments used to modify the flow fields were assembled in the 1970s as part of research into fire spread and avalanche for ecasting in mountainous regions of California. Various heuristic algorithms were developed to approximate down slope drainage flows, terrain blocking and channeling (Geiger, 1965), thermal heat islands (Stern and Malkus, 1953), surface friction retarda tion, capping by an elevated i nversion and so on. These algorithms were based on empirical stud ies in wind tunnels, numerical modeling experiments, and field studies in the Alps, some dating back to the 1930s (Defant, 1933). Later work by Tesche et al., (1986), Kessler et al., (1987) and Douglas and Kessler (1988) integrated the CIT and CTWM mode ling system into a single meteorological model that included algor ithms to blend observationa l data with prognostic meteorological model output. The combined m odel was used extensively for urban-scale ozone studies throughout the U.S. prior to the switch to MM5 as the preferred meteorological model for SIP studies in the mid-1990s.

The CALMET model development incorporat ed the main features of the CTWM and CIT wind model and significantly updat ed the physical parameterizations and improved model input/output (I/O ) schemes (Scire et al., 2000a). Today, CALMET uses the CTWM two-step approach to the computati on of the wind fields. In the first step, an ‘initial-guess’ wind field is c onstructed and then adjusted to approximate the kinematic effects of terrain, slope flows, and terrain blocking. Currently, the gridded MM5 field is used as the initial guess prio r to terrain-perturbation. The second step consists of an objective analysis procedure to blend the MM5 field with observational data to produce a final wind field. This introduction of obser vational data in the second step of the CALMET wind field development is optional. It is also possible to run the model in “no observations” (No-Obs) mode, wh ich involves the use only of MM5 gridded data for the initial guess field followed by fine-scale terr ain adjustments on the scale of the CALMET domain.

Normally, the CALMET computational domain is specified to be at smaller grid spacing than the MM5 dataset used to initia lize the initial guess field. For example,
36/12 km MM5 data sets available for 2000- 2003 over the CENRAP domain have been used to develop the 6 km CALMET grids shown in Figures 5-1 through 5-4.

The current thermal, kinematic, and dyna mic effects parameterized in CALMET, used in the first step of the windf ield development, are as follows:J-14 Kinematic Effects of Terrain : The CTWM algorithms for kinematic effects (Liu and Yocke, 1980) is used to evaluate the influen ce of the terrain on the wind field. The initial guess field winds are used to compute a terrai n-forced vertical velo city, subject to an exponential, stability-dependent decay functi on. The effects of terrain on the horizontal wind components are evaluated by applying a divergence-minimization scheme to the initial guess wind field. The divergence minimization scheme is applied iteratively until the three-dimensional divergence is less than a threshold value.

Slope Flows : The original slope flow algorithm (D efant, 1933) has been upgraded (Scire and Robe, 1997) based on the shooting flow algorithm of Mahrt (1982). This scheme includes both advective-gravity and equilibriu m flow regimes. At night, the slope flow model parameterizes the flow down the sides of the valley walls into the floor of the valley, and during the day, upslope flows are parameterized. The magnitude of the slope flow depends on the local surface sensible heat flux and local terrain gradients. The slope flow wind components are adde d to the wind field adjusted for kinematic effects.

Blocking Effects : The thermodynamic blocking effect s of terrain on the wind flow are parameterized in terms of the local Froude number (Allwine and Whiteman, 1985). If the Froude number at a particular gr id point is less than a crit ical value and the wind has an uphill component, the wind direction is adju sted to be tangent to the terrain.

3.4 Estimation of Region al Haze Contributions

The default procedure for quantifying visibility impacts is described in several documents (IWAQM, 1998; FL AG, 2000). Implementation of these procedures in CALPUFF is described in t he user’s documentation (Scire et al., 2000b). Generally, ‘visibilit y’ may be quantified either by visual range (the greatest distance that a la rge object can be seen) or by the light extinction coefficient, which is a measure of the li ght attenuation per unit distance due to scattering and absorption by gases and particles . Visibility is impaired when light is scattered in and out of t he line of sight an d by light absorbed along the line of sight. The light extinction coefficient (b ext ) considers light extinction by scattering (b scat ) and absorption (b abs ):

b
ext = b
scat + b abs

The scattering components of extinction (b scat ) are represented by light scattering due to air molecules (i.e., Rayleigh scattering, b rayleigh ) and light scattering due to particles, b sp . The absorption components of extinction (b abs ) include light absorption due to gases (b ag ) and particles (b ap ). Furthermore, particle scattering, b sp , can be expressed by its components:

b
sp = b SO4 + b NO3 + b OC + b SOIL + b Coarse

where the chemical species and soot scattering coefficients are given as:J-15 b SO4 = 3 [(NH
4 )
2 SO
4 ] f(RH)

b NO3 = 3 [NH 4 NO 3 ]f(RH)

b OC = 4 [OC]

b SOIL = [Soil]

b Coarse = 0.6 [Coarse Mass]

b
ap = 10 [EC]

The numeric coefficient at the beginning of each equation is the dry scattering or absorption efficiency in meters-squared per gram. The f(RH) te rm is a monthlyaverage relative humidity adjustment factor . The terms in the brackets are the estimated concentrations fro CALPUFF (or other model) in micrograms per cubic meter (μg/m 3 ).

Finally, the total atmospheric extinction is estimated as:

b
ext = b SO4 + b NO3 + b OC + b SOIL + b Coarse + b
ap + b rayleigh

or, substituting in the above terms,

b ext = 3 f(RH) [(NH 4 ) 2 SO 4 ] + 3 f(RH) [NH 4 NO 3 ] + 4[OC] + 1[Soil] + (3-1) + 0.6[Coarse Mass] + 10[EC] + b Ray

This is the so-called IMPROVE exti nction equation currently recommended by EPA (2003). Note that the sulfate (SO 4 ) and nitrate (NO 3 ) components are hygroscopic because their extinction coe fficients depend upon relative humidity. The concentrations, in square brackets, are in μg/m 3 and b ext is in units of Mm -1

. The Rayleigh scattering term (b Ray ) has a default value of 10 Mm -1 , as recommended in EPA guidance for tracking reasonable progress (EPA, 2003a). The effect of relative humidity variabili ty on the extinction coefficients for SO 4 and NO 3 can be estimated in several ways, but following the EPA BART guidelines, the Class I area-specific m onthly f(RH) values shown in Table 6-1 should be used.

Modeled ground level concentrations of each of the above visibility impairing pollutants are used with the IMPROVE equation to deduce the extinction coefficient. The ch ange in visibility (measured in terms of ‘deciviews’) is compared against background c onditions. The delta-deciview, Δ dv, value isJ-16 calculated from the source’s contribution to extinction, b source , and background extinction, b background , as follows:

Δ dv = 10 ln ((b background + b

source )/ b background )

The impact of a source is determined by comparing the Δ dv, or haze index (HI), for estimated natural background conditions with th e impact of the source and without the impact of the source . If the Δ dv value is greater than the 0. 5 dv threshold the source is said to contribute to visibility impairmen t and is thus subject to BART controls.

CALPOST uses a previous IMPROVE f(RH) curve (FLAG, 2000) which differs slightly from the f(RH) now used by IMPROVE and EPA (2003), mainly at high relative humidity. Also, CALPOST sets the maximum RH at 98% by default (although the user can change it), while the EPA’s guidance now caps it at 95% (easily modified in the CALPUFF input file).

For regional haze light extinction calcula tions, use of a plume-simulating model such as CALPUFF is appropriate only when the plume is sufficiently diffuse that it is not visually discernible as a plume per se , but nevertheless its presence could alter the visibility through the background haze. The IWAQM Phase 2 report states that such conditions occur starting 30 to 50 km from a source. This is consistent with the BART guidance recommendation for us ing CALPUFF for source-recepto r distances greater than 50 km. But, CALPUFF is also recommended by EPA as an option that can be considered for shorter transport distances when the pl ume may in fact be discernible from the background haze.

Apart from the chemistry issues discusse d previously, there do not appear to be any major reasons why CALPUFF cannot be us ed for even shorter transport distances than 30 km, as long as the scale of the plume is larger than the scale of the output grid so that the maximum concentrations and the width of the plume are adequately represented and so that the sub-grid details of plume st ructure can be ignored when estimating effects on light extinction. The standard 1-km output grid that has been established for Class I area analyses should serve down to source- receptor distances somewhat under 30 km; how much closer than 30 km will depend on the topography and meteorology of the area and should be evaluated on a case-by-case basis with individual CENRAP State modelers. (For reference, the width of a Gaussian plume, 2 σ y , is roughly 1 km after 10 km of travel distance, assuming Pasqui ll-Gifford dispersi on rates under neutral conditions.)

3.4.1 CALPOST Methods

Calculation of the impact of the simu lated plume particulate matter component concentrations on light extinction is carried out in the CALPOST postprocessor. For BART applications, this processor is of considerable importance.J-17 CALPOST is used to process the CALP UFF outputs, producing tabulations that summarize the results of the simulations, iden tifying for example, the highest and secondhighest hourly-average concentrations at each receptor. When performing visibilityrelated modeling, CALPOST uses concentra tions from CALPUFF to compute light extinction and related measures of visibil ity (deciviews), reporting these for a 24-hour averaging time. The CALPOST processor contai ns several options for evaluating visibility impacts, including the method described in the BART guidance, which uses monthly average relative humidity values. CALPOS T contains implementations of the IWAQMrecommended and FLAG-recomme nded visibility techniques and additional options to evaluate the impact of natura l weather events (fog, rain and snow) on background visibility and visibility impacts from m odeled sources. CALPOST uses Equation 3-1 to calculate the extinction increment due to the source of interest and provides various methods for estimating the background extinction against which the increment is compared in terms of percent or deciviews.

For background extinction, the CALPOST processor contains seven techniques for computing the change in light extinction due to a source or group of sources (i.e., Methods 1 through 7). These ar e usually reported as 24-hour average values, consistent with EPA and FLM guidance. In addition, there are two techniques for computing the 24-hour average change in extinction (i.e., as the ratio of 24-hour average extinctions, or as the average of 24-hour ratios). Method 2 is the current default, recommended by both IWAQM (EPA, 1998) and FLAG (2000) for sour ce-specific. Method 6 is recommended by EPA’s BART guidance (70 FR 39162).

In Method 2, user-specified, speciated mont hly concentration values are used to describe the background. When applied to natural conditions, for which EPA’s default natural conditions concentrations are annual averages, the same component concentrations would have to be used throughout the year (unless potential refinements to those default values resulted in concentrations that vary during the year). Hourly background extinction is then calculated using these conc entrations and hourly, sitespecific f(RH) from a 1993 IWAQM curve or, optionally, the EPA regional haze f(RH)

curve. 2 Again the RH is capped at either 98% (default ) or a user-selected value (most commonly at 95%).

Method 6 is similar to Method 2, except monthly f(RH) values (e.g., EPA’s monthly climatologically representative values ) are used in place of hourly values for calculating both the extinction impact of the source emissions and the background conditions extinction. Hourly source impacts, w ith the effect on extinc tion due to sulfates and nitrates calculated usi ng the monthly-average relati ve humidity in f(RH), are compared against the monthly default natu ral background concen trations. Thus the monthly-averaged relative humidity is appl ied to the hygroscopi c components (i.e., sulfate and nitrate) of both the source impact and the b ackground extinction with Method
6.

2 Note that the hourly-varying natural background extinc tion here is not consistent with that prescribed by the EPA’s natural conditions guidance (EPA, 2003b) , for which a “climatologically-representative” f(RH) that only varies monthly is to be used. Method 6 uses these monthly average humidity values.J-18 3.4.2 POSTUTIL

The POSTUTIL processor allows the cumulative impacts of multiple sources from different simulations to be summed, including computing the difference between two sets of predicted impacts (useful for evalua ting the benefits of BART controls). It also contains a chemistry module to evaluate the equilibrium relati onship between nitric acid and nitrate aerosols. Th is capability allows the potential non-linear effects of ammonia scavenging by background sulfate and nitrate sources to be approximated in the formation of nitrate from an individual sour ce. The processor can compute the impacts of individual sources or groups of sources on sulfur and nitr ogen deposition into aquatic, forest and coastal ecosystems, thereby allo wing changes in depos ition fluxes resulting from changes in emissions to be quantified.

The POSTUTIL processor attempts to overcome the bias introduced when CALPUFF assumes that the full background a mmonia concentration is entrained into each discrete puff. For a single puff, this may be satisfactory, but the model overestimates the production of ammonium nitrate when multiple puffs co-exist and overlap. The POSTUTIL processor re-p artitions the ammoni a and nitric acid concentrations to conform to the ammoni a-limiting processes influencing nitrate formation. Though based on recognized science, this approximate post-processing method is fundamentally dependent on reliable estimates of ambient NH 3 at the Class I receptor of interest.

3.4.3 Refined Extinction and Back ground Visibility Estimates

EPA, the IMPROVE Steering Committee, and the RPOs are evaluating whether refinements are warranted to the methods recommended for calculating extinction and the default estimate of natural background visibi lity. Whether EPA will approve of any changes to the IMPROVE equation is uncertain at this time. Also, the responsibility for incorporating any changes to the algorithm s in CALPUFF (e.g., new f(RH) curves) is unclear. If changes to these methods are recommended by EPA, CENRAP is encouraged to adopt them. However, details of the proc ess for incorporation of any refinements to the IMPROVE equations in the CALPUFF system should be addressed in the State’s or source operators modeling protocol.

3.5 Model Availability

The EPA-approved version of the CALPUF F modeling system is available from Earth Tech, Inc., ( http://www.src.com/calpuff/calpuff1. htm ). The main models (CALMET, CALPUFF, and CALPOST), their GUIs, and many of the processors are available to download. One may also register to receive no tices of model updates. The most recent update to the system (25 May 2005) is a new version of CALMM5 (MM5 V3) that has been added to the Download BETA-Test page. This version of CALMM5 processes MM5 Version 3 output data directly.J-19 Earth Tech offers CALPUFF training cour ses that include a description of the technical formulation of the models, overvie ws of each of the processor program s, and hands-on application of the models to seve ral case study data sets . Attendees of the course receive a training notebook, a wor kbook of case study problems, exercises, and data sets, updates on recent and future mode l enhancements, and the latest (proprietary) versions of the models and Graphical User Interfaces (GUIs). Othe r third-party training courses and materials are also available.

3.6 CALPUFF Evaluation Studies

Tesche (2002, 2003) reviewed results of various CALPUFF evaluation studies and reached the following conclusions:

> There is a paucity of model ev aluation information for CALPUFF at scales of 50 to 200 km and beyond;

> Based on the limited information available, CALPUFF may be able to give unbiased estimates of short-term (i.e., 3-10 hr) concentrations of nonreactive contaminants to within a factor of two (e.g. 200%) out to distances of about 200 km from a source. This level or uncertainty in a 200 km radius around a source is incr eased if one examines CALPUFF’s predictions in a particular modeli ng cell (e.g., one containing a population center) at a specific hour as opposed to considering the question of bias generally over the entire 200 km regi on irrespective of location and time of occurrence;

> For time periods of a day or less, CALPUFF is unable to produce reliable predictions of non-reactive concentrations at a specific location and time;

> What limited experimental data do exist suggest that the accuracy and reliability of the model’s predictions degrade as the distance scale increases;

> While the IWAQM recommendations on the range of applicability of the CALPUFF model (50 to 200 km) rests on very sparse model evaluation information, EPA’s suggestion that th e model can be used for scales beyond 200 km, even with case-by-case approval, is not based on model evaluation data; and

> For chemically reactive pollutants such as SO 2 , NO x , sulfate, nitrate, nitric acid, and other secondary reaction products, the testing of CALPUFF model over extended spatial scales (50 km and beyond) has not been attempted in a rigorous manner.

Scire et al., (2001) report an evaluati on of CALPUFF sulfate, nitrate, light extinction, and sulfur and nitrogen deposition at a Cla ss I areas over a range of source-J-20 receptor distances. In this study, in whic h a large number of sources were modeled simultaneously, sulfate and nitrate predicti ons at the CASTNet monitoring site in Pinedale, Wyoming were evaluated against obs ervations, and light extinction predictions were evaluated using transmissometer m easurements. Wet sulfur and nitrogen predictions were compared to observations at several acid deposition monitoring sites. This study is especially releva nt because it evaluates the performance of the model to predict variables of direct interest in Class I visibility analyses, such as sulfate and nitrate concentrations and light extinction coefficients

More recently, Chang et al., (2003) repor ted an intercomparison of CALPUFF with two other transport and disp ersion models with high resolution fi eld data. CALPUFF predictions for inert SF 6 were compared using two recent mesoscale fi eld datasets: the Dipole Pride 26 (DP26) and the Overla nd Along-wind Dispersion (OLAD). Both fi eld experiments involved instantane ous releases of sulfur hexa fl uoride tracer gas in a mesoscale region with desert basins and m ountains. Tracer concentrations were observed along lines of samplers at distances up to 20 km. CALPUFF predictions were evaluated using the maximum 3-h dosage (concentrati on integrated over time) along a sampling line. At the DP26 sampler array, CALPUFF had mean biases with in 35% and random scatters of about a factor of 3–4. About 50%–60% of the CALPU FF predictions were within a factor of 2 of th e observations. At the OLAD site, the model underpredicted by a factor of 2–3, on average, with random scat ters of a factor of 3–7. Only about 25%– 30% of the CALPUFF predictions of inert SF 6 were within a factor of 2 of observations.

The tracer studies with which CALPUFF tr ansport and diffusion capabilities were evaluated in the IWAQM Phase 2 report were generally over distances greater than 50 km. More recently, model performance has been performed at shorte r distances including a power plant in Illinois in simple terrain at source-receptor distances in arcs ranging from 0.5 km to 50 km from the stack (Strimaitis et al., 1998). Another CALPUFF evaluation study over short-distances is re ported by Morrison et al. (2003). These studies address model performance over s ource-receptor distances from a few hundred meters to 50 km.J-21 Figure 3-1. CALPUFF Modeling Syst em Components. (Scire et al., 2000a)J-22 Table 3-1. Major Features of the CALME T Meteorological Model. (Scire et al., 2000b)

• Boundary Layer Modules of CALMET

- Overland Boundary Layer - Energy Balance Method - Overwater Boundary Layer - Profile Method
- Produces Gridded Fields of: -- Surface Friction Velocity -- Convective Velocity Scale -- Monin-Obukhov Length -- Mixing Height -- PGT Stability Class -- Air Temperature (3-D) -- Precipitation Rate

• Diagnostic Wind Fiel d Module of CALMET - Slope Flows
- Kinematic Terrain Effects - Terrain Blocking Effects - Divergence Minimization
- Produces Gridded Fields of U, V, W Wind Components - Inputs Include Domain-Scale Winds, Observations, and (optionally) Coarse-Grid Prognostic Model Winds - Lambert Conformal Projection CapabilityJ-23 Table 3-2. Major Features of the CALPU FF Dispersion Model (Scire et al., 2000a)

• Source types - Point sources (constant or variable emissions) - Line sources (consta nt or variable emissions) - Volume sources (constant or variable emissions) - Area sources (constant or variable emissions)

• Non-steady-state emissions and meteorological conditions - Gridded 3-D fields of meteorological variables (winds, temperature)

- Spatially-variable fields of mixi ng height, friction velocity, convective velocity scale, Monin-Obukhov le ngth, precipitation rate - Vertically and horizontally-var ying turbulence and dispersion rates - Time-dependent source and emi ssions data for point, area, and volume

sources - Temporal or wind-dependent scal ing factors for emission rates, for all source types

• Interface to the Emission s Production Model (EPM) - Time-varying heat flux and em issions from controlled burns and wildfires

• Efficient sampling functions - Integrated puff formulation - Elongated puff (slug) formulation

• Dispersion coefficient ( σ y , σ z ) options
- Direct measurements of σ v and σ w - Estimated values of σ v and σ w based on similarity theory - Pasquill-Gifford (PG) disp ersion coefficients (rural areas) - McElroy-Pooler (MP) dispersion coefficients (urban areas) - CTDM dispersion coe fficients (neutral/stable)

• Vertical wind shear - Puff splitting - Differential advection and dispersion

• Plume rise - Buoyant and momentum rise - Stack tip effects - Building downwash effects - Partial penetration - Vertical wind shear

• Building downwashJ-24 - Huber-Snyder method - Schulman-Scire method

-

 PRIME method

Table 3-2. Major Features of the CALPUFF Dispersion Model (Concluded).

• Complex terrain - Steering effects in CALMET wind field - Optional puff height adjustment: ISC3 or "plume path coefficient" - Optional enhanced vertical disp ersion (neutral/weakly stable flow in CTDMPLUS)

• Subgrid scale complex terrain (CTSG option) - Dividing streamline, H d , as in CTDMPLUS: - Above H d , material flows over the hill and experiences altered diffusion rates - Below H d , material deflects around the hill, splits, and wraps around the hill

• Dry Deposition - Gases and particulate matter - Three options: - Full treatment of space and time variations of deposition with a resistance model - User-specified diur nal cycles for each pollutant - No dry deposition

• Overwater and coastal interaction effects - Overwater boundary layer parameters - Abrupt change in meteorological conditions, plume dispersion at coastal boundary - Plume fumigation

• Chemical transformation options - Pseudo-first-order chemical mechanism for SO 2 , SO = 4 , NO x , HNO 3 , and NO 3 (MESOPUFF II method) - Pseudo-first-order chemical mechanism for SO 2 , SO = 4 , NO, NO 2 HNO 3 , and NO 3 (RIVAD/ARM3 method) - User-specified diurnal cycles of transformation rates - No chemical conversionJ-25 • Wet Removal - Scavenging coefficient approach

- Removal rate a function of precip itation intensity and precipitation typeJ-26 Table 3-3. Parameter Variations in Box Model Simulations Used to Develop the CALPUFF Sulfate and Nitrate Forma tion Algorithms. (Morris et al., 2003).

Surrogate
Parameter Number of Variations Model Input Parameters And Variations

Season 3 Temperatures of 30, 20 and 10 ° C were used for the, respectively, summer, fall and winter seasons. Diurnally varying clear sk ies solar radiation was assumed for each season corresponding to a latitude
of 40 °

. Background Air Reactivity 4 For the summer season the following four levels of background ozone and VOCs were used: Ozone (ppb) VOC (ppbC)
20 50
50 250
80 500 200 2,000 For fall and winter the ozone concentrations were assumed to be 75% and 50% of the summer levels. Dispersion 2 Two different rates of plume dispersion were used: (1) a stable case with a wi nd speed of 1.5 m/s and; (2) a slightly unstable case w ith a wind speed of 5.0 m/s. Release Time 2 Photochemical box model simulations were performed with release tim es of sunrise and noon. Plume NOx Concentration 3 Initial plume NOx concentrations of 7, 350 and 1400 ppb were used.

# Appendix 9.2B BART Exemption CALPUFF Modeling Input and Output Files

# Appendix 9.2C Subject-to-BART CALPUFF Post Control Modeling Input and Output Files