InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Environmental Engineering » "Air Quality", book edited by Ashok Kumar, ISBN 978-953-307-131-2, Published: August 18, 2010 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 2

Development and Application of a Methodology for Designing a Multi-Objective and Multi-Pollutant Air Quality Monitoring Network for Urban Areas

By Nicolas Mazzeo and Laura Venegas
DOI: 10.5772/9752

Article top

Overview

Location of the city of Buenos Aires and an aerial view of Buenos Aires city.
Figure 1. Location of the city of Buenos Aires and an aerial view of Buenos Aires city.
Estimated emissions of NOx (expressed as NO2) by source category. (*) cars, taxis, trucks.
Figure 2. Estimated emissions of NOx (expressed as NO2) by source category. (*) cars, taxis, trucks.
Estimated emissions of CO by source category. (*) cars, taxis, trucks.
Figure 3. Estimated emissions of CO by source category. (*) cars, taxis, trucks.
Estimated emission of PM10 by source category. (*) cars, taxis, trucks.
Figure 4. Estimated emission of PM10 by source category. (*) cars, taxis, trucks.
Mean annual background concentrations of NO2 in Buenos Aires city.
Figure 5. Mean annual background concentrations of NO2 in Buenos Aires city.
Mean annual background concentrations of CO in Buenos Aires city.
Figure 6. Mean annual background concentrations of CO in Buenos Aires city.
Mean annual background concentrations of PM10 in Buenos Aires city.
Figure 7. Mean annual background concentrations of PM10 in Buenos Aires city.
Horizontal distribution of the total exceedance score (
								
									
										
											E
											S
											=
											
												∑
												
													k
													=
													1
												
												4
											
											E
											
												S
												k
											
										
									
								
								
							) in the city of Buenos Aires.
Figure 8. Horizontal distribution of the total exceedance score ( E S = ∑ k = 1 4 E S k ) in the city of Buenos Aires.
Horizontal distribution of the population factor (PF) in the city of Buenos Aires.
Figure 9. Horizontal distribution of the population factor (PF) in the city of Buenos Aires.
Horizontal distribution of the Total Score (TS) in the urban area.
Figure 10. Horizontal distribution of the Total Score (TS) in the urban area.
Monitoring locations selected considering the total score and the distance constraint (preliminary configuration).
Figure 11. Monitoring locations selected considering the total score and the distance constraint (preliminary configuration).
Sites chosen to locate an air quality monitoring station and pollutants to be measured.
Figure 12. Sites chosen to locate an air quality monitoring station and pollutants to be measured.
“Spatial representativeness” for NO2 monitoring stations (▲).
Figure 13. “Spatial representativeness” for NO2 monitoring stations (▲).
“Spatial representativeness” for CO monitoring stations (●).
Figure 14. “Spatial representativeness” for CO monitoring stations (●).
“Spatial representativeness” for PM10 monitoring stations (▼).
Figure 15. “Spatial representativeness” for PM10 monitoring stations (▼).

Development and application of a methodology for designing a multi-objective and multi-pollutant air quality monitoring network for urban areas

Nicolás A. Mazzeo1 and Laura E. Venegas

1. Introduction

Air pollution has been with us since the first fire was lit, although different aspects have been important at different times. Air pollutants are substances which, when present in the atmosphere under certain conditions, may become injurious to human, animal, plant or microbial life, or to property, or which may interfere with the use and enjoyment of life or property. Air pollution is, however enacted on all geographical and temporal scales, ranging from strictly “here and now” problems related to human health and material damage, over regional phenomena like acidification and forest die back with a time horizon of decades, to global phenomena, which over the next centuries can change the conditions for man and nature over the entire globe.

Three classes of factors determine the amount of pollution at a site: a) the nature of relevant emissions, b) the state of the atmosphere and c) topographical aspects.

In this respect the cities act as sources. Cities are by nature concentrations of humans, materials and activities. They therefore exhibit both the highest levels of pollution and the largest targets of impact. Air pollution problems in urban areas generally are of two types. One is the release of primary pollutants and the other is the formation of secondary pollutants. Since a major source of pollutants is motor vehicles, “hot spots” of high concentrations can occur especially near multilane intersections where the emissions are especially high from idling vehicles. The “hot spots” are exacerbated if high buildings surround the intersection, since the volume of air in which the pollution is contained is severely restricted. The combination of these factors results in high concentrations. These cause effects on health and the environment. Increasingly rigorous legislation, combined with powerful societal pressures, is increasing our need for impartial and authoritative information on the quality of the air we all breathe.

Monitoring is a powerful tool for identifying and tackling air quality problems, but its utility is increased when used, in conjunction with predictive modelling and emission assessment, as part of an integrated approach to air quality management (Rao, 2009).

The monitoring of air pollution level is of significance especially to those residents living in urban areas. Planning and location air quality monitoring networks is an important task for environmental protection authorities, involving: a) ensuring that air quality standards are achieved, b) planning and implementing air quality protection and air pollution control strategies, and c) preventing or responding quickly to air quality deterioration. Therefore, environmental protection authorities need to plan and install air quality monitoring networks effectively and systematically. There are no hard and fast rules for air quality network design, since any decisions made will be determined ultimately by the overall monitoring objectives and resource availability.

Before starting the air quality monitoring network design it is essential to establish what problem has to be solved and what constraints have to be imposed on an “ideal” measuring system. The overall objectives of the monitoring network have to be clearly stated. Some of the specific monitoring objectives can be: to quantify ambient air quality and its variation in space and time; to provide data for air pollution control regulations; to provide real-time data for an alert and warning system; to provide trends for identifying future problems or progress against management/control targets; to provide data for development/validation of management tools.

The goals of this study are: a) to develop an objective procedure to determine the monitoring site locations to detect urban background air pollutant concentrations greater than reference concentrations in an urban area, taking into account the consideration of “protection capability” for areas with higher population density, b) to apply the proposed methodology for designing a multi-pollutant (NO2, CO and PM10,) urban air quality network for Buenos Aires city and c) to evaluate “the spatial representativeness” of mean concentrations measured at each monitoring station. The proposed network design methodology is based on the analysis of the results of atmospheric dispersion models; an exceedance score; a population factor and on the application of the t-Student test for comparison air pollutant mean concentrations at different sites.

2. Introduction to Air Quality Monitoring Network Design

Since one cannot expect to monitor air quality at all locations at all times, selection of sites to give a reliable and realistic picture of air quality becomes a problem in the efficient use of limited resources. The selection of monitoring objectives for optimal allocation of air quality monitoring stations may have to cover several design principles. The required design principles usually consist of the considerations of protection capability for regions with higher population density and significant area with higher economic growth as well as the detection capability of higher pollution concentrations, higher frequency of violation of stipulated standards, and the major industrial/traffic sources in an urban region. Moreover, the cost for siting a pollutant-specific monitoring network would be higher than that for a common monitoring network with respect to several pollutants simultaneously. Thus, for practical reasons, most monitoring networks install different detection instruments together in a common monitoring network that could be viewed as more economic and feasible applications.

Even with a clear set of network objectives, there is no universally accepted methodology for implementing such objectives into the network design, with the approaches used being as varied as the regions being managed. Different methodologies on air quality monitoring network design have been reported in the literature. Among them, statistical methods take advantage of the fact that most air quality measurements are correlated either in time at the same location or in space with other monitors in a network. In this way, networks can be optimized by examining time series correlations from long measurement records or spatial correlations among measurements from many nearby monitors (Munn, 1975, 1981; Elsom, 1978). Various statistical and optimization schemes were applied for designing a representative air quality monitoring network with respect to a pollutant-specific case (Smith & Egan, 1979; Graves et al., 1981; Pickett & Whiting, 1981; Egmond & Onderdelinden, 1981; Handscombre & Elsom, 1982; Husain & Khan, 1983; Nakamori & Sawaragi, 1984; Modak & Lohani, 1985a,b; Liu et al, 1986; Langstaff et al., 1987, Hwang & Chan, 1997). Furthermore, Noll & Mitsutome (1983) developed a method to establish monitor locations based on expected ambient pollutant dosage. This method ranked potential locations by calculating the ratio of a station’s expected dosage over the study area’s total dosage.

It usually happens that an initial monitoring network evolves over time. Therefore after some time a redesign may be required to maximize its capacity to meet modern demands. In this case, it may be desirable the new network maximizes the amount of information it will provide about the environmental field it is being asked to monitor. Equivalently, it should maximally reduce uncertainty about that field. These ideas can be formalized through the use of entropy that quantifies uncertainty and can be used as an objective function. Caselton et al. (1992) used it to rank monitoring sites for possible elimination, an idea extended by Wu & Zidek (1992). Recently, Ainslie et al. (2009) used the entropy-based approach of Le & Zidek (2006) to redesign a monitoring network in Vancouver (Canada) using hourly ozone concentration.

The consideration of multi-pollutant air quality monitoring network design with respect to different objectives was introduced in a series of papers by Modak & Lohani (1985a,b,c). The design principles of a minimum spanning tree algorithm for single or multiple pollutants with respect to one or two objectives was illustrated in these studies. Kainuma et al. (1990) developed a similar procedure to evaluate several types of siting objectives and used a multi-attribute utility function method to determine optimal locations.

Several methods of air quality monitoring design or optimization also include the analysis of atmospheric dispersion models estimations (Hougland & Stephens, 1976; Koda & Seinfeld, 1978; McElroy et al., 1986; Mazzeo & Venegas, 2000 , 2008; Tseng and Chang, 2001; Baldauf et al., 2002; Venegas & Mazzeo; 2003a, 2010). For example, Hougland & Stephens (1976) selected monitoring site locations maximizing coverage factors, such as strength of emission source, distance from the source, and local meteorology for each source included in the study. The basis of this "source oriented" method was to consider for each source and wind direction, the monitor with the largest coverage factor. Koda & Seinfeld (1978) presented a methodology for distributing a number of monitoring stations in a study area in order to obtain the maximum sensitivity of the collected data to achieve the variations in the emissions of the sources of interest. The developed methodology used model estimations of ground level concentrations of pollutants for different meteorological scenarios. McElroy et al. (1986) applied air quality simulation models and population exposure information to produce representative combined patterns and then employed the concept of ‘sphere of influence’ (SOI) developed by Liu et al. (1986) to determine the minimum number of sites required. The monitor’s SOI is defined as the area over which the air quality data for a given station can be considered representative, or can be extrapolated, with known confidence. The site’s SOI can be determined using the covariance structure of the concentrations. Thus, a monitor site’s SOI comprises those neighbouring sites whose variance can be explained by the original site’s variance within a certain degree of confidence.

Tseng & Chang (2001) integrated a series of simulation and optimization techniques for generating better siting alternatives of air quality monitoring stations in an urban environment. The analysis presented used atmospheric dispersion models to estimate air pollution concentrations required in the optimization analysis. Three planning objectives for the minimization of the impacts of the highest concentrations and the highest frequency of violation, as well as the maximization of the highest protection potential of population were emphasized subject to budget, coverage effectiveness (the ratio between effective detection area and total detection area for a monitoring station), spatial correlation, or concentration differentiation constraints. In this case, the concentration differentiation constraints takes into account that the spatial correlation between grids can be high, but the order of magnitude of measured or predicted concentrations between both grids may present significant difference, given the fact that grids are only spatially correlated in terms of concentration pattern.

Baldauf et al. (2002) presented a simple methodology for the selection of a neighbourhood-scale site for meeting either of the following two objectives: to locate monitors at the point of maximum concentration or at a location where a population oriented concentration can be measured. The proposed methodology is based on analyzing middle-scale (from 100 to 500 m) atmospheric dispersion models estimations within the area of interest.

Sarigiannis & Saisana (2008) presented a method for multi-objective optimization of air quality monitoring systems, using both ground-based and satellite remote sensing of the troposphere. This technique used atmospheric turbidity as surrogate for air pollution loading. In their study, Sarigiannis & Saisana (2008) also defined an information function approach combining the values of the violation score, the land-use score, the population density, the density of cultural heritage sites and the cost function. Furthermore, similarities among locations were assessed via the linear correlation coefficient between locations. A gain of information was defined as the product between the correlation coefficient and the information function. The location with the maximum value of the gain information was selected as the best monitoring location.

Elkamel et al. (2008) presented an interactive optimization methodology for allocating the number of sites and the configuration of an air quality monitoring network in a vast area to identify the impact of multiple pollutants. They introduced a mathematical model based on the multiple cell approach to create monthly spatial distributions for the concentrations of the pollutants emitted from different emission sources. These spatial temporal patterns were subject to a heuristic optimization algorithm to identify the optimal configuration of a monitoring network. The objective of the optimization was to provide maximum information about multi-pollutants emitted from each source within a given area.

Pires et al. (2009) applied principal component analysis to identify redundant measurements in air quality monitoring networks. To validate their results, authors used statistical models to estimate air pollutant concentrations at removed monitoring sites using the concentrations measured at the remaining monitoring sites.

Mofarrah & Husain (2010) presented an objective methodology for determining the optimum number of ambient air quality stations in a monitoring network. They developed an objective methodology considering the multiple-criteria, including multiple-pollutants concentration and social factors such as population exposure and the construction cost. The analysis employed atmospheric dispersion model simulations. A multiple-criteria approach in conjunction with the spatial correlation technique was used to develop an optimal air quality monitoring network design. These authors used triangular fuzzy numbers to capture the uncertain (i.e., assigning weights) components in the decision making process. The spatial area coverage of the monitoring station was also determined on the basis of the concept of a sphere of influence.

3. Proposed Methodology

The purpose is to design a multi-pollutant air quality monitoring network for an urban area, considering two objectives: one is the detection of higher pollutant concentrations and the other is the “protection capability” for areas with higher population density. The first one is analysed measuring the potential of a monitoring site to detect violations of reference concentrations in terms of violation scores.

The proposed approach consists of seven steps. The first step is to select the air pollutants of concern and their reference concentration levels for each averaging time less-equal 24h. The values for different intervals of reference concentrations can be chosen based on air quality guideline values for the selected pollutants. Furthermore, weighing factors are defined to penalize the exceedance of higher reference concentrations with regard to exceedance of lower ones.

The second step is to apply atmospheric dispersion models to compute the time series of pollutant concentrations in each grid cell in which the urban area is divided.

In the third step an exceedance score (ESk) of pollutant k is computed for each grid cell. ESk is given by the following equation (Modak & Lohani, 1985b):

ESk=i=1Nkj=1nk(ωj+1-ωj)(Ci,k-CRj,k)Z(CR(j+1),k-CRj,k)
(1)

where Ci,k is a simulated concentration value of pollutant k, Nk is the number of concentration values (Ci,k) of pollutant k, j is the weighing factor corresponding to the reference value CRj,k, nk is the number of reference values for each pollutant, Z is a factor defined by

Z{1Ci,kCRj,k0Ci,kCRj,k
(2)

The fourth step is to evaluate a population factor (PF) for each grid cell, defined by

where P is the number of inhabitants in the grid cell, PT is the total population in the urban area.

In the fifth step the total score (TS) defined by Equation (4) is evaluated for each grid cell.

TS=PFk=1MESk
(4)

where M is the number of pollutants (if one pollutant has more than one averaging time, each of them has to be considered separately).

In the sixth step the grid cells are ranked according to TS values. The location with the maximum TS value is selected as the best monitoring location. All grid squares located nearer than a given distance (D) (for example, 1 km) to the selected one, are discarded for further site selections. The next site locations are determined according to the same procedure. The number of locations is arbitrary (usually limited by the economical constraint). All grid cells with high TS separated more than distance a D are selected for installing a monitoring station. These selected grid cells constitute the preliminary network.

The seventh step is to evaluate if average concentrations of each pollutant at near selected sites are significantly different. Considering one pollutant at a time, and using the t-Student test, if the difference between mean concentrations at a pair of near sites is statistically significant at the 99% confidence level, both sites remain in the network. Otherwise, the site with less TS can be eliminated from the preliminary network. This procedure is repeated considering all sites. The proposed network is obtained in this step.

Furthermore, “the spatial representativeness” of the monitoring sites of the proposed network can be evaluated. Applying the t-Student test to each pollutant mean concentration, “the spatial representativeness” of each monitoring site can be given by all the near grid cells where mean concentrations are not statistically significant different at the 99% confidence level.

4. Application to the city of Buenos Aires

4.1. The city of Buenos Aires and its surroundings

The city of Buenos Aires (34 35’S – 58 26’W) is the capital of Argentina and is located on the west coast of the de la Plata River. It has an extension of 203km2 and 3058309 inhabitants (INDEC, 2008). The city (Figure 1) is surrounded by the Greater Buenos Aires (24 districts) of 3627km2 and 9575955 inhabitants. Both the city of Buenos Aires and the Greater Buenos Aires form the Metropolitan Area of Buenos Aires (MABA), which is considered the third megacity in Latin America following Mexico City (Mexico) and Sao Paulo (Brazil).

media/image5.jpeg

Figure 1.

Location of the city of Buenos Aires and an aerial view of Buenos Aires city.

The MABA is located on a flat terrain with height differences less than 30 m. The de la Plata River is a shallow estuary of 35000km2, approximately. It is 320km length, and its width varies from 38km to 230km in the upper and lower regions, respectively. In front of the city, the width of the river is about 42km. Mean water temperature varies from 12 C in winter to 24 C in summer. The de la Plata River plain has a temperate climate. The city is hot and humid during summer months (December to February), with a mean high temperature of 27 C. Fluctuating temperatures and quickly changing weather conditions characterise autumn and spring seasons. The winter months (June to August) are mild but humid, with a mean minimum temperature of 6 C. The annual average temperature is 18 C in the city, and it varies between 15-16 C in the suburbs. In the city, frost may occur from June to August, but snowfall is extremely rare. The annual rainfall varies between 900mm and 1600mm, influenced by winds that advect humidity from the Atlantic Ocean. Rainfall is heaviest in March. Winds are generally of low intensity. Strong winds are more frequent between September and March, when storms are more frequent. Annual frequency of winds blowing clean air from the river towards the city is 58%.

The air quality in the city has been the subject of several studies carried out during the last years. Some of these studies analysed data obtained from measurement surveys of pollutants in the urban air (Bogo et al., 1999, 2001, 2003; Venegas & Mazzeo, 2000, 2003b; Mazzeo & Venegas, 2002, 2004; Mazzeo et al., 2005; Bocca et al., 2006). Other studies reported results of the application of atmospheric dispersion models (Venegas & Mazzeo, 2005, 2006). In the Greater Buenos Aires, very few air quality measurements have been made (Fagundez et al., 2001, SAyDS, 2002).

4.2. Emission inventory for the city of Buenos Aires

Mazzeo & Venegas (2003) developed a first version of CO and NOx (expressed as NO2) emission inventory for Buenos Aires city. Also Pineda Rojas et al. (2007) presented an emission inventory of these pollutants for the Metropolitan Area of Buenos Aires which includes updated emissions for the city of Buenos Aires. An emission inventory of particulate matter (PM10) for the city of Buenos Aires has been presented by Venegas & Martin (2004). The inventories for the city of Buenos Aires include: a) area sources: residential, commercial, small industries, aircrafts LTO (landing/take-off) cycles at the domestic airport, and road traffic (cars, trucks, taxis, buses) and b) point sources: stacks of three Power Plants. The spatial resolution of the inventories is 1x1 km and a typical hourly variation. The emission factors used in preparing the emission inventories were derived considering: a) monitoring studies undertaken in Buenos Aires (Rideout et al., 2005; b) the EMEP/CORINAIR Atmospheric Inventory Guidebook (EMEP/CORINAIR, 2001); c) the US Environmental Protection Agency’s manual on the Compilation of Air Pollution Emission Factors (EPA, 1995). These factors were applied to fuel consumption, gas supply data and vehicle kilometres travelled within each grid square. Data on traffic flow, fleet composition and bus service frequencies was also available. Aircraft emissions were computed knowing the scheduled hourly flights, the type of aircraft, the information available on LTO (landing/take-off) cycles and emission factors (Romano et al, 1999, EMEP/CORINAIR, 2001). Spatial and temporal dependent NOx (expressed as NO2), CO and PM10 emission distributions in the Buenos Aires Metropolitan Area were obtained.

Figure 2, 3 and 4 show in diagrammatic form the distribution of annual emission of NOx (expressed as NO2), CO and PM10 by source category, for the city of Buenos Aires. Since the Power Plants burn natural gas most of the year and consume fuel oil as much as twenty days in wintertime only, they are responsible for approximately, 51.6% of NOx (expressed as NO2), 0.02% of CO and 2.3% of PM10 total annual emissions in the city. Carbon monoxide and particulate matter are very much associated with emissions from road traffic. Within Buenos Aires city, road traffic is responsible for 43% of NOx (expressed as NO2), 99.43% of CO and 94% of PM10 annual emissions in the city.

media/image6.jpeg

Figure 2.

Estimated emissions of NOx (expressed as NO2) by source category. (*) cars, taxis, trucks.

media/image7.jpeg

Figure 3.

Estimated emissions of CO by source category. (*) cars, taxis, trucks.

media/image8.jpeg

Figure 4.

Estimated emission of PM10 by source category. (*) cars, taxis, trucks.

4.3. Brief description of the atmospheric dispersion models used in this study

In this work, area sources and point sources contributions to air pollutant concentrations in the city are estimated applying the following atmospheric dispersion models: a) DAUMOD model for urban area sources and b) AERMOD model for urban point sources.

The DAUMOD model

The DAUMOD urban atmospheric dispersion model can be used to estimate background concentrations due to area source emissions. This model has been described elsewhere (Mazzeo & Venegas, 1991; Venegas & Mazzeo, 2002, 2006) however a brief description of it is included. In the development of the DAUMOD model (Mazzeo & Venegas, 1991), a semi-infinite volume of air, bounded by the planes z=0 and x=0 is considered. According with Gifford (1970), for steady-state conditions, with the x-axis in the direction of the mean wind and the z-axis vertical, the concentration [C(x,z)] of pollutants emitted from an area source at the surface, satisfies the bi-dimensional diffusion equation :

u(z)C(x,z)x=z[K(z)C(x,z)z]
(5)

because only the vertical diffusion contributes significantly to the process. In Equation (5), u(z) is the mean wind speed and K(z) is the vertical eddy diffusivity for contaminants. If the effluents are emitted continuously from the surface level with source strength (Q) expressed as mass per unit area per unit time, the concentration [C(x,z)] satisfies the lower boundary condition:

K(z) C(x,z)z|z=0=Q
(6)

Another basic assumption is that at a given distance, the vertical extension of the plume of contaminants is h(x), where C(x,h(x)) = 0. Then, there is no transport of mass through the upper limit of the plume, and this can be expressed as:

K(z) C(x,z)z|z=h=0
(7)

The boundary condition C(x,h(x))= 0, can be satisfied assuming that the solution of Equation (5) is given by the following polynomial form:

C(x,z)=C(x,0)j=06Aj(zh)j
(8)

The coefficients Aj have been computed fitting Equation (8) to the results given by the following expression (Pasquill & Smith, 1983):

C(x,z)=C(x,0)exp[4.605(zzm)s]
(9)

where s is a shape factor which depends on atmospheric stability and surface roughness (Gryning et al., 1987) and zm is the height at which C(zm) = 0.01C(0). The height zm is usually considered to be the upper limit of the plume, so we assumed h= zm. Considering different atmospheric stability conditions, the coefficients (A0, A1,….A6) of the polynomial of grade 6 are obtained for each fitting.

Atmospheric stability Expressions for Aj(z0/L)
z0/L < -10-2 Aj = Aj(z0/L=-0.01)
-10-2 ≤ z0/L < -10-4 A0 = 1.0
A1 = -9.254667- 0.8043134 ln(|z0/L|)
A2 = -26.88303107 – 197.989893 [ln(1.2146 |z0/L|)]-1
A3 = -38.00005 + exp[4.16612 - 373.1065 |z0/L|)]
A4 = -84.48740174 – 333.915544 [ln(7.5651|z0/L|)]-1
A5 = -33.25054 + exp[4.13875 - 289.5308 |z0/L|)]
A6 =-14.47563571 – 43.4735075 [ln(14.5776 |z0/L|)]-1
-10-4 ≤ z0/L ≤ 10-4 A0 = 1.0
A1 = 3853.3 (z0/L) - 1.461
A2 = -18740 (z0/L) - 6.797
A3 = 27740 (z0/L) + 26.931
A4 = -16270 (z0/L) - 39.652
A5 = 965 (z0/L) + 27.781
A6 = 1635 (z0/L) – 7.655
10-4< z0/LA0 = 1.0
A1= -0.05478233 - 0.0001021171 [ln[(z0/L)+1]]-1
A2 = -6.55023478 + 0.02035983 [ln(z0/L)]3 +0.00191583[ln(z0/L)]4
A3 = 12.9282233 + exp[2.917612 -1007.8064 (z0/L)]
A4 =-0.59677391 + 0.05583574 [ln(z0/L)]3 +0.00040899[ln(z0/L)]4
A5 = -1.9551195 + exp[3.5211141 -1255.2843 (z0/L)]
A6 = 2.66883478 + 0.00810494[ln(z0/L)]3 -0.00053199[ln(z0/L)]4

Table 1.

Expressions of Aj(z0/L) as functions of (z0/L).

The fittings of polynomial form given by Equation (8) to the results of Equation (9) are excellent; with coefficients of determination of ≈ 1.0 (the reader can find details of these results in Mazzeo & Venegas, 1991). Coefficients Aj (j=0,…6) depend on surface roughness and atmospheric stability. Plotting the values of Aj vs (z0/L) (z0 is the surface roughness length and L is the Monin-Obukhov´s length) the expressions of Aj (z0/L) included in Table 1 have been obtained

The following expressions are considered for the wind speed and the eddy diffusivity (Arya, 1999):

u(z)=u*kv[ln(zz0)+Ψ(zL)]
(10)
K(z)=kvu*(z+z0)φ(zL)
(11)

where u* is the friction velocity, kv is the von Karman´s constant (=0.41), (z/L) functions determine stability correction due to stratification and (z/L) is the dimensionless wind shear (Wieringa, 1980; Gryning et al., 1987):

Ψ(zL)={6.9zLzL01φ1(zL)zL0

φ(zL)={1+6.9zLzL0(122zL)14zL0

(0)= 0 and (0)= 1.

Substituting Equations (8) and (11) into Equation (6) and operating, the following expression for C(x,0) can be obtained:

C(x,0)=Qh(x)|A1|kvu*z0
(12)

C(x,0) can be estimated from Equation (12) knowing the form of h(x). Therefore, considering the equation of pollutant mass continuity expressed by (Pasquill & Smith, 1983):

0xQ dx=zohu(z)C(x,z) dz
(13)

and Equations (8), (10) and (12) along with the boundary condition C(x,h)=0, the following expression can be obtained when Q= constant:

xz0=1|A1|kv2(hz0)2G(z0h;hL)
(14)

The form of G(z0/h; h/L) is not simple (the complete expression is included in Mazzeo & Venegas, 1991), however the values of (h/z0) computed from Equation (14) can be fitted with great accuracy (coefficient of determination ≈1) to potential functions (Mazzeo & Venegas, 1991) given by:

a and b depend on atmospheric stability. The expressions for a(z0/L) and b(z0/L) are included in Table 2.

Atmospheric stabilityExpressions for a(z0 /L) and b(z0 /L)
z0/L < -10-4a= 3.618833 + 0.2369076 ln(|z0/L |)
b= 0.5356147 + 0.0234187 ln [(|z0/L| ) + 0.01]
-10-4 ≤ z0/L ≤ 10-4a= -384.73 (z0 /L) + 1.4
b= -130.0 (z0 /L) + 0.415
10-4 < z0/La= 0.6224632 + 7.37387x10-5 / ln [(z0/L) + 1 ]
b= 0.5065736 - 1.196137 / ln [2802.315 + 9 / (z0 /L) ]

Table 2.

Expressions of a(z0/L) and b(z0/L) as functions of (z0/L).

Substituting Equation (15) in Equation (12), it becomes:

C(x,0)=aQ(x/z0)b|A1|kvu*
(16)

which is the expression for a semi-infinite area source emitting continuously with a uniform strength (Q). The expression for a finite and continuous source located between x= 0 and x= x1, with strength Q, can be derived from Equation (16) by subtracting the concentration due to a continuous semi-infinite source with strength Q lying through x>x1, from Equation (16):

C(x,0)=aQ[xb(xx1)b]|A1|kvz0bu*
(17)

In an urban area, we may assume horizontal distribution of area sources with strength varying according to a typical square grid pattern. Each grid square has a uniform strength Qi (i = 0, 1, 2, …, N). The variation of the concentration with x, for x > xi (i = 0, 1, 2, …, N) is given by:

C(x,0)=a[Qoxb+i=1N(QiQi-1)(xxi)b](|A1|kvzobu*)
(18)

The form of C(x,z) can be obtained substituting Equation (18) into Equation (8),

C(x,z)=a[Qoxb+i=1N(QiQi-1)(xxi)b](|A1|kvzobu*)j=06Aj(zh)j
(19)

A constant wind direction is required for application of Equations (18) and (19). It has been noted from applications of Equation (18) that estimated concentration at any receptor is mainly originated from the emission in the grid square in which the receptor is located. This is because area source distributions in a city are generally fairly smooth and, the contribution of upstream grid squares (from Equation (18)) rapidly reduces with distance to the receptor. The simplification of assuming that the uniform area source strength Qi only varies with x (in the wind direction), suppose to consider the “narrow plume” hypothesis. This assumption has also been included in other simple urban dispersion models (Gifford, 1970, Gifford & Hanna, 1973, Arya, 1999).

The performance of DAUMOD model in estimating background concentrations has been evaluated comparing estimated and observed concentration data from several cities. Results for Bremen (Germany), Frankfurt (Germany) and Nashville (USA) have been reported in Mazzeo & Venegas (1991) and for Copenhagen (Denmark) can be found in Venegas & Mazzeo (2002). The validation of the application of DAUMOD to estimate NOx, CO and PM10 in Buenos Aires City can be found in Mazzeo & Venegas (2004 ), Venegas & Mazzeo (2006) and Venegas & Martin (2004). Results show that the performance of the model in estimating short-term concentrations (hourly and daily) is good and it improves when estimating long averaging time values (monthly and annual).

At present, photochemical transformations involving NO, NO2 and O3 are not included in DAUMOD model. However, output concentrations of NO2 are calculated on the basis of an empirical relationship between NO2 and NOx (Derwent & Middleton, 1996; Dixon et al., 2001; Middleton at al., 2008). The concentration of NO2 is calculated using the polynomial expression (Derwent & Middleton, 1996, CERC, 2003):

[NO2] = 2.166 – [NOx] (1.236 – 3.348 B + 1.933 B2 – 0.326 B3)

where B= log10([NOx]) and [NOx] is hourly-averaged concentration in ppb.

The AERMOD model

AERMOD (EPA, 2004) is a steady-state plume dispersion model for assessment of pollutant concentrations from a variety of sources. AERMOD simulates transport and dispersion from multiple point, area or volume sources based on an up-to-date characterization of planetary boundary layer (PBL). AERMOD’s concentration formulations consider the effects from vertical variations in wind, temperature and turbulence. These profiles are represented by “effective” values constructed by averaging over the layer through which plume material travels directly from the source to receptor. Sources may be located in rural or urban areas, and receptors may be located in simple or complex terrain. The model employs hourly sequential pre-processed meteorological data to estimate concentrations for averaging times from one hour to one year (also multiple years). AERMOD is designed to operate in concert with two pre-processor codes: AERMET processes meteorological data for input to AERMOD, and AERMAP processes terrain elevation data and generates receptor information for input to AERMOD. AERMOD is applicable to primary pollutants and continuous releases of toxic and hazardous waste pollutants. Chemical transformation is treated by simple exponential decay. A more complete description of AERMOD model can be found in EPA (2004) and Cimorelli et al. (2005).

In stable boundary layer (SBL), the concentration distribution is assumed to be Gaussian, both vertically and horizontally. During stable conditions (i.e., stable and neutral stratifications, when L > 0), AERMOD estimates concentrations (Cs) from:

Cs(x,y,z)=Qp u ˜ σzsFy                     xm=-{exp[(zhes2mzieff)2zs2]+exp[(z+hes+2mzieff)2zs2]}
(20)

where, QP is the point source emission rate, ũ is the effective wind speed, zieff is the effective mechanical mixing height, zs is the total vertical dispersion, hes is the plume height (Weil, 1988; Cimorelli et al, 2005) and Fy is the lateral distribution functions.

In the convective boundary layer (CBL), the horizontal distribution is assumed to be Gaussian, but vertical distribution is described with a bi-Gaussian probability density function (Willis & Deardoff, 1981; Briggs, 1993). In CBL the transport and dispersion of a plume is characterized as the superposition of three modelled plumes: the direct plume (from the stack), the indirect plume, and the penetrated plume, where the indirect plume accounts for lofting of a buoyant plume near the top of boundary layer, and the penetrated plume accounts for the portion of a plume that, due to its buoyancy, penetrates above the mixed layer, but can disperse downward and re-enter the mixed layer. In the CBL, plume rise is superposed on the displacements by random convective velocities (Weil et al, 1997).

The total concentration (Cc) in the CBL is found by adding the contribution from three sources: a “direct” source, an “indirect” source and a “penetrated” source (above the CBL top). For the horizontal plume state,

Cc(x,y,z)=Cd(x,y,z)+Cr(x,y,z)+Cp(x,y,z)
(21)

The concentration (Cd) due to the direct plume is given by

Cd(x,y,z)=QPfp2πu˜σyexp[y22σy2]xj=12m=0λjσzj{exp[(zΨdj2mzi)22σzj2]+exp[(z+Ψdj+2mzi)22σzj2]}
(22)

where Ψdj= hs + Δhd + wj¯ x/ũ is the height of the direct source plume, ũ is the effective wind speed, Zi is the mixing height, hd is the plume rise and wj¯ is the vertical velocity. The subscript j is equal to 1 for updrafts and 2 for downdrafts with j defined as the weighting coefficient for each distribution:

λ1=w2¯w2¯w1¯
(23)
λ2=w1¯w2¯w1¯
(24)

Equation (22) uses an image plume to handle ground reflections by assuming a source at z = -hs. All subsequent reflections are handled by sources at z=2zi+hs, -2zi - hs, 4zi + hs, -4zi-hs and so on. The factor fp accounts for the fraction of source material that does not penetrate the top of the CBL.

The concentration (Cr) due to the indirect source is calculated from

Cr(x,y,z)=QPfp2πu˜σyexp[y22σy2]xj=12m=0λjσzj{exp[(z+Ψrj2mzi)22σzj2]+exp[(zΨrj+2mzi)22σzj2]}
(25)

where zj and Fy are the same as defined for the direct source, the plume height Ψrj= hs + Δhr + wj¯ x/ũ (j=1,2) with hr= hd - hi and hi is the indirect source plume rise.

For the penetrated source, the vertical and lateral concentration distributions have a Gaussian form, such that the concentration (Cp) is given by

Cp(x,y,z)=QP(1fp)2πu˜σypσzpexp[y22σyp2]xm={exp[(zhep2mzieff)22σzp2]+exp[(z+hep+2mzieff)22σzp2]}
(26)

where hep is the height of the penetrated plume height and zieff is the height of the upper reflecting surface in a stable layer.

For flow in complex terrain, AERMOD incorporates the concept of a dividing streamline (Snyder et al., 1985), and the plume is modelled as a combination of terrain-following and terrain-impacting states. This approach has been designed to be physically realistic and simple to implement while avoiding the need to distinguish among simple, intermediate and complex terrain. As result, AERMOD removes the need for defining complex regimes; all terrain is handled in a consistent and continuous manner that is simple while still considering the dividing streamline concept in stably-stratified conditions.

The model considers the influence of building wakes and it enhances vertical turbulence to account for the “convective like” boundary layer found in night-time urban areas.

4.4. Annual concentration distributions in the city of Buenos Aires

Hourly urban background concentrations of NO2, CO and PM10 at each grid square (1 x 1 km) in which the city of Buenos Aires has been divided, are estimated applying the two atmospheric dispersion models described above. DAUMOD model is applied to compute area source emissions contribution and AERMOD to estimate point sources contribution. Hourly data registered at the meteorological station of the domestic airport (located in the city) and the emissions inventory for Buenos Aires city (Section 4.2.) have been used in calculations. Figure 5, 6 and 7 show the horizontal distributions of computed mean annual NO2, CO and PM10 concentrations within the city of Buenos Aires. High concentration values can be found downtown and near the main train stations in the city, where there are usually more activity and traffic. The NO2, CO and PM10 concentration patterns show large spatial variability that is mainly related to the distribution of area sources emissions. Higher concentration values are related to higher traffic density area.

media/image36.jpg

Figure 5.

Mean annual background concentrations of NO2 in Buenos Aires city.

media/image37.jpg

Figure 6.

Mean annual background concentrations of CO in Buenos Aires city.

media/image38.jpg

Figure 7.

Mean annual background concentrations of PM10 in Buenos Aires city.

4.5. Application of the methodology for designing the air quality monitoring network

The site selection proposed methodology described in Section 3 has been applied to determine the site locations for monitoring 1-h average NO2 concentration (k=1), 1-h (k=2) and 8-h (k=3) average CO concentrations and 24-h average PM10 concentration (k=4) in the atmosphere of the city of Buenos Aires, greater than reference concentration levels.

Reference concentration levels have been selected considering the Air Quality Guidelines established by the World Health Organization (W.H.O., 2000, 2006). A weighing factor (ωj) of 1 is assigned to the Air Quality Guideline levels. Table 3 presents the reference concentration levels and the corresponding weighing factors for the three pollutants considered. According to the proposed methodology, in this application, M=4 and nk=7.

Weighing factor
(ωj)
NO2
(averaging
time: 1h)
CO
(averaging
time: 1h)
CO
(averaging
time: 8h)
PM10
(averaging
time: 24h)
mg/m3mg/m3mg/m3mg/m3
0.50.101550.025
0.70.142170.035
0.90.182790.045
1.00.2030100.050
1.20.2436120.060
1.50.3045150.075
2.00.4060200.100

Table 3.

Reference concentration levels and weighing factors (j) for NO2 (averaging time: 1h), CO (averaging times: 1h and 8h) and PM10 (averaging time: 24h).

The values of 1-h average NO2 concentrations, 1-h and 8-h average CO concentrations and 24-h average PM10 concentrations (Ci,k) are estimated for each grid square in which the city has been divided applying DAUMOD and AERMOD atmospheric dispersion models.

Using model estimations of (Ci,k) and data on Table 3, the values of the exceedance score ESk (Equation (1)) for each combination (air pollutant, averaging time) (k=1,...,4) are estimated for each grid square. Figure 8 presents the horizontal distribution of the total exceedance score (ES= ES1 + ES2 + ES3 + ES4) in the urban area.

In order to evaluate the population factor (PF) given by Equation (3), the horizontal distribution of population density (inhab/km2) in the city has been elaborated using the information of the last National Census (INDEC, 2008). Figure 9 shows the horizontal distribution of the population factor (PF) in the city of Buenos Aires.

media/image39.jpg

Figure 8.

Horizontal distribution of the total exceedance score ( ES=k=14ESk ) in the city of Buenos Aires.

media/image41.jpg

Figure 9.

Horizontal distribution of the population factor (PF) in the city of Buenos Aires.

In this way, knowing the total exceedance score (ES) and the population factor (PF) for each grid square, the horizontal distribution of the total score (TS= PFxES) (Figure 10) is estimated using Equation (4). A preliminary monitoring network configuration is defined considering a budget constraint that limits the number of monitoring sites to twelve. The location of the 12 stations is obtained maximizing TS and considering a minimum distance of D=1km between two sensors of the same pollutant. Figure 11 shows the monitoring sites of the preliminary network.

media/image42.jpg

Figure 10.

Horizontal distribution of the Total Score (TS) in the urban area.

media/image43.jpg

Figure 11.

Monitoring locations selected considering the total score and the distance constraint (preliminary configuration).

The final configuration of the proposed monitoring network is defined comparing the mean concentration of each pollutant at near preliminary sites, using the t-test. If the difference between average concentrations at two nearby sites is statistically significant at the 99% confidence level, both sites remain in the final network configuration. Otherwise, the site with lower TS can be discarded. Table 4 and 5 show if the difference between average concentrations at a pair of near sites is statistically significant at the 99% (Yes:Y, No: N).

NO2123456789101112CO123456789101112
1YYYYYY1YYYYYY
2YYYY2YYYY
3YNY3YNN
4NYYYY4YYYYY
5YYYYY5YYYYY
6NY6NN
7NY7NY
8Y8Y
9YY9YY
1010
11N11Y

Table 4.

It is indicated if the difference between the average concentrations at two nearby monitoring locations is statistically significant at the 99% confidence level (Y: yes; N: no).

PM10123456789101112
1NYYYYY
2YYYY
3NYY
4NYYYY
5NYYYY
6YN
7NY
8Y
9YY
10
11N

Table 5.

It is indicated if the difference between the average concentrations at two nearby monitoring locations is statistically significant at the 99% confidence level (Y: yes; N: no).

Figure 12 shows the proposed monitoring sites and the pollutants to be measured at each site according to results on Table 4 and 5.

media/image44.jpg

Figure 12.

Sites chosen to locate an air quality monitoring station and pollutants to be measured.

Once the proposed network (Figure 12) is in operation, the environmental authority of the city may be interested to know the horizontal extension of the “spatial representativeness” of mean concentration values registered at each monitoring site. Figure 13, 14 and 15 shows the areas near each monitoring site where the NO2, CO and PM10 mean concentrations cannot be considered statistically significant different at the 99% confidence level, obtained after the application of the t-Student test.

media/image45.jpg

Figure 13.

“Spatial representativeness” for NO2 monitoring stations (▲).

media/image46.jpg

Figure 14.

“Spatial representativeness” for CO monitoring stations (●).

media/image47.jpg

Figure 15.

“Spatial representativeness” for PM10 monitoring stations (▼).

5. Conclusion

A multiple objective and multi-pollutant planning procedure for designing an urban air quality monitoring network is presented in this study. The considered monitoring objectives are to maximize the “detection capability” of higher pollutant concentrations and the “protection capability” for areas with higher population density. The design methodology is based on the analysis of air pollutant concentrations estimated by atmospheric dispersion models. It simultaneously considers an exceedance score and a population factor. A statistical analysis is used for optimization.

The methodology is applied to design a NO2, CO and PM10, monitoring network for the city of Buenos Aires considering a spatial resolution (for the emission inventory and model estimations) of 1 x 1km. Air pollutant concentrations in the city have been estimated using the DAUMOD and AERMOD atmospheric dispersion models, that evaluate the contribution of area and point sources, respectively.

The optimal alternative of the proposed network can be summarized as: a) seven locations for monitoring NO2, CO and PM10; b) two sites for NO2 and CO; c) one location for CO and PM10 and d) one station for NO2 only. The “spatial representativeness” of mean concentrations at monitoring sites varies with each pollutant: a) for NO2, between 1-7km2; b) for CO between 2-11km2 and c) for PM10, between 4-12km2.

It must be noted that the ultimate decision in site selection is left to the air quality monitoring authority.

Future studies could be focused on: a) the evaluation of how sensitive is the proposed methodology for air quality network design to slight changes in the input data (e.g. the weighing factors, the spatial resolution) and b) the inclusion of other optimization objectives (e.g. land use, frequency of violations of air quality standards, protect damage to vulnerable receptors as historic and/or artistic valuable property).

References

1 - B. Ainslie, C. Reuten, D. G. Steyn, N. D. Le , J. V. Zidek, 2009 Application of an entropy-based Bayesian optimization technique to the redesign of an existing monitoring network for single air pollutants. Journal of Environmental Management, 90 2715 2729 .
2 - S. P. Arya, 1999 Air Pollution Meteorology. Oxford University Press. New York
3 - R. Baldauf, R. W. Wiener, D. K. Heist, 2002 Methodology for Siting Ambient Air Monitors at the Neighborhood Scale. Journal of Air & Waste Management Association, 52 1433 1442 .
4 - B. Bocca, S. Caimi, P. Smichowski, D. Gómez, S. Cairoli, 2006 Monitoring Pt and Rh in urban aerosols from Buenos Aires, Argentina. Science of the Total Environment, 358 255 264 .
5 - H. Bogo, R. M. Negri, E. San Román, 1999 Continuous measurement of gaseous pollutants in Buenos Aires city. Atmospheric Environment, 33 2587 2598 .
6 - H. Bogo, D. R. Gómez, S. L. Reich, R. M. Negri, E. San Román, 2001 Traffic pollution in downtown of Buenos Aires City. Atmospheric Environment, 35 1717 1727 .
7 - H. Bogo, M. Otero, P. Castro, M. J. Ozafrán, A. Kreiner, E. J. Calvo, R. M. Negri, 2003 Study of atmospheric particulate matter in Buenos Aires city. Atmospheric Environment, 37 1135 1147 .
8 - G. A. Briggs, 1993 Plume dispersion in convective boundary layer. Part II: Analysis of CONDORS field experiment data. Journal of Applied Meteorology, 32 1388 1425 .
9 - W. F. Caselton, L. Kan, J. V. Zidek, 1992 Quality data networks that minimize entropy. In: Statistics in the Environmental and Earth Sciences. Walden, A. & Guttorp, P. (Eds.), 10 38 , Halsted Press, New York.
10 - CERC 2003 ADMS-Urban. An Urban Air Quality Management System. User Guide. Version 2.0. Cambridge Environmental Research Consultants Ltd., Cambridge.
11 - A. J. Cimorelli, S. G. Perry, A. Venkatram, J. C. Weil, R. J. Paine, R. B. Wilson, R. F. Lee, W. D. Peters, R. W. Brode, 2005 AERMOD: A dispersion model for industrial source applications Part I: General model formulation and boundary layer characterization. Journal of Applied Meteorology, 44 682 693 .
12 - R. G. Derwent, D. R. Middleton, 1996 An empirical function for the ratio 2 NOx. Clean Air, 26 57 62 .
13 - J. Dixon, D. R. Middleton, R. G. Derwent, 2001 Sensitivity of nitrogen dioxide concentrations to oxides of nitrogen controls in the United Kingdom. Atmospheric Environment, 35 3715 3728 .
14 - N. D. V. Egmond, D. Onderdelinden, 1981 Objective analysis of air pollution monitoring network data: spatial interpolation and network density. Atmospheric Environment, 15 1035 1046 .
15 - A. Elkamel, E. Fatehifar, M. Taheri, M. S. Al-Rashidi, A. Lohi, 2008 A heuristic optimization approach for Air Quality Monitoring Network design with the simultaneous consideration of multiple pollutants. Environmental Management, 88 507 516 .
16 - D. M. Elsom, 1978 Spatial correlation analysis of air pollution data in an urban area. Atmospheric Environment, 12 1103 1107 .
17 - EMEP/CORINAIR. 2001 Atmospheric Emission Inventory Guidebook, Third Edition, European Environment Agency, Copenhagen.
18 - EPA. 1995 Compilation of Air Pollution Emission Factors, AP-42, 5th ed., United States Environmental Protection Agency, Research Triangle Park, NC.
19 - EPA. 2004 User’s Guide for the AMS/EPA Regulatory Model-AERMOD, EPA-454/B-03-001. United States Environmental Protection Agency, Research Triangle Park, NC.
20 - L. A. Fagundez, V. L. Fernández, T. H. Marino, I. Martín, D. A. Persano, M. Rivarola y Benítez, I. V. Sadañiowski, J. Codnia, A. Zalts, 2001 Preliminary air pollution monitoring in San Miguel, Buenos Aires. Environmental Monitoring and Assessment, 71 61 70 .
21 - F. A. Gifford, 1970 Atmospheric Diffusion in an Urban Area, NOAA Research Lab. Nº 33. Oak Ridge, N. C.
22 - F. A. Gifford, S. R. Hanna, 1973 Modelling urban air pollution. Atmospheric Environment, 7 131 136 .
23 - R. J. Graves, T. D. Lee, L. F. J. McGinnis, 1981 Air Monitoring Network Design: case study. Journal of Environmental Engineering ASCE, 107 941 955 .
24 - S. E. Gryning, A. A. M. Footslog, J. S. Irwin, B. Sivertsen, 1987 Applied dispersion modelling based on meteorological scaling parameters. Atmospheric Environment, 21 79 89 .
25 - C. M. Handscombe, D. M. Elsom, 1982 Rationalisation of the National Survey of Air Pollution Monitoring Network of the United Kingdom using spatial correlation analysis: a case study of the Greater London area. Atmospheric Environment, 16 1061 1070 .
26 - E. S. Hougland, N. T. Stephens, 1976 Air Pollutant Monitor Siting by Analytical Techniques. Journal of the Air Pollution Control Association, 26 51 53 .
27 - T. Husain, S. M. Khan, 1983 Air Monitoring Network Design Using Fisher’s Information Measures-A Case Study. Atmospheric Environment, 17 2591 2598 .
28 - J. S. Hwang, Ch. Ch. Chan, 1997 Redundant measurements on urban air monitoring networks in air quality reporting. Journal of Air & Waste Management Association, 47 614 619 .
29 - INDEC 2008 Estimaciones de población total por departamento y año calendario. Período 2001-2010. Serie Análisis Demográfico. Nº 34. Instituto Nacional de Estadística y Censos. (www.indec.gov.ar) (in Spanish)
30 - Y. Kainuma, K. Shiozawa, S. Okamoto, 1990 Study of the Optimal Allocation of Ambient Air Monitoring Stations. Atmospheric Environment, 24B 395 406 .
31 - M. Koda, J. H. Seinfeld, 1978 Air monitoring siting by objective. EPA-600/4-7-036, United States Environmental Protection Agency, Las Vegas, Nevada.
32 - J. Langstaff, C. Seigneur, M. K. Liu, 1987 Design of an optimum air monitoring network for exposure assessment. Atmospheric Environment, 21 1393 1410 .
33 - N. D. Le, J. V. Zidek, 2006 Statistical Analysis of Environmental Space-Time Processes. Springer, New York.
34 - M. K. Liu, J. Avrin, R. I. Pollack, J. V. Behar, J. L. McElory, 1986 Methodology for designing air quality monitoring networks. Environmental Monitoring and Assessment, 6 1 11 .
35 - N. A. Mazzeo, L. E. Venegas, 1991 Air pollution model for an urban area. Atmospheric Research, 26 165 179 .
36 - N. A. Mazzeo, L. E. Venegas, 2000 Practical use of the ISCST3 model to select monitoring site locations for air pollution control. International Journal of Environment and Pollution, 14 246 259 .
37 - N. A. Mazzeo, L. E. Venegas, 2002 Estimation of cumulative frequency distribution for carbon monoxide concentration from wind-speed data, in Buenos Aires (Argentina). Water, Air and Soil Pollution, Focus, 2 419 432 .
38 - N. A. Mazzeo, L. E. Venegas, 2003 Carbon monoxide and nitrogen oxides emission inventory for Buenos Aires City (Argentina). Proceedings of the Fourth International Conference on Urban Air Quality- Measurement, Modelling and Management, 159 162 . Prague, Czech Republic, March 2003. University of Hertfordshire, Hatfield.
39 - N. A. Mazzeo, L. E. Venegas, 2004 Some aspects of air pollution in Buenos Aires city. International Journal of Environment and Pollution, 22 365 378 .
40 - N. A. Mazzeo, L. E. Venegas, H. Choren, 2005 Analysis of NO2 O3 and NOx concentrations measured at a green area of Buenos Aires City during wintertime, Atmospheric Environment, 39 3055 3068 .
41 - N. A. Mazzeo, L. E. Venegas, 2008 Design of an air quality surveillance system for Buenos Aires city integrated by a NOx monitoring network and atmospheric dispersion models. Environmental Modeling and Assessment, 13 349 356 .
42 - J. L. McElroy, J. V. Behar, T. C. Meyers, M. K. Liu, 1986 Methodology for designing Air Quality Monitoring Networks: II. Application to Las Vegas, Nevada, for carbon monoxide. Environmental Monitoring and Assessment, 6 13 34 .
43 - D. R. Middleton, A. R. Jones, A. L. Redington, D. J. Thomson, R. S. Sokhi, L. Luhana, B. E. A. Fisher, 2008 Lagrangian modelling of plume chemistry for secondary pollutants in large industrial plumes. Atmospheric Environment, 42 415 427 .
44 - P. M. Modak, B. N. Lohani, 1985a Optimization of ambient Air Quality Monitoring Networks: Part I. Environmental Monitoring and Assessment, 5 1 19 .
45 - P. M. Modak, B. N. Lohani, 1985b Optimization of ambient Air Quality Monitoring Networks: Part II. Environmental Monitoring and Assessment, 5 21 38 .
46 - P. M. Modak, B. N. Lohani, 1985c Optimization of ambient Air Quality Monitoring Networks: Part III. Environmental Monitoring and Assessment, 5 39 53 .
47 - A. Mofarrah, T. Husain, 2010 A Holistic Approach for optimal design of Air Quality Monitoring Network Expansion in an Urban Area. Atmospheric Environment, 44 432 440 .
48 - R. E. Munn, 1975 Suspended particulate concentrations: Spatial correlations in the Detroit-Windsor area. Tellus, Vol. XXVII, 397 405
49 - R. E. Munn, 1981 The Design of Air Quality Monitoring Networks. MacMillan, London.
50 - Y. Nakamori, Y. Sawaragi, 1984 Interactive Design of Urban Level Air Quality Monitoring Network. Atmospheric Environment, 18 793 799 .
51 - K. E. Noll, S. Mitsutome, 1983 Design methodology for optimum dosage air monitoring site selection. Atmospheric Environment, 17 2583 2590 .
52 - F. Pasquill, F. B. Smith, 1983 Atmospheric Diffusion, John Wiley & Sons, New York.
53 - E. E. Pickett, R. G. Whiting, 1981 The design of cost-effective Air Quality Monitoring Networks. Environmental Monitoring and Assessment, 1 59 74 .
54 - A. L. Pineda Rojas, L. E. Venegas, N. A. Mazzeo, 2007 Emission inventory of carbon monoxide and nitrogen oxides for area sources at Buenos Aires Metropolitan Area (Argentina). Proceedings of the Sixth International Conference on Urban Air Quality, Cyprus, March 2007, University of Hertfordshire, Hatfield.
55 - J. C. M. Pires, M. C. Pereira, M. C. M. Alvim-Ferraz, F. G. Martins, 2009 Identification of redundant air quality measurements through the use of principal component analysis. Atmospheric Environment, 43 3837 3842 .
56 - S. T. Rao, 2009 Environmental Monitoring and Modeling Needs in the 21st Century, EM, October, 3 4 .
57 - G. Rideout, D. Gourley, J. Walker, 2005 Measurement of in-service vehicle emissions in Sao Paulo, Santiago and Buenos Aires. ARPEL Environmental Report Nº25, Environmental Services Association of Alberta, Edmonton.
58 - D. Romano, D. Gaudioso, R. De Lauretis, 1999 Aircraft emissions: a comparison of methodologies based on different data availability. Environmental Monitoring and Assessment, 56 51 74 .
59 - D. A. Sarigiannis, M. Saisana, 2008 Multi-objective optimization of air quality monitoring. Environmental Monitoring and Assessment, 136 87 99 .
60 - SAyDS. 2002 Estudio o línea de base de concentración de gases contaminantes en atmósfera en el área de Dock Sud en Argentina. Agencia de Cooperación Internacional del Japón en Argentina y Secretaría de Desarrollo Sustentable y Política Ambiental, Buenos Aires. www.medioambiente.gov.ar/dock_sud/default.htm (in Spanish)
61 - D. G. Smith, B. A. Egan, 1979 Design of monitor networks to meet multiple criteria, Journal of Air and Waste Management Association, 29 710 714 .
62 - W. H. Snyder, R. S. Thompson, R. E. Eskridge, R. E. Lawson, I. P. Castro, J. T. Lee, J. C. R. Hunt, Y. Ogawa, 1985 The structure of the strongly stratified flow over hills: Dividing streamline concept. Journal of Fluid Mechanics, 152 249 288 .
63 - C. C. Tseng, N. B. Chang, 2001 Assessing relocation strategies of urban air quality monitoring stations by GA-based compromise programming. Environment International, 26 523 541 .
64 - L. E. Venegas, P. B. Martin, 2004 Particulate Matter Concentrations in the City of Buenos Aires, Proceedings of the 14th Congreso Argentino de Saneamiento y Medio Ambiente, Buenos Aires, November 2004, AIDIS-Argentina, Buenos Aires (in Spanish).
65 - L. E. Venegas, N. A. Mazzeo, 2000 Carbon monoxide concentrations in a street canyon at Buenos Aires City (Argentina). Environmental Monitoring & Assessment, 65 417 424 .
66 - L. E. Venegas, N. A. Mazzeo, 2002 An evaluation of DAUMOD model in estimating urban background concentration. Water, Air and Soil Pollution: Focus, 2 433 443 .
67 - L. E. Venegas, N. A. Mazzeo, 2003a Design methodology for background air pollution monitoring site selection in an urban area. International Journal of Environment and Pollution, 20 185 195 .
68 - L. E. Venegas, N. A. Mazzeo, 2003b Air quality in an area of Buenos Aires City (Argentina), Proceedings of the III Congresso Interamericano de Qualidade do Ar, Canoas, Brasil, July 2004, Asociación Interamericana de Ingeniería Sanitaria y Ambiental-Asociacion Brasileira de Engenharia Sanitaria e Ambiental, Seçao, Rio Grande do Sul. (in Spanish).
69 - L. E. Venegas, N. A. Mazzeo, 2005 Application of atmospheric dispersion models to evaluate population exposure to NO2 concentration in Buenos Aires. International Journal of Environment and Pollution, 25 224 238 .
70 - L. E. Venegas, N. A. Mazzeo, 2006 Modelling of urban background pollution in Buenos Aires city (Argentina). Environmental Modelling & Software, 21 577 586 .
71 - L. E. Venegas, N. A. Mazzeo, 2010 An ambient air quality monitoring network for Buenos Aires city. International Journal of Environment and Pollution, 40 184 194 .
72 - J. C. Weil, 1988 Dispersion in the convective boundary layer. In: Lectures on Air Pollution Modeling. Venkatram, A. & Wingaard, J.C. (Eds.), 167 227 , American Meteorological Society, Boston.
73 - J. C. Weil, L. A. Corio, R. P. Brower, 1997 A PDF dispersion model for buoyant plumes in the convective boundary layer. Journal of Applied Meteorology, 36 982 1003 .
74 - W.H.O. 2000 Air Quality Guidelines for Europe, World Health Organization Regional Publications, European Series Nº91, Copenhagen.
75 - W.H.O. 2006 Air Quality Guidelines. Global Update 2005. Particulate matter, ozone, nitrogen dioxide and sulphur dioxide, World Health Organization, Geneve.
76 - J. Wieringa, 1980 A revaluation of the Kansas mast influence on measurements of stress and cup anemometer overspeeding. Boundary-Layer Meteorology, 18 411 430 .
77 - G. E. Willis, J. W. Deardorff, 1981 A laboratory study of dispersion in the middle of the convectively mixed layer. Atmospheric Environment, 15 109 117 .
78 - S. Wu, J. V. Zidek, 1992 An entropy based review of selected NADP/NTN network sites for 1983-86. Atmospheric Environment, 26A 2089 2103 .