P.1 

OAR Box 1126 

Prepped by Ollie Stewart 

Document Number: 

169)11-1-382 

Docket Number: 

A-83-22 

Printed 4/18/2008 9:58:29 AM Sheet 169 of 182 
P.2 
p.; 

Final Report 

AN EXAMINATION OF AIR QUALITY 

DATA COMPLETENESS REQUIREMENTS 

SYSAPP-84/076 . 

29 May 1984 

Prepared for 

Monitoring and Reports Branch 

Office of Air Quality Planning and Standards 

U.S. Environmental Protection Agency 

Research Triangle Park, NC 27711 

Prepared by 

A. D. Thrall 

J. L. Baptista 

C. S. Burton 

Systems Applications, Inc, 

101 Lucas Valley Road 

San Rafael, CA 94903 

68306 ef 831t6r 
P.4 

O 

o 

DISCLAIMER 

O 

This report has been reviewed by the Office of Air Quality Planning 

and Standards, U.S. Environmental Protection Agency, and approved for 

publication. Approval does not signify that the contents necessarily 

Q reflect the views and policies of the U.S. Environmental Protection 

Agency, nor does mention of trade names or commercial products constitute 

endorsement or recommendation for use. 

O 

O 

O 

O 

O 

O 

n 

O 
P.5 

ACKNOWLEDGMENTS 

The authors are grateful to Scott Chaplin for his research on current data 

completeness criteria and to Mike Harper, Belle Hudischewskyj, and Gary 

Lundberg for assistance in computing. The authors also thank Bill Hunt, 

Neil Frank, and Tom Curran of the U.S. Environmental Protection Agency for 

their thoughtful comments. 

83146 1 m 
o 

o 

o 

P.6 

CONTENTS 

O 

O 

O 

O 

O 

O 

Acknowl edgments iii 

List of Illustrations ., vi 

List of Tables x 

Notation xiii 

1 INTRODUCTION 1 

Background 1 

Methodology and Organization of Results 3 

2 STATISTICAL PROCEDURES FOR DETERMINING 

THE AIR QUALITY STATUS OF A SITE 6 

Annual or Quarterly (Long-Term) Averages 6 

Short-Term Concentrations. 10 

The Observed Number of Short-Term Exceedances 10 

The Estimated Number of Short-Term Exceedances 15 

Current Practice in Determining the Short-Term 

Attainment Status of a Site 19 

Limiting the Maximum Probability of Type I Error 19 

Summary of Results and Alternative Approaches 27 

3 EMPIRICAL ESTIMATION OF THE PROBABILITY OF 

MISCLASSIFYING ATTAINMENT STATUS 31 

Long-Term Ambient Standards 31 

Short-Term Ambient Standards 43 

O L i (Time t I J f ooooooo»ooooooooooooooooo.>0eo«o»«.oo*eo**o«oo*oo«oo« * J w 

o 

83m6r - i v 
o 

P.7 

o 

o 

o 

o 

o 

STATISTICAL ESTIMATION OF THE STATUS AND 

TRENDS IN NATIONAL AIR QUALITY 

rt LJU cb i I Uil 0 1 DlG3o«o»ooa««ooaoooo«ooft«»ao o 

Variance of the Estimators 

AN EMPIRICAL ASSESSMENT OF THE DATA COMPLETENESS 

NEEDED TO RELIABLY ESTIMATE NATIONAL TRENDS.. 

Sulfur Dioxide. 

Carbon Monoxide 

Nitrogen Dioxid 

UZUMC t « » « o « •« • • 

L6d Q OOOOOOOOO O 

Total Suspended 

Summa ry 

Parti 

oooooooooo o 

culate. 

00000*00000 0 

o 

KcTcrCMCc J oooooooo o 

Appendix A SULFUR DIOXIDE DESCRIPTIVE STATISTICS 

Appendix B CARBON MONOXIDE DESCRIPTIVE STATISTICS 

Appendix C NITROGEN DIOXIDE DESCRIPTIVE STATISTICS 

Appendix D OZONE DESCRIPTIVE STATISTICS 

Appendix E LEAD DESCRIPTIVE STATISTICS 

Appendix F PARTICULATE MATTER DESCRIPTIVE STATISTICS 

Appendix G SUMMARY OF CURRENT DATA COMPLETENESS CRITERIA 

55 

55 

56 

58 

58 

58 

65 

65 

72 

72 

72 

80 

O 

O 

O 

O 

831H6r 1 
P . 8 

O 

o 

ILLUSTRATIONS 

O 

2-1 Binomial approximation to a and 3 in the one-exceedance 

case, r = l 21 

Q 2-2 Binomial approximation to a and 3 in the three-exceedance 

case, r = 3 23 

2-3 Binomial approximation to a and 3 for r = 1 and a 

specification of CXQ = 0.2 26 

u 2-4 Binomial approximation to a and 3 for r = 3 and a 

specification of ag = 0.2 28 

3-1 Lead: the probability of misjudging actual attainment status 

as determined by the mean, the coefficient of variation, 

O and the proportion of data completeness 36 

3-2 TSP, annual geometric mean: the probability of misjudging 

actual attainment status as determined by the natural log 

of the geometric mean, the coefficient of variation, and 

O the proportion of data completeness 37 

3-3 TSP, annual arithmetic mean: the probability of misjudging 

actual attainment status as determined by the mean, the 

coefficient of variation, and the proportion of data 

_. COffl U i ClCllcb b oooooooooooooooooooooooooooooooooooooooooooooooo o O O 

A A 

3-4 Estimated probabilities a and 3 of misjudging attainment 

of the 24-hour sulfur dioxide standard as functions of the 

proporti on of data compl eteness 46 

A A 

O 3-5 Estimated probabilities a and 3 of misjudging attainment 

of the one-hour carbon monoxide standard as functions of 

the proportion of data completeness... 47 

A A 

3-6 Estimated probabilities a and 3 of misjudging the 8-hour 

O carbon monoxide standard as functions of the proportion of 

GuU u COITi p 1 GLGUGS S ooooooooooooooooooooooooooooooooooooooooooo o H O 

83m6r 1 vi 

O 
P.9 

O 

3-7 Estimated probabilities a and 3 of misjudging attainment 

O of the maximum daily one-hour ozone standard as functions 

of the proportion of data completeness 49 

A A 

3-8 Estimated probabilities a and 3 of misjudging attainment 

of the 24-hour TSP standard as functions of the proportion 

\^S O T Qdl u COITi p I CtCllCb b oooooooooooooeoooooooooooooooooooooooooo o *) U 

5-1 Estimated long-term sulfur dioxide air quality parameters 

as functions of the required data completeness.. 60 

Q 5-2 Estimated short-term sulfur dioxide air quality parameters 

as functions of the required data completeness 62 

5-3 Estimated one-hour carbon monoxide air quality parameters 

as functions of the required data completeness 64 

O 

O 

O 

O 

5-4 Estimated 8-hour carbon monoxide air quality parameters as 

functions of the required data completeness 67 

5-5 Estimated long-term nitrogen dioxide air quality parameters 

as functions of the required data completeness 69 

5-6 Estimated annual ozone air quality parameters as functions 

of the required data completeness 71 

5-7 Estimated long-term lead air quality parameters as functions 

O of the required data completeness 74 

5-8 Estimated long-term TSP air quality parameters as functions 

of the required data completeness.. 76 

5-9 Estimated short-term TSP air quality parameters as functions 

of the required data completeness 78 

Appendix A: Sulfur Dioxide 

A-l Histogram of the percentage of data completeness for 

O sites reporting sulfur dioxide concentrations in 1981 

A-2 Scatter plots of ambient air quality statistics versus the 

percentage of total hours recorded in 1981 by sites 

reporting sulfur dioxide concentrations 

A-3 Histogram of expected primary exceedances for sites 

reporting 24-hour sulfur dioxide concentrations in 1981 

0 eaiwr 1 vii 
P.10 

O 

A-4 Scatter plot of expected primary exceedances versus the 

percentage of total days recorded in 1981 by sites 

O reporting 24-hour sulfur dioxide concentrations 

Appendix B: Carbon Monoxide 

B-l Histogram of the percentage of data completeness for 

O sites reporting carbon monoxide concentrations in 1981 

B-2 Histogram of expected primary exceedances for sites 

reporting one-hour carbon monoxide concentrations in 1981 

Q B-3 Scatter plot of the number of expected primary exceedances 

versus the percentage of total hours recorded in 1981 by 

sites reporting one-hour carbon monoxide concentrations 

O 

B-4 Histogram of expected primary exceedances for sites 

reporting 8-hour carbon monoxide concentrations in 1981 

B-5 Scatter plot of the number of expected primary exceedances 

versus the percentage of total hours recorded in 1981 by 

sites reporting 8-hour carbon monoxide concentrations 

O Appendix C: Nitrogen Dioxide 

C-1 Histogram of the percentage of data completeness for 

sites reporting nitrogen dioxide concentrations in 1981 

Q C-2 Scatter plots of ambient air quality statistics versus the 

percentage of total hours recorded in 1981 by sites 

reporting nitrogen dioxide concentrations in 1981 

O 

O 

O 

Appendix D: Ozone 

D-1 Histogram of the percentage of data completeness for 

sites reporting ozone concentrations in 1981 

D-2 Histogram of the number of expected primary exceedances 

for sites reporting ozone concentrations in 1981 

D-3 Scatter plot of the number of expected primary exceedances 

versus the percentage of data completeness for sites 

reporting ozone concentrations in 1981 

83m6r i vm 

O 
P.11 

Appendix E: Lead 

E-1 Histogram of the percentage of data completeness for 

sites reporting lead concentrations in 1981 

E-2 Scatter plots of ambient air quality statistics versus the 

percentage of data completeness for sites reporting lead 

concentrations in 1981 

Appendix F: Total Suspended Particulate 

F-1 Histogram of the percentage of data completeness for 

sites reporting TSP concentrations in 1981 

F-2 Scatter plots of ambient air quality statistics versus the 

percentage of data completeness for sites reporting TSP 

concentrations in 1981 

F-3 Histogram of the number of expected primary exceedances 

for sites reporting TSP concentrations in 1981 

F-4 Scatter plot of the number of expected primary exceedances 

versus the percentage of data completeness for sites 

reporting TSP concentrations in 1981 

83146r - ix 
0 

o 

o 

TABLES 

p.i; 

O 

O 

O 

o 

o 

o 

2-1 The required proportion of data completeness as a function 

of the total number of possible measurements, the coefficient 

of variation, and the prescribed maximum relative error 

occurring with 95 percent probability 8 

2-2 Relative error for data completeness proportions p^ = 0.50 

and Pi = 0.75 9 

2-3 The proportion of data completeness required to limit the 

risk to 3 of not detecting r + 1 exceedances 14 

2-4 The maximum probability of a Type I error and the binomial 

approximation to the probability, as functions of the proportion of data completeness.... 18 

2-5 Binomial approximations to probabilities a and 3 in the 

O i l G ""GXCGGQ O MC C vQbt . ) P "~ Xoo000000«09000000000000«000*0««*>00 0 C U 

2-6 Binomial approximations to probabilities a and 3 in the 

three-exceedance case, r = 3. 22 

2-7 Attainment status as a function of observed exceedances and 

the proportion of data completeness 25 

2-8 Attainment status for the one-exceedance form of the 

SX» a MQc i l U \ \ *"*" iyooooo«©o»oooooo©ooo*oo©ooo©ooooooo»>oooooos»s»o o CO 

3-1 Estimated probability that a site' s attainment status for 

long-term standards i s misclassified as a function of 

minimum data completeness proportion 34 

3-2 Estimated probability that a site' s attainment status for 

short-term standards i s misclassified as a function of 

minimum data completeness proportion... . 51 

O 

O 

831 ^ 1 
<•> 

P.13 

O 

o 

o 

o 

o 

5-1 

5-2 

5-3 

5-4 

5-5 

5-6 

5-7 

Estimated composite annual average of one-hour sulfur dioxide 

and its estimated variance as functions of the required 

percentage of data compl eteness 

Estimated mean annual rate of 24-hour sulfur dioxide 

exceedance and its estimated variance as functions of the 

required percentage of data completeness 

Estimated mean annual rate of one-hour carbon monoxide 

exceedance and its estimated variance as functions of the 

required percentage of data completeness 

Estimated mean annual rate of 8-hour carbon monoxide 

exceedance and its estimated variance as functions of the 

required percentage of data completeness 

Estimated composite annual average of one-hour nitrogen 

dioxide and its estimated variance as functions of the 

required percentage of data completeness 

Estimated mean annual rate of maximum daily one-hour ozone 

exceedance and its estimated variance as functions of the 

required percentage of data completeness 

Estimated composite average of maximum quarterly lead and 

its estimated variance as functions of the required 

percentage of data completeness 

59 

61 

63 

66 

68 

70 

73 

O 

O 

5-8 

5-9 

Estimated composite average of 24-hour TSP and its estimated 

variance as functions of the required percentage of data 

COIli p I GUGllG j J OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO i 

Estimated mean annual rate of 24-hour TSP exceedance and its 

estimated variance as functions of the required percentage 

of data compl eteness 

75 

77 

O 

O 

Appendix A: Sulfur Dioxide 

A-l 

A-2 

Number of sites recording one-hour sulfur dioxide during at 

least 75 percent of the total possible hours in 1981 

Number of sites recording one-hour sulfur dioxide during at 

least 50 percent of the total possible hours in 1981 

83146r 1 x i 

O 
P.M 

o 

A-3 Frequency table of the percentage of total days recorded in 

1981 versus expected primary exceedances for sites reporting 

^ 24-hour sulfur dioxide averages 

Appendix B: Carbon Monoxide 

B-l Number of sites recording carbon monoxide during at least 

O 75 percent of the total possible hours in 1981 

B-2 Number of sites recording carbon monoxide during at least 

50 percent of the total possible hours in 1981 

O B-3 Frequency table of the percentage of total days recorded in 

1981 versus expected primary exceedances for sites reporting 

one-hour carbon monoxide averages 

O 

B-4 Frequency table of the percentage of data completeness in 

1981 versus expected primary exceedances for sites reporting 

8-hour carbon monoxide averages 

Appendix D: Ozone 

D-1 Frequency table of the percentage of data completeness in 

0 1981 versus expected primary exceedances for sites reporting 

maximum daily ozone concentrations 

Appendix F: Total Suspended Particulate 

O F-1 Number of sites recording 24-hour TSP during at least 75 

percent of the scheduled monitoring days in 1981 

F-2 Number of sites recording 24-hour TSP during at least 50 

percent of the scheduled monitoring days in 1981 

F-3 Frequency table of the percentage of total days recorded in 

1981 versus expected primary exceedances for sites reporting 

24-hour TSP averages 

O 

O 

Appendix G: Current Data Completeness Criteria 

G-l National Aerometric Data Bank's criteri a for validit y 

of data 

G-2 California Air Resources Board's criteri a for representaO tiveness of air quality measurements and statistic s 

G-3 Indiana Air Pollution Control Division's criteri a for 

completeness of data 

83146r

 1 x1 1 
P.15 

NOTATION 

ft 

Fixed or controlled parameters 

N = number of possible measurements in a year 

PQ = required data completeness proportion or percentage 

e = specified tolerance for estimating a site-specific long-term 

average concentration 

'• r = the maximum number of times short-term concentrations at a 

site may exceed a reference concentration without violating 

a given short-term air quality standard 

a = probability that the attainment status of a site will be 

> incorrectly determined as nonattainment 

3 = probability that a site's failure to attain an air quality 

standard will go undetected 

A Site-specific quantities (i = site index) 

pn- = proportion of possible measurements actually recorded during 

the year 

A 

Ui = average of the p^N concentrations recorded during the year 

^i = average of all N concentrations that could have been 

recorded during the year 

o. = standard deviation of the p^N concentrations recorded during 

• the year 

Oi = standard deviation of all N concentrations that could have 

been recorded during the year 

0, ki = the number of times a site is reported to have exceeded a 

short-term reference concentration 

831<+6 7 xiii 
P.16 

X. = the estimated number of short-term exceedances 

A. = the actual number of short-term exceedances (equal to ki 

when Pi = 1) 

Sample and population quantities 

/Jr = set of indices i of sites comprising the sample 

V(PQ ) = number of siites in the sample whose data completeness 

proportion equals or exceeds pg 

A A 

U(PQ ) = sample average of Pi across sites with p. > pQ 

P(PQ ) = population mean of pi across sites with p. > PQ 

P = population mean of pi 

"2 "2 

a (pn) = sample average of variance estimates o* across sites 

with p. > pQ 

2 2 

a (pQ) = population mean of variances Oi across sites with p. *» p^ 

2 2 

a = population mean of variances a^ 

A A 

x(pQ) = sample average of estimated site-specific exceedances (A..) 

for those sites with p. > pQ 

A(P 0 ) = population mean of site-specific exceedances (Ai) for those 

sites with p. > pQ 

A = population mean of site-specific exceedances 

8311+6 7 X1 v 
P.17 

EXECUTIVE SUMMARY 

PURPOSE AND SCOPE OF THE STUDY 

The U.S. Environmental Protection Agency (EPA) stores ambient air quality 

data, obtained from a nationwide network of monitoring stations, in the 

National Aerometric Data Bank (NADB). Statistics summarizing the air 

quality at each monitoring station are also stored in the NADB. These 

summary statistics are made available to the government-and public 

provided that certain data validity criteria are met, chief among these 

being requirements for data completeness (given in Appendix G) . 

This data screening policy helps to ensure that NADB summary statistics 

are reliable. Moreover, the policy of using one set of data validity 

criteria regardless of the intended use of the data is necessary since 

the EPA must respond quickly to numerous requests for NADB data. It is 

possible, however, that some information may be unnecessarily lost to the 

general NADB user because of the stringency of the criteria for data 

validity. Nevertheless, stringent criteria may be needed for some 

purposes, notably for determining whether National Ambient Air Quality 

Standards (NAAQS) have been attained. 

The purpose of this report is to examine data reliability as a function of 

data completeness in order to assist the EPA in its review of the adequacy 

and appropriateness of the current NADB data validity criteria. 

We consider two important cases at opposite ends of the spectrum of data 

usage: (1) the determination of NAAQS attainment at a site, which is 

essential to the implementation of the Clean Air Act and which depends 

heavily on complete data, and (2) the estimation of the status of, and 

831t6r 1 xv 
P.18 

trends in, air quality nationwide, for which much less complete data can 

be tolerated since data are pooled from a large number of sites. 

METHODOLOGY 

We take the view in this report that summary statistics based on incomplete data are estimates of the summary statistics that would have been 

obtained from complete data. Moreover, we treat the set of observations 

at a site as a random sample from the entire set of pollutant concentrations, even in the case of regular but intermittent sampling (so that a 

sampling schedule of once every six days produces at most 17 percent 

completeness). 

From this viewpoint, the correct NAAQS attainment status of a site depends 

solely on the complete set of measured pollutant, concentrations, so that 

misclassification is possible only for those sites having incomplete 

data. Similarly, the accuracy of estimated trends in' national air quality 

depends on the data completeness at each site.as well as the set of sites 

used to form the estimates. We have thus set aside other important 

components of data reliability, i.e., measurement error (due to a variety 

of inherent limitations, such as instrument calibration and maintenance or 

recording and reporting practices) or the degree of spatial representativeness of measured concentrations. 

We adopt two forms of analysis: (1) we analyze hypothetical "mostsensitive" scenarios (in which sites barely attain or fail to attain 

short-term air quality standards, or where the hypothetical high variability of pollutant concentrations makes the estimation of long-term 

averages more difficult) and (2) we estimate actual error rates based on 

1981 NADB data (described in Appendixes A through F) . 

The criterion we use to evaluate the effect of data completeness requirements on the estimation of national air quality trends is the variance of 

the estimated national average in a single year. We show that strict data 

831i+6r 1 

xvi 
o 

P.19 

o 

o 

completeness requirements affec t the variance in two distinc t ways. That 

i s , excluding sites with less complete data has two effects : (1) an 

increase in the reliabilit y of those site-specifi c estimates included in 

the national estimate and (2) a reduction in the number of sites included 

in the national estimate. Using 1981 NADB data we evaluate these competing effects on the variance of the estimated national average. 

O 

O 

O 

Q 

O 

G 

MAJOR FINDINGS 

In general, reliable determination of the short-term NAAQS attainment 

status of a site requires more complete data than the other data uses we 

consider. However, as several investigators have noted, increasing the 

data completeness proportion (up to a point) also increases the chance of 

misclassifying a site that, in fact, has one exceedance. This peculiarit y 

is due to the fact that the more complete the data, the greater the chance 

of including the exceedance in the sample and thus deriving an expected 

exceedance rate greater than unity when correcting for the incompleteness 

of the data record. This scenario may be guarded against by requiring 

essentially complete data (more than 95 percent completeness), or by 

modifying the procedure for determining NAAQS attainment status. 

We used 1981 data in two ways to analyze data completeness needs: (1) we 

estimated the expected rate of misclassification for sites that met 

current and alternative completeness requirements in 1981, and (2) for 

long-term lead and particulate standards we estimated "critica l intervals " 

about the standards within which misclassification probabilities exceed 

10 percent, and we show how these intervals narrow as more complete data 

are required. (Critical intervals are defined in Section 3, and should 

not be confused with confidence intervals.) 

O 

Each form of analysis serves a different purpose. The firs t indicates 

rates of misclassification under current completeness requirements but is 

less useful for estimating the consequences of tightening or relaxing 

these requirements. The second form of analysis uses characteristics of 

O 

831<+6r 1 x v n 
G 

P.20 

G 

G 

the 1981 data to estimate the consequences of various data completeness 

requirements on the classification of lead or particulate monitoring sites 

whose long-term averages are near the standard. The second form thus 

focuses on those sites whose long-term average concentrations are near the 

standard without addressing the question of how many such sites there are 

likely to be (which cannot be done reliably in the absence of much more 

data or strong assumptions of uncertain validity). 

G 

G 

O 

The probability of misclassifying an attainment site under current data 

completeness requirements is 1 percent for the 1-hour CO standard, 2 

percent for the 8-hour CO, 24-hour S02, and 24-hour TSP standards, and 

8 percent for the daily maximum ozone standard (see Table 3-2). Whenever 

reliable estimation was possible, we estimated the probability of misclassifying a nonattainment site . For the 24-hour SO2 standard and the 

1-hour CO standard, reliable estimation was not possible due to the small 

number of nonattainment sites (less than 10). The probability is less 

than 1 percent for the 8-hour CO standard and the daily maximum ozone 

standard, and 6 percent for the 24-hour TSP standard. 

O 

O 

o 

We attempted to use the first form of analysis to estimate the effects of 

alternative date completeness requirements on the probability of misjudging the attainment of short-term standards. This was done by restricting 

the analysis to those sites meeting the alternative requirements. For 

continuously monitored pollutants, a relaxation of the 75 percent requirement to 50 percent produced no appreciable difference in the probability 

of misclassification. It is theoretically possible, but impractical, to 

reduce the probability of a false determination of nonattainment to zero 

by requiring 95 percent complete data. For the current 24-hour TSP standard, we examined sites meeting completeness requirements of 5, 12, and 

25 percent, i.e., requirements of 5, 11, and 23 valid samples per quarter, 

respectively. Again, the analysis was insensitive to differences in data 

completeness, although the estimate for 5 and 25 percent completeness 

differed slightly. 

G 

O 

83m6r -

xvm 
P.21 

G 

For long-term standards, the 1981 data indicate that the probability of 

Q misclassifying the status of a site that has actually attained a standard 

and that meets current NADB completeness requirements is 2 percent for TSP 

and less than 1 percent for other pollutants (see Table 3-1). The chance 

of misclassifying the status of a TSP monitoring site that has not 

Q attained the annual standard is estimated to be 9 percent. The estimated 

chance of misclassification is greater for TSP sites that have not 

attained the standard than for those that have attained the standard 

because the nonattainment sites tend to be only slightl y above the 

O standard, whereas the attainment sites tend to be well below the standard. For other long-term standards, there are too few nonattainment 

sites to reliably estimate the rate of misclassifying such sites . These 

results suggest that misciassifications are unlikely under current data 

G completeness requirements. This is true provided that the 1981 data are 

representative of other years, i.e. , so that actual long-term averages are 

well below the standard for al l but a small fraction of sites . 

G As in the case of short-term standards, the misclassification rate was 

approximately the same for sites meeting various data completeness 

requirements. The exception is that the estimated rate of misclassification of TSP sites that have not attained the standard drops from 9 

O percent to 3 percent when 25 percent completeness is required. 

Thus, direc t estimation of expected rates of misclassification generally 

requires more data or stronger assumptions (of uncertain validity ) to 

G distinguish the effects of alternative data completeness requirements. 

Using the second form of analysis, in which we estimated "critica l 

intervals " about a given standard, we analyzed the effec t of differen t 

O completeness requirements on the accuracy of classifications of attainment 

of the quarterly lead standard and the annual TSP standard. We gave 

special attention to these two pollutants since they are typicall y sampled 

only once every six days. For each pollutant we estimated the expected 

G total number of misclassifications (that is , misclassifications of both 

83U6r i xix 

G 
o 

o 

o 

o 

o 

attainment and nonattainment) to be less than 5 percent; the estimated 

percentage for lead was less than that for TSP (see Figures 3-1 through 

3-3). For lead (quarterly standard =1.5 pg/nr), the estimated critical 

interval when 12 samples per quarter are required is -17 percent to +27 

percent, i.e., 17 percent below to 27 percent above the standard. For TSP 

(annual standard = 75 pg/m3), the critical interval when five samples per 

quarter are required is -11 percent to +13 percent. If 12 TSP samples per 

quarter were required, the interval would narrow to -7 percent to +8 

percent. 

If our goal were to contain these critical intervals within ± 10 percent 

of the respective standards, then strengthened requirements should be 

considered for lead and possibly for particulates. Indeed, judging from 

our analysis of TSP, this goal would be achieved under a recently proposed 

revision to standards for particulate matter that includes a requirement 

of 12 valid observations per quarter (EPA, 1984; 49 FR 10434). 

O 

O 

A similar critica l interval analysis of short-term standards is possible, 

although thi s would require additional assumptions about the frequency 

distribution of concentrations at individual sites . In the interests of 

time and project resources we did not conduct thi s analysis, although i t 

may be useful to do so in the future. 

O 

O 

O 

For reliable estimation of national trends, the benefits of relaxing data 

completeness requirements, i.e., broadening the base of sites from which 

to form the national average over successive years, outweigh the disadvantages of including the sites that have less complete data and 

therefore less reliable summary statistics (see Section 5). This is 

because the variability in air quality among sites outweighs the sampling 

variability at individual sites. Although no optimum requirement can be 

identified from the analysis, the current 50 percent data completeness 

requirement used for the annual report on national air quality trends 

could be relaxed without increasing calculated variances. This conclusion 

is based on the distribution of the percentages of data completeness found 

O 

831"+6r 1 XX 
in the NADB records for 1981, and assumes that the distribution would not 

change substantially in future years if the current requirement for 

estimating national trends were relaxed. 

RELATED ISSUES 

We have attempted to be thorough within the scope of this study, however 

we must point out that data completeness requirements need to be evaluated 

in a broader context that includes 

(1) overall data quality; 

(2) the procedures used to determine the attainment status of a 

site, and 

(3) the ambient standards themselves. 

Some alternatives to the current procedure for determining short-term 

attainment status use a greater number of the high concentrations recorded 

at a site (see Section 2). Although these alternative procedures are 

somewhat more complicated, they appear to be less sensitive to incomplete 

data and avoid the counter-intuitive result, noted above, that the chance 

of misclassification increases with increasing data completeness. 

Although this report focuses on data completeness, there are other aspects 

of data reliability, e.g., quality assurance procedures. It is important 

to remember that requiring more complete data, e.g., requiring more 

frequent intermittent sampling, is not the only way to increase data 

reliability, and not necessarily the most cost-effective way. Thus all 

aspects must be considered together if enhanced data reliability is the 

objective. 

Similarly, the question of data reliability belongs to a larger cirele of 

issues concerning effective air quality management. If increased data 

831"46r i XXI 
reliability reduces the chance of misclassification, then the marginal 

costs of improved data reliability should be weighed against the economic 

and public health benefits of fewer misclassifications. If we reduce the 

chance that an attainment site will be misclassified as nonattainment, 

then unnecessary administrative and control costs will also be reduced. 

Likewise, if we reduce the chance that a nonattainment site will be 

misclassified as attainment, then presumably public health is better 

maintained. In evaluating the latter benefit, however, one should bear in 

mind that a margin of safety is incorporated into ambient standards so 

that actual pollutant concentrations at a nonattainment site may yet be 

below levels causing adverse health effects. Moreover, it is possible 

that with sufficient effort devoted to data reliability, our certainty 

about actual air quality could exceed our certainty about air quality 

heaith effects. 

In short, monitoring requirements (data completeness, quality assurance, 

numbers of monitors, frequency of sampling, etc.) and other regulatory 

issues, such as the margin of safety embodied in ambient standards, are 

interdependent. For example, the smaller the margin of safety used in 

deriving an air quality standard, the less we can afford to misjudge the 

attainment status of a site, and the greater the need for highly reliable 

data. 

RECOMMENDATIONS 

One use of the NADB data that we consider is the estimation of the status 

of and trends in air quality nationwide. Tightening the data requirements 

has two competing effects: (1) the sites meeting the stricter requirements (and therefore included in the national average) have more reliable 

air quality summary statistics, but (2) the exclusion of sites with less 

reliable statistics reduces the sample size used to form the national 

average. Our conclusion is that the current 50 percent data completeness 

requirement for determining national trends, which is a relaxation of the 

75 percent NADB requirement, should be further relaxed because the 

83mer l xxii 
o 

o 

advantage of a large sample size outweighs the advantage of reliable 

statistics at individual sites. Unfortunately, the optimum requirement 

cannot be identified from our analysis, so that any specific recommendation will require further consideration. 

G 

O 

G 

O 

© 

G> 

G 

G 

The EPA staff responsible for the estimation and reporting of national air 

quality trends have direct access to NADB data. This is an example of a 

special group of NADB users who can apply data completeness requirements 

according to their particular needs. In considering the release of NADB 

data to the general user, however, it appears wise to continue the policy 

of using one set of data completeness requirements. Special users with 

direct access to NADB data, such as those responsible for the estimation 

of national air quality trends, should continue to determine the degree of 

data completeness needed for their particular application. 

Because of the practicality of using one set of requirements, the requirements should be based on the common application having the greatest need 

for complete data. On this basis, we have focused on the issue of whether 

current NADB criteria provide sufficient protection against misclassifying 

the attainment status of a site with respect to short-term or long-term 

standards. 

Certainly, accurate estimates of national trends require less complete 

data than do accurate classifications of individual sites (compare Figures 

5-1 through 5-9 with Tables 3-1 and 3-2). Of course, we cannot compare 

these data needs quantitatively until we prescribe the accuracies acceptable in each case, but the evidence supporting the statement is clear 

nevertheless. It is also clear that there are other possible data uses, 

not examined in this report, that depend even more heavily on complete 

data than does the determination of attainment status, e.g., inferring the 

maximum 10-minute concentration from hourly data. But, for the purpose 

of establishing appropriate NADB requirements for complete data, the 

determination of attainment status seems to be the most exacting use of 

air quality data that is common to all the pollutants for which we have 

national ambient air quality standards. 

83U6T 1 XXI11 

G 
o 

o 

o 

Direct estimates of the expected misclassification rate for continuously 

monitored pollutants do not change appreciably when the current NADB 

requirement of 75 percent data completeness is relaxed to 50 percent. 

Analyses of the most sensitive scenarios (which could become more common 

than they were in 1981) warn against this relaxation, however. In 

particular, if the procedures for determining the attainment of short-term 

standards remain unchanged, it seems best for the time being to also 

maintain current data completeness requirements for continuously monitored 

pollutants. 

O 

O 

O 

Similarly, in the case of intermittent sampling, e.g., for lead and 

particulate matter, direct estimates of misclassification rates barely 

distinguish between the alternative data completeness requirements we 

consider. In general, the NADB data validity criteria require only five 

samples per quarter, and thus only 20 samples per year. For lead, 

however, 12 samples per quarter are required for attainment determinations 

(EPA, 1979; 43 FR 46266). In either case, it is much more difficult to 

ensure the correct classification of air quality status under an intermittent sampling. Critical interval analyses (Figure 3-1 through 3-3) 

demonstrate the advantages of strengthening data completeness requirements 

for lead and particulates. It is therefore advisable to adopt the 

recently proposed requirement of 12 particulate samples per quarter (EPA, 

1984; 49 FR 10434) and to consider even more frequent sampling of lead. 

O Concerning the issues related to data completeness--less central to the 

scope of this study but no less important for effective air quality 

management--we recommend consideration of two changes in procedures. 

O (1) A greater number of the high concentrations recorded at a site 

should be used in order to provide more robust short-term 

attainment inferences. 

O 

O 

83146T

 l 

XXI v 
P.27 

(2) Data completeness considerations should be incorporated in the 

review, revision, and interpretation of air quality standards, 

along with such related issues as the cost of more complete 

data, other limitations on data reliability, the cost of 

misclassifying attainment status, and the existing margins of 

safety used to set thresholds for ambient concentrations. 

831t6r 1 

xxv 
INTRODUCTION 

The U.S. Environmental Protection Agency (EPA) stores ambient air quality 

data, obtained from a nationwide network of monitoring stations, in the 

National Aerometric Data Bank (NADB). Statistics summarizing the air 

quality at each monitoring station are also stored in the NADB. These 

summary statistics are made available to the government and public 

provided that certain data validity criteria are met, chief among these 

being requirements for data completeness (given in Appendix G) . 

This data screening policy helps to ensure that NADB summary statistics 

are reliable. Moreover, the policy of using one set of data validity 

criteria regardless of the intended use of the data is necessary since the 

EPA must respond quickly to numerous requests for NADB data. It is 

possible, however, that some information may be unnecessarily lost to the 

general NADB user because of the stringency of the data validity 

criteria. Nevertheless, a stringent criteria may be needed for some 

purposes, notably, for determining whether National Ambient Air Quality 

Standards (NAAQS) have been attained. 

The purpose of this report is to examine data reliability as a function of 

data completeness in order to assist the EPA in its review of the adequacy 

and appropriateness of the current NADB data validity criteria. 

BACKGROUND 

Air quality data are provided to the National Aerometric Data Bank (NADB) 

primarily by state and local agencies. Included in the NADB record of air 

quality summary statistics for a monitoring station is a "criteria flag" 

83146r 2 
indicating whether NADB data validity criteria are met. These (and other) 

data validity criteria are presented in Appendix G, but they may be 

roughly summarized as a 75 percent data completeness requirement. (Note 

that the monitoring schedule for lead and particulate matter is once every 

six days. In these cases of intermittent sampling of 24-hour concentrations, NADB criteria require five samples per quarter and thus 20 samples 

per year. 

In responding to typical requests for NADB data, the record of a site 

failing to meet the criteria may be written with question marks to 

indicate questionable data reliability, or the record may not be written 

at all. Thus, whereas authorized users of NADB data have access to all 

records, the records indicating that NADB criteria have been met constitute a common subset of NADB data that is reliable and widely available. 

In this report we consider the effect of excluding air quality data from 

sites that fail to meet certain data completeness requirements. More 

specifically, we consider whether current NADB requirements appear to be 

satisfactory or need to be changed, and, in the latter case, whether data 

completeness requirements should vary with data usage and type of monitoring (intermittent or continuous). 

We consider two important cases at opposite ends of the spectrum of data 

usage: (1) the determination of NAAQS attainment at a site, which is 

essential to the implementation of the Clean Air Act and which makes a 

heavy demand for data completeness, and (2) the estimation of the status 

of and trends in air quality nationwide, for which much less complete data 

can be tolerated (since data are pooled from a large number of sites). 

For each of these two data uses we consider summary statistics relevant to 

short-term and long-term air quality standards. That is, we consider the 

number of exceedances of the short-term standards (having an averaging 

time of 24 hours or'less), and the annual or maximum quarterly average. 

Again, by treating both peak value summary statistics (where the effect of 

831<+6 2 
incomplete data is pronounced) and long-term averages (where data completeness needs are minimal), we hope to provide useful benchmarks for 

other data uses. 

Currently, the attainment of most national ambient air quality standards 

is determined separately for each year of data. However, a three-year 

attainment period is used in the case of ozone, and multi-year periods are 

being considered for other short-term standards. We therefore address the 

data completeness needed for both single-year and multi-year attainment 

standards. 

METHODOLOGY AND ORGANIZATION OF RESULTS 

We take the view in this report that summary statistics based on incomplete data are estimates of the summary statistics that would have been 

obtained from complete data. Moreoever, we treat the set of observations 

at a site as a random sample from the entire set of pollutant concentrations, even in the case of regular but intermittent sampling (so that a 

sampling schedule of once e^/ery six days produces at most 17 percent 

completeness). 

From this viewpoint., the correct NAAQS attainment status of a site depends 

solely on the complete set of measured pollutant concentrations, so that 

misclassification is possible only for those sites having incomplete 

data. Similarly, the accuracy of estimated trends in national air quality 

depends on the data completeness at each site as well as the sample of 

sites used to form the estimates. We have thus set aside other important 

components of data reliability, i.e., measurement error (due to a variety 

of inherent limitations, such as instrument calibration and maintenance or 

recording and reporting practices) or the degree of spatial representativeness of measured concentrations. 

We adopt two forms of analysis. In Section 2 we analyze hypothetical 

"most-sensitive" scenarios (in which sites barely attain or fail to attain 

83m6 2 3 
P.31 

O 

0 

o 

o 

short-term air quality standards, or in which the hypothetical high 

variability of pollutant concentrations makes the estimation of long-term 

averages more difficult). In Section 3 we assess the likely error rates 

based on 1981 NADB data (described in Appendixes A through F) . 

In our treatment of short-term NAAQS, the assigned status of a site is 

based on a quantity known as the "expected number of exceedances," that 

is, attainment status is determined assuming that the ratio of observed 

exceedances to total observations is maintained for the complete data set. 

The criterion we use to judge the degree of data completeness required to 

accurately determine NAAQS attainment is the chance of misclassifying the 

status of a site. The chance is zero for sites having 100 percent 

^ complete data or for sites having no exceedances. Otherwise the chance of 

misclassification depends on both the degree of data completeness and the 

actual number of exceedances. 

Q In Section 2 we consider the effect of data completeness on the chance of 

erroneously classifying the status of a site as nonattainment when in fact 

the reference concentration has been exceeded exactly once in the course 

of the year. We also consider the chance of erroneously classifying the 

Q status as attainment when the site has exactly two exceedances. This 

analysis is also carried out for air quality standards allowing no more 

than three exceedances in three years. 

Q In Section 3 we use 1981 NADB data to estimate the actual (but unknown) 

rates of misclassification. We first estimate the frequency distribution 

among sites of the actual number of exceedances. We then use this 

distribution to estimate the probability of misclassifying sites meeting 

Q various data completeness requirements. 

For long-term (annual or quarterly) ambient standards, the chance of 

misclassifying a site depends not only on the degree of data completeness 

O at the site but also on the empirical frequency distribution of the entire 

year of pollutant concentrations at the site. In Section 2 we assess data 

831<+6 2 

O 
completeness requirements by assuming that the actual average concentration at the site is 10 percent above or below the standard, and that the 

standard deviation of the year of concentrations at the site is twice the 

average. In Section 3 we empirically estimate probabilities of misclassification using 1981 data. 

In Section 4 we introduce a criterion for evaluating the effect of data 

completeness requirements on the estimation of national air quality 

trends, namely the variance of the estimated national average in a single 

year. Excluding sites with less complete data has two effects: (1) the 

reliability of the site-specific estimates included in the national 

estimate is increased, and (2) the number of sites included in the 

national estimate is reduced. In Section 5, using 1981 NADB data, we 

evaluate these competing effects on the variance of the estimated national 

average. 

8311+6 2 
STATISTICAL PROCEDURES FOR DETERMINING 

THE AIR QUALITY STATUS OF A SITE 

ANNUAL OR QUARTERLY (LONG-TERM) AVERAGES 

Suppose we wish to know (up to measurement error) the annual average 

concentration pi of a pollutant at a particular site, with index i, based 

on a complete set of N measurements, but that only ni valid measurements 

were recorded, the remaining mi = N - ni being either invalid or missing. 

Letting Pi denote the proportion, ^-/N, of data completeness, it can be 

A 

shown (Cochran, 1963) that the sample average p. has expected value p, and 

variance (1 - pi )/pi ,(o|/N): 

E u. 

Var 

Pi a. 

Moreover, annual air quality measurements usually entail large numbers of 

both missing and observed values (mi and nM) with relatively few extreme 

values. Under these conditions, formulated more precisely by Cochran 

A 

(1963), who cites Hajek (I960), the sampling distribution of Pi is 

approximately normal. 

2 Thus [(1 - p^)/p^Do./N is a natural measure of the error in estimating the 

annual average (pi) of all N measurements by the average of a proportion 

(Pi) of the N measurements. 

831H6r 3 
P.34 

We may use these results—the expected value, variance, and approximate 

A 

normality of Pi~to determine the minimum acceptable proportion of 

measurements, pg, provided that (1) we specify a relative error tolerance, 

e, such that errors of larger magnitude have, say, only a 5 percent of 

occurring, and (2) we have some idea of the coefficient of variation, 

aJ\i,. Specifically we have 

|Mi " Mil 

Pr l— ' — > E ) = Pr 

I"1 " ui ep. 

(Varp\) 1 / 2 (Var^) 1 '

2 

* 2 

EM. 

.(Var Ui ) 1/2 

$ e 

Wi (Np.,) 

i/2\n 

'i

 l ' > . 

where * is the standard normal cumulative distribution function. The 

above expression may be used to calculate either pg from a prescribed 

probability (e.g., 5 percent), or a confidence bound of a prescribed level 

A 

(e.g., 95 percent) on the absolute relative error (|p. - p.|/p.) 

given the proportion of data completeness (Pi) and the coefficient of 

variation (oi/p.) at a particular site. 

Table 2-1 shows the minimum proportion pg corresponding to several 

combinations of e, a^/\x^, and N. Table 2-2 shows E as a function of 

N, oJvi, and p.. Note that the 75 percent criterion used by the SAROAD 

system implies that partial annual averages are not likely to differ by 

more than appprbximately one percent from complete annual averages if the 

coefficient of variation equals unity (100 percent). The corresponding 

error for the 50 percent completeness criterion commonly used in reporting 

trends in annual averages is approximately two percent (for an individual 

site). 

83146 3 
TABLE 2-1. The required data completeness proportion (pg) as a function of 

the total number of possible measurements (N), the coefficient of variation (oi/pi), and the prescribed maximum relative error (E ) occurring 

with 95 percent probability. 

N = 8760 

(hrs/yr) 

N = 2920 

(3 hrs/yr) 

N = 1095 

(8 hrs/yr) 

N = 365 

(days/yr) 

N = 91 

(days/qtr) 

N = 30 

(days/mo) 

°i/yi 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

.001 

0.99096 

0.99772 

0.99943 

0.99697 

0.99924 

0.99981 

0.99886 

0.99972 

0.99993 

0.99962 

0.99991 

0.99998 

0.99991 

0.99998 

0.99999 

0.99997 

0.99999 

1.00000 

.01 

0.52297 

0.81431 

0.94607 

0.76684 

0.92936 

0.98135 

0.89765 

0.97229 

0.99292 

0.96339 

0.99059 

0.99763 

0.99061 

0.99764 

0.99941 

0.99689 

0.99922 

0.99980 

e 

.05 

0.04201 

0.14923 

0.41233 

0.11626 

0.34479 

0.67793 

0.25971 

0.58390 

0.84878 

0.51278 

0.80805 

0.94394 

0.80848 

0.94409 

0.98541 

0.92756 

0.98085 

0.99514 

.10 

0.01084 

0.04201 

0.14923 

0.03184 

0.11626 

0.34479 

0.08063 

0.25971 

0.58390 

0.20831 

0.51278 

0.80805 

0.51346 

0.80848 

0.94409 

0.76197 

0.92756 

0.98085 

1.00 

0.00011 

0.00044 

0.00175 

0.00033 

0.00131 

0.00523 

0.00088 

0.00350 

0.01384 

0.00262 

0.01041 

0.04040 

0.01044 

0.04050 

0.14446 

0.03102 

0.11351 

0.33871 

83146 

83132 
TABLE 2-2. Relative error (E) for data completeness 

proportions p. = 0.50 and p. = 0.75. 

P.36 

N = 8760 

(hrs/yr) 

N = 2920 

(3 hrs/yr) 

N = 1095 

(8 hrs/yr) 

N = 365 

(days/yr) 

N = 91 

(days/qtr) 

N = 30 

(days/mo) 

a • /u • r

 Mi 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

0.5 

1.0 

2.0 

.50 

0.01047 

0.02094 

0.04188 

0.01814 

0.03627 

0.07254 

0.02961 

0.05923 

0.11846 

0.05129 

0.10259 

0.20518 

0.10273 

0.20546 

0.41092 

0.17892 

0.35784 

0.71568 

.75 

0.00605 

0.01209 

0.02418 

0.01047 

0.02094 

0.04188 

0.01710 

0.03420 

0.06839 

0.02961 

0.05923 

0.11846 

0.05931 

0.11862 

0.23724 

0.10330 

0.20660 

0.41320 

8311+6 

8313 2 
P.37 

O 

O 

o 

In interpreting the minimum proportions given in Table 2-1, note that a 

given maximum relative error E may seem either excessively stringent or 

too lenient depending on the particular application. For example, 10 

percent relative error may be too large if we need to determine the 

attainment status of a site whose annual average is typically within 10 

percent of the annual average ambient standard whereas for a site whose 

annual average is typically less than half the standard , a 10 percent 

relative error specification may be unnecessarily stringent. 

O 

SHORT-TERM CONCENTRATIONS 

The Observed Number of Short-Term Exceedances 

O 

Current national ambient air quality standards (NAAQS) for pollutant 

concentrations averaged over 24 hours or less influence the presentation 

and analysis of high short-term concentrations. The short-term standards 

~ for most pollutants stipulate an ambient concentration that should not be 

exceeded more than once per year, which implies that the second highest 

short-term concentration recorded at a site determines whether the site is 

in attainment of the standard. 

O 

It is possible to treat the second highest value in a year's partial 

sample of short-term concentrations as an estimate of the second highest 

concentration in a complete sample, just as we treated long-term 

Q averages. This turns out to be a difficult approach from a statistical 

viewpoint, however, because of the great variability of the second highest 

concentration and the influence of the unknown tail distribution of shortterm concentrations on the sampling behavior of the second highest value. 

O 

Another, more workable approach, at least from a statistical perspective, 

is to examine the sequence of short-term concentrations at a site over a 

one-year period and to simply count the number of times the ambient 

Q concentration exceeds the standard. If the count is no more than one, the 

short-term standard is attained; otherwise it is not. 

8311+6 3 10 

O 
P.38 

O 

o 

G 

O 

O 

O 

O 

In this approach the sequence of short-term concentrations is converted to 

a 0-1 valued sequence, in which 1 indicates that the concentration exceeds 

the standard and 0 indicates otherwise. The.key question regarding the 

effect of data completeness on determining attainment status is whether an 

exceedance has been overlooked. More specifically, we wish to know the 

probability that a site has two or more exceedances but that no more than 

one exceedance is observed in our partial sample. We may formulate this 

as the following conditional probability: 

Pr (k. < 1 | xi = 2) 

where 

ki = the number of exceedances in the observed sample of size ni 

Xi = the unknown total number of exceedances that occurred during the 

year. 

The reason for focusing on the case X. = 2 is that this is the most 

challenging problem for a fixed partial sample size, ni; we would do well 

to prescribe the chance, say 3, we are willing to take that we fail to 

observe the two exceedances, since for larger values of Xi our chance of 

detecting a violation only improves. More generally, for a standard 

allowing no more than r exceedances we would want to calculate the 

Q probability that r or fewer exceedances are observed in the partial sample 

when there are r + 1 exceedances in the complete sample. 

This conditional probability is given by the hypergeometric distribution: 

O For future reference we note that 

E (ki | x., p.) = pi xi 

O 83m6 3 11 
P.39 

and 

N - Xi 

var (ki i x., P .) =-**pn: pi

( 1 "

 pi} xi 

* Pi U - Pi) Xn. 

(see, for example, Johnson and Kotz, 1969). Returning to the chance that 

not all• r + 1 exceedances are detected, we have 

or 

Noting that 

B*--Pr-(k.< r | x. = r + 1) 

1 - B = Pr (ki = r + 1 | Xi = r + 1) 

'N - r - 1\ //N 

ni ni - 1 nT ' r 

N N - 1 " • N - r 

i " j _ "T / N \ _J_ 

- j N \H - j / " N -

N v N - j ; N - j 

= Pi -TH-J^-PI ) 

83146 3 12 
P.40 

we see that 

i - e - Pl[P, - jj-i-j- d - P,)] ... [pf - ^ d - P,)] 

1

 Pi h -i^-p,-)] --- [PI -•&<! - PI>] 

_ r+1 /r + 1\ 1 „r ,, . . 

" Pi " I 2 JN

 Pi (1

 " Pi' * " • 

+ (-l) r

r!(l) rp j ( l - P i )

r 

- pr1 * 1 ^ + l ) (i)j prj+l a - Pi>j 

where the S^i are Stirling numbers of the first kind (Abramowitz and 

Stegun, 1964). 

Hence 

0 = 1 r+1 

= 1 for r = 1. 

Table 2-3 gives the data completeness proportion, p0, required to reduce 

the risk of failing to detect two exceedances to a specified level, 3. 

For the 75 percent completeness criterion, there is an approximate 44 

percent chance that at least one of the two exceedances will not be 

detected. For the 50 percent completeness criterion the risk is 75 

percent. It would seem then that the determination of short-term attainment status requires either nearly (e.g., 95 percent) complete data, or 

else some method of inference other than merely counting the number of 

8311+6 3 13 
P.41 

TABLE 2-3. The proportion, pg, of data completeness required 

to limit the risk to 3 of not detecting r+ 1 exceedances. 

a 

0.01 

0.02 

0.05 

0.10 

0.20 

0.44 

0.75 

r = 1 

0.995 

0.990 

0.975 

0.949 

0.894 

0.748 

0.500 

PO 

r = 3 

0.997 

0.995 

0.987 

0.974 

0.946 

0.865 

0.707 

8311+6 8 14 
P.42 

exceedances in the partial sample (see EPA, 1979; Pace and Frank, 1982; 

Thrall and Burton, 1983). Table 2-3 also shows the required data 

completeness for a three-exceedance standard, a topic to which we now 

turn. 

The Estimated Number of Short-Term Exceedances 

An alternative procedure for determining whether short-term air quality 

standards have been attained at a site is to extrapolate from the observed 

exceedances, ki, to the total number of exceedances, X-, i.e., to 

estimate Xi by 

^ - k1/p1 . 

From the conditional moments derived for ki, we have 

E (x\ | X., Pi ) = Xi 

.

 X " Pi 

Var (xi | X., p.) =~J— Xi 

i 

so that xi is an unbiased estimate whose variance approaches zero as Pi 

approaches unity. 

For an r-exceedance short-term standard, we might judge a site to be in 

attainment of the standard only if 

^ < r+ | . 

This procedure reduces the risk of failing to detect nonattainment 

(x. > r + 1), since the estimate Xi is an extrapolation of the observed 

number of exceedances (k^) among ^ observed concentrations to the total 

number of exceedances (x.) in N concentrations. But we have also introduced the risk that Xi might be overestimated. 

831<+6 3 15 
P.43 

This procedure for determining whether a short-term air quality standard 

has been attained at a site can be cast as a statistical test of a 

hypothesis. A site is presumed to be in attainment of the standard (i.e., 

1 the null hypothesis is that Xi < r) unless Xi *> r + -^ , in which case it 

is judged to have failed to attain the standard, i.e., we infer that the 

alternative hypotheis, X- *» r + 1, holds. 

A "Type I" error occurs when the site has actually attained the standard 

but we incorrectly infer that it has not. Previously, when we determined 

attainment solely on the basis of the observed number of exceedances, ki, 

this was impossible, since ki can never be greater than the actual number 

of exceedances, Xi. Now, however, the probability of a Type I error is a, 

where 

a = Pr (Xi > r + j | Xi < r j . 

Whatever the proportion of data completeness, Pi, the probability of a 

Type I error will be greatest when the site is just barely in attainment 

of the standard, i ,e., when x, = r, in which case 

P r \ V r + 1 I x i = r) * 

Now i f Pi > 2r/(2r + 1) i t is impossible for Xi to equal or exceed r + j 

since 

w 

<

 (2r)/(2r + 1) 

"

 r + 2 

so that 

a = 0 i f pi > 2r/(2r + 1) . 

83m6 3 16 
P.44 

1 

If 2(r - l)/(2r +1) < p. < 2r/(2r + 1) , Xi will equal or exceed r + j 

only if k- = r. In general, as Pi becomes smaller, fewer observed 

exceedances (k^) are required to produce an estimate (xH) equaling 

1 

or exceeding r + j , and the probability (a) of a Type I error increases, 

as shown in Table 2-4. Of course, when the sample size (ni) falls below 

r, the observed number of exceedances (ki) may be no more than ni and the 

risk of a Type I error is reduced. (We assume throughout this discussion 

that N > (r+1 ) (r + |) so that (r + 1)/N < 2/(2r + 1).) The binomial 

approximation given in Table 2-4 is accurate so long as N is large 

relative to r. 

Thus, basing an attainment determination on Xi rather than ki introduces 

some probability (a) of a Type I error, but the maximum probability (3) of 

a Type II error is substantially reduced for small values of Pi. We have 

= Pr (^i < r + | | xi = r + l ) 

• I (r: Tn;:. K ) 

f o r 2(J ~ 1) < p < 2j _ 

r o r

 2r + 1 p i * 2r + 1 

and for j = 1, 2, 3, ... , r+ 1 

t

1

 (r + 1 \ k n n xr+1-

S i l k ) pi

 (1

-

p

i

} 

Of special interest are the cases r = 1, the one-exceedance form of 

short-term standards, and r = 3, a three-year, three-exceedance form of 

831»+6 3 17 
oo 

• 

TABLE 2-4. The 

functions of the 

Pi Interval 

7TTT < Pi

 < l 

2r + 1 K i 

r

 - D 2 

N

 p i 2r + 1 

1 < „ "i < r 

N

 p i " N N 

# 

na ximum 

proport 

2j 

2r + 1 

• 

probabili 

ion pi of 

, 

, j = 2, 

t y , a 

data 

*3 j • • t 

• 

of a Type I 

completeness 

, r 

Pr 

r 

£ 

k=j 

r 

£ 

k=l 

n. 

l 

£ 

k=l 

• 

error < 

» 

(x. > r 

(W; 

(0(n N: 

(0(S( 

• 

ind the binomial 

4lV r> 

0 

-)/(", ) 

-)/(»i ) 

-V(-i) 

•-

approximation 

• 

to the 

, • 

probabil 

• 

ity , as 

Binomial Approximation 

k=j 

£ ( • 

k=l

 Vl 

0 

;)Pi< i 

:) A ** 

;) A <* 

-»i> r-

k 

-Pi) r"

k 

-Pi) r"

k 

83H+6 
the standards. Tables 2-5 and 2-6 show the binomial approximations to a 

and 3 for r = 1 and r = 3, respectively. The corresponding plots are 

shown in Figures 2-1 and 2-2. 

Current Practice in Determining the 

Short-Term Attainment Status of a Site 

In the previous.section we calculated a and 3, the probabilities of 

misclassifying attainment or nonattainment sites, respectively, if the 

estimated number of exceedances (x ^ is rounded to the nearest integer and 

then compared to the allowed number of exceedances (r). This procedure 

generated the regularly spaced data completeness cutpoints 2j/(2r + 1) or 

j/(r + 0.5) of Table 2-4. 

In fact, the current EPA policy is to round Xi to the nearest tenth and 

then compare the rounded value to r. That is, the status of a site is 

classified as nonattainment only if x., rounded to the nearest tenth, 

exceeds r. This procedure generates the data completeness cutpoints 

j/(r + .05) which are greater than the corresponding cutpoints j/(r + 0.5) 

of Table 2-4. 

For a site with data completeness proportion between j/(r + 0.5) and 

j/(r + .05), rounding to.the nearest tenth rather than to the nearest 

integer thus increases the chance of misclassifying the status of the site 

as nonattainment if in fact the total number of exceedances is r. On the 

other hand, rounding to the nearest tenth decreases the chance of misclassifying the status as attainment if the total number of exceedances is 

actually r+1 . 

Limiting the Maximum Probability of a Type I Error 

The effect of determining attainment status by comparing x- to the fixed 

1 ' 

value r + -p- is to create the regularly spaced data completeness cutpoints 

2j/(2r + 1) of Table 2-4. If Pi = 2r/(2r + 1), the maximum probability ct 

83146 3 19 
TABLE 2-5. Binomial approximations to 

probabilities a and 3 in the one-exceedance 

case, r = 1. 

Pi interval 

Binomial Approximation 

a 3 

2/3 < p. < 1 

0 < Pi < 2/3 

0 

Pi 

1-P i 

(1 - Pi) 1 

831i+6 8 20 
1.0 

0.6r 

0.6 

II 

sz 

CL 

0,4 

(a) a versus p. 

1.0 

0.3 

0.6 

0.4 

0.?-

o.o1 ' L 

1 . 0 

(b) 3 versus Pi 

8311+6 

FIGURE. 2-1. Binomial approximation to a and g in 

the one-exceedance case, r = 1. 

21 
TABLE 2-6. Binomial approximations to probabilities a and 3 in the three-exceedance case, r = 3. 

Pi Interval 

Binomial Approximation 

6/7 < P i 1 1 - P< 

P-0 

4/7 < p. < 6/7 6 Pi 2 (1 - Pi) 2 + 4 p. (1 - Pi )

3

 + (1 - Pi )

4 

2/7 < p. < 4/7 p. 3

 + 3 Pi

2

 (1 - Pi ) 4 Pi (1 - Pi) 3

 + (1 - Pi) 4 

3/N < p. < 2/7 1 - (1 - p.)" (1 - Pi) 

Pi = 2/N 3 Pi 2

 (1 - Pi) + 3 P i (1 - Pi )

2 (1 - Pi) 

Pi = 1/N 3 Pj (1 - Pj)' (1 - Pj) 

8 3 11+6 
1.0 

o.6r 

0.6 

0.4 

0.2 

1.0 

o.e 

0.6L 

1 3 

•*- > 

.o 

0. 4 

0.2 

(a) a versus p-j 

0.0 

0.0 

p, 

(b) 6 versus Pi 

FIGURE 2-2. Binomial approximations to a and 6 in 

the three-exceedance case, r = 3. 

8311+6 

23 
P.51 

of a Type I error will be between 0.60 and 0.67, depending on r. The 

point is that by using this method of determining the short-term attainment status of a site we run a high risk that the status might be erroneously determined as nonattainment. 

It is possible to modify the method so that the probability of a Type I 

error never exceeds a prescribed level, ag, at the cost of increasing the 

probability of a Type II error (in which a site's failure to attain the 

short-term standard goes undetected). The modification is to simply shift 

the data completeness cutpoints to smaller values, TT,, IT^, •••» it. which 

are obtained from ag. The procedure is depicted in Table 2-7. 

2-4, with 2j/(2r + 1) replaced by IT., for j 

J 

The resulting maximum probability of a Type I error is that given in Table 

1, 2, ..., r. The maximum 

probability, 3, of a Type II error is also obtained in this way from 

previous expressions. Below we demonstrate the derivation of TTJ from a 

prescribed level of risk, ag, for the cases r = 1 and r = 3. 

The procedure for the current one-exceedance form of a standard (r = 1) is 

depicted in Table 2-8. We see that a Type I error can only occur if ki = 

1, Pi < ITJ, and Xi = 1. The probability of a Type I error is thus zero if 

P i > ffl. If Pi < IT. then 

Pr(k. = l | Xi = 1) 

V 

so the absolute maximum (i.e., across all sites) occurs at Pi = TTJ. We 

may therefore limit the risk of a Type I error to ag simply by setting IT j 

= ag. Figure 2-3 depicts a and 3 as functions of Pi when the prescribed 

absolute maximum probability of a Type I error is ag = 0.20. 

In the case r - 3, ki will be taken at face value if Pi > v^* i.e., we 

will only declare the status to be nonattinament if we observe k^ > 4 

8311+6 3 24 
P.52 

TABLE 2-7. Attainment status as a 

function of observed exceedances (ki) and 

the proportion (pi) of data completeness. 

* 

0 

1 

1 

2 

2 

o < Pi < 1 

0 < Pi < TTj 

ITj < P i < 1 

0 < Pi < TT2 

TTo < PA < 1 

Assigned 

kn- Pi Interval Status 

A 

N 

A 

N 

A 

r 0 < pn- < irr N 

r irp < Pi < 1 A 

> r + l 0 < pi < 1 N 

A = attainment, N = nonattainment 

TABLE 2-8. Attainment status for the 

one-exceedance form of the standard (r = 1) 

k

i 

CVJ 

O IH iH A 

Pi Interval 

o < P i < 1 

0 < Pi < TT-j 

Ttj < Pi < 1 

o < P i < 1 

Assigned 

Status 

A 

N 

A 

N 

* 

A = attainment, N = nonattainment 

25 

83146 8 
1.0 

0.6 

-i | r "i r n 1 r 

0.6 

sz 

a 

0.4 

0.2 

0. 8 To 0. 2 0. 4 0. 6 0. 5 i 

Pi 

(a) a versus p. 

0 . 0 0. 2 0. 4 0. 6 

Pi 

1 . u 

0. B 

0.6 

0.4 

0.2 

n r, 

\ i"\ A! \ i -' 

-

iii i 

i'i ' 

-

\ 

\ 

\ -

i I- I i \ 

(b) 3 versus p. 

FIGURE 2-3. Binomial approximation to a and B 

for r = 1 and a specification of aQ = 0.2. 

8311+6 26 
P.54 

exceedances. Hence there is no chance of a Type I error. If *£ < Pi < 

113, three observed exceedances are required for a nonattainment determination, and the maximum probability of an incorrect nonattainment determination is approximately p*j. Continuing, we obtain a modified Table 2-6 with 

cutpoints 2j/7 replaced by TTJ for j = 1, 2, 3. Close approximations to 

the required n.j are obtained by solving 

a

0 - IT 3 

_ 2 

a Q - TT2 

+ 3TT2 (1 T T 2 ) 

a„ = 1 - a •i) ' 

For a0 = 0.20, we obtain 

ir 3 = 0.585 

TT2 = 0.287 

TTJ = 0.072 

Figure 2-4 shows the resulting a and 3, 

Summary of Results and Alternative Approaches 

In discussing attainment of air quality standards at an individual site 

during a specified year (or other attainment period), we have adopted the 

view that a complete set of measured ambient concentrations at the site 

for the year would unambiguously define the status of the site to be 

either that of attainment or nonattainment. We have also adopted the view 

that, however incomplete the data may be, we must come to a decision 

regarding whether, on the basis of a partial set of observations, the 

status appears to be either that of attainment or nonattainment, i.e., we 

have excluded an "indeterminate" status. Rather than prescribing a data 

8311+6 3 27 
1.0 

0.B 

0.5 

DJ 

0. 4 

0. 2 

(a) a versus p. 

0. n 1 1 1 1 

1 . " < 

(b) 3 versus p. 

FIGURE 2-4. Binomial approximation to a and 3 

for r = 3 and a specification of an = 0.2. 

8311+6 28 
o 

P.56 

0 

o 

completeness threshhold, pg, above which we may reliably asssign attainment status, we have addressed the chance of coming to an erroneous 

decision as a function of the data completeness proportion, Pi, whatever 

it may be. Moreover, we have focused on situations in which the total 

number of exceedances (x.) is just within or just above that prescribed by 

the standard, i.e., situations in which the chance of misjudging the 

attainment status is greatest. 

O 

O 

0 

Within this context, we examined three procedures for judging attainment 

status: (1) based solely on the observed number of exceedances, k ^ (2) 

based solely on the estimated number of exceedances, Xi = k./p.; and (3) 

based on both ki and the data completeness proportion, Pi. In the first 

procedure there is no possibility of a Type I error, in which the status 

is erroneously determined to be that of nonattainment. In the second 

procedure the maximum probability, a, of a Type I error is greater than 

0.6, but the probability, 3, of a Type II error, in which the status is 

erroneously determined to be that of attainment, is reduced for moderate 

data completeness porportions. In the third procedure a is limited to a 

preassigned level ag at the cost of increasing 3. 

O 

G 

The Type I/Type II error analysis does not give a definitive answer 

regarding which of the three procedures is preferable, because this 

depends on the priorities of the decisionmaker. For example, if a Type I 

error is to be completely avoided, the first procedure is preferable. It 

is hoped that the analysis will nevertheless be useful in considering the 

interplay between the proportion of data completeness and the likelihood 

of an erroneous attainment determination. 

G 

® 

Frank and Curran (1982) present a similar analysis in the case of particulate matter, which is typically sampled once every six days. Their 

results correspond to those we have obtained (although they explicitly 

treat a five-exceedances-per-year, three-year standard, which we treat 

only implicitly in Table 2-4, setting r = 15). Moreover, the theoretical 

probabilities they calculate correspond to those they obtained empirically 

CD 

8311+6 3 29 
from a subset of sites having nearly complete data. Their study emphasizes a comparison of the first two procedures, based on ki and 

Xi, respectively. They show that when the number of exceedances is well 

above the number (r) allowed in the standard, our enhanced ability to 

detect nonattainment using the second procedure is even greater than that 

shown by the present worst-case analysis (in which Xi = r + 1) . 

A different form of theoretical analysis is presented by Davidson and 

Hopke (1983). The authors construct a probability distribution of 

pollutant concentrations at a given site having a characteristic concentration, c, that is exceeded on average once per year. They define the 

status of a site as "attainment" if c is at or below the standard; 

otherwise it is "nonattainment." 

It should be noted that theoretical analysis of the most sensitive 

scenarios is not the only possible way to view the interplay between data 

completeness and procedures for inferring attainment status. An alternative is to empirically estimate the probability of misclassifying a site 

with a given proportion of data completeness. This approach is discussed 

in Section 3. 

There are other procedures for determining attainment. Breiman, Gins, and 

Stone (1978) substantially ameliorate the effects of incomplete data using 

a different estimate, Xi s of the actual number of exceedances, \j. Their 

approach is to fit an exponential distribution to the highest 20 percent 

observed concentrations and to use this fitted tail exponential distribution to estimate the expected number of exceedances among the missing 

concentrations (also see Breiman and Stone, 1983). It might prove 

worthwhile to include a procedure based on Xi in future investigations of 

the effects of incomplete data on the chance of misclassification. 

8311+6 3 30 
EMPIRICAL ESTIMATION OF THE PROBABILITY OF 

MISCLASSIFYING ATTAINMENT STATUS 

LONG-TERM AMBIENT STANDARDS 

Using the statistical framework introduced in Section 2, we consider the 

ith site to have attained a long-term ambient standard if the average of a 

complete set of concentrations, p i s does not exceed the reference 

concentration, C^^gg. Typically, a complete set of concentrations is not 

recorded, i.e., the data completeness proportion (Pi) is less than unity, 

so that the average«of the recorded concentrations (p.) is used instead of 

Pi to classify the site as attainment or nonattainment. 

If Pi is less than unity and the site has in fact attained the standard, 

there is a chance, denoted by a* , that the site will be misclassified as 

nonattainment. Similarly, if the site has in fact failed to attain the 

standard, there is some chance, 3i, that the site will be misclassified as 

attainment. Using the approximate normality of \x- discussed in Section 2, 

if *i < CNAAQS w e hav e 

a

i

 = P r (yi > CNAAQS} 

* 1 

' NPi 

,1/2 

'NAAQS M

i 

= * 

N p i \ ^ i ' CNAAQS 

1 l 

8311+6 10 31 
and if ui > CNAAQ S

 we have 

*1

 = Pr K K CNAAQS} 

'l»MAAnC ~ \-

i 

Setting 

-[(Ar 

Np,- \1/2 in - C 

NAAQS

 Mi 

a. 

NAAQS 

Ci = \ 1 - Pi 

i "NAAQS 

where a- is the standard deviation of the recorded concentrations, we may 

estimate Oi and 3. as 

«i

 = $(? i} if yi * SAAQS 

* i

 = * ( _ « 1 ) if ui > CNAAQS 

We then estimate a(p0) and 3(pQ ), the probability of misclassifying a site 

whose actual status is attainment and nonattainment, respectively, as 

a(pQ) = average Oi, p^ > PQ and ^ < C NAAQS 

3(p0) = average iv p. >.pQ and M. > C NAAQ S 

8311+6 10 32 
P.60 

The calculation of a(pg) and 3(pQ) requires that each site be tentatively 

classified as attainment or nonattainment based on u^ prior to the 

calculation of a.: or 3.. This tentative classification may be erroneous 

for a given site, particularly if ^i is near C NAA Q S , but averaging over a 

large number of sites produces estimates a(pQ ), 3(P0) whose reliability is 

much greater than the corresponding estimates, Oi or 3., at a given site. 

Table 3-1 summarizes the results of these calculations for sulfur dioxide, 

nitrogen dioxide, lead, and total suspended particulate matter. In some 

cases no estimate is given because the number of sites is insufficient for 

reliable estimation. In the cases where reliable estimation is possible, 

we see that a(pg) and e(po) are less than one percent. For sulfur dioxide 

and nitrogen dioxide, data completeness requirements of 50, 75, and 95 

percent were examined. For lead and total suspended particulate, which 

are typically monitored once per six days, we examined data completeness 

requirements of 5 percent (corresponding to the current requirement of 

five recorded concentrations per quarter), 12 percent (which may be 

interpreted as requiring a 75 percent data recovery from a once per six 

day sampling schedule), and 25 percent (75 percent recovery from a once 

per 3 day sampling schedule). 

When interpreting Table 3-1 it is important to distinguish between rates 

of misclassification that are possible under various data completeness 

requirements and rates that are estimated from 1981 data. The distinction 

is especially important for lead and particulate matter since their onceper-six-day sampling schedules limit the data completeness proportion to 

one-sixth. Figures 3-1 through 3-3 illustrate this point for, respectively, 

the quarterly lead standard, the annual TSP standard (based on the annual 

geometric mean concentration), and a hypothetical TSP standard based on 

the annual arithmetic mean concentration. The last figure is included 

because an annual particulate (PM 10) standard based on an arithmetic mean 

is currently under consideration. 

Each figure consists of a set of curves, a set of plotted points, and a 

bar chart. Each curve shows possible misclassification probabilities for 

831i+6r 10 33 
P.61 

TABLE 3-1. Estimated probability that a site's attainment status for longterm standards is misclassified as a function of minimum data completeness 

proportion, pg. The table entries are the estimated percentages of attainment sites (A) meeting a given data completeness requirement, p0, likely to 

be misjudged as nonattainment, or the percentage of nonattainment sites (N) 

likely to be misjudged as attainment. No estimate is made if the number of 

sites is insufficient. (Note: there are no known specific instances of 

misclassified sites--see text for a discussion of the estimation procedure.) 

S02 (annual): 891 sites 

A (99%)c 

N (1%) 

Assigned Status 

Actual p0 = .50 (85%V 

Status A N 

< 1%C 

p0 = .75 (71%) 

A N 

< 1% 

p0 = .95D (8%) 

A N 

< 1% 

N 02 (annual): 324 sites 

Actual 

Status 

A (97%) 

N (3%) 

p0 = .50 (78%) 

A N 

< 1% 

p0 = .75 (54%) p0 = .95 (8%) 

A N A N 

< 1% < 1% 

Pb (maximum quarterly): 696 sites 

Actual p0 = .05 (61%) 

Status A N 

A (97%) 

N (3%) 

< 1% 

p0 = .12 (39%) p0 = .25 (1%) 

A N AN ' 

< 1% 

TSP (annual): 3654 sites 

Actual p0 = .05 (92%)e 

Status A N 

A (87%) 

N (13%) 

2% 

9%1 

p0 = .12 (77%) p0 = .25 (6%) 

A . N A N 

2% 

8% 

< 1% 

3% 

Notes on next page. 

831i+6r 8 34 
P.62 

Notes for Table 3-1 

The number in parentheses is the percentage of sites meeting the required 

data completeness proportion, pg. 

By 0.95 the value y^- = 0.9524 is meant. 

The number in parentheses is the estimated percentage of sites that have 

attained or failed to attain the standard. Estimates are based on reported 

long-term averages. 

This is the estimated a-probability that a site in attainment and meeting 

the data completeness proportion p0 = 0.5 would be misjudged 

as nonattainment. This entry can also be interpreted as an estimate of the 

percentage of such sites expected to be misclassified. 

Note that intermittent monitoring in itself results in less than complete 

data. The data completeness requirements of 0.05, 0.12, and 0.25 correspond 

approximately to five recorded 24-hour measurements per quarter, and 75 

percent data recovery from once-per-six-day and once-per-three-day sampling 

schedules, respectively. 

This is the estimated 3-probability that a site not in attainment of the 

standard and meeting data completeness proportion pg = 0.05 would be 

misjudged as attainment. Alternatively, we can interpret this entry as an 

estimate of the percentage of such sites expected to be misclassified. 

831i+6r 8 35 
P.63 

(A) f n = 1/18, N = 91 ; (1/18) < = p. < (1/12), 50 sites 

U . I B. 

E.5 : 

5.A 

s

 E

-

5

-

"S B.4 -

U ? c] r-

*. B.2 ; 

o 

~ 0.1 r 

E.E 

o 

£ -e. e.-3 

o/y = 0.36 (IOJ tile ) J 

o/y • 0.59 (50% tile ) ' 

o/.w « 0.93 (90% tile ) 

:. ^ 

o 

E.5 

E.4 

E.S -

E.2 

E . : 

e. 7 AH-: 

(B) f n = 1/12, N = 91; (1/12) < = p. < (1/8), 86 sites . 

U IE . 

o/y = 0.36 (10% tile) _' 

alx. = 0.59 (50% tile) \ 

o/y = 0.93 (90% tile) ~ 

J 

i 

i 

E.3 i . 2 2 . 4 3 . 3 

(C) pQ = 1/8, N= 91; (1/8) < = p. < = (1/6), 180 sites . 

E. 

E. 5 

e . 4 :-

o/y = 0.36 (10% tile ) -j 

o/y = 0.59 (50% tile ) ! 

o/y = 0.93 (90% tile ) -I 

>> 

£ 0.1 

.0 0. 0 -SSST 

O 

^•i , r. 

0. 0 

1 . 2 

Figure 3-1. Lead: the probability of misjudging actual attainment status as 

determined by the mean (y), the coefficient of variation {o/v), and the 

proportion of data completeness (p). Continuous curves show probabilities 

determined by fixed combinations of p and o/\i, whereas the plotted points 

are determined by combinations obtained from 1981 data. 

35 
(A) pQ = 1/18, N = 365; (1/18) < = p. < (1/12), 199 sites . % 

g 0. 6 

S B.5 

* 0. 4 

ms 

S 0.3 

% B.2 

Percent 06s 

se. 

I 

"S 0. 0 

o 

-••J . 30. 0 

l n og/l n yg = 0.08 (10% tile ) 

l n og/l n yg = 0.10 (50% tile ) 

l n og/l n y g = 0.15 (90% tile ) 

_L 

60. 0 90. 0 

_L 

120.0 150.0 

(B) pQ = 1/12, N = 365; (1/12) < = pi < (1/8), 432 sites . 

g 0- 6 

Percent Dbs 

se. 

S 0. 5 

mm 

« 0. 4 -

tn 

a 0- 3 

»- 0. 2 -

o 

>> 

£ 0. I 

o 

0 . 0 

°- -0. 0 . 0 

I 

30. 0 

mmH*m& 

I 

In og/l n yg = 0.08 (10% tile ) 

In og/l n yg = 0.10 (50% tile ) 

l n o g /ln'y g = 0.15 (90% tile ) 

t,'":i<<""K * MX X 

_L 

60. 0 90. 0 120.0 150.0 

"SPMtP 

(C) pn = 1/8 , N = 365; (1/8 ) < = p. < = (1/6), 2125 sites, 

c 0- 6 

o 

Percent DDs 

5B. 

In og/ln y g = 0.08 (10% tile ) -

l n og/l n yg = 0.10 (50% tile) _ 

I n og/ln yg = 0.15 (90% tile ) 

•aKOs-^WB^seot;—* act 

120.0 150.0 

Figure 3-2. TSP, annual geometric mean: the probability of misjudging 

actual attainment status as determined by the natural log of the geometric 

mean (yg ). the coefficient of variation (ln ag/ln y g ) , and the proportion 

of data completeness (p). Continuous curves show probabilities determined 

by fixed combinations of p and (In Og/ln yg), whereas the plotted points are 

determined by combinations obtained Trom ly81 data. 

37 
(A) pQ = 1/18, N = 365; (1/18) < = p. < (1/12), 199 site s Percent 009 

SE. 

30. 0 60. 0 90. 0 120. 0 150. 0 

(B) pQ = 1/12, N = 365; (1/12) < = pi < (1/8), 432 sites, 

g 0-6 | 

Percent Obs 

5B. 

0.5 -

<« 0.4 -

o 

0 . 3 

0 . 0 

o 

°- -0 . 0 . 0 30. 0 

o/y = 0.32 (10% tile ) 

o/y = 0.42 (50% tile ) 

o/y = 0.61 (90% tile ) 

60. 0 90. 0 

x vx xw. jftfe&^MMfafe.—;c 

120.0 150.0 

(C) pf 

I

 B6 

S 0.5 

1/8, N = 365; (1/8) < = Pi < (1/6), 2125 sites, Percent Obs 

5B. 

;-<M. i? 60.0 90.0 120.0 150.0 

IBB. 

1 I 

ill 

1 1 

o/y = 0.32 (10% tile) -

o/y = 0.42 (50% tile) 

o/y - 0.61 (90% tile) 

1 1 

Figure 3-3. TSP, annual arithmetic mean: the probability of misjudging 

actual attainment status as determined by the mean (y), the coefficient of 

variation (o/y), and the proportion of data completeness (p). Continuous 

curves show probabilities determined by fixed combinations of p and o/y, 

whereas the plotted points are determined by combinations obtained from 

1981 data. 

38 
P.66 

a range of mean concentrations by combining a given data completeness 

proportion (1/18, 1/12, or 1/8) with a given coefficient of variation 

(obtained from the 10th, 50th, or 90th percentile of the 1981 distribution). The given data completeness proportions approximately correspond 

to 5, 8, and 12 samples per quarter, respectively. The last proportion, 

1/8, may also be interpreted as a 75 percent data recovery rate from a 

once-per-six-day sampling schedule. 

The plotted points in each figure represent sites whose data completeness 

proportions fall with a specified interval. Parts (A), (B), and (C) of 

each figure represent completeness intervals of 1/18 to 1/12, 1/12 to 1/8, 

and 1/8 to 1/6, respectively. For each site the estimated misclassification probability (in part determined by the coefficient of variation 

reported for that site) is plotted against the reported mean concentration. Within each data completeness interval the bar chart shows the 

percentages of sites having a given probability of misclassification. 

The figures show that the data completeness requirements by themselves do 

not limit probabilities of misclassification. It can be seen, however, 

that for the great majority of 1981 sites, reported mean values are 

sufficiently far above or below the standard and coefficients of variation 

are sufficiently small so that the estimated probabilities are effectively 

zero. Thus, misclassifications are not ruled out by these data completeness requirements but they are unlikely, assuming that the means and 

coefficients of variation reported in 1981 are representative of other 

years. 

We can also use the figures to infer the benefits of requiring 12 rather 

than 5 samples per quarter from intermittent samplers. In the case of 

lead we are assessing the benefits of current requirements since 12 

samples per quarter are required to determine attainment status (see Table 

G-l). The bar charts of parts (A) and (C) of Figure 3-1 show that the two 

groups of sites have nearly the same low susceptibility to misclassification. The curves distinguish parts (A) and (C) more clearly, however, and 

are probably more indicative of the future consequences of strengthening 

83146r 1 0 

39 
or relaxing data requirements for all sites. We can express this difference by comparing "critical concentration intervals," i.e., by noting the 

interval of maximum quarterly mean concentrations within which the 

probability of misclassification exceeds 0.10 for a site having a typical 

(i.e., median) coefficient of variation. 

This definition of a "critical interval" bears some explanation. We are 

interested in correctly classifying the attainment status of a site whose 

actual long-term average concentration, y, is estimated by y based on 

incomplete data. For simplicity, we assume the coefficient of variation, 

o/v, to be the median of the coefficients of variation calculated for our 

sample of sites. If y is very close to the reference concentration, 

^NAAQS'

 the n u will be above or below CNAA Q5 with nearly equal probability. Therefore, if y is only slightly above or below CN A A Q $ , the risk 

of misclassification will be close to 50 percent. As y moves away from 

^NAAQS»

 t^1e probability of misclassification decreases until, at some 

point above and below C^A AQ5, the probability drops below 10 percent. The 

critical interval is thus the set of y for which attainment determinations 

are subject to error with nonnegligible probability. 

In mathematical notation we define the critical interval as the union of 

intervals A and B where 

1

 M : p < CNAAQS an d Pr y (p > CNAAQS) > a * 

B = { y : y > C NAAQS and Pru(y < CNAAQS) > a j 

so that the critica l interval is 

A U B = { y : a < Pr^y" > CNAAQ$) < 1 - a } . 

In the framework of testing the statistical hypothesis that y < ^NA Ans 

based on y, the critical interval is the set of y for which the power of 

83m6r io 40 
o 

P.68 

O 

o 

o 

o 

the test, i.e., the probability that y exceeds CNAA g $, is between 

a and 1 - a, i.e., setting a equal to 0.1, between 10 and 90 percent. 

We may also describe the critical interval as the interval within which 

the probabilities of apparent attainment and nonattainment (based on a 

comparison of y and CwA Age) are both between 10 and 90 percent. It is 

thus an interval within which the determination of attainment status is 

unreliable. We may therefore use critical intervals to judge the adequacy 

of completeness requirements--the more stringent the requirement, the 

narrower the critical interval. 

The critical interval for a site recording only five valid lead samples 

per quarter is 1.13 to 2.23 yg/m (25 percent below to 49 percent above 

the standard of 1.5 yg/m3 ), whereas the interval narrows to 1.24 to 1.90 

yg/m (17 percent below to 27 percent above the standard) for a site that 

records 12 valid samples per quarter. 

O 

O 

O 

O 

O 

Figures 3-1 and 3-2 support similar conclusions regarding TSP long-term 

standards. Since long-term average concentrations are well below the 

standard at most sites, the benefits of 12 versus five samples per quarter 

are clearer from the curves depicting possible misclassification 

probabilities than by a direct comparison of sites with different proportions of data completeness. We consider two forms of annual TSP standard, 

both having a reference concentration of 75 yg/m , namely, the current 

standard based on the geometric mean and a hypothetical standard based on 

an arithmetic mean (the latter is relevant to the recently proposed PM-10 

standard). For the current standard, the critical interval for the case 

of five valid samples per quarter is 67 to 85 yg/m (11 percent below to 

13 percent above the standard), which narrows to 70 to 81 yg/nr (7 percent 

below to 8 percent above the standard) for the case of 12 valid samples 

per quarter. These results are the same (to two-digit accuracy) for the 

hypothetical standard concentration represented in Figure 3-3. 

It is instructive to compare the cases of lead and TSP represented by 

Figures 3-1 and 3-2, respectively. The probability of misclassification 

is determined not solely by the proportion of data completeness, but also 

O 831i+6r io 41 
P.69 

O 

by the number of valid observations to be used in judging attainment. 

Q For example, 5.5 percent complete data yields 20 observations with 

which to judge attainment of an annual standard, but only five 

observations for assessing the attainment of a quarterly standard. 

Therefore, a given number of lead samples per quarter with which to judge 

Q attainment of the quarterly lead standard leads to higher 

misclassification probabilities than the same number of TSP samples per 

quarter (i,e. , four times that number per year) with which to judge 

attainment of the annual TSP standard, other factors being equal. Thus, 

Q the critica l interval corresponding to the 12 lead samples per quarter 

required for attainment determinations is relativel y wider (17 percent 

below to 26 percent above the quarterly lead standard) than the interval 

corresponding to the required 5 TSP samples per quarter (10 percent below 

O to 13 percent above the annual TSP standard). 

Therefore, i f our goal were to contain the 0.10 probability critica l 

interval within , say, + 10 percent of the standard, then the requirement 

Q of five TSP samples per quarter nearly achieves the goal whereas the 

requirement of 12 lead samples per quarter does not. This suggests that 

any tightening of current data completeness requirements to improve 

attainment assessments of long-term standards should begin with lead,. 

O 

On the other hand, i f our goal were to restric t the fraction of sites 

vulnerable to misclassification, e,g. , sites for which the misclassification probability exceeds 0,10, i t appears from the 1981 data that any 

O tightening of completeness requirements should be directed at TSP rather 

than lead, since the observed fraction of "vulnerable" sites is greater 

for TSP than for lead. This difference between TSP and lead may be solely 

a consequence of the smaller sample of lead sites , however, 

O 

The two goals thus lead to different conclusions about which pollutant, 

lead or TSP, warrants the most consideration for strengthened data 

completeness requirements. Since the "critica l interval" analysis seems 

O to be more reliable than the "vulnerable sites" analysis, i t seems that 

requirements for lead warrant firs t consideration. This is especially so 

since strengthed requirements for particulates have already been 

O 83ii+6r io

 4 2 
o 

P.70 

O 

considered; proposed revisions to the particulate NAAQS include a 

requirement of 12 valid observations per quarter (EPA, 1984). 

:o 

O 

G 

The recommendation that strengthened completeness requirements be considered for lead is offered only as a point of departure, however, since 

realisti c recommendations must include many practical considerations 

outside the scope of our analysis. For example, we have not attempted to 

distinguish in general between data completeness proportions resulting 

from good data recovery from intermittent sampling, on the one hand, 

versus poor data recovery from continuous sampling on the other. The 

latte r case may be indicative of generally unreliable data. Pending a 

more comprehensive review of data completeness requirements, however, the 

results presented in Table 3-1 and Figures 3-1 through 3-3 should be 

useful because they show that current requirements are adequate to render 

reliable determinations of the attainment of long-term standards at al l 

but a small fraction of sites . 

G 

SHORT-TERM AMBIENT STANDARDS 

O 

G 

G 

G 

Correctly classifying the short-term attainment status of a site having a 

certain data completeness proportion is much more difficul t than classify -

ing the long-term attainment status. That is , the estimated long-term 

average is a more stable statisti c than the observed number of short-term 

exceedances or the estimated total number of exceedances. Because of 

this , i t is not possible to reliably estimate the chance of misclassifying 

the status of most sites . We can, however, estimate the average chance of 

misclassifying sites whose proportion of data completeness meets a 

required proportion. 

This can be done i f we make the reasonable assumption that the actual 

number of short-term exceedances (X.) and the data completeness proportion 

(Pi) are statisticall y independent. That is , we may use an estimate of 

the total number of exceedances at each site to estimate the statistica l 

distributio n of short-term exceedances across al l sites . We then apply 

831<+6r i o 43 

G 
P.71 

this statistical distribution to those sites having data completeness 

proportions, Pi, near a given proportion, p, and thereby estimate the 

fraction of these sites having a given number of exceedances 

(x = 0, 1, 2, ...) . For a one-exceedance standard (for which r = 1 in 

the notation of Section 2), we may therefore estimate the probability of 

misclassified attainment and nonattainment sites, respectively, a and 3 , 

r r 

associated with a given proportion of data completeness. The estimates of 

a and 3 are 

ap = Pr(xi > 1.05 | pi = p, Xi = 1) Pr(xi = 1 | Xi < 1) 

= p-Pr(xi = 1 | xi < 1), i f p <Y-^ -

° '

 1f P>

T^5 ' 

and 

Xmax 

3p = E Pr (

x

i < i- 0 5 I Pi = P. *i

 = *) Pr U -j = X | X. > 2) 

X

- ^ 

^

X

 \ « ' ' 1 

= E d - P) Pr(x. = x | x. > 2), i f p <T±E 

A ~ £-

max \ \ * 1 

= E CU " P)

 + *Pd " P) ] p^(Xi = x | xn. > 2), i f p > y-i ^ 

We then estimate a(pQ) and 3(pQ) as 

A A 

a(p0) = average ap , ' p1 > pQ 

A A 

3(p0) = average 3p , Pi > P0 

In estimating the statistical distribution of the total number of exceedances, instead of using Xi = k./p- we use X-, obtained by fitting an 

83ii+6r

 10 44 
exponential distribution to the upper 10 percent of the distribution of 

concentrations at a site (the exponential tail is fitted using available 

statistics—the 90th percentile and the maximum concentration). This is 

A 

done to avoid a certain lumpiness in the distribution of Xi (e.g., the 

minimum positive estimate obtained from a once per six day sampling 

schedule is 6) , and results in a larger estimate in several cases. The 

formula for x is 

X = k + (N - n + 1)(1 - Y ) exp 

r/CNAAQS " V ln 

max 

1 - y 

where we have suppressed the index i for ease of notation and where 

k = 

N = 

n = 

*Y = 

x^ = 

"Y 

<max 

n 

the number of recorded exceedances 

the number of observations in a complete data set 

the number of recorded concentrations 

0.9, corresponding to the use of the upper 10 percent of the 

concentration distribution 

x0.9» the ^t h percentile of the concentration distribution 

the maximum recorded concentration 

n/(n + 1) 

Provision is made in the calculation of X for cases in which either x > 

"Y 

CNAAQS or xmax = x * In these cases x was set equal to x = k/p. 

It should also be noted that x is rounded to the nearest integer to estimate 

the statistical distribution of the total number of exceedances among all 

sites. 

The results of these calculations are shown in Figures 3-4 through 3-8 and 

A 

in Table 3-2. The figures indicate the chance a that a site randomly 

selected from a group of attainment sites having data completeness 

A 

proportion p would be misclassified as nonattainment, or the chance 3 that 

a nonattainment site would be misclassified. In contrast to the most 

sensitive scenarios treated in Section 2, Figures 3-4 through 3-8 indicate 

83146 T 10 45 
(a) a versus p 

0.06 

0.06'r 

CL 

0.04|-

0.02-

o.c 

(b) 3 versus p 

0.08 

0.06 

-i i 1 r 

CD 

n 

0.04 

0.02 

0. OJ 

FIGURE 3-4. Estimated probabilities a and 3 of 

misjudging attainment of the 24-hour sulfur dioxide 

standard as functions of data completeness proportion, p. 

8311+6 46 
(a) a versus p 

O.OB -

0.06 r 

0.04 

0.02 u 

0.00 

0.2 0.4 0.6 0. 8 1.0 

(b) 3 versus 0 

1 1 r 

0.06 -

0.06 

X) 

0.04 -

0.02 1-

0.0 

1 ;0 

FIGURE 3-5. Estimated probabilities a and 3 of 

misjudging attainment of the one-hour carbon monoxide 

standard as functions of data completeness proportion, p. 

8311+6 47 
• 

0.08 

0.06 -

(a) a versus p 

-1 1 1 1——r 1 n r 

ti 

sz 

D 

0.04 

0.02 -

(b) B versus p 

0.08 

0.06 

QJ 

0.04 

0.02 L 

0.0 

O.S 1.0 

FIGURE 3-6. Estimated probabilities a and 3 of 

misjudging the 8-hour carbon monoxide standard as 

functions of the data completness proportion, p. 

8311+6 48 
P.76 

1.0 

0.8 

0.6 -

x: 

Q 

0.4 -

0.2 

(a) a versus p 

-i 1 1 1 1 1 r 

FIGURE 3-7. Estimated probabilities a and 3 of 

misjudging attainment of the maximum daily one-hour 

ozone standard as functions of the data completeness 

proportion, p. 

8311+6 49 
(a) a versus p 

0.08 

0.06 

o 

SZ' 

a 

0.04 

0.02 

0. ot 

-i 1 1 r -i 1 r 

0.06 -

0.06 -

o 

X) 

0.04 -

0.02 -

0.0 

FIGURE 3-8. Estimated probabilities a and 3 of 

misjudging attainment of the 24-hour TSP standard as 

functions of the data completeness proportion, p. 

8311+6 50 
P.78 

TABLE 3-2. Estimated probability that a site's attainment status for shortterm standards is misclassified as a function of minimum data completeness proportion, pg. The table entries are the estimated percentages of attainment sites (A) meeting a given data completeness requirement, p0 , likely to 

be misjudged as nonattainment, or the percentage of nonattainment sites (N) 

likely to be misjudged as attainment. No estimate is made if the number of 

sites is insufficient. (Note: There are no known specific instances of 

misclassified sites—see text for a discussion of the estimation procedure.) 

S02 (24 hours): 891 sites 

Assigned Status 

Actual p0 = .50 (85%)' 

Status A N 

A (98%) c 

N (1%) 

2%c 

p0 = .75 (71%) p0 = .95° (8%) 

A N A N 

2% 0% 

CO (1 hour): 480, sites 

Actual 

Status 

A (99%) 

N (1%) 

p0 = .50 (79%) 

A N 

< 1% 

p0 = .75 (65%) pg = .95 (22%) 

A N A N 

< 1% 0% 

CO (8 hours): 480 sites 

Actual 

Status 

A (57%) 

N (43%) 

p0 = .50 (79%) 

A N 

2% 

< l%e 

p0 = .75 (65%) 

A N 

2% 

< 1% 

p0 = .95 (26%) 

A N 

0% 

< 1% 

O3 (Daily maximum 1-hour): 723 sites 

Actual 

Status 

A (63%) 

N (37%) 

p0 = .50 (86%) 

A N 

p0 = .75 (70%) p0 = .95 (16%) 

A N A N 

8% 8% 

< 1% < 1% 

0% 

< 1% 

Continued 

831H6T 8 bl 
P.79 

O 

O 

o 

o 

o 

o 

TABLE 3-2. (Concluded) 

TSP (24-hours): 3667 sites 

Assigned Status 

Actual p0 = ,05 (92%)f pQ = ,12 (77%) p0 = ,25 (6%) 

Status A N A N A N 

A (79%) 2% 2% 3% 

N (21%) . . 6% 6% 4% 

a The number in parentheses is the percentage of sites meeting the required 

data completeness proportion, pg. Under current procedures for determining 

the attainment of short-term standards, this completeness requirement rules 

out the possibility of a false nonattainment determination. 

b By 0.95 the value -—^ = 0.9524 is meant. 

c The number in parentheses is the estimated percentage of sites that have 

attained or failed to attain the standard. These estimated percentages are 

based on the tail-exponential method discussed in the text. 

^ This is the estimated a-probability that a site in attainment and meeting 

the data completeness proportion p0 = 0,5 would be misjudged 

as nonattainment. This entry can also be interpreted as an estimate of the 

percentage of such sites expected to be misclassified. 

e This is the estimated 3-probability that a site not in attainment of the 

^ standard and meeting data completeness proportion pg = 0.5 would be 

misjudged as attainment. Alternatively, we can interpret this entry as an 

estimate of the percentage of such sites expected to be misclassified. 

Note that intermittent monitoring in itself results in less than complete 

O data. The data completeness requirements of 0.05, 0,12, and 0,25 correspond 

approximately to five recorded 24-hour measurements per quarter, and 75 

percent data recovery from once-per-six-day and once-per-three-day sampling 

schedules, respectively, 

O 

831i+6r 8 52 

O 
P.80 

that a and 3 are for the most part well below 10 percent. This is 

because most sites are estimated to have no exceedances, in which case 

there is no chance of misclassification, and those sites estimated to have 

more than one exceedance frequently appear to have more than two exceedances. In short, the most difficult cases treated in Section 2 appear to 

be rare in actuality. It should be noted that the data are from 1981, so 

that, for purposes of comparison, ozone is treated as having the same oneyear one-exceedance form of standard as the other pollutants. 

Table 3-2 shows an approximate 1 or 2 percent probability of misclassifying sulfur dioxide and carbon monoxide sites. The a- and 3-probabi-lities 

of misclassification are estimated to be 8 percent and less than 1 

percent, respectively, for ozone sites. For particulate sites a and 3 

are estimated to be 2 and 6 percent, respectively. 

SUMMARY 

The 1981 data indicate that under current data completeness requirements, 

the rate of misjudging the attainment of either long-term or short-term 

standards is under 10 percent, and in several instances the estimated 

misclassification rate is under one percent. The current completeness 

requirements therefore seem to be adequate for generally reliable attainment determinations. This conclusion should be regarded as tentative, 

however, pending development of specific objectives for limiting rates of 

misclassification. 

We have also conducted what must be regarded as a preliminary analysis of 

the benefits of strengthened requirements for data completeness. In 

general, direct estimates of the expected percentages of misclassified 

sites reveals virtually no change resulting from the changes in requirements that we consider. Taking a less direct approach in the case of 

long-term standards, however, we have estimated a critical interval about 

the standard, within which misclassification probabilities exceed 10 

percent, and we have shown how these critical intervals narrow as more 

8311t6 r 10 53 
P.81 

complete data are required. The analysis suggests that strengthened 

requirements be considered for lead and possibly for particulates (the 

recently proposed revisions to the particulate standards include such a 

strengthened requirement). A similar analysis is possible for short-term 

standards, although the analysis would require some assumptions about the 

frequency distribution of concentrations at individual sites. 

In any case, realistic recommendations for strengthened data completeness 

requirements must address the broader issues concerning data reliability 

and cost-effective air quality management. For example, we should 

probably consider the benefits of strengthened quality assurance procedures along with strengthened data completeness requirements. Similarly 

we might consider whether more accurate attainment determinations might be 

achieved under current data completeness requirements by changing the form 

of current short-term standards. We have not attempted to include such 

issues in the present analysis, although these issues certainly deserve 

future consideration. 

831i+6r 10 54 
P.82 

STATISTICAL ESTIMATION OF NATIONAL 

AIR QUALITY STATUS AND TRENDS 

We now consider the estimation of national averages of site-specific 

annual averages, y., and site-specific exceedance rates, x, . We face the 

choice of a minimum data completeness proportion, pg, such that only those 

sites in-our sample with p. > pn will be used.to estimate the national 

average. 

A QUESTION OF BIAS 

Let 

1

 v(Po}

 . 

p(Pg) = vTpgT E W-j 

and 

1

 v(p o}

 . 

X(Pg) = v(pQ) E X i 

A A 

where y. and x. are estimates of the site-specific long-term average, M^, 

and exceedance rate, X-, respectively, and where the summation is over 

sites in.the sample with p. > p Q . If we assume that from site to site the 

air quality parameters yi and \. are independent of the proportion pn- of 

data completeness, and if we regard the set St of sites in our sample as 

having been randomly drawn from the population of sites of interest, then 

A A 

for each pg, y(pQ) and x(pQ) should be unbiased, i.e., we should have 

A 

E y(pQ) = y = E(U i ) 

8311+6 it 55 
P.83 

and 

E x(pQ) = X = E (x.) . 

A A 

This lack of bias follows from the unbiasedness of yi and x., 

A 

E Ui = Ui 

E x i = X i 

as discussed in Section 2. 

A A 

Under these assumptions, y(pQ) and x(pQ) are approximately normally 

distributed unbiased estimators of y and x, respectively. Thus, a natural 

criterion to use in selecting pg is the variance of these estimators. 

VARIANCE OF THE ESTIMATORS 

A A 

In Chapter 2 we obtained expressions for the variance of y, and x., 

conditional on the choice of site, i, and thus conditional on y., x•, and p. 

I I A 1 

Extending the conditional variance to the averages taken over the set trf 

of sites in our sample, we have 

Var(„(pn) | fi) «-i E -s- L4 

V 

2

(P„) ft pi

 N 

and 

1

 V( 'o)l " P i Var (x(p0 ) |^ ) --jr-L. - E -y" 1

 x. 

v (pQ) 1=1

 K i 

8311+6 i+ 56 
P.84 

To obtain the unconditional variance, we must allow for site-to-sit e 

differences in y. and X-. Thus (see, for example, Brillinger, 1966) 

Var (y(pQ)) = E {var Q(pQ) \ $)} + Var{E(y(p0) \ fr)} 

W. 1 - Pl of 

V (Pg) 1 = 1 H l 

v(p0 ) 

MPg) fax 1 

which may be estimated as 

A A 

(Pg) 

Var (M(PQ)) = ~2 " E 

1 - p. a i "i , 1 

v(p0) 

v (pQ) i=l N

+

 2 v'(Pg) 1=1 

E (l^ - u(Pg))' 

Similarly , we have 

X(p0) 1

 " Pi 

v(Pg) 

Var U(p ) = Ej-g-^ - E -5— L^ I ?

 + Var]7TFT £ x i 0

 lv (p0) 1=1

 p i M r[po> 1 = 1 1 

which may be estimated as 

A A 

Var (x(pQ) = - £ 1 ^'-p , 

-

L— E —1— - ^ + 

v(Pg) 

V (Pg) i= l

 H l 

i + 7 ^ S ui-i(po, v 

The variance expressions are similar for y(pQ) and x(p Q ). In both cases, 

the second term may be expected to be proportional to l/v(p Q ) while the 

firs t term depends on the distribution of Pi (as p0 increases to unity , 

(1 - Pi)/Pi approaches zero but l/v(p Q ) increases). 

8311+6 i+ 57 
AN EMPIRICAL ASSESSMENT OF THE DATA COMPLETENESS 

NEEDED TO RELIABLY ESTIMATE NATIONAL TRENDS 

SULFUR DIOXIDE 

Table 5-1 and Figure 5-1 show y(pQ ), the estimated composite annual 

average sulfur dioxide (S0£) concentration, and Var (y(pn)), its estimated 

variance as discussed in Section 4, for several values of required data 

completeness proportion pg. The estimator y(pQ) is fairly constant with 

p0, but its estimated variance increases substantially for pQ > 80 percent, consistent with the discussion in Section 4. 

Similarly, Table 5-2 and Figure 5-2 show a substantial increase in the 

A 

estimated variance of x(pQ) for pQ > 80 percent. In addition, Figure 5-2a 

suggests that increasing the required percentage, pg, may induce a 

positive bias in the estimate of the national average exceedance rate. 

The increase in X(pQ) accelerates for pg > 80 percent. 

In summary, it appears that no more than an 80 percent data completeness 

requirement should be used to estimate the national status and trend in 

SO2 concentrations. 

CARBON MONOXIDE 

Table 5-3 and Figure 5-3 show the estimated exceedance rate, x(pn ), and 

A A *•* 

its estimated variance, Var (x(pQ)), for hourly averaged carbon monoxide 

(CO) concentrations. Interestingly, sites with more complete data have no 

A A A 

observed exceedances, so that both x(pQ) and Var (x(pQ)) are calculated 

to be zero for pg > 94 percent. For smaller values of pg, however, both 

x(p0) and Var (x(pQ)) increase with pQ , with a fairly sharp rise in 

8311+6 9 58 
CO oo 

a 

o o 

TABLE 5-1. Estimated compo site annual average 1 

-o its estimat T3 r

C 

o 

n 

o CL 

< 3-

CVJ 

R So 

3 

.c 

i cu 

o 

c 4 -

o 

Ol +J 

1+^ 

t/

o 

1 

+-

o 

> 

i_> c: zs 4 -

1/1 r

o 

a CL 

<P-

s-

ra 

c > 

(IJ (J r

c 

a 

s_ r

> 

a 

required percentage, pQ, of data comple teness. 1 

o 

o Q-

< 3-

$_ ro 

< > 

o Q . 

< P

O - Q . 

LD rl O iH o ro 

ro ro ro ro vo r-

vo r- <x> ri CM r-

CN ro oo r-t in •«* 

r H 

OOOOO O 

OOOOO O m ^ r-i ^ <Tt m 

•>* oo r» a\ r- cn 

<* w» o oo •« 

cr> vo r^ r- vo r-

CN CN CN CM CS CN 

Ol <Tl O O Cl O 

0\ 0\ O O CTV O 

CT> O^l O O ^ O 

ffl o> O O Ol o 

co CN CN o r» in 

lOffiuiHoxNinHoioiffioor^off i 

r--oor--oooocovoCT>r-inr~CNoor-~cr>o 

r-voorocrioinoovorocMocNoinc H 

OOOCNrHiHVOOinrHOOmfOCNCOr-lf H 

CM<T>oor~voir»in^j

,

"^rcororororocoro 

rHOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOO O 

HMOCNinoinooinvorTJ'ffi'j'criv o 

vorocMinr-rHCT>vor-OiHcnr~-crvroc N 

CNrHVO

I -*CT>OOVOr~VOrOCMOOr-CNO O 

i~-r"~vovoininir>ir>ir)inir»'3''3'ir>inin 

CMCNCNCNCNCMCNCNCNCNCNCNCMCMCMC N 

a\oOToo\oimoiooa\ffiooo o 

cnocnootoioioioooMnooo o 

o\ O 0NO^O>CTlCT»OO^CT(OOOO 

siocnooxnaioioocnoiooo o 

r~CNCNCNr~roi>-r-rHooooooooinooro 

voin-<i

,

cOf-iooor-»incNOO'«j'or-i-icNror~cr>r-i'«j'o 

m m m m m o» oo oo oo ooM^r^vowmi'fnNC N 

o 

o o o 0 0 
(A) y(pQ) versus pQ (B) Var(y(pn )) versus pr 

- + 

o 

c 

iO 

cu 

XJ 

OJ 

-(-> 

rO 

E 

l/l co. 

+ + + + + + ++ + 4 + 

++

+V+ 

10. 20. 30. 10. 53 . 60 . 70. 80. 80. 

Required Percent 

1C0. 

rsi 

m _ 

Varianc e Estimated 

10 

i / i 

.-

u

 + + 

i + + 

i 

+ + + 

i 

+ 

+ 

+ 

+ 

i i ; i 

10. 20. 30. 10. 50 . 60 . 70 . 

Required Percent 

80. 90. 130. 

FIGURE 5-1 . Estimated long-term SO? air qualit y parameters 

as functions of the required data completeness, p„. 

R 1 1 'i G 
1 0 O O 0 

OD 

1—« 

. r 

cn 

CT> 

I—* 

O O Q O 

MU>uj.e».Lncy»CT\-~j^j--Jooooaocovovovovoiovo 

OjbHONlUHOWOU)00l0i^Nl00OK)UU10\v l 

towwMioiouioa)iti>iooH«5i-'vi*bomwiti>o o 

>JUl-slWUIOOOOiM-

,

vlOOIOWWa(l-

, UHDit>.W O 

il>OOOvHD^|vlO«>HIOHHHN)vlHili(X)IOvlO O 

oooovovooooovovovovovooovovoovo o 

OOOOVOVOOOOOVOVOVOVDVOOOVDVOOVD O 

OOOOOOOOOOOOOOOOOOOOO O 

HHHHMHHHHHHIOIOWIOIOWIOMWIOO I 

|OMUIUi^U19ivlOOa)VOON)U)U1IX)OvlU}OaVOt O 

•JOOUimft.WIOWIOOlCTimMVOVOWtOOIOOOO A 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOOt O 

OOOOOOOOOOOOOOOh-'MI-'l-'tOLn O 

HHMHMHtOIOIOWWit»UimOOHWWa)vlH* > 

uw*uiviiotouiy>woouiitkvi*wo-Jowwoov i 

Hitk7j^>JOOfcfcOOmO>W|\J(MOibOOU100 H 

-

TD 

o 

>- > 

-a 

o 

< > 

CU 

- s 

•>• > 

TD 

o 

centage, pn, o

f dat

a completeness. 

o 

O 

0

-

-— 

>

T

O 

^

), a

s functions o

f the require

d per
' 00 UJ 

t — < 

! -F 1

 CT1 

! 

i 

: 

i 

i 

1 

i 

i 

i i 

i 

i 

i 

• 

t 

i 

i 

i 

I 

i 

! 

i ii 

i 

OJ 

m - h 

c 

3 

n <-t

-

mU. 0 SUO - b 

rt- 3 " fD req c e- 5 d data co 

--n 

r—t cn c= PO 

m t n 

i ro 

m i n 

r-t-

3 Ote

J d m 3 " 

o rt-term SO 

ple

3 r o t fD 

3 fD VI to 

V T D 

O 

QJ _lualit

o y p QJ ramet ers 

PO fquire D d Percent PO fquire D d Percent 

sP-0 10. 20. 30, tn 60. 7 ts CD 

C"3*" 

ea' 

r a

-

Q 

sP0 0 

(S3-

10. 50. 50. 70. e0. 93. I 

r 

h 

rh 

V 

¥ 

t-t-tf 4-

+ 

tX 

f^3 

29 

Estimated Exceedance Rate 

0 .1 0.2 0.3 0.1 0.5 0.6 0.7' 0.8 0.9 1.0 

t . l r I I t I I I I 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+

 + 

Estimated Variance 

0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.10 0.15 0.50 i i i i i i I I i i +

 + 

(A) >> > T D 

O 

*-* 

< rsu fD s T D 

O 0 0 < > QJ - 5 >< > 

•a 

O 

< fD "5 

(/! 

c 

l/l 

• a 

o 

• ' 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

I 
•tr 

CT) 

tO KJ 

O O CD UT 

u>£kO\cn^JoooocoaovovoK>vov0vov0V0v o 

LH to I-J oo^HOiootooiouiMjTcriffivK s 

UUIOUOOOvllOIOOOOOUIv J 

OOOOOVOOOVOVOOVOVOV O 

OOOOOVOOOVDVOOVDVDV O 

OOOOOVOOOVOVOOVOVOV O 

OOOOOVOOOVOVOOVOVOV O 

OONNOOOSIT J 

OVOOOVOVOOV O 

OVOOOVOVOOV O 

OVOOOVOVOOV O 

OVOOOVOVOOV O 

oo 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOO O O O OOOOOOOOCD O 

OMJiviuimmmvivisjoomuirnHHooooo o 

COOOHNIOWVIIOWOOUIIONIWIOWOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

H-'MMI-'MI-'l-'l-'tOtOWh-'tOtTOOOOOOOO O 

WWUIK)WUIvliDMO\OlOWIOHHOOOOO O 

oo\owwvi,i>*ioou\oui7j«)Oi«>a)Ooooo o 

TD 

O 

>- > 

*—% 

• a 

O 

. -

< > 

cu 

-s 

."—^ 

5-> > 

•"—'N 

TD 

O 

% * 

'— 

O 

o 

3 

TD 

— - m 

tt 

<-tfD 

3 

fD 

in 

in 

• 

J

T

T

O 
; 

1 

i 

t 

i 

1 

i 

i 

| 

• 

i 

i 

1 

1 

! 

j 

> 

i 

, 

i 

! 

1 

1 

! 

i 

i 

i 

i 

: 

i 

' 

\ 

] 

1 

1 

1 

t 

| 

, 

1 

i 

I 6 d 

CD • P J l 

QJ "Tl in i—t 

rn 

-h CZ C 73 3 m 

n 

rr- t n 

->• i 

o OJ 3 • in o m - h m 

rtrt- ->• 3- 3 fD QteJ d one-hour required dat QJ O 

o 

O O QJ ir qua mpl ete 3 —' fD -"• in rt- tn << 

m 

T D 

TD QJ 

O "5 • QJ 

3 (D 

rt-

fD -5 in 

179 

Estimated Exceedance Rate 

_ 0 . 0000. 0100.0200.0 300.0100.0500.0600.0700.0800.09 0 

•^ t i i t i t . i i i 

<s-

C J " 

^ -. UJ TO s ~ fD • JD 

— J* eg-

-J • 

fD CX ui "O ' fD -5 oi O s

-

fD 

3 

r t- -J 

C i " • X 

" J " LO 

i 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ + 4- . 

ca 

Estimate Variance 

Jl.0000.0000.0010.0010.0020.0020.0030.0030.0010.0010 

S r i I l.l I I I I I I 

s-

fD £-

-Q • 

LZ _" ^ -5 s H fD • Q . 

in 

~0 C3 " fD 60. rcent Q " CD Q " 

T 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ + * 

3a 

>m > TD 

O 

S M ^ 

< CD 

-% m 

c in TD 

O 

0 0 w. ^ < > 01 -J 

^-m. V > 

^ - N • a 

o ~ .. ) ver in 

C 

l/l •a 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

i 
P.92 

Var (x(p0)) beginning at pg = 60 percent. The results suggest that no 

more than 65 percent data completeness is needed to estimate the national 

hourly CO exceedance rate. 

Table 5-4 and Figure 5-4 show results for 8-hour CO exceedances. Both the 

A 

estimated mean exceedance rate, x(pn ), and its estimated variance, 

A A ^ 

Var (x(p0)) sharply increase for pg > 90 percent, suggesting that no more 

than a 90 percent data completeness requirement should be used to estimate the 

national 8-hour CO exceedance rate. 

NITROGEN DIOXIDE 

A 

Table 5-5 and Figure 5-5 show u(pn ), the estimated composite annual 

*J A A 

average nitrogen dioxide (N02) concentration, and Var (u(pQ)), its 

estimated variance, for several values of the required data completeness 

proportion, pg. The estimator p(pQ) is fairly constant with pg, but its 

estimated variance increases substantially for pg > 90 percent. This 

suggests that no more than a 90 percent data completeness requirement 

should be used to estimate the national status and trend in N0£ 

concentrations. 

OZONE 

Table 5-6 and Figure 5-6 show the estimated exceedance rate, x(pQ ), and 

its estimated variance, Var (x(pn)), for daily maximum 1-hour ozone (Oo) 

* W rs ° 

concentrations. Both x(pQ) and Var (x(p0)) increase with pg, with a 

fairly sharp rise in Var (x(pn)) beginning at pg = 85 percent. The 

results suggest that nb more than a 90 percent data completeness 

requirement should be used to estimate the national daily maximum ozone 

exceedance rate. 

83m6 9 65 
o 

o 

o 

TABLE 5-4. Estimated mean annual rate of 

8-hour CO exceedance, x(pn), and its 

estimated variance, Var(x(pQ), as functions 

of the required percentage, PQ, of data 

completeness. 

x(pQ) Var(x(pQ)) 

O 

O 

O 

O 

99.30399 

98.78999 

97.61400 

96.90599 

95.77599 

95.06799 

93.65300 

92.48900 

91.17599 

90.17099 

88.54999 

86.92899 

80.87900 

74.13199 

68.40199 

61.34699 

44.58900 

35.26299 

25.58200 

20.04600 

0.51400 

0.29700 

43.75500 

47.91299 

49.96200 

43.69399 

39.28200 

34.01600 

34.15400 

33.28799 

35.17000 

33.69099 

31.16200 

31.36500 

30.44500 

29.97100 

31.42300 

30.30300 

29.75400 

30.59200 

29.68500 

33.51299 

32.93100 

32.86200 

254.13135 

183.81097 

155.08441 

92.20396 

63.46610 

45.13207 

40.93763 

35.94974 

38.84214 

33.20287 

27.71750 

25.19687 

22.05127 

20.39421 

24.06196 

21.27951 

19.35506 

18.64980 

16.95455 

26.13456 

24.23474 

24.13839 

O 

O 

O 

O 831U6 66 
CO 

OJ 

V 

Ol 

Z9 

QJ 

m 

- h 

c 

3 

O 

rt-

^ t 

<D 

3 

in 

-n 

>-H 

o 

cz 

PO 

m 

t n 

i 

-e» 

. 

o m 

3 " 

fD 

-5 

fD 

JD 

c 

mm. 

-s 

fD 

3 

QJ 

rtfD 

CX 

00 

1 

3 " 

o 

O . C 

a. OJ 

rtQJ 

O 

O 

3 

TD 

-^ fD 

r lfD 

3 

-5 

o 

o 

QJ 

.^. -T" 

JD 

c 

QJ 

mmi 

mm. 

rtfD << 

in 

in 

>. 

TD 

o 

. 

TD 

Oi 

-J 

QJ 

3 

(D 

rtfD 

-s in 

PO 

fD 

JD 

CJ 

- S 

CD 

Q . 

T3 

fD 

-5 

O 

fD 

3 

r+ 

Estimated Exceedance Rate 

10. 20. 30. 10. 50. 

_ j 1 1_ ^ 1 . 

+ 

+ 

TD 

O 

< 

fD •% 

in 

TD 

O 

t 

+

* ' x y 

Estimated Variance 

sP 50. 150. 250. 300. 

PO 

CD 

JD 

C 

_ j . 

-m 

fD 

CL 

-o 

(D 

- i 

o 

fD 3 

rtQJ 

S " 

.£> 

eat / i 

K > -

ca- 70. 8 cato 

+ 

+ 

+ 

+ 

+ 

+ 

++ 

03 

OJ 

-5 

TD 

O 

< 

(D 

in 

c 

in 

• a 

o 

+4 

+ + 
P.95 

TABLE 5-5. Estimated composite annual average 

of one-hour N02 , y(pQ ), and its estimated 

variance, Var(p(p0)), as functions of the 

required percentage, pQ , of data completeness. 

p o 

96.50000 

95.09999 

93.79999 

92.39999 

90.50000 

88.20000 

85.70000 

82.00000 

78.20000 

75.20000 

69.70000 

65.50000 

58.39999 

51.39999 

45.30000 

37.10000 

31.80000 

26.90000 

23.10000 

17.50000 

0.30000 

0.30000 

A 

•v(P0) 

33.44099 

34.75300 

41.21500 

44.69399 

44.60000 

44.89600 

44.36099 

44.71400 

45.92799 

45.62500 

44.94700 

44.63499 

45.01400 

44.45399 

43.99200 

44.30199 

44.23399 

43.65500 

43.77599 

43.56400 

43.36700 

43.24699 

A /•. 

Var(y(p0)) 

14.59878 

6.81544 

11.62743 

10.46519 

6.89711 

5.32763 

4.35964 

3.41974 

3.34880 

3.19546 

2.92656 

2.74853 

2.56410 

2.42393 

2.24541 

2.22146 

2.05578 

1.92911 

1.78446 

1.64982 

1.57100 

1.57575 

68 

63146 
(A) P(P 0 ) versus pQ (B) Var(p(pQ)) versus pQ 

CT) 

VO 

c 

to 

TD 

<V 

+-> 

ro 

E , 

+ + + + + + + + + + + ++++ 

+ 

+ 

10. 20. 30. 10. 50 . 60. 70 . 80 . 90 . 100. 

Required Percent 

in. 

CU -

o 

c 

ro 

• i — 

S_ 

rc 

TD ' 

OJ 

"b. 

+ 

+ 

+ + 

.+ + 

+ + + + + + + + 

" l — 1 1 1 1 1 

2 0 . 30. 10. 50 . 60 . 70 . 80 . 90 . 100. 

Required Percent 

FIGURE 5-5. Estimated long-term N02 air quality parameters 

as functions of the required data completeness, pQ. 

8 3 1 li C 
„P,97 

O 

0 

O 

TABLE 5-6. Estimated mean annual rate of 

maximum daily one-hour ozone exceedance, 

x(pQ ), and its estimated variance, Var(x(pn)), 

as functions of the required percentage, 

P Q , of data completeness. 

O 

O 

O 

O 

O 

p0 

98.89999 

97.79999 

96.39999 

95.29999 

93.89999 

92.89999 

91.50000 

90.09999 

88.29999 

86.89999 

85.00000 

83.00000 

80.29999 

78.00000 

73.79999 

67.20000 

61.20000 

50.89999 

40.00000 

26.80000 

3.70000 

0.40000 

x(P0) 

12.82300 

20.10000 

19.90200 

18.18400 

15.01800 

13.10600 

12.83500 

11.58300 

10.65700 

9.71900 

9.16100 

8.53400 

8.17600 

7.83800 

7.43100 

7.45200 

7.36400 

7.20000 

7.78500 

7.93800 

7.92600 

7.91500 

A A. 

Var(x(pQ)) 

30.20399 

22.94687 

15.08907 

10.62845 

7.01568 

5.04998 

4.27442 

3.30427 

2.67449 

2.18219 

1.82070 

1.53777 

1.32614 

1.16137 

1.01475 

1.00815 

0.89779 

0.82834 

0.93169 

0.95963 

0.96646 

0.96394 

o 

o 

83146 

o 

70 

B B ESD 
IZ 

Estimated Exceedance Rate 

Q) -T| 

i n i—i 

CT 

- h CRZ 

E 5 

unct 

-6. ions o rn 

-h in 

£ . 

JD 

C 

, -B 

fD 

*-— t3n 

. 60 

Perc 

fD • 

3 

r t- -j 

rrs 

CD QJ 

-5 

fD 

JD 

C 

- J , 

-5 

fD 

Q . 

Q . 

Ql 

rtQj 

r> O 

3 

TD 

m-1 

fD 

fD 

a. 

QJ 

3 

3 

C 

QJ 

— j 

O 

M 

o 

3 

fD 

OJ 

_ j . 

-5 

rt- JD 

fD 

3 

fD 

in 

in 

* 

TD 

O 

C 

QJ 

—' —,. rl-

<< 

TD 

QJ 

-5 

QJ 

3 

fD 

rtfD 

-s in 

20. 30. 

+ 

+ 

+ 

+ 

+ . 

•a 

o 

< 

a> 

m 

c 

m 

•O 

O 

"^ 

«?,-

Estimated Variance 

10. 20. 30. 

PO 

fD 

c 

-s 

fD 

CL 

TD 

fD 

- 5 

n fD 

3 

rtc r 10. 20. 3 s

-

j ^ 

S " 

cn 

^ j 

e a -

CD 

IS" 

UD 

s " 

T . 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ + 

^ . 

• 

+ + 

** r'

 w r * i 

+ + 

10. 

03 

at 

TD 

O 

< 

fD 

m 

o 
P.99 

LEAD 

Table 5-7 and Figure 5-7 show u(pn ), the estimated composite maximum 

V A A 

quarterly lead (pb) concentration, and Var (u(pn)), its estimated variA A ^ 

ance. Both y(pQ) and Var (p(pQ)) are fairly constant for pQ < 17 percent 

(corresponding to complete data recovery from a once-per-six-day sampling 

schedule). Only a small number of sites have more than 17 percent data 

completeness. The consequent rise in the estimated variance suggests that 

no more than a 17 percent data completeness requirement should be used to 

estimate the national status and trend in lead concentrations under 

current sampling requirements. 

TOTAL SUSPENDED PARTICULATE 

Table 5-8 and Figure 5-8 show y (Pg), the estimated composite annual 

average concentration of total suspended particulate (TSP), and its 

A A A 

estimated variance, Var (y(pQ)). The estimator u(pQ) is fairly constant 

with pg, and its estimated variance shows no pattern for increasing pg. 

Therefore, it is not possible to judge the effect that an increased data 

A 

completeness requirement would have on the variance of u(p0 ). 

Similarly, Table 5-9 and Figure 5-9 show no appreciable pattern in the 

A 

estimated variance of x(pQ ). Figure 5-9,. however, suggests that increasing the required percentage, pg, may induce a positive bias in the 

estimate of the national average exceedance rate. Because of this 

possible bias, it appears that no more than a 17 percent data completeness 

requirement should be used to estimate the national status and trend in 

TSP concentrations. 

SUMMARY 

In general, we see that the advantages of including a large number of 

sites in the estimation of national air quality trends outweighs the 

831*t6 9 7 2 
CO CO oo 

OOOOOOOOO*.^J0D0BOW*»Ut<7>a>^J-J^ J 

OOOO-'-^-'-'aiWOi'oMOONMWWWW W 

-»-*0ol-»{>.U10>-

,

^0l0l0a)lkCD-*00l0l0i D 

—« -« vl 00 -» Xk C7> -vl —'00lDO>O>lOWO>^l^fc-* —* —* 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

tk^VJIO>4>lkJ>,UII\)-»l\)-UUlUI~JO*ivJ-JOi-vir>ov O 

wo>WK)0)Wa)Mioa)a)K)N)i\)ON)Oo>-'a)ON) 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O ^ ^ i > Ul M - - —- — •OOOOOOOOO-' O — IS} Ul 

IS)UIVJ»04^Ul-J~Jis.^\JlCTivJOiCX5-JvO—• lO * -• Ul 

lOOlOWvllOMlOOOT.IWMJwUiOUI-'UllOWUi W 

TD 

O 

tr > 

.—. TD 

O 

V ^ < > OJ -5 . s -e > 

-.—.. TD 

O 

^ — • 

.—.-

o 

o 

3 TD — J fD r+ fD 

3 fD in in 

. 

O -h c+ 3 " fD -s fD JD 

c 

— J * -5 fD CX -a fD - j 

o fD 

3 

rt-

OJ UD fD V* TD 

CD >* 

O - h CL at 

r l-

OJ 

CD in 

rt-

—*. =1 QJ 

rt- fD a. 

< OJ 

~. mJ. 

Ql 

3 

O fD 

^ 

< a> 

-. mm-.s -c ^ V TD 

3 

- .. 

^-^ 

sr 

Ot in - h 

c 

3 ctions 

=

a 

X 

— 

=

c -

J

C

O

-

r

f

-

—

VC 

. 

C

> Q

Q 

V 

> x

*^

•o 

o 

%

t « Q

3 C

—

rt

tn 
frZ 

OJ 

I/I 

- t i 

c 

3 

O 

rtmJ. 

o 3 

i n 

O 

- h 

rt3 " 

fD 

n 

fD 

JD 

C 

—t . 

-i 

fD 

CX 

CL 

OJ 

r l -

Q) 

o 

o 3 

TD 

• ^ 

fD 

rtfD 

3 

fD 

i n 

i n 

,* 

TD 

o 

. 

- n 

i—, 

CT 

CZ 

PO 

m 

U l 

i 

v i 

. 

m 

in 

rt-

—*• 3 

QJ 

rtfD 

CX 

— J 

O 

3 

UD 

1 

rtfD 

1 

3 

—-J 

fD 

O i 

Q . 

Q l 

— i . 

-i 

JD 

C 

O i 

—' _ J . 

rtt < 

TD 

OJ 

-% Q) 

3 

fD 

fD 

-5 

l/l 

) . 0 0. 1 

pa 

tD 

JD 

c 

fD 

CX 

T3 

fD 

- 5 

o 

fD 

3 

0 . 2 

I 

Estimated Mean 

. 3 0. 1 0. 5 0, 6 0, 

4-

0 . 5 

_ i 

0 . 6 

I 

0.8 0.9 

i 

1.0 

4-

44-

4-

+ 

+ + 4-

• o 

o 

< 

ro 

- i 

m 

• o 

o 

^.00 0 

pa 

fD 

JD 

"5 

fD 

T3 

fD 

-s 

n 

fD 

3 

<-+ 

Estimated Variance 

0,002 0.00 1 0.006 ' 0.008 

-f^- I 4- + 

0.010 

0 0 

QJ 

4-

4-

O 

< 

fD 

T 

m 

c 

in 

• o 

o 

4- 4-

Z d 
P.3 

TABLE 5-8. Estimated composite average of 

24-hour TSP, y(p Q ), and it s estimated 

variance, Var(p(pn )), as functions of the 

required percentage, pn , of data completeness. 

p o 

30.68500 

16.71200 

16.71200 

16.43800 

16.43800 

16.16400 

15.89000 

15.89000 

15.61600 

15.34200 

15.06800 

14.79500 

14.24700 

13.97300 

12.87700 

12.05500 

10.41100 

8.21900 

6.02700 

3.83600 

0.27400 

0.27400 

y(P0) 

51.08199 

53.59200 

56.05000 

56.52800 

56.72600 

56.57700 

57.10899 

57.17000 

57.28699 

57.48299 

57.43500 

57.27200 

57.13400 

57.28999 

57.31000 

57.41100 

57.20499 

56.97199 

56.78200 

56.70900 

56.52599 

56.52599 

A A 

Var(y(pQ)) 

3.17931 

1.35666 

0.90998 

0.69591 

0.52045 

9.22075 

2.12360 

1.13663 

0.75091 

0.56889 

0.44404 

4.18731 

1.42163 

0.86768 

0.65872 

0.56483 

0.45510 

3.04553 

1.37268 

0.87463 

0.66334 

0.66334 

8311+6 75 
9Z 

fD ICI 

cx i 

rtQ . fD 

QJ "5 

rt- 3 

QJ 

—I 

o oo 

O -O 

3 

TD QJ 

(D - 5 

rtfD JD 

3 C 

fD OJ 

i n —' 

i n ->• 

» rtTD 

O TD 

• DJ 

-s 

QJ 

3 

fD 

rtfD 

-S 

in 

sP 

QJ -T I 

in i—i 

CT 

- h <= 

C PO 

3 m 

o r i - t n 

O CO 

3 • 

i n 

o n -ti m 

rtrt- -J3" 3 

fD OJ 

rt-

-i tD 

fD CL 

C —' 

_.. o 

po 

fD 

JD 

fD 

CX 

TJ 

fD 

(-> fD 

rt10. 

I 

20. 

i 

Estimated Mean 

30. 10. 50. 

L_ 

60. 

4-

4-

4-

4-

t 

} o 

< 

ro 

- j 

m 

o 

Estimated Variance 

3.0 1.0 5.0 6.0 7.0 8.0 9.0 10.0 

PO 

CD 

JD 

C 

-5 

fD 

CX 

TD 

fD 

- 5 

r> 

(D 

3 

-• =r J-

4-

4-

" 4-

++4- + 

i 

+ 

i 

4-

4-

1 1 ! 

4-

1 1 1 1 

4-

0 3 

OJ 

- J 

•o 

o 

< 

fD 

m 

T3 

O 

Vd 
-r

-

cn ^ J 

OOWOlOOMNWAJiiUIUIUIVIIUlOlOiOiOiOi O 

MMODOM^OO)iOM>IOOIOiO)(»-

> A<kNlNlO i 

Nl>JWNJ-»-'UlNlvl*,l00»lk-'l0l00M>H>l-'-» O 

JkAOlNjiO-*UlvlW-slUIO)N)(7ioOJi.a)a)NMU l 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O KJ IV) 

Ji^lwUIUIUIUIAWWWO.^UIOvJOivHO'OM* ' 

O5O)00WJkvliO4>.M— VIW(OWVIUIO>JMVJ1N)N ) 

~J-J-»UllO-»-JOiaDUlCT>Ulv000f>000O00OO-Nl 0 

OOOOOOOOOOOOOOOOOOOOO O 

OOOOOOOOOOOOOOOOOOOOO O 

0000000000000000--0000 0 

-.-.-.NjOl-'-* — — -* Ul O O -• — WMO -•-» N X". 

WWOBlOWOM-'-'(»7>IOia)-'(D^*>IOK)rv)U I 

OOOOO-'UHiOlOtO^lOWWWOiUlUIJikOOOW W 

-'fc^OO-'OOOWWOOlWfcOOiOOOvlMlOv J 

uwwooNio-JajNiuiwuiwww^ji.-'uioioia ) 

TD 

o TD 

o 

ai -5 TD 

o 

o 

o =? -a —J fD 

r l-

fD 

3 ro in in 

. 

o 

- t l 

rt-

3-

fD -5 fD JD 

c 

—i. -5 (D CX TD fD -5 

O fD 

3 

rt-

a> 

I Q fD 

mf TD 

O 

<* 

O 

fD in 

rt-

_*. =3 Q J 

rt- fD CX 

< CD -s 

m l . QJ 

3 

n fD 

» 

< at -s 

*-^ > 7 

*—~* TD 

O 

^—* 

^.^ ,. QJ in 

> 

> 

O 

CX - h 

Ol c 

rt- 3 

OJ n 

o 

3 in 
8Z 

QJ 

i n 

- t i 

c 

3 

<-> rt-

— i . 

O 

3 

i n 

O 

- t i 

r l -

ZT 

fD 

"1 

fD 

JD 

SZ 

mm. 

-s 

fD 

CX 

CX 

QJ 

rf 

QJ 

-n 

t—m 

CT 

c: 

PO 

m 

O l 

i 

VO 

• 

m 

in 

rt-

— j . 

3 

OJ 

rtfD 

CL 

i n 

3 " 

O 

-s r l -

r+ 

tD 

~l 

3 

O —I 

O CO 

3 "O 

TD 

—• QJ 

fD -"• 

r l - -5 

fD 

3 JD 

ID C 

in QJ 

i n —< 

TD v< 

o 

• TD 

QJ 

-5 

QJ 3 

fD 

rtfD 

~i 

i n 

X I 

fD 

JD 

fD 

Q L 

tD 

-5 

o 

fD 

3 

Estimated Exceedance Rate 

0.5 1.0 1.5 ' 2.0 • 2.5 

, , 4 1 . _ 

3.0 

+4-4 

4-

4-

f

++rt + -H- + TD 

O 

< 

fD 

-. 

m 

c 

m 

TD 

o 

Estimated 

.5 ' l.l 

Variance 

1.5 2.0 

PO 

fD 

-O 

c 

fD. 

CX 

T3 

fD 

-5 

O 

fD 

3 

TO 

4-

4-

4- 4-

"Hf 4-

o> 

• o 

o 

< 

m 

-. 

m 

c 

m 

"O 

o 

9 d 
P.7 

disadvantages of including sites having a low data completeness proportion, p.j. The broad inference from these 1981 results—which seems likely 

to hold even if more frequent intermittent samples or more complete data 

are required—is that the wide variability between sites argues against 

the exclusion of a large fraction of sites solely on the basis of data 

completeness requirements. Thus, the requirement for estimating national 

trends should be more relaxed than that for determining the attainment 

status of an individual site. 

83146 9 79 
REFERENCES 

P.8 

Abramowitz, M., and I. A. Stegun. 1964. Handbook of Mathematical 

Functions. Dover, New York. 

Brillinger, D. R. 1969. The calculation of cumulants via conditioning. 

Annals Inst. Stat. Math., 21:215-218. 

Breiman, L., 0. Gins, and C. J. Stone. 1978. "Statistical Analysis and 

Interpretation of Peak Air Pollution Measurements." Technology 

Service Corporation, Santa Monica, CA (TSC-PD-A190-10). 

Breiman, L., and C. J. Stone. 1983. "Approximate confidence bounds for 

extreme quantiles." In progress (will be submitted in 1984 to 

J. Amer. Stat. Assoc.). 

Cochran, W. G. 1963. Sampling Techniques, 2nd ed. Wiley, New York. 

Davidson, J. E., and P. K. Hopke., 1983. "Implications of Incomplete 

Sampling in a Statistical Form of the Ambient Air Quality Standard 

for Particulate Matter." Paper presented at the 76th Annual Meeting 

of the Air Pollution Control Association, Atlanta, Georgia. 

EPA. 1979. "Guideline for the Interpretation of Ozone Air Quality 

Standards." Environmental Protection Agency, Research Triangle Park, 

NC (EPA-450/4-79-003). 

EPA. 1981. "U.S. Environmental Protection Agency Intra-Agency Task Force 

Report on Air Quality Indicators." Environmental Protection Agency, 

Research Triangle Park, NC (EPA-450/4-81-015). 

Frank, N. H., and T. C. Curran. 1982. "Statistical Aspects of a 24-Hour 

National Ambient Air Quality Standard for Particulate Matter." Paper 

presented at the 75th Annual Meeting of the Air Pollution Control 

Association, New Orleans, Louisiana, June 1982. 

83146 6 80 
P.9 

Hajek, J. 1960. Limiting distributions in simple random sampling from a 

finite population. Pub. Math. Inst. Hungarian Acad. Sci., 5:361-374. 

Johnson, N. L., and S. Kotz. 1969. Discrete Distributions, Wiley, New 

York. 

83146 6 

81 
P.10 

Appendix A 

SULFUR DIOXIDE DESCRIPTIVE STATISTICS 

83146 5 
P.11 

APPENDIX A 

SULFUR DIOXIDE DESCRIPTIVE STATISTICS 

Using program NA282 of the EPA's Storage and Retrieval of Aerometric Data 

(SAROAD) system, various summary files of 1981 data were obtained. Figure 

A-l shows the distribution of data completeness among air monitoring 

stations. Table A-l shows the number of sites meeting an annual quota of 

75 percent data completeness and meeting a corresponding quarterly quota 

in a specified number of quarters. Table A-2 shows the same information 

for an annual quota of 50 percent. In considering whether a quarterly 

quota should be required so that an even distribution of data would be 

guaranteed throughout the seasons of the year, it is of interest to note 

that roughly 70 percent of the sites that met an annual quota also met a 

corresponding quarterly quota in each of the four quarters of 1981. 

There are low correlations between the percentage of data completeness and 

reported annual averages, standard deviations, and coefficients of 

variation (Figure A-2). 

Histograms of "expected" exceedances (i.e., X. or k./p. in the notation of 

Section 2) show that zero expected exceedances were reported by 96 percent 

of the stations (Figure A-3). The greatest number of expected exceedances 

were reported by the stations that had the most complete data (Figure A-4 

and Table A-3). 

83146 5 
p 

*- 3..' '-MMX 

•-VG ..'.'.-7.0.0/7 

* J . .'.'.O'.'.'.u.'J 

* i •:. u.fjijii 

* is..0.0.0../ 

* 1 b . <J!JXJ:J 

*z\ ytvrj:) 

* 2 4 ..O.O.'/.:' 

•••••2V ..'J.O.'.7J 

" 3«' . tJiUu:] 

* 3 3..'J.O:./.-/' 

* 3 G .lir,,",'! 

•"3'.' ..'J.O ,'J .a 

••''-12..07.'.•.•';/ 

* .15./,'.'.".•'.'J 

*!L I ..'/.O.o.-/ 

* 5 1 .ffO.XX 

:

'

;

 Ij <! . w.0.0.0 

*<y/ .#:•;<)-J 

'•':ufj. <:)£}:>.'• 

••uJ .ij:-J>:J:J 

* G 6 . ,0.<7.A./ 

••••'O'J..OX'A': J 

"Vi ! .!>!!:...) 

-'/{, .aU:.>:) 

••'•'/o ..00.0.7 

- 0 1 ..0.0 A J 

••Li'- ..'JX'.J..' 

*UV.0.O;J J 

«*,J.U'.A)'.O.'-J.,T 

*' J 3.0.0.0.'./ 

* 3 0 ..0.0.0.7 

*VJ,.U.0.J.7 

*l.'*:i.jL'.y.'/ 

•••l>..3..'J.o// 

'••l.ij:j..07.'./ 

* 111 . 0.0./ 

•••1 1 4 . JJJJ.-J 

5 Iff 

•AA A 

+ 

• AA 

+AAA 

+ 

•AAAAAAA 

• AA 

•AAAAAAAA 

•AAAA A 

• A/'A A A 

.|. *\/\/> \t\ 

•AAAAA A 

+,V\A.A,AA 

•AAAAAA A 

•AAAAA A 

• AAA 

• AAA 

•AAAAAA A 

• AAAA 

•AAAAAAAA A 

• AAAAAAAA. 

• AAA 

SYN20L COUNT MEAN ST.DEV . 

A 4/ 2 7G.98 7 23.05 8 

EACH SYMBOL REPRESENTS 1 OBSERVATIO N 

15 2ff 2 5 3.3 3 5 4 0 4 5 5 0 5 5 6 0 6 5 7 0 7 5 

+ A\ * AAAAAAAAAAAA 

:-AAAA AAAAAAAA 

+ AAAAAAAAAAAAAAAA 

• A,". A, A, A A A. A. A A A, A 

• A.AAAAAAAAAAA 

+A,'->AiA,AA.A.'\A1AAA 

+ AA.A.A.A,'\A.A,AAA,A 

• AAAA.AAAAAAAA 

•I'AA.A \ /\A.AAX '\J\* A 

•AAAAAAAAAAA / 

• A 

+ 

+ 

+ 

+ 

5 Iff 

AAAAAAAAA 

AAAAAAAAAAAAAA 

AAAAAA.AAAAAAAAAAAAAAAAAAAAAA 

AAAAA.".A.A.\AAA.AAA.\AAAAAA.AAAAAAAAAAAAAAAAA 

AAAAA.A.AAAA.AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

AAA.\A.-'.AAAA.A.AAAAA.AAAi.'-.A.AAAA.AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

..'.AAAAAAAAAAAAAAAAAAAAAAAAA 

15 2.0 ;:L> 3.7 3 5 4ff 4 5 Uff 5 5 6ff 6 5 7 0 V5 

8 0 

8ff 

FREQUENCY 

INT . 

3 

0 

2 

• 3 

0 

7 

2 

0 

5 

5 

5 

C 

6 

7 

6 

3 

3 

7 

4 

9 

9 

8 

3 

15 

12 

16 

2 1 

2 6 

4 0 

5 1 

6 6 

7 5 

3 8 

1 

0 

0 

0 

SS 

CUM . 

3 

3 

5 

8 

8 

15 

17 

2 5 

3 0 

3 5 

4 0 

4 6 

5 2 

5 9 

6 5 

6 0 

7 1 

7 8 

8 2 

9 1 

1 0 0 

1 0 8 

1 1 1 

1 2 6 

1 3 0 

.15 4 

1 7 5 

2 0 1 

2 4 1 

2 9 2 

3 5 8 

4 3 3 

4 7 1 

4 7 2 

4 7 2 

4 7 2 

4 7 2 

4 V 2 

PERCENTAGE 

INT . 

0 . 6 

0 . 0 

0 . 4 

0 . 6 

0 . 0 

1 . 5 

0 . 4 

1 . 7 

1 . 1 

1 . 1 

1 . 1 

1 . 3 

1 . 3 

1 . 5 

1 . 3 

0 . 6 

0 . 6 

1 . 5 

0 . 0 

1 . 9 

1 . 9 

1 . 7 

0 . 6 

3 . 2 

2 . 5 

3 . 4 

4 . 4 

5 . 5 

0 . 5 

10. 8 

14. 0 

15. 9 

8 . 1 

0 . 2 

ff.0 

0 . 0 

0 . 0 

0 . 0 

CUM . 

0 . 6 

0 . 6 

1 . 1 

1 . 7 

I . V 

3 . 2 

3 . 6 

5 . 3 

6 . 4 

V . 4 

0 . 5 

9 . V 

11. 0 

12. 5 

13. 8 

14. 4 

15. 0 

16. 5 

IV. 4 

19. 3 

21. 2 

22. 9 

23. 5 

26. 7 

29. 2 

32. 6 

3 V . 1 

42. 6 

5 1 . 1 

6 1 . 9 

75. 0 

91. 7 

99. 8 

100. 0 

100. 0 

1 BM. 0 

MJPl.il 

1 0 0 . 0 

Number of Sites 

8314 6 

FIGURE A-la. Histogram of the percentage, p, of data completeness for those sites 

reporting one-hour S0? concentrations in 1981: SLAMS (including NAMS). 
P.13 

LU . Si o — r-->^rioinoro'*co>.i-vcri^r^^tvi'Cco5r^iac7iL'i"vDi^u-)'5-inL-)S:5iS!5-v 5 

c s • • 

< z> mm —. re v L I o i ^ 5> >M •«• :•» cr.—• o a ,M n o s ? S) n a o & .-oi i 5 n Q ^ ^ ^ -r; > 5, 5, 5 . 

I—i_> _«_. « — csjMC\irererereri'<t^S''*Lio<Dor^r«.r».cociCiS!S;^S-,l 3 

Z m-,m-. rm mm m* 

U l 

0 • sSCM(\ji^oii^CT--i'^«ji-5»-^--'-p~-c-pj---<i--<i^oi."noioin-.r>vC\rvj!3L-)aj5).SiS i 

01 hUJ rr —i^-«~is~:3^c^ i — re—i^ifi^re^pj^N—•pjLT^'>*pj>j.'«3-reinr7«r7-LT5a^S7S5>. 

c >-• 

>- • 

o IE -«fcore o — oiM-^j-reS'r'-)— • —•r\j — ^-—• —ioin«t3 ,fM,'nr)Mfl^r)oioMi3ip?j > 

p • 

UJI — ^^ifliononivisiNroo^'-'SinrsSiiflBiosivMco- > ~ n n < n 5 ~« rsi 5) ^ ^ Ss 

O i 3 « — „ ^ mm mm (\jre-*-«-<-<~t\irereC^ J 

re 

•>? 

> • 5 

ixi y ~4 

• < 

I- > 

t/i zz 

LJ 

t/1 

CO 

re D 

cri 

z —, - . < • 

L J re 

C -•.• zz 

S) + 

CO 

i n + 

65 + 

un + 

UD 

ts + 

Ll + 

L l 

IT) + 

«a + 

•s> + 

re i 

+ 5) 

CO 

+ u l 

1^ 

+ S} 

+ LO 

+ i s 

o 

+ JT) 

L l 

+ B> 

+ LO 

•<*• 

+ s» 

O — re 

O 

:D 

^ + 

rj 

C2 CO. 

! CD ^ 

ca ca : 

+ + 

Crj 

CO 

ca 

a 

CO c3 

CJ CD 

ca a 

ca c i 

~2 C3 ca ca 

S - 1 3 J 

C3 13 C3 i i 

+ + + * 

a 

CD 

ca 

C3 

Hi 

C3 

CD 

CD 

CO 

ca 

cn 

ca 

ra 

a 

i :3 

: J3 

- -j^S 

\ ca 

! ca 

. 22 

CD : 

CD : 

' f~~ ; -

! Ca C i CD . 

! CD D3 CD : 

+ + + 

C3 

C3 

CJ 

CO 

C2 C3 

".: ca 

C3 C2 

c i -A) 

ca CJ 

+ + 

ca 

. •_) 

ca ca 

•Ai c3 

r i ca : 

ca ca : 

+ + 

ca : 

ca : 

ca . 

LD : 

ca : 

• CZ: : 

ca 

:a 

i ca ca ra ca a ca 

: ca ca 

+ + 

ca 

ca 

a 

ca ca 

t_': :a 

r.-- .a 

f-• ca 

c: za 

ca ca 

cc Cl' :J 

c ca _a 

CA : j ca 

ca CD cr. ca 

ca CD ca ca 

ca ::, A' -a ca 

•a .j- . . -a 

ca :a ca A , ca 

ca -.-J ca •_:• ;D 

ca CD ca c'i ca 

ca ::• CD CD ia 

CD CD CD CJ A "J 

ca ca ca ca ca ca ca ca 

cn ca ca ca ca CD ca ca 

• • - j c a :.-;: c 1 „ •-• C D 

JJC J jjcj ca La ca ca 

caracacacacracDr- a 

CD -Zi CD CD CA; JJ CCD ^ £3 

CD ca CD ca 

+ + + + + + + + + 

l / l 

CD 

tn 

l » -

o 

s~ 

OJ 

JD 

E 

3 

+ c l 

-\J 

+ -5 

CM 

+ ^5 

^ Ti ** .-, - ^ ' 

'. t ; 

:

-i s to 

Dr, • • 

• csi LT co —. -*• r^ ^a. re a: T . ~ J IJD 

re ' a n —< —• —1 ~ j ^4 .\i y re re :a -r -r *•:<•::*•: : ********* * * * i: * * 

'a? T> T: "D> 

re o T I ?T! LT ra -> 

D'JOT-SN' V 

• . ^ * • * •: < * * * • 

r> -7 & •:; ; 

:-> cS ~i Q r_ 

: a; -z; y. :•> 5 

•• a> c; Zs zz ^ 

r>- o re '.a r i 

"j Tl .71 T-I —• 

* * % rt * 

I/) 

OJ 

CU 

I 

Jcr 

cc 
P.M 

TABLE A-l. Number of sites recording one-hour SO2 during at least 

75 percent of the total possible hours in 1981, by the number of 

quarters in which at least 75 percent of the total hours per quarter 

are recorded. 

Site Type 

SLAMS 

(includes 

Other 

No 

wi 

75% 

NAMS) 

. of 

th 

Quarters 

at Least 

Hours/Quarter 

0 

1 

2 

3 

4 

Total 

0 

1 

2 

3 

4 

Total 

No. of 

> 75% 

0 

0 

17 

79 

244 

340 

0 

0 

8 

48 

117 

173 

Sites 

< 75% 

35 

48 

38 

9 

0 

130 

49 

89 

51 

59 

0 

248 

Total 

35 

48 

55 

88 

244 

470 

49 

89 

59 

107 

117 

421 

TABLE A-2. Number of sites recording one-hour S02 during at 

least 50 percent of the total possible hours in 1981, by the 

number of quarters in which at least 50 percent of the total 

possible hours per quarter are recorded. 

No. of Quarters 

with at Least 

Site Type 50% Hours/Quarter 

SLAMS 0 

(includes NAMS) 1 

2 

3 

4 

Total 

Other 0 

1 

2 

3 

4 

Total 

No. of 

> 50% 

0 

0 

19 

42 

343 

404 

0 

1 

15 

102 

170 

288 

Sites 

< 50% 

11 

30 

25 

0 

0 

66 

24 

48 

59 

2 

0 

133 

Total 

11 

30 

44 

42 

343 

470 

24 

49 

74 

104 

170 

421 

83m6 8 
SLAMS , SYMB0L=S 

OTHER , SYMBOL=0 

N= B91 

COR= .0500 

MEAN ST.DEV. REGRESSION LINE RCS.MS. 

X 70.490 25.722 X= ..0SG73*Y+ 68.328 6G3.36 

Y 25.025 16.639 Y= .03S13*X+ 22.479 276.30 

CM 

o tn 

to 

tu 

140.0 + 

122.5 + 

105.0 + 

37.50 + 

70.03 + 

52.50 + 

35.0.7 + 

,X+. .+* 

0 s 

so 

o b 

oo 

so OS 

0 

OS 

SOS S 

0 

0 s 

s 

ss s 

s 

s 

17.50 + 

.0 0 

000 

00* 

0 ss 0 OS 

S S O S 0 

0 0 0 0* 00 

0 SS 00 s 

s ssoooo 

0 0 0' ss *s 

* * SSS OS OS s 

<

;0 0 s * 

SO 0 S 0 s SOOO 0 OOOO* 0 sssoooo*** 

00 oos 000 SOO *00 S 00 OS 000*00 

SS OOOSOOOOOOO 0 OS OSO 0*S 0 

OS SO 00 OS-00 00 00 *0 0 0*0 

*s 

so 

. 0 

OS 

0 

000 S 

0*SS*0 

00*0000 

SSSSS00*-00S 

* "000 *s 

SS* S0"0*'0 

SOO OOS 0 SSS 00* 0 s 0 00 s 

S OS 

s 

S 0 

S * 

*0 

0 ss 

*S OS 

S *0 

0*S*** 

s* ** 

*S0*00 

S OS** 

S OS s 

o*s*oo 

0 0* 

s 

s 

s 

0 

s 

SOS 

SSSSO 

oss 

s s 

ssssos 

ss s* 

s*s*ss 

s*ss«* 

3 * * * S * 

s**s** 

2 * * -J; rt ft 

rt rt * ft rt rt 

Q * * * * rt 

s**sos 

SS**;--

s*s *s 

ft rt 

ss 

SSS 

ssos 

rtrtrtQ 

rtrtrt 

«*(-)* 

rt rt rt rt 

--s*s 

0** 0 

*s*s 

•••so 

*sss 

SSS 

0.00.0" + 

7.50 0 22.5 0 37.5 0 52.5 0 67.5 0 82.5 0 97.5 0 

.00 0 15.0 0 30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 105. 0 

FIGURE A-2a. Scatter plots of ambient air quality statistics versus 

the percent of total hours recorded in 1981 (p) by sites reporting 

one-hour S09 : annual average versus p. 

83146 
SLAMS , SYMBOL-S 

OTHER , SYMBOL=0 

N= 091 

COR= .071 2 

liEAN 

X 7.v. 4 93 

Y 39.20 5 

ST.DEV. REGRESSION LINE RES.II*;. 

25.732 X= . .vriU.j-3-Y+ 6U.223 GG2.07 

32.733 Y= ..J9.J33:-X+ 32.-^5 10G7.6 

4.C0 

300 

2 00 

17-.7 

1.C0 

0.0 

, + .X, ,+* 

0 

s 

0 

s 

0 

00 

00 

so 

0 

s s 

0 0 

OS ss 

0 so 

00* 0 OS 

s 0 OS 

S 0 0 OS 

* 0 0 00 

00 030 0 

003 S S'.'O-'O 0 0 

*••* S SO 0 'Ju'-* 

-OS 00•-•••SS 00'"'S 

iOO so ••• 30•••• 00 

0 SO 00b005-0* 

00 OSS 0 

0 

s s 

00 s 

rt 

0 0 0 s 

S 0 

0 00 0 S 0 

O'-OSS s 

OSOOOO 0 ••OS * 

ss s soaoo--*osu3*o 

UO 00S SS03S00 «rto f t O•*-"*" ••• 

OSOO ;,; OOS*S ;*S*" ; MMVSO3SS ;'*" 3 

•'SS *0 * S--SS 00 = 3-3S3**S--* 

00 * *0S0 00 OSO 0 * *0"-v 

0 0 0 S 0 0 0 

....+....+....X....+....+ . 

s 

S O S 

0 0 s 

ooso 

0 OS * 

OS o •••"oo 

*SSSS*0** 

s*s***s*o 

rt rt rt rt rt rt rt rt rt 

* * rt * rt .-.. 2 * S 

rt rt s rt rt rt s rt s 

••• * s * •'••• s s s 

s soss s 

7.50.J 22. C.J 37.50 52.50 67.50 32.50 97.50 

700 15.00 30.00 4 5.00 60.00 7 5.00 90.00 105.0 

FIGURE A-2b. Standard deviation versus p. 

831«+6 
P.17 

SLAMS 

OTHER 

N- 891 

COR=-.0330 

SYMBOL-S 

SYMBOL-0 

MEAN ST.DEV. REGRESSION LIME RES.MS 

X 70.493 25.732 X=-1.0273*Y+ 72.173 664.72 

Y 1.G345 .32749 Y = - ..O0l'u6*X+ 1.7091 .60477 

,+....+....+....+ ,X+. , +* 

12.0 + 

C 7 . 3.0 

0 

E 

r, rrr 

3.00 

1 .30 

0.00 

.0 

.S 

.s 

OOOO 

O S O 

000 

-OSOO s 

OS 

so 

00 

0 

OS 

OS 

O-O 

o*. 

0 

0 000 

OS 0 

OS 0* 

0 0 

0 0 00 

so 

SOO 

00 

0 

so 

0 s 

0 0 OS 

0 S O 0*00 S 

0 0 000 * 

s o*o-o •••* * 

00 SOO-O--* -'S 

00 0 0003 *000'000* 

0 000 

00 0 

SO 0 

0 0 SOO 

OS 0 03 00 

S"30ft*0-'0S0 0**00 00SS 

0 

0 

s s 

s s 

00 0 

00 s 

o * ••* 

* s * * 

rt* « n 

*0'- 0 s 

* =•• o •• s s • 

SOS 

0 

OS o••••'• SO'-1 

O S O 

0 

•• a 

0 

' S SSS'-S 

sss *••••• so 

s 

0**3** 0 —S 

!

'.:0"S*S*3 3 

0 S3 SSS* 

00 

. + . .X . + . . 

<0*S! 

*sy 

s * 0 

0 0 * s 

00 ••*o 

**os* o s 

S*0 "* 0 

rtrtrt rt .'.- rtrtrt 

* rt * rt rt rt rt rt 

* rt rt rt rt rt s s 

s** 

S3*00 0 

7.500 22.50 37.50 52.50 G7.50 02.50 97.50 

0.000 15.00 30..'70 4 5.00 60.00 7 5.00 90.00 105.0 

FIGURE A-2c. Coefficient- of variation (standard deviation divided by annual 

average) versus p. 

831U6 
SYMBOL COUNT MEAN ST.DEV. 

A 4 7/T 0.180 1.424 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

E 

* . 6.00000 

* 1 . 2.0000 

* 1 . 0.70.00 

* 2 . 4000 0 

*3 .000r''7' 

*3 . G 00.00 

* 4.20000 

* 4. 0000.7 

* 5 . 1.0.000 

* 6..7.7000 

*6.G0O>)'0' 

ft).20000 

* 7.8.0000 

*3 . 4.0000. 

«-j . i.r rr or/7 

* 9. G 0.0 0 0 

M 0.2.00.7 

-10.0000 

')1.1000 

* 1 2 . 000.7 

'M2.G000 

* 13 . 2.0.00 

*13.i.!0('.'T 

'Ml . 4.0.00 

« 1 5 linf'AI 

•M 5. 6'"''• 7 

'•••• 1 fi . 205 0 

"1G.'V".0 

* 1 y. -I' • '7 

* 1 0 0'-;..a 

*Ki. G0O.7 

* 19.2007 

*19 .8""' J 

"20.4070 

rt ,> 1 . 000/ 7 

*;:i .G.0.00 

ft22.2000 

-2 2.0000 

•/".AA 

•AAA 

•AAA 

• AA 

• A 

+ 

i-

+ 

+ 

+' 

+ 

+ 

•A 

• A 

• A A 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1-

+ 

+ 

+ 

+ 

+ 

•A 

+ 

+ 

• A 

+ 

+ 

» • 

+ 

+ 

+ 

5 

"•.AAA/ 

\ 

< 

5 

1 0 15 

.AAAAAAAAA. 

1 0 15 

2 0 25 

-.AAAAAAAA/' 

2 0 2 5 

3 0 

rt rt A Ar t 

3 0 

3 5 4 0 45 5 0 5 5 

AAAAAAAAAAAAAAAAAAAAAAAA A 

3 5 4 0 45 5 0 55 

G 0 6 5 7 0 7 5 8 0 

AAAAAAAAAAAAAAAAAAAAAAA * 

G 0 6 5 7 0 7 5 8 0 

FREQUENC Y 

INT . 

4 5 4 

4 

3 

2 

1 

0 

13 

0 

0 

0 

0-

0 

1 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

CUM . 

4 54 

4 5 8 

4 61 

4 6 3 

4 6 4 

4 G 4 

4 G 4 

4 04 

4 6 4 

4 0 4 

4 6 4 

4 64 

4 0 5 

4 GG 

4-6 fl 

4 03 

4 G'J 

4 0 n 

4 6 Cl 

4 G 3 

4 0 0 

4 0 (.1 

40 0 

40 0 

4 Oil 

4 00 

4 •'.'> 0 

4 li 0 

4 0 9 

4 0 9 

4 0 9 

4 7 0 

4 /0 

4 7 0 

4 7 0 

4 7 0 

4 7 0 

4 7 0 

PERCENTAG E 

INT . 

96. 6 

0 .9 

0 . 6 

0 .4 

0 .2 

0 . 0 

0 . 0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0 .2 

0 .2 

0 .4 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 .2 

0 . 0 

0 . 0 

0 .2 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

CUM . 

96. 6 

97. 4 

9 8 . 1 

9S . 5 

90. 7 

98. 7 

98. 7 

98. 7 

98. 7 

98. 7 

90. 7 

98. 7 

90. 9 

9 9 . 1 

99. 6 

99. 6 

99. 6 

9 9 . G 

99. G 

99. 6 

99. 6 

9 9 . G 

9 9 . G 

9 9.6 

99. 6 

9 9 . G 

"99.6 

9 9 . G 

99. 8 

99. 8 

99. 0 

100. 0 

1 0 0 .0 

100. 0 

10.7.0 

1 0 0 .0 

1 0 0 . 0 

170. 0 

Number of Sites 

FIGURE A-3a. Histogram of expected primary exceedances, E, for those sites 

reporting 24-hour S02 concentrations in 1981: SLAMS (including NAMS). 

8311+6 
EAC 

E 

+ — 

* . G7.70.7.7 +8IJO 

M .20.0.7.7 + BBBC 

* 1 .0.0.000 +0o 

* 2 . 4.7.0.0.7 H-B 

* 3..0.7.0.0.7 + BB 

* 3 . 0.0.7.0.7 +B 

* 4 . 2 00.7.7 +B 

* 4 .09.00.7 + 

* 5 . 4.7.00.7 + 

* 6.007.70 + 

"6. 0 0.0 077 + 

* 7 .2 0 07 7 + 

* 7 . 07.700 + 

A ;j _ .; .1.1000 + 

'•'9.000.0.7 • 

* 9 . 60.0.70 + 

* 10 . 20.0.7 + 

"10.0.70 0 + 

* 1 1 .4 00 0 + 

'12.0.0.0.7 + 

* 12 . 0.0.0.7 + 

-13.270 0 + 

* 13 . 0.0 0.7 + 

* 14 . 4 0.077 + 

* 15.0.0.0.7 + 

* 15 . 600.7 + 

* 1 G . 2.0.00 + 

* 1 G . 0.70.7 + 

* 17 . 4.7.00 + 

* 1 0 . .70' 7.7 + 

* 18 . 0 0 0.7 + 

* I 9 . 200.7 + 

* 19.000 0 + • 

* 20.400 0 + 

* 21.000 0 + 

*21.000 0 + 

* 2 2 . 2.00.7 • 

* 22.000 0 + 

10 1! 

GUB. 

BEiB'. 

L- '3 B £ C u 3 0 b 

UB 

SY 

H SYMBOL 

20 25 

+ ;- -

'_> 0 L> B 3 0 U •. 

M70 L 

B 

Kill* RE 

3.7 

.- ;. — 

• Li 7 2' 3 B 

COUNT 

42 1 

SLUTS 

MEAN ST.DEV. 

0.067 0.300 

1 OBSERVATION 

35 40 45 50 55 60 65 70 75 80 

f (. ( 1 1 1. + + + + 

BliUBBBCBUBBBBBbDBBBBBU'BBBBBUB BBBBBBBBBBBBBBBBBBBB* 

FREQUENCY 

INT. CUM. 

404 

10 

2 

1 

2 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

404 

414 

416 

417 

419 

4 20 

421 

421 

421 

4.21 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

421 

42 1 

4 2 1 

42 1 

421 

421 

421 

421 

421 

421 

421 

421 

421 

PERCE 

INT. 

96.0 

2.4 

0.5 

0.2 

0.5 

0.2 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

NTACE 

CUM. 

.0 

.3 

,C 

.0 

.5 

,8 

.0 

.0 

.0 

96. 

99, 

93. 

99. 

99. 

99, 

100. 

100. 

100. 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100. 

100. 

100. 

100 . 

100.0 

10.7.0 

10.7.0 

100 .0 

1.70.0 

100 .0 

1.7.7.0 

. 0 

. 0 

. 0 

. 0 

100:1 

1.70 

100 

1 0.7 

1.7.7 

10.7.0 

1.7.7 . .7 

10.7.0 

17 7. 0 

1.70 . 0 

1.7.7.0 

1 0.7.0 

• - + -

10 20 

-- + -

25 30 35 

-- + -

40 

-- + -

45 

-- + -

60 

•- + -

70 

-- + 

80 

8314 6 

FIGURE A-3b. Other sites 
P.20 

SLAMS , SYMB0L=S 

OTHER , SYMBOL=0 

N= 891 

COR= .0813 

MEAN 

70.457 

.12632 

ST.DEV. 

26.073 

1.0679 

REGRESSION LIME RES.MS. 

X= 1.S335-Y+ 70.275 675.55 

Y= .00333*X-.10730 1.134 1 

, +* 

10.0 + 

17.5 + 

15.0 + 

12.5 

10.0 

7.5.7 

5.J7 

2 . 50 + 

OO 

SOOS 

0 

ss 

0 

0 

s s 

0*00* 0** 

0.0 0 y*Q*Q****QQ**Q##*rt**********"***Q*****Q rt** *rtrt***rtrt*rt***:7 . rt rtrtrt 7": rtrtrtrtg y 

+....+....+....+....+....+....+....+....+....+X...+....+....+....•....+ . 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

.000 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE A-4. Scatter plot of the percentage of total days recorded in 1981 

(p) versus expected primary exceedances (E) for sites reporting 24-hour 

average SCL83H*6 
TABLE A-3. Frequency table of the percentage of total days recorded in 1981 (p) 

versus expected primary exceedances (E) for sites reporting 24-hour S0? averages. 

0 - 28X 43X <» 60.5X <- 71X B0X 86X 89.5X 93X <-. 96.5X <- 100X TOTAL 

0 

1 

<- 2 

<- 5 

<» 20 

TOTAL 

87 

0 

0 

1 

0 

88 

88 

fi 

0 

0 

0 

88 

87 

0 

0 

0 

0 

87 

88 

0 

1 (SI 0 

91 

88 

0 

5 

1 

0 

94 

81 7 

0 

1 

1 

0 

83 

77 

fj 

3 

0 

0 

80 

87 

2 

1 

0 

3 

93 

102 

2 

0 (SI 2 

108 

73 

A 

1 

0 

1 

79 

8S8 

e 

12 

7 

6 

891 

83146 

TJ 
Appendix B 

CARBON MONOXIDE DESCRIPTIVE STATISTICS 

83146 5 
P.23 

APPENDIX B 

CARBON MONOXIDE DESCRIPTIVE STATISTICS 

Figure B-l shows the distribution of data completeness among air monitoring stations. Table B-l shows the number of sites meeting an annual quota 

of 75 percent data completeness and meeting a corresponding quarterly 

quota in a specified number of quarters. Table B-2 shows the same 

information for an annual quota of 50 percent. Note that roughly 80 

percent of the sites that met an annual quota also met a corresponding 

quarterly quota in each of the four quarters of 1981. 

Histograms of one-hour expected exceedances show that zero expected 

exceedances were reported by 99 percent of the stations (Figure B-2). A 

scatter diagram of expected one-hour exceedance versus percentage of data 

completeness is presented in Figure B-3 and in Table B-3. 

Similar information for overlapping eight-hour exceedances is shown in 

Figures B-4 and B-5 and in Table B-4. More than 10 percent of the 

stations reported over 40 expected exceedances. No clear relationship 

between expected exceedances and data completeness is evident. 

83146 5 
p 

* 3 . .7000.7 

*G .007''0 7 

* 9 . JU.70.7. J 

*12.0.7A' 0 

* 1 5 . 00/7.7 

* 18..07.'.'..f 

* 2 1 .0.7.'A > 

*'AA' ..707.7 

* 2 7..'J,.,.0 

*.J0.0..'0. 7 

* 3 3 . f)'ifi 7 

* 3 G ..'10.0.7 

* 3 9.0. 0 0.7 

* 4 2 .0.00 0 

* 4 5 .07.'.0.0 

* 4 o .0.00.7 

* 5 1 0"'i-"---T 

* 5 4 ..U.'J.O.'J 

*57.00.0. 7 

* 3 0 . 00.0.7 

* 63..0.0.O.J 

* GO.'0.7.70 

*G 0.0.0 . ::•..' 

* 7 2.000 . j 

* 7 5 . 0'fy.7 

* 7 0.0.0:7; J 

* 3 17 07707 

* U 4 .000. 7 

* 0/ . .0.7.0.7 

*y.O..'/.«'.'0.7 

"•03.0. 7 0.7 

* JO . 0 0.0.7 

*9'J..O.o.O0 

* 102..0.0.7 

* 1.05 .0.0 0 

*lu:j . "'' 7 

* 1 11 .00 0 

* 1 1 4.00 0 

SYMBOL COUNT 

A 36 4 

EACH SYMBOL REPRESENTS 

5 1 0 15 2 0 2 5 3.7 3 5 4 0 

• A 

• A 

• A' 

• AAA 

• AA 

+ 

• A 

•AAAAAA 

• v ^ \AAA A 

• AAA 

• AA 

• AAAA 

• AAA 

•AAAAAA A 

• A 

•AAAAA A 

• AA 

+ 

'•AAA 

• AAA 

• AAA *\AAA 

•AAAAAAA A 

+ A.\A.'\A.\AAA 

•AAAAA A 

• AA.". AAA AAAAA 

•AAAAAA A 

•AAAAA A 

•AAAAAAAAAAAAAA A 

+ AAAAA;*- 1 "• \<VAA 

+AAAAAA.V.AAAAAAAAAAAAAAAAAAAAAAAAA 

MEAN ST.DEV . 

80.35 3 23.13 8 

1 OBSERVATION 

4 5 5 0 5 5 6 0 6 5 7 0 

• AAAAA/-..'.AAAAAAA.V\AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

' + A A AAA;'.,'.. \AAA. '.AAA AAA AA A A AA AA A AA AA AA A A A AAA Afi .AAAAAAAAAAAAAAAAAAAA 

• A."..'.. '.AAA. :-.AAAAAAAA,a ,'^AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+ 

+ 

+ 

+ 

5 10 15 2 0 2 5 3 0 3 5 4 0 4 5 5 0 5 5 6 0 6 5 7 0 

7 5 

7 5 

8 0 

8 0 

FREQUENCY 

INT . 

1 

1 

1 

3 

' 2 

0 

1 

6 

8 

3 

2 

4 

3 

7 . 

1 

6 

2 

0 

3 

3 

7 

8 

9 

6 

11 

7 

6 

15 

12 

3 3 

4 7 

6 1 

7 0 

15 

0 

0 

0 

0 

CUM . 

1 

2 

3 

6 

8 

8 

9 

15 

2 3 

2 6 

2 8 

3 2 

3 5 

4 2 

4 3 

4 9 

5 1 

5 1 

5 4 

5 7 

6 4 

7 2 

8 1 

8 7 

9 0 

1 0 6 

11 1 

1 2 6 

1 3 8 

1 7 1 

2 1 8 

2 7 9 

3 4 9 

3 6 4 

3 6 4 

3 6 4 

3G4 

3 6 4 

PERCENTAGE 

INT . 

0 . 3 

0 . 3 

0 . 3 

0 . 8 

0 . 5 

0 . 0 

0 . 3 

1 . 6 

2 . 2 

0 . 8 

0 . 5 

1 . 1 

0 . 8 

1 . 9 

0 . 3 

1 . 6 

0 . 5 

0 . 0 

0 . 8 

0 . 8 

1 . 9 

2 . 2 

2 . 5 

1 . 6 

3 . 0 

1 . 9 

1 . 6 

4 . 1 

3 . 3 

9 . 1 ' 

12. 9 

16. 8 

19. 2 

4 . 1 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

CUM . 

0 . 3 

0 . 5 

0 . 8 

1 . 6 

2 . 2 

2 . 2 

2 . 5 

4 . 1 

6 . 3 

7 . 1 

7 . 7 

8 . 8 

9 . 6 

11. 5 

11 . 8 

13. 5 

14. 0 

14. 0 

14. 8 

15. 7 

17. 6 

19. 8 

22. 3 

23. 9 

26. 9 

28. 8 

30. 5 

34. 6 

37. 9 

47. 0 

59. 9 

76. 6 

95. 9 

100. 0 

100. 0 

100. 0 

100. 0 

100. 0 

Number of Sites 

83146 

FIGURE B-la. Histogram of the percentage, p, of data completeness 

for those sites reporting CO in 1981: SLAMS (including NAMS). 
• • 

p 

* 3.00000 

* 6 . .0.7.700 

* 9 .07.000 

* 12 . 00.00 

* 1 5 .0.0.7.7 

* 1 o . .7.0.70 

* 2 1 .0.0.70 

* 24 . 7''-'"'f 

* 2 7.00.007 

* 30.0.0.0.7 

* 3 3..0.7:70 

* 3 0.00:70 

* 3 9 f •"' '> '} 

* 4 2 . 0.0.7:7 

* 4 5.0.0 7.7 

.*4 0..7r77.7 

* 5 1 .007.7 

* 5 4.00.0.7 

* 5 7.00.7.7 

*G 0.0.0.7.7 

* 63.00.077 

* 6 6.007.7 

* 6 9.007.7 

.-* 72.0.0.70 

*75 nrc-'/j 

* 7 U . 07.0.7 

* 3 1 .0.00.7 

•34 .0007 

* 3 7.0.07) .7 

* 90.0.7.0.7 

* 9 3.0.0.0.7 

* 9 G . 00'00 

* 9 9 .000.7 

* 102.0.00 

* 1.0 5..0770 

•108.000 

•111. 00.7 

* 1 1 4 . .0.70 

83146 

• • 

EACH SY 

5 10 15 20 

.•Lit} 

+ CISB 

+ CB 

+ B 

+ L;B 

• UL, 

• LiL.nC 

+ .v:. i',:.-.33 

• OOB B 

+ U0C 

+ !',:L:I!D7 

• Coll!" 

+ 17 0 07', B 

• UOL B 

• CB 

+ 

+ B 

+ 

+ B 

+ 

+ 03B 

+ G 

• lA'OB 

+BBBB B 

• CB 

+ 

+I;L.[;B 

• UL.B 

+ B 

+ D7,;.733BBBBBB 

• CO.'UB^iBBBBBBB 

+CLBBB7 B 

• BBLBBiiBB 

+ 

+ 

+ 

+ 

+ 

5 10 15 20 

------

• • • 

SYMBOL COUNT MEAN ST.DEV. 

B 116 58.952 31.631 

•IBOL REPRESENTS 1 OBSERVATION 

25 30 35 40 45 50 55 60 

25 30 35 40 45 50 55 60 

Number of Sites 

FIGURE B-lb. Other sites. 

--..-._ — ... — . . _ . — 

65 

65 

. - . 

• 

70 

70 

— . . 

75 

75 

• 

80 

80 

. 

% 

FREQUENCY 

INT. 

2 

3 

2 

,1 

2 

2 

4 

6 

4 

3 

6 

4 

5 

4 

2 

0 

1 

0 

1 

0 

3 

1 

4 

5 

2 

0 

4 

3 

1 

12 

14 

7 

8 

0 

0 

0 

0 

0 

CUM. 

2 

5 

7 

8 

10 

12 

16 

22 

26 

2 9 

35 

39 

44 

48 

50 

50 

51 

51 

52 

52 

55 

56 

60 

65 

67 

67 

71 

74 

75 

87 

101 

108 

116 

116 

116 

116 

1 16 

116 

• • 

PERCENTAGE 

INT. 

1.7 

2.6 

1 .7 

0.9 

1 .7 

1 .7 

3.4 

5.2 

3.4 

2.6 

5.2 

3.4 

4.3 

3.4 

1 .7 

0.0 

0.9 

0.0 

0.9 

0.0 

2.6 

0.9 

3.4 

4.3 

1 .7 

0.0 

3.4 

2.6 

0.9 

10.3 

12. 1 

6.0 

6.9 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

1 .7 

4.3 

6.0 

6.9 

8.6 

10.3 

13.8 

19.0 

22.4 

25.0 

30.2 

33.6 

37.9 

41 .4 

43.1 

43. 1 

44.0 

44.0 

44.8 

44.8 

47.4 

48.3 

51.7 

56.0 

57.8 

57.8 

61 .2 

63.8 

64.7 

75.0 

87. 1 

93. 1 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 
P.26 

TABLE B-l. Number of sites recording CO during at least 75 percent 

of the total possible hours in 1981, by the number of quarters in 

which at least 75 percent of the total hours per quarter 

are recorded. 

Site Type 

SLAMS 

(inc 

Other 

:ludes 

No. 

wi 

75% 

NAMS) 

of Quarters 

th at Least 

Hours/Quarter 

0 

1 

2 

3 

4 

Total 

0 

1 

2 

3 

4 

Total 

N o. of 

> 75% 

0 

0 

5 

53 

212 

"270" 

0 

0 

0 

17 

33 

"50" 

Sites 

< 75 % 

15 

37 

33 

9 

0 

"94 

17 

35 

8 

6 

0 

"66" 

Total 

15 

37 

38 

62 

212 

"364 

17 

35 

8 

23 

33 

Tl? 

TABLE B-2. Number of sites recording CO during at least 50 percent of 

the total possible hours in 1981, by the number of quarters in which 

at least 50 percent of the total hours per quarter are recorded. 

Site Type 

SLAMS 

(includes 

Other 

No. of 

with 

50% 

NAMS) 

Quarters 

at Least 

Hours/Quarter 

0 

1 

2 

3 

4 

Total 

0 

1 

2 

3 

4 

Total 

N o. of Sites 

> 50% 

0 

0 

10 

42 

261 

"313" 

0 

0 

1 

16 

49 

1J6" 

< 50 % 

7 

23 

20 

1 

0 

Jl 

9 

30 

11 

0 

0 

"50 

Total 

7 

23 

30 

43 

261 

364" 

9 

30 

12 

16. 

49 

116 

831U6 8 
• • 

E 

* . 3.7.7000 

'•'' . 00.000 0 

•

:

 . 9.00.000 

* 1 . 200.0.7 

* 1 . 50.0.70 

* 1 .170.0.00 

•'••2 . 1.0.000 

•••2. 4.70.00 

ft2. /.O.O70 

- 3 .0.0.000 

••'•'3. 30.0.7.7 

•••3.''0.0 0 0 

" 3 .90.00.7 

-4.2000/ 7 

v

-'4 .!"' 00 7 <~f 

* 4 . li .0.0.0.7 

*5.1000 0 

* 5 . 4 00.0.7 

* 5 . 70.0.0.7 

*6.0000 0 

* 6 . 30.00.7 

* 6 . 6 00.0.7 

* 6 . 9.0.000 

* 7 . 2.0.0.7.7 

* 7 . 5000. 7 

* 7 .30.00 0 

''••Xi. 101''f-f 7 

* 3 . 4 0.00.7 

* 3 . 700.0 0 

* 7 . fiiin.00 

* 9 . 3000 0 

* 9 . 600.0 0 

* 9 . 90.00 0 

* 10 . 2.0.00 

* 10 . 5.000 

* 10 . 8.000 

*11.100 0 

* 1 1 . 400 0 

8314 6 

- - - - -

• •••• • 

SYliBOL COUNT MEAN ST . DEV . 

A .7, 4 0..7G 0 .7.09 1 

EACH SYliUOL REPRESENTS 1 OBSERVATIO N 

5 1 0 15 2 0 2 5 3.7 3 5 4 0 45 5.7 5 5 G0 6 5 7 0 7 5 8 0 

•AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA-AAAAAAAAAAAAAAAAAAAAAAAAA 

+ 

+ 

+ 

• A 

+ 

+ A 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ • 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

• A 

+ 

+ 

+ 

+ 

+ 

+ 

• A 

+ 

+ 

•4-

+ 

5 1 0 15 2 0 25 3 0 3 5 4 0 46 5 0 5 5 6 0 6 6 7 0 7 5 8 0 

Number of Sites 

FIGURE B-2a. Histoqram of expected primary exceedances, E, for those 

sites reporting one-hour CO concentrations in 1981: SLAMS (including NAMS) 

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

• 

FREQUENCY 

INT . 

* 30 0 

0 

0 

0 

1 

0 

l 

0 

0 

0 

0 

0 

0 

H 

0 

0 

ll 

kl 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

.7 

0 

0 

0-

^^^^ M 

CUM . 

3 0 0 

3!;.7 

"6 <J I'J 

3 0 0 

3 6 1 

3 0 1 

UCi 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 6 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 0 2 

3 7 2 

3 0 2 

3 0 2 

3 0 3 

3 0 3 

.36 3 

3 0 3 

3 0 7 

3 0 3 

3 6 3 

3 0 4 

3 0 4 

3 0 4 

3 0 4 

3G4 

-

• 

PERCENTAGE 

INT . 

9U. 9 

0 . 0 

0 . 0 

0 . 0 

0 . 3 

0 . 0 

(J . U 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 3 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 3 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

CUM . 

9 8. 9 

9 0 . 9 

9 3. 9 

90. 9 

9 9 . 2 

9 9. 2 

0 0 . D 

9 9. 5 

9 9. 6 

9 9. 5 

9 9. 5 

99. 5 

9 9 . 5 

9 9 . 5 

9 9. 5 

9 9. 5 

99. 5 

9 9. 5 

9 9. 5 

99. 5 

9 9. 5 

9 9. 5 

9 9. 5 

9 9. 5 

9 0. 5 

9 9. 5 

9 9. 7 

9 9 . 7 

9 9. 7 

9 9 . 7 

' 9 9 . 7 

9 9. 7 

99. 7 

1.7.7.0 

10.7. 0 

1.70. 0 

1.0.0 . 0 

100. 0 

• 
300000 

600000 

900000 

20000 

50000 

80000 

10000 

40000 

70000 

00000 

30000 

60000 

90000 

20000 

50000 

80000 

10000 

40000 

70000 

6.00000 

6.30000 

60000 

90000 

20000 

50000 

80000 

10000 

*8.40000 

* 8.70000 

*9.00000 

* 9.30000 

* 9.60000 

•9.90000 

•10.2000 

•10.5000 

•10.8000 

•11.1000 

•11.4000 

* 1. 

*1 . 

•1 . 

• 2, 

• 2, 

• 2. 

•3. 

•3. 

•3. 

•3. 

• 4. 

• 4. 

• 4. 

•5, 

•5. 

•5. 

SYMBOL COUNT 

OTHER B 116 

EACH SYMBOL REPRESENTS 

MEAN ST.DEV. 

0.094 1.012 

1 OBSERVATIONS 

10 

•- + -

15 20 25 30 35 

-- + -

40 

-- + -

45 

-- + -

50 

-- + -

60 

-- + -

65 

-- + -

70 

-- + -

75 80 

-- + 

FREQUENCY 

INT. CUM. 

+ + -

+BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBUBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

+ 

• 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

•- + -

10 

•- + -

15 

-- + -

20 

-- + -

25 

-- + -

30 

-- + -

35 

-- + -

40 45 50 

-- + -

55 

-- + -

60 

-- + -

65 

-- + -

70 

-- + -

75 

-- + 

80 

Number of Sites 

115 

0 

0 

0 

• 0 

0 

0 

0 

0 

.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

l 

0 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

115 

1 15 

115 

1 15 

115 

115 

115 

115 

1 15 

115 

1 15 

1 16 

116 

PERCENTAGE 

INT. CUM. 

99. 1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.9 

0.0 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99.1 

99. 1 

99.1 

99. 1 

99.1 

99. 1 

99. 1 

99. 1 

99. 1 

99.1 

99. 1 

99. 1 

99.1 

99. 1 

99. 1 

99.1 

99. 1 

99.1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

99. 1 

100.0 

100.0 

FIGURE B-2b. Other sites 

83146 
P.29 

SLAMS 

OTHER 

N- 480 

COR=-.0068 

MEAN 

X 75.234 

Y .06792 

SYMBOL=S 

SYMBOL=0 

ST.DEV. 

27.015 

.77943 

REGRESSION LINE RES.MS. 

X=-.23588*Y+ 75.250 731.29 

Y=-196E-6*X+ .08269 .60876 

, + . .X . + . 

12.0 + 

10.5 + 

9.00 + 

7.50 

6.00 

4.50 

3.00 + 

1.50 + 

0.00 Y*0 0** S**0 0****** "*0*S***0* SS* • SSS***S*****SSS**************S Y 

.+....+....+....+....+....+....+....+....+....+....X . ... + ... . + ... . + ... . + . 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

0.000 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE B-3. Scatter plot of the number of expected primary exceedances 

(E) versus the percentage (p) of total hours recorded in 1981 for sites 

reporting one-hour CO concentrations. 

8 311*6 
TABLE B-3. Frequency table of the percentage of total days recorded in 1981 (p) 

versus expected primary exceedances (E) for sites reporting 1-hour CO averages. 

E 

0 

<- 11 

0 -• 26X 

46 

1 

<= 46.5X 

48 

0 

<=• 71% 

50 

1 

<- 83X 

46 

0 

<- 90% 

65 

1 

P 

O 91.5X 

28 

1 

<- 94X 

50 

1 

<- 96X 

48 

0 

<- 97.5X <- 100% 

43 51 1 

0 0 1 

i 

TOTAL 

475 

TOTAL 47 48 51 46 66 29 51 48 43 51 I 480 

83146 
EACH 

SYMBOL COUNT 

A 3G4 

SYMBOL REPRESENTS 

MEAN ST.DEV. 

34.301 112.449 

1 OBSERVATION 

•40.0000 

•80.0000 

•120.000 

•160.000 

•200.000 

•240.000 

•280.000 

•320.000 

•360.000 

•400.000 

•440.000 

•480.000 

•520.000 

•560.000 

•600.000 

•640.000 

•680.000 

•720.000 

•760.000 

•800.000 

•840.000 

•880.000 

•9 20.000 

•960.000 

•1000.00 

•1040.00 

•1080.00 

•1120.00 

•1 160.00 

•1200.00 

•1240.00 

•1280.00 

•13 20.00 

•1360.00 

•1400.00 

•1440.00 

•1480.00 

•1520.00 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+ + + +• + + + ^ + * (. ^ + + + 1. + 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

•AAAAAAAAAAAAAAAAAAA 

•AAAAAAA 

•AAAAAAAA 

• AAA 

AAAA 

A 

AAA 

AA 

AA 

AAA 

FREQUENCY PERCENTAGE 

INT. CUM. INT. CUM. 

+ 

+ 

• A 

10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

309 

19 

7 

a 3 

4 

1 

3 

2 

2 

3 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

3.09 

328 

335 

343 

346 

350 

351 

354 

356 

358 

361 

3G1 

361 

361 

362 

362 

362 

362 

362 

362 

362 

363 

363 

363 

363 

363 

303 

363 

363 

36 3 

363 

363 

363 

363 

364 

364 

364 

364 

84.9 

5.2 

1 .9 

2.2 

0.8 

1 . 1 

0.3 

0.8 

0.5 

0.5 

0.8 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

• 84 

90 

92 

94 

95 

96 

96 

97 

97 

98 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

9 9 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

100 

100 

100 

100 

9 

1 

0 

2 

1 

2 

4 

3 

8 

4 

2 

2 

2 

2 

5 

5 

5 

5 

5 

5 

5 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

0 

0 

0 

0 

Number of Sites 

83146 

FIGURE B-4a. Histogram of expected primary exceedances, E, for those 

sites reporting 8-hour CO concentrations in 1981: SLAMS (including NAMS) 
E 

•40.0000 

•30.0000 

•120.000 

•160.000 

•200.000 

•2 40.000 

•280.000 

•320.000 

•3G0.000 

•400.000 

•440.000 

•480.000 

•520.000 

•560.000 

•600.000 

•640.000 

•600.000 

•720.000 

•760.000 

•800.000 

•840.000 

•880.000 

•920.000 

•960.000 

•1000.00 

•1040.00 

•1000.00 

•1120.00 

•11G0.00 

•1200.00 

•1240.00 

•1200.00 

•1320.00 

•1360.00 

•1400.00 

•1440.00 

•1480.00 

•1520.00 

SYMBOL COUNT MEAN ST.DEV. 

OTHER B 116 28.348 90.300 

EACH SYMBOL REPRESENTS 1 OBSERVATIONS 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

• BBBB 

+ B 

+ B 

+ B 

+ 

+ B 

+ B 

+ B 

+ B 

+ B 

+ 

+ 

+ 

+ B 

+ 

+ -

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

FREQUENCY 

INT. 

103 

4 

• 1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CUM. 

103 

107 

108 

109 

110 

1 10 

11 1 

112 

113 

114 

115 

1 15 

1 15 

115 

116 

116 

116 

116 

1 16 

116 

116 

116 

116 

116 

1 16 

116 

1 16 

1 16 

116 

1 16 

116 

116 

1 16 

116 

116 

116 

116 

116 

PERCENTAGE 

INT. 

88.8 

3.4 

0.9 

0.9 

0.9 

0.0 

0.9 

0.9 

0.9 

0.9 

0.9 

0.0 

0.0 

0.0 

0.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

83.8 

92.2 

93.1 

94.0 

94.8 

94.8 

95.7 

96.6 

97.4 

98.3 

99.1 

99. 1 

99. 1 

99. 1 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100 .0 

100.0 

Number of Sites 

FIGURE B-4b. Other sites , 

83146 
SLAMS 

OTHER 

SYMBOL=S 

SYMBOL-0 

1400 + 

1225 + 

1050 

875. 

700. 

525. 

350. 

175. 

0.00 

N= 480 

COR=-.0376 

MEAN 

X 75.482 

Y 32.862 

OS" 

ST.DEV. 

27.260 

107.46 

,+X. 

REGRESSION LINE RES.MS 

X=-.00953*Y+ 75.795 743.63 

Y=-.14812*X+ 44.043 11555 

Y S O 

+*0 0** S*00 00****S0«*** *S*SO*SS 

0 s 

s 

s 

s 

0 s 

s 

s 

O S s 

ss 

s s ss 

s ss ss 

S OS ss 

o s*s** s**ss 

s s s 

ssoss s s 

osss**os********s**s*s*********s 

i-1... + .... + ....+X.. . + .... + .... + ..• 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

0.00 0 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE B-5. Scatter plot of the number of expected primary exceedances 

(E) versus the percentage (p) of total hours recorded in 1981 for sites 

reporting 8-hour CO concentrations. 

83146 
TABLE B-4. Frequency table of the percentage of total days recorded in 1981 (p) 

versus expected primary exceedances (E) for sites reporting 8-hour CO averages. 

0 25.5 46X <= 70.5X <- 82.5X <- 89X <- 91.5X <- 94.5X <- 96.5X <- 98.5X <- 100X TOTAL 

0 

1 - 5 

<- 10 

<=• 50 

<- 100 

<- 1400 

TOTAL 

36 

0 

1 

3 

0 

6 

46 

31 

a 5 

7 

1 

5 

49 

31 

3 

3 

7 

3 

3 

50 

29 

3 

6 

6 

1 

3 

48 

27 

5 

4 

5 

1 

4 

46 

31 

1 

3 

9 

0 

2 

46 

26 

7 

5 

4 

2 

6 

50 

21 

6 

5 

12 

0 

1 

45 

25 

S 

0 

15 

3 

4 

52 

21 

2 

1 

11 

4 

9 

48 

278 

32 

33 

79 

15 

43 

4 80 

83146 
Appendix C 

NITROGEN DIOXIDE DESCRIPTIVE STATISTICS 

83146 5 
P.36 

APPENDIX C 

NITROGEN DIOXIDE DESCRIPTIVE STATISTICS 

Figure C-1 shows the distribution of data completeness among air monitoring stations. There are low correlations between the percentage of data 

completeness and reported annual averages, standard deviations, and 

coefficients of variation (Figure C-2). 

83146 5 
UJ • IO "3 Si Cl Cl Ol •« s3 CO (^ CSI CSJ VO mm 55 LO -C Cl Cl -CO & LNI O Sl 3 ! Cl r~ — •«• CM LO rv. LO 3 : Sl 3 ! Si S3 

w s: 

< z> ^ —> csi csi csl CM ^ vo Ci 5: —, co LO o c i S; csi -o vo r^ c i ,-NJ o ro mm m _-. ,-v. t-n _ _ CM co ff: 5; S; S; Si 

t - O rt rt rt rt rt rt CM CM .NJ CSJ ,\J ,SJ 01 CO CO KT * * LO LO VO CS CO Cl Cl 3 ; vv Si 3 ; 31 

2 mm rt m« mm m* 

LU 

<J • IO LO Sl Sl Sl 5 } IJ ^ •<» LO LO 51 ^ L l 01 LO Sl LO ^ LO Sj • « ^ LO CJ1 Cl Cl ^ C70 CO CSI CM CO LO Sl 31 Si ^ 

oil -

UJ I I 5 ) ^ — -H^ ^ — CSJCSI—i.SlCSJCSISCSJ — CSl->CSJrt(SICS|i3-—iCSI^fLOLOvOrv^—.OSiSlSlG ^ 

Q_ ,—• .rt »—* 

>- • 

C> Z TM csi "=r vo vo vo o> -=• cn M co >. "0 ro c i ,-sj to c i •>=• ^ —, o LO CO ^ t r VO r-» 3? vo ^ Sj -c LO _O LO LO LO 

z 3) -< *- csi csi csi co co co •*? •«• •=? _o ^ j o \ , s a ; , 5 « , o ^ o CA 3- 'i> A- 3' 3- : ^ 

L I U rt - « TM - , « N N ,-j Cv cs; CM 

=> 

o • 

UJ >— rt — CM CSI S l ^ CO LO LO CO -> < • LO rt vo r o •>? r o LO CO f LO C l CO VO 5 } CSI —i CO VO - i C l •>» mi ^ 3 i 5 § 

CJ Z rt rt rt — —, Cv CM —• 

•_) - I 

t l 

(S + 

CO 

i n + 

ca + 

rv. 

IT) + 

vo 

c^ 

CO 

u>j -- • : 

o ts: 

. H 

to 

CM 

r. 

z -•• 

< • UJ —• 

s : rv 

3 

—, — '

< 

> 

UJ 

in 

CO 

-1 

—7 

;-c 

o 

L O 

L O 

• o 

L O 

:o 

•^. 

fe) + 

+ S3 

CO 

+ Sl 

+ LO 

VO 

+ 53 

o 

+ LO 

LO 

+ 3 

LO 

D L 

t f 

+ 

*—; 

" 

- • 

~ 

'0 

< 

+ 

-

: 

-

• ? 

+f <* 

< < 

+ + + 

^ s :- » '- . 

• - : -" : " • 

..- , > • ..' • - . ' 

-** > — ••-* . 

• . 0 -"5 

• 1 — . T-K 

< < < 

< < < < 

. <<<< ^. < 

<: <<< < < 

<<<<<<< 

+ + + + + + + + 

"^^-?>-.>7*i>rN>-^"7^ ^ 

--: •-! o ~-\ '-: •": ~i •' 

_ i ^ .. • - • •. ' • . ^ . • 

*~» '-) '"- ~: *"• T: .-~" — 

H . . i 

::• — -.:- r- r- n r= n 

—t : J ~ J :- J v") : o :•) -:o 

< 

+ 

^ 

• T i 

"*; 

C s i 

-r 

< <r 

<r* 

* • ' 

<:•• " 

*=r 

+ 

^ 

:~-

_0 

- 7 

<rr' 

<• 

+ 

7 ^ 

7; 

— 7 

:c-

•-

< 

<<r" < < r <•• < < - < 

•"-• -T.-' ^ * -•• '.-' 

< < < t ' < 

+ + + + + 

3- — - ^ &7 >-. 

• - * " ^ •".** '."•> T * 

- ! % — —• .— -• — '- ' TM- -

..1 ,-) -.' J'J 

< <r 

< — • 

f * 

•^ + 

— 

-"> 

. • ^ 

•i.. 

J 

AAA 

<r 

<< 

<" *-«•:•; " 

- r ' 

+ 

^-* 

• - ! 

,~; 

T-. 

? 

cr 

<r 

<T 

<r <r 

T ' <r* 

< < 

+ + 

^ > ^ 

•~\ T-: 

1 

- ^ •- * 

r-J :.") 

:^ "^ 

< 

< 

< < 

< 

c 

<r 

<.' <r' 

*r'' 

e ^ ' 

+ 

^ 

— ."

:D 

"*-

< 

< • -

t< f 

< • 

"' r ' 

+ 

• - S 

0 

" • 

—1 

'•} 

* * • 

t f 

-f 

t^ 

r < 

<r 

« r * 

c 

< 

•=•: " 

+ 

- s 

- ' 

-. 

. 

- • • 

:D 

AAA 

<+* < 

cf 

t f 

^ 

r. ' 

t f 

t : ' 

<. <•' 

-.* 

+ 

-^, 

:; 

J- ' - 1 

. 

>. :o 

<* ^ 

<< 

c - r 

t-.' 

<< <r 

< •»"• 

<* 

+ 

*-% 

^ A. 1 

,~! 

r^ 

•'\ 

< 

< 

< 

< 

< 

<r 

<, <? ' 

t f 

t f 

t f 

^ <f 

t f 

< r J 

<f 

< 

+ 

?-> 

Q 

y 

' - s 

. 

-^ " • • 

< 

< 

< 

<:•' 

•rf 

< 

<f 

< 

<<f < 

< 

<f 

•f 

<C* 

<C»* 

•c' 

T 

T * 

< 

+ 

^ 

."• 

r> 

• j ; 

7 i 

t f ^AA 

<r 

<r 

< • 

< 

<* <*• 

+ 

J - > 

o 

"3 

y . 

•XT-

< 

+ 

>^ 

C; 

t 

r.' 

- , • 

—« 

+ 

— 

c> 

_ * * :f) 

z.-

—• 

+ + + 

•^ ^ ^ 

^> T i . 3 

5 r^

 :I5 

• • * -O ^ • * 

Q ^_ ^^ 

— . *- H — 

+ 5 

co 

+ LO 

rM 

t/i 

OJ 

0 0 

+ si o 

ai 

JD 

+ LO E 

CO 3 

+ S3 

OO 

c 

4 J 

S -

o 

CL 

OJ 

sl/l 

<u 

4-> 

a* 

1/1 

o 

^3 

4-> 

So 

M -

v/1 

(/l 

CU 

c 

OJ 

r— tn 

o-s: 

E <C 

o z 

o 

cn 

m c 

+ J -rfO X I 

TD 3 

4- 13 

o c 

OJ e t 

Oi_ l 

IO 0 0 

c 

(IJ . . 

U 1— 

5- 0 0 

QJ CTi 

Q.r — 

cu c 

MO 

t/1 

c 

0 

• r* 

E ••-> 

<a 

s_ 

ro 

S_ 

on-> 

O 

4 -i 

1/1 

•r~ 

X 

, 

c 

OJ 

0 

c 

0 

0 

CM 

to O 

'~ 

0 

L U 

OC 

ZD 

O 

k r t 

L L . 

z 

s3 

O 

s: 

1 

OJ 

c 0 

a *: x * # 
cr 

X I m 

o 

i 

- i 

CD 

l/l 

•- »- t 3: 3. vo c. to t . --: •:- — c: f. r: vo c c: 3. v j ;- rt c; to r J vr c. vet s 

:

 • L IL-;.- ?•'? A ? A ^ A A J: lv;? % ? A 5:? S ;:;??;?§:: ir ? 

— ^ c: co 

~ >r ~ J-' 

r-7 . » . 

"' -

C

 :~ • 

+ + + + + + + 

Si + 

ro 

i a + 

CJ 

6; + 

to 

CO + 

fO 

O l,f 

tn 

—'• t 

cf t s + 

fD 

(/i 

61 + 

C l 

CO + 

S; + 

u i + 

CO 

6; + 

CO 

co 

ra 

cv BBBBBB c~ B 

— 

CO 

ra 

CO 

CT 

+ + + + + . 

ra co c; co — tr. 

ra cc: co — 

P C D 

C3 C3 

CtJ 

+ + + + + + + + + + + + + + + 

— (o - ct cs c; ct1 — rararac-racoractt! 

C3 CO CV Ctt Ct: :;-: ~ CC Ct. Ct CC. C" C5 to Ct 1 

co cc- ~ c= ct cc- co cp ct ct; 

ct; — CO ra CO w CO CO 

Cd C" Cti ra 03 Cti CO 

CJ c n a n 

ra ra ra 

ra CJ 

+ + + + + + 

co rcra 

ra 

+ >-3 

ro C 

Lt r- t t 

r t r. 

- - O 

CT 

+ 6 } 

to 3 

rc n 

• > 

63 C to 

co 

to 

LO > 

L". —: ro ro 

• - -o: ro 

S S; ^ ^ 63 co — co to rt rt rs; co to to ro ro ro to ~ LO cr to rt LO rc ro cc v ; to c: c i ro to 6; ro -- to —; n 

• O 

. • e 

„ « « rt « rt rt 1 - (-. n 

rt rt « rt rt rt rt 3; o vo o vo vo cc cc v j v j v : v j c- c~ c- to to t-" to to to to ro ro rt •- cr 2 

vo vo vo vo vo o to to c; to to ^ iz; LO rs- cc cr. to ro vo c: to -vj to to c; c i .». o vo LO v j rt vo c i c i to to . o; O 

• • < 

•-< 0 3 

6: 6; 63 6; i5 to o to to 6; 65 « *>. ro to rt rt rt ^ 6; to co ro ^ to « rt o co to rn co ~ ro s ; rt s ; ro — r.i 

—.7-

ia i s tS i s i s rc rs1 ro to c: co v j rs1 to to -vi v j v j to c: ro ir; LO ix ts- v i v j v j vo to v i i s v ; to 6; VM co to . o 

m 

fe &. i^ is tt; ir. '^r co c? vj vi v: vi v! ci o o ci C7 c- to to to to to to to to to ro ro »- OH 

is is is is ,'r:

 :i~- co ci rv vo cc v; LO — cc <.o v.o ro ic; cc v: r; v' co to <a cc- v: lo.- to >— to vo vi co co to ro c > 

0; '.n 

ts&isisiasjoccritoisrotocitovotovoisjcoia'—ovOtocotovi3;toto&:toroci7Sts;toc o . m 

sed 
SLAMS , 5YMD0L=5 

OTi-IER , SYMDOL=0 

N = 324 

COR= .092 4 

MEAN ST.DEV. REGRESSION LINE RLS.KM. 

X 64.23 1 2 S . t Z 7 X= .11430-Y + 5U.27 5 7 3 ! . J u 

Y 43.24 7 22.53 7 Y= . .0744 3:-;<+ 33.43 5 50/.3 3 

1 40. 0 

ta. 

CSJ 

O 

c 

ro 

OJ 

7 0 . 

,+X. 

0 0 

s 

s 

0 0 

0 os s 

SO 0 

003 00 

0 003 

S O S 

; . O.S 0 0 

>-s s 

0 0 00 s 

s so 

s ss 

' s 

0 

S S S 

0 

0 

SO OS 

OS 

S 

I 7.5.0' + 

0.00 0 

. 0 

. 0 

+ 

00 

s 

00 s 

00 

.00 

0 0 

OS s 

0 * 

0 

00 

0 

-0 

0 

0 

S3 

s S S 0 

S S 3 3 OS 

s s 

SS S3 O S SS 

0 SSS S SOS 

s s-s s s* 

s o o ss:; s-- s 

S S 3 3 S3•• 3 VJ 

S S OS 0 S3 S3333 ; 

S S S 3 S * 3 3 3 3 

S SOSO S S SS 5 -'-3 

0 0 SS S 

S3 0003 

0 

3 SOS S 

0 

0 

3

s

 03333 

0 3 00 

0 03 S 

000 -

0 

s 

7.50.7 22.5.7 37.5 5 52.5 0 67.5 0 82.5 0 97.5 0 

0.00 0 15.0 0 30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 \Z. 

FIGURE C-2a. Scatter plots of ambient air quality statistics versus 

the percentage of total hours recorded in 1981 (p) by sites reporting 

N 02 : annual average versus p. 
SLAMS 

OTHER 

, 5YMDOL=S 

, SYMBOL=0 

N= 324 

COP.= .004 4 

MEAN ST.DEV. 

X 64.231 23.027 

Y 27.535 13.335 

REGRESSION LINE RE S.MS. 

X= ...-.I';'J 2* Y+ 63.335 737.36 

Y= .0.32 19 •••;:+ 2/..145 193.39 

SS 

OO 

+ 0 

.0 

70 

0 

0 s o oos 

0 0 0 so s 

0 0 OS S 03 

0*0 * 0 3* 

CO S SOOS 3 3 

os;o * o os 

s S 0 s 

0 00 S 0 00 

00 0 

0 . 00 0 

3 0 0 

s 

so 

OS 

SS 

0 

S OS 

30 

S O S 

s s 

0 ••S * 0 

0 0-- 0 

S 000S 

3 

S c r* c r* 

JO o o 

0 SS 5 S SSO S 

0 0 333S '••3 30 

S SSS 333333-30SS 

S S O S 33303U33S 

0> o O o 0) J IJ j _, j, .. ._. j 

SO 00 3 3 S3 JJ 3-A:;

S 3" SO 03S J' 

3 S 

0 0 

0-00 

3"0. 3 

30 00 

0 0 

7.300 

..•....•....•....•....•..•.+....+..X.+... . 

22.30 37.20 52.50 67.50 n? nc; 

13.00 30.00 45.00 60.00 75.00 

UL.J S 97.50 

90.00 105.0 

FIGURE C-2b. Standard deviation versus p. 
P.41 

SLAMS , SYMDOL=S 

Oi/liiR' , 3Yi-iBOL=0 

N= 324 

COR--.143 0 

iiEAN ST.DEV. REGRESSION LINE RES.HI, 

•' 6-1.231 23.027. X = -2.3 . 1 39* Y + 77,;06 /IJ.I. 

Y ,J/30 4 .20351 -Y = -..:.J103*X + .74277 .0414C 

X.. + . 

0 

0 

s 

0 

0 

0 

0 0 

0 

OJ 0 

JO s 

3 S 30033 0 

0 0 30 003 3 - J - 0 

3 S O S3 

* 3 0 0 S3 

0 0 3 0 0 

0 ' 0 

0 0 S S 3 S S 

SOOO 0 0 0 *3S S S O 

0 0 S * S 0* 0 03S 3 3 S*0 5 

0 SJOS * 3 30 3 33 SS 330 033 3 

S SSO 0 * S 3 33 SS 3"3"3-;-'-"--t00 

3 - 0 SS O SSOS 3 3-0'""33 • M J33 3-:S-• •• 33 

S 0J3 0 S 0 - ViSJ : .3-3-33S S 

3 0 3 SO 03 30 33--303S 

S 3 SO 

0 

1 .y. 

0 . 

. . + ... . + ... . + ... . + ... . + ... . + ... . + ... . + ..X . + . .. . + . . . . + .... + . 

!2.3 0 37.3 0 52.3 0 67.3 0 32.3 0 9/.3 0 

30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 105. 0 

FIGURE C-2c. Coefficient of variation (standard deviation 

divided by annual average) versus p. 
P.42 

Appendix D 

OZONE DESCRIPTIVE STATISTICS 

831A6 5 
P.43 

APPENDIX D 

OZONE DESCRIPTIVE STATISTICS 

Figure D-1 shows the distribution of data completeness among air monitoring stations. 

Histograms of "expected" exceedances (i.e., x. or k./p. in the notation of 

Section 2) show that roughly 85 percent of the stations reported seven or 

fewer expected exceedances in 1981 (Figure D-2). (These statistics are 

purely for descriptive and comparative.purposes and do not directly 

pertain to the attainment of the ambient standard, since this is 

determined from three years of data, not from a single year.) Somewhat 

greater expected exceedances were reported by stations having nearly 

complete data; 4 percent of the stations reported zero expected exceedances (Figure D-3 and Table D-1). Judging by Table D-1 and the low 

correlation given in Figure D-3, no clear relationship is evident between 

expected exceedancs and data completeness. 

83146 5 
*3.00000 

* 6.00000 

* 9.00000 

*12 

*15 

*10 

*21 

*24 

*27 

*30 

*33 

*36 

*39 

*42 

*45 

*48 

*51 

*54 

*57 

*G0 

*63 

*66 

*69 

*72 

•75 

*78 

*81 

*34 

*37 

*90 

*93 

*96 

*99 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

00.00 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

. 0000 

* 102.000 

* 105.000 

* 108.000 

* 1l1.000 

* 114.000 

SYMBOL COUNT 

A 548 

EACH SYMBOL REPRESENTS 

MEAN ST.DEV. 

81.767 18.528 

1 OBSERVATION 

10 15 10 25 30 35 40 45 50 55 60 65 70 75 80 

•AAAAAAA 

• A 

• A 

•AAA 

• AAAA 

• A 

• AA 

• AA 

• AAA 

• AAA 

•AAAAAAAAA 

• AAA 

•AAAAAAAA 

•AAAAAAA 

•AAAAAAAAA 

•AAAAA 

•AAAAAAAAAA 

•AAAAAAAAAAAAAA 

•AAAAAAAAA 

•AAAAAAAAAAA 

+AAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAA 

10 

— »•-

15 20 25 30 35 

-- + -

40 50 

- - + -

55 60 

• - + -

65 

• - + -

70 

-- + -

75 

- - + 

80 

FREQUENCY PERCENTAGE 

INT. CUM. INT. CUM. 

0 

0 

0 

0 

7 

1 

1 

3 

4 

1 

Z 

2 

3 

3 

9 

3 

8 

7 

9 

5 

10 

14 

9 

11 

12 

26 

27 

31 

55 

51 

62 

73 

78 

21 

0 

0 

0 

0 

0 

0 

0 

0 

7 

8 

9 

12 

16 

17 

19 

21 -

24 

27 

36 

39 

47 

54 

63 

68 

78 

92 

101 

1 12 

124 

150 

177 

208 

263 

314 

376 

449 

527 

548 

548 

540 

5 40 

5 48 

0.0 

0.0 

0.0 

0.0 

1 .3 

0.2 

0.2 

0.5 

0.7 

0. 

0. 

0. 

0. 

0, 

1 . 

0. 

1 . 

1 . 

1 . 

0. 

1 . 

2. 

1 . 

2. 

2. 

4, 

4. 

5. 

10.0 

9.3 

11.3 

13.3 

14.2 

3.8 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 .3 

1 , 

1 , 

2, 

2 

3 

3. 

3 

4 

4 

6 

7, 

8 

9 

11 

12.4 

14.2 

16.8 

18.4 

20.4 

22.6 

4 

3 

38.0 

4 3.0 

27. 

32. 

57.3 

68.6 

81". 9 

96.2 

100.0 

100.0 

100.0 

100.0 

100.0 

Number of Sites 

H ? ] i,(, 

FIGURE D-la. Histogram of the percentage, p, of data completeness for those sites 

reporting 0-, concentrations in 1981: SLAMS (including NAMS). 
SYMBOL COUNT 

B 175 

EACH SYMGOL REPRESENTS 

MEAN ST.DEV. 

64.906 27.499 

1 OBSERVATION 

10 1! 

*3.00000 

*6.00000 

* 9.00000 

* 12.0000 

* 15.0000 

•18.000.9 

*2 1 .0000 

*24.0000 

* 2 7.0000 

* 30.0000 

•33.0000 

'36. 

'39. 

'42. 

'45, 

'48. 

'51 . 

•57. 

*G0. 

*63. 

•66. 

•09. 

•72. 

•75. 

•73. 

•81 . 

•84. 

•87. 

*90. 

•93. 

•96. 

•99. 

0000 

0000 

0000 

0000 

0000 

00.00 

•54.0000 

00 00 

0000 

0000 

0000 

.0000 

0000 

0000 

0.000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

'102.000 

'105.000 

1108.000 

'111 . 

•114'. 

+ + + •» 

+ B 

•BB 

+ B 

+ 

+BBBBB 

+ BBB 

• ED 

• BB 

• BBBB 

+BBBCBB 

•BBBBBBBB 

• BBBB 

+ B3B 

•BBBBBB 

• BB 

+BBBBB 

• BBBB 

+ B 

• BB 

+ B R 

•BBBBBBBB 

+BBBBB 

• BB 

•BBBBBBBB 

•BBBBBBBB 

+BBBBB 

•BBBBBBBBBB 

•BBBBBBB 

•BBBBBB 

•BBBBBBBBBBB 

•BBBBBBBBBBBBBB 

+GBB3BBBBBBBBB 

+0BBBBGB3BB 

+BBBBB 

20 

- - + -

25 30 35 40 45 

-- + -

50 55 60 65 70 75 

-- + -

80 

000 

000 

10 

• - + -

15 

-- + -

20 30 

-- + -

35 

•- + -

40 

-- + -

45 

-- + -

50 60 

-- + -

65 

-- + -

70 75 

-- + 

80 

Number of Sites 

FREQUENCY 

INT. CUM. 

1 

2 

1 

0 

5 

3 

2 

2 

4 

6 

8 

4 

3 

6 

2 

5 

4 

1 

2 

2 

8 

5 

2 

8 

8 

5 

10 

7 

6 

11 

14 

13 

10 

5 

0 

0 

0 

0 

1 

3 

4 

4 

9 

12 

i a 

16 

20 

26 

34 

38 

41 

47 

49 

54 

58 

59 

61 

63 

71 

7 6 

70 

OG 

94 

99 

109 

116 

122 

133 

147 

160 

170 

175 

175 

175 

175 

175 

PERCENTAGE 

INT. CUM. 

0.6 

1 . 1 

0.6 

0.0 

2.9 

1 

1 

1 

2 

3 

4 

2 

1 

3 

1 

2 

2 

0 

1 

1 

4 

2 

1 

4 

4 

2 

5 

4 

3.4 

6.3 

8.0 

7.4 

5.7 

2.9 

0.0 

0.0 

0.0 

0.0 

0. 

1 . 

2. 

2. 

5. 

6. 

g .0 

9. 

11 . 

14. 

19. 

21 . 

23. 

26.9 

20. 

30. 

33. 

33. 

34. 

36. 

40. 

'4 3. 

4 4. 

49. 

5 3. 

56. 

62. 

66. 

69. 

76. 

04. 

91 . 

97. 

1 00 . 

100.0 

100.0 

100. a 

100.0 

.0 

9 

1 

7 

o 

.0 

.6 

.4 

,6 

. 1 

,7 

.6 

.3 

.3 

.7 

.0 

.0 

4 

1 

,0 

FIGURE D-lb. Other sites. 

831'if, 
E 

• 7.00000 

•14.0000 

•21.0000 

•28.0000 

•35.0000 

•42.0000 

•49.0000 

•56.0000 

•63.0000 

•70.0000 

•77.0000 

•84.0000 

•91.0000 

•98.0000 

•105.000 

•112.000 

•119.000 

•126.000 

•13 3.000 

•140.0100 

•14 7.000 

•15 4.000 

•161.000 

•16 8.000 

•175.000 

•182.000 

•189.000 

•196.000 

•203.000 

•210.000 

•217.000 

•224.000 

•231.000 

•238.000 

•245.000 

•252.000 

•259.000 

•266.000 

SYMBOL COUNT MEAN ST.DEV. 

A 540 7.445 22.398 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAA 

• AAAA 

•AAAAAA 

• AA 

• A 

•AA 

• AAAA 

+ 

• A 

+ 

• A 

• AAA 

• AAA 

+ 

• AA 

• AAA 

+ 

• AA 

• AAA 

• A 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

FREQUENCY 

INT. 

463 

35 

12 

4 

6 

2 

1 

2 

4 

0 

1 

0 

1 

3 

3 

0 

2 

3 

0 

2 

3 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CUM. 

463 

49 8 

510 

514 

520 

522 

523 

525 

529 

529 

530 

530 

531 

534 

537 

537 

539 

542 

542 

544 

547 

540 

540 

540 

548 

540 

540 

548 

540 

540 

540 

540 

540 

548 

548 

548 

540 

548 

PERCENTAGE 

INT. 

34.5 

6.4 

2.2 

0.7 

1 . 1 

0.4 

0.2 

0.4 

0.7 

0.0 

0.2 

0.0 

0.2 

0.5 

0.5 

0.0 

0.4 

0.5 

0.0 

0.4 

0.5 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

8 4.5 

90.9 

93. 1 

93.8 

94.9 

95.3 

95.4 

95.8 

96.5 

96.5 

96.7 

96.7 

96.9 

97.4 

98.0 

98.0 

98.4 

98.9 

93.9 

99.3 

99.0 

100.0 

100.0 

100.0 

100.0 

100.0 

1 00 .0 

100. 0 

100.0 

100 .0 

100.0 

100 . 0 

100. 0 

100 .0 

100.0 

100.0 

100.0 

100. 0 

Number of Sites 

83146 

FIGURE D-2a. Histogram of the number of expected primary exceedances 

(E) for sites reporting 03 concentrations in 1981: SLAMS (including NAMS) 
SYMBOL COUNT MEAN ST.DEV. 

B 175 9.388 35.657 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

E 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+ + + + + + + + + + + _ + + + + + + 

•7.00000 +B3BBBBB3BBBBBBBBB5BBBBBBB3BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

•14.0000 +BBBBB3BBB 

•21 .0000 +BBBBB 

•28.0000 +BB 

•35.0000 +BBB 

*4 2.00>30 +B 

•49.0000 +BB 

•56.0000 + 

•63.0000 + 

•70.0000 + 

•7 7.000.7 + 

•84.0000 +B 

•91.0000 + 

•90.0000 + 

•105.000 + 

•112.000 + 

•119.000 + 

•126.000 +B 

•133.000 + 

•140.000 + 

•147.000 + 

•15 4.000 + 

•161. .000 + 

•160.0.0.7 + 

•175.00.0 + 

•182.000 + 

•189.000 +B 

•196.000 + 

•203.000 + 

•210.000 + 

•217.000 +B 

•224.0.00 + ' 

•231.000 + 

•238.000 +B 

•245. .00.0 + 

•252.000 + 

•259.000 +B 

•266.000 • 

+ + + + + + + + + + + + + + + + + 

5 1.0 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

Number of Sites 

FIGURE D-2b. Other sites. 

FREQUENCY 

INT. 

147 

9 

5 

2 

3 

1 

2 

. 0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

CUM. 

147 

156 

161 

163 

166 

167 

169 

169 

169 

169 

169 

170 

170 

170 

170 

170 

170 

171 

171 

171 

171 

171 

171 

171 

171 

171 

172 

172 

172 

172 

173 

173 

17.3 

174 

174 

174 

175 

175 

PERCENTAGE 

INT. 

84.0 

5. 1 

2.9 

1.. 1 

1 .7 

0.6 

1 . 1 

0.0 

0.0 

0.0 

0.0 

0.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.6 

0.0 

0.0 

0.0 

0.6 

0.0 

0.0 

0.6 

0..U 

0.0 

0. 0 

0.0 

CUM. 

84.0 

89. 1 

92.0 

93. 1 

94.9 

95.4 

96.6 

96.6 

96.6 

96.6 

96.6 

97. 1 

97. 1 

97. 1 

97. 1 

97. 1 

9 7. 1 

9 7.7 

97.7 

97.7 

97.7 

97.7 

97.7 

97.7 

97.7 

97.7 

90.3 

90.3 

90.3 

98.3 

90.9 

9 0.9 

90.9 

99.4 

9 9.4 

9 9.4 

100.0 

1 00 . 0 
SLAMS 

OTHER 

K- 723 

COR- .0073 

SYMBOL=S 

SYMBOL=0 

, •. 

MEAN 

77.685 

7.9154 

ST.DEV. 

22.238 

26.214 

REGRESSION LINE RES.MS. 

X= .00621*Y+ 77.636 495.19 

Y= .00863*X+ 7.2452 688.08 

. . + + ..X. + + + 

280 + 

245 + 

210 

175 + 

140 + 

105 + 

70. 

35. 

0.0 

0 

OOS 

OO 

S 

0 S 

SO 

0 S 

S 

S 

s 

s s 

s 

SSS 

S 

s 

s 

so 

s 

S S S 

s 

s s 

s s 

s s 

s s s s 

0 s s SOS ss 0 s osossss ss s s 

Y o ssooo s s * s s**os ssssss*s*s*s*sss Y 

+0 0 0 **o*o*o s** **s o****ss*s*ss* *****s*s*********************s + 

.....................+....+....+....+....+....+....+.X..+....+....+....+ . 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

0.000 15.00 30.00 45.00 60.00 75.00 . 90.00 105.0 

FIGURE D-3. Scatter plot of the percentage of data completeness in 1981 

(p) versus the number of estimated primary exceedances (£) for sites 

reporting Og concentrations. 

831U6 
TABLE D-1. Frequency table of the percentage of total days recorded in 1981 (p) 

versus expected primary exceedances (E) for sites reporting 1-hour ozone averages. 

0 41X <= 62.5X <• 75X <- 81.5X <- 86X 89X 92X 95X <- 97.5X <- 100X TOTAL 

0 

1 - 2 

<- 5 

<"255 

TU i AL 

51 

at 

9 

12 

72 

41 

e 

5 

20 

71 

42 

O 

10 

14 

75 

43 

ia 

7 

10 

70 

31. 

1 2 

16 

9 

73 

29 

1 2 

13 

9 

63 

34 

12 

16 

17 

79 

37 

14 

16 

1 1 

78 

25 

12 

8 

20 

65 

30 

13 

15 

19 

77 

368 

99 

115 

141 

723 

83146 

s 
P.50 

Appendix E 

LEAD DESCRIPTIVE STATISTICS 

83146 5 
P.51 

APPENDIX E 

LEAD DESCRIPTIVE STATISTICS 

Figure E-1 shows the distribution of data completeness among air monitoring stations. For SLAMS, the mode of the data completeness distribution 

is in the interval of 15-18 percent, corresponding to nearly complete data 

recovery from a once-per-six-day sampling schedule. There are low 

correlations between the percentage of data completeness and reported 

annual averages, standard deviations, and coefficients of variation 

(Figure E-2). 

83146 5 
'3 

*6. 

l 9. 

•12 

00000 

00000 

00000 

0000 

•15.000.0 

•13.0000 

. 0000 

. 0000 

. 00.00 

. 000.0 

. 0000 

, 000.7 

•39.0000 

•42.0000 

. 0.0.00 

•21 . 

•24. 

•27. 

•30. 

•33. 

•36. 

•45 

*48 

*51 

•54 

•57 

•60 

•63 

•66 

•69 

•72 

•75 

•70 

*31 

*G4 

•87 

•90 

0000 

0000 

0000 

0000 

0000 

0000 

.0000 

0.0.0.0 

0000 

0.000 

0000 

0.000 

0000 

0.000 

.0.00.0 

"93.0.0.00 

•96.000.0 

* 9 9.0.0.00 

* 102.0.00 

•10!}. 000 

•100.00.0 

•111.000 

•114 .£'00 

SYMBOL COUNT 

A 171 

EACH SYMBOL REPRESENTS 

MEAN ST.DEV. 

11.917 10.059 

1 OBSERVATION 

10 1! 55 

-- + -

60 65 70 75 

-- + -

80 

-- + 

20 2 5 30 35 40 45 50 

+ + + ,. + + ^ + + + + . 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

• AA 

•AAA 

•AAAAA 

•AAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 

•AAA 

•A 

•A 

• AA 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

FREQUENCY 

INT. CUM. 

•A 

+ 

10 20 30 35 40 45 60 70 75 80 

49 

2 

3 

5 

19 

05 

3 

0 

0 

1 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

49 

51 

54 

59 

78 

163 

1G6 

166 

166 

167 

168 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

170 

171 

171 

171 

171 

171 

171 

171 

PERCENTAGE 

INT. CUM. 

28. 

1 . 

1 . 

2. 

11 . 

49. 

1 . 

0.0 

0.0 

0.6 

0. 

1 . 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.6 

0.0 

0 . 0 

0.0 

0. 

0. 

0.0 

28.7 

29.8 

3 1 .C7 

34.5 

45.6 

95.3 

97.1 

97. 1 

97. 1 

97 

90 

9 9 

99 

99 

99 

99.4 

99.4 

99.4 

99.4 

99.4 

93.4 

9 9..' 

99.* 

99, 

99. 

99. 

99. 

99. 

99.4 

9 3./; 

99.4 

100.0 

1 00 . 0 

1 00 . 0 

1 00 . .0 

1 00 . 0 

1 00 . 0 

1 00 . 0 

Number of Sites 

FIGURE E-la. 

reporting Pb 

Histogram of the percentage, p, of data completeness for those sites 

in 1981: SLAMS (including NAMS). 

P, 3 1 4 <> 
p 

* 3.00000 

•6.0000.0 

* 9 .0000.1 

•12. 0.000 

•15.0000 

•13.0005 

•21.0000 

•24.0000 

•27.0000 

•30.0000 

•33.0.000 

•36.0000 

•39.0.O00 

•42.0.000 

•45 .0.0.00 

•4 8.0.000 

•51 .0.0.00 

*54 .Cnna 

•5 7.0.000 

•60.00.00 

•6 3.0.000 

•6 6.00.00 

•69.0.080 

•72 , "ff'00 

•7 5.0000. 

•78 .0.0.00 

•81 .0.000 

* S1.0000 

* 8 7.0000 

•90. .000.0 

•9 3..00.00 

•9G..O000 

'•99.0000 

•10 2.000 

•105..000 

• 1 fJ 3 .000 

•111. .00.0 

•114. .000 

5 

SYM30L COUNT MEAN ST.DF.V. 

B 525 7.614 7.023 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+ B BBBB EBBBBBBBBBBB^BBBBCSBBBQBBBBEB BEBBB BBBBBB BBSBBBBBBBBB'L-BBBBBBBBBBBBBBBBBBBBBB * 

• EBECBL'B 3BBBBBQBB 

• BBBBBBBISBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB5BBBBBBBEBBBB* 

•BBBBBBB 

•BBCBBBB 

•BBBBBBB 

+ 

+ 

+ 

+ B 

+ B 

+ B 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5 

iGBBBBBUBP-iCBCBBBBCGUBBBBB 

2B BBBB BBBB ;;• BEE BBBBBBBBBBIJBB 

3BBBEBBBB3'i QBEBBBBBBBBBBBBBBBBBBBBBBBBEBDBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB3B* 

10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

FREQUENCY 

INT. 

212 

16 

95 

33 

34 

132 

0 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

• 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CUM. 

212 

228 

323 

356 

390 

522 

622 

52 2 

5 22 

523 

524 

52 5 

525 

5 2 5 

5 2 5 

525 

625 

525 

5 2 5 

525 

52 5 

515 

5 2 5 

52 5 

5 25 

52 5 

525 

5 25 

52 5 

525 

525 

525 

525 

525 

525 

5 2 5 

5 2 5 

525 

PERCENTAGE 

I N T . 

40.4 

3.0 

18. 1 

6.3 

6.5 

25. 1 

0.0 

0.0 

0.0 

0.2 

0.2 

0.2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

40.4 

43.4 

61 .5 

67.8 

74.3 

99.4 

99.4 

99.4 

99.4 

99.G 

9 9 . U 

100.0 

100.0 

100.0 

1 00 .0 

100.0 

1 0 0 . (? 

100. a 

1 00 . 0 

100.0 

1 00 . 0 

100.0 

100.0 

1 00 . 0 

1 00 . 0 

1 00 . 0 

100.0 

100.0 

100 .0 

1 00 . 0 

1 00 . 0 

1 00 . 0 

1 00 . 0 

1 00 . 0 

100.0 

1 00 . 0 

100.0 

1 00 . 0 

Number of Sites 

FIGURE E-lb. Other sites, 

H ''• ] 'i f, 
SLAMS 

OTHER 

N- 696 

COR=-.0080 

MEAN 

X 8.6711 

Y .44578 

SYMBOL-S 

5YMBOL-0 

ST.DEV. 

8.0850 

.67639 

REGRESSION LINE RES.MS. 

X=-.09538*Y+ 8.7136 65.457 

Y=-668E-6*X+ .45157 .45814 

cn 

JD 

a. 

c 

to 

01 

12. 0 + 

10. 5 + 

9.0 0 + 

7.50 + 

6.00 + 

4.50 

3.00 

1.5 0 

0.0 0 

SSO 

* S O 

S O * 

S* 0 ss * 

*0 0. 00 s** 

** *oo*o*****s s s 

•*00**0****** S OS * ' S Y 

••00*00*0*0** + 

+ ....+X.. . + ... . + ... . + ... . + ... . + ....•... . + ..'..•... . + ... . + ... . + . .. . + ... . + . 

7.500 22.50 37.50 .52.5 0 67.50 82.50 97.50 

000 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE E-2a. Scatter plots of ambient air quality statistics versus the 

percent (p) of data completeness for sites reporting Pb concentrations in 

1981: maximum quarterly average versus p. 

83146 
SLAMS , SYMBOL'S 

OTHER , SYMBOL=0 

N- 670 

COR=-.0068 

MEAN ST.DEV. REGRESSION LINE RES.MS. 

X B.9734 Q.0O.J1 X = - .08 i 29*Y + 8.9991 65.545 

Y .31610 .676o0 Y = -559E-5*:<+ .32120 .45B72 

10.0 0 + 

3.75 0 + 

7.50 0 + 

S 6.25.7 + 

0 

. 0~0 + 

. X . . . . + . 

3 . 75.3 +S 

2.50 0 + S 

0 c 

s 

S oss 

1 .23 0 +* 0 0 0 , 

•0 0 OS O 

<' SO 0 SS:-* 

• ooo o i •••""•••••••••s 

Y^OO-'O —S*-A.. S 

0.J20 + •••••oo 0'---v0--0'A •*•••'• 

7.5.J.3 22.5 0 37.5 0 52.5 0 67.5 0 82.5 0 97.5 0 

0.05 0 15.0 0 30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 105. 0 

FIGURE E-2b. Standard diviation versus p, 
SVMBOL=S 

SYMBOL=0 

SLAMS 

OTHER 

N= 670 

COR=-.0252 

MEAN ST.DEV. REASSIG N LINE P.CS.AS. 

X 8.9734 B..US.71 X = -. S5225vi Y-i- 9.3;;?.9 65.5J7 

Y .54311 .31213 Y = -971 L-5-;!+ ,CJ5182 .09754 

C 

0 

E 

F 

F 

V 

A 

R 

3.2 

2.3 

2. 4 

1.6 

1 .2 

.00 

.40 

0.0 

0 0 

0 

S - s 

0 

0 

3 0 *S 

••* 0 s SO 

* * 0 SS 0 

* 0 s o 

* 00 s S 00 

00 0 * * 0 

•'•' OOOCO OCO-'S 

* o oo*•'•-••••"S 

* SCO - *••••• •••••••vs 

* OOOOOO 0"'•'•'* 

"••0 * 00 0 "*•** * 

•••o OOGGO*** * 

* oooooo o * •••••••• 

4

0 OO-'-'O o ••••••* 

*S 00 0 0 0'* 

0 0 ss. 

- 0 0 

00 0 

0 

0 

s 

S 

s 

0 

7.50 0 22.5 0 . 37.5 0 52.5 0 67.5 0 82.5 0 97.5 0 

0 0 0 15.0 0 30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 105. 0 

FIGURE E-2c. Coefficient of variation (standard deviation divided by 

maximum quarterly average) versus p. 
P.57 

Appendix F 

TSP DESCRIPTIVE STATISTICS 

83146 5 
© 

© 

P.58 

APPENDIX F 

e> TSP DESCRIPTIVE STATISTICS 

CD 

©• 

0 

O 

05 

© 

Figure F-1 shows the distribution of data completeness among air monitoring stations. The mode of the distribution is the interval of 15-18 

percent data completeness, corresponding to nearly complete data recovery 

from a once-per-six-day sampling schedule. Table F-1 shows the number of 

sites meeting an annual quota of 75 percent data completeness and meeting 

a corresponding quarterly quota in a specified number of quarters. Table 

F-2 shows the same information for an annual quota of 50 percent. Note 

that roughly 70 percent of the sites that met an annual quota also met a 

corresponding quarterly quota in each of the four quarters of 1981. 

There are low correlations between the percentage of data completeness and 

reported annual averages, standard deviations, and coefficients of 

variation (Figure F2). 

Histograms of "expected" exceedances (i.e., X- or k./p. in the notation of 

Section 2) show that roughly 95 percent of the stations reported eight or 

fewer expected exceedances in 1981 (Figure F-3). The greatest number of 

expected exceedances were reported by the stations having relatively less 

complete data; 91 percent of the stations reported zero expected exceedances (Figure F-4 and Table F-3). Judging by Table F-3 and the low 

correlation given in Figure F-4, no clear relationship is evident between 

expected exceedances and data completeness. 

© 

83146 5 

® 
SYMBOL COUNT MEAN ST.DEV. 

A 2295 15.722 9.718 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

FREQUENCY PERCENTAGE 

P 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 INT. CUM. INT. CUM. 

+ + + + + + + + + + + + + + + + + 

* 3.00000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 40 

•6.00000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 76 

•9.00000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 86 

•12.0000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 165 

* 15.0000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA• 53 4 

•18.0000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA* 12 40 

•2 1.0000 +AAAAAAAAAAAAAAAAA 17 

•2 4.0000 +AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 4 7 

•2 7.0000 +AAAAAAAAA 9 

* 30.0000 +AAAAAAAA 8 

•3 3.0000 AAAAAAAAAAAAAAAAAAAAAAAAAAAA 2 7 

* 36.0000 + 0 

•39.0000 +AA 2 

•4 2.0000 +AA 2 

* 45.0000 +A 1 

•4 8.0000 +AA 2 

•51.0000 +AA 2 

•54.0000 +A 1 

•57.0000 +A 1 

•60.0000 + 0 

•63.0000 +A 1 

•66.0000 +A 1 

•69.0000 +AAAAAA 6 

•72.0000 +AAA 3 

•75.0000 + 0 

•78.0000 +A 1 

•81.0000 +AAAA 4 

•84.0000 +A 1 

•87.0000 + 0 

•90.0000 +AA 2 

•93.0000 +AAAA 4 

•96.0000 +AAAAAAA 7 

•99.0000 +AAA 3 

•102.000 +AA 2 

•105.000 + 0 

•108.000 + 0 

•111.000 + 0 

•114.000 + 0 

+ + + + + + + + + + + + + + + + + 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

Number of Sites 

FIGURE F-la. Histogram of the percentage, p, of data completeness for those 

sites reporting TSP concentrations in 1981: SLAMS (including NAMS). 

40 

116 

202 

367 

901 

2141 

21G 9 

2205 

2214 

2222 

2249 

2249 

2251 

2253 

2254 

2256 

2258 

2259 

2260 

2260 

2261 

2262 

2263 

2271 

2271 

2272 

2276 

2277 

2277 

2 279 

2283 

2290 

2293 

2295 

2295 

2295 

2295 

2295 

1 

3 

3 

7 

23 

54 

0 

2 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

3 

7 

2 

3 

0 

7 

0 

4 

3 

2 

0 

1 

1 

0 

1 

1 

0 

0 

0 

0 

0 

3 

1 

0 

0 

2 

0 

0 

1 CSI 3 

1 

1 

0 

0 

0 

0 

1 

5 

8 

16 

39 

93 

Q 4 

96 

96 

96 

98 

98 

98 

98 

98 

98 

98 

98 

98 

98 

98 

98 

93 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

100 

100 

100 

100 

100 

7 

1 

8 

0 

3 

3 

a 

1 

5 

8 

0 

0 

1 

2 

2 

3 

4 

4 

5 

5 

5 

6 

8 

0 

0 

0 

2 

2 CSI 3 

5 

8 

9 

0 

0 

0 

0 

0 

8 3 I 'if. 
* 3.00000 

* 6.00000 

* 9.00000 

•12.0000 

•15.0000 

•18.0000 

•21.0000 

•24.0000 

. 0000 

. 0000 

. 0000 

. 0000 

'27. 

'30. 

"33. 

'36. 

"39. 

"42. 

'45, 

'48. 

'51 

'54. 

. 0000 

. 0000 

, 0000 

. 0000 

, 0000 

. 0000 

*57.0000 

•60.0000 

•63.0.00.0 

* 6 6.0000 

•63.0000 

•72.0000 

. 000.0 

.0000 

. 0000 

. 0000 

•3 7.0000 

•90.0000 

•9 3.0000 

•96.0000 

•99.0000 

•102.000 

•105.000 

•108.000 

•111.000 

* 114.000 

•75. 

•78. 

•31 

•84. 

SYMBOL COUNT 

B 1372 

EACH SYMBOL REPRESENTS 

MEAN ST.DEV. 

15.845 15.563 

1 OBSERVATION 

10 

+ 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

• BBBB 

+ BBB 

• BBBB 

+ BBB 

+ B 

+ BBB 

15 

•- + -

20 30 35 40 

-- + -

45 50 55 60 65 70 75 80 

FREQUENCY PERCENTAGE 

INT. CUM. INT. CUM. 

BBB3BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

BBB3BBBBBBBBBBBBBBBCBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

BBBBB3BBBBBBB3BBBBBBB0BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

BBBBBBBBBBBBBBBBBBBBBEBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB* 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB3BBBBBBBBBBBBBBB* 

BBBBBBBBBBBB 

BBBBBBBBBBBBB 

B3BBBBBBBB 

BBBBBBBBBBBBBBBBBBBBBBBBBB 

SBB 

JB 

+ 3BB 

+BBBBE 

+ B 

+BBBBE 

+ B 

• BB 

• BBBB 

+BEBBE 

• BBBB 

+ BU 

• BB 

+BBBBE 

+BBEBE 

+BBBBE 

• BB 

+ B 

• BBBB 

•BBBBBBBB 

• B 

JBBB 

SBBB 

!BBB 

10 

- + -

15 

• - + -

20 25 30 

-- + -

35 

-- + -

40 

-- + -

45 

-- + -

50 

-- + -

55 60 

-- + -

65 

-- + -

70 

-- + -

75 

-- + 

80 

87 

145 

135 

119 

259 

468 

16 

17 

14 

3 

30 

3 

1 

3 

0 

3 

6 

1 

5 

1 

2 

4 

5 

4 

2 

2 

7 

7 

7 

2 

1 

4 

8 

1 

0 

0 

0 

0 

87 

232 

367 

4 36 

745 

1213 

1229 

1246 

1260 

1263 

1293 

1296 

1297 

1300 

1300 

1303 

1309 

1310 

1315 

1316 

1318 

1322 

1327 

1331 

1333 

1335 

1342 

1349 

1356 

1358 

1359 

1363 

1371 

1372 

1372 

1372 

1372 

1372 

6.3 

10.6 

9.8 

8, 

18. 

34. 

1 . 

1 , 

• 1 . 

0. 

2. 

0. 

0. 1 

0.2 

0.0 

0.2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 1 

0. 1 

0.3 

0.6 

0. 1 

0.0 

0.0 

0.0 

0.0 

6. 

16. 

26. 

35. 

54. 

88. 

89. 

90.8 

91 .8 

92. 1 

94.2 

94.5 

94.5 

94.8 

94.8 

95.0 

95.4 

95.5 

95.8 

95.9 

96. 1 

96.4 

96.7 

97.0 

97.2 

97.3 

97.8 

98.3 

98.8 

99.0 

99. 1 

99.3 

99.9 

100.0 

100.0 

100.0 

100.0 

100.0 

Number of Sites 

FIGURE F-lb. Other sites. 
P.61 

TABLE F-1. Number of sites recording 24-hour TSP during at least 

75 percent of the scheduled monitoring days in 1981, by the number of 

quarters in which at least 75 percent of the scheduled days per quarter 

are recorded. 

Site Type 

SLAMS 

(includes 

Other 

No. 

wi 

75% 

NAMS) 

of Quarters 

th at Least 

Days/Quarter 

0 

1 

2 

3 

4 

Total 

0 

1 

2 

3 

4 

Total 

No. of 

> 75% 

0 

2 

69 

452 

1335 

1858 

0 

0 

70 

185 

590 

845 

Sites 

< 75% 

79 

148 

145 

65 

0 

437 

160 

204 

119 

44 

0 

527 

Total 

79 

150 

214 

517 

1335 

2295 

160 

204 

189 

229 

590 

1372 

TABLE F-2. Number of sites recording 24-hour TSP during at least 

50 percent of the scheduled monitoring days in 1981, by the number 

of quarters in which at least 50 percent of the scheduled days 

per quarter are recorded. 

Site Type 

SLAMS 

(includes 

Other 

No. of Quarters 

with at Least 

50% Days/Quarter 

0 

NAMS) 1 

2 

3 

4 

Total 

0 

1 

2 

3 

4 

Total 

No. of 

> 50% 

0 

0 

52 

235 

1825 

2112 

0 

3 

62 

175 

785 

1025 

Sites 

< 50% 

32 

89 

61 

1 

0 

183 

75 

186 

83 

3 

0 

347 

Total 

32 

89 

113 

236 

1825 

2295 

75 

189 

145 

178 

785 

1372 

83146 8 
P.62 

SLAMS 

OTHER 

SYMBOL=S 

SYMBOL=0 

N= 3654 

COR=-.0122 

MEAN 

15.821 

56.526 

ST.DEV. 

12.220 

20.827 

REGRESSION 

X=-.00714*Y+ 

Y=-.02074*X+ 

LINE 

16.22 5 

56.35 4 

RES.MS. 

149.3 4 

433.3 1 

400 + 

350 + 

300 

,X+. 

tn 250 

c 

to 

tv 

i-

+ J 

tu 

B 

o 

Ol 

200 +' 

0 

S 

s 

0 s 

150 + 

s 0 

s s 

0 00*0 

0 0 

100 +s *os**s * 

* * * * 2 '* " 

0 0 * *>'•'*** 

* rt ••": rt rt rt rt * 

) rt * V; rt * rt -.* * 

} * rt rt rt rt rt * rt 

)**;.-Qrtrtrt * 

Q V * * 0 •'-•* 

*0 00 

0 

^... . + . . . 

Y0* 

5 0. +0" 

. 0" 

. 0 " 

0 . 0 + 

**S 

S* 

S* 

rt * rt 

rtrtrt 

rtrt * 

rtrtrt 

rtrtrt 

rtrtrt 

rtrtrt 

rtrtrt 

rtrtrt 

* s 

s 

s 

s s 

s 

ssss 

*s* 

OS 

SOS 

s**s **os s*os*oo 

*s»*« *s***os 

*0***o 0*S 0 

S 0 000 00 

00 00 0** c 

O S 0 

s s s 

S SOS S 0 SOS 0 

0 SS S 0 

0 S *00 S 00 0 00 

00 0 SO 00*000 

00 *0 0 000 00 

00 0 00 0 

0 

S SO 

s 

s oos 

SSSS * S 

so *0 

S 00 

0 0 

0 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

0.000 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE F-2a. Scatter plots of ambient air quality statistics versus the 

percent of data completeness (p) for sites reporting TSP concentrations 

in 1981: annual geometric mean versus p. 

831^6 
P.63 

SLAMS 

OTHER 

SYMBOL=S 

SYMBOL=0 

N- 3654 

COR= .079 7 

MEAN 

15.82 1 

1.56.70 

IT.DEV. REGRESSION LINE RES.MS. 

12.22 0 X= 4.3S43*Y+ 8.951 1 148.4 2 

.2220 0 Y= .00145-X + 1.544 1 .0489 9 

5 . 0 + 

4 . 5 

4 . 0 

3 . ! 

3 . 0 0 

s 

2 . 5 + 0 

2 . 0 + 

1.5 Y 

1 . 0 

*....+ . 

0 

.•X . 

so 

0 

O O 0 

0 

0 s 

COO* 

oos s s 

* 0 * s s 

• 0 0 * OS s 

0 00 * *s s 

SQS00*3"*0S 

0 -•'•' SO*" * '•' * * * * 

0*S-0 ' 

SSO 

S 0 

so 

0 

0 0 

00 

000 

0 

0 

00 rt rt * * A 3 0 * * S * 

* rt Q« * v.- ** *** rt QOS s oooos 

* rt rt rt rt rt rt rt rtrtrt* rtrtrt g Q g * S 

2 rt rt rt rt * rt rt rt * rt rt g 2 * S 

V: rt 

OSO 

rt * rt rt rt rt rt s 

rt rt * rt rt -rt rt i 

* rt rt * * * * 

* OS 

S*SO**OS 0 0 

**OOOS**S SS * 0 

0 

S 0' 

SSO 

s 

0 

. s 

0 

0 

0 0 

S 0 

0 0 

0 0 0 

0 

0 00 0 

*OSOS so 

SS s 

s 

0 

oooo 

OOSOOO 

0 S 00 

S s 

0 

0 

* 0 

OSOS *0SS 

S SS *0S 

SS SS S 

SO 0 

,x+. 

7.50 0 22.5 0 37.5 0 52.5 0 67.5 0 82.5 0 97.5 0 

0 0 0 15.0 0 30.0 0 45.0 0 60.0 0 75.0 0 90.0 0 105. 0 

FIGURE F-2b. Geometric standard deviation versus p. 
2.4 

2.1 

1.8 

C 

0 

E 

F 

F 

V 

A 

R 

1.5 

1.2 

SLAMS 

OTHER 

N= 3654 

COR= .0082 

MEAN 

X 15.821 

Y .03262 

SYMBOL=S 

SYMBOL=0 

ST.DEV. 

12.220 

.03931 

REGRESSION LINE RES.MS. 

X= 2.5264*Y+ 15.739 149.35 

Y= 26GE-7*X+ .03219 .00159 

.90 

.60 

.30 

0.0 

0 S 

0 OS 

0* 0 0 s 

o*s OO***** OS O 000 0 

V/Qrt rt „*r t rtrtrt * X rt * rt * * * Q * Q * * ;VQ * QQ 

<s*s***o S s* .*rt ST****!! ; * ^*Art'Jt**lCts 

0 00 

OOOO *0 *00 00*000 00*000 0 OS *0 

'"SS 0 0 SOSS S OS S SO 0 SSS*S*OS* 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

0.000 15-.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE F-2c. Coefficient of variation (geometric stnadard deviation 

divided by the annual geometric mean) versus p. 

9. 3 1 L, e. 
SYMBOL COUNT MEAN ST.DEV. 

A 2295 1.489 9.438 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

*8.00000 

•16.0000 

•24.0000 

•32.0000 

•40.0000 

•4 8.0000 

•56.0000 

•6 4.0000 

* 7 2 .00.00 

*80.0000 

* 3 3.0000 

* 9 6.0000 

•104.000 

•112. 000 

•120.000 

•128.000 

•136.000 

•144 .000 

•15 2.000 

•160.000 

•163.000 

•17 6.000 

•184.000 

•19 2.000 

* 200.000 

•203.000 

•216.000 

•224.000 

•232.000 

•2 40.000 

* 2'! 0 . 000 

•256.000 

•264.000 

•272.000 

•2 30.000 

•288.000 

•296.000 

•304.000 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

+AAAAAAAAAAAAAAAA 

•AAAAAAA 

•AA 

• AAA 

•AAAA 

•AA 

+ 

• AA 

•AA 

•AA 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

•A 

+ 

+ 

+ 

•A 

+ 

+ 

+ 

+ 

+ 

FREQUENCY 

INT 

196 

57 

16 

7 CSI 3 

4 (SI 0 

2 

2 CSI 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

CUM. 

2196 

225 3 

2269 

2276 

2278 

2281 

2285 

2287 

2287 

22&9 

2291 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2293 

2294 

2294 

2294 

2294 

2295 

2295 

2295 

2295 

2295 

2295 

PERCENTAGE 

INT. 

95.7 

2.5 

0.7 

0.3 

0. 1 

0. 1 

0.2 

0. 1 

0.0 

0. 1 

0. 1 

0. 1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

95.7 

98.2 

9e.9 

99.2 

99.3 

99.4 

99.6 

99.7 

99.7 

99.7 

99.8 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

99.9 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

Number of Sites 

FIGURE F-3a. Histogram of the number of estimated primary exceedances (E) 

for sites reporting TSP concentrations in 1981: SLAMS (including NAMS). 

831M6 
SYMBOL COUNT MEAN ST.DEV. 

B 1372 1.483 8.144 

EACH SYMBOL REPRESENTS 1 OBSERVATION 

E 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

+ + + + + + + + —.-- + - + + + + + + + + 

•8.00000 +BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB*1302 1302 

•16.0000 +BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

•24.0000 +BBBBBB3BBB 

•32.0000 +BBBBBB 

•40.0000 +BBBB 

•4 8.0000 +BBB 

•56.0000 +B 

•6 4.00.90 + 

•7 2.0000 +BB 

•80.0000 +B 

•8 8.0000 +B 

•9 6.0.000 + 

•104.000 + 

•112..CO0 +B 

•120.000 +B 

•128.000 + 

•136.000 +B 

•144.000 + 

•152.000 + 

•160.000 + 

•168.000 + 

M7G.Jf-.O0 + 

•18 4.000 + 

•192..0.00 + 

•2 0.0.000 + 

* 203.000 + 

•216.000 + 

•2 2 4.0.00 + 

•2 32.000 + 

•24 0.000 + 

•24 0.000 + 

•250.000 + 

•264.. 100 + 

•27 2.000 + 

* 280.000 + 

* 23U.000 + 

•20 (3..0.00 + 

•30 4.000 + 

+ + + + + + + + + + + + + + + + + 

5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 

Number of Sites 

FIGURE F-3b. Other sites. 

FREQUENCY 

INT. 

302 

39 

10 

6 

4 

3 

1 

0 (SI 1 

1 

0 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CUM. 

1302 

1341 

1351 

1357 

1361 

1364 

1365 

1365 

1367 

1368 

1369 

1369 

1369 

1370 

137 1 

1371 

1372 

1372 

1372 

13 72 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

1372 

PERCENTAGE 

INT. 

94.9 

2.8 

0.7 

0.4 

0.3 

0.2 

0. 1 

0.0 

0. 1 

0. 1 

0. 1 

0.0 

0.0 

0. 1 • 

0. 1 

0.0 

0. 1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CUM. 

94.9 

97.7 

98.5 

98.9 

99.2 

99.4 

99.5 

99.5 

99.6 

99.7 

99.8 

99.8 

99.8 

99.9 

99.9 

99.9 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100 .0 

100.0 

1.00.0 

100.0 

100 .0 

100.0 

100 .0 

100 .0 

100.0 

100.0 

100.0 

100.0 

1 00 .0 

1 00.0 

100.0 

H ~s I 'if, 
P.67 

SLAMS 

OTHER 

N= 3667 

COR= .0125 

SYMBOL=S 

SYMBOL=0 

MEAN 

15.765 

1.4867 

ST.DEV. 

12.233 

8.9745 

REGRESSION LINE RES.MS. 

X= .01703*Y+ 15.741 149.66 

Y= ,00919*X+ 1.3418 80.551 

280 

245 + 

210 

175 

140 

105 

70. 

35. 

0.0 

,X+. 

0 

0 

S S 

0 OS * 

s ss 

D 0 OO S 0 

0 so so 

0 osss** 

0** S***Q 

ss 

so 00 

s s 

s s 

****** sss* S 000 OS S S 0 0 0 

y******«*ft********«r*****Q**0S0 S0*0*0 *o *o*os*o*oo* *0*000 *S**S*OS* Y 

+....+....X....+....+....+....+....+....+....+....+....+....+....+....+. 

7.500 22.50 37.50 52.50 67.50 82.50 97.50 

000 15.00 30.00 45.00 60.00 75.00 90.00 105.0 

FIGURE F-4. Scatter plot of the percentage of data completeness (p) 

versus the number of estimated exceedances (E) for sites reporting 

TSP concentrations in 1981. 

83146 
TABLE F-3. Frequency table of the percentage of total days recorded in 1981 (p) 

versus expected primary exceedances (E) for sites reporting 24-hour TSP averages, 

E 

0 

i - 10 

<= 20 

<» 260 

TOTAL 

0 - 6% < 

339 

0 

2 

7 

348 

= 10.5X < 

326 

7 

16 

6 

355 

C= 13X < 

308 

28 

3 

8 

347 

- 14.5X <- 15% < 

308 

27 

8 

5 

348 

231 

13 

2 

2 

248 

- 15.5X < 

367 

23 

10 

3 

403 

C- 16X < 

426 

23 

6 

5 

460 

- 16.5X <- 17X < 

542 

21 

12 

8 

583 

219 

17 

10 

3 

249 

- 100X 

257 

52 

9 

8 

326 

TOTAL 

3323 

21 1 

78 

55 

3667 

8314 6 

m 
P.69 

Appendix G 

SUMMARY OF CURRENT DATA COMPLETENESS CRITERIA 

83U6 5 
APPENDIX G 

SUMMARY OF CURRENT DATA COMPLETENESS CRITERIA 

The present data completeness criteria used to establish valid air quality 

summary statistics vary by use and organization. Tables G-l through G-3 

summarize criteria for NADB, the California Air Resources Board, and the 

Indiana Air Pollution Control Division. In addition, the EPA Trends 

Report (EPA, 1983) requires 50 percent or 4380 hourly observations per 

year (or in some cases 4000 hours per year is used). 

The "EPA fact book" in the Trends Report, which summarizes air quality 

levels in Standard Metropolitan Statistical Areas, uses the same criteria 

as in the Trends Report but notes if the volume of annual data is less. 

Summaries from intermittent measurements based on less than 30 days or 

less than the NADB criteria are also noted. 

REFERENCES 

CARB. 1983. California Air Quality Data, Summary of 1982 Air Quality 

Data, Gaseous and Particulate Pollutants. California Air Resources 

Board, Sacramento, CA. 

EPA. 1977. Aeros Manual Series Volume III: Summary and Retrieval. U.S, 

Environmental Protection Agency, Research Triangle Park, NC (EPA450/2-76-009a). 

EPA. 1979a. "Implementation Plans for Lead National Ambient Air Quality 

Standards." 43 Federal Register, 46266. 

83m6 5 
o 

o 

o 

o 

0 

o 

o 

o 

o 

P.71 

EPAo 1979b. Guideline for Lead Monitoring in the Vicinity of Point 

Sources. U.S. Environmental Protection Agency (EPA-450/4-01-006). 

EPA. 1983. National Air Quality and Emission Trends Report, 1981. U.S. 

Environmental Protection Agency, Research Triangle Park, NC. 

EPA. 1984. "Proposed Revisions to the National Ambient Air Quality 

O Standards for Particulate Matter." 49 Federal Register, 10434. 

O 831»t6 5 

uiumH 
P.72 

TABLE G-l. National Aerometric Data Bank's criteri a 

for validit y of data. (Source: EPA, 1977) 

Continuous Sampling (1 , 2, and 4-hour data) 

Quarterly statistics--75% or 1642 hours 

Annual statistics--75% or 6570 hours 

Intermittent Sampling (24-hour data) 

Quarterly statistics—fiv e samples per quarter 

( I f one month has no values, at least two 

values in other months) 

Annual statistics—four valid quarters 

ir 

Note that 12 samples per quarter are required to 

demonstrate attainment of the NAAQS for lead 

(EPA, 1979a, b). This requirement is also under 

consideration for particulates. (EPA, 1984). 

831H6r 8 
o 

o 

o 

o 

o 

o 

. 0 

o 

o 

o 

TABLE G-2. California Air Resources Board's criteri a for representativeness 

of air quality measurements and statistics . (Source: CARB, 1983) 

Representative 

Calendar 

Statistic 

Year | 

( 

Quarter \ 

\ 

/ 

Month < 

( 

Sampling Time 

Period 

4 

representative 

calendar 

quarters 

24 hours 

Or 

< 24 hours 

24 hours 

Or 

< 24 hours 

Day 

Mean of 

M Hour 

Period 

1-hour 

Or 

2-hour; or 

3-hour; or 

24-hour 

24 

8 

6 

3 

2 

1 

Basis of Statistic or 

Requirement 

Number of 

Representative Hours 

or Days Required 

3 representative months 

Based on a daily 

statistic; or 

Based on hourly statistics 

4 or more 24-hour samples 

Based on a daily 

statistic; or 

Based on all hourly 

statistics; or 

Based on all 2-hour 

statistics; or 

Based on all 3-hour 

statistics 

Mo more than 2 consecutive 

samples missing in each 

1/3 day (hours 0 thru 7, 

8 thru 15, 16 thru 23) 

69 or more calendar 

days 

1,638 or more hours 

23 or more calendar 

days 

548 or more hours 

274 or more 2-hour 

statistics 

183 or more 3-hour 

statistics 

6 .or more hours in 

each 1/3 day 

11 or more samples 

7 or more samples 

22 but not more than 26 hours of sampling 

Number of Samples Needed 

22 or more hourly samples 

7 or more hourly samples 

5 or more hourly samples 

4 hourly samples 

3 hourly samples 

2 hourly samples 

30 minutes or more of sampling* 

* ARB Quality Assurance Manual Section 2.0.23. 

O 

I3146 

uHl 
P.74 

TABLE G-3. Indiana Air Pollution Control Division's 

criteria for completeness of data. (Source: Buch, 1983) 

Time Interval Minimum Number of Valid Observations 

A. Continuous 

Hour 

Day 

Month 

Quarter 

Year 

Three-Hour Running Average 

Eight-Hour Running Average 

At least forty-five (45) minutes* 

At least eighteen (18) hours 

At least twenty-one (21) days 

Three (3) months 

At least nine (9) months with 

at least two (2) months per quarter 

Three (3) consecutive hours 

Six.(6) hours out of eight (8) 

consecutive hours 

B. Six-Day Sampling Schedule 

Quarter 

Year 

At least ten (10) samples unless one 

(1) month has zero (0) valid samples, 

then the remaining two (2) months must 

have at least four (4) samples each. 

At least nine (9) months with at least 

two (2) months per quarter and at least 

forty-five (45) samples. 

C. Daily Sampling Schedule 

Quarter 

Year 

At least sixty (60) samples unless one 

(1) month has zero (0) valid samples, 

then each remaining month must have at 

least twenty-four (24) samples. 

At least nine (9) months witn at least 

two (2) months per quarter and at least 

two-hundred seventy-three (273) samples 

*If an hour includes the daily zero/span sequence, the hour is 

considered valid if at least fifty percent (50.0%) of the hour 

is valid. 

83146