The use of wind fields in a land use regression model to predict air pollution concentrations for health exposure studies

A methodology is developed to include wind flow effects in land use regression (LUR) models for predicting nitrogen dioxide (NO₂) concentrations for health exposure studies. NO₂ is widely used in health studies as an indicator of traffic-generated air pollution in urban areas. Incorporation of high-resolution interpolated observed wind direction from a network of 38 weather stations in a LUR model improved NO₂ concentration estimates in densely populated, high traffic and industrial/business areas in Toronto-Hamilton urban airshed (THUA) of Ontario, Canada. These small-area variations in air pollution concentrations that are probably more important for health exposure studies may not be detected by sparse continuous air pollution monitoring network or conventional interpolation methods. Observed wind fields were also compared with wind fields generated by Global Environmental Multiscale-High resolution Model Application Project (GEM-HiMAP) to explore the feasibility of using regional weather forecasting model simulated wind fields in LUR models when observed data are either sparse or not available. While GEM-HiMAP predicted wind fields well at large scales, it was unable to resolve wind flow patterns at smaller scales. These results suggest caution and careful evaluation of regional weather forecasting model simulated wind fields before incorporating into human exposure models for health studies. This study has demonstrated that wind fields may be integrated into the land use regression framework. Such integration has a discernable influence on both the overall model prediction and perhaps more importantly for health effects assessment on the relative spatial distribution of traffic pollution throughout the THUA. Methodology developed in this study may be applied in other large urban areas across the world.     

Empirical Evaluation of Ambient Ozone Interpolation Procedures to Support Exposure Models

ABSTRACT

An intercomparison study has been performed with six empirical ozone interpolation procedures to predict hourly concentrations in ambient air between monitoring stations. The objective of the study is to use monitoring network data to empirically identify an improved procedure to estimate ozone concentrations at subject exposure points. Four of the procedures in the study are currently used in human exposure models (nearest monitors daily mean and maximum, regression estimate used in the U.S. Environmental Protection Agency's (EPA) pNEM, and inverse distance weighting), and two are being evaluated for this purpose (kriging in space and kriging in space and time). The study focused on spatial estimation during June 1-June 5, 1996, with relatively high observed ozone levels over Houston, Texas. The study evaluated these procedures at three types of locations with monitors of varying proximity. Results from the empirical evaluation indicate that kriging in space and time provides excellent estimates of ozone concentrations within a monitoring network, while the more often used techniques failed to capture observed pollutant concentrations. Improved estimation of pollutant concentrations within the region, and thus at subject locations, should result in improved exposure modeling.     

Temporal and spatial statistical analysis of ambient air quality of Assam (India)

ABSTRACT

Present paper represents the spatio-temporal variation of air quality and performances of geostatistical tools for the identification of pollutants zone in various districts of Assam (India). Geographic Information System (GIS) and geostatistical analysis were utilized to estimate the spatio-temporal variations (2015–2017) of gaseous and particulate air pollutants. Data of 23 fixed monitoring stations were collected from the Central Pollution Control Board (CPCB). It was observed that SO₂ and NO_x concentrations are the major pollutants to the deterioration of air quality in Assam State. Exploratory data analysis was considered for the determination of spatial and temporal patterns of air pollutants. Air Quality index (AQI) was calculated based on the air pollutants and particulate matter. Radial Basis Function (RBF) interpolation techniques were used to analyze the spatial and temporal variation of air quality in Assam. Cross-validation is applied to evaluate the accuracy of interpolation methods in terms of Root Mean Square Error (RMSE), Mean Absolute Percentage Error (MAPE), Nash–Sutcliffe Equation (NSE) and Accuracy Factor (ACFT). In 2015, the high value of AQI portrayed in the central and northeast of the state. In 2016, the central and entire east of the study area was recorded the highest value of AQI. In 2017, it was observed that mostly the central part of the state recorded the high value of AQI. The spatio-temporal variation trend of air pollutants provides sound scientific basis for its management and control. This information of air pollution congregation would be valuable for urban planners and decision architects to efficiently administer air quality for health and environmental purposes.     

A holistic approach for optimal design of air quality monitoring network expansion in an urban area

This paper presents an objective methodology for determining the optimum number of ambient air quality stations in a monitoring network. The methodology integrates the multiple-criteria method with the spatial correlation technique. The pollutant concentration and population exposure data are used in this methodology in different ways. In the first stage, the Fuzzy Analytic Hierarchy Process (FAHP) with triangular fuzzy numbers (TFNs) is used to identify the most desirable monitoring locations. The network configuration is then determined on the basis of the concept of sphere of influences (SOIs). The SOIs are dictated by a predetermined cutoff value (r_c) in the spatial correlation coefficients (r) between the pollutant concentrations at the monitoring stations identified from first step and the corresponding concentrations at neighboring locations in the region. Finally, the optimal station locations are ranked by using combined utility scores gained from the first and second steps. The expansion of air quality monitoring network of Riyadh city in Saudi Arabia is used as a case study to demonstrate the proposed methodology.     

Development and field validation of a community-engaged particulate matter air quality monitoring network in Imperial,California, USA

The Imperial County Community Air Monitoring Network was developed as part of a community-engaged research study to provide real-time particulate matter (PM) air quality information at a high spatial resolution in Imperial County, California. The network augmented the few existing regulatory monitors and increased monitoring near susceptible populations. Monitors were both calibrated and field validated, a key component of evaluating the quality of the data produced by the community monitoring network. This paper examines the performance of a customized version of the low-cost Dylos optical particle counter used in the community air monitors compared with both PM_2.5 and PM₁₀ (particulate matter with aerodynamic diameters <2.5 and <10 μm, respectively) federal equivalent method (FEM) beta-attenuation monitors (BAMs) and federal reference method (FRM) gravimetric filters at a collocation site in the study area. A conversion equation was developed that estimates particle mass concentrations from the native Dylos particle counts, taking into account relative humidity. The R² for converted hourly averaged Dylos mass measurements versus a PM_2.5 BAM was 0.79 and that versus a PM₁₀ BAM was 0.78. The performance of the conversion equation was evaluated at six other sites with collocated PM_2.5 environmental beta-attenuation monitors (EBAMs) located throughout Imperial County. The agreement of the Dylos with the EBAMs was moderate to high (R² = 0.35–0.81).

Implications: The performance of low-cost air quality sensors in community networks is currently not well documented. This paper provides a methodology for quantifying the performance of a next-generation Dylos PM sensor used in the Imperial County Community Air Monitoring Network. This air quality network provides data at a much finer spatial and temporal resolution than has previously been possible with government monitoring efforts. Once calibrated and validated, these high-resolution data may provide more information on susceptible populations, assist in the identification of air pollution hotspots, and increase community awareness of air pollution.     

A Preliminary Study of Particulate Emissions from Small Wood Stoves

This paper summarizes the methodology developed to analyze alternative oxidant control strategies of the 1979 Air Quality Plan for the San Francisco Bay Area. The analysis of alternative oxidant control strategies is a complex task, particularly when a grid-based photochemical model is the primary analysis tool. To handle quantitatively spatial and temporal variations in emissions under both existing and projected future conditions, as well as to simulate the effects of a wide variety of control strategies, a system of computer-based models was assembled. The models projected and distributed a number of variables in space and time: population, employment, housing, land use, transportation, emissions, and air quality. Given time and budget constraints, an approach to maximizing the information return from a limited number of model runs was developed. The system was applied in three sequences to determine (1) what future air quality would be if no further controls were implemented, (2) the degree of hydrocarbon and NO_x emission control necessary to attain the oxidant standard, and (3) the effectiveness of alternative stationary source, mobile source, transportation and land use control strategies in contributing to attainment and maintenance of the oxidant standard.

A number of significant modeling assumptions had to be developed in order properly to interpret the modeled results in the context of the oxidant standard. In particular, a Larsen-type analysis was used to relate modeled atmospheric conditions to “worst case” conditions, and a proportional assumption was made to compensate model results for an imperfect validation. The specification of initial and boundary conditions for future year simulations was found to be a problem in need of further research.     

Air quality mapping using GIS and economic evaluation of health impact for Mumbai City,India

Mumbai, a highly populated city in India, has been selected for air quality mapping and assessment of health impact using monitored air quality data. Air quality monitoring networks in Mumbai are operated by National Environment Engineering Research Institute (NEERI), Maharashtra Pollution Control Board (MPCB), and Brihanmumbai Municipal Corporation (BMC). A monitoring station represents air quality at a particular location, while we need spatial variation for air quality management. Here, air quality monitored data of NEERI and BMC were spatially interpolated using various inbuilt interpolation techniques of ArcGIS. Inverse distance weighting (IDW), Kriging (spherical and Gaussian), and spline techniques have been applied for spatial interpolation for this study. The interpolated results of air pollutants sulfur dioxide (SO₂), nitrogen dioxide (NO₂) and suspended particulate matter (SPM) were compared with air quality data of MPCB in the same region. Comparison of results showed good agreement for predicted values using IDW and Kriging with observed data. Subsequently, health impact assessment of a ward was carried out based on total population of the ward and air quality monitored data within the ward. Finally, health cost within a ward was estimated on the basis of exposed population. This study helps to estimate the valuation of health damage due to air pollution.

Implications: Operating more air quality monitoring stations for measurement of air quality is highly resource intensive in terms of time and cost. The appropriate spatial interpolation techniques can be used to estimate concentration where air quality monitoring stations are not available. Further, health impact assessment for the population of the city and estimation of economic cost of health damage due to ambient air quality can help to make rational control strategies for environmental management. The total health cost for Mumbai city for the year 2012, with a population of 12.4 million, was estimated as USD8000 million.     

Characterizing the spatiotemporal variability of PM2.5 in Cusco,Peru using kriging with external drift

Advancing the understanding of spatiotemporal aspects of air pollution in the urban environment is an area where improved methods can be of great benefit to exposure assessment and policy support. This paper explores the potential of a technique known as kriging with external drift (KED) to provide high resolution maps of fine particulate matter for a downtown region of Cusco, Peru. There were three stages in this research. The first was to conduct a pilot level monitoring campaign to investigate ambient, regional, and street-level air pollutant concentrations for particulate matter (PM_2.5, PM₁₀) and carbon monoxide (CO) in the Province of Cusco. The second was to compile observations within a geographic information system (GIS) in order to characterize the proximal effect of the local transportation network, elevation, and land use classifications on PM_2.5. Third, regression, ordinary kriging and kriging with external drift were used to model PM_2.5 for three select time periods during a 24-h day. Statistical evaluations indicate kriging with external drift resulted in the strongest models explaining 64% of variability seen with morning particle concentrations, 25% for afternoon particles, and 53% in evening particles. These models capture spatial and temporal variability for air pollution in Cusco. These variations seem to be influenced, to varying degrees, by elevation, meteorological conditions, spatial location, and transportation characteristics. In conclusion, combining GIS, meteorological data and geostatistics proved to be a complementary suite of tools for incorporating spatiotemporal analysis into the air quality assessment.     

A land use regression for predicting fine particulate matter concentrations in the New York City region

We developed regression equations to predict fine particulate matter (PM_2.5) at air monitoring locations in the New York City region using data on nearby traffic and land use patterns. Three-year averages (1999–2001) of PM_2.5 at US Environmental Protection Agency (EPA) monitors in the 28 counties including and surrounding New York City were calculated using daily data from the EPA's Air Quality Subsystem. As the secondary contribution to PM_2.5 concentrations is lowest in the winter, we also calculated and modeled average winter 2000 PM_2.5 to conduct a preliminary evaluation of model sensitivity to source contribution. Candidate predictor variables included traffic, land use, census and emissions data from local, state and national sources and were tabulated for a series of circular buffer regions at varying distances around the monitors using a geographic information system. In total, more than 25 variables at 5 different buffer distances were considered for inclusion in the model. Before evaluating the variables we removed several samples from the modeling for validation. For comparison and validation purposes we computed both a model using data for the full 28-county region as well as a more urbanized 9-county region. We found that traffic within a buffer of 300 or 500 m explains the greatest proportion of variance (37–44%) in all 3 models. Measures of urbanization, specifically population density, explain a significant amount of the residual variation (7–18%) after including a traffic variable. Finally, a measure of industrial land use further improves the 28-county and 9-county models based on the 3-yr annual averages, explaining an additional 4% and 11% of the variation, respectively, while vegetative land use improves the winter model explaining an additional 6%. The final models predicted well at validation locations. In total, the final land use regression models explain between 61% and 64% of the variation in PM_2.5.     

Improved space–time forecasting of next day ozone concentrations in the eastern US

There is an urgent need to provide accurate air quality information and forecasts to the general public and environmental health decision-makers. This paper develops a hierarchical space–time model for daily 8-h maximum ozone concentration (O₃) data covering much of the eastern United States. The model combines observed data and forecast output from a computer simulation model known as the Eta Community Multi-scale Air Quality (CMAQ) forecast model in a very flexible, yet computationally fast way, so that the next day forecasts can be computed in real-time operational mode. The model adjusts for spatio-temporal biases in the Eta CMAQ forecasts and avoids a change of support problem often encountered in data fusion settings where real data have been observed at point level monitoring sites, but the forecasts from the computer model are provided at grid cell levels. The model is validated with a large amount of set-aside data and is shown to provide much improved forecasts of daily O₃ concentrations in the eastern United States.     

Sources of variability in ambient air toxics monitoring data

Under the Clean Air Act Amendments, the United States Environmental Protection Agency is required to regulate emissions of 188 hazardous air pollutants. The EPA, Office of Air Quality Planning and Standards is currently conducting a National-scale Air Toxics Assessment with a goal to identify air toxics which are of greatest concern, in terms of contribution to population inhalation risk. The results will be used to set priorities for the collection of additional air toxics emissions and monitoring data. Expanded ambient air toxics monitoring will take the form of a national air toxics monitoring network. With all monitoring data, however, comes uncertainty in the form of environmental variability (spatial and temporal) and monitoring error (sample collection and laboratory analysis). With this in mind, existing data from the Urban Air Toxics Monitoring Program (UATMP) were analyzed to obtain a general understanding of these sources of variability and then provide recommendations for managing the data uncertainties of a national network. The results indicate that environmental variability, in particular temporal, comprises most of the overall variability observed in the UATMP data. However, at lower ambient levels (on the order of 0.1–0.5 ppbv or lower) environmental variability tends to dissipate and monitoring error takes over, most notably analytical error. Overall, the results suggest that common techniques in ambient air toxics monitoring for carbonyls and volatile organic compounds may satisfy many of the primary objectives of a national air toxics monitoring network.     

A Comparative Validation of the RAM and PTMTP Models for Short-Term SO2 Concentrations in Two Urban Areas

A comprehensive and comparative model validation of two EPA models for short-term SO₂ concentrations was performed. The two models tested were RAM (Urban version) and PTMTP (Terrain version). Both are multiple source, multiple receptor gaussian plume models, recommended in the EPA Guideline On Air Quality Models. ¹ The principal difference between the two models is in their use of empirical dispersion coefficients. It was because of the potential for markedly different predicted maximum SO₂ concentrations, and the absence of any testing data on the RAM model, that the validation analysis was undertaken. The current study utilized a full year of air quality data from monitoring sites in two Indiana cities, Michigan City and Indianapolis. Cumulative frequency distributions for each site and model were prepared and comparisons made. The results indicate that the RAM (Urban) model was highly inaccurate in predicting maximum short-term SO₂ concentrations. The PTMTP model, although conservative in its estimates, produces results which more closely resemble the distribution of observed SO₂ concentrations. The body of information presented in this paper is directed to environmental scientists responsible for air quality modeling, and to those persons who set policy on the use of models in air quality studies.     

Regulation of Indirect Sources

The Division of Air Pollution Control, Illinois Environmental Protection Agency, has conducted an ambient air quality monitoring project focusing on carbon monoxide levels in and around several indirect sources. An analysis of the data indicates that highway-type pollutant emissions have the greatest impact on receptors in the vicinity of indirect sources. This implies that the principal, localized constraint on the siting of indirect sources will be the carbon monoxide generated on public roadways servicing those indirect sources. Clearly, adequate procedures must be developed to link such highway-type emissions to pollutant concentrations. An area-source model and a line-source model were tested using the data generated during the monitoring project. Favorable results were achieved using the line-source model. The proper siting of indirect sources involves the allocation of roadway capacity by the governmental units responsible for transportation network design, working in conjunction with regional planning bodies. A regulatory structure is suggested which emphasizes a regional approach, and an example of an air quality allocation scheme is given. The methodology is applicable to all automotive air pollutants although, in general, localized sensitivity is lost for N0₂ and photochemical oxi-dants.     

Comparison of land-use regression models between Great Britain and the Netherlands

Land-use regression models have increasingly been applied for air pollution mapping at typically the city level. Though models generally predict spatial variability well, the structure of models differs widely between studies. The observed differences in the models may be due to artefacts of data and methodology or underlying differences in source or dispersion characteristics. If the former, more standardised methods using common data sets could be beneficial. We compared land-use regression models for NO₂ and PM_10, developed with a consistent protocol in Great Britain (GB) and the Netherlands (NL).Models were constructed on the basis of 2001 annual mean concentrations from the national air quality networks. Predictor variables used for modelling related to traffic, population, land use and topography. Four sets of models were developed for each country. First, predictor variables derived from data sets common to both countries were used in a pooled analysis, including an indicator for country and interaction terms between country and the identified predictor variables. Second, the common data sets were used to develop individual baseline models for each country. Third, the country-specific baseline models were applied after calibration in the other country to explore transferability. The fourth model was developed using the best possible predictor variables for each country.A common model for GB and NL explained NO₂ concentrations well (adjusted R² 0.64), with no significant differences in intercept and slopes between the two countries. The country-specific model developed on common variables for NL but not GB improved the prediction.The performance of models based upon common data was only slightly worse than models optimised with local data. Models transferred to the other country performed substantially worse than the country-specific models. In conclusion, care is needed both in transferring models across different study areas, and in developing large inter-regional LUR models.     

Predictions of U.K. regulated power station contributions to regional air pollution and deposition: a model comparison exercise

Contributions of the emissions from a U.K. regulated fossil-fuel power station to regional air pollution and deposition are estimated using four air quality modeling systems for the year 2003. The modeling systems vary in complexity and emphasis in the way they treat atmospheric and chemical processes, and include the Community Multiscale Air Quality (CMAQ) modeling system in its versions 4.6 and 4.7, a nested modeling system that combines long- and short-range impacts (referred to as TRACK-ADMS [Trajectory Model with Atmospheric Chemical Kinetics-Atmospheric Dispersion Modelling System]), and the Fine Resolution Atmospheric Multi-pollutant Exchange (FRAME) model. An evaluation of the baseline calculations against U.K. monitoring network data is performed. The CMAQ modeling system version 4.6 data set is selected as the reference data set for the model footprint comparison. The annual mean air concentration and total deposition footprints are summarized for each modeling system. The footprints of the power station emissions can account for a significant fraction of the local impacts for some species (e.g., more than 50% for SO2 air concentration and non-sea-salt sulfur deposition close to the source) for 2003. The spatial correlation and the coefficient of variation of the root mean square error (CVRMSE) are calculated between each model footprint and that calculated by the CMAQ modeling system version 4.6. The correlation coefficient quantifies model agreement in terms of spatial patterns, and the CVRMSE measures the magnitude of the difference between model footprints. Possible reasons for the differences between model results are discussed. Finally, implications and recommendations for the regulatory assessment of the impact of major industrial sources using regional air quality modeling systems are discussed in the light of results from this case study.     

Determination of air quality zones in Turkey

In this study, the particulate matter (with an aerodynamic diameter <10 μm; PM₁₀) profile of Turkey with data from the air quality monitoring stations located throughout the country was used. The number of stations (119) was reduced to 55 after a missing data treatment for statistical analyses. First, a classification method was developed based on ongoing national and international (European Commission directives) legislations to categorize air zones into six groups, from a “Very Clear Air Zone” to a “Polluted Air Zone.” Then, a Geographic Information System (GIS)-based interpolation technique and statistical analyses (correlation analysis and factor analysis) were used to generate PM₁₀ pollution profiles of the annual heating time and nonheating time periods. Finally, the coherent air pollution management zones of Turkey, based on air quality criteria and measured data using a GIS-based model supported by statistical analyses, were suggested. Based on the analysis, four hot spots were identified: (i) the eastern part of the Black Sea region; (ii) the northeastern part of inland Anatolia; (iii) the western part of Northeastern Anatolia; and (vi) the eastern part of Turkey. The possible reasons for the elevated PM₁₀ levels are discussed using topographic, climatologic, land use, and energy utilization parameters. Finally, the suggested air zones were compared with the administrative air zones, which were newly developed by the Turkish Ministry of Environment and Forestry, to evaluate the level of agreement between the two.

Implications: The evaluation of air quality profiles of specific regions is important in the development and/or application of an effective air quality management strategy. Factor analysis (FA), together with correlation analysis (CA), provides useful information to classify air pollution management areas over regional networks that have historical time-series air quality data. In this study, which relied on a FA- and CA-based methodology, the coherent air pollution management zones of Turkey after using a GIS-based model were suggested. Policy makers and scientist can use these suggested zones to construct better air quality management strategies.     

Effects of Land Use Data on Dry Deposition in a Regional Photochemical Model for Eastern Texas

ABSTRACT

Land use data are among the inputs used to determine dry deposition velocities for photochemical grid models such as the Comprehensive Air Quality Model with extensions (CAMx) that is currently used for attainment demonstrations and air quality planning by the state of Texas. The sensitivity of dry deposition and O₃ mixing ratios to land use classification was investigated by comparing predictions based on default U.S. Geological Survey (USGS) land use data to predictions based on recently compiled land use data that were collected to improve biogenic emissions estimates. Dry deposition of O₃ decreased throughout much of eastern Texas, especially in urban areas, with the new land use data. Predicted 1-hr averaged O₃ mixing ratios with the new land use data were as much as 11 ppbv greater and 6 ppbv less than predictions based on USGS land use data during the late afternoon. In addition, the area with peak O₃ mixing ratios in excess of 100 ppbv increased significantly in urban areas when deposition velocities were calculated based on the new land use data. Finally, more detailed data on land use within urban areas resulted in peak changes in O₃ mixing ratios of ~2 ppbv. These results indicate the importance of establishing accurate, internally consistent land use data for photochemical modeling in urban areas in Texas. They also indicate the need for field validation of deposition rates in areas experiencing changing land use patterns, such as during urban reforestation programs or residential and commercial development.     

National estimates of outdoor air toxics concentrations

The Clean Air Act identifies 189 hazardous air pollutants (HAPs), or "air toxics," associated with a wide range of adverse human health effects. The U.S. Environmental Protection Agency has conducted a modeling study with the Assessment System for Population Exposure Nationwide (ASPEN) to gain a greater understanding of the spatial distribution of concentrations of these HAPs resulting from contributions of multiple emission sources. The study estimates year 1990 long-term outdoor concentrations of 148 air toxics for each census tract in the continental United States, utilizing a Gaussian air dispersion modeling approach. Ratios of median national modeled concentrations to estimated emissions indicate that emission totals without consideration of emission source type can be a misleading indicator of air quality. The results also indicate priorities for improvements in modeling methodology and emissions identification. Model performance evaluation suggests a tendency for underprediction of observed concentrations, which is likely due, at least in part, to a number of limitations of the Gaussian modeling formulation. Emissions estimates for HAPs have a high degree of uncertainty and contribute to discrepancies between modeled and monitored concentration estimates. The model's ranking of concentrations among monitoring sites is reasonably good for most of the gaseous HAPs evaluated, with ranking accuracy ranging from 66 to 100%.     

Spatial surrogates for the disaggregation of CORINAIR emission inventories

CORINAIR atmospheric emission inventories are frequently used input data for air quality models with a domain situated in Europe. In CORINAIR emission inventories, sources are broken down over 11 major source categories. This paper presents spatial surrogates for the disaggregation of CORINAIR atmospheric emission inventories for input of air pollutants and particulate matter to grid or polygon based air quality model domains inside Europe. The basis for the disaggregation model was the CLC2000 land cover data to which statistical weights were added. Weights were population census data for residential emissions, employment statistics for agricultural and industrial area emissions, livestock statistics for ammonia emissions and annual aircraft movements for emissions realized by air transport. Additional road and off-road network information was used to disaggregate emissions realized by traffic. A comparison of top down produced emission estimates with spatially resolved national emission data for The Netherlands and the United Kingdom gave confidence in the present spatial surrogates as a tool for the top down production of atmospheric emission maps. Explained variance at a spatial resolution of 5 km was >70% for CO, NMVOC and NO_x, >60% for PM10 and almost 50% for SO₂.     

An Evaluation of CO2, Measurements as an Indicator of Air Pollution

Measurements of ambient carbon dioxide (CO₂), made at the Continuous Air Monitoring Program station in downtown Cincinnati, Ohio, and at a rural location near Cincinnati are presented and evaluated to determine the significance of CO₂ data in urban air quality monitoring programs. Through analysis of rural CO₂ data and evaluation of combustion-sources by means of a diffusion model, it is demonstrated that the variation of urban CO₂ concentrations around the prevailing atmosphere background level results from combustion-and noncombustion {natural) sources. The concentration from natural sources can be substantial and in fact override the combustion sources. Because it is not yet practical to predict the contribution of natural sources to urban CO₂ concentrations, data obtained for-this gas have only limited utility as an index of air quality. Significant statistical relationships between CO₂ data and air quality measurements for summer months are shown to result from similar meteorological effects rather than similar sources. A seasonal and spatial variation of ih-ese relationships is postulated and subsequently demonstrated by analysis of CO₂ and air quality measurements from New Orleans, Louisiana, and SU Louis, Missouri.