EPA-450/3-77-006
March 1977
EMISSION DENSITY
ZONING
U.S. ENVIRONMENTAL PROTECTION AGENCY
Office of Air and Waste Management
Office of Air Quality Planning and Standards
Research Triangle Park, North Carolina 27711
-------
EPA-450/3-77-006
EMISSION DENSITY ZONING
by
Frank H. Benesh, Phillip D. McLellan,
Michael T. Mills, and Robert Patterson
GCA/Technology Division of GCA Corporation
Bedford, Massachusetts
Contract No. 68-02-1376
Task Order No. 26
EPA Project Officer: John L. Robson
Prepared for
ENVIRONMENTAL PROTECTION AGENCY
Office of Air and Waste Management
Office of Air Quality Planning and Standards
Research Triangle Park, North Carolina 27711
March 1977
-------
This report is issued by the Environmental Protection Agency to report
technical data of interest to a limited number of readers. Copies are
available free of charge to Federal employees, current contractors and
grantees, and nonprofit organizations - in limited quantities - from the
Library Services Office (MD-35), Research Triangle Park, North Carolina
27711; or, for a fee, from the National Technical Information Service,
5285 Port Royal Road, Springfield, Virginia 22161.
This report was furnished to the Environmental Protection Agency by
GCA/Technology Division of GCA Corporation, Bedford, Massachusetts,
in fulfillment of Contract No. 68-02-1376, Task Order No. 26. The contents
of this report are reproduced herein as received from GCA/Technology
Division of GCA Corporation. The opinions, findings, and conclusions
expressed are those of the author and not necessarily those of the Environ-
mental Protection Agency. Mention of company or product names is
not to be considered as an endorsement by the Environmental Protection
Agency.
Publication No. EPA-450/3-77-006
11
-------
ABSTRACT
Emission density zoning is an air quality control strategy whereby
all parcels of land within an air quality maintenance area, excluding
public rights of way and reserved open space, are assigned maximum legal
emission allowances for particulate matter and sulfur oxides, expressed
in terms of mass' of pollutant per time period per lot area. The subject
of this publication is the development of criteria for (a) partitioning
an air quality maintenance area into component areas, (b) the setting of
emission density limits for all land uses in each component area and,
separately, all major emission sources, and (c) the revision of the emis-
sion density limits as conditions change. The criteria are then trans-
lated into a generalizable decision model, a linear program, for setting
the emission density limits. The criteria and model are tested in two
air quality maintenance areas, Baltimore and Louisville.
iii
-------
CONTENTS
Page
Abstract ill
List of Figures viii
List of Tables xi
Sections
I Introduction and Summary 1
Introduction 1
Organization of the Report 5
Summary of Project 6
Findings and Limitations 7
II Criteria for Partitioning AQMAs 12
Introduction 12
Air Quality Goals 13
Land Use Goals 20
Administrative Goals 22
III Criteria for Setting EDZ Limits 24
Economic Goals 25
Air Quality Goals 27
Land Use and Environmental Planning Goals 31
-------
CONTENTS (continued)
Sections
Administrative Goals 33
IV Criteria for Revising EDZ Limits 35
Criteria for Revising Limits Without an Emission Rights
Market 36
Criteria for Revising Limits With an Emission Rights
Market 37
V Specification of Decision Model 40
Maximization of Emissions (Model No. 1) 43
Minimization of Control Costs (Model No. 2) 48
Alternative Formulation of Control Cost Minimization
(Model No. 3) 49
Land Value (Model No. 4) 50
Variants of Model No. 1 52
VI Guidelines and Model Implementation Considerations 56
Resource Requirements 57
Linear Programming Considerations 58
Guideline for Executing the EDZ Model 62
VII An Example of the Emissions Maximization Formulation in
the Baltimore AQMA 95
Data Compilation 95
Execution of MPS 100
Optimal Solution to Sulfur Oxides Maximization 103
Particulate Matter Example 115
VIII An Example of the Emissions Maximization Formulation in
the Louisville AQMA 119
vi
-------
CONTENTS (continued)
Sections Page
Data Compilation 119
Execution of MPS 123
Particulate Matter Example 123
Sulfur Oxides Example
IX References 140
Appendixes
A An Introduction to Linear Programming 144
B Provision for Short-Term Standards in Air Quality Constraint 148
C The SO -PM Feasibility Constraint 160
D Specification and Utilization of Control Cost Parameters 169
E The Application of the Emission Density Setting Model in an
Attainment Situation . 178
vii
-------
LIST OF FIGURES
No. Page
1 Offset Gridding 15
2 Receptor Locations Outside an AQMA 19
3 Test City Example Output 72
4 ROWS Section of SOX Emissions Maximization Program 98
5 COLUMNS Section of SOX Emissions Maximization Program 99 .
6 RHS Section of the SOX Emissions Maximization Program 101
7 BOUNDS Section of SOX Emissions Maximization Program 102
8 MPS Control Program 104
9 ROWS Section of SO Optimal Solution 105
10 COLUMNS Section of SO Optimal Solution 106
11 Slack SO Levels (micrograms per cubic meter) 109
X
12 Status of Bounds on Emission Density Limits, Current
Light Industry 110
13 Status of Bounds on Emission Density Limits, Future
Light Industry 111
14 Status of Bounds on Emission Density Limits, Current
Heavy Industry 113
15 Status of Bounds on Emission Density Limits, Future
Heavy Industry 114
16 Slack TSP Levels in Baltimore (micrograms per cubic meter) 116
17 Status of Bounds on Emission Limits, Current Light Industry 117
viii
-------
LIST OF FIGURES (continued)
No. Page
18 Louisville AQMA Component Area Grid System 120
19 Numbering of Louisville AQMA Component Areas 121
20 Slack SOX Levels (micrograms per cubic meter) 126
21 Projected SOX Concentrations (microgram per cubic
meter) 127
22 SOX Levels Attributed to EDZ Sources (microgram
per cubic meter) 128
23 SOX Levels Attributed to non-EDZ Sources
(micrograms per cubic meter) 129
24 Status of Bounds on Emission Limits, Current Light Industry 130
25 Status of Bounds on Emission Limits, Future Light Industry 131
26 Status of Bounds on Emission Limits, Current Heavy Industry 132
27 Status of Bounds on Emission Limits, Future Heavy Industry 133
28 SOX Emission Density, Base Year (tons per day per
square kilometer) 134
i
29 SOX Emission Density, Design Year (tons per day
per square kilometer) 135
30 SOX Emission Density, All Sources at Upper Limit
(tons per day per square kilometer) 136
31 Construction of Feasible Region 146
32 Determination of Optimal Solution 146
33 Location of TSP Monitoring Stations and Associated Geometric
Standard Deviation 154
34 Hypothetical Relationship Between SO and PM Emissions 161
X
35 Plot of Baltimore Heavy Industrial Point Source PM Emission
Rates Versus SOX Emission Rates (tons/year) 164
ix
-------
LIST OF FIGURES (continued)
No. Page
36 Plot of Baltimore Heavy Industrial Point Source PM Rates
Versus SOX Emissions Rates, Sources With Emissions Less
Than 100 Tons/Year 165
37 Plot of Baltimore Light Industrial Point Source PM Emission
Rates Versus SOX Controlled Emission Rates (tons/year) 166
38 Plot of Baltimore Commercial Point Source PM Emission Rates
Versus SOX Controlled Emissions Rates, Excluding Coal
Burning Sources (tons/year) 167
39 Scheduel of Emission Reductions Versus Control Costs for
Prototypical Heavy Industrial Point Source 174
40 Plot of Linear Approximation of the Control Cost - Emission
Reduction Schedule 175
-------
LIST Of TABLES
No^ Page
1 Definitions of Variables Used in Section V 44
2 Control Cost and Emission Schedule, Example of Prototypical
Heavy Industrial Point Source 51
3 Problem Characteristics and MPS Timing Estimates, Emissions
Maximization Linear Programming 60
4 Requisite Data for Preparing an Emissions Maximization
Linear Program 62
5 Test City Base Year Emission Inventory, SO (tons/day) 67
X
6 Test City Land Use Inventory (square kilometers) 68
7 Calculation of the Right-Hand Side 69
8 1974 TSP Concentrations Recorded at Monitoring Stations
in the Baltimore AQMA 150
9 1974 SO Concentrations in the Baltimore AQMA 156
x
10 1980 to 1990 Projections of Growth in Manufacturing
Earnings ($000) 170
11 Emission Reductions and Associated Control Costs 173
xi
-------
SECTION I
INTRODUCTION AND SUMMARY
INTRODUCTION
Emission density zoning is an air quality control regulation whereby
a maximum legal rate of air pollutant emissions is prescribed for certain
parcels of land. The limits are expressed in terms of mass of pollutant
per time period per lot area. The maximum legal rate, or emission density
limit, would vary with the location and the planned use of the land parcel.
The concept of an emission density strategy is not new. An early
discussion of its potential was published by the National Air Pollution
Control Administration in 1968.29 Subsequently, John Roberts and Edward
Croke of Argonne National Laboratory contributed greatly to a definition
of an EDZ control strategy, its,potential, and its limitations.30*31
Emission density related control strategies, similar but not the same
as those contemplated in this report, have been implemented in Chicago,
Illinois and Louisville, Kentucky.26'32
In concept, it is very similar to the typical use of density and
bulk regulations in land use zoning. Each parcel of land is assigned
a permitted use; e.g., residential, commercial, and the permitted uses
have associated maximum densities. For example, a residential zone may
have a maximum density of two dwellings per acre; i.e., a minimum lot
size of 1/2 acre. Similar restrictions are placed on commercial and
office zones, frequently in terms of a floor area ratio (FAR), which
specified the maximum amount of building floor area per lot area. For
-------
example, a four-story building covering one-half of the lot would have
a FAR of two.
Land-use zoning is usually applied to all areas of the city, both
developed and undeveloped. Developed; i.e., currently built on, land
can expand or redevelop up to the maximum permitted density and use for
its zoning class. One land use may have differing maximum permitted
densities in various areas of the city; e.g., for residences, one may
find areas zoned for a minimum lot size of 2 acres and other areas in
the same city zoned for 1/2-acre lots. Emission density zoning could
also be implemented in this manner, such that all areas of the region
are controlled and that similar land uses have different maximum emis-
sion densities, depending on their location within the city.
This report describes the criteria for setting the emission density
limits in a metropolitan region, a model for setting the limits, and
examples of the application of the criteria and model. The region
considered is an air quality maintenance area1 because this is convenient
for air quality planning purposes (though other regions could be em-
ployed) . The model sets emission density limits for a grid network of
squares called component areas. It is assumed that either all the land
use within a particular category in a component area will have identical
emission density limits, or that a local jurisdiction will reassign
emission density limits reflecting local planning goals but not increas-
ing the total permitted emissions within the component area.
In a manner similar to land-use zoning, the emission density zoning
model described in this report is based on a land use plan. It adopts
the amount and configuration of land use for a specified future date,
called a design year, and sets a pattern of emission density limits that
achieve an adopted objective. The emission density limits are suf-
ficiently high on every parcel of land to permit development as contemplated
in the land use plan. However, in some areas they would be so stringent as
-------
to permit only certain types of the broad category of development out-
lined in the land-use plan, to require the application of emission control
techniques, or to require development of a lesser density; i.e., the same
kind of development on a larger parcel of land. In other areas, the po-
tential development will not be so restricted.
This dependence on a region's land-use plan, as well as specific
implementation concerns, requires the substantial involvement of city
and regional planners as well as air pollution control engineers in the
setting of emission density limits and the implementation of an emission
density zoning ordinance.
It is possible (and, in fact, may be desirable) that emission density
zoning be implemented with the provision for transferable emission rights.
This would be similar to transferable development rights in conventional
zoning. A market for transferable emission rights is attractive in that
it will allow the pattern of set emission limits to change with changes in
growth patterns, yet still maintain the National Ambient Air Quality
Standards (NAAQSs). Transferable emission rights may be necessary to ensure
the legal and political acceptability of emission density zoning. The En-
vironmental Protection Agency's proposed Emission Offset policy (41 FR
55524 and 55558) almost certainly will stimulate development of such a mar-
ket. The required new institutional mechanisms, however, are not considered
in this report.
Emission density zoning, as contemplated in this report, would in-
clude the following procedures:
• Each compartment area and, separately, each major stationary
source, will be assigned emission density limits by the
regional authority that are the maximum densities of
annual* emissions of sulfur oxides (SOX) and particulate
matter (PM).
Annual emission density limits, to be effectively enforced, probably
will require translation into daily limits.
-------
• Except for major stationary sources, all land of the same
land use category in a component area, as specified in
land use plans and/or zoning ordinances, will be assigned
the same emission density.
• Every parcel of land in the region, excluding public
rights-of-way and reserved open space, will have a fixed
emission density limit.
• Local jurisdictions may reallocate but not increase emis-
sion limits within a component area as local plans, zoning,
or planning goals change. Local jurisdictions may reallo-
cate emissions among component areas under their authority
subject to the constraints discussed in Section IV of this
report.
• A market for transferable emission rights among industrial,
and perhaps commercial, sources may be instituted according
to the procedures described in Section IV of this report.
• Variances may be granted by a regional authority subject to
the limitations described in Sections III and IV.
Brail26 has defined four types of emission quota strategies, viz.
• Emission Allocation Planning - (jurisdictional emission
quotas) under which the air pollution control agency
allocates permissible emission quotas by local govern-
ment jurisdiction within its region. No attempt is made
directly to apply the emission quota to new stationary
sources.
• Floating Zone Emission Quotas - in which a limitation
on pollutants generated within an area of a specified
size which can be drawn about any specific location
within the metropolitan area.
• District Emission Quotas, a step down from emission
allocation planning, a strategy which limits the amount
of pollutants to be emitted during some time period from
a planning district within a jurisdiction.
• Emission Density Zoning in which the emissions of a
pollutant be limited to prescribed levels for a selected
unit area.
The subject of this report does not fit neatly into any of the above
categories. The output of the decision model does provide unique emission
-------
density limits for every major point source; thus it may be called emis-
sion density zoning. It also provides emission density limits for every
area source. The area source emission density limits are uniform for all
land use of a particular category within a component area; i.e., grid
square. Thus, the emission density limits for a specific land use within
a component area can be summarized as an emission quota for that land
use, an approach similar to district emission quotas. Moreover, in the
contemplated implementation scenario, the total emissions in a municipality
are treated as a jurisdictional emission quota which the municipality
would be free (and encouraged) to reassign. Thus the subject of this
report is a decision model for setting emission density limits, the con-
templated implementation of which may include elements of emission allo-
cation planning and possibly district emission quotas. If desired, the
model can be used as a means of setting jurisdictional emission quotas
or district emission quotas without assigning emission density limits.
ORGANIZATION OF THE REPORT
The remainder of this section is devoted to a summary of the report.
The following three sections of this report describe the criteria for
partioning a region into component areas, setting emission density limits,
and revising the emission density limits at a later date as conditions
change. Next, Section V describes the specification of a mathematical
model, a linear program, for setting the emission density limits. Sec-
tion VI describes a step-by-step procedure for executing the model.
Sections VII and VIII illustrate the application of the criteria and model
in two example regions, Baltimore and Louisville.
There are several appendixes. The first is a brief introduction to
linear programming and should be read first by the reader unfamiliar with
the technique. The second and third discuss two refinements of the model,
the consideration of short-term air quality standards and the consideration
of a functional relationship between sulfur oxide and particulate matter
emission rates. Appendix D discusses the estimation of control costs for
-------
future emission sources, which is necessary for an alternate formulation
of the emission limit setting model. Appendix E discusses the application
of the model if attainment is a problem (i.e., air quality standards are
currently being violated).
SUMMARY OF PROJECT
Criteria for partioning a region into component areas for the pur-
poses of setting emission density limits were developed. This partioning
of a region very closely approximates the process of gridding an air
quality control region (AQCR) for the purposes of dispersion modeling.
More attention, however, is given to minimizing the number of large grids
spanning jurisdictional boundaries.
Criteria for setting emission density limits in component areas were
developed and translated into a generalized mathematical model. The model
was tested in two air quality maintenance areas (AQMAs), Baltimore, Mary-
land and Louisville, Kentucky. Criteria were also developed for revising
the emission density limits due to unforeseen growth or changes in commu-
nity land use plans.
The model for setting emission density limits is in the form of a
linear program. The simplest formulation of the model is one that maxim-
izes the aggregate regional emissions of air pollutants subject to the
constraints that the NAAQSs are not violated and the emission density limit
set for each source be in a feasible range. The maximization of emissions
is assumed to be loosely associated with minimizing the conflict between
maintaining the NAAQSs and regional economic growth. The opportunities
for other formulations of the model are presented; among them is the mini-
mization of control costs of firms complying with emission density limits.
-------
The model for setting emission density limits (EDLS) is based on and
assumes the validity of the land use element of the comprehensive plan for
the region. The emission density zoning (EDZ) model does not alter the
location of development within the broad categories of the plan. The emis-
sion density limits set by the model do place additional constraints on the
nature of development at any one location. Within a broad category of
development; e.g., light industry, the pattern of emission density limits
identifies those locations in the region planned for such development that
can more easily accommodate development with high emission rates and those
areas that are more suited for low emission rates. In other words, the
model identifies the optimal pattern of "carrying capacities" that maximizes
the total emissions in the region (while maintaining the NAAQSs).
FINDINGS AND LIMITATIONS
The emissions maximization formulation of the model was tested in
two AQMAs. In each AQMA it was executed twice, once for particulate matter
and once for sulfur oxides, the two pollutants that are principally from
stationary sources. (Pollutants such as carbon monoxide, which is gener-
ated primarily by mobile sources; e.g., the automobile, are not well-
suited to control by emission density zoning.) Overall, the testing of
the model was successful.
The model, however, assumes attainment of the NAAQSs; i.e., the air
quality in the region is currently below the standards. If this is not
the case, control strategies must first be developed to restrict the
emissions of existing sources so that attainment of the NAAQSs is achieved.
Moreover, the model assumes that if the region develops according to the
land use plan and all future development employs the best available emis-
sion control technology while building facilities at the lowest develop-
ment density economically feasible, the NAAQSs can still be maintained.
The model operates, then, on the middle ground between the extremes of
future development without additional controls (which is assumed to violate
the NAAQSs) and requiring all future development to adopt the most stringent
-------
controls (which is assumed not to violate the NAAQSs). The model can be
revised to accommodate situations outside this range.
Several problems may confront a user of the model. First, the size
and complexity of the linear program increases rapidly with the number
of emission sources; i.e., component areas. For metropolitan areas
roughly the size of'Baltimore or larger, the computer resources required
to solve the linear program can become quite large; i.e., over an hour
of computer time. In such cases it is possible to solve a series of
smaller problems. One can solve the problem for a sample; e.g., 10 per-
cent, of the receptors (locations at which air quality is calculated), and
then use the solution to that problem as a starting point to solve the
linear program with a full complement of receptors (or a dense field of
receptors near those receptors projected to approach the NAAQSs).
Second, the user of the model must be able to specify the maximum
emission density that can be expected in each land use category as well
as the minimum emission density that is economically feasible and still
compatible with development according to the land use plan. This is a
critical data item; one can expect to devote a significant amount of man-
power (several man-months) to compiling it. Many emission density limits
will be set at one of these limits.
Third, the general pattern of emission density limits that can be
expected implies that the land use plan must have addressed urban sprawl.
In the model tests that were conducted, emission sources on the periphery
of the region were set at their upper limits and sources at the center
of the region were set at their lower limits. This pattern means that:
• Capacity expansion of existing sources at the periphery
is favored over those in the center of the region.
Centrally located sources cannot expand their capacity,
and hence increase emissions, without reducing their
emissions from current activities. Suburban sources
are not so limited.
-------
• The permitted emission densities on undeveloped land,
within a specified category, are higher on the periphery
of the region and lower'in the center. This would seem
to encourage development away from the center of the
region unless there are more overwhelming locational
advantages at the central location.
While such a pattern of emission density limits does not affect the amount
of land planned for specific uses, it confers on peripheral locations a
decided advantage. The congruence of this phenomena with other planning
goals must be examined. A more balanced allocation of emission limits may
be desirable, facilitating additional development at central locations at
the expense of a lesser aggregate amount of permitted emissions in the
region. It is this attempt at balance that the control cost and land
value formulations of the decision model address.
Implementation Considerations
There are several considerations that will be faced by an agency
implementing EDZ. One should first note that EDZ cannot be used to
control natural and transported particulate matter, and emissions from
nontraditional sources. Background pollution, auto exhaust, reentrained
road dust, fugitive dust from construction, demolition, and quarrying and
dust from unpaved areas and other urban activities fall into this difficult-
to-control category. The contribution of background and nontraditional
sources to each receptor must be projected for the design year. If the
EDLs are met by every controlled source, it is still possible to violate
the NAAQSs if the noncontrolled EDZ sources actual emissions are greater
than was projected.
Implemention of EDZ will require substantial intergovernmental coopera-
tion. A close working relationship between the staffs of the regional air
pollution control agency and the regional planning agency is necessary —
more so than in other aspects of air quality maintenance planning. More-
over, EDZ must be integrated into the regional planning process and must
influence the transportation plan, water resource management plan, and land
-------
use element of the regional plan. The EDZ limit setting model computes an
"optimal" set of EDLs without considering their effect on mobile source
pollutants, water pollution, housing quality, core area revitalization,
conversion of agricultural land, and other regional goals. This means
that the regional plan must have already dealt with these issues. The EDZ
model will then optimize the potential emission density pattern within the
limits of the regional plan. The 208 and transportation plans must be
based on, and be compatible with the regional plan.
Cooperation among and within regional agencies is not all that is
required. EDZ also requires the involvement and commitment of local
governmental units. Ideally, the municipalities in the region will
internally reallocate the regionally-specified EDLs so that they will
reflect local planning goals, and will include the implicit municipal
emission quota as a constraint in subsequent revisions of the local land
use plan. At a minimum, the municipality must include compliance with
the EDL as an additional criterion when issuing building and occupancy
permits. For example, the EDZ ordinance might require a certificate
signed by a professional engineer stating an estimated average daily
emission rate to accompany all permit applications. The municipality
could then check compliance with the maximum permitted emission density
as it now does for compliance with the maximum permitted development
density. If the municipality is unwilling to cooperate in administering
EDZ, a scenario like the one postulated earlier may be unfeasible. An
unexamined issue is how one might require or force a municipality to
administer EDZ.
If municipality-based administration of EDZ is unfeasible, a regional
or state agency will have to administer it. This is less desirable,
because it is likely that only sources above a certain threshold size
could be checked for compliance. This limits the effectiveness of EDZ.
A possible loss of flexibility in a regionally-administered system, and
a possible lack of integration with local planning, also appear to make
a regional system less desirable.
10
-------
Another issue is how to ensure compliance with the EDL once a
facility is operating. Most of the physical changes in a facility that
would alter its emission rate can be controlled through the building per-
mit process. There are, however, at least two other actions that alter
emission densities and may be more difficult to control. An industrial
facility may desire to increase its hours of operation. It may then
not be in compliance with its EDL. More traditional emission control
for point sources, administered by the regional agency, are thus still
necessary in conjunction with EDZ. Also, a facility may change its land
holdings, thereby changing its emission density. Some mechanism must
be designed either to prevent a source from selling vacant land used only
to satisfy EDZ, or to reduce the EDL on the land that was transferred to
a new owner. A transferable emission rights market may alleviate this
problem.
The legal basis for emission density zoning is unclear. While there
appear to be no substantial issues raised if emission density zoning is
implemented as a maintenance strategy,26 a specific implementation
strategy has not been reviewed. For example, because emission density
zoning is patterned after land use zoning, should emission density
limits for the same land use category in a component area differ among
parcels? In conventional zoning, the regulations must be uniform for
each class and kind of building within a district.
Finally, the political considerations inherent in assigning emis-
sion density limits should be noted. First, as previously mentioned, the
issue of balance between the center city of a region and suburban loca-
tions must be addressed by the implementing agency. Second, the assign-
ment of emission densities to component areas and, hence, the allocation
of emissions to local jurisdictions constrains the zoning prerogative of
the jurisdiction. Though the local jurisdiction may be permitted to
reallocate emission rights among districts within its boundaries, the
total emissions permitted from within a grid cell is essentially fixed.
Third, emission density zoning, by assigning limits according to a regional
land use plan, can become a substantial tool for implementing the plan.
11
-------
SECTION II
CRITERIA FOR PARTITIONING AQMAs
INTRODUCTION
Emission density zoning (EDZ) requires the partitioning of an air
quality maintenance area into smaller units. Such divisions are common
in air quality studies, most notably in connection with air quality model-
ing. When an array of square or rectangular units is used, one refers
to this array as a "grid system" comprised of "grid cells" or "grid squares."
The problem in partitioning any area for EDZ is to choose the sizes
of grid cells, and the spatial arrangement of grid cells and receptors,
to account for various concerns of EDZ. So that these concerns may be
better defined, the overall goals of EDZ are broken down into three
categories with particular relevance to the gridding problem. These
categories are maintenance of air quality standards, adequate considera-1
tion of existing and future land use patterns, and administration of EDZ.
There is, of course, some overlap among these categories, but the dis-
tinctions provide a useful framework for discussion of the gridding problem.
In the following discussion the specific concerns within each of these
categories are presented, along with the specific gridding objectives
suggested by them. Criteria for meeting each objective will then be
presented.
Attention should first be drawn to EPA's Guidelines for Air Quality
Maintenance Planning and Analysis, hereafter referred to as the Guidelines.
12 3
Volumes 7, 8, and 13 are particularly relevant to the gridding concerns
of EDZ, and the user should be familiar with the techniques described in
12
-------
them. Volume 7 presents methods for projecting pollutant emissions by
counties, and Volume 13 describes how these emissions can be allocated
to subcounty areas. Volume 8 presents a semiautomated (computerized)
procedure for constructing an area source emissions inventory grid, based
on population densities. A manual overlay gridding technique is also
presented in Volume 13.
As noted previously, an EDZ analysis assumes the"existence of studies
such as those described in the Guidelines. If the grid system for the
air quality analysis and the concurrent emission inventory are prepared
in a manner compatible with the requirements of EDZ, a substantial
amount of effort will be saved.
AIR QUALITY GOALS
The need for air quality maintenance planning is based on the likeli-
hood that growth and development will lead to a violation of air quality
standards. The key to maintenance planning with EDZ is to identify
existing and potential problem areas and to set emission density limits
in these areas sufficient to maintain air quality. Analyses of emission
inventories and air quality monitoring data and estimates of growth in
emissions are sources of this information. It is also necessary to relate
emissions to ambient concentrations through the use of a diffusion model
4
such as the Air Quality Display Model (AQDM). Thus any grid system must
meet the requirements of the diffusion model used to determine this rela-
£
tionship. For example, AQDM requires that all area source grid cells
be squares, and that all receptor locations have positive coordinate values.
*
AQDM is the only model currently capable of producing the source-receptor
file (described in Section III)~used as an input to the decision model.
However, the Climatological Dispersion Model (CDM)5 is being modified to
generate such a file, and may be recommended by EPA for future applica-
tions. CDM currently imposes certain requirements on the grid system,
and these requirements are likely to remain in force for any revised ver-
sions. These requirements are as follows: large grid cells must be in-
tegral multiples of the smallest grid cell in the array; the grid system
outline must be rectangular; and, all grid cells must be squares. These
13
-------
Grid cell side lengths should be between 1 and 8 km, as this size
range ?llows for resolution of most data normally used for air
quality planning. Furthermore, grid cells should be of sizes that can
be easily aggregated or subdivided. The easiest way to obtain a grid sys-
tem that meets these requirements is to use side lengths of 1, 2, 4, or
8 km for all grid cells. The best method for preparing such a grid sys-
tem is to overlay an array of 8 km x 8 km cells on the study area, and
then to subdivide these cells according to the criteria presented below.
It is recommended that the Universal Transverse Mercator (UTM)
coordinate system be used for grid registration since this system is
widely used for other air quality planning purposes.
These considerations can be summarized in the first criterion for
partitioning an air quality maintenance area:
1. Grid cells shall be squares with side lengths between
1 and 8 km, based in the UTM coordinate system. Addi-
tionally, they shall be of sizes that are easily aggre-
gated and subdivided. In order to comply with this
criterion, side lengths of 1, 2, 4, and 8 km are strongly
recommended.
The grid system is the means for spatially allocating existing and
projected area source emissions for modeling. Thus, grid sizes must be
chosen that give adequate resolution of emissions. "Adequate resolution"
simply means that the grid size must be small enough to prevent under-
estimating concentrations, but not smaller than the scale of the data and
requirements also mean that any offset between grid cell boundaries (see
Figure 1) must be an integral multiple of the side length of the small-
est grid cell. In fact, the user should avoid offsets whenever possible,
as they can lead to the use of small, inefficient dummy cells to "fill
out!1 the rectangular boundary. If it is necessary for the user to fill
out an existing grid system with zero-emission dummy cells to make it
rectangular, the added grid cells should be as large as possible; e.g.,
even larger than the 8 km x 8 km limit imposed in the criteria. CDM's
requirements strongly suggest use of the 1-2-4-8 gridding system
described in the text.
14
-------
Figure 1. Offset gridding
15
-------
modeling efficiencies will allow. Thus the higher the area source emis-
sion density, the smaller the grid size required to provide resolution.
Emission densities correspond to the intensity and type of develop-
ment in an area. Industrial land uses typically emit more pollutant
per hectare than residential or commercial uses, the high-rise residential
development typically has a higher emission density than single-family
houses on large lots.
In preparing the grid the user need not know actual emission den-
sities. Instead, he may use other values as proxies, such as the values
used to allocate emissions in maintenance planning analyses.8 The most
common such value is population. Emission densities for several land
uses are considered to be directly related to population densities.
One obvious exception to this relationship is industry. Industrial
emissions are more typically related to employment, process volumes, or
value added. However, since industrial emission densities are among the
highest of any land use type, gridding can be simplified by applying the
smaller grid sizes to heavy industrial areas. The user should keep in
mind that a large part of industrial emissions may be point source emis-
sions, which need not be reflected in the grid.
The following criteria are based on these considerations and allow
the user to prepare a grid system that reflects likely emission densities.
2. Areas having an existing or projected population density
greater than 2,000 persons per square kilometer should be
contained in grid cells no larger than 1 km x 1 km.
3. Areas having an existing or projected population density
between 1,000 and 2,000 persons per square kilometer should
be contained in grid cells no larger than 2 km x 2 km.
4. Areas having an existing or projected population density
between 200 and 1,000 persons per square kilometer should
be contained in grid cells no larger than 4 km x 4 km.
16
-------
5. Areas containing heavy industrial sources* that are not in-
dividually itemized in the emission density limit setting process
should be contained in grid cells no larger than 2 km x 2 km.
The relationship between emission density and grid size also holds true
for concentrations. High-concentration areas (i.e., those areas where
concentrations are near the air quality standard) may be identified from
isopleth maps based on monitored or projected concentrations at selected
receptor points. Resolution of emissions is especially important in
these areas, as great care must be taken in locating any new sources
there. Therefore:
6. Areas in which existing or projected pollutant concen-
trations are 90 percent or more of the applicable air
quality standard shall be contained in grid cells no
larger than 2 km x 2 km.
Problems can still arise in areas having lower concentrations. It
is possible that a given grid cell may have an existing or predicted con-
centration low enough to allow for a substantial amount of development.
However, the spatial distribution of emissions within the cell may lead
to a "hot spot" problem if all the development is clustered in a small
area. Thus, if the cell is larger than 1 km x 1 km, it should be sub-
divided to reflect emission patterns more accurately. This requires that
emissions be allocated to the grid cells. Then:
7. If the projected emission density of a grid cell in the
system produced by criteria 1 through 6 is higher than
the median value for the air quality maintenance area,
and if by subdividing the cell into four equal parts
and reapportioning emissions to the resulting
cells, 40 percent or more of the total emissions for the
original cell can be assigned to a smaller cell, the sub-
division shall be incorporated into the grid system. This
criterion clearly cannot apply to 1 km x 1 km grid cells.
For example, defined as Standard Industrial Classification (SIC)
categories 26 through 33.9
17
-------
The "hot spot" situation just described highlights the importance
of another concern, the location of receptor points. Clearly a strategy
like EDZ requires detailed information on emissions and concentrations
throughout the air quality maintenance area. Therefore:
8. A model receptor point shall be located within each grid
cell, preferably at the center. If an existing monitoring
site is located within a grid cell, that site can be sub-
stituted for the receptor point to lower modeling costs.
One other consideration for receptor location must be mentioned.
Because one of the goals of EDZ is maximizing the use of the assimilative
capacity of the air, the decision model may tend to set higher emission
density limits in grid cells near the outer edge of the air quality
maintenance area. This happens because emission sources located in those
areas have a smaller cumulative impact on air quality within the air
quality maintenance area than those located in the center of the grid
system. This is an undesirable outcome, as it may cause violations of the
NAAQS if sufficient industrial development is planned for along the bound-
ary of the AQMA, and regardless of the type, of development, it may cause
a violation of the allowable PSD increment in the county adjacent to the
AQMA. To be able to control this possibility, the receptor grid should
be extended into the adjacent areas, as shown in Figure 2. Therefore:
9. In AQMA's where there is or may be substantial industrial
development near the AQMA boundary, or where the adjacent
area is a Class II or Class I PSD area, a row of receptor
points shall be added at regular intervals around the out-
line of the grid system at a distance of not less than
4 km. It is recommended that these receptors be placed
at what would be the center points of 8 km x 8 km grid
cells, as illustrated in Figure 2.
Another air quality goal related to model receptor location is pro-
tection of sensitive receptors. The issues involved in identifying these
receptors and setting air quality standards for them are presented in
Prevention of Significant Deterioration.
18
-------
8 km —*|
i r
~r
I
\
X
X
-AQMA GRID
BOUNbARY
•"
X
^
RECEPTORS
Figure 2. Receptor locations outside an AQMA
19
-------
Section III. If the user decides to protect sensitive receptors, receptor
points representing them must be included in model runs. However,
increasing the number of receptor points increases the cost and complexity
of modeling. Therefore:
10. A receptor point shall be placed at the location of each
sensitive receptor designated as "protected" by the
implementing authority. Receptor points located under
previous criteria may be substituted if they are located
within 1 km of the sensitive receptor.
A final air quality goal with relevance to both gridding and re-
ceptor location is prevention of significant deterioration. As currently
proposed, significant deterioration regulations will effectively require
air quality standards more stringent than the NAAQSs to be met in desig-
nated areas. The standard will depend largely on the area's classifica-
tion. Therefore:
11. A single grid cell shall not contain areas with different
classification for significant deterioration requirements,
unless at least one receptor point is located within each
area.
It is quite likely that by following the above criteria, so many
receptors are identified that in a large AQMA the computer resources
necessary to execute the linear program for setting emission density
limits may become prohibitive. Techniques for modifying the receptor
grid to overcome this problem are discussed in Section VI.
LAND USE GOALS
From the discussion of air quality goals it is clear that land use
patterns must be considered in the gridding process. Land use patterns
serve as a proxy for emissions, and projections of future emissions must
be based on future land uses. Furthermore, land use plans represent
community preferences about development. For EDZ to be meaningful it must
be implemented in a way that reinforces and is consistent with the land
use plan; the grid system should reflect the patterns in that plan.
20
-------
For the purposes of setting and administering emission limitations,
the ideal grid would have a single land use in each cell. Otherwise, a
single-number emissions limit must be translated into emissions limits
for a number of source types, unless the community has no particular
land use preferences in a given grid cell. Given the scale of most urban
development and the size range of grid cells, isolating land uses with
the grid is not practical. Ranking of land uses by their air quality
impacts might be possible, but in most cases emissions result from a
combination of activities. Considering these combinations, rather than
individual land uses, in the gridding process will give better informa-
tion about emissions, at least in urbanized areas.
However, it may be desirable in some cases to consider plans more
explicitly, and to grid so as to separate incompatible land uses. An
example might be a planned industrial area bordered by low density housing.
While undue expansion of the industrial area can be effectively prevented
with a low emissions ceiling for a large grid cell containing both uses, .
location of facilities within the cell is not directly addressed. If
the cell were subdivided so that the two areas were placed in separate
grids, limits could be set that would directly reinforce the desired
development patterns.
"Incompatible land uses" must be identified from their differential
air quality impacts, arid not simply from the kinds of activities they re-
present if such a strategy is to be successful. Further, such a decision
is most appropriately made by local jurisdictions. One obvious category
of land uses to be included is sensitive environmental areas. If such
areas are not otherwise protected, they can be separated from other land
uses by grid boundaries, and appropriate emission limits can be set to
prevent development of them.
Another feature of some plans is their consideration of natural
boundaries to development, such as water bodies. Such boundaries can
simplify gridding based on land use patterns. Grid cell subdivision is
not necessary if an effective boundary already exists.
21
-------
12. When the user judges it necessary he may separate in-
compatible land uses with grid boundaries. The deci-
sion should be based on applicability of other criteria,
likelihood of conflicts over land use that may be re-
solved in a way detrimental to air quality goals, and
the presence of natural boundaries that effectively
limit the spread of a given land use.
ADMINISTRATIVE GOALS
EDZ can be implemented and enforced through any of a number of agency
structures, and will require substantial interagency cooperation. Careful
consideration of grid placement can help reduce the potential for such
conflicts. In the California situation cited by Brail, et al.,26 local
jurisdictions were concerned with whether emissions limits assigned by
the State Air Resources Board would be set for local jurisdictions or for
grid cells. In such a situation it would clearly be desirable for each
grid cell to be wholly contained within a single jurisdiction. Again,
given the size range of grid cells this is neither efficient nor feasible.
Rather, it may be desirable to identify areas of high conflict potential
and to subdivide grid cells in these areas to fit jurisdictional boundaries
as closely as possible.
For example, such a situation was encountered in the Louisville test
of these criteria. The Indiana Port Commission has planned construction
of an industrial area and port facility on the Clark County, Indiana shore
of the Ohio River, located within the Louisville Interstate air quality
maintenance area. The Kentucky shore and an island in the River opposite
the proposed site are used for recreational purposes, and Kentucky is
opposed to development of the facility.10 It would clearly be desirable
to have grid boundaries conform to state boundaries in this case, but the
grid scale makes it unlikely. Emission density limits can be set for any
grid cell spanning the state boundaries, based on both states' regional
land use plans. Once these limits have been set, allocation of emissions
to each portion of the grid cell may be done on the same basis. Disputes
22
-------
between the states might then be settled by some negotiating procedure.
The point is simply to facilitate such negotiations by subdividing grid
cells in such areas to the smallest allowable size.
13. If a single grid cell spans jurisdictional boundaries
and the potential for conflicts between jurisdictions
over land uses in the grid cell is high, the cell shall
be subdivided until each grid cell is contained as
nearly as possible within a single jurisdiction.
A second general administrative goal is to simplify the process of
regridding in response to developments in the air quality maintenance
area. First, the need to regrid will be minimized by meeting the grid-
ding criteria in the first gridding effort. Second, the criteria must
also be observed for any regridding effort, and hence provide guidance
as to when regridding is necessary; e.g., for obtaining adequate resolu-
tion of emissions. Finally, selection of grid sizes that can be aggre-
gated or subdivided easily will make regridding simpler.
14. Grid boundaries shall be redrawn whenever changing condi-
tions cause any of the criteria to be violated.
In any administrative situation minimizing costs is important. Since
the criteria are designed to simplify the process of implementing EDZ,
meeting them should help minimize costs. Since cost (especially modeling
costs) is partly dependent on the number of grid cells, unnecessary
subdivisions should be avoided. For example, the gridding technique in
Volume 13 of the Guidelines requires that the user subdivide grid cells
until each portion of the air quality maintenance area for which alloca-
tion data are available is contained as nearly as possible within a single
grid cell. This may be unnecessary for achieving adequate resolution
of emissions for EDZ. Therefore:
15. A cluster of adjoining grid cells may be aggregated if
the resulting cells do not violate any of the previous
criteria.
23
-------
SECTION III
CRITERIA FOR SETTING EDZ LIMITS
The relatively long-term nature of air quality maintenance requires
that those measures used to ensure maintenance be as compatible as
possible with overall community goals. Air quality is oftentimes not as
tangible'an objective as other public desires and is thus open to per-
haps closer scrutiny than more tangible desires such as more jobs. It
is for these reasons that considerations other than air quality must
also be included in setting EDZ limits.
The following discussion has two objectives. First, the criteria
that planners are likely to require are identified. Secondly, these
criteria are expressed in terms compatible with the decision model used
to set emission density limits. Four groups of goals from which cri-
teria can be developed have been identified and are discussed in turn.
These are economic goals, air quality goals, land use and environ-
mental planning goals, and administrative goals.
The discussion assumes the reader's familiarity with the basic ter-
minology of linear programming, as this is the technique used to set emis-
*
sion density limits. Some criteria have thus been formulated as objec-
tives or constraints.
*
The reader unfamiliar with linear programming is referred to Appendix A,
which is a brief introduction to the subject.
24
-------
ECONOMIC GOALS
The economic goals of any region are, in the general sense, well
known and often stated. In a word, prosperity is sought. Economic
prosperity is typically identified with high levels of employment, stable
prices and tax levels, and a diversity of economic activity. EDZ is likely
to influence the location, density, amount, and type of industrial activity
in the region. Therefore, criteria in this area must be established.
Emission density limits can affect urban development either by pre-
cluding certain types of activity altogether or by necessitating controls
for new establishments. Therefore, accounting for these impacts will
require two kinds of information. First, the user must have some idea
of the amount of growth expected to occur during the period covered by
the maintenance plan, and how the growth will be distributed across the
air quality maintenance area. The best sources of this information are
the land use element of the comprehensive regional plan(s) applicable to
the air quality maintenance area, and local plans and zoning. These plans
and regulations reflect community preferences and expectations about
growth, and the emission density limits should reflect the amount and
location of growth called for in them. To accomplish this, the user
must define a value for the projected maximum emissions (without EDZ)
from a given site for the last year of the planning period. These val-
ues can then be used in the decision model in various ways, depending
on the formulation the user selects.
In most instances the user cannot determine exactly which industry or
other development will locate where in a region. For example, most plan-
ning and zoning schemes include only one or two industrial land use cate-
gories; e.g., heavy and light, or general and permitted. Therefore, the
user should select those categories within each land use category that
are expected to show the most growth as prototypes for calculations of
future emissions.
25
-------
With this information, emission density limits can then be set by
maximizing the aggregate emission loading in the region, given the pro-
jected land-use growth. Emission density limits would thus be set, while
keeping within the bounds of minimum and maximum feasible rates that will
allow the development specified in the land-use plan, in a manner that
allows the greatest amount of emissions consistent with the NAAQs. One
assumes that facilitating the maximum emission loading in a region is
loosely synonomous with facilitating the maximum amount of economic growth
consistent with the NAAQSs. Note, however, that in some areas whose eco-
nomic base consists of little high emission density industry, it may not
be true that maximizing potential emission loading will help economic
growth.
The second kind of information that may be required is control costs.
These costs vary by industry, process, and control efficiency, and should
also be estimated for prototype industries. They are most useful when
expressed in their incremental form; i.e., the cost per unit of emission
reduction. If resources permit this information to be obtained, emis-
sion density limits can be set that will reflect the cost of controls for
the projected level of activity. Such costs would be expected to vary
with the location of activity within an air quality maintenance area,
simply because the degree of control needed for a given source would vary
from location to location. Thus emission limits can be established that
represent the least control cost development pattern from an air quality
maintenance perspective, given a desired or projected level of economic
and demographic activity.
Economic goals can thus be summarized in the following criterion.
1. Emission density limits shall reflect the costs of emission
controls associated with the desired level and location of
economic activity for a community, or shall reflect the maxi-
mum amount of emissions in the region consistent with the
NAAQSs.
26
-------
AIR QUALITY GOALS
For the purposes of this report EDZ has as its primary goal the
maintenance of NAAQSs for particulate matter and sulfur oxides. Emis-
sion density limits must be set low enough so that predicted concentra-
tions of these pollutants do not exceed the standards. This means that
the user must know the relative air quality impact of each existing and
potential source at every receptor point in the air quality maintenance
area. This impact is determined by the location of the source relative
to the receptor and by the variables that influence the atmospheric dis-
persion of pollutants.
However, emissions from inventoried sources in an air quality main-
tenance area are not the only components of concentrations at a given
receptor point. Background concentrations resulting from sources outside
the air quality maintenance area are one component; fugitive emissions that
are not inventoried within the air quality maintenance area are another.
These values must be accounted for before emission density limits can be
set.
Existing background and fugitive components are often assumed to be
accounted for by "statistical" calibration of the model. Calibration is
accomplished using simple linear regression, so that predicted concentra-
tions are adjusted by the formula
/>,
X = A + BX (1)
where X is the predicted concentration.
The statistical validity of the calibration process is discussed in two
monographs by Brier.16,17
27
-------
• * ('-I. t
12
Calibration is discussed elsewhere in more detail, and the user
should refer to these discussions if he is unfamiliar with the technique.
The intercept A should include the background and fugitive concentration
components, and can simply be subtracted from the standard being used to
set emission limits for the pollutant in question.
Future background and fugitive components may be larger or smaller than
existing values. The user must attempt to estimate the changes in these
values and to adjust the concentration standards accordingly. Projected
increases would have to be subtracted from the standards, but decreases
might be held in reserve by the implementing agency as a safety factor.*
In any case a safety factor or reserve should be available for several
reasons. Unusual meteorological conditions or unanticipated increases in
background or fugitive emissions could lead to violations of the standards.
Error in the modeling process could lead to underprediction of concentra-
tions. Finally, changes in the land use plan are likely to occur over
time. In order to avoid rerunning the decision model for small changes,
the implementing authority should have a reserve concentration increment
available so that variances can be granted. This point is discussed in
more detail in Section IV.
In summary, the primary air quality constraint can be expressed as:
.J.k,* £ (Ej,k,£,m aj,k,m> lXi>£ - b^ - fi)£ - pi>a - r^
- s for all i,A
At this time there is no approved methodology for inventorying fugitive
emission let alone projecting them. As methodologies are developed, they
should be used. Background concentrations are best estimated by regional
air quality planners familiar with the area and a specific methodology is
not given here. To the extent that these components of the emission inven
tory are omitted, the calibration intercept may account for them. If sub-
stantial sources are uninventoried, it may be prudent to project the cali-
bration intercept to the design year by the population forecast to account
for growth in these sources.
The variable definitions shown in Table 1 are employed.
28
-------
This states that the emissions at a source and the transfer coefficient
between that source and a specified receptor, summed over all sources,
shall not exceed a certain critical value. The critical value is cal-
culated by subtracting the concentration at the specified receptor due
to emission sources not in the EDZ limit setting model from a desired
air quality standard. There would be a similar constraint equation for
each receptor.
The transfer coefficient, t. . , „ is derived from executing AQDM.
!> J >K» *•
This relates annual emissions to annual mean concentrations (i.e., long-
term standards). It is sometimes referred to as x/Q or the coupling coef-
ficient. A method for considering short-term standards is discussed in
Appendix B.
There may be state or local air quality standards in an air quality
maintenance area that are more stringent than the NAAQSs. These can simply
be entered in place of the NAAQSs in the formulation if they are expressed
in the same way for the same pollutants. State and local ambient air .
quality standards may be expressed for different averaging times and with
different definitions of violations; e.g., a 4-hour particulate concentra-
tion standard not to be exceeded more than once a month. State and local
agencies should have sufficient experience with such standards to be able
to make them compatible with the model. Another air quality goal to which
EDZ must be applied is prevention of significant deterioration. This ef-
fectively means that more stringent air quality standards must be met in
13
designated areas. The standards as currently proposed are expressed in
terms of an allowable concentration increment for each designated area of
a certain class. The implementing agency might want to phase the use of
this increment over time to prevent its being rapidly used up in a way that
limits development opportunities. Thus, the standard on which emission
limits are based would represent that portion of the allowable increment
to be used during the period in question.
29
-------
Yet another potential air quality goal that may or may not be codified
is protection of sensitive receptors, namely persons, plants, or materials
whose susceptibility to air pollution warrant their protection with more
stringent air quality standards. Previously, a criterion was presented
requiring model receptor points to be located at any sensitive receptor
that the user would want to protect. Concentration limits sufficient
to protect these receptors could then be assigned to these points in the
process of setting emission limits. The process would involve several
steps. First, the user would have to define which receptors would be
protected, as the number of possibilities could be large. Then
concentration limits sufficient to protect the designated receptors
would have to be established. These could be based on information about
pollutant effects, or, if such information is difficult or costly to ob-
tain, on some arbitrarily defined "safe" value; e.g., the significant
deterioration increment for a Class I area. Whatever the method used, the
standard selected would simply be used in place of the NAAQS for the des-
ignated sensitive receptors. In deciding which sensitive receptors to
include, the user should keep in mind that increasing the number of re-
ceptor points increases the cost and complexity of modeling, and that any
tightening of air quality standards will impose additional costs on the
community.
2. Emission density limits shall be sufficient to maintain the
NAAQSs and any standards based on significant deterioration
requirements and community air quality objectives. The stan-
dard shall be adjusted to account for background concentra-
tions, fugitive emissions, and a safety margin selected by
the implementing authority.
A final technical issue involved in setting limits is the question
of stack heights. Clearly, the user will not know exact stack heights
for future sources and must thus make certain assumptions. The usual
practice is to assume a standard stack height for all area sources of a
given type (i.e., industrial and nonindustrial types), which applies to
existing and future sources. For future point sources the following
criterion should apply:
30
-------
A conservative stack height for all future point sources
shall be assumed for modeling purposes, based on state
and local minimum-stack-height regulations and on exist-
ing sources of the same category in the base year emis-
sion inventory.
LAND USE AND ENVIRONMENTAL PLANNING GOALS
The comprehensive planning process normally involves consideration
of several sets of goals, which, when taken together, define the desired
development pattern for a community. These goals include those on eco-
nomic development and housing, public service, environmental quality,
and transportation.
Because EDZ is a factor in land-use decisions, and especially in
industrial location decisions, it will affect travel patterns in an air
quality maintenance area. It is -necessary to account for this effect
in setting emission density limits. Land use plans normally account in
some way for travel patterns, especially those resulting from work trips.
Future land uses shown in plans are determined in part by access to
transportation. There may be a transportation control plan that imposes
a limit on travel. Emission density limits can be set to account for
these concerns if they are reflected in the land use plan.
EDZ affects land use patterns because different land uses have dif-
ferent emission characteristics. Emission density limits must thus be
set for each land use type in each grid cell. Furthermore, the land use
types for which limits will be set must be those defined in planning and
zoning documents, as these are the categories for which future locational
information is available. This allows the user to tailor the emission
limits to conform to the land use pattern presented in the plan.
Tailoring of emission limits can be accomplished in two steps. First,
in each formulation of the limit setting model, the emission limit is
weighted by the land area available to each source category in a grid cell,
31
-------
based on the regional land use plan. Secondly, the emission limits must
reflect the emission characteristics of the given land uses. For example,
in the emissions maximization model (Model No. 1 of Section V), if emis-
sion limits were based on planned areas alone, a residential area might
receive a higher limit than an industrial area in the same grid cell.
This situation can be prevented by setting minimum emissions for each
land use type in the grid cell based on the planned or projected activity
levels for that land use. These minimum emissions would correspond to the
minimum feasible emissions constraint used in Model No. 1.
Furthermore, certain land uses might be prevented from locating in
an area by setting maximum limits low enough to make control costs pro-
hibitive. If, for example, a residential area is to be protected against
industry, a limit might be set low enough to preclude any industrial loca-
tion there. If no development is to be allowed in a given area in order
to protect certain resources, emissions from that area would not be al-
lowed to increase at all.
EDZ can also be used to aid in the attainment of environmental goals
expressed in the land use plan. For example, water quality management
might be aided, as development types and densities relate to waste loads
entering streams. Water quality planning information might be used to
place a limit on the amount of development occurring in various locations
in a basin over a given period, depending on the phasing of construction
of treatment facilities. EDZ limits could then be set to reflect the cor-
responding level of development.
These concerns can be summarized in the following criteria:
4. Emission density limits shall be set for individual land
use types within a grid cell. The land use types selected
shall be those for which future location data are available;
i.e., existing planning and zoning categories. Since the
numbers and definition of these categories are likely to
vary from municipality to municipality, for uniformity and
modeling efficiency they should be placed within broader
groups. A suggested grouping follows:
32
-------
rural/open space
low density residential
high density residential
commercial/office
institutional/service
light industry
heavy industry.
5. Emission density limits shall reflect development patterns
congruent with the land use and environmental planning
goals of the community, which are expressed as land use
or development plans or as policy statements. This may
be accomplished by setting limits that bar certain sources;
i.e., land use types, from a given area, or encourage their
location in other area.
ADMINISTRATIVE GOALS
Some administrative concerns of EDZ have been mentioned previously,
such as reserving a concentration safety margin to avoid model reruns.
A major administrative goal that was mentioned previously in relation to
gridding is conflict minimization. The implementing agency will have to
balance some competing interests in the process of setting limits. By
meeting the other criteria presented in this section, the agency can
hopefully enter the goals of these interests into the decision model.
Some issues will remain to be resolved. One goal that should be ad-
dressed directly is minimum interference with growth patterns. This goal
is partly met by using emission limits to reinforce the land use plan.
However, it is possible that the decision model could set very different
emission density limits for similar land uses in adjoining grid cells.
The user should guard against this by requiring that a cluster of parcels
of a given land use type has the same emission density limit over its
whole area.
This might be expressed as a model constraint in the following way:
33
-------
Ej,k,*,m - V-.k,i,. - ° I for certain A and J; (3)*
or as a criterion.
6. When continguous areas of the same land use type are
contained within separate grid cells, emission density
limits for the two areas shall be equal.
The user must manually identify each j and A. If any grid cell has no
area in a particular land use class, it cannot be identified in this con-
straint. That is, where a. , = 0 or a... . = 0, the constraint is
meaningless. j'k'm J+A*k»m
34
-------
SECTION IV
CRITERIA FOR REVISING EDZ LIMITS
Once emission density zoning limits have been set (i.e., emissions
quotas have been allocated to grid cells) and implemented, several
conditions may arise that require the revision of the prescribed emission
limits. These conditions can be categorized as
• those arising from unexpected community growth
• those arising from significant changes in community
goals (e.g., transportation and land use plan revisions)
• those arising from unexpected changes in the techno-
logical, administrative, or economic precepts of the
emission density limit setting decision model; e.g.,
changes in fuel prices, changes in state or federal
regulations, or unplanned closings or openings of major
emission sources.
The differences in the scope and influence of these conditions suggest
the application of different methods of revising emission density limits.
In general, small changes can be accommodated through variance proceedings;
only in the event of major changes in these conditions should the limits
be revised globally.
In addition, whether a market is instituted for the exchange of emis-
sion rights will significantly limit the frequency and extent of emission
limit revisions. A market for emission rights will depend on a stable in-
ventory of emission limits to be orderly and fair.
Criteria for revising limits with and without an emission rights
market are outlined separately below.
35
-------
CRITERIA FOR REVISING LIMITS WITHOUT AN EMISSION RIGHTS MARKET
In defining the original set of emission density limits, one will
use the best information available on existing emission sources, commu-
nity growth plans, control costs, etc. However, it is almost certain
that this information is imperfect. One should expect unforeseen and
undesired hardships imposed on individual firms. For an emission den-
sity zoning control strategy to be equitable, and perhaps legal, one
must be able to ameliorate these hardships by revising the emission
density limits in individual cases; i.e., by granting of variances.
To be able to grant a variance without violating the air quality stan-
dards, a reserve is required. Therefore,
1. The emission density limits shall be set with a reserve
concentration at each receptor.
2. To provide for orderly community growth and unanticipated
individual hardship situations, a quasi-judicial air pol-
lution control board or similar institution shall be em-
powered to grant variances to emission density limits
from time to time. Variances could be granted either
to individuals for specific sites or to municipalities
for allocation among sites in component areas. An in-
dividual would need the endorsement of the appropriate
municipality. In each calendar year, up to 10 percent
of the remaining reserve at each receptor could be al-
located in such a manner.
3. Changes in conditions that lower emission levels will be
treated on an individual basis, with the reduction in
emissions in a particular component area accruing to the
air quality reserve at each receptor.
A revision to the transportation and land use plans may alter the
configuration of the emission sources. If the revised plans are based on
new demographic and economic projections, the total quality of projected
emissions may change as well. If these changes are significant, it will
be necessary to revise the emission density limits globally.
36
-------
4. After a significant revision in the transportation and land
use plans, the emission density limits will be redefined ac-
cording to the criteria for setting emission density limits.
5. Existing emission densities at the time of a plan revision
will be used as a lower bound on the redefined emission den-
sity limits.
Aside from individual variances, unplanned closings of major sources,
and plan revisions, other unexpected changes in the technological, ad-
ministrative, or economic assumptions of the emission density limit set-
ting process may occur. For example,
• a New Source Performance Standard that was promulgated
or projected to be promulgated may be rescinded, or
• (in the event the control cost minimization objective
function is employed) changes may occur in the costs of
various control devices or in the relative prices of fuels.
Such unexpected changes may alter the optimal mix of emission density
limits. If the changes are of a large magnitude, it would be desirable
to revise the limits globally. If the changes are small or affect only a
small number of sources, it would be less disruptive to grant variances
to the appropriate sources. While the threshold between the two possible
reactions would have to reflect the experience of the control agency,
6. Sensitivity analysis should be used when a change in the
assumptions that were used in setting the limits occurs.
If the sensitivity analysis indicates a significant, say
more than 20 percent, change in the optimal solution for
more than 10 percent of the sources, then the limits should ,
be reset according to criteria 4 and 5. Otherwise, the
sources could be treated on an individual bases according
to criterion 2.
CRITERIA FOR REVISING LIMITS WITH AN EMISSION RIGHTS MARKET
The ability of a private or public emission rights market to provide
for orderly community growth, changes in fuel prices, and variances
37
-------
from original limits in individual cases markedly limits the revisions
in emission limits that would have to be implemented by a control agency.
In fact, only when a major revision is made in the land use plan would
a control agency have to consider the global redefinition of the EDZ
limits. That redefinition would be accomplished according to the cri-
teria outlined previously. Otherwise, one could depend on the market
to provide for the requisite changes in emission limits, subject to the
following criteria.
1. Emission rights may be transferred between parcels within
a component area on a one-to-one basis.
2. Emission rights may be transferred between parcels in
different component areas according to the following
relationship:
the smallest over all i (A)
where Ae is the change in emissions in the new component area
Ae is the change in emissions in the old compartment area
X
t , is the transfer coefficient from source x to
x' receptor i
t ^ is the transfer coefficient from source y to
receptor i.
This is based on the identity,
lAey fcy,i I " lAex tXfll
In other words, there is some increase in emissions Ae
in component area y, whose resulting increase in con-
centration at receptor i will exactly balance the reduc-
tion in concentration at receptor i due to the decrease in
emissions at x. If, over this set of emission increases,
S(Ae ), the smallest value of Ae , the amount, AS ,
v y ' y' y
is permitted to be transferred to component area y,
then concentrations at all receptors will not increase
due to the transfer. The corresponding Aex is sub-
tracted from the old component area.
38
-------
3. After a transfer between component areas according to
criterion 2, the incremental decrease in concentrations
at receptors other than the binding receptor; i.e., the
receptor i for which Ae is the smallest, would accrue
to the control agency as a reserve.
For each receptor i, except for the binding receptor(s),
this reserve would be
Ri =
A. To facilitate orderly community growth, portions of the
reserve at each receptor (the sum of the initial reserve
that was set aside when emission density limits were first
set and any increment due to transactions of emission
rights between component areas) should be auctioned in a
public market at appropriate time intervals. For example,
50 percent of the remaining reserve could be sold every
10 years or after a major change in the transportation or
land use plans.
A control agency could publish a matrix of minimum values of the
tv j/t . ratios. A source wishing to buy emission rights then could
x, i y,i
search for high values in the appropriate row (or column) to find the
potentially best bargin. On the other hand, a source wishing to sell
emission rights could search for high values in the appropriate column
(or row) to find the potentially best prospective purchaser.
39
-------
SECTION V
SPECIFICATION OF DECISION MODEL
The selection of a decision model for setting emission density limits
in an AQMA must rest with the individual agency implementing an EDZ main-
tenance strategy. The choice of a particular model will represent the
resolution of the tradeoff between a specific agency's concerns and avail-
able resources. Four candidate formulations are described in the follow-
ing sections. These formulations differ in their objective functions;
i.e., exactly what is being maximized or minimized. The possibility of
other objective functions is discussed. Next, several variants of the
first formulation, emissions maximization is described. While each
variant does maximize emissions, they each treat the constraint set
slightly differently.
The simplest model is the maximization of emissions. This formula-
tion is the most inexpensive model formulation to implement in terms of
both data requirements and computer time. Given the community goals
expressed in the land use element of the regional plan, it will allocate
emission rights to component areas in a pattern that maximizes the aggre-
gate emission loading in a region and still maintains ambient air quality
standards. The assumption implicit in this formulation is that the
maximization of emissions is closely related to the maximization of eco-
nomic activity permitted by the NAAQSs. Within the limits of the upper
*
and lower bounds on emission density limits , the maximization of emissions
implicity values emissions from different land use categories in the same
That is, the zero cost and technological constraints.
For example, light industry and commercial land use.
40
-------
component area Identically. That is, if both sources are not at their
upper or lower bound, several alternate optimal solutions are possible
since the transfer coefficients are identical.
More significantly, in the case of one binding receptor (which was
the case in the Baltimore TSP test run and was closely approached in the
other tests), the emissions maximization formulation implicitly values
emissions from different areas of the region, weighted by the reciprocal
of their transfer coefficients as having equal value. An example will
clarify this condition. Assume two emission sources, whose emission
density limits, in tons per year per acre, are E.. and £„. Their transfer
coefficients (X/Q) to a certain receptor are, respectively, 0.5 and 0.1
ym/m3/ton/year and their areas of the emission sources are both 2 acres.
Mathematically, the linear program may be stated as
Maximize emissions, E = 2 E..+ 2 E- (7)
such that (0.5) (2) EX + (0.1) (2) E2 <. S
where S is the desired NAAQS. In the maximization of emissions then, the
linear program will implicitly value 5 tons of emissions from source 2 as
equivalent to 1 ton of emissions from source 1 for the purposes of meeting
the air quality constraint. The linear program would assign a higher
emission density limit to the second source regardless of the control
cost schedules (or other costs) of the two emission sources. The remain-
ing three decision models attempt to value the emissions from various
sources according to the location and type of emission source.
The premise behind the control cost minimization formulations (the
second and third models) is that for two emission sources with equal
transfer coefficients (x/Q) to a critical receptor, the emission source
with the higher incremental control costs should be required to meet a
less stringent standard. The goal of this formulation is to set a pattern
of emission density limits such that, if the region is developed to the
41
-------
levels specified in the regional land use plan, the individual source's
incremental control cost of reducing the ambient air concentration at a
critical receptor would be identical with that of all other sources. With
suitable constraints, the minimization of control costs is equivalent to
maximizing emissions weighted by control costs. It is an attempt at a
cost effective approach to setting EDZ limits.
The final model is the maximization of emissions weighted by land
value (model no. 4). There are several reasons why it may be preferred
to the second and third models. First, land value represents the upper
bound on control costs, and in some situations may be easier to obtain
than other control cost data. An individual source is faced with several
options in complying with an EDZ ordinance. One, of course, is the instal-
lation of control equipment in addition to that required by other existing
regulations. Another possibility is the reduction of the emission density.
This means, in effect, that some of the site must remain comparatively
underutilized (not developed to the firm's desired emission generating
capacity). The cost of the underutilized land, therefore, is one form of
control cost facing the firm. However, because there appears to be no
clear theoretical justification for using maximum control cost — land
value — as opposed to average or minimum control cost in the model, some
users may prefer the second or third model. Second, because land costs
tend to be higher in the center of a metropolitan area, this formulation
will not exhibit the tendency of the emissions maximization formulation
to set high emission density limits on the urban fringe. (Model No. 1
does this because the critical receptors tend to be in the center of the
metropolitan area.) Even though the land use plan may restrict intensive
development on the urban fringe, the emissions maximization model, as
opposed to model no. 4, will encourage near suburban industrial develop-
ment to the detriment of expansion of existing core city industrial
facilities and redevelopment of core city industrial sites. Third, model
no. 4 may reduce the potential desruption of existing land use patterns
caused by the introduction of EDZ. The relative value of a parcel is
taken as a partial indicator of that parcel's rightful share of the bundle
42
-------
of emission rights available to the region. High value property would
receive higher emission limits than low value property, other things
being equal. Model no. 4, therefore, would perturb the already existing
value ranking of land as determined by market forces the least, and pre-
sumably change the location and timing of future development the least.
This is done at the expense of the total emissions in the region.
Aside from the four models described below, others are possible. The
user is limited only by the resources and ingenuity of the implementing
agency. The models would be basically the maximization of emissions
weighted by some parameter; e.g., control cost, land value. Within this
format, other formulations can easily be specified if they more adequately
reflect the goals of the community and if the requisite data are available.
The four models are presented below. Variable definitions are shown on
Table 1.
MAXIMIZATION OF EMISSIONS (MODEL NO. 1)
Maximizing emissions is relatively the most straightforward formula-
tion of-an objective function. The most specious relationships to specify
are the SO -PM feasibility constraints, which ensure that a very stringent
X
emission density limit for one pollutant and a relatively lax unusable limit
for the other pollutant are not set in the same cell. The linear program
may have to be run several times as this approximation is fine tuned. The
necessity of this constraint and the estimation of its parameters are dis-
cussed in Appendix C. The formulation of Model No. 1 is as follows:
Maximize emissions = E Z Z £ E. . . _ a.« i, _ (8)
j k Am j'k'*'m j'k'm
43
-------
Table 1. DEFINITIONS OF VARIABLES USED IN SECTION V
SUBSCRIPTS
i = receptor index
j = grid square or point source index
k = land use class index
Si = pollutant species index
1. particulate matter
2. sulfur oxides
m = current and growth index
1. existing sources
2. future sources
n = 1,2 . . . for each potential emission rate (tons/year)
PARAMETERS
ai k m = t*ie lan(* area of emission source j,k,m for existing sources
' ' (m = 1) , obtained from base year land use inventory, and
for future sources (m = 2) , obtained from land use plan
^1 i ~ t*ie Pr°Jected concentration at receptor i of pollutant i
' due to background
c, . = the incremental control cost incurred to reduce emissions
' 'm from land use category k. Obtained from execution of IPP.a
c. = the cost of land for land use type k in grid square j; or for
3» point source j
i
c = the control cost incurred at source j to emit at the emis-
' sion rate indexed by n, i.e., at emission rate e
J j*1*11
ei o n = for each source j, the set of possible emission rates,
J' ' indexed by n
e. , ^ = the maximum feasible future emission density after application
of existing and planned regulations for existing sources
(m = 1), obtained from AQMP projection; for future sources
(m = 2) , obtained from floor area ratio and emission factors
e . .• . = the minimum feasible design year emission density (defined
j,fc,J6,m Qnly for m . 2)
. = the base year emission density (defined only for m = 1) , ad-
J. . >m justed for current and pending emission control regulations.
E. , ,, = the emission density limit, the variable that is set by the
j,K,x,,m settlng model
44
-------
Table 1 (continued). DEFINITIONS OF VARIABLES USED IN SECTION V
f. » = the projected concentration due to fugitive emissions for
' receptor i (usually defined only for £ = 1)
s. ,, = the projected concentration at receptor i of pollutant H
' due to mobile sources
r. ., = safety margin or reserve
p. ^ = the projected concentration at receptor i of pollutant £
' due to power plants and other sources not considered in the
emission limit model
t. . ....... = the transfer coefficient, the incremental air quality impact
' at receptor i of source j, land use class k, for pollutant £.
Obtained from AQDM
p^, v, = coefficients relating the minimum possible levels of particu-
late matter emissions for a given sulfur oxide emission level,
and the converse
X. „ = the reverse calibrated standard, pollutant £, receptor i
1, X,
Z. „ = a zero-one variable
j > ^j n
The Air Quality Implementation Planning Program (IPP) determines control
costs for point sources given emission limits. See TRW Systems Corp.
Air Quality Implementation Planning Program. 2 volumes. Environmental
Protection Agency, Research Triangle Park, North Carolina. November,
1970.
45
-------
subject to:
nlr quality constraint: I Z t.
£ (E
for all
* r
i,£
,
attainment constraint:
m | for all j,k,£; for m = 1
proximity constraint:
- E.j+A kj£ - 0 | for certain A and j
zero cost constraint:
for all j,k,£,m
technological constraint: E
> e°in „ | for all j,k,£,m
— j,K,x,,m
SO -PM feasibility
constraint: E.
-Pk
1 vk
for
PM in
each first term
for £ = SOX in
each second term
for all j,k,m
The objective function is the sum of the total emissions in the
region; i.e., the sum of the products of each source's emission density
limit and its area. The first constraint is the air quality constraint.
The sum of the products of emissions at each source and its transfer
coefficient represents the ambient concentration at the specified receptor
due to sources controlled by EDZ. This must be less than the NAAQS after
subtracting the concentration due to mobile sources, fugitive dust, power
plants, background, etc. There is a separate constraint equation for each
receptor and pollutant.
46
-------
The attainment constraint indicates that an existing source shall not
have an emission density limit assigned to it that is less than its exist-
ing emission density (after considering compliance schedules and other
planned emission controls). This places the maintenance burden on future
sources and capacity expansion of existing sources. The possibility of not
doing so is discussed in Appendix E. There is a separate constraint
equation for each existing source.
The proximity constraint allows one to require emission density limits
in adjoining component areas to be identical. It permits one to set
identical emission density limits for all the grid cells in a political
subdivision; e.g., municipality, township, to facilitate the administration
of EDZ.
The zero cost constraint specifies an upper bound on the emission
density limit for each source. The upper bound is the maximum emission
density that may be expected for that source in the absence of EDZ. It
implies that no control equipment beyond what is already required by
other regulations would be installed; hence, it is the point where the
firm incurs no additional costs due to EDZ.
It is necessary to include an upper bound on all sources; otherwise,
the linear program will assign an excessively high and unusable emission
density limit at a few sources that have the smallest transfer coefficients
•ff
to the critical receptors. In similar manner, the technological constraint
specifies the minimum emission density for each source that is technologi-
cally or economically feasible.
*
This occurred by accident in preparing the Louisville TSP test case.
Several sources had their upper limits omitted. The linear program maxi-
mized emissions by setting most sources at their lower limit and assigning
an emission density limit of 2456 tons/day/km^ at one source that had no
upper limit! In the same run, several sources also had their lower limits
omitted by mistake. Their emission density limits were set at zero.
47
-------
MINIMIZATION OF CONTROL COSTS (MODEL NO. 2)
The minimization of control costs is formulated as follows:
Minimize C = I. I. I c, „ 1 a. . eTT 0
k*m k'£'m j j'k'm \J>k>*>m
subject to:
air quality constraint: Z Z t Z (E
| for all i.X,
attainment constraint: E. , . > e. , I for all j,k,&; for m = 1
j,K,x.,m j,K.,x,,m
proximity constraint: E. , 0 - E..A , g = 0 | for certain A and j
~
zero cost constraint: E..0,m
technological constraint: EJ i. p > e™^ | for all j,k,Jl,m
It should be noted that this is equivalent to maximizing emissions
weighted by control costs, viz.,
Maximize E' = Z Z Z c,, „ „ Z a. ,, m E. . „ ,iriN
k H m j J.k*^1" (10)
subject to the same constraints.
The major difficulty in using this formulation, as well as the next
one (Model No. 3), is the specification of the control cost parameter,
GV o ™» fc^e incremental annualized cost per ton per day of an emission
iiC y X* y ni
48
-------
reduction at each source. A method for specifying the parameter is dis-
cussed in Appendix D. However, the estimation of control costs of future*
emission sources is tenuous at best. If, in a given application, they
cannot be estimated with some confidence it may be prudent to use the
simpler emissions maximization formulation.
ALTERNATIVE FORMULATION OF CONTROL COST MINIMIZATION (MODEL NO. 3)
Rather than forcing control cost to be a linear parameter, one can
recognize its nonlinear and discrete nature. Through the use of inte-
ger programming, one can directly minimize control costs, control costs
being defined as a set of costs associated with each of the possible
emission levels. The formulation is as follows:
Min C' = Z Z Z c. 4 n Z. £ n (11)
j H n
subject to:
= 1 I for
single emission level
per source
for
air quality constraint: Z Z t. . „ e. • Z. i
j i»J»* j>*»n j,£,n I
j n ......
- Xi,£ " bi,£ ~ fi,* " ri,£ - Pi,£ - Si,£ I
for all i,£
Z is a zero-one variable which is constrained to add to one for each
source. Therefore, for the set of Z's for a single source and pollutant,
in the solution to the linear program all Z's except one will be zero.
49 '
-------
The emission density associated with the nonzero Z is the limit for that
source and pollutant. In the case of the sample point source for Balti-
more discussed in Appendix D, the set of emission rates (not emission
densities) and associated control costs are shown in Table 2.
The advantage of this formulation is that it allows the emission
density limits to be set at specified levels that are related to certain
control devices, as was not the case in the prior control cost formulation.
In most contexts, we believe this is not that great an advantage and does
not outweigh the higher per run cost associated with integer programming
(as much as 10 times that of linear programming).
LAND VALUE (MODEL NO. 4)
The minimization of the value of additional land required to be pur-
chased in order that the same level of economic activity be maintained
is not feasible, viz.,
max
min Z E E c^
j k 4
(12)
since one cannot take the reciprocal of the decision variable, E, in the
objective function in an ordinary linear program.
An alternative formulation is the minimization of the cost of excess
land that must be owned by a firm to comply with EDZ, viz.,
... _ „ _ r land max land
jk'm Cj'k aJ.k.* eJ.k»*.* " 'j.fc
max
50
-------
Table 2. CONTROL COST AND EMISSION SCHEDULE,
EXAMPLE OF PROTOTYPICAL HEAVY
INDUSTRIAL POINT SOURCE
Cost
dollars
0
6,876
13,752
20,628
22,769
24,910
27,051
39,122
51,193
63,264
75,410
87,556
99,702
122,203
135,116
157,617
170,530
193,031
205,944
299,791
393,638
487 ,485
581,322
675,179
Emissions
(64 J> >
J»i»n'
tons per year
6,852
5,952
5,052
4,152
3,927
3,702
3,477
3,081
2,685
2,289
1,959
1,629
1,299
1,224
939
864
579
504
219
193
166
140
113
86
Source: Derived from data in
Appendix D
51
-------
Equivalently, one can maximize emissions weighted by the cost of
land, viz. ,
maximize E E E E c.&? a. , E. , „ ,,,^
J' J.k.m "j.k.Jl.m (14)
subject to:
air quality constraint: E E t. . , E (E * a )
j k i.J.k,* m j,k,£,m j,k,ni
| for all i,£
i
attainment constraint: E. , . > e. u „ „ | for all j,k,«,; for m = 1
j,K.,x,,m — j,K,x,,ni
proximity constraint: E, . - E . - 0 | for certain A and j
K.x.m ' J
zero cost constraint: E. , « <_ e!?8? « _ | for all j,k,£,m
technological constraint: E-i k £ m — ^^ Si \ for a11 J»k»&»m
and perhaps the SOx~PM feasibility constraint.
VARIANTS OF MODEL NO. 1
Variant 1A - Municipal Aggregation
The proximity constraint allows one to require the linear program to
set emission density limits such that the limits in certain specified
component areas (grid squares) are equal. Where this is done to a
limited extent, model no. 1 is the most straightforward approach. However,
one may desire to use the proximity constraint more extensively. For ex"
ample, it may be desirable for emission density limits to be assigned by
municipality; i.e., a municipality would be assigned one emission density
limit which could then be allocated among land use categories. If there are
-------
100 municipalities in a urban area and several component areas in each
municipality, using the proximity constraint of model no. 1 becomes a
cumbersome process.
An alternative approach is to define a composite transfer coefficient
for the municipality. After calculating the transfer coefficients for each
component area, one would determine the average transfer coefficient for
each municipality, or preferably an average of the component area transfer
coefficients weighted by the area of that land use in the component area.
Specifically,
aj,k>m
a
j km J
for all !,«, (15)
where the summation is taken over all component areas in the municipality
and the area, a. , is the area of the land use in that portion of the
J ,tc,m,
component area that is also in the municipality. Model 1 can then be used
directly, where j, instead of representing a grid square index, becomes a
municipality index. Note that in this variation, the proximity constraint
could be redefined to allow one to set identical emission density limits
in adjacent municipalities.
The use of this variant is most suited to an urban area that has many
small suburban municipalities; e.g., Boston or Minneapolis-St. Paul.
Even in these areas, it would be advisable to divide the core city into
several sectors or to retain the component area gridding in the core city.
Variant IB - Census Tract Component Areas
While the use of grid squares as the basic geographical unit for com-
piling data is well suited for the requirements of dispersion modeling,
the use of a square grid is both an awkward basis for compiling the emis-
sion and land use inventories as well as for administering a set of emission
53
-------
density limits. If an approach such as variant 1A, municipal aggregation,
i.a conteii,, uted, even the accuracy of the dispersion model output is com-
promised when the average municipal transfer coefficient is calculated.
In addition, an unwieldly data item should be necessary; viz., the amount
of each grid square that is in each municipality.
An alternate approach is the use of census tracts as the component
area. This allows one to compile a more accurate emission and land use
inventory in less time, as well as simplifying the use of variant 1A since
census tracts generally follow political boundries. The only difficulty:
with this approach is the question of dispersion model accuracy. We are
unaware of the answer to this question. It should be noted that AQDM
has been successfully run several times in this manner with the coordinates
of the centroid of the census tract used as the center of a pseudo grid
28
square.
Variant 1C - Constant Emission Density Limits
The logical extrapolation of 1A would be to set one emission density
limit for each land use category in the region. That is, the proximity
constraint would be used for every component area. The advantages of this
*
approach are significantly smaller linear programming costs and, later,
a much simpler task of administration. It will definately allow less
total emissions in the region, and therefore potentially constrain growth
to a greater extent. •
Though possible to simply add a proximity constraints for each grid
square and land use, it is equivalent to reformulate model no. 1 without
the grid square and point source index j, viz.,
Maximize emissions = E L E E, „ E a. , (16)
kt£,m k**'m j j'k'm
About 10 percent of the costs of executing model no. 1.
54
-------
subject to
air quality constraint: £ £ £ t.
4 ' k m j
E,
k,
~ P
i,£ i,
for all i,X,
and corresponding attainment, zero cost, and technological constraints.
One problem with this approach is the definition of the upper and
lower bounds. For future sources (m=2) it presents no problem. There are
however, at least three options on how current sources (m°l) may be
treated
1. Use the same upper and lower bounds for current sources
as future sources. While this is the simplest approach,
it will certainly lead to many nonconforming uses where
the emission density limit for a current source is less
than its present emission density. This will be espe-
cially true of large point sources.
2. Set a constant emission density limit for all area
sources, as in equation (15), but treat each point source
separately as in the original formulation, equation (8).
This should lead to only a moderate increase in program-
ming costs and eliminate most of the nonconforming uses.
3. Set a constant emission density limit only for future
area sources. Existing area sources and point sources
would be treated individually, with unique upper and
lower bounds. This will eliminate all nonconforming
uses, but it will only eliminate about half of the
linear programming costs and cause the administration
of the emission density limits to be almost as complex
as the original emissions maximization, Model No. 1.
55
-------
SECTION VI
GUIDELINES AND MODEL IMPLEMENTATION CONSIDERATIONS
This section first reviews emission inventory considerations and
technical considerations of the linear program. Subsequently, a guide
for executing an EDZ limit setting model is presented and illustrated
using test city of the AQDM manual."*
All of the formulations of the emission density limit setting model
require significant amounts of information on base year and projected de-
sign year emissions. These data must be summarized in a way that usually
would not be available if the emission inventory is prepared without the
requirements of emission density zoning in mind. However, if sufficient
foresight is used when the emission inventory is prepared, the inventory
can meet the requirements of the emission density limit setting model
while incurring a negligible increment in cost.
The following considerations should be observed.
• The grid square network used for the subcounty allocation
of emissions should also observe the criteria in Section II
for partioning an air quality maintenance area. In
practice, the grid system that one would devise using
CAASE^ will require only minor modification. Most of
the modifications would be implemented in the absence
of EDZ in the interest of accurate dispersion modeling.
The exception to this would be modifications resulting
from the criterion of subdividing grid squares spanning
jurisdictional boundaries.
56
-------
The emission inventory must be prepared using the source
classes for which emission density limits are contemplated.
Most emission inventories use only three stationary area
source classes* i.e., residential, commercial, and indus-
trial; it is probable the emission density zoning regula-
tion would utilize a more detailed classification scheme.
Since the emission density limit setting model also re-
quires the base year and planned land use of each source
category in each grid square, this imposes no additional
data requirements. The land use data can be conveniently .,
used for the subcounty allocation of the emission inventory.
In summary, the gridding of the AQMA must meet the additional criteria
discussed in Section II, the land use and emission inventories (both cur-
rent and projected) must have the same classification scheme, their geo-
graphical basis must be the gridded component areas, and the classification
scheme must be the same as (or more detailed than) the categories contem-
plated for the administration of EDZ.
RESOURCE REQUIREMENTS
If the emission inventory is prepared with the requirements of EDZ
in mind, it will require a negligible increase in resources committed
to its preparation.
If one attempts to use the emission density limit setting model with
a pre-existing inventory, the inventory must be modified to reflect the
cited considerations. This will likely require the splitting of grid
squares and a disaggregation of the emission source classes. The methods
one is prone to use will probably lead to a more inaccurate inventory than
if the emission inventory had originally been compiled with EDZ in mind.
It will also require a moderate amount of time and resources; e.g., it
required 1 man-month and an hour of computer time to modify the inventory
for the Baltimore example discussed in the subsequent section. Finally,
in modifying an existing inventory, care must be exercised that the
resulting emission inventory is compatible with the land use inventory.
57
-------
In summary, the compilation of data for input into the linear programs
will be simpler and involve fewer problems if the emission inventory is
developed being cognizant of the requirements of emission density zoning.
The land use inventory and plan will, in most cases, already exist.
The resources required to summarize it on the geographical basis of the
component areas will vary widely, depending on the form and scale of the
inventory and plan. A typical situation is where the inventory is shown
on a series of Prismicolor maps. One common approach is to prepare an
acetate overlay of the grid system and, by planimeter or random dot grid,
determine the amount of land use in each grid. It should also be noted
that if the land use inventory has been summarized by census tract, con-
sideration should be given to using the census tract as the component
area.
Beside the land use inventory and the emission inventory, sufficient
data must be compiled to determine the upper and lower bounds on each
emission density limit. Of the three items, this can require the largest
resource commitment. It is, as noted earlier, absolutely necessary and
its accuracy is, in addition, critical, as many of the emission density
limits may be set at either their upper or lower limits. It will require
several man-months to acquire the necessary information with sufficient
accuracy.
LINEAR PROGRAMMING CONSIDERATIONS
Considering the size; i.e., number of sources, of a typical air qual-
ity maintenance area, it is apparent that for all but the smallest the
linear program to set emission density limits is a significant computer
problem. The computer time required to solve a linear program is diffi-
cult to predict even if the time required for a similar problem is known.
There are, however, several rules of thumb describing the importance of
various parameters of the problem in determining solution time, viz.,
58
-------
"The most important factor is the number of ordinary
function constraints. In fact, computation time tends
to be roughly proportional to the cube of this number,
so that doubling this number may multiply the computation
time by a factor of approximately 8. By contrast, the
number of variables is a relatively minor factor. Thus
doubling the number of variables probably will not even
double the computation time. A third factor of some
importance is the density of the table of constraint
coefficients (i.e., the proportion of the coefficients
that are not zero) because this affects the computation
time per iteration. One common rule of thumb for the number
of iterations is that it tends to be roughly twice the num-
ber of functional constraints."^5
Setting aside Model No. 3, which was not tested here, it should also
be noted that the EDZ problem is an atypical linear program. The matrix
has a density of almost 100 percent. That is, almost every single vari-
able — each of which is an emission 'density limit for a particular
pollutant and source — appears with a nonzero coefficient — the trans-
fer coefficient — in every one of the air quality constraints. (The
attainment, zero cost, and technological constraints are handled as
upper and lower bound statements, and are not included in the matrix.
The feasibility and proximity constraints, if used, would reduce the
density.) A more typical linear programming problem has a density of
20 percent to 40 percent. In other words, more variables are omitted
from each constraint than are left in. General purpose linear program
computer application packages, such as MPS,18>19 are designed for the
more typical problem and can be expected to be inefficient in solving
20
a problem such as that posed by EDZ. A linear program code more
appropriate than MPS for this type of problem could not be located.
The extent of the inefficiency using MPS instead of a special purpose
code could not be estimated.
With these caveats in mind, the Central Processing Unit (CPU) time
and Execute Channel Program (EXCP) count required for the execution step
of MPS are presented in Table 3 for several emission maximization
linear programs. The information processing requirements increase
59
-------
Table 3. PROBLEM CHARACTERISTICS AND MPS TIMING ESTIMATES, EMISSIONS
MAXIMIZATION LINEAR PROGRAMMING
Application
Problem
characteristics
Machine
Rows (Including object)
Columns
Elements
Dens Ity
Iterations
Invert calls
Execute step CPU time
EXCPS
Core
Test run of MPS
on teat city of the
AQDM Manual4
IBM 370/158
18
43
774
71.66
0
0
2.91 seconds
68
25 6K
County example
Air Quality Analysis
Workshop21.b
IBM 370/195
237
121
28,914
34.07
65
1
20 seconds
2,092
400K
Sulfur oxide
example,
Baltimore
IBM 370/158
53
2256
119,568
99.74
1,514
20
20 minutes
23,200
44 8K
Partlculate matter
example,
Baltimore
IBM 370/158
53
2256
119,568
99.74
2,253
22
17.8 minutes
18,950
448K
Sulfur oxide
example #1
Louisville
IBM 370/158
103
629
64,787
86.04
442
25
*•- 3 minutes
~ 1000
448K
Particulate matter
example,
Louisville0
IBM 370/158
103
629
64,787
86.04
419
4
~2 minutes
~ 1000
44 8K
Sulfur oxide
example #2
Louisville0
IBM 370/158
257
629
161,653
71.07
404
5
12 minutes
8556
448K
*Ihe input matrices were all 100 percent dense, however the introduction of a slack variable in every constraint affects the final
density.
Dual problem solved.
CInfe«slble solution reached.
-------
dramatically with the size of the problem; one who is formulating a
linear program for a large air quality maintenance area should exercise
caution so that the problem does not exceed the resources of the user.
If a problem is too large, one effectively has two options for simplify-
ing it, viz., reducing the number of rows (receptors) or the number of
columns (sources). (Reducing the number of feasibility or proximity
constraints would have the same effect.)
To reduce the number of columns, one essentially has to (a) aggre-
gate or remove land use categories or (b) aggregate component areas. A
substantial reduction in the number of columns can only be accomplished
by the former method. Some of the more homogeneous or low emission
density source categories can be removed from the model. For example,
the emission density limits for the residential source categories can be
set a priori and removed as variables from the linear program. Their
contribution to ambient air quality at each receptor would have to be
subtracted from the standard in a manner similar to background or mobile
source contributions. Alternatively, one can collapse or aggregate the
land use categories. For example, the commercial and institutional
categories could be combined. This is feasible either if they have
approximately the same upper and lower bounds on feasible emission den-
sities; or if component areas tend to have one or the other category of
land use, or neither, but not both.
An alternative to reducing the number of columns is reducing the
number of receptors. This is less preferable, because one does not know
where critical receptors will be located. The sensitivity of the optimal
solution to the number and location of receptors has not been completely
investigated. The Louisville sulfur oxides test run was executed for
two receptor sets. A different solution was reached, as discussed in
Section VIII.
61
-------
GUIDELINE FOR EXECUTING THE EDZ MODEL
Presented below are step-by-step instructions for preparing input
for and executing an emissions maximization linear program for one pollu-
*
tant using MPS. Familiarity with MPS is assumed. Data, as described in
Table 4, are assumed to have been compiled:
Table 4. REQUISITE DATA FOR PREPARING AN EMISSIONS
MAXIMIZATION LINEAR PROGRAM
Base year and design year emission inventory for all emission
sources. The stationary area source emission categories must be
identical to those contemplated for regulation by emission density
zoning. The subcounty allocation must be done on a grid system
that is compatible with the component area criteria of Section II.
Base year land use inventory and design year land use plan that
is compatible with the emission inventory in terms of grid network
and subcounty allocation. The land use should be in terms of net
land area, not gross land area; i.e., it should be in terms of lot
area.
Lot area and land use classification of all major point sources.
Meteorological data, air quality, and stack height data necessary
for executing AQDM for the base year and design year.
This includes estimates of stack heights of future sources.
Maximum permitted floor area ratios and minimum feasible floor
area ratios (for nonindustrial development) for each land use in
each component area — obtained from the land use plan, zoning
ordinance, and discussions with the regional planning agency.
The guideline is divided into three sections, viz., right hand side,
columns, and bounds, corresponding to discrete data elements required for
the MPS input.
See the Mathematical Programming System/360 Application Description
for an overview of MPS.
62
-------
Right Hand Side Section
1. Execute AQDM with the calibration option for the base
year emission inventory and the air quality monitoring
stations.
2. Determine the desired ambient air quality standard, in
terms of an annual arithmetic mean for each receptor
to be used in the EDZ model. These may be the long-
term NAAQSs, a controlling long-term standard reflecting
the short-term NAAQS,* or ambient air quality levels
reflecting PSD or local regulations.
3. Adopt a suitable safety margin or reserve at each
receptor. Subtract this from the standard. The reserve
should be based, in part, on the validity of the AQDM
calibration.
4. Reverse calibrate the standard at each receptor in the
EDZ model', i.e., subtract the intercept from the results
of step 3, then divide by the slope.
5. Execute AQDM with the source-receptor file option for the
design year emission inventory.+
6. Using the source-receptor file and the design year emis-
sion inventory, calculate the contribution of each source
not in the EDZ model to each receptor. Sources not in
the EDZ model are those sources in the design year inven-
tory for which emission density limits are not being set.
They may include:
• background
• fugitive emissions
• fugitive dust emissions
« power plants
• mobile sources
Sum the contributions for each receptor.
See Appendix B.
The version of AQDM utilized in this study produces an uncalibrated
source-receptor file. If the version used produces a calibrated file,
Step No. 4 is not necessary. The version of AQDM utilized also produces
a unit emission source receptor file. If this is not the case, the trans-
fer coefficients must be divided by the emission rates.
63
-------
Subtract the vector of non-EDZ source contributions
(step no. 6) from the vector of reverse calibrated
standards (step no. 4). The result is a column vector
that is the right hand side (RHS) of the air quality
constraint.
Columns Section
8. For each plant that is treated as a major source in EDZ;
i.e., it remains in the source-receptor matrix, and has
more than one paint, develop a source-receptor contribu-
tion for each plant by averaging the individual contribu-
tions for each point within the plant weighted by the
emissions from each point.
9. Reduce the source-receptor matrix by deleting each row
that is a non-EDZ source. Take the transpose of this
matrix. AQDM provides, for each receptor, the contribu-
tion from each source; MPS requires, for each source, the
contribution at each receptor.
10. Using the land use inventory, land use plan, and point
source lot area, compute the amount of each component
area taken up by:
a. current land use, by category, excluding point
source land use, and
b. planned growth in land use, by category, excluding
point source land use.
11. Using the results of steps no. 9 and 10 construct a vector
for each emission source', i.e., current area source,
future area source, point source, that consists of:
• the area of that source
• the product of the area of the source and its trans-
fer coefficient (obtained from the transposed source
receptor file) to receptor no. 1
• the product of the area of the source and its
transfer coefficient to receptor no. 2
(and so on for each receptor)
(and so on for each receptor)
(and so on for each receptor)
64
-------
• the product of the area of the source and its trans-
fer coefficient to the ifch receptor.
This is the columns input to MPS; i.e., the constraint
matrix.
12. Determine, as discussed in Appendix C, whether an
SOx-PM feasibility constraint is necessary. If
one is required, the computed parameters are added
to the end of the column vectors and RHS vector.
13. For each proximity constraint required, a 1 and_-l
must be added to the ends of column vectors of the
appropriate sources, and zeros must be added to the
end of the column vectors of all other sources and
the RHS vector.
Bounds Section
14. All necessary input to MPS has now been prepared except
for the BOUNDS section; i.e., the minimum and maximum
feasible values for each emission density limit. The
importance of these parameters is critical; many of the
emission density limits will be set at one or the other
of these values. For existing sources, the lower bound
is their current emission density which is computed by
dividing their total annual emissions by their land
area. For point sources, the lower bound is the smallest
of their current emission density and their projected
design year emission density excluding capacity expansion.
15. The emission density upper bound on existing sources
should be defined from the AQMP projections.^'^ Note
that this is projected growth for existing sources.
16. The emission density upper bound for nonindustrial future
sources is obtained by multiplying the maximum permitted
floor area ratio of the land use in a specific component
area and the land use based emission factor, 2 after
applying existing control regulations.
17. The emission density lower bound for nonindustrial future
sources is obtained by multiplying the minimum feasible
floor area ratio and the most controlled land use based
emission factor.
18. Upper and lower bounds for industrial future sources are
obtained by considering the dirtiest and cleanest indus-
try that can be expected to locate in each industrial land
use category. ....,__ .
65
-------
19. The data must now be formatted for input to MPS. Consult
the following sections for examples of input data, MPS
control programs, .and sample output.
Example of Guideline
An example of the guideline is illustrated for sulfur oxides using
A
test city of the AQDM manual as a sample problem. The test city emis-
sion inventory was modified to show emission classes, as shown in
Table 5. An assumed land use inventory is presented in Table 6.
1. The AQDM execution of test city, as shown on page C-9 of
the AQDM manual, yields a calibration of
y = 0.2 + 0.5132 x
2. Applying the average geometric standard deviation of
1.93, one calculates the annual mean equivalent to
the short-term standard, viz.,
m = 80.9 = 365 (1.93)0*5 ln i'93"2'62
(see equation 25). Therefore, the long-term standard, 80 yg/m3,
is controlling.
3. For the purposes of this example, no safety margin is
assumed.
4. Reverse calibrating,
155.49 = (80 - 0.02)70.5132
5. As shown in the AQDM manual.
6. Source No. 2, assumed to be a power plant, is the only
non-EDZ source. The product 6f its projected emission
rate by its uncalibrated transfer coefficient (x/Q) to
each receptor yielded the second column on Table 7.
7. As shown in Table 7, the vector of non-EDZ source
contributions (in this case, due solely to source no. 2)
is subtracted from the reverse calibrated standard,
yielding the third column in Table 7, the right-hand side.
8. Of the three remaining point sources; i.e., source no. 1,
3 and 4, none have multiple emission points.
66
-------
Table 5. TEST CITY BASE YEAR EMISSION INVENTORY,
SO (tons/day)
Source
no.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
Light
residential
-
-
-
-
0.130 '
0.120
-
0.140
0.094
0.250
-
-
0.252
0.250
- .
0.200
0.260
0.230
-
0.270
Medium
residential
-
-
'-
-
-
-
-
-
-
-
1.244
-
-
-
0.750
0.300
-
-
0.510
0.530
- .
High-rise
residential
-
-
-
-
-
-
-
-
-
-
0.622
-
-
-
-
-
-
-
-
-
-
Conerclal
-
-
-
-
-
-
-
0.017
-
1.244
-
0.378
-
-
-
-
-
-
.
-
Light
Industry
-
-
-
10.500
-
-
0.250
-
-
-
-
0.208
-
-
-
-
-
-
-
-
H««vy
Industry
130.000
-
21.120
-
-
-
0.250
-
-
-
-
0.312
-
-
-
-
-
-
-
-
-
Other
-
150.51
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
67
-------
Table 6. TEST CITY LAND USE INVENTORY (square kilometers)
Score*
Bo.
1
2
3
4
5
6
7
I
9
10
11
12
13
14
15
It
17
18
19
20
21
Base year land use Inventory
Light
residential
15.00
15.54
22.50
15.27
14.06
10.38
14.06
7.38
13.24
13.24
9.00
Medium
residential
44.55
17.50
10.30
17.50
12.25
High-rise
residential
5.46
Cooerclal
0.80
13.64
5.39
Light
industry
0.0809
11.67
7.50
Heavy
industry
0.8094
0.0202
4.94
11.25
Transportation,
open space,
and utilities
1.4165
2.50
2.44
4.99
2.48
2.50
2.50
34.76
6.25
4.49
2.50
7.50
5.70
2.49
2.49
7.50
12. 75
16.00
Total,
excluding
vacant
17.50
17.98
21.60
24.98
18.57
16.56
98.41
25.00
20.26
16.56
25.00
23.58
15.73
15.73
25.00
25. 80
25.00
Vacant
7.50
6.21
3.32
0.0
6.43
8.44
1.59
0.0
4.74
8.44
0.0
0.0
9.27
9.27
0.00
0.00
0.00
Planned land use, design year
Light
residential
-
6.21
-
-
6.11
8.44
-
-
3.12
8.44
-
-
9.27
9.27
-
-
-
Hediua
residential
-
-
-
-
-
-
1.11
-
-
-
-
-
-
-
-
-
-
High-rise
residential
-
-
-
-
-
-
0.14
-
-
-
-
-
-
-
-
-
-
Comercial
-
-
-
-
0.32
-
0.34
-
1.62
-
-
-
-
-
-
-
-
Light
industry
-
-
2.33
-
-
-
-
-
-
-
-
-
-
-
-
-
-
Heavy
Industry
0.0
-
-
0.99
-
-
-
-
-
-
-
-
-
-
-
-
-
-
Ideal.
excltt£l3£
poinc S3UZZKC
25 .CC
21.13
24.«
24.98
25. OC
25.90
vn.se
15.-3O
25. 30
25.90
25.90
23.58
25.08
25.03
25.00
25.08
2v ae
ON
oo
-------
Table 7. CALCULATION OF THE RIGHT-HAND SIDE
Standard
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155.49
155 . 49
Concentration
due to source no. 2
2.58
3.08
4.28
5.58
5.58
1.78
2.18
5.88
8.98
1.58
4.18
11.18
0.88
1.18
1.58
1.58
4.08
Right-hand
side
152.9
152.4
151.2
149.9
149.9
153.7
153.3
149.6
146.5
153.9
151.3
144.3
154.6
154.3
153.9
153.9
151.4
69
-------
9. The transfer coefficients from source no. 2 to each
receptor are deleted from the source receptor matrix.
The transpose is taken.
10. The land use inventory, Table 6, already presents
the information in the required format.
11. For example, the first source has an area of 0.8094
square kilometers, and its transfer coefficients to
the first two receptors are 0.1607 and 0.5067 pg/m3
per ton/day, so the beginning of the first vector
would be
• 0.8094
• 0.1301
• 0.4101
12. The SO -PM feasibility constraint is assumed to be
unnecessary.
13. The proximity constraint is assumed to be unnecessary.
14. The lower bounds for the current sources are their cur-
rent emission density. For example, the first source
has a current emission rate of 130 tons per day and
an area of 0.8094 square kilometers. Its lower bound
would be 160.6 tons per day per square kilometer.
15. For example, based on the AQMP analysis, assume source
no. 1 plans to increase its capacity by one-third.
Its upper limit would therefore be 213.6 tons per day
per square kilometer.
16. For example, assume the maximum permitted development
density for light residential land use is 1/2 acre
lots, or two dwelling units per acre. Multiplying
by the appropriate land use based emission factor,22
assume this product yields an emission density of
0.0197 tons per day per square kilometer.
17. Assume for light residential land use, the minimum
expected development density is 1-1/2 acre lots. If
the emission factor remains the same, the lower bound
would be one-third of the upper bound calculated in
step 17, or 0.00657 tons per day per square kilometer.
70
-------
18. For example, assume the range of expected floor area
ratios for light industrial development is 0.30 to
0.36, that light industrial sources are currently
required to have a 75 percent control efficiency and
the maximum feasible control efficiency is 90 percent.
Also, based on an analysis of the expected industrial
development, assume that the average expected emission
factor is 0.0674 tons per day per square foot of floor
area (658 kg/day/m2). Therefore the upper and lower
bounds would be computed as follows:
(0.3) (0.9) (0.0674) = 0.01821
(0.36) (0.75) (0.0674) = 0.02604
19. The format requirements for MPS are illustrated in
Section VII.
Based on the data developed in this sample program, a MPS input data file
was prepared and executed. The output of this execution is shown in
Figure 3 which includes a listing of the MPS control program. The system
used in this example consists of 8-digits, viz.,
Column Index
1-3 Source number, sources 1 to 4 are point sources
Source number, sources 5 to 21 are area sources
4-6 Unused, all zeros
7 Land use, 2 = light residential
3 = medium residential
4 = high rise residential
5 = commercial
6 = light industry
7 = heavy industry
8 Current or future source, 1 = current
2 = future
The interpretation of the output is explained in Section VII.
71
-------
//GC35EKHR JOL> (GCA02. JOB 34* UNIT=?314,SPACF=tTRK.(tPRI.10>..CONTIC)
IEFIS531 SUBSTITUTION JCL - UN IT = ?314.SPACE=
XXSCRATCH? OD UNIT=?314,SPACE=ITRK,ftPRI,10)..CONTIG)
I!rF6531 SUORTITUTION JCL - UN I T = ?31 4 . SP»CE= « TRK , ( 10 . 1 0 > . . CONT IG )
XXSCRATCH3 DO UMI T=23 1 4 .SPACE =( TRK . »CPR I .1 0 ) . ,CONT 1C)
IFF653I SUBSTITUTION JCL - UN IT=231 4 . SPACE =» TPK . I 10 . 1 0 ) . . CONT IG )
XXICRAfCHI DD UNIT=23I»,SPACE = ( TRK.(CPRI.10)..COVTIG)
IEPf.S3I SunsTITUTION JCL - UN IT =231 A . SPACE= (TRK . ( 1 O . 1 0 ) . .CONT IG )
XXSYSMLCP ID l/NIT = 2314.SPACE=(TRK.(tPRI .10) . .CONT IG ) ,O" i="" sp="<" new.="" pass="" )="" ief653i="" substitution="" jcl="" -="" un="" 1t-23i="" 4="" .="" space^ttrk="" (="" 1o.="" 1="" 0="" .cont="" 1g=""> .DI SP= (NEW .PASS >
XXSYSP9INT DD SYSOUT=A
XXSYSUDUMP OD SYSOUT=A
XXSYSABS DD SYSOUT=A
XXSYSPUNCH OD SYSOUT=R
//MPSI .SYSIN DO »
IFF2361 ALLOC. FOR GC35EKHR MPSI
IEF237I 243 ALLOCATED TO STF.PLIB
IEF237I 340 ALLOCATED TO SCRATCHl
IFF237I 341 ALLOCATED TO SCRATCH2
IEF237I 342 ALLOCATED TO SCRATCH3
IEFP37I 343 ALLOCATED TO SCRATCH4
IEF237I 344 ALLOCATED TO SYSMLCP
IEF237I 6O2 ALLOCATED TO SYSPRINT
IEF?37I 603 ALLOCATED TO SYSUDUMP
IEF237I 60« ALLOCATED TO SYSABS
IEFP37I 73? ALLOCATED TO SYSPUNCH
1EF?37I 70? ALLOCATED TO SYSIN
1EFI4?! - STEP WAS EXECUTED - COND CODE OOOO
IFF2B5I SYS 1.MPS.SYSTEM KEPT
IEFP8SI WL SER NOS= RESPAG.
IEF285I SYS77I2-5.T09I51 I .RVS77.GCSSEKHR.ROOOOOO 1 DELETED
IF.F28!il VOL SER NOS= PPDO03.
I = F?8-SI SYS77125.TOQ151 1 . RV677 .GCS5EKHR . HOOOO 00 ? DELETED
IFF28S! VOL SFO NOS= KEM002.
IFFaSSI SVS77]?S.TO"151 1 ,RV877.GCa5FKHR.ROOOOO03 DELETED
IFF285I VOL SER NOS= BSD003.
IEF2B5I SYS7712S.T091S11.RV877.GCaSEKHR.R0000004 DELETED
IEF2R5I VOL SER NOS= BSD002.
1EF?85I SYS77125.TO91511.RV877.GC35EKHR.ROOOOOOS PASSED
IEF285I VOL SER NOS= MDS444.
1EF373I STTP /MPSI / START 77125.09IS
1EF374I STEP /MPSI / STOP 77125.0916 CPU OM IN 01.63SEC STCR VIRT 25fiK
XXMPS2 EXEC PC,M=EXECUTOR.CONO = ( 4 ,LT .MPSI )
XXSTEPL1B DD DSN=SYS1.MPS.SYSTEM,DISP=SHR
XXSYS«LCP OD UNIT=2314.DSN=*.MPSI.SYSMLCP.DISP-(OLD.DTLETE)
XXPROHFILF OD UNIT=2314,SPACF=(CYL.(tPRI?)..CONTIG)
Figure 3. Test city example output
72
-------
IF.F6531 SUBSTITUTION JCL - UN 1 T=231 4 . SPACE=(C YL .( 5 ) . . CONT 16 )
XXMATRIXl ID UNIT=231«.SPACE«=(CYL.(EPOI2>..CONT?G>
IEF6531 SUBSTITUTION JCL - UN IT=231« . SPACE=(C YL .( 5 > . .CONT IGI
XXETA1 OD UNIT=2314.SPACE=«CVL.ftPRI2»..CONTIG>
IEF653I SUBSTITUTION JCL - UN IT=231 4 . SP ACE=(CYL .< 5 ) . . CONT 1 G >
XXSCRATCH1 DO UNI T=23I4 . SPACE = < CVL. « CPR I 2 ) . .CONTI G)
IFF653I SUBSTITUTION JCL - UNIT=23I 4 . SPACE=»CYL. < 5 ) • t CONT IG )
XXSCRATCH2 DO UNI T=2314 ,SPACE=(CYL. ( EPf» I 2) . .CONTI G»
IEF653I SUBSTITUTION JCL - UNIT=23I4. SPACEMCVL .< 5 > ..CONT IG)
XXSYSPRINT OD SYSOUT = A
XXSVSUDUMP DO SVSOUT=A
XXSVSABS DD SVSOUT=A
XXSVSPUNCH DD SYSOUT=B
//MPS2.SYSIN DD *
IEF236I ALLOC. FOR GC35FKHR MPS2
IEF237I 243 ALLOCATED TO STEPLI8
IFF237I 34« ALLOCATED TO SYSMLCP
1EF237I 3« 1 ALLOCATED TO PROBFILE
IEF237I 34(S ALLOCATED TO MATRIX!
IEF237I 340 ALLOCATED TO ETA1
IEF237I 343 ALLOCATED TO SCRATCHI
IEF237I 341 ALLOCATED TO SCRATCH2
IEF237I 602 ALLOCATED TO SYSPRINT
IEF237I 603 ALLOCATED TO SYSUDUMP
IEF237I 604 ALLOCATED TO SVSABS
IEF237I 732 ALLOCATED TO SYSPUNCH
IEF?37I TO2 ALLOCATED TO SYSIN
IEFI42I - STEP WAS EXECUTED - COND CODE 0000
1EF285I SYS I .MPS. SYSTEM KEPT
IEF20S1 VOL SER NOS= RESPAG.
IRF2B51 SYS77125.TO9151 1.RV877.GC35EKHR.ROO00005 DELETED
IEF2B5I VOL SER NOS= MDS444.
IEF285I SYS7712S.T09I51 I .RV877.GC95EKHR.ROO0001 I DELETED
lEF?^! VDL SFR NOS= KEM002.
IEF?B«iI SYS77125.T09I51 I . RV877. GC3SEKHR.HOOOOOI 2 DELETED
IEF?B5I VOL SER NOS= PPD002.
IEF28SI SYS77I2S.TOOI5I1.RV877.GCS5EXHR.ROOOOOI3 DELETED
IEF285I VOL SER NOS= PPDO03.
IFF285I SYS77125.TO915I l.RV877.GC3)5EKHH.ROOOObl4 DELETED
IEF285I VOL SFR NOS= BSD002.
IEF2S5I SYS77I25.TO9I5II.RVB77.GC85EKMR.ROOOOOI5 DELETED
IEF28SI VOL SER NOS= KEMO02.
1EF373I STEP XMPS? / START 77125. O9I6
IEF374I STFP /MPS2 / STOP 77125.0917 CPU OMIN 02.9ISEC STOR VIRT 2S6K
IEF3751 JOB /GCaSEKHRX START 77125.0915
IEF376I JOB XGCSSEKHR/ STOP 77125.0917 CPU OMIN 04..54SEC
Figure 3 (continued). Test city example output
73
-------
CONTROL PROGRAM COMPILFR - MPS/3K) V."-«l 1
0001 • PROGRAM
0002 INITIALZ
OO65 MOVE(XDATA.«TESTCITY«)
0066 MOVE
OO67 CONVERT!'SUMMARY•)
O068 BCOOUT
0069 SETUP»'BOUND*."EM1SL«)
O0.7O MOVE(XOBJ.*EMISS'>
O071 MOVEfXRHS.*SO?«)
DOT? PRIMAL
OO73 SOLUTION
OO74 RANGE
0075 EXIT
0076 PEND
EXTCUTOR. MPS/360 V?-MU
CONVERT TESTCITY TO PBF ILF
TIME = 0.01
SUMMARY
I- ROWS SFCT1ON.
0 MINOR CTROR(S) - 0 MAJOR ERROR(S).
2- COLUMNS SECTION.
O MINOR BRROR(S) - 0 MAJOR ERROR(S).
3- RHS'S SECTION.
SO?
0 MINOR ^RROR(S) - 0 MAJOR ERROR IS).
5- BOUNDS RfTCTION.
EMISL
0 MINOR E(?POR(S) - 0 MAJOR ERRORfSI.
Figure 3 (continued). Test city example output
74
-------
FXFCUTOR. MPS/360 V2-M11
NUMBER OF ELEMeiTTS BY COLUMN ORDER
PACE
77X135
19
26
33
4O
47
54
O010O071 .
OO80O021.
01200061.
OI6OOO31 .
00600022.
011000*2.
.. .18
...18
...ta
.. .18
...18
...18
O03OOO71....18
OO9O0021....16
012OOO71 ....13
01700021 ....18
OO7OOO62....1S
01 100052.. ..I 8
OOAOO061 .
OO9OOOS1 .
O130OO2I .
0180 OO 21 .
00 70 007? .
01300022.
. .. 18
... 16
...18
...18
...18
...18
OOSOOO21..
0100O02I..
OI3OOO5I ..
O19O0031..
OO9OOO22..
0130005?..
..18
..18
..18
. .18
..ta
..18
OO60O021 ..
O110O031..
0140O021 ..
02OOO031 ..
OO90OOS2..
OI40O022..
..18
..18
..18
..18
..IB
..18
OO70OO6I.... 18
Ol 1OOO41.. .. 18
O15O0031.. ..18
O210OO21.. ..18
01000022. ...18
O17OOO22....I8
00700071 .. .. 18
Ol lOOOSi .. .. 18
016OOO21 .. .. 18
OO5OOO22.. .. 18
Ol 1OOO32-. .. 18
O1OOOO22.. ..18
Figure 3 (continued). Test city example output
-------
EXECUTOR. MPS/36O V2-MI1 PAGE 3 - 77X125
NUMBER OF ELEMENTS BY ROW ORDER. EXCLUDING BUS'S. INCLUDING SLACK ELEMENT
1 N EMISS ....43 L ROOI ....43 L ROO2 ....43 L ROO3 ....43 L ROO'4 ....43 L RODS .,..43 L ROO6 ....43
8 L ROO7 ....43 L RO08 ....43 L ROO9 ....43 L R010 ....43 L RO11 ....43 L R012 ....43 L RO13 ....43
15 L RO14 ....43 L ROiS ....43 L ROI6 ....43 L RO17 ....43
PROBLEM STATISTICS - 18 ROWS. 60 VARIABLES. T74 ELEMENTS. DENSITY = 71.66
THESE STATISTICS IMCLUDE ONE SLACK VARIABLE FOR EACH ROW.
0 MINOR ERRORS. O MAJOR ERRORS.
Figure 3 (continued). Test city example output
-------
EXECUTOR.
8CDOUT - USING PBF ILE
TIME = 0.24
MPS/360 V2-MI1
NAME
ROW?
N BOSS
L ROOI
L ROO2
L RO03
L BOO*
«. 0005
L "0 O6
L POO7
L Rooe
L ROO9
L R01O
L ROM
L BO I 2
L ROI3
L RO 1 4
L BO 1 5
L R016
U RO17
COLUMNS
OOIOO07I
OO 100071
00100071
OOIOOO71
OO10OO71
00100071
O010OO71
00100071
OOIOO071
00300071
00300071
OO30O071
00300071
O03O0071
00300071
0030007)
00300071
003OOO71
00400061
004OOO6I
00400O6I
004OOO61
00400061
O0400061
00400061
OO400061
00400061
00500021
OOSOOO21
00500021
0050OO?!
00500021
00500021
EXECUTOR
TESTCITY
EHISS
R002
BOO 4
R0a6
R008
B010
ROI2
BO14
R016
EMI SS
B002
R004
R006
BOOB
R01O
B012
BOI4
B016
EHISS
BOO 2
R004
ROO6
BOOS
B010
P012
RO14
BOI6
EMISS
BOO?
B004
B006
B00.8
BO10
MPS/36O V2-M11
.8O940
.41010
. 2463O
.O648O
.18700
.06850
.1O20O
.07357
.O6475
.02023
.OO209
.O1018
.oooes
.OO3II
.OD064
.00230
.OO091
.00183
.08094
.Ol 301
.04932
.OO476
.02028
.00376
.01086
.00598
.0071 1
15.00000
9.09SOO
3.63200
2.322OO
3.36600
1.72200
R001
RO03
ROO5
B007
ROO9
RO1 I
ROI3
B01S
B017
B001
ROO3
BOOS
ROOT
B009
RO1 1
H013
•R015
R017
B001
R003
BOOS
RO07
ROO9
BO 11
R013
BOIS
R017
BOOl
R003
RODS
ROO7
RO09
BO 11
.13010
.451 1O
.1957O
.0676O
.teiae
.1O81O
.059T3
.O709O
.06749
.00199
.00323
.01221
.00112
.OO5O3
.00230
.00062
.00149
.0017*
.OO84T
.03519
.02473
.OO569
.01883
-O1O07
.00367
.00736
.OO647
7.5930O
4.55700
2.877OO
1.886OO
2.51300
1.886OO
Figure 3 (continued). Test city example output
77
-------
EXECUTOR.
MPS/3f-0
OOS00021
O050O02I
005000?!
00600021
OO600021
00600021
00600021
OO600021
00600021
00600021
00600021
0060O021
007000ISI
00700061
00700061
00700061
00700OA1
007O006I
007O0061
00700061
007OO061
O0700071
OO7O0071
00700O71
00700071
0070OO71
00700071
O0700071
00700071
OO7OOO71
OOBOOO21
00800021
008OOO21
00800021
00600021
OOBO0021
00800021
00800021
008000?!
OO9O0021
00900021
009O0021
O09OOO21
O09O0021
009OOO21
00900021
00900021
0090O02I
009O0051
00900051
OO9OOO51
O0900051
OO900OSI
00900051
R012
ROI 4
R016
EM1SS
ROO2
R004
RO06
RO08
ROIO
H012
R014
RO16
EMISS
0002
ROO4
R006
R008
0010
R012
ROI 4
R016
EMISS
R002
R004
R006
R008
R01O
R012
RO14
R016
EMISS
ROO2
R004
R006
ROOB
ROIO
0012
R014
ROI ft
EMISS
ROO2
POO*
ROO6
ROOfl
ROIO
R012
ROI4
ROl*)
EMISS
R002
R004
ROOC.
RODS
ROIO
I.8390O
1.31100
1.25400
I5.5400O
7.87200
4.728OO
1.24300
3.606OO
.19900
.95800
.4123O
.2430O
1 .67000
.876OO
7.109OO
.68600
2.92400
.54130
1 .5660O.
.86220
1.0260O
- 4.938OO
.79400
3.0OBOO
.29O30
1.23700
.22910
.6626O
.3649O
.43400
22.50000
2.324OO
I 1.3200O
.98100
3.4610O
7.1 IOOO
2.781OO
1.00800
2.03400
15.270OO
I. 18300
2.421OO
.50690
1.20200
.47O30
1.72500
.464 1O
1. I 12OO
. - .8O36O
.06228
.I274O
.O266R
.06324
.02475
R013
R015
RO17
9O01
R003
R005
R007
RO09
R01 1
ROI3
R015
RO17
ROO1
R003
R.O05
R007
R009
R01 1
no 13
RO15
R017
ROO1
ROO3
RO05
ROO7
R009
R011
RO13
ROI5
ROI7
R001
RO03
RO05
RO07
R009
ROI I
"013
R015
R017
ROO1
RO03
ROO5
RQ07
RO09
ROI 1
0013
RO15
0017
R001
R003
R005
R007
ROO9
ROI 1
1 .36800
1.1970O
1 .2660O
2.498OO
8.66000
3.7570O
1.29700
3.4B9OO
2.076OO
1.1*700
I .1MOO
1.30000
1 .??000
5.0 C30O
3.5t>*OO
.82020
2.715OO
1.45100
.5297O
.5192O
.93330
.5165O
2.147OO
I.50800
.34710
1.149OO
.61420
.22420
.45030
.395OO
I.766OO
3.501OO
13.58OOO
1 .24200
S.589OO
2.56300
.6930O
1.65400
1.958OO
.9039O
I.56300
7.606OO
-6397O
2.32100
1.I59OO
.29310
.6794O
1.3680O
.O4757
.08229
.40030
.03367
.122^0
.06099
Figure 3 (continued). Test city example output
78
-------
MPS/36O V7-MII
0090O051
00900051
OO90OO51
01000021
010O0021
O1OOOO21
0100OO?!
010OOO21
Ol 0000?!
010000?!
01OOOO?!
O10OOO21
01 IOOO31
01 1OOO?I
01 100031
01 1OOO31
01100031
01 10O031
01 1OOO31
01 10O031
01 1OO031
01 10OO41
01 100041
01 10OO41
01 100041
01100041
Ol 100041
01 1O0041
01100041
Ol 10OO4I
01 IOO051
01 100051
01100051
01100051
01 10OO51
O11OOO51
01 100O5I
01 100OFS1
01 100051
01200061
01200O61
01200061
O12OO061
O12OOO61
012OOO61
012000*1
01200061
0120O061
012OO071
01200O7I
012O0071
01200071
OI20O071
012OOO71
R012
9014
R016
EMISS
0002
R004
R006
ROO6
WO10
3012
R014
ROlft
EMISS
ROD?
0004
R006
R008
ROIO
R012
R014
ROI6
EMISS
R002
ROO4
R006
RO08
ROIO
R012
R014
RO16
EMISS
R002
R004
RO06
Rooe
ftOlO
RO12
P014
R016
EMISS
R002
ROD 4
R006
ROO8
ROIO
ROI2
Pt>14
R016
EMISS
RO02
POO 4
ROO6
ROOS
ROIO
.09080
.02443
.osaso
14.06000
4.137OO
2.51200
7.27300
4.O61OO
2.599OO
2.813OO
?.O85OO
2. 1 15OO
44.55000
I5.56OOO
6.6S10O
7.190OO
13. 6300 O
3.S8100
10.06000
5.15400
6.OOOOO
5.45500
1.90500
.81440
.88150
1.669OO
.43850
1.232OO
.IS311O
.73470
13.64OOO
4.7620O
2.03600
2.20400
4. 173OO
1.09600
3.079OO
1.578OO
1.83700
- 7.50OOO
1.25OOO
1. 12OOO
.59700
3.272OO
3.3300O
1.6870O
.34880
.66780
11.25000
1.8750O
3.93OOO
.8955O
4.90600
. 499SO
ROI3
R015
ROI7
R001
ROO3
ROOS
ROO7
R009
R01 1
R013
R015
RO17
ROOI
RO03
ROOS
RO07
RO09
R01 1
R013
R015
R017
ROOI
R003
ROOS
R007
ROO
-------
FXECUTOP.
MPS/360 V2-W1I
01200O71
012OOO71
01200071
013000?!
013000?!
01300021
01300031
01300021
01300021
0130OO?!
OI30OO2I
013OOO?!
013000S1
0130OO51
Ol 300051
013000S1
01300O51
01300051
013O0051
01300O51
01300051
O14OOO?!
01AOOO?!
O140OO21
OI400O21
OI4OOO21
OI4O0021
014000?!
014000?!
0140O021
01500031
01500031
OI500O31
01500031
01500O3I
01500O31
01500031
015OO031
01500031
016OOO?!
OI6OO02I
016000?]
OI6O002I
016OOO21
OI600O?!
O160O021
0160OO?!
016O0021
01600031
016OOO31
01600031
01600031
01600031
01600O3I
ROI 2
RO1A
R016
EMISS
R002
ROOA
R006
ROOB
ROIO
ROI 2
R014
B016
RMI51S
R002
0004
ROO6
»OO«
OO1O
ROI 2
R014
R016
EM1SS
ROO2
ROO4
ROO/S
ROO8
ROIO
ROI?
POI4
R016
EMISS
ROO?
R004
R006
ROOS
ROIO
R012
ROJ4
ROI 6
EMISS
ROO2
ROO4
ROO6
ROOB
ROIO
ROI?
RO14
R016
EMISS
ROO2
ROO4
ROO6
ROOB
ROIO
2.83100
.97310
I.302OO
1O. 36000
I.3O6OO
2.T74OO
.747TO
1.98700
.41S4O
1.Q19OO
.39T7O
.05640
5.36500
.6774O
1.43800
.38770
I. 03100
.2154O
.99510
.2062O
.4959O
- I4.06OOO
2.83600
1 . 1150O
5.81100
2.51300
7.27300
3.212OO
4.19500
3.76500
- 17.50OOO
2.03400
4.17600
.2730O
6.09QOO
1 .3930O
9.81400
1.034OO
2.72700
7.37500
.86O7O
1.58000
.48970
1.66200
.44S4O
3.218OO
.32670
.59300
10.50000
1.22500
2.25000
.69720
2.36700
.63840
ROI 3
R015
9017
R001
R003
ROOS
ROO7
RO09
ROI 1
R013
R015
ROI7
RO01
ROO3
ROOS
ROO7
ROO9
ROI I
RO13
RO15
R017
ROO I
ROO3
ROOS
ROO7
R009
ROI 1
ROI3
ROI 5
ROI 7
R001
ROO3
ROOS
ROOT
R009
R011
ROI 3
RO1S
R017
ROOI
QOO3
ROOS
ROOT
R009
ROI 1
R013
ROI5
R017
ROOt
ROOS
ROOS
RO07
ROO9
ROI I
.3A340
.117300
1.409OO
.8 1520
2.0S40O
4.2950O
.694 TO
5.3720O
.«876O
.38730
.S7?20
I.4230O
.4??TO
1.O6SOO
2.22700
.36O2O
2.78500
.5I21O
.20080
.296 TO
.737TO
3.9TOOO
1.54400
I .2ISOO
2.69000
2.018OO
4.05600
2.6O200
3.9ISOO
2.81300
2.03400
3.7S20O
2.979OO
1.93700
4.3?eOO
6.81800
.7735O
1.40700
4.4OIOO
.71460
1-2070O
1.73TOO
.6018O
2.57100
1.188OO
.24930
.43660
1 . 15OOO
1.01700
1.718OO
2.501OO
.85660
3.660OO
1.69200
Figure 3 (continued). Test city example output
80
-------
EXECUTOR.
MPS/360 V2-M11
01600031
OI6OOO31
01600031
017000?!
01700O21
01700021
OI7OOO21
017000?!
01700O21
017OOO21
01 7OOOPI
0170OO21
01600021
OIBOOO21
01800031
018OOO?!
0180OO?!
01S00021
oisooo?!
OI80O021
0180OO21
01900O31
OI900O3I
01900031
01900O31
O1900O31
OI900O31
OI9OOO31
OI900O31
O1900O31
02000031
02000O31
02OOOO31
0200OO31
0?OOOO31
02000031
02000031
07000031
O200OO31
O2100021
02100021
0210O02I
021000?!
07100021
O2I00021
02100O21
OP10002I
02100021
005O0022
005O002?
OOSOOO2?
OO50OO22
00500022
OO50O022
R012
RO14
R016
FMISS
R008
R004
ROO6
nods
R01O
RO1?
R01A
R016
6MISS
ROD?
R004
ROOft
ROOS
RO1 0
R01?
RO14
ROlft
EMISS
RO02
ROO4
R006
ROO8
RO10
RO12
R014
HOlft
EMISS
R002
RO04
R006
ROOS
ROtO
R012
R014
RO16
EMISS
R002
ROD 4
0006
a oo e
RO10
RO12
R014
R016
EMISS
ROO2
BO04
R006
R008
ROIO
4.5S2OO
.46520
.8442O
13.24000
2.20600
2.83600
3.73800
1.0480O
5.469OO
I.9O100
«. 79900
3.8240O
I3.24OOO
2.7S70O
1.239OO
3.362OO
1.45300
3.53600
2.3640O
6.8480O
4.80200
17.50000
3.442OO
2.912OO
3.39500
3.53000
3.46500
3.08700
3.3S5OO
t 1.64OOO
12.25000
2.04200
2.62500
1.47OOO
3.5340O
1.47000
3.67SOO
1.597OO
6.2400O
9.0000O
.74970
1.38600
.72000
1.90400
.72000
3.06000
.783OO
1.52800
7.SOOOO
4.S4900
1.8I6OO
I.1610O
1.6830O
. .»61OO
R013
R015
"017
R001
R003
ROOS
ROO7
ROO9
'JO 11
ROI3
R015
R017
9001
R003
ROOS
RO07
R009
R011
ROI3
R015
R.OI7
R001
R003
ROOS
ROOT
RO09
RO1 1
ROI3
R01S
R017
ROOI
ROO3
ROOS
ROO7
R009
R011
R013
R015
t«OI7
ROOI
R003
ROOS
R007
R009
ROI 1
R013
R015
0017
ROOt
R003
ROOS
R007
ROO9
RO11
.35*90
.62I6O
1 .6370O
2.7600O
1.2360O
.70940
1. R6000
1 .14100
2.36IOO
7.997OO
4.800OO
3.024OO
?.«06OO
.7346O
2.2290O
1 .05IOO
2.3290O
2.54000
8.800OO
3.82400
2.82600
3.6400O
1 .6380O
3.332OO
1.92OOO
5.O94OO
2.I74OO
9.056OO
6.34600
I .884OO
2.45000
2.227OO
1 .93IOO
2.52700
5.063OO
.94200
2.40200
8.16700
.69210
I .260OO
1 .63600
.92610
2.14300
2.040OO
.69210
.882OO
4.00OOO
3.79700
2.27900
I.A39OO
.9428O
1.256OO
.94280
Figure 3 (continued). Test city example output
81
-------
MPS/3M) V2-M11
00500022
00500022
O050OO22
00600022
00600022
O060OO22
00600022
OO600O22
OO60O022
006OOO?2
00600022
006OOO22
OO7000I42
O0700O62
00700062
O070006?
00700062
O0700062
0070O062
00700062
007O0062
0070O072
O0700072
OO70O072
007OOO72
OO700072
O07OO072
007OOO72
00700072
00700072
00900022
O09OOO22
OO9OO022
OO90OOP2
009OO022
009O0022
00900O22
0090OO22
009O0022
00900052
00900052
009O0052
009O0052
00900O52
00900052
0090O05?
009O0052
00900052
010OOO22
01000O22
010O0022
01000O?2
01000022
01000O??
ROI2
ROI4
R016
EMISS
R002
«004
ROO6
R008
0010
R012
R014
R016
EMISS
R002
R004
ROO6
R008
ROIO
R012
R014
R016
EMISS
R002
R004
R006
RO08
P01O
ROI2
R0t4
9016
EMISS
R002
ROO4
R006
R008
ROIO
RO1?
R014
R016
EMISS
R002
R004
R006
ROO8
ROIO
R012
R014
R016
EMISS
R002
R004
ROO6
R008
ROIO
.91950
.65550
.627OO
6.21000
3. 147OO
1.88900
.49680
1.44IOO
.47940
.78250
.56450
.4968O
- 2.33300
.37510
1.42100
.13720
.5846O
. 1O83O
.31310
.17240
.2O510
.98750
.15880
.6017O
.0581O
.24750
.04582
. 1325O
.07298
.08680
- 6.1O7OO
.47330
.9686O
.20280
.48060
.18810
.69010
. 1 857O
.4446O
.32140
.02491
*O5O97
.01067
.02529
.00990
.03632
.O0977
.0234O
- 8.43800
2.46700
1.50700
4.364OO
2.43700
1.559OO
R013
ROIS
R017
R001
RO03
R005
R007
RO09
P011
R013
R015
0017
ROOI
R003
RO05
R007
R009
R01 I
R013
ROIS
R017
R001
R003
RODS
ROO7
ROO9
ROII
ROI3
ROIS
RO17
ROD I
R003
RODS
ROO7
R009
R01 1
R013
ROIS
RO17
R001
R003
ROO5
R007
R009
R011
RO13
901 S
ROI7
R001
R003
RODS
ROO7
0009
ROII
.6840O
.S9850
.63300
.99860
3.46100
1.5O2OO
.5185O
1.39500
.82970
.4S83O
.34400
.51980
.244OO
1.01400
.7I87O
.16400
.54890
.29020
.10890
.2128O
.I 8660
.10330
.42940
.30170
.06942
.22980
.12280
.O4463
.09006
.O79OO
.36150
.62340
3.04300
.25590
.93130
.46350
.11730
.27180
.54780
.01903
.03291
.16010
.01347
.04901
.08439
.02931
.01430
.02800
3.492OO
1 .4680O
1.2I2OO
1.78SOO
1.92700
2.237OO
Figure 3 (continued). Test city example output
82
-------
EXECUTOR.
MPS/360 V2-N11
0100O022
0100002?
01OOO022
01 100032
01 100032
01100032
01 100O32
01 100032
01 100032
01 100O32
01 100032
Ol 1OOO32
Ol 10004?
Ol 100042
01 10004?
01 10004?
01 10004?
01 10004?
01 100042
01 IOOO42
01 10O042
01 IOOO1?
01 IOOOS2
01 1OOO52
01 10OOS2
01 100 OS2
01 100OS?
01 100052
01 1OOO52
01 100052
O130OO??
01300O22
013O0022
01300O22
0130002?
01300022
013000P2
01300022
0130002?
01300O52
01300052
01300052
013000S2
01300052
OI30005?
O13000S?
O1300O52
01300052
01400022
0140002?
014OOO2?
014OOO??
0140O022
0140002?
RO12
1014
R016
EMISS
ROD?
POO*
RO06
ROO8
R010
R01?
001 4
9016
EMISS
RO02
ROO4
ROO6
ROOfl
R010
ROI2
ROI4
R016
EMISS
P002
0004
ROO6
ROO8
ROIO
R012
R014
R016
EM1SS
RO02
R004
RO06
POOS
ROIO
R012
RO14
ROI6
EMISS
RO02
R004
R006
BOOS
ROIO
RO12
RO14
R016
EMISS
RO02
RO04
RO06
ROO8
ROIO
1.6BSOO
I.2S1OO
1.26900
1 .1 I4OO
. 3890O
.16630
.18000
.34090.
.O895T
.2515O
.12B9O
. 15O1O
. I364O
.04763
.02036
.02204
.04173
.01097
.03080
.01578
.Ol 837
.34090
.11900
.OS09O
.05509
. IO43O
.02741
.07698
.03944
.O4592
3. 1 15OO
.39190
.83200
.22430
.59650
. 1 2460
.5757O
.11930
.28690
1.615OO
.20320
.43140
. 1 1 630
.30930
.06460
.29850
.06185
.I487O
8.43800
1.70200
.66910
3.487OO
1.50800
4.36400
R013
R01S
R017
ROOI
R003
POOS
ROO7
P009
P01 1
R0»3
ROI5
R017
ROOt
R003
POOS
ROO7
R009
ROM
P013
R015
R017
ROOI
R003
ROO5
ROO7
POO9
ROM
ROI3
RO15
RO17
ROOI
R003
POOS
POO 7
ROOQ
ROI I
ROI3
R015
RO17
ROOI
POO3
R005
ROO7
ROO9
ROI 1
R013
POI5
R017
ROOI
POO 3
RO05
RO07
ROO9
not i
I.IS9OO
1.347OO
1.23400
.251 1O
.27540
.16810
.25470
.2713O
.25950
.08666
.1396O
.14110
.03O74
.03378
.07058
.03118
.03321
.0317T
.O1O6O
.OI7O9
.01778
.076*4
.08427
.0514*
.07793
.08301
.079*0
.026*9
.04271
.04319
.244SO
.61610
I.288OO
.2084O
I.61100
.29620
.1162O
.17160
.42680
.12680
.31940
.66800
.1O80O
.835*0
.15360
.O6O24
.08899
.22130
?.382OO
.9265O
.729OO
I.614OO
I.21100
2.434OO
Figure 3 (continued). Test city example output
83
-------
f?XECUTOB.
MPS/360 V2-M11
RHS
01400022
01400O22
01*00022
01 700022
OIT00022
017OOO22
017O0022
01700022
017O0022
01700022
01700022
017OOO?2
O18OOO22
0180O022
0180O022
01800022
OIB00022
01800O22
01800022
01600022
01800022
S02
S02
SO2
SO?
SO 2
SO 2
SO2
SO 2
SO 2
R012
0014
R016
EM1SS
R002
R004
ROO6
RO08
R010
R012
R014
ROI6
EMISS
ROO2
R004
R006
R008
R01O
R012
R014
R016
ROOI
RO03
R005
R007
R009
R01 I
R013
R015
0017
1.92700
2.517OO
2.2S9OO
- 9.265OO
1.54400
.50960
2.61600
.73380
3.8280O
1.33000
6. 15900
2.67700
9.265OO
1.93000
.86720
2.35300
1.0170O
2.4760O
1.65500
4.79400
3.36100
152.9000O
151.20000
I49.400OO
153.300OO
146.50000
151.30000
154.60000
153.90OOO
151.40000
ROI3
R015
POI7
ROOI
R003
ROO5
ROOT
R009
RO1 1
R013
ROI5
R017
ROOI
R003
R005
R007
R009
R011
ROI3
R015
R017
R002
ROO4
R006
R008
R010
R012
RO14
ROI6
1 .56100
2.349OO
I .6B800
1 .93200
.86540
.49660
I.3020O
.79860
1 .6530O
5.59800
3.3600O
2.11 700
1 .824OO
1 .544OO
.51420
1 .5600O
.7356O
1 .6310O
1 .77800
6.16000
2.67TOO
152.400OO
149.90OOO
153.70000
149.600OO
153.90OOO
1 44.3OOOO
154.30000
153.30000
BOUNOS
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
EMISL
EMISL
EMISL
EMISL
FMISL
EMISL
EM JSL
EMISL
EMISL
EMISL
EMISL
EXISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMTSL
EMISL
EMISL
PMISL
OO100071
0010O071
00300071
00300071
O04O0061
OO400O61
00500021
005O0021
OO60OO21
00600021
OO70OO61
00700061
007O007I
OO700O71
008OO021
00000021
00900021
00900021
00900051
009OOOS1
01000021
01000021
160.60000
213.60OOO
1044.OOOOO
2088.00000
129.700OO
194.SOOOO
.00737
.00958
.00657
.00854
.01821
.O2367
.O43O4
.OS595
.00529
. 0068R
.00520
.00677
.01745
.02269
.0151 1
.01964
Figure 3 (continued). Test city example output
84
-------
FXECUTOH.
V?-M1I
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
CO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
LO
UP
FMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EM1SL
FMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
FMISL
EMISL
EMISL
EMISL
EMISL
EMISL
EMISL
01100031
01100031
01100041
01100041
011O0051
01100051
O120O061
012O006I
01200071
012O0071
013OOO21
013OOO21
01300051
01300051
01400021
01400021
01500031
O1500O31
01600021
01600021
OI60003I
0160OO31
01700021
01700021
OIBOOO21
OIS00021
0190O03I
01900031
02000031
02000031
02100021,
02100021
00500022
OO50O022
00600022
006OOO22
00700062
OO700062
00700O72
00700072
OO9 00022
0090002?
OO9OO052
OO9OOO52
01000022
01OO0022
01100032
0110003?
01IO0042
O110O042
01100052
OI1000S2
01300082
0130002?
.02374
.03086
.09693
.12600
.O77S4
. IOO8O
.02357
.03064
.02353
.03O58
.02063
.02682
. 05967
.07757
.01511
.01964
.03643
.04 736
.02305
.02997
.02429
.03158
.01 669
.O2I70
.Ol»77
.01 920
.02477
.03280
.03670
.047RI
.02555
.03322
.00657
.01970
.00657
.01970
. 0 1 82 1
.02604
.02663
.06155
.00657
.01970
.01152
.02836
.00657
.01970
.03643
.03940
.06494
.15750
.05195
.12600
.00657
.01970
Figure 3 (continued). Test city example output
85
-------
EXECUTOR.
MPS/360 V2-N1I
LO EMISL
UP EMISL
LO EMISL
UP EMISL
LO EMISL
UP EMISL
LO EMISL
UP EMISL
ENOATA
OI30OO52
0130OOS2
OI4OO022
0140O022
01700022
017OOO22
01000022
OI80OO22
.03998
. 09696
.OO657
.01970
.O0657
.0197O
.00657
.01970
EXECUTOR.
SETUP PBPILE
TIME ° 0.28
BOUND = EMISL
MPS/360 V2-M11
MATRIX1 ASSIGNED TO MATRIX!
ETA1 ASSIGNED TO ETA1
SCRATCH1 ASSIGNED TO SCRATCHl
SCRATCH2 ASSIGNED TO SCRATCH2
MAXIMUM PRICING NOT REQUIRED - MAXIMUM POSSIBLE
NO CYCLING
POOLS
H.REG-BITS MAP
WORK REGIONS
MATRIX BUFFERS
ETA BUFFERS
ROWS (LOG.VAR.t
COLUMNS (STR.VAR.)
774 ELEMENTS - DENSITY = 71.66 -
BE1
9
3
5
R SIZE
168
3504
7152
CORE
92
1512
IOS12
35760
TOTAL NORMAL .FREE.
18
42
17 1
0 0
FIXED BOUNDED
0 0
0 42
3 MATRIX RECORDS (WITHOUT RHS'S)
PRIMAL OBJ = EM1SS
TIME = 0.3O MINS.
SCALE = .
SCALE RESET TO 1.00000
OPTIMAL SOLUTION
RHS = SO2
PRICING 7
FXECUTOR. MPS/360 V2-M11
SOLUTION (OPTIMAL)
TIME = O.31 MINS. ITERATION NUMBER = 0
...NAME...
FUNCTIONAL
RESTRAINTS
BOUNDS....
...ACTIVITY...
242.15117-
DEFINED AS
EMISS
SO?
EMISL
Figure 3 (continued). Test city example output
86
-------
EXECUTOR.
SECTION I - BOWS
MPS/360 V2-M11
NUMBER ...ROW.. AT
.ACTIVITY...
SLACK ACTIVITY
1
2
3
4
S
6
7
8
<»
to
1 1
12
1.1
14
15
16
17
IS
ewiss
ROOI
RO02
RO03
RO04
POOS
R006
9007
ROO8
ROO9
R010
R011
R012
R013
R014
R01S
R016
ROI7
BS
BS
BS
BS
BS
BS
BS
BS
BS
BS
RS
BS
BS
BS
BS
BS
SS
BS
242.15117-
34.81524
97.17264
112.48956
85.86729
74.34O54
18.32188
19.60387
53.47604
55.87774
16.92212
32.34656
31.85386
15.89226
20.3751 1
2I.5302S
81.17437
21.53262
242.15117
118.08476
55.22736
38.7IO44
64.O3271
75.05946
135.37812
133.69613
96.12396
9O. 62226
136.97788
118.95344
112.44114
138.70774
133.92489
132.36975
I3P. 12563
129.86738
.LOWER LIMIT.
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
NONE
.UPPER LIMIT. .DUAL ACTIVITY
152.
152,
151
149
140,
153
153,
149,
146
153,
151,
144,
154
154,
193
183
151
NONE
,«OOOO
,4OOOO
,2OOOO
,9OOOO
,40000
,TOOOO
,30000
,60000
.SOOOO
,90000
,30000
,30000
.6OOOO
,30000
,90000
.30OOO
,40000
i.ooooo
Figure 3 (continued). Test city example output
87
-------
M°S/36O V2-M11
SECTION 2 - COLUMNS
r*\ji* r>C rf
19
?0
?1
22
?3
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
4fl
49
SO
51
52
S3
54
55
56
57
SB
59
60
. v-ni_ vimn .
O010007I
O0300071
OO400O61
O0500O21
OO600021
00700O61
00700O71
OOBOO021
O090OO21
OO90OO51
OIOOO021
01100031
Ol IOOO41
O110OO51
OI20O061
OT2OOO71
O13OOO21
013OO051
OI4QOO2I
O15OOO31
O16OOO21
OI6OOO3I
O17OOO21
01800021
0 190OO31
OPOOOO31
021OO021
00500022
00600022
00700062
OO7OOO72
OO9OOO22
O090O052
010OOO22
OI10O032
O110OO42
O I1OO052
013OO022
O1300O52
OI4OOO22
01700022
01BOO022
n i •
UL
UL
in.
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
Ut-
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
213.60000
2O88.00OOO
194.50000
.O09S8
.OOR54
.02367
.05595
.00688
.00677
.02269
.01964
. 03O86
.12600
.10080
.03064
. O3OS8
.02682
.07757
.01964
.O4736
.02997
.03158
.O2170
.01920
.03220
.O4781
.03322
.01970
.01970
.02604
.06155
.0197O
.02836
.01970
.03940
.15750
.126OO
.0197O
.09696
.01970
.OI970
.01970
INPUT COST..
.80940-
.O2O23-
.08094-
15.OOOOO-
15.54000-
1 1.670OO-
4.93800-
22.5000O-
15.27000-
.BO360-
14.06000-
44.55OOO-
5.45500-
13.64000-
7.5OOOO-
I 1.25OOO-
10.38000-
5.385OO-
14.06000-
17.5OOOO-
7.37500-
1O.50000-
13.240OO-
13.24OOO-
17.50OOO-
12.25OOO-
9.0000O-
7.SO 000-
6.21000-
2.33300-
.98750-
6.1O700-
i32140-
8.43800-
1.11400-
.13640-
.3409O-
3.1150O-
1.61500-
8.43800-
9.2650O-
9.26500-
O LIHIT.
I6O.60OOO
1044.00000
129.70000
.00737
.00657
.01821
.04304
.00529
.OO52O
.01745
.01511
.02374
.09693
.07754
.O2357
.O2353
.02063
.05967
.01511
.03643
.023O5
.O2429
.01669
.01477
.•O2477
.03678
.O2555
.00657
.O0657
.01821
.02663
.00657
.01152
.00657
.03643
.O6494
.OS19S
•OO657
.03998
.00657
.00657
. 00657
.UPPER LIMIT.
213.60000
2088.00000
194.5000O
.00958
.00854
.02367
.05595
.O06B8
.00677
.02269
.01964
.03086
.12600
.IO08O
.O3064
.03 058
.O2682
.07757
.01964
.04736
.02997
.03158
.0217O
.01920
.O322O
.04781
.O3322
.O197O
.01970
.02604
.06155
.O197O
.02836
.01970
.03940
. 157SO
.1260O
.01970
.09696
.01970
.O197O
.01970
COST.
.80940-
.02093-
.08094-
IS.OOOOO-
15.S400O-
11.67000-
4.938OO-
22.SOOOO-
15.27OOO-
.8036O-
14.060OO-
44.55000-
5.4S5OO-
13.64OOO-
7.50OOO-
I 1.25000-
I0.38OOO-
5.38500-
14.06000-
17.5OOOO-
7.375OO-
10.50000-
13.2400 O-
13.24OOO-
17.500OO-
12.25OOO-
9 . OO OO 0—
7.500OO-
6.2IOOO-
2.33300-
.98750-
6.1O7OO-
.32140-
A.4380O-
1 .114OO-
.13640-
.3*090-
3. I1SOO-
1.6150O-
8.438OO-
9.26500-
9.26SO0-
Figure 3 (continued). Test city example output
88
-------
RANGE
TIMfi =
EXECUTOP.
MPS/360 V2-MI1
0.32 KINS. ITERATION NUMBER
...NAME...
FUNCTIONAL
RESTRAINTS
ROUNDS....
...ACTIVITY...
242.15117-
OEFINED AS
EMISS
SD2
EM ISC
Figure 3 (continued). Test city example output
89
-------
VO
O
EXECUTOR. MPS/36O V2-M11
SECTION 2 - COLUMNS AT LIMIT LEVEL
.COLUMN. AT ...
19 OO10O071 UL
20 OO3O0071 UL
21 OO4O0061 UL
22 OO5O0021 UL
23 OO6O0021 UL
24 OO70O061 UL
25 OO7OOO71 UL
26 O08OOO21 UL
27 OO9OOO2I UL
28 OO90O051 UL
29 010OOO21 UL
30 01100031 UL
31 01100041 UL
32 0110OO5I UL
33 01200061 UL
34 O12OOO71 UL
PAGE 20 -
213.59991
2O87. 99921
194.49992
.00958
.00854
.02367
.05595
.O0688
.OO677
.02269.
.01964
.O308A
.126OO
.10080
.O3O64
.030-58
.80940-
.02023-
.08O94-
1S.OOOOO-
15.54000-
I 1.6700O-
4.9380O—
22.SOOOO-
15.2700O-
.80360-
I4.O600O-
44.55003-
5.45500-
13.64000-
7.5OOOO-
11.25001-
• • U.*J"C. ^ l_ 1 >" * i •
..UPPER LIMIT.
160.59995
213.59993
1O43.9996O
2087.99954
129.69996
194.49994
.OO737
.00958
.OO657
.OO854
.O1821
.02367
.O43O4
.O5595
.00529
.OO688
.00520
.OO677
.01745
.02269
.0151 1
.01964
.O2374
.03O86
. O9693
.12600
.07754
.10080
.02357
.O3064
.02393
.O305*
l_LJ Ht T «\_ 1 I V I 1 I
UPPER ACTIVITY
INFINITY-
299.41330
INFIN1TY-
8235.37408
INFINITY-
1294.54006
INFINITV-
6.O7985
INFINITY-
4.47856
INF1NITY-
7.65435
INFINITY-
IS. 08596
INFINITY-
5.534O8
INFINITY-
9.87522
1NF1N1TY-
187.53O70
INFINITY-
13.36925
1NF1NITY-
3.54679
INF1NITY-
28.84294
INFINITY-
1 1.58417
INFINITY-
21.57657
INFINITY-
14.39454
* • * \J ri i f t.w -» • • •
...UNIT COST..
.80940
.8O94O-
,O?O23
.O2O23-
.08094
,03094-
15.OOOOO
IS.OOOOO-
I5.54OOO
I5.54OOO-
1 1.67OOO
. 1 1.67OOO-
4.93800
4. 9380 O-
22.500OO
22.50000-
15.27000
15.270OO-
.80360
.80360-
14.O6OOO
14.O6OOO-
44.S50O3
44.55OO3-
5. 4550 0
5.45500-
13.640OO
13.64OOO-
7.500OO
7.5OOOO-
11. 2500 I
11.25001-
• •ur-r-t-r* \JLJ .a I ••
..LOWER COST..
1NFIN1TY-
INFINITY-
1NF1NITY—
INFINITY-
*
INFINITY-
INFINITY-
*
1NFINITY-
INFINITV-
INFIN1TY-
INFINITY-
INFINITV-
INFINITY-
INFINITY-
INF1NITV-
INFIN1TY-
INF1NITY-
i_i FI i it nu
PROCESS.
NONE
RO03
NONE
RODS
NONE
ROO3
NONE
R002
NONE
ROO3
NONE
R003
NONE
R003
NONE
RO05
NONE
R005
NONE
ROO5
NONE
R002
NONE
R003
NONE
R003
NONE
ROO3
NONE
R009
NONE
f* 1
AT
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
Figure 3 (continued). Test city example output
-------
EXECUTOR. HPSX36O V2-HII
PAGE 21 - 77/125
VO
35 O13O0021 UL
36 O130O051 UL
37 01400021 Ul.
38 OI5OOO31 UL
39 016O0021 UL
40 O160O031 UL
01700021 Ut.
42 OISOOO21 UL
43 01900031 UL
44 O2OOO031 UL
45 O210O02I UL
46 OO500022 UL
47 OO60OO22 UL
48 OO70O062 UL
49 OO7O0072 UL
SO O 090 0022 UL
51 OO9OOO52 UL
026B2
07757
O1964
O4736
O2997
031S8
O2I7O
01920
O322O
04781
O3322
O197O
O197O
02604
06155
0197O
O2836
• • r*r~\j i ^_u ^> • • » » * h
..l
10.380OO-
5.38500-
14.O6OOO-
1 7. SO 000-
7.375OO-
10.50001-
13.24OOO-
13.24000-
17.SOOOO-
12.25000-
9.OOOOO-
7.5OOOO-
6.21 OOO-
2.3330O-
.9875O-
6.1O700-
.32140-
,U Wt." k. •""*•» *
JPPER LIMIT. l
.02O63
.O2682
.05967
.O7757
.01511
.01964
.O3643
.O4736
.02305
.02997
.O2429
.03158
.O1669
.021 7O
.O1477
.O1920
.02477
.03220
.04 781
.02555
.03322
.OO657
.01970
.00657
.01970
.01821
.O26O4
.O2663
.06155
.00657
.01970
.01 132
.O283*
JPPER ACTIVITY ..
INFINITY-
16.89619
INFINITY-
32.61697
INFINITY-
IB. 85338
INFINITY-
ID. 36 464
INFINITY-
32.10157
INFINITY-
22.56383
INFINITY-
IS. 24216
INFINITY-
IS. 06121
INFINITY-
10.66693
INFINITY-
IS. 84 798
INFINITY-
3 O. 75579
INF1N1TY-
1 2.I6O24
INFINITY-
1 1.20445
INFINITY-
38.202OO
1NFINITY-
90.21 162
INFINITY-
24.68596
IMFINITT-
468.SS680
..UNIT COST.. .,
10.38OOO
1O.38OOO-
5.38SOO
5.385OO-
14.0600O
14.06000-
17.5OOOO
17.50OOO-
7.37500
7.37500-
IO.SO001
10.50001-
I3.24OOO
13.24000-
13.24OOO
13.240OO-
17.50000
I7.5OOOO-
12.2SOOO
12.25000-
9.OOOOO
9.OOOOO-
7.50OOO
7.SOOOO-
6.21OOO
6.21000-
2.33300
2. 3330 O-
. 98 75 0
.9875O-
6.I07OO
6. 10700-
.32140
.32140-
.LOMER COST..
INFINITY-
•
INFINITY-
•
INFINITY-
INFINITY-
•
INFINITY-
•
INFINITY-
•
INFINITY-
•
INFINITY-
INFINITY-
•
INFINITY-
»
INFINITY-
•
INFINITY-
•
INFINITY-
•
INFINITV-
INFINITY-
INFIN1TY-
INFINITV-
PBOCESS.
NONE
ROO9
r*JME
R009
NONE
R010
NONE
RO03
NONE
R003
NOff
R003
NONE
ROI4
NONE
ROI5
M3NE
ROO3
NONE
ROO3
NONE
9003
NONE
R002
ROO3
NONE
ROO3
NO»«
ROO3
NONE
R005
NOME
no 05
AT
UL
UL
UL
UL
UL
OL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
Figure 3 (continued). Test city example output
-------
N)
EXECUTOR. MPSX36O V?-«ll
NUMBER .COLUMN. AT ...ACTIVITY.
52 OIOOO022 UL
53 O1IO0032 UL
S4 01100042 UL
55 O1IOOOS2 UL
56 O1300022 UL
57 OI3O0052 UL
S8 01400022 UL
59 01700022 UL
6O O18OOO2S UL
PAGE
22 - 77X125
VITY... ..INPUT COST.. ..LOWER LIMIT. LOWER ACTIVITY ..
..UPPER LIMIT. UPPER ACTIVITY ..
.01970 8.438OO-
.O394O 1.1140O-
.I575O . 1364O-
. I26OO .34O90-
.OI970 3.11500-
.09696 1.61500-
.OI97O 8.4380O-
.O197O 9.26500-
.01970 9.26500-
.00657
.0197O
.03643
.O394O
.O6494
.1S75O
.05195
. T260O
.O0657
.01970
.03998
.O9696
. OO6S7
.0197O
.OO6S7
.01970
.00657
.O197O
INFINITY-
22.40614
INFINITY-
14O.6OOI1
INFINITY-
1148.15377
INFIN1TY-
459.48789
INF1NITY-
56. 27 J 83
1NFINITY-
108.57464
INFINITV-
3.1.40781
INFINITY-
21.76427
INFINITY-
2 1 .SO 029
.UNIT COST..
.UNIT COST..
B. 43800
8. 4380 O-
1. 1140O
1.1I40O-
.1364O
.13640-
.34O9O
.34O9O-
3.11500
3.115OO-
1.615OO
1.6I5OO-
8.43800
8. 4380 O-
9.26500
9.265OO-
9.265OO
9. 26 SO 0-
.UPPER COST.. LIMITING AT
.LOWER COST.. PROCESS. AT
. NONE
INFINITY- ROO2 UL
. NONE
INFINITY- ROO3 UL
NONE
INFINITY- ROO3 Ut
INFINITY- ROO3 UL
. NONE
INFINITY- ROO9 UL
NONE
INFINITY- ROO9 UL
NONE
INFINITY- ROIO UL
. NONE
INFINITY- DO 14 UL
NONE
INFINITY- HOIS UL
Figure 3 (continued). Test city example output
-------
CXIFCUTOR. NPS/36O V2-N11
SECTION 3 - ROWS AT INTERMEDIATE LEVEL
23 - 77/125
VO
LO
CR ...OOW.. AT ...ACTIVITY... SLACK ACTIVITY ..LOWE" LIMIT. LOWER ACTIVITY .
..UPPER LIMIT. UPPER ACTIVITY .
2 ROOI BS 34.81525
3 R002 BS 97.17264
4 RO03 BS 112.46056
5 ROD* BS 85.86728
6 RODS BS 74.34O55
7 R006 BS 18.32188
8 R007 BS 19.60387
9 RO08 BS 53.47604
10 ROO9 BS 55.87773
11 R010 BS 16.92212
12 ROll BS 3?. 34656
13 R012 BS 31.P5886
14 R013 BS 15.89226
15 RCM4 BS 20.37511
16 R015 BS 21.53O24
17 R016 BS 21.17438
118.08475 NONF
152.S9999
55.22736 NONE
152.39999
38.71043 NONE
151.2000O
64.O3271 NONE
149.89999
75.O5945 NONE
149.39999
13S. 37811 NONE
153.7OOOO
133.69612 NONE
153.29999
96.12395 NONE
149.59999
9O. 62225 NONE
146.4999S
136.97787 NONE
153.89999
118.95343 NONE
151.29999
112.44113 NONE
144.29999
138.70773 NONE
154.59999
133.92488 NONE
154.29999
132.36975 NONE
153.89999
132.12561 NONE
153. 20999
34.76538
34.81525
97.11289
97.17264
1 1 2.44411
112.48956
82.67134
85.86728
61.59331
74.34055
18.28894
18.32188
19.6O3I1
19.6O387
53.45291
53.476O4
55.84799
55.87773
16.88917
16.92212
32.29O71
32.34656
31.75159
31.85886
15.81873
15.89226
2O. 29421
?O. 37511
2 1.44934
21.53O24
21.08789
21.1 7438
..UNIT COST.. ..UPPER COST.. LIMITING
..UNIT CO.ST.. ..LOWER COST.. PROCESS.
1.97524
INFINITY
1.64871
INFINITY
1.79428
INFINITY
1.64112
INFINITY
1.65684
INFINITY
1.93318
INFINITY
4.37377
INFINITY
2.29216
INFINITY
1.78317
INFINITY
1.93318
INFINITY
2.41951
INFINITY
1.78317
INFINITY
1.65 SOS
INFINITY
1.50430
INFINITY
1.50406
INFINITY
1.50344
INFINITY
00500022
NONE
OOSOOO22
NONE
OO6OOO22
NONE
OO400O61
NONE
OO3OOO71
NONE
OIOOOO21
NONE
Ol IOOO32
NO ft
O120OO61
NONE
012OOO61
NONE
O14OOO21
NONE
O2 00 OO 31
NONE
01 5O OO 31
NONE
O17OOO22
NONE
01700022
NONE
01800022
NONE
O19OOO31
NONE
AT
AT
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
Figure 3 (continued). Test city example output
-------
EXECUTOR. MPS/36O V2-M11 PAGE ?» - 77/125
NUN3ER ...»0».. AT ...ACTIVITY... SLACK ACTIVITY ..LOWER LIMIT. LOWER ACTIVITY ...UNIT COST.. . .UOPER COST.. LIMITING AT
..UPPER LIMIT. UPPER ACTIVITY ...UNIT COST.. ..LOVER COST.. PROCESS. AT
18 R017 OS 21.5326? 129.86737 HONE 21.4*25* l.»999« O2OOOO31 UL
151.39999 21.53262 INFINITY NONE
Figure 3 (continued). Test city example output
-------
SECTION VII
AN EXAMPLE OF THE EMISSIONS MAXIMIZATION
FORMULATION IN THE BALTIMORE AQMA
The emission density limit setting model was executed in two AQMAs,
Baltimore, Maryland and Louisville, Kentucky-Indiana. The emissions
maximization formulation was executed for both pollutants in the Baltimore
AQMA and the Louisville AQMA. The Baltimore example is discussed in this
section; the Louisville examples are discussed in Section VIII.
DATA COMPILATION
The emission inventory and land use inventory available for the
Baltimore AQMA were not compatible; in addition, neither geographic system
for which the data were compiled was consistent with the component area
criteria of Section II. Consequently, a new grid system was prepared.
The emission inventory was reallocated to the new grid system; the land
use inventory and plan were compiled by grid square using acetate overlays.
During the latter stages of the preparation of this report it was dis-
covered that an error was made in the partitioning of industrial area
source emissions in Baltimore into the light and heavy industrial land
use classes. For each grid square the inventoried industrial emissions
were partitioned into light and heavy industrial emissions by the light
and heavy industrial land use, weighted by the mean emission densities
of light and heavy industrial point sources. Through a programming error
the mean plant emission densities were calculated for general commercial
and light industrial rather than for light industrial and heavy industrial
sources. Since the emission densities for the general commercial plants
were much lower than for the light industrial plants, the heavy industrial
95
-------
Point source lot areas were obtained by consulting real estate tax assess-
ment records; estimates of floor area ratios were developed from discus-
sions with the staff of regional and local planning agencies.
After the required data had been developed the preparation of the
input proceeded as described in Section VI. A total of 2256 emission
sources in 524 component areas were identified. There were seven land
use categories, consistent with those specified in the regional plan
and local zoning ordinances, viz:
• low density residential
• medium density residential
e high density residential
• local commercial
• general commercial
• light industrial
• heavy industrial
There was an average of 4.3 emission sources, area and point, in each
component area.
All power plants and incinerators were classified an non-EDZ sources,
as well as all mobile sources and the miscellaneous emission sources on
NEDS area source card 5. All point sources in the emission inventory
that were larger than 25 tons per year were treated as point sources in
the EDZ model; those less than 25 tons were assigned to the appropriate
area source. A 10-meter stack height was assumed for all nonindustrial
sources. Based on an examination of the point source inventory, a
65-meter height was assumed for all industrial area sources.
area sources were assigned a disproportionally high percentage of the to-
tal industrial area source emissions, leaving effectively a single indus-
trial land use, namely, heavy industrial. For the Louisville simulation
this error was corrected. This error does not invalidate the usefulness
of the test; however, it is less representative of Baltimore and more of
a similar city with a light industrial base that has a very low mean emis-
sion density.
96
-------
To facilitate the manipulation of the emission and land use inven-
tories by computer, an eight-digit identifier was constructed for each
source. It consisted of the following:
digits 1-3 the component area index number
4-6 point source identification number,
(zero if an area source)
7 land use category (values of 2 through 8 and
corresponding to the seven land uses
noted above)
8 current-future emission source identi-
fier (1 = current, 2 = future growth)
18
The input to the MPS execution step consists of four discrete sets of data:
the ROWS section, the COLUMNS section, the RHS section, and the BOUNDS
section.
The ROWS section is shown in Figure 4. The rows section is pre-
ceded by a name card with a user specified problem name. Shown is EDZBAL
(emission density zoning - Baltimore). Following the name card is the rows
card indicating the beginning of the ROWS section. Following that is a
record for each row indicating the type of row and the user specified name
of the row. The first row is the objective function, called EMISS (emis-
sions). The N indicates it is the objective function. Subsequent records
are included for each receptor. The L indicates that the constraint is a
"less than or equal to" type. The name of each row is the receptor number.
The COLUMNS section, shown in Figure 5, contains, for each source,
a vector consisting of first, the area of the source (the objective func-
tion parameter) and second, for each receptor, the product of the area of
the source and its transfer coefficient to that receptor. Thus the first
source shown in Figure 5 is the current low density residential land
97
-------
NAME EDZBAL
ROWS
N EMISS
L RQOl
L R2
L R523
L R524
Figure 4. ROWS section of SOX emissions
maximization program
98
-------
COLUMNS
00100021 EMISS
00100021 R001
00100021 R002
-0.5440E 00
0.3076E 00
0.2950E 00
00100021 R524
00100032 EMISS
00100032 R001
0.1690E-01
-0.2400E 00
0.1357E 00
02101761 R524
03501861 EMISS
03601861 R001
0.3730E-02
-0.1919E 02
0.5622E 00
Figure 5. COLUMNS section of SOX emissions
maximization program
99
-------
use in component area 1. It has an area of 0.544 square kilometers.
3
With a transfer coefficient of 0.5654 yg/m /ton/day, the value of the
element in the matrix is 0.3076 (0.5654 x 0.5440). Since there are 524
receptors and 2256 emission sources, there are 1,184,400 ((524 + 1) *
2256) records in the COLUMNS section.
The RHS section, shown in Figure 6, contains the right-hand side
constraint value (i.e., upper limit on air quality) for each row (i.e.,
receptor). The vector is given a user specified name, in this case, S02.
The BOUNDS section, shown in Figure 7, contains the lower and
upper bounds on emission density for each source; this vector is also
given a user specified name, in this case EMISL. For example, the first
source shown, source No. 12200181, has a lower limit of 0.122 and an
upper limit of 2.249. The BOUNDS section contains two records for each
source, or 4512 records. It is followed by an ENDATA card.
The SO -PM feasibility constraint was found to be unnecessary based
on the analysis in Appendix C. If it had been required, it would have
been added to the columns section. The proximity constraint was also not
employed to simplify the testing. If used, it too would have been added
to the columns section.
EXECUTION OF MPS
The problem, as compiled above, contained more matrix elements than
had been expected. We had estimated that there would be fewer emission
sources per component area. To reduce the problem to a more manageable
size, the number of rows (i.e., receptors) was reduced. This was done
*
The linear program actually run was a minimization. Hence the minus sign
for the area, which in effect led to the minimization of the negative of
emissions, equivalent to maximizing emissions.
100
-------
RHS
S02 R001 0.1215E+3
S02 R002 0.1214E+3
S02 R003 0.1213E+3
S02 R050 0.1180E+3
Figure 6. RHS section of the SOX emissions
maximization program
101
-------
BOUNDS
LO EMISL
UP EMISL
LO EMISL
UP EMISL
12200181 0.1220E+00
12200181 0.2249E401
03500361 0.7936E-02
03500361 0.1090E+00
UP EMISL
ENDATA
49400082 0.3048E+01
Figure 7. BOUNDS section of SOX emissions
maximization program
102
-------
by using an option in MPS to use a subset of the constraint equations.
It is not recommended that this option be normally employed, for it is
relatively inefficient. The number of constraint rows; i.e., receptors,
was reduced to a set of 52.
Because it was unnecessary to include the SOX-PM feasibility con-
straints, the maximization of particulate matter and sulfur oxide emis-
sions could be solved independently. Therefore two separate linear
problems could be executed, one for each pollutant.
The control program for MPS used in this example is shown in Fig-
ure 8. The CONVERT command reads the input data, converts it to
binary, and stores it on a temporary disk file. In this example, it is
also using a masking option to use only those rows that end in the digit
five. SETUP takes certain sections of the problem and places them in a
working file. CRASH is a procedure to obtain a starting point that con-
tains fewer infeasibilities. It is an attempt at solving the linear
program more efficiently. PRIMAL is the command that solves the linear
program. SOLUTION outputs the optimal solution (see the example output
in Figures 9 and 10) and RANGE is a postoptimal procedure that calculates
how sensitive the final solution is to changes in the right-hand side of
the air quality constraint and to changes in the land area of each source.
OPTIMAL SOLUTION TO SULFUR OXIDES MAXIMIZATION
After presenting information on the characteristics of the problem
and an iteration log, the output of MPS presents the optimal solution to
the problem and the sensitivity of the solution to changes in the param-
eters. The solution is described in two sections. The first is the rows
*
Through the use of the masking option on the CONVERT command.
+The number of rows was thus 53, 52 constraints and the objective function.
103
-------
PROGRAM
INITIALZ
MOVE(XDATA,'EDZBAL')
MOVE(XPBNAME,'PBFILE')
CONVERT('SUMMARY','RMASKS','R**5', 'EMISS',' ')
SETUP('BOUND'.EMISL1)
MO VE(XOBJ,'EMISS')
MOVE(XRHS,'S02')
CRASH
PRIMAL
SOLUTION
RANGE
EXIT
PEND
Figure 8. MPS control program
104
-------
EXECUTOR.
SECTION 1 - ROWS
MPS/36O V2-HII
77/01O
o
Ul
NUMBER
2
3
: 4
5.
6
7
8
9
10
11
12
13
14
15
16
17
! 20
! 21
22
23
24
26
26
i 27
i 28
29
30
31
32
— if-
35
36
37
38
2.9
40
42
; 43
. . .ROM. .
R005
RO1S
R02S
R035
R045
R0o5
R065
R075
RO85
R095
R105
R115
R125
R135
R145
R155
R16S
R175
R18S
R195
R205
R215
^225
R235
R245
R2S5
R265
R275
R285
R295
R305
R315
— R32S
R33S
R345
R355
~R36S
R375
R385
R395
R4O5
R415
R425
R435
R445
R45S
R475
AT
BS
BS
BS
. UL
UL
BS
BS
.-BS_
BS
BS
BS
BS
BS
UL
BS
UL
BS
BS
BS
UL
BS
BS
UL
BS
UL
BS
BS
bS
BS
BS
BS
BS
bS
UL
__BS_ _
UL
BS
US
BS
as
BS
BS
BS
BS
BS
BS
. . .ACT IV1TY. . .
0*3.66417-
26.46435
44.7399O
59.57467
119.90000
118.50000
101 .92672
1 14.33388
48.7740
-------
EXECUTOR•
SECTION 2 - CULUMNS
MHS/360 V2-M11
PAGE
5O -
77/010
NUMBER
54
55
56
57
• 58
59
60 .
61
62
63
64
65
66
67
66
69
70
71
72
73
74
75"
76
77
78
79
80
81
82
— 83""
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
9*
100
101
102
.COLUMN.
" '00 10*0031
00200021
O020O022
00200051
"00300021
00300O31
00300051
00400021
00400022
O04000S1
00400071
005OO021
00500022
00300051
O05O0071
00600021
00600051
00600071
O0600O82
O070O021
OO700O22
00700051
00700071
O070008i
00800021
00800022
00800051
OO800071
"6090O02 1 '.'
00900022
01 O00021
010OOO22
"01 SOOT/51
01 100021
01100022
01100031
Ol 10O051
01 100061
01 100072
O1200021
0120002?
01200031
01200052
O12OOO61
O120OO71
01200072
Ol JOOO21
O13OO022
A r . . . AC
„ .jJi
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
UL
.06_900
.12900
.06900
.069OO
.28200
.06900
.12900
.28200
.06900
.06900
.28200
.3O600
.06900
.06900
.2U200
.3O600
.06900
.28200
.30600
3.0480O
.06900
.26200
.30600
3.04800
" ~.0~6900 '"
.06900
.28200
.30600
— .069W
.06900
.06900
.06900
V2"8200 ~"
.06900
.06900
. 12*OO
.28200
.50dOO
.3OoOO
.06900
.06*00
.12900
.28200
.5O600
.30*600 " -"'
.30600
.06*00
.06*00
:>UT COST.. ..LOWER LIMIT. ..UPPE* LIMIT. .REDUCED COST.
.5*400- .0145? .06900 .=-2o59-
.24000-
2. 1760O-
3.80800-
.72OOO-
.59840-
.24000-
.2-000-
1.08600-
13.60000-
. 4800O- -
.09600-
1.36000-
2.99200-
.24000-
.24000-
1.25100-
.O96OO-
.09600-
.04800-
2.720OO-
2.17oO6-
.4*000-
.04600-
.04600-
1.36000-
2.72000-
. 24OOO-
.0360O-
1.3600O-
2.99200-
2.1760O-
4.216OO-
.T20T)0-
1.63200-
3.80800-
.24OOO-
".24000-
2.496OO-
.46COO-
1.63200-
I 14400-
.48000-
.41600-
' .04600-
,43200-
.81600-
8.7040O-
.03458
.00767
.02300
.00361
.O1387
.03375
.Ol 167
.01645
.02300
.00604
.00625
.02300
.02300
.02667
.00625
.01463
.04200
.01354
.09500
.00728
."52300
.01208
.O2600
.09500
.01669
.02300
.01O42
.01667
".0*300
.02300
.01466
.O2300
.03633
.00619
.02300
.01500
.00042
.00168
.02600
.01060
V02"300
.02569
.04200
.01442
" V0~2"600
.02600
.02300
.02300
. 1?900
.O69OO
.069OO
.282OO
.O69OO
. 12900
.28200
. 06900
.06900
.28200
.30600
.06900
'""- V06900 "
.28200
.30600
. O6900
" .28200 "" "
.30600
3.04800
.0690O
. 069OO
.2820O
.30600
3.04800
.0690O
.06900
.28200
.30600*
i'06$"00
.06900
.06900
.06900
.28200
. 069OO
. 069OO
.12900
I 50800
.oObOO
.06900
I 12900
.28200
.3O80O
.30600
.30600
.06900
.06900
.1O820-
1 .6*739-
2.97054-
.56164-
.45371-
.18197-
.18197-
.8O390-
10.04a62-
.35465-
.07099-
1.00678-
2.21484-
.17766-
.17780-
.89957-
.O6903-
.O69I4-
.03457-
1.86742-
I033O1-
.03301-
.93136- "
1.86278-
.16436-
.O247O-
2ll0994-
1 .42907-
2.76CJ95-
.07882-
1 .04305-
2.43357-
.15338-
.13338-
1.59520-
.30776-
1.O4886-
.09255-
.30851-
.26737-
.03093-
.27961-
.490O9-
5.31220-
Figure 10.
COLUMNS section of SO optimal solution
X
-------
sections, shown in part in Figure 9, which lists for each row an
internal variable number, the row name, its status (BS indicates the row
is in the basis and not binding, UL indicates the row is binding and at
its upper limit), the row activity (the ambient air quality at each re-
ceptor), the slack activity, the lower and upper limits as assigned by
the user, and the dual activity. Of most interest is whether a row is
at the upper limit, and the level of slack activity; i.e., unused tol-
erance for SO or the difference between the NAAQS and the projected con-
X
centration at the receptor. As can be seen in Figure 9, nine receptors
are at their upper limit. The lower limit for all receptors is zero;
i.e., ambient air quality is constrained to be nonnegative.
The second section of the solution output is the columns section, a
portion of which is shown in Figure 10. For each column (source), the
internal variable number, the column name (source index number), its
status, the column activity (assigned emission density limit in units of
2 2
tons/km /day), its objective function parameter (area in units of km ),
The dual activity is the value of the dual variable in the dual problem.
The dual can be interpreted as minimizing the marginal value of the air
quality standards used in the primal right-hand side, that is, minimizing
the rate at which total emissions may be increased by increasing the
standards. The value of the dual variables are sometimes called shadow
prices, and can be thought of as the maximum unit price one should be willing
to pay to increase the allocation of that resource; i.e., the air quality
standard. The shadow prices provide a convenient ranking of the importance
of each receptor. If one receptor is found to have a very high shadow
price, any non-EDZ sources (e.g., power plants) that have a significant
contribution to the ambient air quality at that receptor should be examined
for the possibility of additional emission controls, thereby increasing
the right-hand side value for the receptor. Also, if the receptor is a
sensitive receptor (thus with a lower total concentration than the NAAQS),
the desirability of maintaining its status should be reassessed.
107
-------
the lower and upper limits assigned by the user, and the reduced cost
of the associated variable. Of most interest is the activity column,
which indicates the emission density limit assigned to the source.
Interpretation of Solution
Tables of numbers such as those shown in Figures 9 and 10 are diffi-
cult to interpret, especially in a large problem such as the one at
hand. For this reason, the output of MPS was routed to a disk file.
The more interesting values were retrieved and displayed graphically using
27
SYMAP. The most interesting parameter in the rows section of the so-
lution output is the slack air quality, or the difference between the
standard and projected ambient concentrations at each receptor. These
values are plotted as a contour map shown in Figure 11. Immediately
apparent from the figure is that all the binding receptors are grouped
together in the center of the city or to the east over the Chesapeake Bay.
One would suspect, therefore, that emissions would be allocated by the
linear program to sources most distant (and therefore with smaller transfer
coefficients) from this grouping of binding receptors. This is obvious
from inspecting the spatial distribution of the emission sources that
are either at their upper or lower limit, or are unconstrained; i.e., at
neither limit. Figures 12 and 13 show, respectively, the status of each
grid square component area for current and future light industrial
sources. In both figures, those component areas on the periphery of the
region tend to be binding at the upper limit while those sources at the
*
The reduced cost is the rate of increase in the objective function, total
emissions, per unit increase in the emission density limit. At optimality,
an increase in a source currently at its upper limit will cause the total
emissions to increase, an increase in a source at its lower limit will
cause the total emissions to decrease, and an increase in a source in the
basis will cause no change in total emissions. See the section on sensi-
tivity analysis for more detail on these relationships.
108
-------
Figure 11. Slack SO levels (micrograms per cubic meter)
A
109
-------
LEGEND.
L AT LOWER LIMIT
* AT UPPER LIMIT
Figure 12. Status of bounds on emission density limits, current
light industry
110
-------
LEGEND.
L AT LOWER LIMIT
• AT UPPER LIMIT
Figure 13. Status of bounds on emission density limits, future
light industry
111
-------
center of the region tend to be binding at the lower limit. The same
pattern is evident in Figures 14 and 15 which show the same information
for heavy industrial sources.
Three critical insights to the emission density maximization for-
mulation as applied to sulfur oxides in Baltimore can be stated, viz:
• The critical (i.e. , binding) receptors tend to be in the
center of the region.
• The emission density limits for most sources will be set
at either the upper or lower limit.
• The emission density limits for future industrial sources,
which have identical bounds, will be higher on the periphery
of the region. There is then a possible sprawl inducing
potential inherent in simply maximizing emissions.
Sensitivity Analysis
The remainder of the output provided by MPS is a sensitivity analysis
of the optimal solution. The problem is, however, so constrained that
much of the analysis is of little interest. First, the maximum ambient
air quality at each of the binding nine receptors can be changed only
3
slightly (less than 0.3 yg/m ) before the optimal solution changes; i.e.,
an emission source not at an upper or lower limit in the current solution
reaches that limit.
The second part of the sensitivity analysis indicates how much the
assigned emission density limit at binding sources can change if the land
area of that source were changed. It also indicates how the value of the
objective function will change. In general, the model is fairly insensi-
tive to changes in the area of an emission source. The area of most
sources could change by a factor of two or more before the emission den-
sity limit on another source is affected. The emission density limit of
a source whose area is changed is greatly altered, of course.
112
-------
LEGENDi
L AT LOWER LIMIT
• AT UPPER LIMIT
Figure 14. Status of bounds on emission density limits, current
heavy industry
113
-------
LEGENDi
L AT LOWER LIMIT
• AT UPPER LIMIT
Figure 15. Status of bounds on emission density limits, future
heavy industry
114
-------
The third part of the sensitivity analysis indicates how the maximum
ambient air quality at unbinding receptors can change without altering the
3
optimal solution. This analysis shows that a slight change (0.1 yg/m )
in the ambient air quality of a nonbinding receptor will cause an emis-
sion source to reach an upper or lower limit.
The final section of the sensitivity analysis indicates how much the
area of an emission source not at an upper or lower limit may change
without altering the optimal solution. The emission density limit for
2
such sources can change only slightly; i.e., less than 0.2 tons/day/km ,
before another emission source currently at its upper or lower limit is
forced to change its emission rate.
PARTICULATE MATTER EXAMPLE
In a similar manner, a linear program was executed to maximize par-
ticulate matter emissions in the Baltimore AQMA. Baltimore currently is
recording violations of the annual particulate NAAQS. To allow the
linear program in this example to reach a feasible solution, attainment
was assumed. Total suspended particulate levels were assumed to be
20 yg/m3 less than they are in the base year.
The optimal solution to the problem was very similar to that for the
sulfur oxide example. However, in this case there was only one critical
receptor, located in Baltimore City. A contour plot of slack TSP level
is shown in Figure 16.
A plot of the status of the emission density limits is shown on
Figure 17 for current light industry. In this example most of the
sources are at their lower limits and only a few sources at the northern
edge of the region are at their upper limits. In contrast with the sulfur
In other words, a uniform rollback of 20 percent of the emissions at all
current sources was assumed.
115
-------
Figure 16. Slack TSP levels in Baltimore (micrograms per cubic
meter)
116
-------
LEGENDi
L AT LOWER LIMIT
• AT UPPER LIMIT
Figure 17. Status of bounds on emission limits, current light
industry
117
-------
oxide example, there was relatively less room for growth in emissions in
the particulate matter example. Almost all sources had their emission
limits set at the specified lower limit to allow for all the development
contemplated in the land use plan. The particulate matter example does,
however, support the observations made at the conclusion of the preceding
interpretation section. In this case:
• There is one critical receptor at the center of the region.
• Sources at the periphery of the city tend to have their
limits set at the specified upper bound and those nearer
the center of the city (and the binding receptor) are set
at the lower bound.
118
-------
SECTION VIII
AN EXAMPLE OF THE EMISSIONS MAXIMIZATION
FORMULATION IN THE LOUISVILLE AQMA
DATA COMPILATION
The emission inventory had been prepared on a census track basis.
Though it would have been possible to use census tracts as component
areas; to fully test the criteria of Section II it was decided to prepare
a new grid system and reallocate the emission inventory to this grid
system. The new grid system is shown in Figures 18 and 19.
Each county in the AQMA; i.e., Clark, Floyd, and Jefferson, had a
prepared land use plan. However, the plan was presented on a different
scale map in each county, respectively, 1 in. = 5000 ft, 1 in. = 500 ft,
and 1 in. = 6500 ft. Clark and Floyd counties also had land use inven-
tories compiled on colored base maps — of differing scales, 1 in. = 2000 ft,
and 1 in. = 500 ft. Jefferson County's land use inventory was compiled
by census tract. Jefferson County's land use inventory was reallocated
to the grid square system by measuring the percentage of each census
tract in a grid square. The remaining two land use inventories and the
land use plans were transferred to the grid system by preparing an acetate
overlay of the appropriate scale and measuring the amount of each land
use category in a particular grid. This process took 2 man-months.
Point source lot areas were obtained by consulting real estate tax
assessment records. This required 1 man-week.
119
-------
Figure 18. Louisville AQMA component area grid system
120
-------
Figure 19. Numbering of Louisville AQMA component areas
121
-------
To conserve resources, light and medium density residential land use
was treated as a non-EDZ source, their emission density limit was es-
tablished a priori and their projected contribution to each receptor was
subtracted from the right hand side. Their emission density limits were
set two geometric standard deviations larger than the geometric mean emis-
sion density of each class. The remaining land use was partioned into
three categories, viz.,
• Commercial, consisting of high rise residential,
commercial, office, and institutional land uses.
• Light industrial
• Heavy industrial
All power plants and incinerators were classified as non-EDZ sources,
as well as all mobile sources and the miscellaneous emission sources on
NEDS area source card 5. All point sources in the emission inventory
that were larger than 25 tons per year were treated as point sources;
those less than 25 tons were assigned to the appropriate area source. The
same stack heights used in Baltimore; i.e., 10 and 65 meters, were assumed.
Following these procedures, a total of 629 emission sources in 248
component areas were identified; thus, there were an average of 2.5 emis-
sion sources in each component area. This is 59 percent of the average
in Baltimore and can be attributed to the deletion of the two residential
classes and the aggregation of high rise residential, commercial, office
and institutional into one category.
AQDM was executed to determine the calibration coefficients. However,
the calibration was poor; consequently, the linear programs were executed
on uncalibrated data. It is thought that the poor calibration can be
attributed, at least in part, to questionable monitoring data.
The lower bound for existing sources was set at their current emis-
sion density. Point sources, whose emission density was expected to
122
-------
decline in the desing year (1990) due to current regulations, had their
lower limit set at this projected emission density (which does not include
capacity expansion). The upper bound for existing sources was set at
their projected 1990 emission density, based on the projected emission
inventory. Sufficiently detailed floor area ratios could not be readily
obtained for future nonindustrial sources. Consequently, their upper and
lower bounds were determined in the same manner as future industrial
sources. The upper and lower bounds for all future sources, therefore,
were set by estimating the smallest and largest emission density that
could be expected in each land use class. This was determined by analyzing
the distribution of existing emission densities in each class, which was
found to approximate a log-normal distribution. The mean was selected
as the lower bound and the sum of the mean and two standard deviations
was selected as the upper bound.
The input to the MPS execution step was identical in format to
Baltimore. The feasibility and proximity constraints were also not used.
EXECUTION OF MPS
MPS was executed three times, once for particulate matter and twice
for sulfur oxides. The sulfur oxide maximization was run once with
40 percent of the total receptors and once with a full compliment of
receptors. The control program was similar to that illustrated in the
Baltimore example.
PARTICULATE MATTER EXAMPLE
The linear program reached an infeasible solution in the particulate
matter example. That is, with all sources at their lower limits, the
NAAQS was still violated. Emissions at this point included the projected
1990 emissions from non-EDZ sources (power plants, mobile sources, etc.),
the existing emissions from current sources, no capacity expansion at
123
-------
existing sources, and complete development of the vacant land according
to the land use plan with a moderately stringent emission density limit.
Assuming the validity of the dispersion model and the emission inventory,
several options present themselves, viz.,
& Apply additional control strategies at non-EDZ sources.
o Require reductions in the emission density of existing
sources covered by EDZ.
• Rerun the EDZ model with a smaller lower limit for
future sources.
All three are viable strategies in the present case. Non-EDZ sources
generally account for half the particulate concentration at most receptors.
There are several large point sources included in the emission density
limit setting model with significant contributions to receptors projected
to violate the NAAQS. Finally, a more stringent lower limit on future
sources could be considered.
SULFUR OXIDES EXAMPLE
Optimal Solution - 40 Percent Sample
The sulfur oxide example of the emissions maximization formulations was
first run with a 40 percent random sample of total receptor set. There-
fore, 102 receptors were utilized. The optimal solution to the emissions
maximization had a total of 238.8 tons per day of sulfur oxide emissions.
Two receptors, R083 and R087, were at their upper limit; i.e., at the
NAAQS. Thirty-six sources were set at their upper limit, two at neither
bound, and the remainder at their lower limit.
Interpretation of Solution - 40 Percent Sample
The slack air quality is shown in Figure 20. As was the case in
Baltimore, the two binding receptors, along with four others that are
close to the standard are in the core of the metropolitan region. (This
124
-------
is the area with the symbol '.' in Figure 20. The opposite of slack
air quality is shown in Figure 21, which is the total projected sulfur
oxide concentrations. It is interesting to note the contribution to
ambient air quality from EDZ sources and non-EDZ sources. This is shown
in Figures 22 and 23, respectively. The values represented in Figure 23
are, of course, the amount that was subtracted from the NAAQS to obtain
the right-hand side.
Reviewing the pattern of emission density limits set at the upper
limit, the lower limit, or unconstrained, Figures 24, 25, 26, and 27, the
same pattern is discerned as was the case in Baltimore; sources at the
core of the region tend to have their emission density limit set at their
lower bound, those on the periphery tend to have theirs set at the upper
limit. In addition, note that for some component areas, a current source
is set at the upper limit while a future source of the same category is
set at the lower limit. Since the transfer coefficients for two such
area sources would be the same, the emission density limits could be
traded off between current and future sources and still remain at the
optimal solution — so long as total emissions from the component area for
that land use category remained unchanged.
The emission densities in the Louisville AQMA for the base year, the
design year, and the design year with all sources at their upper limit are
shown in Figures 28, 29, and 30.
Sensitivity Analysis - 40 Percent Solution
/
As discussed in the Baltimore example, the sensitivity analysis is
presented in four parts, viz., rows (receptors) at the limit level,
columns (sources) at the limit level, rows at an intermediate level,
and columns at an intermediate level.
The first part indicates how the optimal solution would change if the
value of right-hand side were increased or decreased. For example, at the
125
-------
Figure 20. Slack SO levels (micrograms per cubic meter)
X
126
-------
Figure 21. Projected SO concentrations (microgram per cubic meter)
X
127
-------
Figure 22. SOX levels attributed to EDZ sources
(micrograms per cubic meter)
128
-------
..r
Figure 23. SOX levels attributed to non-EDZ sources
(micrograms per cubic meter)
129
-------
JJJJJJJJJJJJJ
.U.LLLLILU.
LU.LLU.U.U.
LLLLLU.LLU.
LLLLLXLULLL
tLU.LLU.LLU.
LLLLLLU.LLLL
U.U.LLLLLLULL
LLLLLLLLLLU.L LLI
LLLU.LLLLLLLLL U.IU.L N. tXMM •
ULtLLLLLLLULLL ^J j
LLIU.H.LLLLLLL, LL^ttL *• *M LLL «MT
LLIULULLLLLLLL LLt I4.L •» «•» LLL «••
I.LI LLLLLLLLt-LL * . /
txl U.LLLLIL1.I.U. ^ • ~^*+^ 97
LEGEND*
L AT LOWER LIMIT
« AT UPPE«
Figure 24. Status of bounds on emission limits, current light industry
130
-------
- —
rf..
j psssss
V ^ralwrnmB
t ' • iwpnftDMifiarmn
• fmBsmmnmmn
J Ben«RBfwnanen LLI
I raeieeeeBaoea LLI
* ' BPPOBttBBMBHP LL
I " BafUlMMBtMiBQR LL
1 LLLLLLLLLLLLL LLI
1 " — r LLLLLLLLLLLLL IL1
* — ••!••!. — — y LLLLLLLLLLLLL LLI
t LLLLLLLLLLLLL LLI
\ LLLLLLLLLLtLL
1 LLLLLLLLLLLLL LLI
1% LLLLLLLLLLLLL LLI
& LLLLLLLLLLLLL LL
\ LLLLLLLLLLLLL LLI
\ LLLLLLLLLLLLL Lt
.^LLLLLLL
^~ LLLLLLL
v. f *^ LLLLLLL *
1 • I LLLLLLL LL LLL Uf
• - ^^ LL LLL LLL LLL/LL
I j LL LLL LIL LUt LL
> *J -v /
L :('
i J
1 i L
/•
^
LLLLLL
LLL LLLLLL
/ . LLL .LLLLLL
9 LL LLL LLLI LL
I—— ••"/ LLLLLL
I UU LLL * ' LLL LLLLLL
I UU LLL LLL LLLLLL
1 LLLLLLLLLLLLL LL
1 LLLLLLLLLLLLL LL
1 LLLLLLLLLLLLL LL
L
/
y
^^r
+^r
^ s~
ALL X
-LLL /
XLL/
ALL/
4.L< LLLLKLL
JLLL LLLLLLL
.LLL LLLLUL
.LUL LLLLUL
A.m-
.L|L
3 \
/ ^LLLLLLLL LLLLLL
' XLLLLLLLL LLLLLL
*tt^^LLL LLLLLL
LLLLCkLL LLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLMM.
LLLLLLLLLLLLX^
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL '
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLILLL LLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LLLLLLLLLLLLL
LL LLL LLLLLL
LL LLL LLI LIL
LLLLLLLLLLL
1 LLLLLLLLLLLLL LLLLLLLLLLLLL ^
i LLLLLLLLLLLLL LL
LLIL^rWUJL LL
»* LLLLLLLLLLLLL LL
X LLLLLLLLLLLLL LL
* ^ LLLLLLLLLLLLL LL
/ x'
» X""
L. '^
LEGEND:
L AT LOWER LIMIT
»-AT UWtfi LIMIT
B IN BASIS ( UNCONSTRAINED 1
LLLLL LLLLLL S
.LLLLLLL LLL f
XflMi^^^U. ^T
LLLLLLLLLLIT* ^m* .
LLLLLLLLLCL
LLLLLLLLLLL
\.
k
I
I
jl
Figure 25. Status of bounds on emission limits, future light industry
131
-------
LEGEND:
L AT LOWER LIMIT
»AT UPPER LIMIT
Figure 26. Status of bounds on emission limits, current heavy industry
132
-------
LLU.l-LU.LLLk
LLLLLLLLLLLL
LLLLLLLLLLLL
LLLLLLLLLLLL
LLLLLLLLLLLL
LLLLLLLLLLLL
LLLLi
. JLLLLLLLLLL
"LLLLLLUO.LLLLL
LtLLLLLLLLLLLL
J/
.-.J
LEGENDt
L AT LOWER LIMIT
• AT UPPER LIMIT
Figure 27. Status of bounds on emission limits, future heavy industry
133
-------
r
~n
.•.•U'UJLUVfc"llllll"'««
EfmmimmmiimmmiE;
•^;3i
:>S
/:
:S:
1
b.;
-J
MINIMUM
MAXIMUM
.00000
.00055
.00055
.00554
.00554
.02770
.02770
.27702
Figure 28. SOX emission density, base year (tons year
day per square kilometer)
134
-------
MINIMUM
MAXIMUM
.00000
.00035
.00055
.00554
.00554
.02770
.03770
.27702
.27702
.83106
Figure 29. SOX emission density, design year (tons per
day per square kilometer)
135
-------
I -
MINIMUM .00000
MAXIMUM .oooss
,00055
.00554
.00554
.02770
.02770
.27702
XX
XX
.27702
,83106
HH
MH
.83106
1.20000
Figure 30. SOX emission density, all sources at upper limit
(tons per day per square kilometer)
136
-------
upper limit the first receptor, R083, has an upper bound in the problem
of 61.98 micrograms per cubic meter. (The concentration at this receptor
due to non-EDZ sources is thus 80 - 61.98, or 18.02.) . The upper bound
on this receptor could be changed to any value between 61.97 and 62.03.
The total emissions; i.e., the objective function, would decrease or
increase, respectively, at a rate of 0.66 tons per day per unit change
in the receptor's upper bound. At the values 61.97 or 62.03, emission
source 9404021 (a light industrial point source) currently in the basis
would have its emission density limit set at, respectively its lower or
upper bound. For changes in the receptor's upper limit beyond 61.97 or
62.03, the rate of change of the objective function would be different
and other emission sources besides 9404021 would be affected. Similar
information is given for the other receptor at its upper limit.
The second part provides similar information for each column at the
limit level. For example, the first source given is 200012 (a commercial
area source, future development) is at its upper bound, 0.18790 tons per
square kilometer per day, in the optimal solution. It has a planned area
of 1.94 square kilometers. The sensitivity analysis indicates that if
its area were reduced to 1.15 square kilometers its emission density limit
would be set at 0.42 tons per square kilometer. Total regional emissions
would decline 0.80 tons per day per unit decrease in this source's emission
density limit. Any change in this source will affect the emission density
limit of source 9404021, currently in the basis. If source 200012's emis-
sion density limit were set below 0.42 tons per square kilometer, source
9404021's emission density limit would be set at its upper bound. Similar
information is presented for all of the emission sources.
The third part indicates the sensitivity of the solution to changes
in the ambient concentrations at receptors not at their upper limit. The
first receptor is at 33.47 micrograms per cubic meter in the optimal
solution. Its upper bound is 72.6 micrograms per cubic meter. However,
the upper bound can be changed down to 33.40 before an emission source
whose emission density limit currently at its upper bound will be forced
to change.
137
-------
The last part provides information on sources not at their upper or
lower bound. For example, source 9494921 is currently in the basis with
an emission density limit of 18.1 tons per square kilometer. Its upper
and lower bounds are 18.4 and 15.5, respectively. It can however, change
only between 18.2 and 18.0 before another source, currently at its upper
or lower bound, will be affected.
The examples given above are typical of the sensitivity of the solu-
tion to the right hand side and bounds of all the receptors and sources.
As was the case in Baltimore, it is highly sensitive.
Solution - 100 Percent Sample
The emissions maximization was also executed for the full complement
of receptors to determine the sensitivity of the optimal solution to dif-
fering receptor sets. With the full receptor set, the linear program
failed to reach an optimal solution. That is, with all sources set at
their lower limit, the NAAQS is still exceeded. Two receptors, not in-
cluded in the 40 percent execution, showed projected ambient concentra-
tions of, respectively, 10 and 6 yg/m3 in excess of the standard. Several
receptors in the vicinity of these receptors; i.e., 1 kilometer distant,
that were included in the 40 percent execution showed projected concentra-
tions of 15 to 20 yg/m3 less than the standard. Inspection of the emis-
sion inventory indicated several point sources with low stack heights in
the immediate vicinity of these two receptors.
Unlike the infeasible particulate matter example, in which nearly half
or the receptors over a broad area exceeded the NAAQS, the infeasible
100 percent sulfur oxide example is due to a hot spot with most receptors
well below the NAAQS. The appropriate action would be to prepare a sum-
mary of each source's contribution to the two receptors and review the upper
and lower bounds on the emission density for each significant source.
Some form of source specific emission limits may well be appropriate.
In conclusion, several observations can be made.
138
-------
The same pattern of emission density limits set in
Baltimore by the emissions maximization formulation
were set by this model in Louisville; i.e., more
stringent emission density limits in the center of
the region and less stringent ones on the periphery.
It is possible that even though the NAAQS is currently
being attained, the standard may be violated due to
planned growth even if all growth is at the specified
lower emission density limit. In such cases, the
linear program will reach an infeasible solution. The
emission density of any one or all three classes of
sources must then be reduced; viz., non-EDZ sources,
existing EDZ sources, and planned future development.
The model can be sensitive to the number and location of
receptors. In the extreme, an infeasible solution may
be reached with a full complement of receptors when a
feasible solution is reached with a subset of receptors.
If the receptor locations are valid, an analysis of
source contributions to the receptors exceeding the
standard should be used to further define the problem.
To conserve resources, it may be appropriate to run the
model with a subset of receptors. When a satisfactory
optimal solution has been achieved, the model must then be
rerun with all receptors. The preliminary optimal solu-
tion can serve as a starting point for the rerun with all
receptors.
139
-------
SECTION IX
REFERENCES
1. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 7: Projecting County Emissions. U.S. Environmental Pro-
tection Agency, Office of Air Quality Planning and Standards,
Research Triangle Park, North Carolina. January 1975.
EPA 450/4-74-008.
2. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 8: Computer-Assisted Area Source Emissions Gridding Pro-
cedure. U.S. Environmental Protection Agency, Office of Air
Quality Planning and Standards, Research Triangle Park, North
Carolina. September 1974. EPA-450/4-74-009.
3. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 13: Allocating Projected Emissions to Sub-County Areas.
U.S. Environmental Protection Agency, Office of Air Quality
Planning and Standards, Research Triangle Park, North Carolina.
November 1974. EPA-450/4-74-014.
4. Air Quality Display Model. TRW Systems, Inc. Prepared for U.S.
Department of Health, Education and Welfare, Public Health Service,
National Air Pollution Control Administration, Washington, D.C.
November 1969. 2 vols: PB-198~299 and PB-198~300.
5. Busse, A. D. and J. R. Zimmerman. User's Guide for the Climato-
logical Dispersion Model. Prepared for U.S. Environmental Pro-
tection Agency, December 1973. EPA-R4-73-024.
6. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 13. Op. cit. p. 133.
7. Ibid. See also: Guide For Compiling A Comprehensive Emission
Inventory (Second Edition.) U.S. Environmental
Protection Agency, Office of Air Quality Plan-
ning and Standards, Research Triangle Park,
North Carolina, p. 7-10 (APTD-1135)
140
-------
8. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 13. Op. Cit. Chapter 2.
9. Kennedy, A. S., et al. Air Pollution/Land Use Planning Project,
Phase II Final Report, Volume III: An Economic Comparison of Point
Source Controls and Emission Density Zoning for Air Quality Manage-
ment. U.S. Environmental Protection Agency, Office of Air Quality
Planning and Standards, Research Triangle Park, North Carolina.
May 1973. EPA-450/3-74-028c.
10. Pertinent Excerpts from the Preliminary Final Environmental Impact
Statement. Clark Maritime Center, Jeffersonville, Clark County,
Indiana, Ohio River Mile 597. U.S. Army Corps of Engineers,
Louisville District, Kentucky. No date given.
11. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 13. Op. cit. p. 134.
12. Air Quality Display Model. TRW Systems, Inc. Op. cit. p. 2-1
to 2-3.
13. Brown, Barbara D. and C. S. Lipaj. Implications of a Prevention
of Significant Deterioration Policy on State Growth Management.
(Presented at 69th Annual Meeting of the Air Pollution Control
Association. Portland, Oregon. 1976.) p. 3-4-
14. Larsen, R. I. A Mathematical Model for Relating Air Quality
Measurement to Air Quality Standards. U.S. Environmental Pro-
tection Agency, National Environmental Research Center. Re-
search Triangle Park, North Carolina. November, 1971. AP-89.
15. Larsen, R. I. An Air Quality Data Analysis System for Inter-
relating Effects, Standards, and Needed Source Reductions -
Part 2. J Air Pollut Contr Assoc. 24(6). Pittsburgh, Pennsylvania.
June, 1974.
16. Brier, G. W. Validity of the Air Quality Display Model Calibra-
tion Procedure. U.S. Environmental Protection Agency, National
Environmental Research Center. Research Triangle Park, North
Carolina. PB-218 716. March, 1973.
17. Brier, G. W. Statistical Questions Relating to the Validation
of Air Quality Simulation Models. U.S. Environmental Protection
Agency, National Environmental Research Center, Research Tri-
angle Park, North Carolina. EPA-650/4-75-010. March 1975.
18. IBM Corporation. Mathematical Programming System/360, Version 2,
Linear and Separable Programming — Users Manual (360A-CO-14x).
White Plains, New York, GH20-0476-2. October, 1971.
141
-------
19. IBM Corporation. Mathematical Programming System/360 — Applica-
tion Description Manual. White Plains, New York, 6H20-0136-4.
November 1969.
pUFfa; ••- T . • I1 .
20. Personal Communication, Laurel Slate, IBM Corporation, White Plains,
New York.
21. Cirrillo, R. R. et al. Air Quality Analysis Workshop: Volume 1 -
Manual. U.S. Environmental Protection Agency, Office of Air
Quality Planning and Standards, Research Triangle Park, North
Carolina. EPA-450/3-75-080a. November 1975.
22. Benesh, F. H. Growth Effects of Major Land Use Projects -
Volume II: Compilation of Land Use Based Emission Factors.
U.S. Environmental Protection Agency. Office of Air Quality
Planning and Standards, Research Triangle Park, North Carolina,
EPA-450/3-76-012b. September 1976.
23. Baumol, W. J. Economic Theory and Operations Analysis.
Prentice-Hall, Inc., Englewood Cliffs, New Jersey. 1972.
24. Vajda, S. Readings in Linear Programming. John Wiley and Sons,
Inc., New York, 1958.
25. Killer, F. S. and G. J. Lieberman. Operations Research. Holder-
Day, Inc., San Francisco, California. 1974.
26. Brail, R. K. et al. Emission Density and Allocation Procedures
for Maintaining Air Quality. U.S. Environmental Protection
Agency, Office of Air Quality Planning and Standards, Research
Triangle Park, North Carolina. EPA-450/3-75-079. June, 1975.
27. Laboratory for Computer Graphics and Spatial Analysis, SYMAP
Manual. Graduate School of Design, Harvard University, Cambridge,
Massachusetts, October, 1971.
28. Personal Communication with Mr. Dave Lueck, Division of Air Pollu-
tion, Kentucky Department of Natural Resources and Environmental
Protection, Frankfort, Kentucky.
29. Williams, J. D., et al. Air Pollution Emissions Related to Land
Area - A Basis for a Preventive Air Pollution Control Program.
National Air Pollution Control Association, Durham, North Carolina.
APTD-68-11. July 1968.
30. Roberts, J. J. and E. J. Croke. Land Use as an Organizational
Basis for Regional Air Resource Management. Proc 2nd Int Clean
Air Congress, edited by H. M. Englund and W. T. Berry, Academic
Press, New York. p. 1247. 1970.
142
-------
31. Roberts, J. J., E. J. Croke and S. Booras. A Critical Review of
the Effect of Air Pollution Control Regulations on Land Use Plan-
ning. J Air Pollut Contr Assoc. (25)5. Pittsburgh, Pennsylvania.
1975.
32. Guidelines for Air Quality Maintenance Planning and Analysis.
Volume 3: Control Strategies. U.S. Environmental Protection
Agency, Office of Air Quality Performance Standards, Research
Triangle Park, North Carolina. July 1974. EPA 450/4-74-003.
33. 1972 OBERS Projections. U.S. Water Resources Council, Washington,
D.C. April 1974.
143
-------
APPENDIX A
AN INTRODUCTION TO LINEAR PROGRAMMING
\
Linear programming is a mathematical technique whereby a linear
function of several variables is maximized (or minimized) when the vari-
ables are subject to various constraints. The linear function that is
maximized is referred to as the objective function; it is usually some
kind of economic or policy objective such as profit, output, or, in the
problem at hand, the maximization of the emissions of air pollutants or
the minimization of costs imposed on sources of emissions. For example,
the objective function of a linear program that maximizes emissions
would be:
maximize emissions, E' «= a, E, + a. E0 . . . . + a E (17)
1122 n n
where E, , E .... and E are the variables (the emission density limits
12 n
in tons per day per acre) and a , a , and a are the constants (the lot
area of each emission source). All the E's are of the first power, which
is necessary to characterize the relationship as a linear equation.
Each of the variables, E , E .... and E , are subject to some
12 n
constraint. In the problem at hand, the air quality at a specified recep-
tor must not exceed a specified level. Specifically:
t a E + t a E + t a E < S (18)
111 222 n n n —
where t., t., .... t are the transfer coefficients relating unit emissions
12 n
at each source to the ambient concentration at the receptor. The sum,
144
-------
over all sources, of t a E must be less than or equal to the
n n n
standard.
It is useful to consider the graphic solution to a simple problem
with two emission sources and two receptors. Assume the maximum ambient
concentrations desired at receptor 1 is 12 jig/m^ and at the second receptor,
o
8 mg/m . The areas and transfer coefficients are as follows:
Source 1
Source 2
Area, acres
6
7
Transfer coefficie
To receptor no. 1
0.33
0.43
yg/m
at, ,
tons /year
To receptor no.
2
0.33
0.14
The maximization problem would then be formulated as:
Maximize: 6 E + 7 E (19)
Subject to: 2 E. + 3 E <_ 12
2 E1 + E. < 8
j. z —
V E2 i °
In Figure 1 the two constraining equations are shown. All points to the
left of the inequality:
2 El + 3 E2 - 12
(20)
are possible (feasible) combinations of E and E . Similarly, points to
the left of the second constraint(inequality are feasible.
The hatched region in Figure 31 is the set of points that are feasible
under all constraints; it is the feasible region. Clearly, not all of
these points will maximize emissions.
145
-------
(0,8)
2E, + E2 < 8
= (0,4)
2E, + 3E,
-------
The family of objective functions is shown in Figure 32. Each line
is the set of points with constant amounts of aggregate emissions. Con-
sider the objective function for which total emissions, E, equals 18.
This line contains many feasible points, all of which give us 18 tons
per year of emissions. But the line E = 24 is better, because it contains
points with 24 tons per year. Best of all is the line E = 32. Here there
is only one point in the feasible region, D, and this is the optimal
solution. The optimal solution is to set an emission density limit of
3 tons per year per acre for emission source No. 1 and 2 tons per year per
acre for emission source No. 2. At any other combination of emission
limits, either the desired ambient air quality levels of 12 and 8 micro-
grams per cubic meter would be violated or the sum of emissions over the
two sources would be less than 32 tons per year.
For a more detailed and comprehensive treatment of the subject, the
reader is referred to any standard text on the subject, such as:
• W. Baumol. Economic Theory and Operations Analysis23
S. Vajda, Readings in Linear Programming*1*
F. Hiller, G. Lieberman, Operations Research25
These references also note more advanced texts.
147
-------
APPENDIX B
PROVISION FOR SHORT-TERM STANDARDS
IN AIR QUALITY CONSTRAINT
The air quality constraint in the emission density limit setting
model, equation (8), viz.,
'i.J.k.i (EJ,k,*,m «Jikfm> 1 Xi>
s, „ | for all !,«, (21)
1, X.
states that the sum, over all sources, of the product of the emissions
from each source and a transfer coefficient between that source and a
specified receptor i shall not exceed the air quality standard at that
receptor. The transfer coefficients, derived through the execution of
AQDM, relate the effect of annual emissions at a specific source to the
ambient concentration at a specified receptor i. Thus, through the use
of this relatively straightforward constraint, it is possible to set a
pattern of emission density limits such that the annual ambient air qual-
ity standards, expressed as arithmetic means, are not violated.
One problem with this formulation, however, is that the annual stand-
ard for particulate matter, 75 yg/m^, is an annual geometric mean. For
the receptors in the EDZ model that are monitoring stations, it is possible
*
The variable definitions of Section V are employed.
148
-------
to convert the geometric mean to an expected arithmetic mean that is
equivalent to the annual particulate matter standard, if one assumes a
log-normal distribution, viz,
2 *
m = m exp (0.5 In s ) (22)
where m = arithmetic mean
m = geometric mean
s = geometric standard deviation
O
For example, the geometric standard deviation at the St. John's
College monitoring station in the Baltimore AQMA is 1.5 (see Table 8),
Substituting in Equation (22), one obtains
m = 81.4 - 75 exp (0.5 In2 1.5)
(23)
Ignoring the effects of calibration, one would substitute 81.4 in
place of X- o f°r the. *• corresponding to the St. John's College receptor
i, x,
and the H corresponding to particulate matter.
The question remains, however, of how to treat the majority of the
receptors in the EDZ model that are not monitoring stations and therefore
have an unknown geometric standard deviation. If a geographic pattern
existed among the observed geometric standard deviations, it would be
For the convenience of the reader familiar with A Mathematical Model for
Relating Air Quality Measurements to Air Quality Standards,14 EPA Publica-
tion AP-89, the indicated equations in this appendix correspond to rear-
rangements or combinations of the following equations in AP-89.
Equation number AP-89
22
24
25
26
149
-------
Table 8. 1974 TSP CONCENTRATIONS RECORDED AT MONITORING STATIONS IN THE
BALTIMORE AQMA
Monitoring stations
St. John's College
Anne Arundel Co.
Overlook School
Anne Arundel Co.
Arundel High School
Anne Arundel Co.
Glen Burnie Public Works
Anne Arundel Co.
Route 176 & 170
Anne Arundel Co.
Harwood High School
Anne Arundel Co.
Fire Dept. HQ
Baltimore
Sun Ave. Fire Dept.
Baltimore
NE Police Station
Baltimore
Number of
observations
53
53
56
62
53
48
107
79
100
Maximum
24-hr
135
184
169
196
182
137
282
324
175
Second
highest
126
136
115
162
153
119
253
317
143
Arithmetic
mean
51
66
52
80
65
47
117
148
60
Geometric
mean
47
60
47
65
58
43
102
134
53
Geometric
standard
deviation
1.5
1.5
1.6
2.7
1.6
1.5
2.0
1.6
1.6
-------
Table 8 (continued). 1974 TSP CONCENTRATIONS RECORDED AT MONITORING STATIONS IN
THE BALTIMORE AQMA
Monitoring stations
NW Police Station
Baltimore
SE Police Station
Baltimore
SW Police Station
Baltimore
3001 Eastern Blvd.
Baltimore Co.
Fort McHenry
Baltimore
1030 Linwood Ave.
Baltimore
Sewage Disposal Plant
Baltimore
Ft. Howard VA Hospital
Baltimore Co.
Garrison Police Barracks
Baltimore Co.
Number of
observations
91
73
140
53
54
58
52
114
47
Maximum
24 -hr
198
309
309
172
228
265
291
182
123
Second
highest
151
287
294
113
221
263
266
154
122
Arithmetic
mean
74
120
96
57
114
108
126
55
55
Geometric
mean
68
106
86
49
103
95
113
50
49
Geometric
standard
deviation
1.5
1.7
1.6
2.1
1.6
1.7
1.6
1.6
1.7
-------
Table 8 (continued). 1974 TSP CONCENTRATIONS RECORDED AT MONITORING STATIONS IN
THE BALTIMORE AQMA
Monitoring stations
Catonsville Library
Baltimore Co.
County Office Bldg.
Harford Co.
Dundalk
Baltimore Co.
Essex
Baltimore Co.
Whiteford
Harford Co.
Simpsonville
Howard Co.
Goucher College
Baltimore Co.
Carrol Co. Hospital
Carrol Co.
Number of
observations
47
53
108
53
58
61
76
52
Maximum
24 -hr
141
130
203
158
96
145
135
227
Second
highest
91
122
176
155
86
116
110
122
Arithmetic
mean
51
56
84
69
39
51
54
50
Geometric
mean
46
51
76
63
35
46
48
43
Geometric
standard
deviation
1.6
1.5
1.6
1.5
1.6
1.6
1.7
1.7
Ui
S3
Source: SAROAD data system.
-------
possible to use that pattern to assign deviations to the receptor grid
in EDZ. A plot of the location of the monitoring stations and geometric
standard deviations is shown in Figure 33. In our view a convincing
pattern does not exist. If a pattern does not exist, an alternative
is to use the mean of the recorded data, in this case 1.67. Substituting
in equation (22), one calculates a value for the expected arithmetic mean
3
of 85.54 yg/m . Therefore,
every receptor in the grid.
3
of 85.54 yg/m . Therefore, this value would be used as the standard at
In this example, the observed geometric standard deviation, obtained
*
from SAROAD printouts, was utilized. In the actual execution of the
decision model it is more appropriate to compute the geometric standard
deviation using the following formulae.
Sg = (c/m> (2A)
where c is the maximum 24-hour observations, and
z is the z deviate of the normal distribution corresponding
to a frequency calculated by
f = 1 - ±~^- (25)
where n = the number of valid 24-hour observations.
The use of formulae 24 and 25 allows the calculation of a geometric
standard deviation from the top half of a distribution, which most closely
fits the frequencies associated with the standard.
*
Storage and Retrieval of Aerometric Data.
153
-------
Ai.r
CARROLL COUNTY
HARTFORD COUNTY
Al.B
ANNE ARUNOEL COUNTY
Ai.e
Figure 33. Location of TSP monitoring stations and associated
geometric standard deviation
154
-------
The assumption of a log-normal distribution also allows the consid-
eration of the 24-hour standards in the emission density limit setting
model. The ability to make this assumption is critical, in that there
is no other way to consider the short-term standard and maintain a rela-
tively simple EDZ model which would be easily implemented.
The treatment of the short-term standards would proceed in a manner
analogous to the estimation of an arithmetic mean equivalent to the geo-
metric mean particulate standard. For example, in the case of sulfur
dioxide, the average geometric standard deviation measured in the Baltimore
AQMA is 2.76 (see Table 9). Assuming log normality, the second highest
24-hour concentration would be related to the arithmetic mean as follows.
z - 0.5 In s ,„,..
C2ndmax=msg g (26)
The appropriate z deviate would be 2.62. The 24-hour sulfur dioxide
standard is 365 yg/m3. Accordingly, rearranging Equation (26) to solve
for m,
m = 365 s °'5 ln sg - 2'62 (27)
O
Substituting 2.76 for s , one obtains an arithmetic mean of
o
42.7 yg/m . Since this is more stringent than the annual standard of
o
80 yg/m , it would be used in Equation (21) for x. „ where £ corresponds to
»
sulfur oxides.
The use of the second highest 24-hour concentration in a year as the
"design value" means that the standards will be violated (second highest
concentration exceeding standard) on the average of every other year.
As an alternative, it has been suggested-^ that the annual maximum con-
centration be used as the design value. If this is done, the ambient
short-term standards would probably be violated in 1 year out of 8. A
Z deviate of 2.94 would be used instead of the 2.62 used above.
155
-------
Table 9. 1974 SO CONCENTRATIONS IN THE BALTIMORE AQMA
Monitoring stations
St. John's College
Anne Arundel Co .
Arundel High School
Anne Arundel Co.
Glen Burnie Public Works
Anne Arundel Co .
Route 176 & 170
Anne Arundel Co.
Harwood High School
Anne Arundel Co.
Calvert & 22nd Sts.
Baltimore
3001 Eastern Blvd.
Baltimore
Number of
observations
51
50
53
55
49
51
43
Maximum
24 -hr
152
150
162
121
61
120
65
Second
highest
74
92
118
64
58
109
33
Arithmetic
mean
20
23
30
18
13
26
8
Geometric
mean
9
13
17
11
8
13
5
Geometric
standard
deviation
3.4
3.2
3.1
2.9
2.9
3.6
2.3
Ui
-------
Table 9 (continued). 1974 SO CONCENTRATIONS IN THE BALTIMORE AQMA
Monitoring stations
Fort McHenry
Baltimore
1030 Linwood Ave.
Baltimore
Sewage Disposal Plant
Baltimore
Ft. Howard VA Hospital
Baltimore Co.
Garrison Police Barracks
Baltimore Co.
Catonsville Library
Baltimore Co.
County Office Bldg.
Harford Co.
Number of
observations
54
55
52
108
44
45
61
Maximum
24 -hr
169
106
284
76
41
49
80
Second
highest
166
102
262
68
31
26
65
Arithmetic
mean
77
24
70
18
9
8
19
Geometric
mean
64
13
50
11
6
5
12
Geometric
standard
deviation
2.0
3.1
2.4
2.7
2.5
2.3
2.6
-------
Table 9 (continued). 1974 S() CONCENTRATIONS IN THE BALTIMORE AQMA
X
Monitoring stations
Dundalk
Baltimore Co.
Essex
Baltimore Co.
Whiteford
Harford Co.
Simpsonville
Howard Co.
Goucher College
Baltimore Co.
Carrol Co. Hospital
Carrol Co.
Number of
observations
108
56
58
59
73
48
Maximum
24 -hr
169
110
35
54
117
54
Second
highest
156
98
28
39
60
48
Arithmetic
mean
33
29
7
12
18
14
Geometric
mean
16
20
5
9
11
9
Geometric
standard
deviation
3.7
2.6
2.3
2.3
2.7
2.6
Ul
00
Source: SAROAD data system.
-------
In a similar fashion, the annual arithmetic mean equivalent to the
24-hour particulate standard can be computed. Using Equation (27), one
obtains
m = 150 s °'5 ln Sg - 2'62 (28)
O
Substituting 1.67 (the average in Baltimore) for s , one calculates a
O
value of 44.6 for the expected arithmetic mean equivalent to a 24-hour
standard of 150 yg/m3. As this also is more stringent than the expected
arithmetic mean, 85.5, equivalent to the annual standard, the value of
44.6 yg/m3 would be used in equation (21) for x. „ where K, corresponds to
i, x,
particulates.
159
-------
APPENDIX C
THE SO -PM FEASIBILITY CONSTRAINT
x
The sulfur oxides emission density-particulate matter emission density
feasibility constraint is based on the premise that sulfur oxide emis-
sion densities and particulate matter emission densities are not inde-
pendent. The utilization of this constraint in the decision model will
avoid the assignment of a certain emission density to a component area
for one pollutant and the subsequent assignment of an emission density
for the other pollutant that is so stringent that the emission density
for the former pollutant could never be utilized.
Assuming that the respective emission rates are not independent, more-
over, that they are linearly related, one could expect to obtain a scatter-
gram such as that shown in Figure 34. The regression line estimating this
relationship, PM = a + b SO , is shown as a solid line. The intercept,
-\ A *
a, and the slope, b, are estimated values; one can compute a confidence
interval for the regression line. In particular,
(29)*
For a given sulfur oxide emission rate, X , one would then constrain
the particulate matter emission rate, Y , to be greater than or equal to
*
See, for example, J. Johnston, Econometric Methods, McGraw-Hill. 1972.
160
-------
(0
2
O
to
I
UJ
s
0.
-PM
SOX EMISSIONS
S°X,
Figure 34. Hypothetical relationship between
SO and PM emissions
x
161
-------
the lower boundary. However, since X appears underneath the radical, this
cannot be easily employed in a linear program. It is an expedient approx-
imation, though of no statistical meaning, to use the confidence interval
for the mean of Y given the mean of X. Specifically,
1 '
e/2
(30)
The dotted lines in Figure 34 show what the extremes of this interval
may look like. The constraint equation would then be
j,k,m
(3o)
which reduces to
a -
iff
j,k,m
(31)
denoting the slope term in parentheses as y, , and the term in brackets
a k
as y, , this can be expressed as
K.
b a
j,k,£=l,m ~ yk Ej,k,Jl=2,m - pk
(33)
which is the first SO -PM feasibility constraint in the emissions maximi-
X
zation decision model described in Section V.
In the execution of the model, one would specify the confidence level
of the interval by selecting a specific t, say 80 percent.
See Johnson, op. cit.
162
-------
As was stated previously, the necessity of this constraint depends on
the existence of a relationship between sulfur oxide and particulate matter
emissions, after the application of existing regulations. If such a
relationship does not exist, then this constraint is superfluous. The
setting, in a specific component area, of a relatively high sulfur oxide
emission density and low particulate matter emission density, or the
converse, would only encourage those sources that have such emission
characteristics to locate in that component area.
The existence of such a relationship would have to be tested for each
time EDZ is implemented in an AQMA. Whether it existed would depend on
the characteristics of the emission sources in an AQMA of a particular
land use category and the homogeneity of these sources. One should note
that this independence may exist only for some of the land use categories.
An inspection of the point source emission inventory in Baltimore and
Louisville would suggest that, in these AQMAs, emission rates for the
two pollutants are independent in the case of the heavy industrial land
use category. Sulfur oxide and particulate matter emission rates for
this category are plotted in Figure 35. Each asterisk represents one
point source and a number represents the position of multiple sources.
The lack of a linear relationship between sulfur oxide and particulate
matter emission rates is adequately documented by the correlation coeffi-
cient. Figure 36 shows the relationship between emission rates for all
heavy industrial point sources emitting less than 100 tons per year.
No linear relationship is evident. (In fact, it could be more accurately
described as a hyperbolic or reciprocal relationship, rather than linear.)
Light industrial point source emission rates are plotted in Figure 37.
Again, there is little linearity evident. Figure 38 shows the relation-
ship between emission rates for noncoal burning commercial point sources.
2
The R is 0.7, which would suggest a relationship exists. The regression
line has been drawn in as a solid line; however, it appears to have been
greatly influenced by a half dozen relatively large sources. A more
appropriate regression line is shown as a dashed line.
163
-------
SCATTfRSR*"! 0
45000.00
40500.00
36000.00
31900. OU
27000.00
Z250B.OO
11000.00
13500.00
9000.00
HOT. 00
0.0
P 100*4) PM (ACROSS) SOX
361.35 1084.05 1B06.TS 252*. 45 3252.15 3*74.85 46*7.5
* :
•
-.-..-- - - -
*
i* •
»999*7 55 * • •
r •• •••• ' • --
•
i
i
i *
5 5420.25 6142.95 6*65.65
«
I
I
I
1
' ' '— •'•" " ' «
I
-'- ' I
1
+
• r
i
i
t
i
»
4SOOO.OO
40SOXT.OO
36000.00
" **W**"W
27000.00
i
------ — -- • - --j -• - -
i
t
i
• "• - •-• "i
»
•' •" ' ~ " . ~ 7 '
I
I
I
' - " T~
I
*
. - •-• I
I
I
I
I
I
9 »»
0.0 722.70 1445.40 2161.10 2890.30 1M1.50 4336.20 50S8.90 5781.60 6504.10 7227.00
STATISTICS..
CORRELATION IRI- 0.00819 R SQUARED - O.OOOOT SIGNIPICANCt
STO ERR OF EST - 1343.95656 INTERCEPT (A) - 72.00005 SLOPE
-------
SCATTERSRAN Of IDOMNI P«
J.OO 15.00
2S.OO
(ACROSS) SOX
35.00 65.00
75.00
•5.00
20.00
lO'.OO
•* *• * *~
*
•
*
2 2
3
6
* «
13
14
14 .
I*
44
12
I
I* ...
I
»Z . i
[•
13
iV
[ 3 — • * " ~—
12 2
• 2
19
19 2
13
[*2
•a
13
19 2
13 2
13 3
48
19
1922
192 44 4
19*962256 2 t
»99«93 "28 *
o.o 10.06
1
r
i
i
i
i
i
t
i
i
I
i
i
i
i
i
i
i
i
' t ' "
i
i
i
i
i
I
i
i
i
i
i
,
i
i
i
i * .
i
i •
2 I •
2 I
I
I
2 2 I
2 t 2
« 61
9 I
* 24 13
22 2 I •
20.00 30.00 40.00 50.00 60.00
.
•
.:
7
i
4
..
... ... _
2
1
I
I
• I
' . *
70.00 10.00 90.00 100.0
90.00
STATISTICS..
CORRELATION CHI-
STD ERR OF EST -
PIOTTEO VALUES -
0.0673*
14.499S5
1075
It SQUARED
INTERCEPT (A) -
EXCLUDED VALUES-
0.00453
6.20879
>2
SIUUFICANCI
SLOPE IB)
. MISSING VALUES
0.019*3
0.10946
. _ 0
Figure 36. Plot of Baltimore heavy industrial point source PM
rates versus SOX emissions rates, sources with emis-
sions less than 100 tons/year (note: '9' indicates
9 or more observations)
165
-------
Of COOK*I M (ACKOSII SOX _
67.41 101.35 317. JS 471.II 607.0S 741.9« 876.89 IOII.T5 1144.6» IMI.tS.
84O.OO
756.00
672.00
588.00 "
504.00
420.00
'
336.00
252.00
.
168.00
84.00 <
1
0.0 •
1
1
*
.
3
3
3
» *
2 4
3
9» 2
94» 2 39
994927 4
.0 134.90 269.80 404.70
I I
t
I
I
I
I
I
I . .
I ._'.
I
I
I
I
I
I .
I
I
. I . — _. ._. .
i . .. .T. :
i
i
i
t
i
i
i
i
i
i
i . .
i
.1 ... . ...
i
i
i
i
i
i
2 T *
I
I
I
* I .
I . 5.
539.60 674.50 809.40 944.30 1079.20 1214.10.. 134.9.
_..
.".
.
10
M4.00
J«4.00
J36_.00
84.00
STATISTICS..
CORRELATION («)-
STD ERR OF EST -
PLOTTED VALUES -
0.61210
53.08410
1007
R SQUARED
INTERCEPT (A) -
EXCLUDED VALUBS-
0.37467
5.66194
.0
SIGNIFICANCE
SLOPE (ft)
MISCINC VALUES -
0.00001
0.26243
0
Figure 37. Plot of Baltimore light industrial point source PM
emission rates versus SOX controlled emission rates
(tons/year) (note: '9' indicates 9 or more observations)
166
-------
SCATTERGMN Of IOOHMI
m.os
is.oe
22.50
17.10
JJupo
*.»S 10.15 M.T9 «».»» 42. S» 7».«5
27.10
55.60
20 125.10 U».00
STATISTICS..
CORRELATION -
STO ERR OF EST -
PLOTTED VALUES -
0.83675
1.67332
i»»
0 SQUARED
INTERCEPT (A) -
EXCLUDED VALUE5-
O.TOOIS
0.94229
0
SIGNIFICANCE
SLOPE IB)
MISSING VALUES -
_
o.ooooi
0.14413
0
Figure 38. Plot of Baltimore commercial point source PM emission
rates versus SOX controlled emissions rates, excluding
coal burning sources (tons/year) (note: '9' indicates
9 or more observations)
167
-------
In summary, a relationship between emission rates should be tested
for on a sample from the base year emission inventory, stratified
according to the land use categories contemplated for emission density
regulation. If one exists, it is prudent to include a feasibility con-
straint, as discussed above, in the decision model.
168
-------
APPENDIX D
SPECIFICATION AND UTILIZATION OF CONTROL COST PARAMETERS
Both formulations of the control cost minimization objective func-
tion (i.e., model nos. 2 and 3 of Section V) require an estimate of
incremental control cost functions for future emission sources. The
ability to estimate these functions is limited by two problems — the
need to forecast adequately the nature of future emission sources and the
need to estimate control cost functions that can be used in the linear
programs. The discussion below is limited to industrial sources, as the
control cost linear programs would likely be applied only to that source
category.
At best, one can identify those areas that are planned for heavy and
light industry. Within each of these categories, regional employment pro-
*
jections would give some indication of the most likely SIC code categories,
probably, at best, on the two-digit level. One could either use the
industry with the largest employment growth or, within the top 5 or 10
categories, the heaviest polluting industrial category as a prototype
heavy (or light) industry. Using this prototype, one would endeavor to
estimate the incremental control costs.
For example, in Baltimore the U.S. Water Resource Council projections,^
shown below in Table 10, indicate that primary metals (SIC code 33) and
chemicals (SIC code 28) are the two industrial categories that will experience
*
Standard Industrial Classification.
169
-------
the largest absolute growth. One could use the control costs associated
with the primary metals industry as the "c" parameters for future heavy
industrial sources. On the other hand, most of the growth in SIC code 33
may represent the utilization of existing capacity and the capacity ex-
pansion of existing sources. In such a situation, it would be appropriate
to use another industrial class.
Table 10. 1980 TO 1990 PROJECTIONS OF GROWTH
IN MANUFACTURING EARNINGS ($000)
Heavy
26
27
28
29
30
31
32
33
industry
28,000
55,400
65,500
1,400
85,100
20
22
23
24,
34,
35
36
371
37
21,
Light
25
19
ex 371
30-32,
industry
12,300
11,000
12,700
12,500
65,300
69,700
26,800
28,000
10,700
38, 39 58,700
Source: Reference 33.
For the purposes of this example, we will assume that it is more
appropriate to use SIC code 28, which is Chemicals and Allied Products.
However, projections more detailed than two-digit SIC codes are needed
to develop a control cost parameter. Unless more detailed employment
projections are available, it would be consistent with AQMP emission pro-
jection methodologies to assume the characteristics of future sources to
be similar to those of existing sources. Therefore, one would review the
emission characteristics of existing heavy industrial point sources to
develop the incremental cost of emission reductions. Given the specific
characteristics of several prototypical sources, the control cost element
of the Implementation Planning Program (IPP) package could be utilized to
estimate the costs of an emissions reduction.
170
-------
One would
1. Develop projections of regional industrial growth by at least
two-digit SIC code. The OBERS projections or equivalent or
superior projections available locally can be used.
2. For each two-digit SIC code, estimate the percentage of future
growth occurring from:
a. the utilization of existing industrial capacity
b. capacity expansion of existing sources
c. construction of new industrial facilities.
This information can be derived from the regional economic
projections and a survey of existing point sources conducted
for the AQMP emission projection.
3. Rank order the existing point sources in each two-digit SIC
code by their emission density after application of existing
and projected regulations.
4. For both the heavy and light industrial categories, deter-
mine the probable two-digit SIC codes for new industrial
facilities, based on the results of Step 2.
5. Identify prototypical point sources in the heavy and light
industrial categories by selecting the existing point source
with the highest emission density in the fastest growing
two-digit SIC code. If, in either the heavy or light in-
dustrial category, growth in more than one two-digit SIC
code is expected to contribute significantly to the con-
struction of new industrial facilities, identify a proto-
typical point source for each one.
6. For each of the two or more prototypical point sources and
for the expected capacity expansion at existing sources,
execute the control cost segment of the IPP package and
identify the annualized control costs for possible levels
of emission reduction beyond existing and projected
regulations.
The costs of specific levels of emission reduction cannot be directly
used in equation (9). There is, however, a more complex formulation of
a control cost minimization that employs integer programming, discussed
in Section V, in which they can be used.
171
-------
Before reviewing the specification and use of the cost parameter, it
is convenient to develop a specific example. In the Baltimore AQMA, it
was assumed that SIC code 28 was the fastest growing heavy industrial
category. The first SIC code 28 point source in the Baltimore inventory
is the Glidden Durkee paint plant, SIC code 2816. The particulate matter
emission sources are summarized in Table 11. Table 11 also shows the
various levels of emission reductions available by installing control
equipment and the annualized cost of that equipment. The possibility of
tandem devices is ignored. Based on this information, a schedule of
control costs and emission reductions was prepared and is shown in
Figure 39.
To use this information to specify the cost parameter, it is necessary
to develop a linear approximation. A possible one is shown in Figure 40.
Three line segments are shown with slopes of 13.86, 68.93, and 3528.08.
Since the control cost — emission reduction schedule cannot be adequately
approximated with a single line segment, one is faced with several cost
parameters that could be used in the objective function.
The only feasible way to accomplish this, given the control cost
minimization formulation in equation (9), is through successive solu-
tions to the linear program.
For example, in the case of three future heavy industrial sources of
the same area, one would first execute the following linear program.
Max E' = 13.86 a-,^ EI + 13.86 &2 E2 + 13.86 a.j E3 (34)
subject to:
air quality constraint: £ t. . a. E; £ S. | for all i
.: i>J J J 1
172
-------
Table 11. EMISSION REDUCTIONS AND ASSOCIATED CONTROL COSTS
u>
Source and control device
Calcination kilns (three)
Uncontrolled
Low efficiency cyclone
Medium efficiency cyclone
Lou efficiency scrubber
Medium efficiency scrubber
High efficiency scrubber
High temperature fabric filter
Rotary dryers (three)
Uncontrolled
Low- efficiency scrubber
-Low temperature fabric filter
Ball mills (three)
Uncontrolled
Low temperature fabric filter
Boilers (five)
-Distillate oil
Natural gas
Ore sulfation (eight)
Uncontrolled
Total uncontrolled emissions
Device
efficiency
0
60
75
80
90
98
99
0
80
99
0
99
Particulate
emissions,
tons/year
1,500
600
375
300
150
30
15
333
66
3
400
it
29
2.5
1
6,852
Capital
cost, $
0
15,700
15,700
60,000
60,000
60,000
114,000
0
11,100
6,900
0
6,200
0
0
Annual
operating
cost, $
0
5,042
7,183
24,508
38,379
58,211
31,113
0
11,956
11,340
0
11,347
0
93,847
Total annualized
cost, ignoring
tax effects, $
0
6,876
9,017
31,158
45,388
65,220
44,431
0
13,252
12,146
0
12,071
0
93,847 .
Emission
reduction
0
900
1,125
1,200
1,350
1,470
1,485
0
267
330
0
396
0
26.5
Average
control cost
per ton
reduction, S
0
7.64
8.02
26.27
33.62
44.37
29.92
0
44.63
36.81
0
30.48
0
3,541.40
-------
1700,000
8600,000 -
8500,000 -
o 8400,000 -
o
K
u 8300,000 -
8200,000 -
8100,000 -
1000
2000 3000 4000 5000
EMISSION REDUCTION, ton/year
6000 7000
Figure 39. Schedule of emission reductions versus control costs
for prototypical heavy industrial point source
174
-------
8675,179 -
o
u
o
a:
I-
o
o
$205,944 -
$63,264 -
4563
EMISSION REDUCTION, tons/year
6633 6766
Figure 40. Plot of linear approximation of the
control cost - emission reduction
schedule
175
-------
technological constraint: a E > 0
J J ~
zero cost constraint: a. E. <^ 4563
for all j
If, in any source j, the upper bound on E. was binding, it may be
possible to achieve a further emission reduction at that source. However,
the incremental cost would no longer be 13.86; it would be 68.93. There-
fore, one would rerun the program, substituting the parameters of the
second line segment. In this instance, suppose the upper bound on the
third E. was binding; one would then rerun the linear program as follows.
Max Ef = 13.86
+ 13.86
+ 68.93
E3 + 63264
(35)
subject to:
air quality: £ t. . a. E. < S. | for all i
• 1>3 J 3 — •*•
technological
constraints
al El —
a2 E2 -
a3 E3 > 0
al El ^- 4563
zero cost
constraints
a
—
< 6633 - 4563 = 2070
The technological and zero cost constraints are specified for the nodes of
a line segment on Figure 40, plot of linear approximation of the control
cost-emission reduction schedule.
176
-------
While one could use ordinary linear programming and rerun the model
until all the upper and lower bounds were not binding or until the opti-
mum solution did not change, it is expedient to use a spatial formula-
tion of linear programming that will do this all at once; viz., separ-
able convex programming.18
177
-------
APPENDIX E
THE APPLICATION OF THE EMISSION DENSITY SETTING MODEL IN AN
ATTAINMENT SITUATION
The models described in Section V and illustrated in Sections VII
and VIII assume that all current emission sources, after the applica-
tion of current or pending emission control regulations, will not be
required to further reduce their emission rates under emission density
*
zoning. Emission density zoning is thus a technique for managing
growth, both in terms of new emission sources and expansion of existing
ones, as contemplated by air quality maintenance planning. Attainment
of the NAAQSs is assumed.
In regions where the NAAQSs have not been attained or in those re-
gions where the current ambient air quality levels are so close to the
NAAQSs so as to preclude growth without reductions in emission levels of
current sources, the models described in Section V are not applicable.
The former situation will be apparent; the latter may not be so obvious.
It is possible that the linear program will be unable to find a feasible
solution; that is, if the region experiences the growth contemplated
in the land-use plan, the NAAQSs will be violated, even if all new sources
emit at minimum levels and current sources do not expand. To achieve
the NAAQSs it is necessary that current sources reduce their emission
rates beyond that contemplated by existing emission control regulations.
*
Specifically, note the attainment constraint of Model No. 1 on page V-3.
This is what happened in the Louisville particulate matter test case.
178
-------
It is possible to implement emission density zoning in such a
situation, though EDZ is more commonly thought of as a maintenance
strategy than an attainment strategy. One has two options. The first
is to impose additional emission control regulations on existing sources.
If sufficient reductions are imposed it will then be possible for the
linear program to find a feasible solution.
Alternatively, one can impose emission density limits that require
existing sources to reduce their emission rates. To set such emission
density limits in a model such as described in Section V, the lower
bound on the emission density limit of existing sources should be set at
a lower value. As formulated in model no. 1, the attainment constraint,
viz.,
Ej,k,£,mle5?I?*,m I for a11 J.k.*; for » - 1 (36)
restricts the emission density for current sources to be greater or equal
bfl SG '
to e , their current value adjusted to show compliance with applicable
regulations. This base or floor for current sources must be respecified
at a lower yet still feasible value.
It is unlikely that one would set the minimum emission levels for
current sources as low as those of future sources. With older plant and
equipment, current sources would find it more difficult to meet such
standards as new development. Also, the options of locating on a site
with a higher emission density or purchasing additional land are not as
easily available to the current emission source as they are to a future
source.
Within the range between the current emission density and the lower
limit that is used for future sources, one must specify a lower bound.
The most straightforward way of doing this is to:
179
-------
a. Catalogue the possible emission reductions for each
current source category. This would basically be
various types of control equipment, process modifi-
cation or fuel switching.
b. Convert the maximum emission reduction to a minimum
emission density. This would have to be done separately
for each component area and point source, as densities
of development will differ.
c. Use this value as the e ase parameter of equation (36).
This method has two shortcomings. The first is that the cost of
control is not considered if the emissions maximization formulation of
the model is utilized. The second is that it is assumed that emissions
can be abated to any point down to the minimum level. Neither short-
coming is particularly serious, the former can be resolved by using the
control cost formulation of the model, the latter is not serious since
most sources will have emission density limits set at their lower or
upper bounds.
However, combining the attainment and maintenance problems into one
model sidesteps the issue of how much emissions growth will be permitted
at new sources and, concurrently, how much emissions reduction will be
required at existing sources. If it is not the intent or policy to
treat the two classes of sources identically, it is more satisfactory to
separate the two problems and treat the trade-off between growth and
current sources openly.
180
-------
TECHNICAL REPORT DATA
(Please read Inunctions on the reverse before completing)
. REPORT NO.
EPA-450/3-77-006
2.
3. RECIPIENT'S ACCESSION-NO.
4. TITLE AND SUBTITLE
EMISSION DENSITY ZONING
5. REPORT DATE
March 1977
6. PERFORMING ORGANIZATION CODE
7.AUTH0R(s)
Fran|< H< Benesh> phil
Mills, Robert M. Patterson
) Michael T.
8. PERFORMING ORGANIZATION REPORT NO.
GCA-TR-77-01-G
9. PERFORMING ORGANIZATION NAME AND ADDRESS
GCA Corporation; GCA/Technology Division
Burlington Road
Bedford, Massachusetts 01730
10. PROGRAM ELEMENT NO.
11. CONTRACT/GRANT NO.
Contract 68-02-1376, TO 26
12. SPONSORING AGENCY NAME AND ADDRESS
Environmental Protection Agency; Office of Air Quality
Planning and Standards; Strategies and Air Standards
Division; Land Use Planning Office
Research Triangle Park, North Carolina
13. TYPE OF REPORT AND PERIOD COVERED
Final
14. SPONSORING AGENCY CODE
15. SUPPLEMENTARY NOTES
16. ABSTRACT
Emission density zoning is an air quality control strategy whereby all parcels of
land within an air quality maintenance area, excluding public rights of way and
reserved open space, are assigned maximum legal emission allowances for particulate
matter and sulfur oxides, expressed in terms of mass of pollutant per time period per
lot area. The subject of this publication is the development of criteria for
(a) partitioning an air quality maintenance area into component areas, (b) the
setting of emission density limits for all land uses in each component area and,
separately, all major emission sources, and (c) the revision of the emission density
limits as conditions change. The criteria are then translated into a generalizable
decision model, a linear program, for setting the emission density limits. The
criteria and model are tested in two air quality maintenance areas, Baltimore and
Louisville.
17.
KEY WORDS AND DOCUMENT ANALYSIS
DESCRIPTORS
b.lDENTIFIERS/OPEN ENDED TERMS C. COSATI Field/Group
Area emission allocations
Flux Density
Land Development
Land Use
Linear Programming
Particulates
Planning
Regional Planning
Sulfur Oxides
Urban Planning
Zoning
Air Pollution Control
Air Quality Maintenance
Emission Density Zoning
18. DISTRIBUTION STATEMENT
Unlimited
19. SECURITY CLASS (ThisReport)
UNCLASSIFIED
21. NO. OF PAGES
192
20. SECURITY CLASS (Thispage)
UNCLASSIFIED
22. PRICE
EPA Form 2220-1 (9-73)
181
------- |