Next Article in Journal
Effects of Δ9-Tetrahydrocannibinol (THC) on Obesity at Different Stages of Life: A Literature Review
Next Article in Special Issue
Estimation of the Seasonal Inhaled Deposited Dose of Particulate Matter in the Respiratory System of Urban Individuals Living in an Eastern Mediterranean City
Previous Article in Journal
Precarious Work as Risk Factor for 5-Year Increase in Depressive Symptoms
Previous Article in Special Issue
Using Multisource Data to Assess PM2.5 Exposure and Spatial Analysis of Lung Cancer in Guangzhou, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sampling Low Air Pollution Concentrations at a Neighborhood Scale in a Desert U.S. Metropolis with Volatile Weather Patterns

by
Nathan Lothrop
1,*,
Nicolas Lopez-Galvez
1,2,
Robert A. Canales
3,
Mary Kay O’Rourke
1,
Stefano Guerra
1,4 and
Paloma Beamer
2,4
1
Mel and Enid Zuckerman College of Public Health, University of Arizona, Tucson, AZ 85721, USA
2
School of Public Health, San Diego State University, San Diego, CA 92182, USA
3
Program in Applied Mathematics, University of Arizona, Tucson, AZ 85721, USA
4
Asthma and Airway Disease Research Center, College of Medicine, University of Arizona, Tucson, AZ 85721, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(6), 3173; https://doi.org/10.3390/ijerph19063173
Submission received: 14 January 2022 / Revised: 3 March 2022 / Accepted: 4 March 2022 / Published: 8 March 2022
(This article belongs to the Special Issue Frontiers in Air Pollution and Human Health)

Abstract

:
Background: Neighborhood-scale air pollution sampling methods have been used in a range of settings but not in low air pollution airsheds with extreme weather events such as volatile precipitation patterns and extreme summer heat and aridity—all of which will become increasingly common with climate change. The desert U.S. metropolis of Tucson, AZ, has historically low air pollution and a climate marked by volatile weather, presenting a unique opportunity. Methods: We adapted neighborhood-scale air pollution sampling methods to measure ambient NO2, NOx, and PM2.5 and PM10 in Tucson, AZ. Results: The air pollution concentrations in this location were well below regulatory guidelines and those of other locations using the same methods. While NO2 and NOx were reliably measured, PM2.5 measurements were moderately correlated with those from a collocated reference monitor (r = 0.41, p = 0.13), potentially because of a combination of differences in inlet heights, oversampling of acutely high PM2.5 events, and/or pump operation beyond temperature specifications. Conclusion: As the climate changes, sampling methods should be reevaluated for accuracy and precision, especially those that do not operate continuously. This is even more critical for low-pollution airsheds, as studies on low air pollution concentrations will help determine how such ambient exposures relate to health outcomes.

1. Introduction

Exposure to ambient air pollutants has long been implicated in numerous health effects [1,2,3], especially respiratory outcomes; however many long-running birth cohorts lack exposure data, necessitating modeling with other air pollution measures [4,5]. Similar to other birth cohorts, the Tucson Children’s Respiratory Study has followed subjects since birth in 1980–1984 but has no air pollution exposures. To take advantage of this long-running health dataset, we sought to measure air pollution at a neighborhood scale to ultimately model exposures for health effect analyses.
The Tucson, Arizona, metropolitan area (Tucson), in the desert southwestern U.S. is home to 1 million people [6] and is defined by a car-centric, low-density, suburban morphology common in the western U.S. While the ambient concentrations of most air pollutants are historically low, Tucson is prone to high naturally occurring, windblown dust concentrations and a climate marked by volatile precipitation patterns, extreme aridity and heat, and dust storms [7,8,9]. With climate change, meteorology around the world will become more extreme [10,11,12,13], mimicking aspects of Tucson’s current climate. There are few studies assessing measurements in arid environments [14,15] and none, to our knowledge, in low-pollution airsheds. Low-pollution locations offer the chance to study threshold effects of pollutant concentrations below regulatory levels—a potential preview of health effects if air pollution is reduced. This presents an opportunity to examine measurement methods at a neighborhood scale. Past studies have used a variety of sampling methods to measure pollutants with intra-city spatial precision, ranging from passive, low-cost samplers (e.g., Ogawa badge samplers) to active monitors (e.g., chemiluminescence monitors) [16,17,18,19,20,21,22,23]. Given the wide range of air sampling approaches, we emulated well-documented protocols implemented in a variety of settings: the European Study of Cohorts for Air Pollution Effects (ESCAPE).
ESCAPE methods for all aspects of monitoring, including location designation, setup, sample analysis, and later model development, have been used in 36 locations among 16 countries throughout Europe, ranging from bucolic to densely urban, from Stockholm in Sweden to the Grecian isle of Crete [24,25]. More recently, these methods were implemented in Taipei, Taiwan [26], in the coastal cities of Perth [27,28] and Brisbane [29] in Australia, and in Durban, South Africa [30]. However, these methods have not been adapted in locations with low air pollution concentrations and more extreme meteorological patterns. In this paper, we illustrate how we adopted and adapted ESCAPE air pollution measurement methods for ambient nitrogen dioxide (NO2), oxides of nitrogen (NOx), and particulate matter < 2.5 µm (PM2.5) and <10 µm (PM10) in diameter, and then assess and compare our findings to ESCAPE and other relevant studies.

2. Materials and Methods

The Tucson, AZ, metropolitan region (Tucson), with an area of 5275 km2, is located in the eastern half of Pima County in southern Arizona, approximately 100 km north of the U.S.–Mexico border (Figure 1). It is defined by the City of Tucson, smaller surrounding jurisdictions, and unincorporated residential development. While the region has a mean density of 190 persons/km2, the City of Tucson (pop. = 535,677) has 809 persons/km2 [6]. Residential development largely comprises detached, single-family dwellings, with little mixed land use. NO2, NOx, and PM2.5 are byproducts of combustion from industry and vehicle travel, as well as from reactions in the atmosphere. While Tucson lacks heavy industry and large point sources of pollution, Tucson has an extensive street network for personal vehicle travel as the primary mode of transportation, which generates area source pollution, much like other urban areas in the western U.S. Nevertheless, these pollutants have been perpetually below relevant ambient air quality standards [7].
Increased concentrations of PM10 in the study area often come from a combination of high winds, limited precipitation, and increased aridity [7]. Together, these factors result in less PM10 and PM2.5 being scoured from the air, less hygroscopic particulate growth [31];,and less dust-arresting vegetation cover [32]. In the early summer, daytime high temperatures regularly reach 43 °C, while snowfall can be expected at least once over winter [8]. While Tucson receives an average total of 30 cm of precipitation, nearly all rainfall comes in the form of violent, but short-lived, summer monsoons or steady, gentler winter rains [8]. In the future, Tucson’s annual precipitation is expected to be reduced even further compared to that in other parts the western U.S. [33]. In addition, dust storms or “haboobs” are more frequent, particularly during the summer monsoon season [9].
Following methods from ESCAPE [24,25,34], we aimed to sample pollutants at 40 sites within the Tucson, AZ, metropolitan area to ensure sufficient geographic representation while balancing the feasibility of sampler set up and material costs. Sites were divided into three distinct types, defined by their expected sources of air pollution: regional background, urban background, and street (see Supplementary Materials for details). Sites were recruited to ensure that half of the locations of each site type measured NO2 and NOx and the other half, NO2, NOx, PM2.5, and PM10. In Tucson, there are few building canyons. As a result, unlike ESCAPE, in which sites could include balconies, our study utilized only ground-level locations.
Sites were selected from residential areas where two birth cohorts lived by overlaying cohort addresses with vehicle traffic maps using ArcGIS ArcMap 10.4 (ESRI, Redlands, CA, USA). We iteratively recruited via email through multiple listservs at the Mel and Enid Zuckerman College of Public Health and the Asthma and Airways Disease Research Center at the University of Arizona, as well as with word-of-mouth recruitment among existing participants and study personnel. As recruitment progressed, later emails solicited residents in specific areas where site representation was lacking. As a result, sampling sites were private residences.
Similar to ESCAPE, we sampled during 14-day contiguous periods, in which we concurrently measured NO2 and NOx at five sites and NO2, NOx, PM2.5, and PM10 at another five sites [24,34]. We scheduled sampling at the participant’s convenience, and so the proportion of site types matched the whole study and that sampling setup types were geographically representative. Each site was sampled for three periods, with each period in a different season. In ESCAPE, sampling seasons were divided into four months of “winter,” four months of “summer,” and four months of “intermediate” season (i.e., spring and fall) [24,34]. However, Tucson has historically had five seasons: a wet and cool winter (December–February), a windy and temperate spring (March–April), a dry and hot pre-monsoon (May–June), a wet and hot monsoon (July–August), and a windy and temperate fall (September–November) [8].
To adapt, we defined summer as May through August, winter as November through February, and the remaining months as the intermediate season. Sampling in summer was distributed evenly among pre-monsoon and monsoon months. During each sampling period, another sampler setup measuring NO2, NOx, PM2.5, and PM10 was collocated at a U.S. Environmental Protection Agency (USEPA) National Core regulatory monitoring station in the center of the city, which is designed to accurately measure low concentrations of NO2, NOx, and PM2.5 using trace-level instrumentation [7]. In doing so, we were able to assess study measurement method accuracy against reference methods and calculate temporally corrected annual mean concentrations for each site for comparison to other ESCAPE-based studies (see Supplementary Materials for details). Sampling began in September 2015 and concluded one year later.
Following ESCAPE methods, we measured NO2 and NOx using Ogawa Badge Samplers (Ogawa & Company, USA, Inc., Pompano Beach, FL, USA) at a height of 2 m, which were analyzed at the University of Arizona (see Supplemental Materials for details). Particulate matter was sampled using two 10 L/min–rated Harvard impactors (Air Diagnostics and Engineering, Inc., Naples, ME, USA): one for PM2.5 and one for PM10. Unlike ESCAPE, which used a customized single pump to run both impactors, we used separate SKC Leland Legacy pumps (SKC Inc., Eighty Four, PA, USA), which were calibrated immediately before and after sampling at the site using a Mesa Labs DryCal Defender 520 (Mesa Laboratories, Inc., Butler, NJ, USA). The pumps were programmed to sample for a contiguous 15 min period out of every two hours [34].
Impactors were suspended 2 m above the ground and 1 m apart, such that the inlets faced downward to prevent rain from soaking the filters. To adapt to Tucson’s extreme summer heat (daytime high temperatures regularly reach 43 °C) and avoid pump overheating and damage (rated to 45 °C), we ventilated pump cases with ten 2.5-cm-diameter holes, equally spaced along the perimeter of the vertical sides of the case, along the top. A constantly running PC fan was placed inside the container with the pumps to ensure continuous air flow. This container was surrounded by a frame of PVC pipe, and reflective Mylar-foil-covered bubble-pocket insulation was suspended between the pipes to prevent direct sun exposure while also allowing air to circulate around the pump container (Figure 2). In addition, to avoid the need of a PC to program pumps from overheating during summer months, samplers could be set up or taken down only from 5:30 to 10:00 a.m.
To collect PM, we used Whatman 37-mm-wide Teflon filters with a 0.2 micron pore size and support ring (Whatman, Inc., Maidstone, UK), because the ESCAPE-suggested Andersen filters were unavailable. Weighing filters in Tucson, AZ, where the relative humidity is <20% most of the year, required a static de-ionizer to remove static charge prior to each time the filter was weighed. All filters were weighed three times using a Mettler Toledo XP2U Ultra Micro Balance accurate to 0.1 µg (Mettler Toledo, Columbus, OH, USA), and then, the weights were averaged, provided the smallest and largest weights were within 10% of each other. Before weighing, all filters were conditioned for ≥48 h in a climate-controlled chamber. In addition, two lab blank filters were weighed, and their mean difference was used to correct the filter mass post-sampling for changes due to humidity. Unlike the ESCAPE approach, which loaded filters into impactors in the field, we loaded impactors in the laboratory and sealed them in plastic bags until set up at the site to reduce incidental contamination of filters from wind-blow dust when loading the impactors and to speed up sampler installation at the site under extreme heat conditions. Nearly all samples met ESCAPE quality assurance and control metrics (see Supplementary Materials for details).

3. Results

Meteorological measures were taken at the Tucson International Airport, 16 km south of the collocated sampler at the reference monitor site (Table 1 and Figure 3). Winter months were characterized by small amounts of rainfall and cool temperatures, as well as increased wind speeds in November and late January. Sampling periods in spring were relatively dry, with a low dew point and high winds. There was no recorded precipitation during sampling periods in February, April, and May. In the pre-monsoon, temperatures and aridity peaked but were tempered by early monsoon rains that began in June. July was the wettest month, with substantial precipitation variability by the hour. Unexpectedly, August had little rain but was humid and had the second highest mean wind speeds. Fall sampling periods saw the end of the monsoon rains and high summer temperatures.
At the sampling site collocated with a reference monitor, the study monitor geometric mean (standard deviation) concentrations were NO2 = 3.91 (1.86) ppb, NOx = 7.24 (2.14) ppb, and PM2.5 = 4.79 (1.36) μg/m3, and all concentrations were below the relevant USEPA National Ambient Air Quality Standards (NO2 = 53 ppb; PM2.5 = 12 μg/m3). Measurements from our monitor and the reference setup were moderately to highly correlated (Spearman correlations), with NO2 r = 0.80 (p = 0.003) and NOx r = 0.88 (p = 0.002), while PM2.5 had an r = 0.41 (p = 0.13) (Figure 4). Correlations were nearly identical after including four additional measurement periods during November and December during which no other sites were sampled to avoid holidays as per ESCAPE protocols.
Compared to the reference monitor, the study setup tended to undermeasure NO2 and NOx (mean percent differences of −45% and −15%, respectively), with the absolute percent difference decreasing as TAPS concentrations increased (NO2 r = 0.66, p = 0.02 and NOx r = 0.58, p = 0.05). The pre-monsoon months of May and June had the greatest average percent differences: NO2 = −52% and NOx = −34%. For PM2.5, the TAPS monitor overestimated on average by 37% but, with absolute percent difference, negatively correlated with reference monitor concentrations (r = −0.69, p = 0.02). TAPS-monitored PM2.5 concentrations during July through September were an average of 92% greater than the regulatory monitor’s, compared to 18% greater the rest of the time.
Throughout the study area, we measured NO2 and NOx at 4 regional background, 23 urban background, and 12 street sites, and PM at 4 regional background, 11 urban background, and 4 street sites (Figure 3). Due to the study design, street sites tended to have limited secure space for the PM sampler setup, so these locations often accommodated only NO2 and NOx samplers. NO2 and NOx concentrations were significantly different between seasons, with the highest in winter and the lowest in summer (Wilcoxon signed-rank test, paired samples, Bonferroni correction, all p < 0.001), while PM concentrations were not (all p > 0.22) (Figure 5). During all seasons, street sites had higher NO2 and NOx concentrations, followed by urban and then regional backgrounds, though this was not significant (Kruskal–Wallis test). PM2.5 concentrations exhibited this same relationship in winter and intermediate seasons, but this was only significant in winter (p = 0.05). PM10 also showed this trend only in winter, and it was not significant. In summer, there were no trends by site type for PM measures.
NO2 and NOx measurements were highly correlated with one another (r = 0.94, p < 0.0001), with similar correlations among site types. With regard to the seasons, the correlation was the greatest in the intermediate season (r = 0.89, p < 0.0001), followed by that in winter (r = 0.84, p < 0.0001) and summer (r = 0.71, p < 0.0001). NO2/NOx ratios were not significantly different by site type or season (Figure 6).
When comparing annual means that were temporally corrected and those that were not, the former were almost always higher (Table 2). These two annual means were highly correlated for NO2 (r = 0.94), NOx (r = 0.95), and PM10 (r = 0.96) but less so for PM2.5 (r = 0.65). As shown previously, the NO2, NOx, and PM2.5 annual means at street sites were the highest, followed by those at urban and then regional backgrounds; however these differences were not significant. Meanwhile, PM10 had lower concentrations at urban background sites, but this was not significant.

4. Discussion

In this study, we adopted and adapted ESCAPE methods previously used in Europe, coastal cities Australian and South Africa, and densely populated Taipei, Taiwan, to the suburban, low-pollution Tucson metropolitan area in the desert southwestern U.S. While ESCAPE methods with updates for varying seasonality and extreme heat were effective in Tucson’s climatological extremes, measurements of fine particulate matter using an ESCAPE monitor were poorly correlated with those from the collocated reference monitor.
One reason may be that the reference monitor inlet was 4 m off the ground, which could result in the TAPS monitor collecting a greater number of larger particles nearer the ground, resulting in consistently higher PM2.5 concentrations. However, there was no relationship between the TAPS PM2.5/10 ratio and the reference monitor PM2.5 concentrations. Alternatively, sampling 15 min of every 2 h could have resulted in oversampling short-lived high PM2.5 sources (e.g., high-wind events or traffic emissions). This was performed as per ESCAPE protocols so that filters would not be overloaded over the 2-week sampling period.
In a past Tucson-based study, PM samplers were run continuously for a shorter sampling period using a Harvard impactor for lower flow rate (e.g., 4 L/min for 1 week) [35]. While PM2.5 was not reported, PM10 concentrations were higher than those in TAPS, indicating that this approach may not have overloaded filters. In our study, the absolute percent difference between the collocated monitors was negatively related to the reference monitor’s concentration, indicating the sampling method may lose accuracy in low-pollution airsheds, which are critical in demonstrating potential concentration thresholds in health effects.
In comparison to other desert locations, albeit coastal ones, Tucson’s mean PM2.5 and PM10 concentrations are about an order of magnitude lower (42.0 and 161 µg/m3 in Chile and 41.6 and 138 µg/m3 in Kuwait, respectively [14,15]). In Kuwait, using a Harvard Impactor (Air Diagnostics and Engineering, Inc., Naples, ME, USA) modified for high-dust settings for 24 h samples, authors found that PM2.5 measures had better agreement with the collocated regulatory monitor compared to PM10, likely due to the greater spatial distribution of that size fraction [15]. Further, their impactor setup consistently under measured PM2.5, with the greatest differences during dust storms, as opposed to ours which over measured compared to the regulatory monitor. Unfortunately, we do not have regulatory PM10 concentrations at the site to compare.
Interestingly, there were no dust storms or other abnormal or historic weather events, as defined by the NOAA Storm Events Database, during our study. However, the hottest months had the greatest PM2.5 percent difference, suggesting that the pumps may be less reliable when approaching their specified maximum operating temperature of 45 °C. While the maximum ambient temperature during sampling was 43 °C (May), temperatures inside cases could reach 45 °C if the ventilation fan failed. While we found no evidence of fan breakdown, it is possible that certain areas experienced power outages during periods of extreme heat. Unfortunately, there are no instrument diagnostics available for more information on whether temperature exceedances caused equipment issues. However, there was no linear trend between temperature and correlation. While pump diaphragms did need to be replaced periodically, this did not affect any of the reference monitor pumps.
It is notable that in Tucson and other southwestern U.S. cities, only 20–50% of PM variation can be explained by meteorological variables, and the remainder may come from unpredictable sources such as nearby vehicle traffic, forest fires, or transport from other locations [36]. While poor sampling agreement at the collocated site may be an isolated incident, it could systemically result in PM2.5 exposure misclassification later on during air pollution modeling for both measurements at samplers throughout the study area and the background setup at the monitoring station used for temporal correction of all other sites. We are unaware of any collocated PM sampling results for locations using ESCAPE methods to compare our findings. In contrast, NO2 and NOx measurements from the continuously sampling Ogawa badges were highly correlated with collocated regulatory monitors.
For continuous Ogawa sampling, our findings are similar not only to other ESCAPE locations in Europe [24] but also to two other studies using the sampler at reference sites in El Paso, TX, another city on the U.S.–Mexico border in the western U.S., where NO2 concentrations were also highly correlated (r = 0.91 and 0.94) [23,37]. Interestingly, while these El Paso–based studies had concentrations more than double those of TAPS, Perth, Australia, had low pollutant concentrations and correlations nearly identical to Tucson’s [27]. Together, this indicates that Ogawa methods are accurate against reference monitors, even when sampling at relatively low concentrations in arid conditions.
Across seasons and site types in TAPS, measurements and annual mean concentrations were often less than half of the relevant National Ambient Air Quality Standards (NAAQS). Concentrations were highest in winter, much like in El Paso, TX [23,37], likely owing to thermal inversions and stagnation and a swell in anthropogenic sources from Tucson’s seasonal winter visitors [7]. Interestingly, NO2/NOx ratios were not different by season or site type, suggesting that sources do not vary drastically over time. While street sites would be expected to have higher ratios owing to diesel vehicles, which emit a higher proportion of NO2 to NOx compared to gasoline vehicles [38], this was only seen in the winter and intermediate seasons. In summer, increased sunlight may be increasing the conversion rate of diesel-based NO2 to NO [36,39]. However, the ratios are comparable to those of ESCAPE, indicating similar pollutant sources (and/or sinks), despite Europe having a larger percentage of diesel vehicles [40,41].
Similar to the case in ESCAPE, street sites generally had the highest annual mean concentrations, while rural background sites had the lowest. Differences in our study were not significant, unlike in other ESCAPE study locations [24,26,34]. Tucson’s lack of large point sources may be one explanation, combined with differences in fleet characteristics and fuel types [40,41]. In addition, Europe and Taipei are generally more densely developed, resulting in more pollution sources in a smaller area, often with street canyons reducing dilution [42,43]. For example, the City of Tucson, the most densely developed portion of the study area, has 890 people/km2 and no building canyons, while Oslo, Norway, has 1532 people/km2. Compared to Oslo, Tucson’s concentrations are slightly higher, yet they are far lower than those for cities in southern Europe such as Barcelona, Spain (37.4 μg/m3), with more persons at a much higher density (1,620,343 people; 16,000 people/km2). Relatedly, stop-and-go driving, which results in greater PM2.5 emissions per mile, is more common in denser areas [44,45]. Interestingly, PM means and ranges in TAPS were nearly identical to those in the coastal metro area of Perth, Australia (1.97 M persons; 307 persons/km2). This may be explained by Perth’s lower population density compared to Europe and the City of Tucson, as well as vehicle fleet characteristics and wind-blown dust sources [28].
It is worth noting that only 4/19 (20%) of TAPS PM sites were street sites, compared to no less than 35% in ESCAPE, likely because only certain residents had the requisite security, electricity, space, and willingness to host a PM sampler setup. Similarly, just 12/39 (30%) of TAPS NO2 and NOx sites were street sites, versus about 50% in ESCAPE. As a result, our study may have under sampled areas that could meet PM site requirements, specifically in areas with no secure outdoor space, more density, and likely higher pollution concentrations. However, Perth, Australia, which had about 50% street sites, had pollutant concentrations nearly identical to those found in TAPS [27], suggesting more street sites would not have substantively expanded the distribution of Tucson’s concentrations.

5. Conclusions

In conclusion, our study illustrates how ESCAPE methods can be generally adopted and adapted to a desert metropolis with air pollution concentrations below regulatory levels with health impacts [46]. While most methods were robust to Tucson’s volatile precipitation patterns and extreme heat and aridity, study PM2.5 measurements were poorly correlated with those from a collocated reference monitor. This may have been because of differences in inlet heights; oversampling of high PM2.5 events because of non-continuous sampling; the pump operating beyond its temperature specifications; and/or other unknown variables.
While our study has notable strengths based on the uniqueness of Tucson’s volatile climate and low-pollution airshed and the use of well-documented methods coupled with collocated regulatory monitoring, it is not without weaknesses. Our pumps were unable to produce diagnostic output that might have confirmed our suspicions of mechanical issues resulting from high ambient temperatures, potentially leading to poor correlations between collocated study and regulatory monitors. Meteorology data obtained geographically closer to the regulatory monitor could have shed more light on factors that may have influenced poor correlation for PM2.5 with less potential for highly spatially variable weather patterns. PM10 concentrations from a regulatory monitor were also unavailable at the regulatory monitor site, blunting our ability to investigate the relationships between study and regulatory setups.
As climate change makes meteorological patterns more extreme, much like Tucson’s is today [10,11,12,13], it will be critical to reevaluate air pollution monitoring methods. Future studies should take advantage of collocating with regulatory monitor setups to help ensure that study samplers are field validated against “gold standard” methods, especially those meant to more accurately measure low concentrations. This is even more important for low-pollution airsheds, where there is less room for error. Without this assurance, sampling bias in different meteorological or pollution scenarios may compound in later air pollution and health effects models. Low air pollution concentrations with intra-city spatial precision are rare [47], yet they are critical in assessing the impacts of reductions in air pollution.

Supplementary Materials

The following supporting information on site selection and characterization, as well as sampling quality control and assurance can be downloaded at https://www.mdpi.com/article/10.3390/ijerph19063173/s1: Supplementary Materials.

Author Contributions

Conceptualization, N.L., N.L.-G., and P.B.; data curation, N.L.; formal analysis, N.L. and P.B.; funding acquisition, P.B.; investigation, N.L., R.A.C., M.K.O., and P.B.; methodology, N.L., N.L.-G., R.A.C., M.K.O., and P.B.; project administration, N.L. and N.L.-G.; resources, N.L. and P.B.; software, N.L.; supervision, N.L.; visualization, N.L.; writing—original draft, N.L. and N.L.-G.; writing—review and editing, N.L., N.L.-G., R.A.C., M.K.O., S.G., and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Institutes of Health (grants ES006694 and HL103970) and a 2015 Arizona Technology and Research Initiative Fund grant. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

The University of Arizona Human Subjects Research Review Board determined that the project did not constitute human subject research.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon reasonable request by contacting N.L. at [email protected].

Acknowledgments

We would like to thank the people who participated in our study and were gracious enough to donate their time and patience. We gratefully acknowledge Oscar Infante, Deidre Kruckenberg, and Rebecca Monroy for their help in developing and completing this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dockery, D.W.; Speizer, F.E.; Stram, D.O.; Ware, J.H.; Spengler, J.D.; Ferris, B.G. Effects of inhalable particles on respiratory health of children. Am. Rev. Respir. Dis. 1989, 139, 587–594. [Google Scholar] [CrossRef] [PubMed]
  2. Speizer, F.E. Asthma and Persistent Wheeze in the Harvard Six Cities Study. Chest 1990, 98, 191S–195S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ware, J.H.; Dockery, D.W.; Spiro, A.; Speizer, F.E.; Ferris, B.G. Passive smoking, gas cooking, and respiratory health of children living in six cities. Am. Rev. Respir. Dis. 1984, 129, 366–374. [Google Scholar] [PubMed]
  4. Jacquemin, B.; Siroux, V.; Sanchez, M.; Carsin, A.-E.; Schikowski, T.; Adam, M.; Bellisario, V.; Buschka, A.; Bono, R.; Brunekreef, B.; et al. Ambient Air Pollution and Adult Asthma Incidence in Six European Cohorts (ESCAPE). Environ. Health Perspect. 2015, 123, 613–621. [Google Scholar] [CrossRef]
  5. Gulliver, J.; Elliott, P.; Henderson, J.; Hansell, A.L.; Vienneau, D.; Cai, Y.; McCrea, A.; Garwood, K.; Boyd, A.; Neal, L.; et al. Local- and regional-scale air pollution modelling (PM10) and exposure assessment for pregnancy trimesters, infancy, and childhood to age 15 years: Avon Longitudinal Study of Parents And Children (ALSPAC). Environ. Int. 2018, 113, 10–19. [Google Scholar] [CrossRef] [PubMed]
  6. US Census Bureau. 2013–2017 American Community Survey 5-Year Estimates. Available online: https://factfinder.census.gov/faces/nav/jsf/pages/community_facts.xhtml?src=bkmk (accessed on 5 April 2019).
  7. Pima County Department of Environmental Quality. Air Quality Summary Report for Pima County, Arizona; Pima County: Tucson, AZ, USA, 2016. [Google Scholar]
  8. US Environmental Protection Agency. Technical Air Pollution Resources. Available online: http://www.epa.gov/ttn/airs/aqsdatamart (accessed on 20 September 2016).
  9. Tong, D.Q.; Wang, J.X.L.; Gill, T.E.; Lei, H.; Wang, B. Intensified dust storm activity and Valley fever infection in the southwestern United States. Geophys. Res. Lett. 2017, 44, 4304–4312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Christensen, J.H.; Kanikicharla, K.K.; Aldrian, E.; An, S.I.; Albuquerque Cavalcanti, I.F.; de Castro, M.; Dong, W.; Goswami, P.; Hall, A.; Kanyanga, J.K.; et al. Climate phenomena and their relevance for future regional climate change. In Climate Change 2013: The Physical Science Basis; Stocker, T.F., Qin, D.M., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  11. Donat, M.G.; Lowry, A.L.; Alexander, L.V.; O’Gorman, P.A.; Maher, N. More extreme precipitation in the worldâ €TM s dry and wet regions. Nat. Clim. Chang. 2016, 6, 508–513. [Google Scholar] [CrossRef]
  12. Prein, A.F.; Rasmussen, R.M.; Ikeda, K.; Liu, C.; Clark, M.P.; Holland, G.J. The future intensification of hourly precipitation extremes. Nat. Clim. Chang. 2017, 7, 48–52. [Google Scholar] [CrossRef]
  13. Masson-Delmotte, V.; Zhai, P.; Pirani, A.; Connors, S.L.; Péan, C.; Berger, S.; Caud, N.; Chen, Y.; Goldfarb, L.; Gomis, M.I.; et al. IPCC, 2021: Summary for Policymakers; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  14. Jorquera, H.; Barraza, F. Source apportionment of PM10 and PM2.5 in a desert region in northern Chile. Sci. Total Environ. 2013, 444, 327–335. [Google Scholar] [CrossRef]
  15. Alahmad, B.; Al-Hemoud, A.; Kang, C.M.; Almarri, F.; Kommula, V.; Wolfson, J.M.; Bernstein, A.S.; Garshick, E.; Schwartz, J.; Koutrakis, P. A two-year assessment of particulate air pollution and sources in Kuwait. Environ. Pollut. 2021, 282, 117016. [Google Scholar] [CrossRef]
  16. Hoek, G.; Beelen, R.; de Hoogh, K.; Vienneau, D.; Gulliver, J.; Fischer, P.; Briggs, D. A review of land-use regression models to assess spatial variation of outdoor air pollution. Atmos. Environ. 2008, 42, 7561–7578. [Google Scholar] [CrossRef]
  17. Jerrett, M.; Arain, A.; Kanaroglou, P.; Beckerman, B.; Potoglou, D.; Sahsuvaroglu, T.; Morrison, J.; Giovis, C. A review and evaluation of intraurban air pollution exposure models. J. Expo. Sci. Environ. Epidemiol. 2005, 15, 185–204. [Google Scholar] [CrossRef] [PubMed]
  18. De Hoogh, K.; Korek, M.; Vienneau, D.; Keuken, M.; Kukkonen, J.; Nieuwenhuijsen, M.J.; Badaloni, C.; Beelen, R.; Bolignano, A.; Cesaroni, G.; et al. Comparing land use regression and dispersion modelling to assess residential exposure to ambient air pollution for epidemiological studies. Environ. Int. 2014, 73, 382–392. [Google Scholar] [CrossRef]
  19. de Hoogh, K.; Chen, J.; Gulliver, J.; Hoffmann, B.; Hertel, O.; Ketzel, M.; Bauwelinck, M.; van Donkelaar, A.; Hvidtfeldt, U.A.; Katsouyanni, K.; et al. Spatial PM2.5, NO2, O3 and BC models for Western Europe–Evaluation of spatiotemporal stability. Environ. Int. 2018, 120, 81–92. [Google Scholar] [CrossRef] [PubMed]
  20. Gehring, U.; Cyrys, J.; Sedlmeir, G.; Brunekreef, B.; Bellander, T.; Fischer, P.; Bauer, C.P.; Reinhardt, D.; Wichmann, H.E.; Heinrich, J. Traffic-related air pollution and respiratory health during the first 2 years of life. Eur. Respir. J. 2002, 19, 690–698. [Google Scholar] [CrossRef] [Green Version]
  21. Amini, H.; Taghavi-shahri, S.M.; Henderson, S.B.; Nadda, K.; Nabizadeh, R.; Yunesian, M. Land use regression models to estimate the annual and seasonal spatial variability of sulfur dioxide and particulate matter in Tehran, Iran. Sci. Total Environ. 2014, 489, 343–353. [Google Scholar] [CrossRef]
  22. Meng, X.; Fu, Q.; Ma, Z.; Chen, L.; Zou, B.; Zhang, Y.; Xue, W.; Wang, J.; Wang, D.; Kan, H.; et al. Estimating ground-level PM10 in a Chinese city by combining satellite data, meteorological information and a land use regression model. Environ. Pollut. 2016, 208, 177–184. [Google Scholar] [CrossRef]
  23. Gonzales, M.; Myers, O.; Smith, L.; Olvera, H.A.; Mukerjee, S.; Li, W.; Pingitore, N.; Amaya, M.; Burchiel, S.; Berwick, M. Evaluation of land use regression models for NO2 in El Paso, Texas, USA. Sci. Total Environ. 2012, 432, 135–142. [Google Scholar] [CrossRef] [Green Version]
  24. Cyrys, J.; Eeftens, M.; Heinrich, J.; Ampe, C.; Armengaud, A.; Beelen, R.; Bellander, T.; Beregszaszi, T.; Birk, M.; Cesaroni, G.; et al. Variation of NO2 and NOx concentrations between and within 36 European study areas: Results from the ESCAPE study. Atmos. Environ. 2012, 62, 374–390. [Google Scholar] [CrossRef]
  25. Brunekreef, B. Study Manual for the European Study of Cohorts for Air Pollution Effects; Universiteit Utrecht: Utrecht, The Netherlands, 2008. [Google Scholar]
  26. Lee, J.H.; Wu, C.F.; Hoek, G.; de Hoogh, K.; Beelen, R.; Brunekreef, B.; Chan, C.C. Land use regression models for estimating individual NOx and NO2 exposures in a metropolis with a high density of traffic roads and population. Sci. Total Environ. 2014, 472, 1163–1171. [Google Scholar] [CrossRef]
  27. Dirgawati, M.; Barnes, R.; Wheeler, A.J.; Arnold, A.-L.L.; McCaul, K.A.; Stuart, A.L.; Blake, D.; Hinwood, A.; Yeap, B.B.; Heyworth, J.S. Development of Land Use Regression models for predicting exposure to NO2 and NOx in Metropolitan Perth, Western Australia. Environ. Model. Softw. 2015, 74, 258–267. [Google Scholar] [CrossRef]
  28. Dirgawati, M.; Heyworth, J.S.; Wheeler, A.J.; McCaul, K.A.; Blake, D.; Boeyen, J.; Cope, M.; Yeap, B.B.; Nieuwenhuijsen, M.; Brunekreef, B.; et al. Development of Land Use Regression models for particulate matter and associated components in a low air pollutant concentration airshed. Atmos. Environ. 2016, 144, 69–78. [Google Scholar] [CrossRef]
  29. Rahman, M.M.; Yeganeh, B.; Clifford, S.; Knibbs, L.D.; Morawska, L. Development of a land use regression model for daily NO2 and NOx concentrations in the Brisbane metropolitan area, Australia. Environ. Model. Softw. 2017, 95, 168–179. [Google Scholar] [CrossRef] [Green Version]
  30. Muttoo, S.; Ramsay, L.; Brunekreef, B.; Beelen, R.; Meliefste, K.; Naidoo, R.N. Land use regression modelling estimating nitrogen oxides exposure in industrial south Durban, South Africa. Sci. Total Environ. 2018, 610–611, 1439–1447. [Google Scholar] [CrossRef] [PubMed]
  31. Seinfeld, J.H.; Pandis, S.N. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 1st ed.; Wiley-Interscience: New York, NY, USA, 1998. [Google Scholar]
  32. Abhijith, K.V.; Kumar, P.; Gallagher, J.; Mcnabola, A.; Baldauf, R.; Pilla, F.; Broderick, B.; Di, S.; Pulvirenti, B. Air pollution abatement performances of green infrastructure in open road and built-up street canyon environments–A review. Atmos. Environ. 2017, 162, 71–86. [Google Scholar] [CrossRef]
  33. Gonzalez, P.; Wang, F.; Notaro, M.; Vimont, D.J.; Williams, J.W. Disproportionate magnitude of climate change in United States national parks. Environ. Res. Lett. 2018, 13, 104001. [Google Scholar] [CrossRef] [Green Version]
  34. Eeftens, M.; Tsai, M.Y.; Ampe, C.; Anwander, B.; Beelen, R.; Bellander, T.; Cesaroni, G.; Cirach, M.; Cyrys, J.; de Hoogh, K.; et al. Spatial variation of PM2.5, PM10, PM2.5 absorbance and PMcoarse concentrations between and within 20 European study areas and the relationship with NO2–Results of the ESCAPE project. Atmos. Environ. 2012, 62, 303–317. [Google Scholar] [CrossRef] [Green Version]
  35. Quackenboss, J.; Lebowitz, M.; Crutchfield, C. Indoor-outdoor relationships for particulate matter: Exposure classifications and health effects. Environ. Int. 1989, 15, 353–360. [Google Scholar] [CrossRef]
  36. Wise, E.K.; Comrie, A.C. Meteorologically adjusted urban air quality trends in the Southwestern United States. Atmos. Environ. 2005, 39, 2969–2980. [Google Scholar] [CrossRef]
  37. Sather, M.E.; Slonecker, E.T.; Mathew, J.; Daughtrey, H.; Williams, D.D. Evaluation of ogawa passive sampling devices as an alternative measurement method for the nitrogen dioxide annual standard in El Paso, Texas. Environ. Monit. Assess. 2007, 124, 211–221. [Google Scholar] [CrossRef]
  38. Keuken, M.; Roemer, M.; van den Elshout, S. Trend analysis of urban NO2 concentrations and the importance of direct NO2 emissions versus ozone/NOx equilibrium. Atmos. Environ. 2009, 43, 4780–4783. [Google Scholar] [CrossRef]
  39. Diem, J.E.; Comrie, A.C. Air Quality, Climate, and Policy: A Case Study of Ozone Pollution in Tucson, Arizona. Prof. Geogr. 2001, 53, 469–491. [Google Scholar] [CrossRef]
  40. US Energy Information Administration (EIA). Annual Energy Outlook 2015; EIA: Washington, DC, USA, 2015.
  41. ACEA. Fuel Types of New Passenger Cars in the EU. Available online: https://www.acea.be/statistics/tag/category/share-of-diesel-in-new-passenger-cars (accessed on 2 July 2019).
  42. Boogaard, H.; Kos, G.P.A.; Weijers, E.P.; Janssen, N.A.H.; Fischer, P.H.; van der Zee, S.C.; de Hartog, J.J.; Hoek, G. Contrast in air pollution components between major streets and background locations: Particulate matter mass, black carbon, elemental composition, nitrogen oxide and ultrafine particle number. Atmos. Environ. 2011, 45, 650–658. [Google Scholar] [CrossRef]
  43. Boarnet, M.G.; Houston, D.; Edwards, R.; Princevac, M.; Ferguson, G.; Pan, H.; Bartolome, C. Fine particulate concentrations on sidewalks in five Southern California cities. Atmos. Environ. 2011, 45, 4025–4033. [Google Scholar] [CrossRef]
  44. Kim, W.G.; Kim, C.K.; Lee, J.T.; Kim, J.S.; Yun, C.W.; Yook, S.J. Fine particle emission characteristics of a light-duty diesel vehicle according to vehicle acceleration and road grade. Transp. Res. Part D Transp. Environ. 2017, 53, 428–439. [Google Scholar] [CrossRef]
  45. Imhof, D.; Weingartner, E.; Ordóñez, C.; Gehrig, R.; Hill, M.; Buchmann, B.; Baltensperger, U. Real-world emission factors of fine and ultrafine aerosol particles for different traffic situations in Switzerland. Environ. Sci. Technol. 2005, 39, 8341–8350. [Google Scholar] [CrossRef]
  46. Beamer, P.I.; Furlong, M.; Lothrop, N.; Guerra, S.; Billheimer, D.; Stern, D.A.; Zhai, J.; Halonen, M.; Wright, A.L.; Martinez, F.D. CC16 Levels into Adult Life are Associated with Nitrogen Dioxide Exposure at Birth. Am. J. Respir. Crit. Care Med. 2019, 200, 600–607. [Google Scholar] [CrossRef]
  47. Papadogeorgou, G.; Kioumourtzoglou, M.A.; Braun, D.; Zanobetti, A. Low Levels of Air Pollution and Health: Effect Estimates, Methodological Challenges, and Future Directions. Curr. Environ. Health Rep. 2019, 6, 105–115. [Google Scholar] [CrossRef]
Figure 1. The Tucson metropolitan area (starred) in the bolded state of Arizona (AZ).
Figure 1. The Tucson metropolitan area (starred) in the bolded state of Arizona (AZ).
Ijerph 19 03173 g001
Figure 2. A typical sampler setup for NO2, NOx, PM2.5, and PM10. The Ogawa Badge Sampler (a) is attached near the top of the central post, while Harvard impactors hang, inlet down, from each end of the horizontal plank (b). Leland Legacy pumps are in a pump case concealed by UV-reflective foil insulation.
Figure 2. A typical sampler setup for NO2, NOx, PM2.5, and PM10. The Ogawa Badge Sampler (a) is attached near the top of the central post, while Harvard impactors hang, inlet down, from each end of the horizontal plank (b). Leland Legacy pumps are in a pump case concealed by UV-reflective foil insulation.
Ijerph 19 03173 g002
Figure 3. Monitoring sites in the Tucson Air Pollution Study area, including the regulatory monitoring site and meteorology station.
Figure 3. Monitoring sites in the Tucson Air Pollution Study area, including the regulatory monitoring site and meteorology station.
Ijerph 19 03173 g003
Figure 4. Comparison of measured concentrations of (a) NO2, (b) NOx, and (c) PM2.5 by month for the collocated Tucson Air Pollution Study (TAPS) and reference (Ref.) monitors. Points are labeled with month. The solid line is the identity line; the dashed line, the linear regression line; and the gray shaded area, the 95% confidence interval. Months 1a and 1b are the first and second two-week periods in January, respectively.
Figure 4. Comparison of measured concentrations of (a) NO2, (b) NOx, and (c) PM2.5 by month for the collocated Tucson Air Pollution Study (TAPS) and reference (Ref.) monitors. Points are labeled with month. The solid line is the identity line; the dashed line, the linear regression line; and the gray shaded area, the 95% confidence interval. Months 1a and 1b are the first and second two-week periods in January, respectively.
Ijerph 19 03173 g004
Figure 5. Pollutant measures by season and site type for (a) NO2, (b) NOx, (c) PM2.5, and (d) PM10. In boxplots, the bold horizontal line represents the median, lower and upper ends of the box represent the first and third quartiles, the upper whisker extends to 1.5 × IQR above the box,. and the lower whisker extends to 1.5 × IQR below the box. Reg.: regional background; Urb.: urban background.
Figure 5. Pollutant measures by season and site type for (a) NO2, (b) NOx, (c) PM2.5, and (d) PM10. In boxplots, the bold horizontal line represents the median, lower and upper ends of the box represent the first and third quartiles, the upper whisker extends to 1.5 × IQR above the box,. and the lower whisker extends to 1.5 × IQR below the box. Reg.: regional background; Urb.: urban background.
Ijerph 19 03173 g005
Figure 6. NO2/NOx ratios by season and site type. In boxplots: the bold horizontal line is the median, lower and upper ends of the box are the first and third quartiles, the upper whisker extends to 1.5 × IQR above the box, and the lower whisker extends to 1.5 × IQR below the box. Reg.: regional background; Urb.: urban background.
Figure 6. NO2/NOx ratios by season and site type. In boxplots: the bold horizontal line is the median, lower and upper ends of the box are the first and third quartiles, the upper whisker extends to 1.5 × IQR above the box, and the lower whisker extends to 1.5 × IQR below the box. Reg.: regional background; Urb.: urban background.
Ijerph 19 03173 g006
Table 1. Hourly meteorological measures by sampling period month.
Table 1. Hourly meteorological measures by sampling period month.
Precip. (mm)Sea-Level Pressure (hPa)Dew Point (°C)Temp. (°C)Wind Speed (m/s)Visibility (km)
Sampling MonthTotalMean (SD)Mean (SD)Mean (SD)Mean (SD)Mean (SD)Mean (SD)
113.300.009 (0.09)1016 (3.81)−1.65 (4.60)14.6 (6.40)3.42 (2.09)16.1 (0.41)
1a1.020.002 (0.04)1020 (4.00)−2.22 (3.60)9.95 (5.81)2.60 (1.64)16.1 (0.28)
1b4.570.01 (0.16)1020 (5.76)−7.51 (3.72)13.5 (8.00)3.40 (2.28)16.1 (0.30)
20-1015 (3.30)−6.70 (2.64)18.5 (6.55)3.20 (1.74)16.1 (0)
36.450.02 (0.22)1013 (4.17)−5.18 (6.05)19.1 (6.50)3.23 (1.95)16.1 (0.35)
40-1009 (4.21)−5.05 (4.17)21.3 (5.83)3.89 (2.41)16.1 (0)
50-1008 (3.65)−4.60 (3.42)28.0 (7.06)3.25 (1.90)16.1 (0)
656.10.12 (1.00)1009 (2.96)14.0 (2.89)31.9 (5.00)3.24 (1.51)16.0 (0.80)
780.60.20 (1.25)1009 (3.42)19.8 (0.52)29.8 (4.76)3.10 (1.78)15.9 (1.04)
87.970.02 (0.16)1012 (3.02)16.3 (0.81)25.6 (5.00)3.49 (2.19)16.1 (0.27)
938.40.11 (0.92)1009 (2.53)12.0 (4.30)27.8 (4.95)3.26 (2.01)15.9 (1.44)
1040.40.11 (0.74)1011 (3.26)9.68 (3.36)20.4 (5.01)3.26 (2.18)15.9 (1.00)
SD: standard deviation; months 1a and 1b are the first and second two-week periods in January, respectively.
Table 2. Summary of measured and temporally corrected annual mean pollutant concentrations (ppb) by study site type.
Table 2. Summary of measured and temporally corrected annual mean pollutant concentrations (ppb) by study site type.
NO2NOx PM2.5PM10
Site TypenGM (GSD)RangeGM (GSD)RangenGM (GSD)RangeGM (GSD)Range
MeasuredReg.43.42 (1.92)1.35–6.286.32 (1.60)3.84–10.343.38 (1.41)2.06–4.4518.3 (1.24)13.6–21.9
Urb.234.57 (1.53)1.98–11.38.78 (1.48)3.89–21.9114.34 (1.24)3.23–6.7317.1 (1.61)5.35–30.1
Street126.22 (1.58)2.33–10.711.0 (1.56)4.79–23.044.73 (1.10)4.20–5.1818.5 (1.57)9.65–26.6
CorrectedReg.44.09 (2.05)1.61–7.518.00 (1.50)4.64–11.843.27 (1.21)2.67–4.1719.8 (1.19)16.2–24.0
Urb.235.11 (1.79)1.36–13.010.7 (1.79)1.56–30.2113.92 (1.28)2.52–5.7218.8 (1.57)7.25–34.8
Street127.39 (1.72)2.80–15.814.3 (1.65)6.38–33.644.83 (1.19)4.18–6.1921.0 (1.55)11.1–29.1
GM: geometric mean; GSD: geometric standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lothrop, N.; Lopez-Galvez, N.; Canales, R.A.; O’Rourke, M.K.; Guerra, S.; Beamer, P. Sampling Low Air Pollution Concentrations at a Neighborhood Scale in a Desert U.S. Metropolis with Volatile Weather Patterns. Int. J. Environ. Res. Public Health 2022, 19, 3173. https://doi.org/10.3390/ijerph19063173

AMA Style

Lothrop N, Lopez-Galvez N, Canales RA, O’Rourke MK, Guerra S, Beamer P. Sampling Low Air Pollution Concentrations at a Neighborhood Scale in a Desert U.S. Metropolis with Volatile Weather Patterns. International Journal of Environmental Research and Public Health. 2022; 19(6):3173. https://doi.org/10.3390/ijerph19063173

Chicago/Turabian Style

Lothrop, Nathan, Nicolas Lopez-Galvez, Robert A. Canales, Mary Kay O’Rourke, Stefano Guerra, and Paloma Beamer. 2022. "Sampling Low Air Pollution Concentrations at a Neighborhood Scale in a Desert U.S. Metropolis with Volatile Weather Patterns" International Journal of Environmental Research and Public Health 19, no. 6: 3173. https://doi.org/10.3390/ijerph19063173

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop