# Air Quality Analysis

Modeling Protocol

Statoil OCS

Permitting Project

Prepared for:

Environmental Protection Agency,

Region IV

Atlanta, Georgia

On behalf of:

StatOilGulf Services, LLC

Houston, Texas

Prepared by:

ENVIRON Interna tional Corporation

Baton Rouge

, Louisiana

Date:

July 20 12

Project Number:

26

- 27189 A

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A i 

# Contents

Page

1 Introduction 1

2 Modeling Meth odology 2

2.1 Dispersion Model Selection 2

2.2 Meteorological Data 2

2.3 Source Characterization 2

2.3.1 Stack Parameters 3

2.4 Building Downwash 4

2.5 Receptors 4

3 Modeling Analysis 5

3.1 Modeling Scenarios 5

3.2 Pre

- Construction Monitoring 5

3.3 Preliminary and Full Impacts Analysis 6

4 Additional Impacts Analysis 7

4.1 Visibility 7

5 Class I Area Impacts Analysis 8

5.1 Modeling Methodology 8

5.1.1 Dispersion Model 8

5.1.2 Modeling Procedures 9

5.1.3 Modeling Domain 9

5.1.4 Source Characterization 10

5.1.5 Chemical Transformations 10

5.1.6 Building Downwash 10

5.1.7 Receptor Grids 10

5.1.8 Meteorological Data 11

5.2 Modeling Analysis 11

5.2.1 Increment Consumption 11

5.2.2 AQRVs 12

6 References 15

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

ii Project No. 26

- 27189A 

## List of Tables

Table 1: Significant Monitoring Concentrations

Table 2: Applicable PSD Class II SILs, NAAQS and Increments

Table 3 : Class I Area Distances from Statoil Lease Blocks

Table 4: “Q/D” Analysis

Table 5: Class I Area SILs and Increments

T able 6 : Sulfur and Nitrogen Deposition Thresholds

List of Figures

Figure 1: Proposed Modeling Domain

– Statoil OCS

List of Appendices

Appendix A: Proposed CALPU FF Control Parameters

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 1 

# 1 Introduction

StatOilGulf Services, LLC ( Statoil ) is applying for a Title V and Prevention of Significant

Deterioration (PSD) permit from the Unites States Environmental Protection Agency (EPA)

Region IV for the proposed drilling of e xploration wells in the Eastern Gulf of Mexico (GOM)

beyond the 25

- mile seaward boundary of any state. The proposed exploration will consist of the

operation of deepwater drilling rigs and an associated fleet on the Outer Continental Shelf

(OCS) of the GO M. As such, the application will be made under the OCS permitting rules (40

CFR 55). This document presents the proposed modeling protocol for the Air Quality Impacts

Assessment (AQIA) required as part of the permitting process.

The modeling protocol pre sented here is based on EPA guidance from the Revision to the

Guideline on Air Quality Models (November 2005). In addition, the modeling protocol reflects

guidance provided by EPA and FWS during the recent Shell Exploration and Production

Company permitti ng process, and relevant discussions with Stan Krivo (EPA), Kelly Fortin

(EPA), and Jill Webster (FWS)

.

All regulated PSD pollutants will be subject to PSD review. The P SD regulated pollutants

include particulate matter (PM), carbon monoxide (CO), nitroge n dioxide (NO 2 

# ), sulfur dioxide

(SO 2 

# ) and ozone (O 3 

# ). The AQIA for this project will address all pollutants that result in a

significant emissions increase with respect to PSD. At this point, it is assumed the project will

only be significant for NOx . At this time, Statoil will not address CO 2 

# in the AQIA because there

are no existing or proposed ambient standards for CO 2 

# or other greenhouse gases.

Statoil is requesting flexibility to conduct exploratory drilling at multiple locations in the Eastern

GOM. Figure 1 shows the location of Statoil’s lease blocks (Desoto Canyon lease blocks 143,

187, 188, 230

, 231 , 625, 669, 670, 671, 715, 716, 759, 760, and 804 )

. Several of these lease

blocks have been added to the list after the May 2, 2012 pre

- application meeting, after recent

auctions. The added lease blocks are all south of the original list of lease blocks, farther from

the Louisiana/Mississippi/Alabama coastline.

Statoil also proposes the use of multiple rigs, namely the Transocean Discoverer Americ as

( Transocean ), or the Maersk Developer (Maersk) drilling rigs. The Maersk is a semi

submersible drilling rig on two pontoons, yet is dynamically positioned. T he Transocean rig is a

dynamic ally positioned drilling vessel with a single hull

. Neither requ ires towing vessels or

anchor handling vessels. Statoil will utilize the worst

- case emissions and modeling parameters

from these rigs in the AQIA to ensure all possible scenarios are addressed.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

2 Project No. 26

- 27189A 

# 2 Modeling Methodology

EPA established a policy in the 2011 Shell E&P permit for De Soto Canyon, whereby the

“property boundary” of OCS sources is set at 25 nautical miles from the respective State’s

seaward boundary.



# 2.1 Dispersion Model Selection

Because the Statoil lease blocks are more than 50 km from the line 2 5 nautical miles from the

States’ seaward boundary, AERMOD cannot be used (40 CFR 51 Appendix W, 4.2.2(b)). For

this situation, the CALPUFF modeling system is appropriate.

CALPUFF version 5.8, the regulatory default version, will be used for both the Clas s I AQRV

and Increment analyses and the Class II NAAQS and Increment analyses . Using the

conservative assumption that the maximum 1

- hour emission rates are the same as the

maximum 24

- hour emission rates, a single CALPUFF run can be used to assess complian ce at

both Class I and Class II receptors for all averaging periods ≤ 24 hours. A separate CALPUFF

run will be used to asses s annual averages, for both Class I and Class II receptors.



# 2.2 Meteorological Data

Statoil proposes to use the three

- year meteorological dataset (2001

– 2003) for

CALMET/CALPUFF that covers the Gulf of Mexico (GOM) region, which was prepared for

VISTAS BART modeling 1 

# . Specifically, Statoil proposes to use Domain 2, which is centered on

Florida but extends far enough to the west to reach 50 km wes t of the Breton Wilderness.

The CALMET/CALPUFF input data files were constructed using a combination of on

- shore

surface and upper

- air data from the National Weather Service (NWS), offshore buoy data from

the National Data Buoy Center (NDBC), and model o utput from the MM5 model. Ozone data

were also compiled as part of this dataset

, and were obtained from the EPA Air Quality System

(AQS). Statoil will use the most recent version developed by the FLM using the regulatory

version of CALMET.



# 2.3 Source Chara cterization

The primary emission sources consist of engines on the Rigs and mobile equipment, such as

supply ships and life boats. OCS permitting rules consider direct emissions from sources

located within 25 miles of the facility. These rules require St atoil to consider impacts from the

associated fleet of vessels that support the drilling unit.

Maersk Developer emission sources include the following:

• 

# Main Generator Engines (8 )

1 

## http://www.vistas

- sesarm.org/BART/calpuff.asp

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 3



## • 

# Standby Generator Engine (1)

• 

# Cementing Unit Engines (2)

• 

# Life Boats (4)

• 

# Resc ue Boat (1 )

• 

# Fuel Storage Vessels

• 

# Fugitives (diesel fuel system)

• 

# Support Vessels (e.g. crew boats)

Discoverer Americas emission sources include the following:

• 

# Main Generator Engines (6 )

• 

# Emergency Generator Engine (1)

• 

# Life Boats (6 )

• 

# Rescue Boat (1)

• 

# Emergency Air Compressor Engine s ( 2 )

• 

# Fuel Storage Vessels

• 

# Fugitives (diesel fuel system)

• 

# Support Vessels (e.g. crew boats)



# 2.3.1 Stack Parameters

Stationary sources, such as the diesel engines, may have vertical

, horizontal , or angled

discharge orientation. For stacks with a vertical discharge, Statoil will use actual stack

parameter data (i.e., stack diameter, stack exit velocity and stack exit temperature). For stacks

with a horizontal or angled discharge, Statoil will adjust the stack parameters to remove any

moment um rise, but maintain thermal rise , by using a stack exit velocity of 0.001 m/s. Statoil

will use the ac tual stack exit temperature, diameter, and stack height of each point source.

Statoil proposes to model mobile sources (crew boats, offshore support ve ssels, etc.) as a

series of volume source s. These will be located along the line between the main drilling vessel

and a point 25 miles away, in the direction of the B reton Wilderness Area

. Ten (10) volume

sources will be used to represent the combined em issions from all relevant mobile sources.

Volume sources are defined by four parameters:

• 

# An emission rate

• 

# A release height (called an “effected emission height”)

• 

# An initial horizontal dimension called the initial sigma

- Y

• 

# An initial vertical dimension call ed the initial sigma

- Z

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

4 Project No. 26

- 27189A 

# According to Section 1.2 in the User’s Guide for the Industrial Source Complex (ISC3)

Dispersion Models, Volume II

- Description of Model Algorithms (EPA

- 454/95

- 003b, September

1995), the effective release height should be based on t he height of the emissions including any

plume rise for buoyancy of the emissions. The plume heights for the fleet emissions will be

estimated using SCREEN3, which accounts for the mechanical and buoyant lift from the ship’s

stacks. 2 

# The stack characteri stics of the main propulsion engines for each ship will be used with

the SCREEN3 algorithms to define the plume height and sigma Z values for the fleet emissions

.

Initial dispe rsion for volume sources is characterized by two parameters, σ

y 

# (sigma Y) and σ

z 

# (sigma Z). For the fleet, the sigma Y value for each volume source will be determined by

dividing the physical horizontal dimension of the volume by 2.15 as recommended in th e ISC

User’s Guide.

The base elevation of all sources will be referenced to sea level. For example, stationary

sources located on the deck of the drilling unit will have a base elevation corresponding to the

height of the deck above the operating wate r line. Stack parameters for all emission sources

are currently under development and will be presented in detail in the permit application.



# 2.4 Building Downwash

Given that the nearest receptors are located more than 175 km from the source, building

downwash effects will not significantly affect the modeled results. The FLMs typically do not

request downwash be included, and neither “downwash” nor “BPIP” are mentioned in FLAG

(2010).



# 2.5 Receptor s

Statoil proposes to use a single domain for both the Class II and Class I modeling: a subset of

the VISTAS CALMET Domain 2 shown in Figure 1 . Class II receptors have been placed

25 nautical miles from the seaward boundary of the four States, and are shown in red. Class I

receptors were obtained from the National Park Service (see Section 5).

2 

## SCREEN3 is an alternate EPA

- approved screening model that uses the same plume rise equations as

ISCST3 and AERMOD but also provides a printout of plume rise calculations.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 5 

# 3 Class II Modeling Analysis

Maximum modeled impacts will be compared to the applicable ambient air quality standards and

PSD Class II increments. This approach is presented below.



# 3.1 Modeling Scenarios

Statoil is requesting the fle xibility to operate either the Maersk Developer or the Transocean

Discoverer Americas drilling units at any given time. All the drilling units can operate under

varying engine load scenarios and with various fleet vessel support operating scenarios. T he

modeling analysis will include two scenarios to adequately address worst

- case impacts:

• 

# Maersk rig operating under a range of emissions and release parameters (including partial

loads) and with associated support vessels

• 

# Transocean rig operating under a range of emissions and release parameters (including

partial loads) and with associated support vessels

Statoil will estimate maximum hourly (lb/hr) and annual (ton/year) emission rates for each

drilling vessel. This leads to a total of four modeling scen arios: short

- term and long

- term for

each of the two proposed vessels.

Statoil will be requesting permits for portable major sources to allow the Maersk or the

Transocean and their associated fleet to operate in any of Statoil’s lease blocks located withi n

the Eastern GOM. As such, Statoil will represent worst

- case dispersion conditions as follows:

Each lease block has dimensions of approximately five kilometers by five kilometers. The

closest point to shore is the northwest corner of the Desoto Canyon 1 43 lease block, which lies

at 28.851 o 

# N and 87.424 o 

# W . For purposes of this modeling analysis, Statoil will model all

scenarios from this northwest corner

. S tatoil expects impacts to be the same for all locations

within the lease blocks. This location h as been selected to represent the lease blocks as a

whole because it is the closest location to shore and the nearest Class I area (Breton).



# 3.2 Pre

- Construction Monitoring

Pre

- construction ambient monitoring may be required for any criteria pollutant that is proposed

to be emitted above the significant emission rates (or 100 tons per year of volatile organic

compounds [VOCs]). Statoil proposes to emit nitrogen oxides (NO

x 

# ) in significant amounts.

If the predicted concentrations from direct emission sources (i.e., the highest modeled

concentration for the applicable averaging time, caused by the proposed significant emissions

increase) are above the significant monitoring concentrations (SMCs) and EPA Region IV

determines that ambient monitoring is required, Statoil will satisfy such a requirement by using

existing ambient monitoring data (refer to Section 3.6 of this modeling protocol for details).

SMCs are listed in Table 1. No monitoring is required if modeling results of the significance

analysis are bel ow the defined SMCs for a given pollutant and time period.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

6 Project No. 26

- 27189A 

# Table 1. Significant Monitoring Concentrations

Pollutant Averaging Time Air Quality Concentration ( μ g/m 3 

## )

NO 2 

## Annual 14



# 3.3 Preliminary and Full Impacts Analysis

PSD modeling analysis involves two phases: a preliminary analysis (referred to as a Significant

Impact Analysis) and, if required, a Full Impact Analysis. The preliminary analysis estimates

am bient concentrations resulting from the proposed project for pollutants that trigger PSD

review. Because Statoil proposes to emit NO

x 

# in excess of 40 TPY, a PM 2.5 

# review is also

triggered.

The results of the preliminary analysis determine whether a full i mpact analysis (facility plus

background and the concentrations from regional sources) for a particular pollutant is

necessary. If the ambient concentrations from the preliminary analysis are greater than the

PSD Significant Impact Levels (SILs), then a F ull Impact Assessment is required.

The EPA has not identified a regulatory

- default model for secondary PM 2.5 

# . Therefore, only

concentrations resulting from proposed primary PM 2.5 

# emissions will be compared to the PM 2.5



# SILs.

If the preliminary analysis indicates that either a Class I or Class II full impact analysis is

necessary, a separate modeling procedure will developed in conjunction with the appropriate

agency representative (EPA or FLM).

Table 2 lists the applicable Class II SIL, NAAQS and Class II PSD increments.

Table 2. Applicable PSD Class II SILs, NAAQS and Increments

Pollutant Averaging

Period Class II Significant

Impact Level

( μ g/m 3 

## ) NAAQS

( μ g/m 3 

## ) Class II PSD

Increment

( μ g/m 3 

## )

NO 2 

## 1

- hour 7.5 (1) 

## 188 ( 2 ) 

## N/A

Annual 1 100 25

PM 2/5 

## 24

- hour 1.2 15 (3) 

## 9

Annual 0.3 35 4

1 

## Interim SIL, established June 28, 2010 in the memo “Guidance Concerning the Implementation of the 1

- hour NO 2



## NAAQS for the Prevention of Significant Deterioration Program”

. 2 

## To attain this standard, the 5

- year averag e of the 98th percentile of daily maximum 1

- hour concentrations must not

exceed 188 μg/m 3



## , effectiv e April 12, 2010

. 3 

## To attain this standard, the 5

- year average of the 98th percentile of 24

- hour concentrations must not exceed

35 μg /m 3



## .

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 7 

# 4 Additional Impa cts Analysis

40 CFR 52.21 requires an Additional Impacts Analysis for compounds under PSD review. The

Additional Impacts Analysis includes the analysis of the project’s impact on the visibility of the

surrounding area, growth and development in the surroun ding area, and soil and vegetation in

the surrounding area. However, since the site is located offshore and there is no development

or vegetation within 30 km of the project, the additional impacts analysis will be restricted to

visibility impacts only.



# 4.1 Visibility

The project’s potential impact on visibility surrounding the facility will be investigated according

to the procedures in EPA’s Workbook for Plume Visual Impact Screening and Analysis

(Revised). 3 

# This visibility analysis is distinct from the an alysis required for a Class I area. A

Level I screening analysis will be performed using the latest version of VISCREEN. If a Level I

screening analysis does not pass the appropriate thresholds, a Level II analysis will be

implemented by replacing the de fault meteorology with data obtained from the representative

meteorological station ( Tampa West buoy [WBAN# 42039 ] and Panama City NWS station

[WBAN# 03882 ]).

3 

## United States Environmental Protection Agency (USEPA). 1992. Workbook for Plume Vi sual Impact Screening and

Analysis. EPA/454/R

- 92

- 023.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

8 Project No. 26

- 27189A 

# 5 Class I Area Impacts Analysis

PSD guidance requires an analysis of potential impacts to air qual ity related values (AQRVs) in

Federal Class I areas within 100 km (62.1 miles) of the facility subject to PSD review. However,

for most applications, the Federal Land Managers (FLMs) request analyses of AQRV impacts

for additional Class I areas within 200

– 300 km of the site. The locations of all nearby Class I

areas in relation to the facility are shown i n Figure 1. As shown in Table 3 , only the Breton

Wilderness Area is located within 300 km of the facility.

Table 3. Class I Area Distances from Statoi l Lease Blocks

Class I Area Distance 1



## (km) (miles)

Breton Wilderness Area (LA) 175 109

Bradwell Bay Wilderness Area (FL) 307 191

Saint Marks Wilderness Area (FL) 313 194

Chassahowitzka Wilderness Area (FL) 461 286

Notes:

1 

## Distance to northwest co rner of Desoto Canyon 143 lease block (28. 85064 °N, 87. 42378 °W)



# The Breton, Chassahowitzka and Saint Marks Wilderness Areas are managed by the Fish and

Wildlife Service (FWS). The Bradwell Bay Wilderness Area is managed by the Forest Service

(FS). S tatoil will share this modeling protocol with the FWS and FS to obtain their concurrence

with the proposed approach. Statoil has already made initial contact with the FWS to discuss

the proposed project and AQRV analysis. 4



# 5.1 Modeling Methodology

5.1.1 Dispersion Model

On April 15, 2003, EPA adopted the CALPUFF modeling system as the EPA’s preferred model

for long

- range transport assessments and for evaluating potential impacts to Class I areas by

including CALPUFF in Appendix A of the Guideline (November 2005). F eatures of the

CALPUFF modeling system include the ability to consider: secondary aerosol formation;

gaseous and particle deposition; wet and dry deposition processes; complex three

- dimensional

wind regimes; and the effects of humidity on regional visibili ty. Statoil proposes to use the

CALPUFF modeling system for modeling Class I area impacts , using the following elements in

this analysis:

• 

# CALPUFF version 5.8, level 070623

4 

## Bart Brashers of ENVIRON had phone discussions with Ms. Jill Webster of FWS in early March, 2012.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 9



## • 

# CALMET version 5.8, level 070623

• 

# CALPOST version 6.221, level 080724 (AQRV analysis only)

• 

# CALPOST version 5.6394, level 070622 (PSD increment analysis only)



# 5.1.2 Modeling Procedures

Statoil will follow the recommendations of the Interagency Workgroup on Air Quality Modeling

(IWAQM) and the FLM Air Quality Related Values Workgroup (FLAG), outl ined in the FLAG

Phase I Report (December 2000 , Revised October 2010 ). EPA endorsed these procedures in

advance of the IWAQM Phase II draft report (December 1998), and reiterated the endorsement

in the April 15, 2003 Federal Register notice (Volume 68, Nu mber 72) that adopted CALPUFF

as a Guideline model. Proposed CALPUFF model selections are presented in Appendix A.

Statoil proposes to model worst

- case emission rates from the facility at the northwest corner of

the Desoto Canyon Block 143 as described in Section 3.1. This location is the closest to any of

the Class I areas.

Based on initial estimates of the maximum 24

- hour emissions of NO

x 

# , SO 2 

# , H 2 

# SO 4 

# , and PM 10 

# for

each drilling vessel and associated fleet, an initial “Q/D” analysis yields values in ex cess of 10

(Table 4). This indicates that the magnitude of emission do not “pass” the initial screening test

recommended by FLAG (2010) and the full Class I analysis is required.

Table 4. “Q/D” Analysis

Class I Area Distance 1



## (km) Q/D

Maersk Transocea n

Breton Wilderness Area (LA) 175 14 20

Notes:

1 

## Distance to northwest corner of Desoto Canyon 143 lease block (28. 85064 °N, 87. 42378 °W)



# 5.1.3 Modeling Domain

The size of the proposed modeling domain is shown in Figure 1

, and extend s at least 50 km

beyond each Class I area and the facility to allow for the possible recirculation of puffs beyond

the facility and Class I areas. Due to the large size of the proposed domain, it uses the Lambert

conformal projection to represent the locations of modeled source s, grid cells and receptors.

Consistent with grid spacing in the VISTAS CALMET Domain 2 dataset, Statoil proposes to use

a grid spacing interval of 4 km. The vertical domain is also composed of terrain

- following grid

cells. The number and size of these c ells are selected to constrain the boundary layer in which

dispersion and chemical transformations occur. Consistent with the CALMET meteorological

datasets, Statoil proposes to use 10 vertical grid cells with cell face heights set at 0; 20; 40; 80;

160; 320; 640; 1,200; 2,000; 3,000; and 4,000 meters above the modeled ground surface.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

10 Project No. 26

- 27189A 

# 5.1.4 Source Characterization

Source characterization for CALPUFF will be identical to that developed for the Class II

analysis . Source locations will be in Lambert Conformal Coor dinate (LCC) system, using the

reference coordinate system of the VISTAS CALMET Domain 2 data grid.

The CALPUFF modeling will consider emissions of NO

x 

# , SO 2 

# , and PM for the evaluation of

AQRVs at the Class I areas. 5 

# However, the CALPUFF modeling will ca lculate only

concentrations of NO

x 

# and primary PM 2.5 

# for the evaluation of PSD increment consumption at

the Class I areas.



# 5.1.5 Chemical Transformations

5.1.5.1 Ammonia

The NO

x 

# chemistry in CALPUFF depends on the ambient ammonia concentration to establish

the equilibri um between gaseous nitric acid and ammonium nitrate. However, ambient

ammonia concentrations are not explicitly simulated by CALPUFF and the user must select an

appropriate background concentration. The IWAQM Phase II recommendations suggest typical

ammo nia concentrations are: 10 parts per billion (ppb) for grasslands; 0.5 ppb for forests; and 1

ppb for arid lands during warmer weather. Given the low background concentration over the

Gulf of Mexico and as requested by the FLM, Statoil proposes a value of 3 ppb.



# 5.1.5.2 Ozone

Reaction rates in the CALPUFF chemistry algorithms are also influenced by background ozone

concentrations. For this analysis, Statoil proposes to use hourly ozone data to represent

background concentrations. During hours of missing ozone da ta, a constant background ozone

concentration of 65 ppb will be substituted as representative of the average background

concentration for this region.



# 5.1.6 Building Downwash

Statoil will not include buildings and structures located on the drilling rig and the a ssociated fleet

vessels. Building downwash is the effect of nearby structures on the flow of emissions from

their respective sources. Because of the relatively large distance between the emissions and

the Class I receptors, building downwash will not be included

.



# 5.1.7 Receptor Grids

The CALPUFF dispersion model simulations assess AQRVs within each Class I area at discrete

receptors obtained from the National Park Service. 6 

# In addition to these discrete receptor

locations, Statoil may use a receptor grid with 4

- km spacing throughout the CALPUFF modeling

domain for AQRV predictions. Statoil may use the 4

- km mesh size receptors to construct plots

showing the spatial variation of the calculated parameters throughout the modeling domain.

5 

## There are no Class I increments for CO, and CO does not contribute to visibility impairment.

6 

## http://www2.nature.nps.gov/air/Maps/Receptors/index.cfm

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 11 

# Statoil will use these pl ots for diagnostic purposes, as well as to develop the figures presented

in the permit application. Comparisons with AQRV criteria will be based on the discrete

receptor locations and elevations rather than the domain

- wide gridded receptors.



# 5.1.8 Meteorologica l Data

Preprocessed VISTAS

- BART 4 km Domain 2 CALMET data will be used for modeling. This

data was provided by the North Carolina Department of Environment and Natural Resources

(NCDENR). Data from 109 surface stations, 10 upper air stations, 9 overwater stations and

96 precipitation stations was used to process the 2001 CALMET data. Data from 100 surface

stations, 10 upper air stations, 9 overwater stations and 92 precipitation stations was used to

process the 2002 CALMET data. Data from 116 surface stat ions, 10 upper air stations, 9

overwater stations and 103 precipitation stations was used to process the 2003 CALMET data.

This dataset is the most recent version developed by the FLM using the currently approved

regulatory version of CALMET (Version 5.8 ). For the current analysis

, the 4

- km mesh size

simulations were used.

CALMET (Version 5.8), the meteorological preprocessor component of the CALPUFF system,

was used to combine the MM5 simulation data, surface observations, terrain elevations, and

land u se data into the format required by the dispersion modeling component CALPUFF. In

addition to specifying the three

- dimensional wind field, CALMET also estimates the boundary

layer parameters used to characterize diffusion and deposition by the dispersion model.

Selected hours of the three

- year CALMET/MM5 three

- dimensional data set were examined by

extracting data from the CALMET output files and plotting the meteorological fields with the

CALDESK software package. Wind vector plots were examined for diff erent times of year,

different times of day, and for all 10 vertical levels.



# 5.2 Modeling Analysis

The impacts on air quality and air AQRVs estimated by CALPUFF will be compared against

applicable threshold levels. The two types of threshold levels are: (1) no t to be exceeded levels

(i.e., PSD increments) or (2) levels of potential concern (levels that, when exceeded, raise

concerns and should be evaluated for their significance and reliability).



# 5.2.1 Increment Consumption

As part of the USEPA’s PSD policy, Class I area allowable concentration increments have been

established. The cumulative air quality impact of all new sources in a common region are

required to be below the allowable PSD Class I increments at Class I areas. A new source that

causes the cumulative impact of all new sources to exceed a PSD Class I allowable increment

may not be allowed to be built.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

12 Project No. 26

- 27189A 

# In 1996, the USEPA published a Federal Register notice of proposed Class I area SIL

thresholds for single projects. 7 

# These proposed single project SIL s are defined as being

approximately 4% of the PSD Class I area increment. If a project’s impact is below the Class I

area single project proposed SIL thresholds, then its impacts can be assumed to be

insignificant.

If estimated concentrations at Class I areas exceed the proposed single source SIL then Statoil

will perform a cumulative PSD Class I increment compliance assessment. Table 5 summarizes

the PSD Class I area increments and the single project proposed Class I area SILs that are

applicable to t he proposed project. If maximum modeled concentrations exceed the SIL for any

pollutant averaging period, then a cumulative analysis for the pollutant would be required. In

this case, Statoil would consult with the FLM and EPA to develop the emissions inv entory.

The EPA has not identified a regulatory

- default model capable of calculating secondary PM 2.5



# . Therefore, for the SILs comparison that might trigger a cumulative PM2.5 ana lysis, only primary

emissions of PM 2.5 

# will be used. For all other AQRVs (evaluated by the FLMs), both primary

and sec ondary PM 2.5 

# will be included.

Table 5 . Class I Area SILs and Increments

Pollutant Averaging Period Significance Level

( μ g/m 3 

## ) Class I Increment

( μ g/m 3 

## )

NO 2 

## Annual 0.1 2.5

PM 2.5 

## 24

- hour 0.07 2

Annual 0.06 1



# 5.2.2 AQRVs

5.2.2.1 Visibility

The FLAG workgroup recommends procedures for est imating the visibility impacts due to

proposed new sources at Class I areas using refined CA LMET/CALPUFF modeling (FLAG,

201 0). The FLAG visibility metric is the estimated maximum 24

- hour change in extinction ( ∆ b

ext 

# )

over clean natural visibility conditio ns (Natural Conditions) at the Class I area. The

interpretation of the FLAG thresholds for extinction change over natural background is as

follows:

• 

# If the source’s visibility impact is < 0.4% on all days, the source is considered insignificant

and the FLM will not object to the permit.

7 

## Federal Register: July 23, 1996 (Volume 61, Number 142)

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 13



## • 

# If the source’s visibility impact is < 5% on all days, the FLM will likely not object to the

permit.

• 

# If there are days with the source’s visibility impact > 10%, the FLM may object to the

permit.

• 

# If there are days in which t he source’s visibility impact are above 5% the frequency,

magnitude and duration of the visibility impacts are used to make a significance

determination.

If a source exceeds a specific threshold at a Class I area, then the frequency, magnitude and

duration of the impacts and the reliability and accuracy of the modeling are examined to

interpret the modeling results. More recent interpretation of the FLAG procedures for evaluating

the visibility impacts estimated by the CALMET/CALPUFF modeling system has al lowed the

introduction of extenuating circumstances that account for natural obscuration of visibility and

modeling artifacts introduced by the CALMET/CALPUFF modeling system and the combination

of the two.

The FLAG (2010) procedures employ the IMPROVE ext inction equation to calculate b

ext 

# (Method 8, invoked with MVISCHECK=1 in CALPOST). This updated equation for extinction

uses monthly relatively humidity adjustment factors with relative humidity capped at 95%. It

uses annual background aerosol concentra tions recommended by the FLMs for the Class I area

of concern, and assesses the visibility using the 98 th 

# percentile modeled values at each

receptor. In order to use Method 8, CALPOST Version 6.221 (Level 080724) w ill be used to

post

- process the CALPUFF output files.

For FLAG 2010, a Project’s extinction is calculated using the revised IMPROVE reconstructed

mass extinction equation as follows:

b project 

# = 2.2 × f S 

# (RH) × [Small Sulfate] + 4.8 × f L 

# (RH) × [Large Sulfate]

+ 2.4 × f S 

# (RH) × [Small Nitrate] + 5.1 × f L 

# (RH) × [Large Nitrate]

+ 2.8 × [Small Organic Mass] + 6.1 × [Large Organic Mass]

+ 10 × [Elemental Carbon]

+ 1 × [Fine Soil]

+ 0.6 × [Coarse Mass]

+ 1.7 × f SS 

# (RH) × [Sea Salt]

+ Rayleigh Scattering (Site Specific)

+ 0.33 × [NO 2 

# (ppb)] {or as: 0.1755 × [NO 2 

# ( μ g/m 3 

# )]}

Where:

[ ] indicates concentrations in μg/m 3



# f S 

# (RH) = Relative humidity adjustment factor for small sulfate and nitrate

f L 

# (RH) = Relative humidity adjustment factor for large sulfate and nitrate

f SS 

# (RH) = Relative humidity adjustment factor for sea salt

For Total Sulfate < 20 μg/m 3 

# : [Large Sulfate] = ([Total Sulfate] / 20 μg/m 3 

# ) × [Total Sulfate]

For Total Sulfate ≥ 20 μg/m 3 

# : [Large Sulfate] = [Total Sulfate]

And:

[Small Sulfate] = [Total Sulfate]

– [Large Sulfate]

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

14 Project No. 26

- 27189A 

# To calculate large and small nitrate and organic mass, substitute ({Large, Small, Total}

{Nitrate, Organic Mass}) for Sulfate.



# 5.2.2.2 Deposition

Acid deposition impacts are represented by total sulfur and total nitrogen deposition at the Class

I areas. The FLAG procedures r equire estimation of total sulfur deposition from the CALPUFF

estimated wet and dry SO 2 

# plus sulfate (SO 4 

# ) deposition. For nitrogen, wet and dry deposition

from all of the nitrogen species modeled species are included (NO 2 

# , nitric acid, and particulate

ni trate).

The Forest Service (FS) has developed sulfur and deposition thresholds “...below which a land

manager can recommend that a permit be issued” (USDA FS, 1989). Although these values

vary for different locations, the lowest “green line” sulfur and nitr ogen deposition values are 3

kilograms per hectare per year (kg/ha/yr). More recently, the NPS has posted a document

“Guidance on Nitrogen and Sulfur Deposition Analysis Thresholds” on their website. 8 

# The NPS

Deposition Analysis Thresholds (DATs) for nit rogen and sulfur deposition are as follows:

• 

# East DAT: 0.010 kg/ha/yr

• 

# West DAT: 0.005 kg/ha/yr

“East” and “West” refer to Class I areas east and west of the Mississippi River. The east DAT is

applicable t o the proposed Project. Table 6 lists the sulfur and nitrogen deposition thresholds

that the proposed project deposition estimates would be compared against.

Note that because the proposed sources would not emit more than 40 TPY of SO 2 

# , Statoil will

not present the sulfur deposition results in the appli cation.

T able 6 . Sulfur and Nitrogen Deposition Thresholds

Threshold Source Average Sulfur Deposition

(kg/ha/yr) Average Nitrogen Deposition

(kg/ha/yr)

FS Threshold 3.0 3.0

NPS DAT 0.01 0.01

8 

## http://www2.nature.nps.gov/air/Pubs/pdf/flag/nsDATGuidance.pdf

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 2 7189A 15 

# 6 References

FLAG. 2000. “Federal Land Managers’ Air Qual ity Related Values Workgroup (FLAG)”: Phase

I Report. USDI

– National Park Service, Air Resources Division, Denver, CO.

FLAG. 2010. “Federal Land Managers’ Air Quality Related Values Workgroup (FLAG)”: Phase

I Report

– Revised. USDI

– National Park Service, Air Resources Division, Denver, CO.

IWAQM. 1998. Interagency Workgroup on Air Quality Modeling (IWAQM) Phase 2 “Summary

Report and Recommendations for Modeling Long

- Range Transport and Impacts on

Regional Visibility.” EPA

- 454/R

- 98

- 019. U.S. Environmental Protection Agency, Office of

Air Quality Planning and Standards, Research Triangle Part, NC.

USEPA. 2003a. Revision to the Guideline on Air Quality Models: Adoption of a Preferred Long

Range Transport Model and Other Revisions; Final Rule. 40 CFR Part 51, Federal Register

, Vol. 68, No. 72, Tuesday April 15, 2003, Rule and regulations pp 18440

- 18482.

USEPA. 2003b.“Guidance for Estimating Natural Visibility Conditions Under the Regional Haze

Rule”, U.S. Environmental Protection Agency, Resea rch Triangle Park, NC.

USEPA. 1992. “Workbook for Plume Visual Impact Screening and Analysis,” EPA/454/R

- 92

-

023.

USEPA. 1985. “Guidelines for Determination of Good Engineering Practice Stack Height

(Technical Support Document for Stack Height Regulatio ns (Revised)),” EPA

- 450/4

- 809

-

023R.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 17 

# Figures

## 18 Project No. 26

- 27189A 100 km

200 km

300 km 

## Full Vistas Domain

Closest Lease Block Statoil OCS CALPUFF Domain

Brigantine Wilderness

Chassahowitzka Wilderness Saint Marks Wilderness

Breton Wilderness

700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800

LCC (RPO Projection) Easting, km -1600 -1500 -1400 -1300 -1200 -1100 -1000

-900 -800CC o t g,



# Figure 1. Proposed Modeling Domain

– Statoil OCS

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 19 

# Appendix A

Proposed CALPUFF Control Parameters

## 20 Statoil OCS EGOM, Transocean Discoverer CALPUFF run, 18 6x108x4km mesh, VISTAS 2003 4km domain 2.

VISTAS CALMET re

- processed by USFS (Tim Allen) . Hourly ozone, 3 ppb NH3.

No gridded receptors, only Class I receptors and Class II "25 mile" receptors.

---------------- Run title (3 lines)

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

CALPUFF MODEL CONTROL FILE

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

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

INPUT GROUP: 0

-- Input and Output File Names

------------- - Default Name Type File Name

------------ ---- --------- CALMET.DAT input * METDAT = *

or

ISCMET.DAT input * ISCDAT = *

or

PLMMET.DAT input * PLMDAT = *

or

PROFILE.DAT input * PRFDAT = *

SURFACE.DAT input * SFCDAT = *

RESTARTB.DAT input * RSTARTB= *

-------------------------------------------------------------------------------- CALPUFF.LST output ! PUFLST = statoil.20 03.out !

CONC.DAT output ! CONDAT = statoil.2003.con !

DFLX.DAT output ! DFDAT = statoil.2003.dry !

WFLX.DAT output ! WFDAT = statoil.2003.wet !

VISB.DAT output ! VISDAT = statoil.2003.vis !

TK2D.DAT output * T2DD AT = *

RHO2D.DAT output * RHODAT = *

RESTARTE.DAT output * RSTARTE= *

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

Emission Files

-------------- PTEMARB.DAT input * PTDAT = *

VOLEMARB.DAT input * VOLDAT = *

BAEMARB.DAT input * ARDAT = *

LNEMARB.DAT input * LNDAT = *

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

Other Files

- ---------- OZONE.DAT input ! OZDAT = ozone/ozone.2001

- 2003.dat !

VD.DAT input * VDDAT = *

CHEM.DAT input * CHEMDAT= *

H2O2.DAT input * H2O2DAT= *

HILL.DAT input * HILDAT= *

HILLRCT.DAT input * RCTDAT= *

COASTLN.DAT input * CSTDAT= *

FLUXBDY.DAT input * BDYDAT= *

BCON.DAT input * BCNDAT= *

DEBUG.DAT output * DEBUG = *

MASSFL X.DAT output * FLXDAT= *

MASSBAL.DAT output ! BALDAT= statoil.2003.bal !

FOG.DAT output * FOGDAT= *

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

All file names will be converted to lower case if LCFILES = T

Otherwise, if LCFILES = F, file names will be converted to UPPER CASE

T = lower case ! LCFILES = T !

F = UPPER CASE

NOTE: (1) file/path names can be up to 70 characters in length

Provision for multiple input files

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

Number of CALMET.DAT files for run (NMETDAT)

Default: 1 ! NMETDAT = 13 !

Number of PTEMARB.DAT files for run (NPTDAT)

Default: 0 ! NPTDAT = 0 !

Number of BAEMARB.DAT files for run (NARDAT)

Default: 0 ! NARDAT = 0 !

Number of VOLEMARB.DAT files for run (NVOLDAT)

Default: 0 ! NVOLDAT = 0 !

!END!

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 21 -------------

Subgroup (0a)

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

The following CALMET.DAT filenames are processed in sequence if NMETDAT>1

Default Name Type File Name

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

! METDAT = /home/data/calm et/vistas/2002/subdomain2_4km/met2002

- dom2

- dec.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- jan.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- feb.dat ! !END!

! METDAT = /home/data/calme t/vistas/2003/subdomain2_4km/met2003

- dom2

- mar.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- apr.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- may.dat ! !END!

! METDAT = /home/data/calmet /vistas/2003/subdomain2_4km/met2003

- dom2

- jun.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- jul.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- aug.dat ! !END!

! METDAT = /home/data/calmet/ vistas/2003/subdomain2_4km/met2003

- dom2

- sep.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- oct.dat ! !END!

! METDAT = /home/data/calmet/vistas/2003/subdomain2_4km/met2003

- dom2

- nov.dat ! !END!

! METDAT = /home/data/calmet/v istas/2003/subdomain2_4km/met2003

- dom2

- dec.dat ! !END!

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

INPUT GROUP: 1

-- General run control parameters

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

Option to run all periods found

in the met. f ile (METRUN) Default: 0 ! METRUN = 0 !

METRUN = 0

- Run period explicitly defined below

METRUN = 1

- Run all periods in met. file

Starting date: Year (IBYR)

-- No default ! IBYR = 2002 ! * start 2 days early *

(used only if Month (IBMO)

-- No default ! IBMO = 12 !

METRUN = 0) Day (IBDY)

-- No default ! IBDY = 30 !

Hour (IBHR)

-- No default ! IBHR = 01 !

Base time zone (XBTZ)

-- No default ! XBTZ = 5.0 !

PST = 8., MST = 7.

CST = 6., EST = 5.

Length of run (hours) (IRLG)

-- No default ! IRLG = 8784 ! * 364+2 days*

Number of chemical species (NSPEC)

Defaul t: 5 ! NSPEC = 9 !

Number of chemical species

to be emitted (NSE) Default: 3 ! NSE = 9 !

Flag to stop run after

SETUP phase (ITEST) Default: 2 ! ITEST = 2 !

(Used to allow chec king

of the model inputs, files, etc.)

ITEST = 1

- STOPS program after SETUP phase

ITEST = 2

- Continues with execution of program

after SETUP

Restart Configuration:

Control flag (MRESTART) Default: 0 ! MRESTART = 0 !

0 = Do not read or write a restart file

1 = Read a restart file at the beginning of

the run

2 = Write a restart file during run

3 = Read a restart file at beginning of run

and write a restart file during run

Number of periods in Restart

output cycle (NRESPD) Default: 0 ! NRESPD = 0 !

0 = File written only at last period

>0 = File updated e very NRESPD periods

Meteorological Data Format (METFM)

Default: 1 ! METFM = 1 !

METFM = 1

- CALMET binary file (CALMET.MET)

METFM = 2

- ISC ASCII file (ISCMET.MET)

## 22 METFM = 3

- AUSPLUME ASCII file (PLMMET.MET)

METFM = 4

- CTDM plus tower file (PROFILE.DAT) and

surface parameters file (SURFACE.DAT)

METFM = 5

- AERMET tower file (PROFILE.DAT) and

surface para meters file (SURFACE.DAT)

Meteorological Profile Data Format (MPRFFM)

(used only for METFM = 1, 2, 3)

Default: 1 ! MPRFFM = 1 !

MPRFFM = 1

- CTDM plus tower file (PROFILE.DAT)

MPRFFM = 2

- AERMET tower file (PROFILE.DAT)

PG sigma

- y is adjusted by the factor (AVET/PGTIME)**0.2

Averaging Time (minutes) (AVET)

Default: 60.0 ! AVET = 60. !

PG Averaging Time (minutes) (PG TIME)

Default: 60.0 ! PGTIME = 60. !

!END!

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

INPUT GROUP: 2

-- Technical options

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

Vertical distribution used in the

near field (MGAUSS) Default: 1 ! MGAUSS = 1 !

0 = uniform

1 = Gaussian

Terrain adjustment method

(MCTADJ) Default: 3 ! MCTADJ = 3 !

0 = no adjustment

1 = ISC

- type of terrain adjustment

2 = simple, CALPUFF

- type of terrain

adjustment

3 = partial plume path adjustment

Subgrid

- scale complex terrain

flag (MCTSG) Default: 0 ! MCTSG = 0 !

0 = not modeled

1 = modeled

Near

- field puffs modeled as

elongated 0 (MSLUG) Default: 0 ! MSLUG = 0 !

0 = no

1 = yes (slug model used)

Transitional plume rise modeled ?

(MTRANS ) Default: 1 ! MTRANS = 1 !

0 = no (i.e., final rise only)

1 = yes (i.e., transitional rise computed)

Stack tip downwash? (MTIP) Default: 1 ! MTIP = 1 !

0 = no (i.e., no s tack tip downwash)

1 = yes (i.e., use stack tip downwash)

Method used to simulate building

downwash? (MBDW) Default: 1 ! MBDW = 2 !

1 = ISC method

2 = PRIME method

Vertical wind shear mod eled above

stack top? (MSHEAR) Default: 0 ! MSHEAR = 0 !

0 = no (i.e., vertical wind shear not modeled)

1 = yes (i.e., vertical wind shear modeled)

Puff splitting allowed? (MSPLIT) Default: 0 ! M SPLIT = 0 !

0 = no (i.e., puffs not split)

1 = yes (i.e., puffs are split)

Chemical mechanism flag (MCHEM) Default: 1 ! MCHEM = 1 !

0 = chemical transformation not modeled

1 = transformation rates comput ed internally (MESOPUFF II scheme)

2 = user

- specified transformation rates used

3 = transformation rates computed internally (RIVAD/ARM3 scheme)

4 = secondary organic aerosol formation computed (MESOPUFF II scheme for OH)

Aque ous phase transformation flag (MAQCHEM)

(Used only if MCHEM = 1, or 3) Default: 0 ! MAQCHEM = 0 !

0 = aqueous phase transformation not modeled

1 = transformation rates adjusted for aqueous phase reactions

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 23 Wet remova l modeled ? (MWET) Default: 1 ! MWET = 1 !

0 = no

1 = yes

Dry deposition modeled ? (MDRY) Default: 1 ! MDRY = 1 !

0 = no

1 = yes

(dry deposition method specified for each species in Input Group 3)

Gravitational settling (plume tilt)

modeled ? (MTILT) Default: 0 ! MTILT = 0 !

0 = no

1 = yes

(puff center falls at the gravitational

settling velocity for 1 particle sp ecies)

Restrictions:

- MDRY = 1

- NSPEC = 1 (must be particle species as well)

- sg = 0 GEOMETRIC STANDARD DEVIATION in Group 8 is

set to zero for a single particle diameter

Method u sed to compute dispersion

coefficients (MDISP) Default: 3 ! MDISP = 3 !

1 = dispersion coefficients computed from measured values

of turbulence, sigma v, sigma w

2 = dispersion coefficients from int ernally calculated

sigma v, sigma w using micrometeorological variables

(u*, w*, L, etc.)

3 = PG dispersion coefficients for RURAL areas (computed using

the ISCST multi

- segment approximation) and MP coefficients in

urban areas

4 = same as 3 except PG coefficients computed using

the MESOPUFF II eqns.

5 = CTDM sigmas used for stable and neutral conditions.

For unstable conditions, sigmas are computed as in

MDISP = 3, described above. MDISP = 5 assumes that

measured values are read

Sigma

- v/sigma

- theta, sigma

- w measurements used? (MTURBVW)

(Used only if MDISP = 1 or 5) Default: 3 ! MTURBVW = 3 !

1 = use sigma

- v or sigma

- theta measurements

from PROFILE.DAT to compute sigma

- y

(valid for METFM = 1, 2, 3, 4, 5)

2 = use sigma

- w measurements

from PROFILE.DAT to compute sigma

- z

(valid for METFM = 1, 2, 3, 4, 5)

3 = use both sigma

- (v/theta) and sigma

- w

from PROFILE.DAT to compute sigma

- y and sigma

- z

(valid for METFM = 1, 2, 3, 4, 5)

4 = use sigma

- theta measurements

from PLMMET.DAT to compute sigma

- y

(v alid only if METFM = 3)

Back

- up method used to compute dispersion

when measured turbulence data are

missing (MDISP2) Default: 3 ! MDISP2 = 3 !

(used only if MDISP = 1 or 5)

2 = dispersion coefficients from internally calculated

sigma v, sigma w using micrometeorological variables

(u*, w*, L, etc.)

3 = PG dispersion coefficients for RURAL areas (computed using

the ISCST multi

- segment approximation) and MP coe fficients in

urban areas

4 = same as 3 except PG coefficients computed using

the MESOPUFF II eqns.

[DIAGNOSTIC FEATURE]

Method used for Lagrangian timescale for Sigma

- y

(used only if MDISP=1,2 or MDISP2=1,2)

(MTAULY) Default: 0 ! MTAULY = 0 !

0 = Draxler default 617.284 (s)

1 = Computed as Lag. Length / (.75 q)

-- after SCIPUFF

10 < Direct user input (s)

-- e.g., 306.9

[DIAGNOSTI C FEATURE]

Method used for Advective

- Decay timescale for Turbulence

(used only if MDISP=2 or MDISP2=2)

## 24 (MTAUADV) Default: 0 ! MTAUADV = 0 !

0 = No turbulence advection

1 = Computed (OPTION N OT IMPLEMENTED)

10 < Direct user input (s)

-- e.g., 800

Method used to compute turbulence sigma

- v &

sigma

- w using micrometeorological variables

(Used only if MDISP = 2 or MDISP2 = 2)

(MCTURB) Defau lt: 1 ! MCTURB = 1 !

1 = Standard CALPUFF subroutines

2 = AERMOD subroutines

PG sigma

- y,z adj. for roughness? Default: 0 ! MROUGH = 0 !

(MROUGH)

0 = no

1 = yes

Partial plume penetration of Default: 1 ! MPARTL = 1 !

elevated inversion?

(MPARTL)

0 = no

1 = yes

Strength of temperature inversion Default: 0 ! MTINV = 0 !

provided in PROFILE.DAT extended records?

(MTINV)

0 = no (computed from measured/default gradients)

1 = yes

PDF used for dispersion under convective conditions?

Default: 0 ! MPDF = 0 !

(MPDF)

0 = no

1 = yes

Sub

- Grid TI BL module used for shore line?

Default: 0 ! MSGTIBL = 0 !

(MSGTIBL)

0 = no

1 = yes

Boundary conditions (concentration) modeled?

Default: 0 ! MBCON = 0 !

(MBCON)

0 = no

1 = yes, using formatted BCON.DAT file

2 = yes, using unformatted CONC.DAT file

Note: MBCON > 0 requires that the last species modeled

be 'BCON'. Mass is placed in species BCON when

generating boundary condition puffs so that clean

air entering the modeling domain can be simulated

in the same way as polluted air. Specify zero

emission of species BCON for all regular sources.

Individual source contributions saved?

Default: 0 ! MSOURCE = 0 !

(MSOURCE)

0 = no

1 = yes

Analyses of fogging and icing impacts due to emissions from

arrays of mechanical ly

- forced cooling towers can be performed

using CALPUFF in conjunction with a cooling tower emissions

processor (CTEMISS) and its associated postprocessors. Hourly

emissions of water vapor and temperature from each cooling tower

cell a re computed for the current cell configuration and ambient

conditions by CTEMISS. CALPUFF models the dispersion of these

emissions and provides cloud information in a specialized format

for further analysis. Output to FOG.DAT is provided in either

'plume mode' or 'receptor mode' format.

Configure for FOG Model output?

Default: 0 ! MFOG = 0 !

(MFOG)

0 = no

1 = yes

- report results in PLUME Mode format

2 = yes

- report results in RECEPTOR Mode format

Test options specified to see if

they conform to regulatory

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 25 values? (MREG) Default: 1 ! MREG = 1 !

0 = NO checks are made

1 = Technical options m ust conform to USEPA

Long Range Transport (LRT) guidance

METFM 1 or 2

AVET 60. (min)

PGTIME 60. (min)

MGAUSS 1

MCTADJ 3

MTRANS 1

MTIP 1

MCHEM 1 or 3 (if modeling SOx, NOx)

MWET 1

MDRY 1

MDISP 2 or 3

MPDF 0 if MDISP=3

1 if MDISP=2

MROUGH 0

MPARTL 1

SYTDEP 550. (m)

MHFTSZ 0

SVMIN 0.5 (m /s)

!END!

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

INPUT GROUP: 3a, 3b

-- Species list

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

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

Subgroup (3a)

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

The following species are modeled:

! CSPEC = SO2 ! !END!

! CSPEC = SO4 ! !END!

! CSPEC = NOX ! !END!

! CSPEC = HNO3 ! !END!

! CSPEC = NO3 ! !END!

! CSPEC = EC ! !END!

! CSPEC = SOA ! !END!

! CSPEC = PMC ! !END!

! CSPEC = PMF ! !END!

Dry OUTPUT GROUP

SPECIES MODELED EMITTED DEPOSITED NUMBER

NAME (0=NO, 1=YES) (0=NO, 1=YES) (0=NO, (0=NONE,

(Limit: 12 1=COMPUTED

- GAS 1=1st CGRUP,

Characters 2=COMPUTED

- PARTICLE 2=2nd CGRUP,

in length ) 3=USER

- SPECIFIED) 3= etc.)

! SO2 = 1, 1, 1, 0 !

! SO4 = 1, 1, 2, 0 !

! NOX = 1, 1, 1, 0 !

! HNO3 = 1, 0, 1, 0 !

! NO3 = 1, 1, 2, 0 !

! EC = 1, 1, 2, 0 !

! SOA = 1, 1, 2, 0 !

! PMC = 1, 1, 2, 0 !

! PMF = 1, 1, 2, 0 !

!END!

Note: The last species in (3a) must be 'BCON' when using the

boundary condition option (MBCON > 0). Species BCON should

typically be modeled as inert (no chem transformation or

removal).

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

Subgroup (3b)

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

The following names are used for Species

- Groups in which results

for certain species are combined (added) prior to output. The

CGRUP name will be used as the species name in output files.

Use this feature to mod el specific particle

- size distributions

by treating each size

- range as a separate species.

Order must be consistent with 3(a) above.

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

## 26 INPUT GROUP: 4

-- Map Projection and Grid control parameters

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

Projection for all (X,Y):

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

Map projection

(PMAP) Default: UTM ! PMAP = LCC !

UTM : Universal Transverse Mercator

TTM : Tangen tial Transverse Mercator

LCC : Lambert Conformal Conic

PS : Polar Stereographic

EM : Equatorial Mercator

LAZA : Lambert Azimuthal Equal Area

False Easting and Northing (km) at the projection origin

(Used only if PMAP= TTM, LCC, or LAZA)

(FEAST) Default=0.0 ! FEAST = 0.000 !

(FNORTH) Default=0.0 ! FNORTH = 0.000 !

UTM zone (1 to 60)

(Used only if PMAP=UTM)

(IUTMZN) N o Default ! IUTMZN = 0 !

Hemisphere for UTM projection?

(Used only if PMAP=UTM)

(UTMHEM) Default: N ! UTMHEM = N !

N : Northern hemisphere projection

S : Southern hemisphere projectio n

Latitude and Longitude (decimal degrees) of projection origin

(Used only if PMAP= TTM, LCC, PS, EM, or LAZA)

(RLAT0) No Default ! RLAT0 = 40.0N !

(RLON0) No Default ! RLON0 = 97.0W !

TTM : RLON0 identifies central (true N/S) meridian of projection

RLAT0 selected for convenience

LCC : RLON0 identifies central (true N/S) meridian of projection

RLAT0 selected for convenience

PS : RLON0 identifies central (grid N/S) meridian of projection

RLAT0 selected for convenience

EM : RLON0 identifies central meridian of projection

RLAT0 is REPLACED by 0.0N (Equator)

LAZA: RLON0 identifie s longitude of tangent

- point of mapping plane

RLAT0 identifies latitude of tangent

- point of mapping plane

Matching parallel(s) of latitude (decimal degrees) for projection

(Used only if PMAP= LCC or PS)

(XLAT1) No Default ! XLAT1 = 33.0N !

(XLAT2) No Default ! XLAT2 = 45.0N !

LCC : Projection cone slices through Earth's surface at XLAT1 and XLAT2

PS : Projection plane slices through Earth at XLAT1

(XLAT2 is not used)

---------- Note: Latitudes and longitudes should be positive, and include a

letter N,S,E, or W indicating north or south latitude, and

east or west longitude. For example,

35. 9 N Latitude = 35.9N

118.7 E Longitude = 118.7E

Datum

- region

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

The Datum

- Region for the coordinates is identified by a character

string. Many mapping products currently available use the model of the

Ea rth known as the World Geodetic System 1984 (WGS

- 84). Other local

models may be in use, and their selection in CALMET will make its output

consistent with local mapping products. The list of Datum

- Regions with

official transformation param eters is provided by the National Imagery and

Mapping Agency (NIMA).

NIMA Datum

- Regions(Examples)

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

WGS

- 84 WGS

- 84 Reference Ellipsoid and Geoid, Global c overage (WGS84)

NAS

- C NORTH AMERICAN 1927 Clarke 1866 Spheroid, MEAN FOR CONUS (NAD27)

NAR

- C NORTH AMERICAN 1983 GRS 80 Spheroid, MEAN FOR CONUS (NAD83)

NWS

- 84 NWS 6370KM Radius, Sphere

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 27 ESR

- S ESRI REFERENCE 6371KM Radius, Sphere

Datum

- region for output coordinates

(DATUM) Default: WGS

- G ! DATUM = NWS

- 84 !

*** Same as MM5 ***

METEOROLOGICAL Grid:

Rectangular grid defined for project ion PMAP,

with X the Easting and Y the Northing coordinate

No. X grid cells (NX) No default ! NX = 263 !

No. Y grid cells (NY) No default ! NY = 206 !

No. vertical layers (NZ) No default ! NZ = 10 !

Grid spacing (DGRIDKM) No default ! DGRIDKM = 4 !

Units: km

Cell face heights

(ZFACE(nz+1)) No defaults

Units: m

! ZFACE =0.,20.,40.,80.,160.,320.,640.,1200.,2000.,3000.,4000. !

Reference Coordinates

of SOUTHWEST corner of

grid cell(1, 1):

X coordinate (XORIGKM) No default ! XORIGKM = 721.995 !

Y coordinate (YORIGKM) No default ! YORIGKM =

- 1598.000 !

Units: km

COMPUTATIONAL Grid:

The computational grid is identical to or a subset of the MET. grid.

The lower lef t (LL) corner of the computational grid is at grid point

(IBCOMP, JBCOMP) of the MET. grid. The upper right (UR) corner of the

computational grid is at grid point (IECOMP, JECOMP) of the MET. grid.

The grid spacing of the computational grid is the same as the MET. grid.

X index of LL corner (IBCOMP) No default ! IBCOMP = 1 !

(1 <= IBCOMP <= NX)

Y index of LL corner (JBCOMP) No default ! JBCOMP = 56 !

(1 <= JBCOMP <= N Y)

X index of UR corner (IECOMP) No default ! IECOMP = 186 !

(1 <= IECOMP <= NX)

Y index of UR corner (JECOMP) No default ! JECOMP = 163 !

(1 <= JECOMP <= NY)

SAMPLING Grid (GRIDDED R ECEPTORS):

The lower left (LL) corner of the sampling grid is at grid point

(IBSAMP, JBSAMP) of the MET. grid. The upper right (UR) corner of the

sampling grid is at grid point (IESAMP, JESAMP) of the MET. grid.

The sampling grid must be identical to or a subset of the computational

grid. It may be a nested grid inside the computational grid.

The grid spacing of the sampling grid is DGRIDKM/MESHDN.

Logical flag indicating if gridded

receptors are used (LSAMP) Default: T ! LSAMP = F !

(T=yes, F=no)

X index of LL corner (IBSAMP) No default ! IBSAMP = 1 !

(IBCOMP <= IBSAMP <= IECOMP)

Y index of LL corner (JBSAMP) No default ! JBSAMP = 56 !

(JBCOMP <= JBSAMP <= JECOMP)

X index of UR corner (IESAMP) No default ! IESAMP = 186 !

(IBCOMP <= IESAMP <= IECOMP)

Y index of UR corner (JESAMP) No default ! JESAMP = 163 !

(JBCOMP <= JESAMP <= JECOMP )

Nesting factor of the sampling

grid (MESHDN) Default: 1 ! MESHDN = 1 !

(MESHDN is an integer >= 1)

!END!

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

## 28 INPUT GROUP: 5

-- Output Options

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

* *

FILE DEFAULT VALUE VALUE THIS RUN

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

C oncentrations (ICON) 1 ! ICON = 1 !

Dry Fluxes (IDRY) 1 ! IDRY = 1 !

Wet Fluxes (IWET) 1 ! IWET = 1 !

2D Temperature (IT2D) 0 ! IT2D = 0 !

2D Density (IRHO) 0 ! IRHO = 0 !

Relative Humidity (IVIS) 1 ! IVIS = 1 !

(RH file is required for

VISIBILITY analyses)

Use data compression option in output file?

(LCOMPRS) Default: T ! LCOMPRS = T !

*

0 = Do not create file, 1 = create file

QA PLOT FILE OUTPUT OPTION:

Create a standard series of output files (e.g.

l ocations of sources, receptors, grids ...)

suitable for plotting?

(IQAPLOT) Default: 1 ! IQAPLOT = 1 !

0 = no

1 = yes

DIAGNOSTIC MASS FLUX OUTPUT OPTIONS:

Mass flux across specifi ed boundaries

for selected species reported hourly?

(IMFLX) Default: 0 ! IMFLX = 0 !

0 = no

1 = yes (FLUXBDY.DAT and MASSFLX.DAT filenames

are specified in Input Group 0)

Mass balance for each species

reported hourly?

(IMBAL) Default: 0 ! IMBAL = 1 !

0 = no

1 = yes (MASSBAL.DAT filename is

specified in Input Group 0)

LINE PRINTER OUT PUT OPTIONS:

Print concentrations (ICPRT) Default: 0 ! ICPRT = 0 !

Print dry fluxes (IDPRT) Default: 0 ! IDPRT = 0 !

Print wet fluxes (IWPRT) Default: 0 ! IWPRT = 0 !

(0 = Do no t print, 1 = Print)

Concentration print interval

(ICFRQ) in hours Default: 1 ! ICFRQ = 24 !

Dry flux print interval

(IDFRQ) in hours Default: 1 ! IDFRQ = 1 !

Wet flux p rint interval

(IWFRQ) in hours Default: 1 ! IWFRQ = 1 !

Units for Line Printer Output

(IPRTU) Default: 1 ! IPRTU = 1 !

for for

Concentration Deposition

1 = g/m**3 g/m**2/s

2 = mg/m**3 mg/m**2/s

3 = ug/m**3 ug/m**2/s

4 = ng/m**3 ng/m**2/s

5 = Odour Units

Mes sages tracking progress of run

written to the screen ?

(IMESG) Default: 2 ! IMESG = 2 !

0 = no

1 = yes (advection step, puff ID)

2 = yes (2003JJJHH, # old puffs, # emitted puffs)

SPECIES (or GROUP for combined species) LIST FOR OUTPUT OPTIONS

---- CONCENTRATIONS

---- ------ DRY FLUXES

------ ------ WET FLUXES

------ -- MASS FLUX

--

SPECIES

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 29 /GROUP PRINTED? SAVED ON DISK? PRINTED? SAVED ON DISK? PRINTED? SAVED ON DISK? SAVED ON DISK?

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

! SO2 = 0, 1, 0, 1, 0, 1, 0 !

! NOX = 0, 1, 0, 1, 0, 1, 0 !

! HNO3 = 0, 1, 0, 1, 0, 1, 0 !

! SO4 = 0, 1, 0, 1, 0, 1, 0 !

! NO3 = 0, 1, 0, 1, 0, 1, 0 !

! EC = 0, 1, 0, 1, 0, 1, 0 !

! SOA = 0, 1, 0, 1, 0, 1, 0 !

! PMC = 0, 1, 0, 1, 0, 1, 0 !

! PMF = 0, 1, 0, 1, 0, 1, 0 !

Not e: Species BCON (for MBCON > 0) does not need to be saved on disk.

OPTIONS FOR PRINTING "DEBUG" QUANTITIES (much output)

Logical for debug output

(LDEBUG) Default: F ! LDEBUG = F !

First puff to track

(IPFDEB) Default: 1 ! IPFDEB = 1 !

Number of puffs to track

(NPFDEB) Default: 1 ! NPFDEB = 1 !

Met. period to start output

(NN1) Default: 1 ! NN1 = 1 !

Met. period to end output

(NN2) Default: 10 ! NN2 = 10 !

!END!

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

INPUT GROUP: 6a, 6b, & 6c

-- Subgrid scale complex terrain inputs

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

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

Subgroup (6a)

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

Number of terrain features (NHILL) Default: 0 ! NHILL = 0 !

Number of special com plex terrain

receptors (NCTREC) Default: 0 ! NCTREC = 0 !

Terrain and CTSG Receptor data for

CTSG hills input in CTDM format ?

(MHILL) No Default ! MHILL = 2 !

1 = Hill and Receptor data created

by CTDM processors & read from

HILL.DAT and HILLRCT.DAT files

2 = Hill data created by OPTHILL &

input below in Subgroup (6b);

Receptor data in Subgroup (6c)

Factor to convert horizontal dimensions Default: 1.0 ! XHILL2M = 0. !

to meters (MHILL=1)

Factor to convert vertical dimensions Default: 1.0 ! ZHILL2M = 0. !

to meters (MHILL=1)

X

- origin of CTDM system relative t o No Default ! XCTDMKM = 0.0E00 !

CALPUFF coordinate system, in Kilometers (MHILL=1)

Y

- origin of CTDM system relative to No Default ! YCTDMKM = 0.0E00 !

CALPUFF coordinate system, in Kilometers (MHILL=1)

! END !

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

Subgroup (6b)

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

1 **

HILL information

HILL XC YC THETAH ZGRID RELIEF EXPO 1 EXPO 2 SCALE 1 SCALE 2 AMAX1

AMAX2

NO. (km) (km) (deg.) (m) (m) (m) (m) (m) (m) (m)

(m)

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

----

## 30 ---------------

Subgroup (6c)

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

COMPLEX TERRAIN REC EPTOR INFORMATION

XRCT YRCT ZRCT XHH

(km) (km) (m)

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

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

1

Description of Complex Terr ain Variables:

XC, YC = Coordinates of center of hill

THETAH = Orientation of major axis of hill (clockwise from

North)

ZGRID = Height of the 0 of the grid above mean sea

level

RELIEF = Height of the crest of the hill above the grid elevation

EXPO 1 = Hill

- shape exponent for the major axis

EXPO 2 = Hill

- shape exponent for the major axis

SCALE 1 = Horizontal length scale along the major axis

SCALE 2 = Horizontal length scale along the minor axis

AMAX = Maximum allowed axis length for the major axis

BMAX = Maximum allowed axis length for the major axis

XRCT, YRCT = Coordinates of the complex terra in receptors

ZRCT = Height of the ground (MSL) at the complex terrain

Receptor

XHH = Hill number associated with each complex terrain receptor

(NOTE: MUST BE ENTERED AS A REAL NUMBER)

* *

NOTE: DATA for each hill and CTSG receptor are treated as a separate

input subgroup and therefore must end with an input group terminator.

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

INPUT GROUP: 7

-- Chemical parameters for dry deposition of gases

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

SPECIES DIFFUSIVITY ALPHA STAR REACTIVITY MESOPHYLL RESISTANCE HENRY'S LAW

COEFFICIENT

NAME (cm**2/s) (s/cm) (dimensionless)

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

--

! SO2 = 0.1509, 1000., 8., 0., 0.04 !

! NOX = 0.1 656, 1., 8., 5., 3.5 !

! HNO3 = 0.1628, 1., 18., 0., 0.00000008 !

!END!

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

INPUT GROUP: 8

-- S ize parameters for dry deposition of particles

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

For SINGLE SPECIES, the mean and standard deviation are used to

compute a deposition velocity for NINT (see group 9) size

- ranges,

and these are then averaged to obtain a mean dep osition velocity.

For GROUPED SPECIES, the size distribution should be explicitly

specified (by the 'species' in the group), and the standard deviation

for each should be entered as 0. The model will then use the

deposition velocity f or the stated mean diameter.

SPECIES GEOMETRIC MASS MEAN GEOMETRIC STANDARD

NAME DIAMETER DEVIATION

(microns) (microns)

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

! SO4 = 0.48, 2.0 !

! NO3 = 0.48, 2.0 !

! EC = 0.48, 2.0 !

! SOA = 0.48, 2 .0 !

! PMC = 7.0, 2.0 !

! PMF = 0.48, 2.0 !

!END!

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

INPUT GROUP: 9

-- Miscellaneous dry depositio n parameters

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

Reference cuticle resistance (s/cm)

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 31 (RCUTR) Default: 30 ! RCUTR = 30.0 !

Reference ground resistance (s/cm)

(RGR) Default: 10 ! RGR = 10.0 !

Reference pollutant reactivity

(REACTR) Default: 8 ! REACTR = 8.0 !

Number of particle

- size intervals used to

evaluate effective particle deposition velocity

(NINT) Default: 9 ! NINT = 9 !

Vegetation state in unirrigated areas

(IVEG) Default: 1 ! IVEG = 1 !

IVEG=1 for active and unstressed vegetation

IVEG=2 for active and stressed vegetation

IVEG=3 for inactive vegetation

!END!

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

INPUT GROUP: 10

-- Wet Deposition Parameters

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

Sc avenging Coefficient

-- Units: (sec)**(

- 1)

Pollutant Liquid Precip. Frozen Precip.

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

! SO2 = 3.0E

- 05, 0.0E00 !

! SO4 = 1.0E

- 04, 3.0E

- 05 !

! NOX = 0.0E00, 0.0E00 !

! HNO3 = 6.0E

- 05, 0.0E00 !

! NO3 = 1.0E

- 04, 3.0E

- 05 !

! EC = 1.0E

- 04, 3.0E

- 05 !

! SOA = 1.0E

- 04, 3.0E

- 05 !

! PMC = 1.0E

- 04, 3.0E

- 05 !

! PMF = 1.0E

- 04, 3.0E

- 05 !

!END!

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

INPUT GROUP: 11

- - Chemistry Parameters

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

Ozone data input option (MOZ) Default: 1 ! MOZ = 1 !

(Used only if MCHEM = 1, 3, or 4)

0 = use a monthly background ozone value

1 = read hourly ozone concentrations from

the OZONE.DAT data file

Monthly ozone concentrations

(Used only if MCHEM = 1, 3, or 4 and

MOZ = 0 or MOZ = 1 and all hourly O3 data missing)

(BCKO3) in ppb Default: 12*80. ! BCKO3 = 12*65. !

M onthly ammonia concentrations

(Used only if MCHEM = 1, or 3)

(BCKNH3) in ppb Default: 12*10. ! BCKNH3 = 12*3.0 !

Nighttime SO2 loss rate (RNITE1)

in percent/hour Default: 0.2 ! RNITE1 = .2 !

Nighttime NOx loss rate (RNITE2)

in percent/hour Default: 2.0 ! RNITE2 = 2.0 !

Nighttime HNO3 formation rate (RNITE3)

in percent/hour Default: 2.0 ! RNITE3 = 2.0 !

H2 O2 data input option (MH2O2) Default: 1 ! MH2O2 = 0 !

(Used only if MAQCHEM = 1)

0 = use a monthly background H2O2 value

1 = read hourly H2O2 concentrations from

the H2O2.DAT data file

Monthly H2O2 co ncentrations

(Used only if MQACHEM = 1 and

MH2O2 = 0 or MH2O2 = 1 and all hourly H2O2 data missing)

(BCKH2O2) in ppb Default: 12*1. ! BCKH2O2 = 12*1.0 !

--- Data for SECONDARY ORGANIC AEROSOL (SOA) Option

( used only if MCHEM = 4)

## 32 The SOA module uses monthly values of:

Fine particulate concentration in ug/m^3 (BCKPMF)

Organic fraction of fine particulate (OFRAC)

VOC / NOX ratio (after reaction) (VCNX)

to ch aracterize the air mass when computing

the formation of SOA from VOC emissions.

Typical values for several distinct air mass types are:

Month 1 2 3 4 5 6 7 8 9 10 11 12

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Clean Continental

BCKPMF 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1.

OFRAC .15 .15 .20 .20 .20 .20 .20 .20 .20 .20 .20 .15

VCNX 50. 50. 50. 50. 50. 50 . 50. 50. 50. 50. 50. 50.

Clean Marine (surface)

BCKPMF .5 .5 .5 .5 .5 .5 .5 .5 .5 .5 .5 .5

OFRAC .25 .25 .30 .30 .30 .30 .30 .30 .30 .30 .30 .25

VCNX 50. 50. 50. 50. 50. 50. 5 0. 50. 50. 50. 50. 50.

Urban

- low biogenic (controls present)

BCKPMF 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30.

OFRAC .20 .20 .25 .25 .25 .25 .25 .25 .20 .20 .20 .20

VCNX 4. 4. 4. 4 . 4. 4. 4. 4. 4. 4. 4. 4.

Urban

- high biogenic (controls present)

BCKPMF 60. 60. 60. 60. 60. 60. 60. 60. 60. 60. 60. 60.

OFRAC .25 .25 .30 .30 .30 .55 .55 .55 .35 .35 .35 .25

VCNX 15 . 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15.

Regional Plume

BCKPMF 20. 20. 20. 20. 20. 20. 20. 20. 20. 20. 20. 20.

OFRAC .20 .20 .25 .35 .25 .40 .40 .40 .30 .30 .30 .20

VCNX 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15.

Urban

- no controls present

BCKPMF 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100.

OFRAC .30 .30 .35 .35 .35 .55 .55 .55 .35 .35 .35 .30

VCNX 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.

Default: Clean Continental

! BCKPMF = 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00 !

! OFRAC = 0.15, 0.15, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0. 20 !

! VCNX = 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0 !

!END!

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

INPUT GROUP: 12

-- Misc. Dispersion and Computational Parameters

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

Horizontal size of puff (m) beyond which

time

- dependent dispersion equations (Heffter)

are used to determine sigma

- y and

sigma

- z (SYTDEP) Default: 550. ! SYTDEP = 550.0 !

Switch for using Heffte r equation for sigma z

as above (0 = Not use Heffter; 1 = use Heffter

(MHFTSZ) Default: 0 ! MHFTSZ = 0 !

Stability class used to determine plume

growth rates for puffs above the bounda ry

layer (JSUP) Default: 5 ! JSUP = 5 !

Vertical dispersion constant for stable

conditions (k1 in Eqn. 2.7

- 3) (CONK1) Default: 0.01 ! CONK1 = .01 !

Vertical dispersion constant for neutral/

unstable conditions (k2 in Eqn. 2.7

- 4)

(CONK2) Default: 0.1 ! CONK2 = .1 !

Factor for determining Transition

- point from

Schulman

- Scire to Huber

- Snyder Building Downwash

scheme (SS used for Hs < Hb + TBD * HL)

(TBD) Default: 0.5 ! TBD = .5 !

TBD < 0 ==> always use Huber

- Snyder

TBD = 1.5 ==> always use Schulman

- Scire

TBD = 0.5 ==> ISC Transition

- point

Range of land use categories for which

urban dispersion is assumed

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 33 (IURB1, IURB2) Default: 10 ! IURB1 = 10 !

19 ! IURB2 = 19 !

Site characterization parameters f or single

- point Met data files

---------

(needed for METFM = 2,3,4,5)

Land use category for modeling domain

(ILANDUIN) Default: 20 ! ILANDUIN = 20 !

Roughness length (m) for modeling domain

(Z0IN) Default: 0.25 ! Z0IN = .25 !

Leaf area index for modeling domain

(XLAIIN) Default: 3.0 ! XLAIIN = 3.0 !

Elevation above sea level (m)

(ELE VIN) Default: 0.0 ! ELEVIN = .0 !

Latitude (degrees) for met location

(XLATIN) Default:

- 999. ! XLATIN =

- 999.0 !

Longitude (degrees) for met location

(XLON IN) Default:

- 999. ! XLONIN =

- 999.0 !

Specialized information for interpreting single

- point Met data files

-----

Anemometer height (m) (Used only if METFM = 2,3)

(ANEMHT) Default: 10. ! ANEMHT = 10.0 !

Form of lateral turbulance data in PROFILE.DAT file

(Used only if METFM = 4,5 or MTURBVW = 1 or 3)

(ISIGMAV) Default: 1 ! ISIGMAV = 1 !

0 = read s igma

- theta

1 = read sigma

- v

Choice of mixing heights (Used only if METFM = 4)

(IMIXCTDM) Default: 0 ! IMIXCTDM = 0 !

0 = read PREDICTED mixing heights

1 = read OBSERVE D mixing heights

Maximum length of a slug (met. grid units)

(XMXLEN) Default: 1.0 ! XMXLEN = 1.0 !

Maximum travel distance of a puff/slug (in

grid units) during one sampling step

(XSAMLEN) Default: 1.0 ! XSAMLEN = 1.0 !

Maximum Number of slugs/puffs release from

one source during one time step

(MXNEW) Default: 99 ! MXNEW = 99 !

Maximu m Number of sampling steps for

one puff/slug during one time step

(MXSAM) Default: 99 ! MXSAM = 99 !

Number of iterations used when computing

the transport wind for a sampling s tep

that includes gradual rise (for CALMET

and PROFILE winds)

(NCOUNT) Default: 2 ! NCOUNT = 2 !

Minimum sigma y for a new puff/slug (m)

(SYMIN) De fault: 1.0 ! SYMIN = 1.0 !

Minimum sigma z for a new puff/slug (m)

(SZMIN) Default: 1.0 ! SZMIN = 1.0 !

Default minimum turbulence velocities sigma

- v and sigma

- w

for each stability class over land and over water (m/s)

(SVMIN(12) and SWMIN(12))

---------- LAND

---------- --------- WATER

----------

Stab Class : A B C D E F A B C D E F

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

Default SVMIN : .50, .50, .50, .50, .50, .50, .37, .37, .37, .37, .37, .37

Default SWMIN : .20, .12, .08, .06, .03, .016, .20, .12, .08, .06, .03, .016

* SVMIN = 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.370, 0.370, 0.370, 0.370, 0.370, 0.370 *

! SVMIN = 12* 0.5 ! MREG = 1 requirement

## 34 ! SWMIN = 0.200, 0.120, 0.080, 0.060, 0.030, 0.016, 0.200, 0.120, 0.080, 0.060, 0.030, 0.016!

Diver gence criterion for dw/dz across puff

used to initiate adjustment for horizontal

convergence (1/s)

Partial adjustment starts at CDIV(1), and

full adjustment is reached at CDIV(2)

(CDIV(2)) Default: 0.0,0.0 ! CDIV = 0.0, 0.0 !

Minimum wind speed (m/s) allowed for

non

- calm conditions. Also used as minimum

speed returned when using power

- law

extrapolation toward surface

(WSCALM) Default: 0.5 ! WSCALM = .5 !

Maximum mixing height (m)

(XMAXZI) Default: 3000. ! XMAXZI = 3000.0 !

Minimum mixing height (m)

(XMINZI) Default: 50. ! XMINZI = 50.0 !

Default wind speed classes

--

5 upper bounds (m/s) are entered;

the 6th class has no upper limit

(WSCAT(5)) Default :

ISC RURAL : 1.54, 3. 09, 5.14, 8.23, 10.8 (10.8+)

Wind Speed Class : 1 2 3 4 5

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

! WSCAT = 1.54, 3.09, 5.14, 8.2 3, 10.80 !

Default wind speed profile power

- law

exponents for stabilities 1

- 6

(PLX0(6)) Default : ISC RURAL values

ISC RURAL : .07, .07, .10, .15, .35, .55

ISC URBAN : .15, .15, .20, .25, .30, .30

Stability Class : A B C D E F

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

! PLX0 = 0.07, 0.07, 0.10, 0.15, 0.35, 0.55 !

Default potential temperature gradient

for stable classes E, F (degK/m)

(PTG0(2)) Default: 0.020, 0.035

! PTG0 = 0.020, 0.035 !

Default plume path coefficients for

each stability class (used when option

for partial plume height terrain adjustment

is selected

-- MCTADJ=3)

(PPC(6)) Stability Class : A B C D E F

Default PPC : .50, .50, .50, .50, .35, .35

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

! PPC = 0.50, 0.50, 0.50, 0.50, 0.35, 0.35 !

Slug

- to

- puff transition criterion factor

equal to sigma

- y/length of slug

(SL2PF) Default: 10. ! SL2PF = 10.0!

Puff

- splitting control variables

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

VERTICAL SPLIT

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

Number of puffs that result every time a puff

is split

- nsplit=2 means that 1 puff splits

into 2

(NSPLIT) Default: 3 ! NSPLIT = 3 !

Time(s) of a day when split puffs are elig ible to

be split once again; this is typically set once

per day, around sunset before nocturnal shear develops.

24 values: 0 is midnight (00:00) and 23 is 11 PM (23:00)

0=do not re

- split 1=eligible for re

- split

(IRESPL IT(24)) Default: Hour 17 = 1

! IRESPLIT = 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0 !

Split is allowed only if last hour's mixing

height (m) exceeds a minimum value

(ZISPLIT) Default: 100. ! ZISPLIT = 100.0 !

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 35 Split is allowed only if ratio of last hour's

mixing ht to the maximum mixing ht experienced

by the puff is less than a maximum value (this

postpones a split until a nocturnal layer de velops)

(ROLDMAX) Default: 0.25 ! ROLDMAX = 0.25 !

HORIZONTAL SPLIT

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

Number of puffs that result every time a puff

is split

- nsplith=5 means that 1 puff splits

int o 5

(NSPLITH) Default: 5 ! NSPLITH = 5 !

Minimum sigma

- y (Grid Cells Units) of puff

before it may be split

(SYSPLITH) Default: 1.0 ! SYSPLITH = 1.0 !

Minimum puff elongation rate (SYSPLITH/hr) due to

wind shear, before it may be split

(SHSPLITH) Default: 2. ! SHSPLITH = 2.0 !

Minimum concentration (g/m^3) of each

species in puff before it may be split

Enter array of NSPEC values; if a single value is

entered, it will be used for ALL species

(CNSPLITH) Default: 1.0E

- 07 ! CNSPLITH = 1.0E

- 07 !

Integration control variables

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

Fractional convergence criterion for numerical SLUG

sampling integration

(EPSSLUG) Default: 1.0e

- 04 ! EPSSLUG = 1.0E

- 04 !

Fractional convergence criterion for numerical AREA

source int egration

(EPSAREA) Default: 1.0e

- 06 ! EPSAREA = 1.0E

- 06 !

Trajectory step

- length (m) used for numerical rise

integration

(DSRISE) Default: 1.0 ! DSRISE = 1.0 !

Boundary Condition (BC) Puff control variables

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

Minimum height (m) to which BC puffs are mixed as they are emitted

(MBCON=2 ONLY). Actual height is reset to the current mixing height

at the release point if greater than this minimum.

(HTMINBC) Default: 500. ! HTMINBC = 500. !

Search radius (km) about a receptor for sampling nearest BC puff.

BC puffs are typically emitted with a spacing of one grid cell

length, so the search radius should be greater than DGRIDKM.

(RSAMPBC) Default: 10. ! RSAMPBC = 10 !

Near

- Surface depletion adjustment to concentration profile used when

sampling BC puffs?

( MDEPBC) Default: 1 ! MDEPBC = 1 !

0 = Concentration is NOT adjusted for depletion

1 = Adjust Concentration for depletion

!END!

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

INPUT GROUPS: 13a, 13b, 13c, 13d

-- Point source parameters

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

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

Subgroup (13a)

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

Number of point sources with

parameters provided below (NPT1) No default ! NPT1 = 16 !

Units used for point source

emissions below (IPTU) Default: 1 ! IPTU = 1 !

1 = g/s

2 = kg/hr

3 = lb/hr

4 = tons/yr

5 = Odour Unit * m**3/s (vol . flux of odour compound)

6 = Odour Unit * m**3/min

7 = metric tons/yr

## 36 Number of source

- species

combinations with variable

emissions scaling factors

provided below in (13d) (NSPT1) Default: 0 ! NS PT1 = 0 !

Number of point sources with

variable emission parameters

provided in external file (NPT2) No default ! NPT2 = 0 !

(If NPT2 > 0, these point

source emissions are read from

the file: PTEMARB.DAT)

!END!

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

Subgroup (13b)

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

POINT SOURCE: CONSTANT DATA

b c

Source X UTM Y UTM Stack Base Stack Exit Exit Bldg. Emission

No. Coordinate Coordinate Height Elevation Diameter Vel. Temp. Dwash Rates

(km) (km) (m) (m) (m) (m/s) (deg. K)

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

*** Max short

- term emission rates ***

X Y Ht Elev Diam Vel Temp DW SO2 SO4 NOX HNO3 NO3 EC SOA PMF

PMC

! SRCNAM = GEN1 !

! X = 941.931,

- 1186.229,28.93,0,0.81,0.001,623.15,0,2.43E

- 01,7.83E

- 03,8 .99E+00,0,4.06E

- 04,2.42E

- 01,7.34E

- 02,0,0

!

! END !

! SRCNAM = GEN2 !

! X = 941.931,

- 1186.228,28.93,0,0.81,0.001,623.15,0,2.43E

- 01,7.83E

- 03,8.99E+00,0,4.06E

- 04,2.42E

- 01,7.34E

- 02,0,0

!

! END !

! SRCNAM = GEN3 !

! X = 941.931,

- 1186.229,26.82,0,0.81,0.001,623. 15,0,2.43E

- 01,7.83E

- 03,8.99E+00,0,4.06E

- 04,2.42E

- 01,7.34E

- 02,0,0

!

! END !

! SRCNAM = GEN4 !

! X = 941.931,

- 1186.228,26.82,0,0.81,0.001,623.15,0,2.43E

- 01,7.83E

- 03,8.99E+00,0,4.06E

- 04,2.42E

- 01,7.34E

- 02,0,0

!

! END !

! SRCNAM = GEN5 !

! X = 941.931,

- 1186.300 ,28.93,0,0.81,0.001,623.15,0,0,0,0,0,0,0,0,0,0 !

! END !

! SRCNAM = GEN6 !

! X = 941.931,

- 1186.302,28.93,0,0.81,0.001,623.15,0,0,0,0,0,0,0,0,0,0 !

! END !

! SRCNAM = GEN7 !

! X = 941.931,

- 1186.300,26.82,0,0.81,0.001,623.15,0,0,0,0,0,0,0,0,0,0 !

! END !

! S RCNAM = GEN8 !

! X = 941.931,

- 1186.302,26.82,0,0.81,0.001,623.15,0,0,0,0,0,0,0,0,0,0 !

! END !

! SRCNAM = EGEN1 !

! X = 942.005,

- 1186.225,9.45,0,0.41,0.001,888.15,0,1.16E

- 01,3.75E

- 03,7.71E+00,0,1.94E

- 04,1.16E

- 01,3.52E

- 02,0,0

!

! END !

! SRCNAM = LB1 !

! X = 942.024,

- 1186.260,19.70,0,0.89,0.001,830.15,0,1.32E

- 03,2.34E

- 04,1.13E

- 01,0,1.21E

- 05,7.25E

- 03,2.19E

- 03,0,0

!

! END !

! SRCNAM = LB2 !

! X = 942.024,

- 1186.270,19.70,0,0.89,0.001,830.15,0,1.32E

- 03,2.34E

- 04,1.13E

- 01,0,1.21E

- 05,7.25E

- 03,2.19E

- 03,0,0

!

! END !

! SRCNAM = LB3 !

! X = 941.939,

- 1186.256,19.70,0,0.89,0.001,830.15,0,1.32E

- 03,2.34E

- 04,1.13E

- 01,0,1.21E

- 05,7.25E

- 03,2.19E

- 03,0,0

!

! END !

! SRCNAM = LB4 !

! X = 941.939,

- 1186.267,19.70,0,0.89,0.001,830.15,0,1.32E

- 03,2.34E

- 04,1.13E

- 01,0,1.21E

- 05,7.25E

- 03, 2.19E

- 03,0,0

!

! END !

! SRCNAM = MOB1 !

! X = 941.927,

- 1186.264,19.70,0,1.85,0.001,773.15,0,7.44E

- 03,3.23E

- 04,2.53E

- 01,0,1.67E

- 05,1.00E

- 02,3.03E

- 03,0,0

!

! END !

! SRCNAM = CMU1 !

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 37 ! X = 941.979,

- 1186.304,9.45,0,0.20,0.001,794.15,0,2.28E

- 02,4.96E

- 04,7.46E

- 01,0,2.57E

- 05,1.54E

- 02,4.65E

- 03,0,0

!

! END !

! SRCNAM = CMU2 !

! X = 941.981,

- 1186.304,9.45,0,0.20,0.001,794.15,0,2.28E

- 02,4.96E

- 04,7.46E

- 01,0,2.57E

- 05,1.54E

- 02,4.65E

- 03,0,0

!

! END !

--------

a

Data for each source are treated as a separate inp ut subgroup

and therefore must end with an input group terminator.

SRCNAM is a 12

- character name for a source

(No default)

X is an array holding the source data listed by the column headings

(No default)

SIGYZI is an array holding the initial sigma

- y and sigma

- z (m)

(Default: 0.,0.)

FMFAC is a vertical momentum flux factor (0. or 1.0) used to represent

the effect of rain

- caps or other physical configurations that

reduce momentum rise associated with the actual exit velocity.

(Default: 1.0

-- full momentum used)

ZPLTFM is the platform height (m) for sources influenced by an isolated

structure that has a significant open area be tween the surface

and the bulk of the structure, such as an offshore oil platform.

The Base Elevation is that of the surface (ground or ocean),

and the Stack Height is the release height above the Base (not

above the platform). Building heights entered in Subgroup 13c

must be those of the buildings on the platform, measured from

the platform deck. ZPLTFM is used only with MBDW=1 (ISC

downwash method) for sources wi th building downwash.

(Default: 0.0)

b

0. = No building downwash modeled

1. = Downwash modeled for buildings resting on the surface

2. = Downwash modeled for buildings raised above the surface (ZPLTFM > 0.)

NOTE: must be entered as a REAL number (i.e., with decimal point)

c

An emission rate must be entered for every pollutant modeled.

Enter emission rate of zero for secondary pollutants that are

modeled, but not emitted. Units are specified by IPTU

(e.g. 1 for g/s).

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

Subgroup (13c)

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

BUILDING DIMENSION DATA FOR SOURCES SUBJECT TO DOWNWASH

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

Source a

No. Effective building height, width, length and X/Y offset (in meters)

every 10 degrees. LENGTH, XBADJ, and YBADJ are only needed for

MBDW=2 (PRIME downwash option)

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

--------

a

Building height, width, length, and X/Y offset from the source are treated

as a separate input subgroup for each source and therefore must end with

an input group terminator.

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

Subgroup (13d)

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

a

POINT SOURCE: VARIABLE EMISSIONS DATA

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

Use this subgroup to describe temporal variations in the emission

rates given in 13b. Factors entered multiply the rates in 13b.

Skip sources here that have constant emissions. For more elaborate

variation in source parameters, use PTEMARB.DAT and NPT2 > 0.

IVARY determines the type of var iation, and is source

- specific:

(IVARY) Default: 0

0 = Constant

1 = Diurnal cycle (24 scaling factors: hours 1

- 24)

2 = Monthly cycle (12 scaling factors: months 1

- 12)

## 38 3 = Hour & Season (4 groups of 24 hourly scaling factors,

where first group is DEC

- JAN

- FEB)

4 = Speed & Stab. (6 groups of 6 scaling factors, where

firs t group is Stability Class A,

and the speed classes have upper

bounds (m/s) defined in Group 12

5 = Temperature (12 scaling factors, where temperature

classes have upper bounds (C) of:

0, 5, 10, 15, 20, 25, 30, 35, 40,

45, 50, 50+)

--------

a

Data for each species are treated as a separate input subgrou p

and therefore must end with an input group terminator.

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

INPUT GROUPS: 14a, 14b, 14c, 14d

-- Area source parameters

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

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

Subgr oup (14a)

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

Number of polygon area sources with

parameters specified below (NAR1) No default ! NAR1 = 0 !

Units used for area source

emissions below (IARU) Default: 1 ! IARU = 1 !

1 = g/m**2/s

2 = kg/m**2/hr

3 = lb/m**2/hr

4 = tons/m**2/yr

5 = Odour Unit * m/s (vol. flux/m**2 of odour compound)

6 = Odour Unit * m/min

7 = metric tons/m**2/yr

Number of source

- species

combinations with variable

emissions scaling factors

provided below in (14d) (NSAR1) Default: 0 ! NSAR1 = 0 !

Number of buoyant polygon area sources

with variable location and emission

parameters (NAR2) No default ! NAR2 = 0 !

(If NAR2 > 0, ALL parameter data for

these sources are read from the file: BAEMARB.DAT)

!END!

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

Subgroup (14b)

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

a

AREA SOURCE: CONSTANT DATA

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

b

Source Effect. Base Initial Emission

No. Height Elevation Sigma z Rates

(m) (m) (m)

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

--------

a

Data for each source are treated as a separate input subgroup

and therefore must end with an input group termi nator.

b

An emission rate must be entered for every pollutant modeled.

Enter emission rate of zero for secondary pollutants that are

modeled, but not emitted. Units are specified by IARU

(e.g. 1 for g/m**2/s).

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

Subg roup (14c)

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

COORDINATES (UTM

- km) FOR EACH VERTEX(4) OF EACH POLYGON

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

Source a

No. Ordered lis t of X followed by list of Y, grouped by source

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

--------

a

Data for each source are treated as a separate input subgroup

and therefore must end with an input group termin ator.

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 39 ---------------

Subgroup (14d)

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

a

AREA SOURCE: VARIABLE EMISSIONS DATA

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

Use this subgroup to describe temporal variations i n the emission

rates given in 14b. Factors entered multiply the rates in 14b.

Skip sources here that have constant emissions. For more elaborate

variation in source parameters, use BAEMARB.DAT and NAR2 > 0.

IVARY determines the type of variation, and is source

- specific:

(IVARY) Default: 0

0 = Constant

1 = Diurnal cycle (24 scaling factors: hours 1

- 24)

2 = Monthly cycle (12 scaling factors: months 1

- 12)

3 = Hour & Season (4 groups of 24 hourly scaling factors,

where first group is DEC

- JAN

- FEB)

4 = Speed & Stab. (6 groups of 6 scaling factors, where

first group is Stability Class A,

and the speed classes have upper

bounds (m/s) defined in Group 12

5 = Temperature (12 scaling factors, where temperature

classes have upper bounds (C) of:

0, 5, 10, 15, 20, 25, 30, 35, 40,

45, 50, 50+)

--------

a

Data for each species are treated as a separate input s ubgroup

and therefore must end with an input group terminator.

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

INPUT GROUPS: 15a, 15b, 15c

-- Line source parameters

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

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

Subgroup (15a)

--------------- Number of buoyant line sources

with variable location and emission

parameters (NLN2) No default ! NLN2 = 0 !

(If NLN2 > 0, ALL parameter data for

these sources are read from the file: LNEMARB.DAT)

Number of buoyant line sources (NLINES) No default ! NLINES = 0 !

Units used for line source

emissions below (ILNU) Default: 1 ! ILNU = 1 !

1 = g/s

2 = kg/hr

3 = lb/hr

4 = tons/yr

5 = Odour Unit * m**3/s (vol. flux of odour compound)

6 = Odour Unit * m**3/min

7 = metric tons/yr

Number of source

- species

com binations with variable

emissions scaling factors

provided below in (15c) (NSLN1) Default: 0 ! NSLN1 = 0 !

Maximum number of segments used to model

each line (MXNSEG) Default: 7 ! MXNSEG = 7 !

The following variables are required only if NLINES > 0. They are

used in the buoyant line source plume rise calculations.

Number of distances at which Default: 6 ! NLRISE = 6 !

transitional rise is computed

Average building length (XL) No default ! XL = .0 !

(in meters)

Average building height (HBL) No default ! HBL = .0 !

(in meters)

Average building width (WBL) No default ! WBL = .0 !

(in meters)

## 40 Average line source width (WML) No default ! WML = .0 !

(in meters)

Average separation between buildings (DXL) No default ! DXL = .0 !

(in meters)

Average buoyancy parameter (FPRIMEL) No default ! FPRIMEL = .0 !

(in m**4/s**3)

!END!

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

Subgroup (15b)

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

BUOYANT LINE SOURCE: CONSTANT DATA

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

a

Source Beg. X Beg. Y End. X End. Y Release Base Emission

No. Coordinate Coordinate Coordinate Coordinate Height Elevation Rates

(km) (km) (km) (km) (m) (m)

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

a

Data for each source are treated as a separate input subgroup

and there fore must end with an input group terminator.

b

An emission rate must be entered for every pollutant modeled.

Enter emission rate of zero for secondary pollutants that are

modeled, but not emitted. Units are specified by ILNTU

(e .g. 1 for g/s).

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

Subgroup (15c)

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

a

BUOYANT LINE SOURCE: VARIABLE EMISSIONS DATA

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

Use this subgroup to describe temporal variations in the emission

rates given in 15b. Factors entered multiply the rates in 15b.

Skip sources here that have constant emissions.

IVARY determines the type of variation, and is source

- specific:

(IVARY) Default: 0

0 = Constant

1 = Diurnal cycle (24 scaling factors: hours 1

- 24)

2 = Monthly cycle (12 scaling factors: months 1

- 12)

3 = Hour & Season (4 groups of 2 4 hourly scaling factors,

where first group is DEC

- JAN

- FEB)

4 = Speed & Stab. (6 groups of 6 scaling factors, where

first group is Stability Class A,

and the speed classes have upper

bounds (m/s) defined in Group 12

5 = Temperature (12 scaling factors, where temperature

classes have upper bounds (C) of:

0, 5, 10, 15, 20, 25, 30, 35, 40,

45, 50, 50+)

--------

a

Data for each species are treated as a separate input subgroup

and therefore must end with an input gro up terminator.

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

INPUT GROUPS: 16a, 16b, 16c

-- Volume source parameters

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

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

Subgroup (16a)

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

Number of volume source s with

parameters provided in 16b,c (NVL1) No default ! NVL1 = 10 !

Units used for volume source

emissions below in 16b (IVLU) Default: 1 ! IVLU = 1 !

1 = g/s

2 = kg/hr

3 = lb/hr

4 = tons/yr

5 = Odour Unit * m**3/s (vol. flux of odour compound)

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 41 6 = Odour Unit * m**3/min

7 = metric tons/yr

Number of source

- species

combinations with variable

em issions scaling factors

provided below in (16c) (NSVL1) Default: 0 ! NSVL1 = 0 !

Number of volume sources with

variable location and emission

parameters (NVL2) No default ! NVL2 = 0 !

(If N VL2 > 0, ALL parameter data for

these sources are read from the VOLEMARB.DAT file(s) )

!END!

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

Subgroup (16b)

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

a

VOLUME SOURCE: CONSTANT DATA

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

b

X UTM Y UTM Effect. Base Initial Initial Emission

Coordinate Coordinate Height Elevation Sigma y Sigma z Rates

(km) (km) (m) (m) (m) (m)

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

*** Max short

- term emission rates ***

Xutm Yutm Hgt Elev

! SRCNAM = VOL01 !

! X = 901.865,

- 1163.096,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL02 !

! X = 906.319,

- 1165.671,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL03 !

! X = 910.773,

- 1168.245,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL04 !

! X = 915.228,

- 1170.820,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL05 !

! X = 919.682,

- 1173.394,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL06 !

! X = 924.136,

- 1175.969,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRC NAM = VOL07 !

! X = 928.591,

- 1178.543,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL08 !

! X = 933.045,

- 1181.118,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL09 !

! X = 937.500,

- 1183.692,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

! SRCNAM = VOL10 !

! X = 941.954,

- 1186.267,19.41,0,2393,13.7,8.13E

- 02,1.42E

- 03,2.66E+00,0,7.35E

- 05,4.39E

- 02,1.33E

- 02,0,0 !

! END !

--------

a

Data for each source are treated as a separate input subgroup

and therefore must end with an input group terminator.

b

An emission rate must be entered for every pollutant modeled.

Enter emission rate of zero for secondary pollutants that are

modeled, but not emitted. Units are specified by IVLU

(e.g. 1 for g/s).

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

Subgroup (16c)

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

a

VOLUME SOURCE: VARIABLE EMISS IONS DATA

---------------------------------------- Use this subgroup to describe temporal variations in the emission

rates given in 16b. Factors entered multiply the rates in 16b.

Skip sources here that have constant emissions. F or more elaborate

variation in source parameters, use VOLEMARB.DAT and NVL2 > 0.

## 42 IVARY determines the type of variation, and is source

- specific:

(IVARY) Default: 0

0 = Constant

1 = Diurnal cycle (24 scaling factors: hours 1

- 24)

2 = Monthly cycle (12 scaling factors: months 1

- 12)

3 = Hour & Season (4 groups of 24 hourly scaling factors,

where first group is DE C

- JAN

- FEB)

4 = Speed & Stab. (6 groups of 6 scaling factors, where

first group is Stability Class A,

and the speed classes have upper

bounds (m/s) defined in Group 12

5 = Temperature (12 scaling factors, where temperature

classes have upper bounds (C) of:

0, 5, 10, 15, 20, 25, 30, 35, 40,

45, 50, 50+)

--------

a

Data for each species are treated as a separate input subgroup

and therefore must end with an input group terminator.

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

INPUT GROUPS: 17a & 17b

-- Non

- gridded (discrete) receptor information

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

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

Subgroup (17a)

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

Number of non

- gridded receptors (NREC) No default ! NREC = 1086 !

!END!

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

Subg roup (17b)

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

a

NON

- GRIDDED (DISCRETE) RECEPTOR DATA

------------------------------------ X UTM Y UTM Ground Height b

Receptor C oordinate Coordinate Elevation Above Ground

No. (km) (km) (m) (m)

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

00001 ! X = 919.692,

- 1083.836, 0.0, 0 ! !END! al_rec_3plus2 5nm

- 0001

00002 ! X = 918.754,

- 1084.200, 0.0, 0 ! !END! al_rec_3plus25nm

- 0002

00003 ! X = 917.826,

- 1084.587, 0.0, 0 ! !END! al_rec_3plus25nm

- 0003

00004 ! X = 916.890,

- 1084.956, 0.0, 0 ! !END! al_rec_3plus25nm

- 0004

00005 ! X = 915.947,

- 1085.306, 0.0, 0 ! !END! al_rec_3plus25nm

- 0005

...(snip

– reduced for brevity) ...

00066 ! X = 856.056,

- 1095.963, 0.0, 0 ! !END! al_rec_3plus25nm

- 0066

00067 ! X = 855.063,

- 1095.805, 0. 0, 0 ! !END! al_rec_3plus25nm

- 0067

00068 ! X = 854.074,

- 1095.624, 0.0, 0 ! !END! al_rec_3plus25nm

- 0068

00069 ! X = 853.090,

- 1095.417, 0.0, 0 ! !END! al_rec_3plus25nm

- 0069

00070 ! X = 852.110,

- 1095.189, 0.0, 0 ! !END! al_rec_3plus25nm

- 0070

00071 ! X = 851.137,

- 1094.936, 0.0, 0 ! !END! al_rec_3plus25nm

- 0071

00072 ! X = 777.710,

- 1118.013, 0.0, 0 ! !END! bretwild

- 0001

00073 ! X = 779.971,

- 1115.939, 0.0, 0 ! !END! bretwild

- 0002

00074 ! X = 780.697,

- 1114.937, 0.0, 0 ! !END! bretwild

- 0003

00075 ! X = 781.422,

- 1113.936, 0.0, 0 ! !END! bretwild

- 0004

00076 ! X = 785.607,

- 1106.067, 0.0, 0 ! !END! bretwild

- 0005

00077 ! X = 789.227,

- 1101.058, 0.0, 0 ! !END! bretwild

- 0006

00078 ! X = 789.783,

- 1098.197, 0.0, 0 ! !END! bretwild

- 0007

00079 ! X = 791.230,

- 1096.193, 0.0, 0 ! !END! bretwild

- 0008

00080 ! X = 791.146,

- 1095.264, 0.0, 0 ! !END! bretwild

- 0009

00081 ! X = 791.785,

- 1093.333, 0.0, 0 ! !END! bretwild

- 0010

00082 ! X = 791.701,

- 1092.404, 0.0, 0 ! !END! bretwild

- 0011

00083 ! X = 792.339,

- 1090.473, 0.0, 0 ! !END! bretwild

- 0012

00084 ! X = 792.256,

- 1089.543, 0.0, 0 ! !END! bretwild

- 0013

00085 ! X = 792.172,

- 1088.614, 0.0, 0 ! !END! bretwild

- 0014

00086 ! X = 792.088,

- 1087.685, 0.0, 0 ! !END! bretwild

- 0015

00087 ! X = 792.00 4,

- 1086.756, 0.0, 0 ! !END! bretwild

- 0016

00088 ! X = 791.920,

- 1085.826, 0.0, 0 ! !END! bretwild

- 0017

00089 ! X = 791.753,

- 1083.968, 0.0, 0 ! !END! bretwild

- 0018

00090 ! X = 792.559,

- 1083.896, 0.0, 0 ! !END! bretwild

- 0019

00091 ! X = 792.475,

- 1082.967, 0.0, 0 ! !END! bretwild

- 0020

00092 ! X = 791.585,

- 1082.110, 0.0, 0 ! !END! bretwild

- 0021

00093 ! X = 792.391,

- 1082.037, 0.0, 0 ! !END! bretwild

- 0022

0 0094 ! X = 791.502,

- 1081.181, 0.0, 0 ! !END! bretwild

- 0023

00095 ! X = 792.307,

- 1081.109, 0.0, 0 ! !END! bretwild

- 0024

00096 ! X = 791.418,

- 1080.252, 0.0, 0 ! !END! bretwild

- 0025

00097 ! X = 791.334,

- 107 9.323, 0.0, 0 ! !END! bretwild

- 0026

00098 ! X = 790.446,

- 1078.467, 0.0, 0 ! !END! bretwild

- 0027

## StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26

- 27189A 43

00099 ! X = 791.250,

- 1078.394, 0.0, 0 ! !END! bretwild

- 0028

00100 ! X = 790.362,

- 1077.538, 0.0, 0 ! !E ND! bretwild

- 0029

00101 ! X = 791.167,

- 1077.465, 0.0, 0 ! !END! bretwild

- 0030

00102 ! X = 790.279,

- 1076.609, 0.0, 0 ! !END! bretwild

- 0031

00103 ! X = 790.195,

- 1075.680, 0.0, 0 ! !END! bretwild

- 0032

00104 ! X = 790.111,

- 1074.751, 0.0, 0 ! !END! bretwild

- 0033

00105 ! X = 789.223,

- 1073.894, 0.0, 0 ! !END! bretwild

- 0034

00106 ! X = 789.140,

- 1072.965, 0.0, 0 ! !END! bretwild

- 0035

00107 ! X = 788.252,

- 1072.109, 0.0, 0 ! !END! bretwild

- 0036

00108 ! X = 788.169,

- 1071.180, 0.0, 0 ! !END! bretwild

- 0037

00109 ! X = 787.281,

- 1070.324, 0.0, 0 ! !END! bretwild

- 0038

00110 ! X = 786.393,

- 1069.467, 0.0, 0 ! !END! bret wild

- 0039

00111 ! X = 785.506,

- 1068.610, 0.0, 0 ! !END! bretwild

- 0040

00112 ! X = 1422.679,

- 1375.246, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0001

00113 ! X = 1422.015,

- 1374.476, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0002

0 0114 ! X = 1421.364,

- 1373.694, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0003

00115 ! X = 1420.726,

- 1372.902, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0004

00116 ! X = 1420.101,

- 1372.099, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0005

...(s nip

– reduced for brevity) ...

00840 ! X = 926.373,

- 1093.508, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0729

00841 ! X = 925.408,

- 1093.793, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0730

00842 ! X = 924.438,

- 1094.058, 0.0, 0 ! ! END! fl_rec_9plus25nm

- 0731

00843 ! X = 923.462,

- 1094.303, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0732

00844 ! X = 922.520,

- 1094.651, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0733

00845 ! X = 921.589,

- 1095.031, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0734

00846 ! X = 920.655,

- 1095.406, 0.0, 0 ! !END! fl_rec_9plus25nm

- 0735

00847 ! X = 850.291,

- 1095.005, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0001

00848 ! X = 849.328,

- 1095.268, 0.0, 0 ! !END! l ams_rec_3plus25nm

- 0002

00849 ! X = 848.432,

- 1095.725, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0003

00850 ! X = 847.527,

- 1096.164, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0004

00851 ! X = 846.614,

- 1096.585, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0005

00852 ! X = 845.690,

- 1096.983, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0006

...(snip

– reduced for brevity) ...

01082 ! X = 728.239,

- 1250.971, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0236

01083 ! X = 727.305,

- 1250.588, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0237

01084 ! X = 726.381,

- 1250.182, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0238

01085 ! X = 725.465,

- 1249.757, 0. 0, 0 ! !END! lams_rec_3plus25nm

- 0239

01086 ! X = 724.559,

- 1249.311, 0.0, 0 ! !END! lams_rec_3plus25nm

- 0240

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

a

Data for each receptor are treated as a separate input subgroup

and therefore must end with an input group terminator.

b

Receptor height above ground is optional. If no value is entered,

the receptor is placed on the ground.