Lagged Association of Ambient Outdoor Air Pollutants with Asthma-Related Emergency Department Visits within the Pittsburgh Region

Abstract

Asthma affects millions of people globally and is especially concerning in populations living with poor air quality. This study examines the association of ambient outdoor air pollutants on asthma-related emergency department (ED) visits in children and adults throughout the Pittsburgh region. A time-stratified case-crossover design is used to analyze the lagged effects of fine particulate matter (PM_2.5) and gaseous pollutants, e.g., ozone (O₃), sulfur dioxide (SO₂), nitrogen dioxide (NO₂), and carbon monoxide (CO) on asthma-related ED visits (n = 6682). Single-, double-, and multi-pollutant models are adjusted for temperature and analyzed using conditional logistic regression. In children, all models show an association between O₃ and increased ED visits at lag day 1 (OR: 1.12, 95% CI, 1.03–1.22, p < 0.05) for the double-pollutant model (OR: 1.10, 95% CI: 1.01-1.20, p < 0.01). In adults, the single-pollutant model shows associations between CO and increased ED visits at lag day 5 (OR: 1.13, 95% CI, 1.00–1.28, p < 0.05) and average lag days 0–5 (OR: 1.22, 95% CI: 1.00–1.49, p < 0.05), and for NO₂ at lag day 5 (OR: 1.04, 95% CI: 1.00–1.07, p < 0.05). These results show an association between air pollution and asthma morbidity in the Pittsburgh region and underscore the need for mitigation efforts to improve public health outcomes.

1. Introduction

Asthma is a common chronic lung disease with varying phenotypes, some of which may be worsened by environmental factors such as air pollution. Air pollution is a complex mixture of gases and particulate matter (PM) that varies in concentration across the United States. The variability in concentration and composition is due to differing weather patterns and pollution sources [1]. Air pollution also varies from day to day and season to season within a region [2]. Numerous studies have associated air pollutants with adverse health outcomes, such as asthma exacerbations [3,4,5,6,7,8,9]. Short-term exposure to PM_2.5 (PM smaller than 2.5 μm in diameter), O₃, and other gaseous pollutants such as SO₂, NO₂, and CO (common products associated with the burning of fossil fuels and industrial emissions) have been shown to trigger asthma exacerbations and result in increased emergency department (ED) admissions in both children and adults [10,11] O₃ is unique from the other pollutants since it is not a primary pollutant. Rather, it is a secondary pollutant formed by photochemical reactions between sunlight and pollutant precursors, such as volatile organic compounds and nitrogen oxides. In some instances, the literature on health effects from air pollution is inconsistent and thus warrants further investigation [12]. The U.S. Environmental Protection Agency (EPA) sets the National Ambient Air Quality Standards for criteria pollutants based in part on controlled human exposure studies assessing airway hyper-responsiveness (AHR) [13]. These controlled exposure studies help understand how the lung responds to air pollutants.

Pittsburgh, located in southwestern Pennsylvania, has a unique topography because its metropolis is built within the hilly Appalachian Mountains [14]. According to the American Lung Association’s 2020 annual air quality report, Pittsburgh ranks poorly in the country for year-round particle pollution and O₃ [15], thus making this region a good location to study pollution-related health effects. A previously-published study in Pittsburgh examined the potential influence of gender and ethnicity on asthma-related ED visits between 2002 and 2005 and found O₃ and PM_2.5 to be significantly associated with an increased likelihood of ED visits for both Black and White Americans [14]. This prior study used a two-pollutant model and did not assess differential effects based on age. Other studies have found similar associations between ED visits and O₃ [16,17,18]; however, few of these studies included five pollutant models. The primary objective of this study was to examine the association between ambient air pollution and the risk of asthma-related ED visits in both children and adults using a five-pollutant model.

2. Materials and Methods

A time-stratified case-crossover study design with conditional logistic regression was used to examine the short-term associations of fine particulate matter (PM_2.5) and gaseous pollutants with asthma ED visits. The case-crossover design is commonly chosen to investigate transient effects on the risk of acute health events [19]. Specifically, this method is frequently used in epidemiology studies to examine the short-term effects of air pollution on respiratory disease events [20]. This design uses defined cases only and compares an individual’s exposure experience just prior to the event with exposure at other times (the referent periods). An advantage of the case-crossover design is that, since each case serves as its own control, confounding by risk factors that are invariant (such as age, sex, and race) or slowly-changing (such as tobacco smoke exposure and socio-economic status) are controlled by the design.

Data on ED visits from Allegheny Health Network (AHN), a six-hospital regional health system, with a primary discharge diagnosis of asthma from 1 July 2008 to 30 June 2013 were included in our analysis. All have hospital-based emergency departments staffed by emergency physicians with up-to-date diagnostic and treatment capabilities. The visits to each of these EDs ranged from 25,000 to 55,000 per year. Data for each visit included: name of hospital; pseudo ID; dates of admission and discharge; primary discharge diagnosis and up to three secondary diagnoses; zip code of residence; age, gender, and race; and disposition from ED (whether or not an individual was subsequently admitted as an in-patient or discharged to home). Specifically included in these analyses were visits with a primary discharge diagnosis of asthma, defined as ICD-9 Codes 493.XX. Analyses were limited to Allegheny County residents, ages five years and older. 2.1. Environmental Exposure Measures For this study, ambient pollutant data were downloaded from the Environmental Protection Agency (EPA) website and included concentrations and air quality index (AQI) values. A reference monitor in urban Pittsburgh (air quality system (AQS) ID: 420030038) measured PM_2.5, O₃, NO₂, and CO. SO₂ measurements were obtained from a separate reference monitor in Avalon, PA (AQS ID:420030002). This monitor was used as an indicator of regional SO₂ air pollution, since SO₂ was not measured at the Pittsburgh location. These monitoring sites are part of the EPA AQS that is used to monitor compliance with the Clean Air Act.

These reference monitors were selected because they measured the pollutants of interest, provided temporal data, and were centrally located between the six hospitals where asthma ED visit data were collected (Figure 1). To address missing values from the monitors due to technical errors in their functioning, we imputed missing values using the mean of the nearest valid values (one before and one after) to fill gaps in the data. The validation of this imputation has been deemed an acceptable process and has been scrutinized in many studies [15,16,21,22,23,24,25].

Meteorological variables, including the daily minimum and maximum air temperatures, were obtained from the Pittsburgh International Airport for dates between July 2008 and June 2013. These data were downloaded from the National Oceanic Atmospheric Association (NOAA) Climate Data online website. Average daily temperature was defined as the average of the minimum and maximum air temperatures for that day. Similar to the pollutant monitoring, missing values were imputed from the nearest days’ values. 2.2. Statistical Analysis For our analyses, the level of ambient air pollution at the time of day zero (day of admission) as well as lagged days just before an individual visited the ED were compared with levels at referent times. This approach is commonly used to evaluate the acute effects of transient air pollution exposures by comparing outcome risks in the same individual at different times. We used a time-stratified approach to select referent periods in 7-day intervals to minimize confounding factors such as day of the week, which could be associated with adverse pollution events. With this approach, the overall time period is divided into strata and exposure in the time period just prior to the event, and exposures in multiple reference periods are compared within that stratum of time. For this study, we used a 28-day strata and referent periods of 7, 14, and 21 days (either before or after the ED event) within each stratum. To assure independence of events within strata, a washout period of 7 days was used to remove recurrent events (ED visits) for an individual.

We evaluated the effect of exposure to pollutants on the day of the ED visit and prior days; lag 0 represents exposure on the same day as the ED visit. We also examined the effects of cumulative days of exposure to calculate average exposures over several days. The average of lag days 0–5 was calculated as the mean of lag days 0, 1, 2, 3, 4, and 5, respectively. Pollutant effects were examined with models containing just one single day or an average of lag days 0–5. We controlled for temperature effects in all models by including average temperature at the same lag(s) as the pollutant(s). Because case-crossover is so tightly controlled by the referent definition of “within 28 days”, cubic splines was not used. In addition to the single-pollutant models for PM_2.5, O₃, NO₂, SO₂, and CO, we conducted two-pollutant models for PM_2.5 and O₃ adjusted for temperature. Lastly, we conducted multi-pollutant models with all five pollutants. The results are presented as odds ratios (ORs) and the associated 95% confidence intervals for every 10 μg/m³ increase in PM_2.5; 10 ppb for O₃, NO₂, and SO₂; and 1 ppm for CO. Only significant findings are included within the text; however, the full data set may be found within Tables S1–S14.

Data analyses were conducted using the case-crossover tool (C-CAT) developed by Apex Epidemiology Research in collaboration with the New York State Department of Health for use with SAS (Abraham JH and Bateson TF, 2016). C-CAT is public domain software that provides an easy-to-use interface for SAS code to implement the time-stratified case-crossover analysis. Separate analyses were conducted for children (ages 5–17 years) and for adults (ages 18 years and older). 3. Results

The study population distribution by age and sex is shown in Figure 2. The majority of ED cases within our sample population were adults (n = 6682). Specifically, 87% were adults (n = 5842) and 13% were children (n = 840). Overall, there were higher rates of ED visits among women compared to men.

3.1. Air Pollution Data

The time scale for each pollutant recorded varies and is dependent upon EPA reporting requirements. The odds of an asthma ED visit related to an increase in exposure per increase of 10 (unit dependent on pollutant as referenced in the methods section) for a single day was significant at different lags for different pollutants. The strongest associations between pollutant levels and ED visits were observed for O₃ within single-, double-, and five-pollutant models for children, in addition to NO₂ and CO within single-pollutant models for adults. Adjusted odds ratios for ED visits for asthma, according to levels of ambient air pollution for each pollutant, may be found in Tables S1–S14. The average levels for the five year period were: PM_2.5: 11.41 µg/m³, SD ± 5.95; O₃: 38.5 ppb SD ± 15.9; SO₂: 9.25 ppb, SD ± 10.89; NO₂: 22.7 ppb, SD ± 9.4; and CO: 0.51 ppm, S.D. ± 0.27; and average temperature was 11.7 degrees Celsius SD ± 18.1.

3.2. Children Ages 5–17 Years

Statistically significant effects of O₃ were noted for lag day 1 in the single- and two-pollutant models (OR: 1.12, 95% CI: 1.02–1.22, p < 0.01) in the single pollutant model; and two-pollutant model (OR: 1.10, 95% CI: 1.01–1.20, p < 0.01) (Figure 3). This same effect was evident in the multi-pollutant model adjusting for PM_2.5, temperature, and the other pollutants (Supplementary Section S1, Table S14). No significant positive associations were observed between SO₂, NO₂, and CO and ED visits. See Supplementary Section S1, Table S4 for complete data. Of note, PM_2.5 in the two-pollutant model had a marginally significant protective effect, which requires further study.

3.3. Adults Aged 18 Years and Older

Statistically significant positive effects of CO were noted in adult ED visits for asthma on lag day 5 (OR: 1.13, 95% CI: 1.00–1.28, p < 0.01) and average lag days 0–5 (OR: 1.22, 95% CI: 1.00–1.49, p < 0.01) in the single-pollutant model. Similarly, a statistically-significant positive effect of NO₂ was seen in adult ED visits for asthma for lag day 5 (OR: 1.04, 95% CI: 1.00–1.07, p < 0.01) in the single-pollutant model (Figure 4). No significant associations were observed for SO₂ and O₃. See Supplementary Section S1, Tables S12 and S13, for complete data on CO and NO₂, respectively. In the five-pollutant model, no statistically significant positive effects were seen for any of the pollutants at any of the lags examined (Supplementary Section S1, Table S14).

3.4. Summary In children aged 5–17 years, there were statistically significant increases in asthma-related ED visits for O₃ on lag day 1 (one day before the visit) in all models. Typically, this was found because of the delay in seeking care. In adults 18 years and older, there was no statistically significant effect of ozone on asthma-related ED visits in any of the models. Instead, we found statistically significant effects of CO on lag day 5 and average lag days 0–5 and NO₂ on lag day 5 in the single-pollutant model only. 4. Discussion

Our analysis showed significant associations between daily pollution levels and asthma-related ED visits for both children and adults during the five-year study period. Higher ambient levels of O₃ were associated with increased ED visits in children, though not in adults. Increased ambient levels of NO₂ and CO were associated with increased odds of an ED visit in adults. Additionally, our study demonstrates differences in asthma-related ED visits based upon sex. This finding is in line with a prior report that post-pubescent women have been shown to have poorer asthma outcomes compared to their male counterparts [26].

A strength of our study was the extended lagged analyses through day 5. Some previous studies examined lag days 1–3 and found no association between pollution levels and ED visits. Significant associations observed at extended lag periods in the current study (e.g., lag day 5) suggest the potential for underlying biological mechanisms as potential contributors to this delayed response. Another possibility is the potential delay in asthmatics seeking ED treatment.

It was unexpected that our study did not show an association between PM_2.5 and ED visits in the single-pollutant model. The association of asthma exacerbations with elevations in short-term exposures to PM_2.5 are well established [27,28,29]. Studies suggest PM_2.5 may activate pathways associated with oxidative stress, which can increase airway hyper-responsiveness [30]. An interesting feature of PM_2.5 is its ability to serve as a transport vessel for airborne gaseous particles to travel deep into the bronchial airway. Therefore, it is important to consider how PM_2.5 concentrations change when combined with gaseous pollutants, as we did in our analysis.

Although the sample size used for the analysis in children was smaller than for adults, there were statistically significant increases in asthma-related ED visits for O₃ noted on lag day 1 (one day before the visit). Typically, this lag occurs because of the delay in seeking care. This finding was noted to overcome the collinearity among pollutants after adjustment for both temperature and PM_2.5. This finding is consistent with other studies showing an association between O₃ and asthma-related ED visits in children. Another study showed that after adjusting for seasonal variation, high levels of O₃ and SO₂ were associated with asthma exacerbations in children [31]. Another study demonstrated a positive relationship between O₃ levels and asthma-related ED visits, with associations being the strongest during the warm season [16,32,33]. In contrast to our findings in children, we did not find a positive association between O₃ and asthma-related ED visits in adults. Consistent with our findings, an Australian study showed increased air pollution affected ED visits for children but not adults [34]. There are several explanations for these discrepant findings in children and adults. Toxicological studies have shown that children tend to breathe in more air through their mouths, as opposed to adults, who breathe primarily through their nasal passages, which help to filter the air before reaching the lungs [35]. A systematic review of 27 epidemiological studies concluded that children may be at higher risk from O₃ because of their immature immune systems, increased durations of time spent outside, and increased air exchange relative to body mass, as compared to adults; therefore, higher exposure levels may be why we see effects from O₃ in children when compared with adults [8].

Lag day 5 was implicated in the adult analyses of NO₂ and CO. This lag may be related to the potential occupational exposures experienced by adults. Additional studies are needed to better understand the significance of this extended lag effect. Some studies suggest delayed effects might exist between various pollutants and asthma outcomes related to delayed physiologic responses [8]. Our finding of an association between NO₂ exposure, which is a known marker of traffic-related pollution that has been associated with lung inflammation, and asthma-related ED visits in adults is consistent with several previous studies. NO₂ was not implicated in asthma-related ED visits in children in our study.

Lastly, the primary effects of high outdoor CO exposures are related to hypoxia, which results in confusion, headache, and nausea [36]. CO might be a marker for other noxious combustion products, such as the burning of wood, coal, gas, and tobacco [36]. CO pollution has been associated with decreased lung function in adults with asthma [37,38,39,40]. CO was not linked with asthma ED visits in children within our study; however, a prior study reported increased odds of school-based health clinic visits related to high ambient levels of CO [41].

5. Study Limitations

As with all epidemiological studies, limitations exist for this study. For example, environmental factors such as elevation could not be completely controlled for and may have influenced the results. O₃ tends to be increased at higher elevations [42]. Furthermore, we used two EPA-grade reference monitors from different locations within Allegheny County to perform the analyses. Because of this small air monitor network, we relied on regional pollution data rather than hyper-local conditions. Therefore, we were unable to capture local-scale spatial variations that may have occurred for various pollutants. The more distant monitor that provided SO₂ concentrations may not have been representative of the entire cohorts’ ambient exposure, since the participating hospitals were clustered near the Pittsburgh reference monitor. Furthermore, due to the limited nature of our dataset, confounders such as socioeconomic status, type of health insurance, severity of asthma, and use of controller therapy were not examined. Lastly, ED visits for asthma are only one piece of a larger picture that describes asthma burden. Future studies should consider multiple endpoints (i.e., ED visits, school-based clinic visits, and outpatient physician office visits) to obtain a more comprehensive picture of how asthma control is exacerbated by air pollutants.

6. Conclusions There is an association between O₃ exposure in children and NO₂ and CO exposure in adults and asthma-related ED visits within the greater Pittsburgh area. Public health intervention(s) aimed at mitigating the effects of air pollutants targeted to the entire population may have significant benefits for children and adults with asthma, as well as the public as a whole.

Figure 1. Map of Allegheny County, PA, depicting locations of regulatory air monitors and hospitals where asthma-related Emergency Department (ED) visits were recorded between 2008 and 2013.

Figure 2. Cases of asthma-related ED visits between 2008 and 2013 separated by sex and age.

Figure 3.(A,B) Odds ratios of an asthma-related ED visit and significant lag days for children ages 5-17 using a two-pollutant model (single day lags and average lag, adjusted for apparent temperature (at same lag), which includes PM2.5 and O3, adjusted for temperature.) Referent time period is 28 days. * p-value < 0.05.

Figure 4.(A,B) Odds ratios of an asthma-related ED visit and significant lag days for adults aged 18 and older using a single-pollutant model for CO and NO2, adjusted for temperature. Referent time-period is 28 days. * p-value < 0.05.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/17/22/8619/s1, Table S1: Children/Teens Age 5-17: Asthma ED Visits for PM2.5 (per 10 µg/m3). Single Pollutant, Single day lags and average lag, Table S2: Children/Teens Age 5-17: Asthma ED Visits for Ozone (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S3: Children/Teens Age 5-17: Asthma ED Visits for PM2.5 (per 10 µg/m3) and Ozone (per 10 ppb). Two Pollutants, Single day lags and average lag, Table S4: Children/Teens Age 5-17: Asthma ED Visits for Carbon Monoxide (CO) (per 10 ppm). Single Pollutant, Single day lags and average lag, Table S5: Children/Teens Age 5-17: Asthma ED Visits for NO2 (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S6: Children/Teens Age 5-17: Asthma ED Visits for SO2 (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S7: Children/Teens Age 5-17: Asthma ED Visits for PM2.5 (per 10 µg/m3), Ozone (per 10 ppb), CO (ppm), NO2 (per 10 ppb), SO2 (per 10 ppb). Multiple Pollutant models, Single day lags and average lag, Table S8: Adults 18+: Asthma ED Visits for PM2.5 (per 10 µg/m3). Single Pollutant, Single day lags and average lag, Table S9: Adults 18+: Asthma ED Visits for Ozone (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S10: Adults 18+: Asthma ED Visits for PM2.5 (per 10 µg/m3) and Ozone (per 10 ppb). Two Pollutants, Single day lags and average lag, Table S11: Adults 18+: Asthma ED Visits for Carbon Monoxide (per 10 ppm). Single Pollutant, Single day lags and average lag, Table S12: Adults 18+: Asthma ED Visits. NO2 (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S13: Adults 18+: Asthma ED Visits for SO2 (per 10 ppb). Single Pollutant, Single day lags and average lag, Table S14: Adults 18+: Asthma ED Visits for PM2.5 (per 10 µg/m3), Ozone (per 10 ppb), CO (ppm), NO2 (per 10 ppb), SO2 (per 10 ppb). Multiple Pollutant models, Single day lags and average lag.

Author Contributions

Project Conceptualization: D.G., A.V., A.A.P., and E.T.; Statistical Methodology: E.T. and J.R.R.; Other Methodology: D.G., A.V., and A.A.P.; Statistical Analysis: E.T., J.R.R., and B.M.B.-H.; Writing-original manuscript preparation and figures: B.M.B.-H.; Writing-review and editing, B.M.B.-H., E.T., J.R.R., D.G., A.V., and A.A.P.; Electronic medical record data obtained by A.V.; Project supervision: D.G. and E.T.; Air Quality data obtained by A.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Heinz Endowments; grant numbers E2157, E6462 and E6450.

Acknowledgments

Thank you to Avery Hill for assisting with the creation of the manuscript figures.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

1. Hamanaka, R.B.; Mutlu, G.M. Particulate Matter Air Pollution: Effects on the Cardiovascular System. Front. Endocrinol. Lausanne 2018, 9, 680.

2. Chen, R.; Peng, R.D.; Meng, X.; Zhou, Z.; Chen, B.; Kan, H. Seasonal variation in the acute effect of particulate air pollution on mortality in the China Air Pollution and Health Effects Study (CAPES). Sci. Total Environ. 2013, 450, 259-265.

3. Kuo, C.Y.; Chan, C.K.; Wu, C.Y.; Phan, D.; Van Chan, C.L. The short-term effects of ambient air pollutants on childhood asthma hospitalization in Taiwan: A national study. Int. J. Environ. Res. Public Health 2019, 16, 203.

4. Anderson, H.R.; de Leon, A.P.; Bland, J.M.; Bower, J.S.; Emberlin, J.; Strachan, D.P. Air pollution, pollens, and daily admissions for asthma in London 1987-92. Thorax 1998, 53, 842-848.

5. Peel, J.L.; Tolbert, P.E.; Klein, M.; Metzger, K.B.; Flanders, W.D.; Todd, K.; Mulholland, J.A.; Ryan, P.B.; Frumkin, H. Ambient air pollution and respiratory emergency department visits. Epidemiology 2005, 16, 164-174.

6. Strickland, M.J.; Darrow, L.A.; Klein, M.; Flanders, W.D.; Sarnat, J.A.; Waller, L.A.; Sarnat, S.E.; Mulholland, J.A.; Tolbert, P.E. Short-term associations between ambient air pollutants and pediatric asthma emergency department visits. Am. J. Respir. Crit. Care Med. 2010, 182, 307-316.

7. Barnett, G.A.; Williams, M.G.; Schwartz, J.; Neller, H.A.; Best, L.T.; Petroeschevsky, L.A.; Simpson, W.R. Air pollution and child respiratory health: A case-crossover study in Australia and New Zealand. Am. J. Respir. Crit. Care Med. 2005, 171, 1272-1278.

8. Zheng, X.Y.; Ding, H.; Jiang, L.N.; Chen, S.W.; Zheng, J.P.; Qiu, M.; Zhou, Y.X.; Chen, Q.; Guan, W.J. Association between Air pollutants and asthma emergency room visits and hospital admissions in time series studies: A systematic review and meta-Analysis. PLoS ONE 2015, 10, e0138146.

9. Son, J.Y.; Lee, J.T.; Park, Y.H.; Bell, M.L. Short-term effects of air pollution on hospital admissions in Korea. Epidemiology 2013, 24, 545-554.

10. Malig, B.J.; Pearson, D.L.; Chang, Y.B.; Broadwin, R.; Basu, R.; Green, R.S.; Ostro, B. A time-stratified case-crossover study of ambient ozone exposure and emergency department visits for specific respiratory diagnoses in California (2005-2008). Environ. Health Perspect. 2016, 124, 745-753.

11. Pereira, A.A.; Pollard, S.L.; Locke, R.; Romero, K.; Lima, J.J.; Hansel, N.N.; Checkley, W.; GASP Study Investigators. Association between exhaled carbon monoxide and asthma outcomes in Peruvian children. Respir. Med. 2018, 145, 212-216.

12. Lin, M.; Chen, Y.; Burnett, R.T.; Villeneuve, P.J.; Krewski, D. Effect of short-term exposure to gaseous pollution on asthma hospitalisation in children: A bi-directional case-crossover analysis. J. Epidemiol. Community Health 2003, 57, 50-55.

13. Goodman, J.E.; Kennedy, E.M.; Seeley, M. Do individuals with asthma experience airway hyper-responsiveness after exposure to nitrogen dioxide? Regul. Toxicol. Pharmacol. 2017, 89, 279-287.

14. Glad, J.A.; Brink, L.L.; Talbott, E.O.; Lee, P.C.; Xu, X.; Saul, M.; Rager, J. The relationship of ambient ozone and PM2.5 levels and asthma emergency department visits: Possible influence of gender and ethnicity. Arch. Environ. Occup. Heal. 2012, 67, 103-108.

15. Pennsylvania-State of the Air | American Lung Association. Available online: http://www.stateoftheair.org/city-rankings/states/pennsylvania/ (accessed on 6 April 2020).

16. Gharibi, H.; Entwistle, M.R.; Ha, S.; Gonzalez, M.; Brown, P.; Schweizer, D.; Cisneros, R. Ozone pollution and emergency department visits in the Central Valley, California, USA, during June to September of 2015: A time-stratified case-crossover analysis. J. Asthma 2019, 56, 1037-1048.

17. Fauroux, B.; Sampil, M.; Quénel, P.; Lemoullec, Y. Ozone: A trigger for hospital pediatric asthma emergency room visits. Pediatr. Pulmonol. 2000, 30, 41-46.

18. Gleason, J.A.; Fagliano, J.A. Associations of daily pediatric asthma emergency department visits with air pollution in Newark, NJ: Utilizing time-series and case-crossover study designs. J. Asthma 2015, 52, 815-822.

19. Jaakkola, J.J.K. Case-crossover design in air pollution epidemiology. Eur. Respir. J. 2003, 21, 81-85.

20. Maclure, M. The Case-Crossover Design: A Method for Studying Transient Effects on the Risk of Acute Events. Am. J. Epidemiol. 1991, 133, 144-153.

21. Shmool, J.L.; Michanowicz, D.R.; Cambal, L.; Tunno, B.; Howell, J.; Gillooly, S.; Roper, C.; Tripathy, S.; Chubb, L.G.; Eisl, H.M.; et al. Saturation sampling for spatial variation in multiple air pollutants across an inversion-prone metropolitan area of complex terrain. Environ. Heal. A Glob. Access Sci. Source 2014, 13, 1-16.

22. Larsen, L.C.; Shah, M. A context-intensive approach to imputation of missing values in data sets from networks of environmental monitors. J. Air Waste Manage. Assoc. 2016, 66, 38-52.

23. Junninen, H.; Niska, H.; Tuppurainen, K.; Ruuskanen, J.; Kolehmainen, M. Methods for imputation of missing values in air quality data sets. Atmos. Environ. 2004, 38, 2895-2907.

24. Plaia, A.; Bondì, A.L. Single imputation method of missing values in environmental pollution data sets. Atmos. Environ. 2006, 40, 7316-7330.

25. Le, H.Q.; Batterman, S.A.; Wahl, R.L. Reproducibility and imputation of air toxics data. J. Environ. Monit. 2007, 9, 1358-1372.

26. Teague, W.G.; Phillips, B.R.; Fahy, J.V.; Wenzel, S.E.; Fitzpatrick, A.M.; Moore, W.C.; Hastie, A.T.; Bleecker, E.R.; Meyers, D.A.; Peters, S.P.; et al. Baseline Features of the Severe Asthma Research Program (SARP III) Cohort: Differences with Age. J. Allergy Clin. Immunol. Pract. 2018, 6, 545-554.e4.

27. Hopke, P.K.; Croft, D.; Zhang, W.; Lin, S.; Masiol, M.; Squizzato, S.; Thurston, S.W.; van Wijngaarden, E.; Utell, M.J.; Rich, D.Q. Changes in the acute response of respiratory diseases to PM 2.5 in New York State from 2005 to 2016. Sci. Total Environ. 2019, 677, 328-339.

28. Anenberg, S.C.; Henze, D.K.; Tinney, V.; Kinney, P.L.; Raich, W.; Fann, N.; Malley, C.S.; Roman, H.; Lamsal, L.; Duncan, B.; et al. Estimates of the global burden of ambient PM2:5, ozone, and NO2 on asthma incidence and emergency room visits. Environ. Health Perspect. 2018, 126, 107004.

29. Krall, R.J.; Mulholland, A.J.; Russell, G.A.; Balachandran, S.; Winquist, A.; Tolbert, E.P.; Waller, A.L.; Sarnat, E.S. Associations between source-specific fine particulate matter and emergency department visits for respiratory disease in four U.S. cities. Environ. Health Perspect. 2017, 125, 97-103.

30. Esposito, S.; Tenconi, R.; Lelii, M.; Preti, V.; Nazzari, E.; Consolo, S.; Patria, M.F. Possible molecular mechanisms linking air pollution and asthma in children. BMC Pulm. Med. 2014, 14, 31.

31. Kuo, C.Y.; Pan, R.H.; Chan, C.K.; Wu, C.Y.; Phan, D.V.; Chan, C.L. Application of a time-stratified case-crossover design to explore the effects of air pollution and season on childhood asthma hospitalization in cities of differing urban patterns: Big data analytics of government open data. Int. J. Environ. Res. Public Health 2018, 15, 647.

32. Entwistle, M.R.; Gharibi, H.; Tavallali, P.; Cisneros, R.; Schweizer, D.; Brown, P.; Ha, S. Ozone pollution and asthma emergency department visits in Fresno, CA, USA, during the warm season (June-September) of the years 2005 to 2015: A time-stratified case-crossover analysis. Air Qual. Atmos. Health 2019, 1-12.

33. Analitis, A.; De'Donato, F.; Scortichini, M.; Lanki, T.; Basagana, X.; Ballester, F.; Astrom, C.; Paldy, A.; Pascal, M.; Gasparrini, A.; et al. Synergistic effects of ambient temperature and air pollution on health in europe: Results from the PHASE project. Int. J. Environ. Res. Public Health 2018, 15, 1856.

34. Chen, K.; Glonek, G.; Hansen, A.; Williams, S.; Tuke, J.; Salter, A.; Bi, P. The effects of air pollution on asthma hospital admissions in Adelaide, South Australia, 2003-2013: Time-series and case-crossover analyses. Clin. Exp. Allergy 2016, 46, 1416-1430.

35. Bateson, T.F.; Schwartz, J. Children's response to air pollutants. J. Toxicol. Environ. Health Part A Curr. Issues 2008, 71, 238-243.

36. Mirowsky, J.E.; Dailey, L.A.; Devlin, R.B. Differential expression of pro-inflammatory and oxidative stress mediators induced by nitrogen dioxide and ozone in primary human bronchial epithelial cells. Inhal. Toxicol. 2016, 28, 374-382.

37. Epidemiology, Pathophysiology, Clinical Evaluation, and Treatment of Carbon Monoxide Poisoning in Child, Infant, and Fetus. Available online: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5530151/ (accessed on 21 April 2020).

38. Ryter, S.W.; Choi, A.M.K. Carbon Monoxide in Exhaled Breath Testing and Therapeutics. J. Breath Res. 2013, 7, 017111.

39. Canova, C.; Torresan, S.; Simonato, L.; Scapellato, M.L.; Tessari, R.; Visentin, A.; Lotti, M.; Maestrelli, P. Carbon monoxide pollution is associated with decreased lung function in asthmatic adults. Eur. Respir. J. 2010, 35, 266-272.

40. Halpin, D.M.; Meltzer, E.O.; Pisternick-Ruf, W.; Moroni-Zentgraf, P.; Engel, M.; Zaremba-Pechmann, L.; Casale, T.; FitzGerald, J.M. Peak expiratory flow as an endpoint for clinical trials in asthma: A comparison with FEV1. Respir. Res. 2019, 20, 159.

41. Evans, K.A.; Halterman, J.S.; Hopke, P.K.; Fagnano, M.; Rich, D.Q. Increased ultrafine particles and carbon monoxide concentrations are associated with asthma exacerbation among urban children. Environ. Res. 2014, 129, 11-19.

42. Chevalier, A.; Gheusi, F.; Delmas, R.; Ordóñez, C.; Sarrat, C.; Zbinden, R.; Thouret, V.; Athier, G.; Cousin, J.M. Atmospheric Chemistry and Physics Influence of altitude on ozone levels and variability in the lower troposphere: A ground-based study for western Europe over the period 2001-2004. Atmos. Chem. Phys 2007, 7, 4326.