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 International Corporation

Baton Rouge, Louisiana

Date:

July 2012

Project Number:

26-27189A

StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26-27189A i

Contents

Page

1 Introduction 1

2 Modeling Methodology 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

Table 6: Sulfur and Nitrogen Deposition Thresholds

List of Figures

Figure 1: Proposed Modeling Domain – Statoil OCS

List of Appendices

Appendix A: Proposed CALPUFF Control Parameters
StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26-27189A 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 exploration 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 GOM. 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 presented 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 permitting 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 PSD regulated pollutants 

include particulate matter (PM), carbon monoxide (CO), nitrogen dioxide (NO2), sulfur dioxide 

(SO2) and ozone (O3). 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 CO2 in the AQIA because there 

are no existing or proposed ambient standards for CO2 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 Americas

(Transocean), or the Maersk Developer (Maersk) drilling rigs. The Maersk is a semisubmersible drilling rig on two pontoons, yet is dynamically positioned. The Transocean rig is a

dynamically positioned drilling vessel with a single hull. Neither requires 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 25 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 Class 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 compliance at 

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

run will be used to assess 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 modeling1

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

Florida but extends far enough to the west to reach 50 km west 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 output 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 Characterization

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 Statoil 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-27189A 3

• Standby Generator Engine (1)

• Cementing Unit Engines (2)

• Life Boats (4)

• Rescue 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 Engines (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 

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

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

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

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

and a point 25 miles away, in the direction of the Breton Wilderness Area. Ten (10) volume 

sources will be used to represent the combined emissions 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 called 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 the 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 characteristics 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 dispersion 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 the 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 water 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 Receptors

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-27189A 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 flexibility 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. The 

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 scenarios: 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 within 

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 143 lease block, which lies 

at 28.851 o

N and 87.424 oW. For purposes of this modeling analysis, Statoil will model all 

scenarios from this northwest corner. Statoil expects impacts to be the same for all locations 

within the lease blocks. This location has 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 (NOx) 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 below 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/m3)

NO2 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 

ambient concentrations resulting from the proposed project for pollutants that trigger PSD 

review. Because Statoil proposes to emit NOx in excess of 40 TPY, a PM2.5 review is also 

triggered.

The results of the preliminary analysis determine whether a full impact 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 Full Impact Assessment is required. 

The EPA has not identified a regulatory-default model for secondary PM2.5. Therefore, only 

concentrations resulting from proposed primary PM2.5 emissions will be compared to the PM2.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/m3)

NAAQS

(μg/m3)

Class II PSD 

Increment

(μg/m3)

NO2 1-hour 7.5(1) 188(2) N/A

Annual 1 100 25

PM2/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 NO2

NAAQS for the Prevention of Significant Deterioration Program”. 2 To attain this standard, the 5-year average of the 98th percentile of daily maximum 1-hour concentrations must not 

exceed 188 μg/m3

, effective 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/m3.
StatOilGulf Services, LLC

Air Quality Analysis Modeling Protocol

Project No. 26-27189A 7

4 Additional Impacts 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 surrounding 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 analysis 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 default 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 Visual 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 quality 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 in 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 Statoil Lease Blocks

Class I Area

Distance1

(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 corner 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). Statoil 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). Features 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 visibility. 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-27189A 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), outlined 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, Number 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 NOx, SO2, H2SO4, and PM10 for 

each drilling vessel and associated fleet, an initial “Q/D” analysis yields values in excess 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

Distance1

(km)

Q/D

Maersk Transocean

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 extends 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 sources, 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 cells 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 Coordinate (LCC) system, using the 

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

The CALPUFF modeling will consider emissions of NOx, SO2, and PM for the evaluation of 

AQRVs at the Class I areas.5

 However, the CALPUFF modeling will calculate only 

concentrations of NOx and primary PM2.5 for the evaluation of PSD increment consumption at 

the Class I areas.

5.1.5 Chemical Transformations

5.1.5.1 Ammonia

The NOx chemistry in CALPUFF depends on the ambient ammonia concentration to establish 

the equilibrium 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 

ammonia 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 data, 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 associated 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-27189A 11

Statoil will use these plots 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 Meteorological 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 stations, 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 use 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 different 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) not 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 SILs 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 the 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 inventory.

The EPA has not identified a regulatory-default model capable of calculating secondary PM2.5. 

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

emissions of PM2.5 will be used. For all other AQRVs (evaluated by the FLMs), both primary 

and secondary PM2.5 will be included.

Table 5. Class I Area SILs and Increments

Pollutant Averaging Period Significance Level 

(μg/m3)

Class I Increment 

(μg/m3)

NO2 Annual 0.1 2.5

PM2.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 estimating the visibility impacts due to 

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

2010). The FLAG visibility metric is the estimated maximum 24-hour change in extinction (∆bext) 

over clean natural visibility conditions (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-27189A 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 the 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 allowed 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 extinction equation to calculate bext

(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 concentrations recommended by the FLMs for the Class I area 

of concern, and assesses the visibility using the 98th percentile modeled values at each 

receptor. In order to use Method 8, CALPOST Version 6.221 (Level 080724) will 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:

bproject = 2.2 × fS(RH) × [Small Sulfate] + 4.8 × fL(RH) × [Large Sulfate] 

+ 2.4 × fS(RH) × [Small Nitrate] + 5.1 × fL(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 × fSS(RH) × [Sea Salt] 

+ Rayleigh Scattering (Site Specific) 

+ 0.33 × [NO2 (ppb)] {or as: 0.1755 × [NO2 (μg/m3)]} 

Where: 

[ ] indicates concentrations in μg/m3 

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

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

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

For Total Sulfate < 20 μg/m3: 

[Large Sulfate] = ([Total Sulfate] / 20 μg/m3) × [Total Sulfate] 

For Total Sulfate ≥ 20 μg/m3: 

[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 require estimation of total sulfur deposition from the CALPUFFestimated wet and dry SO2 plus sulfate (SO4) deposition. For nitrogen, wet and dry deposition 

from all of the nitrogen species modeled species are included (NO2, nitric acid, and particulate 

nitrate).

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 nitrogen 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 nitrogen 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 to 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 SO2, Statoil will 

not present the sulfur deposition results in the application.

Table 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-27189A 15

6 References

FLAG. 2000. “Federal Land Managers’ Air Quality 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, Research 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 Regulations (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

-800 CC 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, 186x108x4km 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.2003.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 * T2DDAT = *

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 = *

MASSFLX.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/calmet/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/calmet/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/vistas/2003/subdomain2_4km/met2003-dom2-dec.dat ! !END!

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

INPUT GROUP: 1 -- General run control parameters

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

 Option to run all periods found

 in the met. file (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)

 Default: 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 checking

 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 every 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 parameters 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) (PGTIME)

 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 stack 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 modeled 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 ! MSPLIT = 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 computed 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)

 Aqueous 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 removal 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 species)

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

 (valid 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 coefficients 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

 [DIAGNOSTIC 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 NOT 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) Default: 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 TIBL 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 mechanically-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 are 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 must 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 model 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 : Tangential 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) No Default ! IUTMZN = 0 !

 Hemisphere for UTM projection?

 (Used only if PMAP=UTM)

 (UTMHEM) Default: N ! UTMHEM = N !

 N : Northern hemisphere projection

 S : Southern hemisphere projection

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

 Earth 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 parameters is provided by the National Imagery and

 Mapping Agency (NIMA).

 NIMA Datum - Regions(Examples)

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

 WGS-84 WGS-84 Reference Ellipsoid and Geoid, Global coverage (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 projection 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 left (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 <= NY)

 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 RECEPTORS):

 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

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

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

 locations of sources, receptors, grids ...)

 suitable for plotting?

 (IQAPLOT) Default: 1 ! IQAPLOT = 1 !

 0 = no

 1 = yes

 DIAGNOSTIC MASS FLUX OUTPUT OPTIONS:

 Mass flux across specified 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 OUTPUT 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 not 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 print 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

 Messages 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 !

 Note: 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 complex 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 to 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 RECEPTOR INFORMATION

 XRCT YRCT ZRCT XHH

(km) (km) (m)

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

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

1

 Description of Complex Terrain 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 terrain 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.1656, 1., 8., 5., 3.5 !

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

!END!

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

INPUT GROUP: 8 -- Size 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 deposition 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 for 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 deposition 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

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

 

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

 Monthly 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 !

 H2O2 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 concentrations

 (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 characterize 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. 50. 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 Heffter 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 boundary

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

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

 Latitude (degrees) for met location

 (XLATIN) Default: -999. ! XLATIN = -999.0 !

 Longitude (degrees) for met location

 (XLONIN) 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 sigma-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 OBSERVED 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 !

 Maximum 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 step

 that includes gradual rise (for CALMET

 and PROFILE winds)

 (NCOUNT) Default: 2 ! NCOUNT = 2 !

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

 (SYMIN) Default: 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!

 Divergence 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.23, 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 eligible 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

 (IRESPLIT(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 develops)

 (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

 into 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 integration

 (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 ! NSPT1 = 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 !

! SRCNAM = 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 input 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 between 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 with 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 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)
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

 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: 14a, 14b, 14c, 14d -- Area source parameters

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

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

Subgroup (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 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 IARU 

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

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

Subgroup (14c)

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

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

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

Source a

No. Ordered list 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 terminator.
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 in 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 subgroup

 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

 combinations 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 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 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 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 subgroup

 and therefore must end with an input group terminator.

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

INPUT GROUPS: 16a, 16b, 16c -- Volume source parameters

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

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

Subgroup (16a)

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

 Number of volume sources 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

 emissions 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 NVL2 > 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 !

! SRCNAM = 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 EMISSIONS 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. For 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 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 group terminator.

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

INPUT GROUPS: 17a & 17b -- Non-gridded (discrete) receptor information

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

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

Subgroup (17a)

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

 Number of non-gridded receptors (NREC) No default ! NREC = 1086 !

!END!

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

Subgroup (17b)

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

 a

 NON-GRIDDED (DISCRETE) RECEPTOR DATA

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

 X UTM Y UTM Ground Height b

Receptor Coordinate Coordinate Elevation Above Ground

 No. (km) (km) (m) (m)

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

00001 ! X = 919.692, -1083.836, 0.0, 0 ! !END! al_rec_3plus25nm-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.004, -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

00094 ! 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, -1079.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 ! !END! 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! bretwild-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

00114 ! 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

...(snip – 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! lams_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.