Skip to main content
Advertisement
  • Loading metrics

Ambient particulate matter pollution and adult hospital admissions for pneumonia in urban China: A national time series analysis for 2014 through 2017

  • Yaohua Tian,

    Roles Conceptualization, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Maternal and Child Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China, Ministry of Education Key Laboratory of Environment and Health, and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China, Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Hui Liu,

    Roles Formal analysis, Visualization, Writing – review & editing

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Medical Informatics Center, Peking University, Beijing, China

  • Yiqun Wu,

    Roles Validation, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Yaqin Si,

    Roles Data curation, Investigation, Resources, Writing – review & editing

    Affiliation Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Man Li,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Yao Wu,

    Roles Software, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Xiaowen Wang,

    Roles Validation, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Mengying Wang,

    Roles Methodology, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Libo Chen,

    Roles Data curation, Investigation, Resources, Writing – review & editing

    Affiliation Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Chen Wei,

    Roles Data curation, Investigation, Resources, Supervision, Writing – review & editing

    Affiliation Beijing HealthCom Data Technology Co. Ltd, Beijing, China

  • Tao Wu,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China

  • Pei Gao ,

    Contributed equally to this work with: Pei Gao, Yonghua Hu

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – review & editing

    yhhu@bjmu.edu.cn (YH); peigao@bjmu.edu.cn (PG)

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Key Laboratory of Molecular Cardiovascular (Peking University), Ministry of Education, Beijing, China

  • Yonghua Hu

    Contributed equally to this work with: Pei Gao, Yonghua Hu

    Roles Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review & editing

    yhhu@bjmu.edu.cn (YH); peigao@bjmu.edu.cn (PG)

    Affiliations Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China, Medical Informatics Center, Peking University, Beijing, China

Abstract

Background

The effects of ambient particulate matter (PM) pollution on pneumonia in adults are inconclusive, and few scientific data on a national scale have been generated in low- or middle-income countries, despite their much higher PM concentrations. We aimed to examine the association between PM levels and hospital admissions for pneumonia in Chinese adults.

Methods and findings

A nationwide time series study was conducted in China between 2014 and 2017. Information on daily hospital admissions for pneumonia for 2014–2017 was collected from the database of Urban Employee Basic Medical Insurance (UEBMI), which covers 282.93 million adults. Associations of PM concentrations and hospital admissions for pneumonia were estimated for each city using a quasi-Poisson regression model controlling for time trend, temperature, relative humidity, day of the week, and public holiday and then pooled by random-effects meta-analysis. Meta-regression models were used to investigate potential effect modifiers, including cities’ annual-average air pollutants concentrations, temperature, relative humidity, gross domestic product (GDP) per capita, and coverage rates by the UEBMI. More than 4.2 million pneumonia admissions were identified in 184 Chinese cities during the study period. Short-term elevations in PM concentrations were associated with increased pneumonia admissions. At the national level, a 10-ΞΌg/m3 increase in 3-day moving average (lag 0–2) concentrations of PM2.5 (PM ≀2.5 ΞΌm in aerodynamic diameter) and PM10 (PM ≀10 ΞΌm in aerodynamic diameter) was associated with 0.31% (95% confidence interval [CI] 0.15%–0.46%, P < 0.001) and 0.19% (0.11%–0.30%, P < 0.001) increases in hospital admissions for pneumonia, respectively. The effects of PM10 were stronger in cities with higher temperatures (percentage increase, 0.031%; 95% CI 0.003%–0.058%; P = 0.026) and relative humidity (percentage increase, 0.011%; 95% CI 0%–0.022%; P = 0.045), as well as in the elderly (percentage increase, 0.10% [95% CI 0.02%–0.19%] for people aged 18–64 years versus 0.32% [95% CI 0.22%–0.39%] for people aged β‰₯75 years; P < 0.001). The main limitation of the present study was the unavailability of data on individual exposure to PM pollution.

Conclusions

Our findings suggest that there are significant short-term associations between ambient PM levels and increased hospital admissions for pneumonia in Chinese adults. These findings support the rationale that further limiting PM concentrations in China may be an effective strategy to reduce pneumonia-related hospital admissions.

Author summary

Why was this study done?

  • Epidemiological studies have reported associations between short-term exposure to ambient particulate matter (PM) pollution and the risk of pneumonia.
  • Previous studies have been primarily conducted in high-income countries, and the findings remain inconclusive.
  • Few scientific data on a national scale have been generated in low- or middle-income countries, despite their much higher PM concentrations.

What did the researchers do and find?

  • We conducted a nationwide time series analysis using data on more than 4.2 million hospital admissions for pneumonia in 184 cities in China between 2014 and 2017 to estimate city-specific, national, and regional average associations between ambient PM pollution and pneumonia hospitalizations.
  • Our results suggested that short-term increases in PM2.5 and PM10 were associated with increased hospital admissions for pneumonia. The effects of PM10 were stronger in cities with higher temperatures and relative humidity, as well as in the elderly.

What do these findings mean?

  • To our knowledge, this is the first study in China to investigate the short-term associations of PM levels with hospital admissions for pneumonia on a national scale.
  • Our findings support the rationale for further limiting PM concentrations in low- and middle-income countries.

Introduction

Pneumonia is a major cause of mortality and morbidity worldwide [1, 2]. In China, an estimated 2.5 million pneumonia cases occur annually, and 5% of these individuals die of pneumonia-related illness [3]. Furthermore, pneumonia is closely related to complications such as pleurisy, lung abscess, and septicemia, as well as to cardiovascular disease [4, 5]. The risks of pneumonia and pneumonia-related mortality increase with age [6]; thus, the incidence of pneumonia is projected to increase as a result of global population aging [7].

Air pollution has emerged as a significant public health problem worldwide [8]. Despite plausible hypotheses linking air pollution with pneumonia [9–12], limited evidence of the association is available, especially in adults. Some studies have started to report some evidence of an association between ambient particulate matter (PM) pollution and respiratory infections, including pneumonia. A time series study suggested an association between fine particulate matter (PM2.5, PM ≀2.5 ΞΌm in aerodynamic diameter) and increased admissions for respiratory tract infections in Medicare enrollees (aged >65 years) in 204 United States urban counties [13]. Similarly, another study reported significant effects of PM10 (PM ≀10 ΞΌm in aerodynamic diameter) on pneumonia admissions among 36 US citizens aged 65 years or older [14]. Strosnider and colleagues assessed age-specific acute effects of PM2.5 on respiratory emergency department visits in 17 US states and reported that short-term exposure to PM2.5 was associated with increased emergency department visits for pneumonia in people aged 19–65 years, but not in people aged <19 or β‰₯65 years [15]. A recent case-crossover study from the US demonstrated significant associations between short-term elevations in PM2.5 and greater healthcare utilization for acute lower respiratory infection in both children and adults [16]. Host and colleagues found increased admissions for lower respiratory tract infections following PM2.5 exposure in six French cities [17], and two time series studies in Hong Kong provided further evidence on the acute effects of PM on respiratory tract infections/pneumonia [18, 19]. However, these studies were conducted primarily in high-income countries/cities where characteristics of air pollution and socioeconomic status differ from those in low- or middle-income countries. The effects of short-term exposure to PM pollution on pneumonia in low- or middle-income counties require further investigation.

China, the largest low- or middle-income country, has one of the highest PM2.5 concentrations worldwide [20]. However, only a few studies have explored the associations between PM pollution and pneumonia, and were subjected to single city/hospital study and small sample size [21, 22]. In this study, we examined the short-term associations between concentrations of ambient PM pollution and daily hospital admissions for pneumonia in adults in China between 2014 and 2017.

Methods

Study sites

A total of 184 Chinese cities were finally included in this study. Fig 1 shows the locations of the 184 cities, representing a geographic distribution across China. Cities were selected according to the availability of daily hospital admission and air pollution data. Cities with ≀1-year hospital admission records were excluded given the fact that a long enough study period in a time series analysis is needed to ensure credible precision and power. Individuals’ detailed disease diagnosis (International Classification of Diseases [ICD] code and text of disease diagnosis) was required to identify pneumonia admissions. Cities with no information on disease diagnosis recorded in the database were also excluded (e.g., Beijing and Shanghai). In addition, we applied additional rules to control the quality of the original data. We regressed the city-specific average numbers of pneumonia admission per day on the number of enrollers of the city. Cities with an absolute value of standardized residuals >3 were considered as the outliers and were removed from the analysis. China has gradually introduced PM2.5 in the national air quality monitoring network and publicized real-time monitoring data since 2013. In this study, the total study period encompassed 2014–2017, and the years of analysis differed by city based on the availability of air pollution data. Of the 184 cities, 78 cities had 3-year PM data records and 106 cities had 4-year PM data records. The years of the study period for each city were presented in S1 Table.

thumbnail
Fig 1. Locations of the 184 Chinese cities included in the study.

The base map was obtained from Natural Earth (https://naturalearthdata.com).

https://doi.org/10.1371/journal.pmed.1003010.g001

Study population

In China, there are three main health insurance schemes: the Urban Employee Basic Medical Insurance (UEBMI), the Urban Residence Basic Medical Insurance, and the New Rural Cooperative Medical Scheme, covering more than 92% of the population [23]. Private medical insurance covers little in China and is generally supplementary to the basic schemes. Daily pneumonia admission data were obtained from the UEBMI database. The UEBMI covers urban employees and retired employees (aged β‰₯18 years). All employers in urban areas, including government agencies and institutions, state-owned enterprises, private businesses, social organizations, and other private entities and their employees (retirees included), are obligated to enroll in UEBMI [24]. At the end of 2016, 282.93 million individuals were recorded in the database. The original data source was medical claims data for urban employees, which included sex, age, date of admission, and cause of admission. Information regarding the number of individuals enrolled in the database, city residents, and coverage rates of these cities for the UEBMI database were published previously [25]. We have legal access to this national insurance information through the big data platform of national medical insurance data for risk prevention and control, established by the Ministry of Human Resources and Social Security of the People’s Republic of China and Beijing HealthCom Data Technology. Hospitalizations with a principal diagnosis of pneumonia (ICD-10 codes J12–J18) were identified from the database. Because the health information was primarily the city-specific daily count of pneumonia admissions, i.e., summarized data (overall, and by age and sex subgroups) without any individual identifiers, this study was exempted from Institutional Review Board approval by the Ethics Committee of Peking University Health Science Center, Beijing, China. The need for informed consent was also waived by the Institutional Review Board. This study is reported as per the STROBE guideline (S1 STROBE checklist).

Environmental data

Data on hourly PM2.5 and PM10 measurements in each city were obtained from the National Air Pollution Monitoring System. Each city has 1–17 monitoring stations. China has issued detailed standards on placement of monitoring stations and the air pollution monitoring process [26]; the monitoring stations cannot be located in the direct vicinity of apparent emission sources; thus, measurements reflect the general urban background concentrations. These measurements have therefore been used extensively to represent population exposure [27–29]. Data on sulfur dioxide (SO2), nitrogen dioxide (NO2), carbon monoxide (CO), and ozone (O3) concentrations were also collected from the same source. We calculated 24-hour averages by averaging hourly air pollution measurements within a day. In each city, the measurements from all monitoring stations were averaged to represent population exposure [27, 28, 30]. Rates of missing data were 0.66% for PM2.5 and 0.15% for PM10. Days with missing monitoring measurements were excluded from the analysis. Data on hourly atmospheric temperature and relative humidity were provided by the China Meteorological Data Sharing Service System. Each city has 1–3 meteorological monitoring stations, all located in urban areas. Rates of missing temperature and relative humidity data were only 0.25%. Daily (24-hour) mean temperature and relative humidity were calculated by averaging all valid monitoring measurements in each city. More detailed information about health, air pollution, and meteorological data sources is given in our previous publications [24, 25].

Statistical analysis

We used a common two-stage method to examine the associations between PM and pneumonia admissions [13, 24, 27, 30–32]. The method and the model used in this study were designed before the analyses were conducted, and the prespecified statistical analysis plan is present in S1 Appendix. In stage 1, we assessed the city-specific associations between PM and pneumonia admissions using a quasi-Poisson regression model that allows for over-dispersed admission counts. Several covariates were included in the model, following previous published studies [13, 27, 32, 33]: (1) a natural cubic spline for time with 7 degrees of freedom (df) per year to adjust for seasonality and long-term trends, (2) natural cubic splines for 3-day moving average temperature with 6 df and relative humidity with 3 df to control for nonlinear and delayed effects of meteorology, and (3) indicators of day of the week and public holiday. The selection of df values was based on the parameter used in several recent large national studies [27, 28, 32, 33]. The model was as shown below: where E(Yt) is the expected number of pneumonia admissions on day t; Ξ² indicates the log-relative risk of pneumonia associated with a unit increase of PM concentrations; ns() indicates a natural cubic spline function; and temperature and relative humidity indicate 3-day moving averages. This model was in line with those used in previous studies [13, 24, 27, 30, 31] and was tested in sensitivity analyses, described in detail below. We also conducted unadjusted analyses with only PM and a natural cubic spline for time, with 7 df per year in the regression model. In stage 2, we conducted random-effects meta-analyses to combine the city-specific effect estimates. To explore regional heterogeneity of the associations, we divided the 184 cities into six geographical regions, i.e., East, Middle-south, Southwest, Northwest, North, and Northeast regions [27]. We used single lags of 0, 1, and 2 days to assess the temporal association between PM and pneumonia. Considering that single-day lag models may underestimate PM effects, we also estimated associations with 3-day (lag 0–2) moving average PM concentrations.

We plotted the national-average concentration-response curves of the associations of PM with pneumonia admissions using the approach of the Air Pollution and Healthβ€”A European Approach (APHEA) project [27, 32, 34]. Briefly, we replaced the linear term of PM in the main model with a B-spline function with two knots at 60 and 150 ΞΌg/m3 for PM2.5 and 100 and 200 ΞΌg/m3 for PM10, based on the distribution of PM concentrations in each city. We then estimated five regression coefficients of the spline function and the 5 Γ— 5 variance-covariance matrix in each city. Finally, we fitted multivariable random-effect models to combine the city-specific components of spline estimates. Variance-covariance matrix is a matrix whose element in the i, j position is the covariance between the i-th and j-th elements of a random vector. The regression coefficients derived from the main model with a B-spline function for PM.

We examined effect modification of the relationship between short-term ambient PM concentrations and pneumonia hospital admissions in analyses by sex, age, and region (northern and southern China). The regional division followed the Huai River-Qinling Mountains line [27, 30]. We assessed the differences in the estimates using a Z-test [35]. We also fitted meta-regression models to investigate several city characteristics as potential effect modifiers, including cities’ annual-average air pollutants concentrations, air temperature and relative humidity, gross domestic product (GDP) per capita, and coverage rates by the UEBMI [27, 30]. Cities’ annual-average air pollutants concentrations, temperature, and relative humidity were calculated from daily measurements during the study period. City-specific relative risks and their confidence intervals (CIs) as the outcome were meta-regressed on each continuous variable of city characteristics. We fitted the random-effects meta-regression model using city-level estimations from the first stage as yi = Ξ²xi + ΞΌi + ei, where yi is the city-specific estimation, xi is the city-level characteristic variable, ΞΌi is the city-specific random effect characterized by the between-city variance, and ei is a normal error term [36]. Random-effects meta-regression can be considered as an extension to the random-effects meta-analysis (i.e., DerSimonian and Laird model) that includes study-level (i.e., city-level in our case) covariates. The model is implemented by the metareg function in STATA. The algorithm for random-effects meta-regression firstly estimated the between-study (between-city) variance and then estimated the coefficients, Ξ², by weighted least squares by using the weights 1/(Οƒi2 + Ο„2), where Οƒi2 is the standard error of the estimated effect in study (city) i, and Ο„2 is the between-study (between-city) variance, advocated by Thompson and Sharp [37].

Sensitivity analysis

We conducted several sensitivity analyses: (1) We tested the potential confounding effects of gaseous pollutants in two-pollutant models; (2) we separately examined the associations of PM with pneumonia in cities with only 3- or 4-year data to compare with the primary analysis; (3) we examined the associations between PM and pneumonia in cities with different levels of population coverage by the UEBMI (<20% and β‰₯20%); (4) we examined the associations excluding cities with ≀2 monitoring stations; (5) we conducted a sensitivity analysis with additional adjustment of the hospitalization of influenza; and (6) to examine whether we appropriately specified the regression model used in this analysis, we used alternative df for time (6–10 per year) or used penalized spline functions for time and meteorological variables in the model.

Statistical analyses were conducted in R version 3.2.2 (R Development Core Team 2008) and STATA software, version 12 (StataCorp, College Station, TX), and a two-sided P < 0.05 was considered statistically significant. The effect estimates are presented as percentage increases and 95% CIs in pneumonia admissions per 10-ΞΌg/m3 increase in PM concentrations.

Results

There were 4.2 million hospital admissions for pneumonia in the 184 Chinese cities between 2014 and 2017 in our study. Table 1 presents summary statistics on citywide pneumonia admissions, PM concentrations, and weather conditions at an annual-average level during the study period. The mean number of daily hospital admissions for pneumonia was 26 (range: 1 to 376). The national-average PM2.5 concentration was 50 ΞΌg/m3 (range: 15 to 102 ΞΌg/m3), and the mean PM10 concentration was 89 ΞΌg/m3 (range: 28 to 193 ΞΌg/m3). The average citywide annual-mean temperature was 14Β°C (range: 0 to 24Β°C). City-specific characteristics, including the number of individuals enrolled in the database, city residents, coverage rates of the population by UEBMI database, and annual-average PM2.5 and PM10 concentrations, temperature, relative humidity, and daily pneumonia hospitalizations are presented in S1 Table.

thumbnail
Table 1. Summary statistics of citywide annual-mean hospital admissions for pneumonia, air pollutants, and weather conditions in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.t001

Table 2 shows the Spearman correlation coefficient values for the environmental variables. At the national level, daily PM2.5 and PM10 concentrations were positively correlated with SO2 (PM2.5: r = 0.56, P < 0.001; PM10: r = 0.59, P < 0.001), NO2 (PM2.5: r = 0.64, P < 0.001; PM10: r = 0.64, P < 0.001), and CO (PM2.5: r = 0.61, P < 0.001; PM10: r = 0.57, P < 0.001) concentrations. O3 concentrations were not correlated with PM2.5 (r = βˆ’0.02, P = 0.452) or PM10 (r = 0.05, P = 0.215). There were inverse and weak correlations between PM2.5 and PM10 concentrations and temperature (PM2.5: r = βˆ’0.26, P < 0.001; PM10: r = βˆ’0.26, P < 0.001) and relative humidity (PM2.5: r = βˆ’0.08, P < 0.001; PM10: r = βˆ’0.28, P < 0.001).

thumbnail
Table 2. Spearman correlation coefficients among the environmental variables in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.t002

City-specific estimates of the associations between 3-day moving average (lag 0–2) concentrations of PM2.5 and PM10 and hospital admissions for pneumonia are presented in S2 Table. Table 3 shows the national-average percentage increases in daily hospital admissions for pneumonia associated with a 10-ΞΌg/m3 increase in PM2.5 and PM10 concentrations at lags 0, 1, 2, and 0–2 days. Overall, we observed similar lag patterns for PM2.5 and PM10. For single-day lags, the strongest associations occurred at lag 0. The associations weakened considerably at lag days 1 and 2. Lag 0–2 generated the highest estimates. On average, a 10-ΞΌg/m3 increase in PM2.5 and PM10 concentrations at lag 0–2 corresponded to a 0.31% (0.15% to 0.46%, P < 0.001) and 0.19% (0.11% to 0.30%, P < 0.001) increase in pneumonia admissions, respectively.

thumbnail
Table 3. National-average percentage increase with 95% CI in daily hospital admissions for pneumonia associated with a 10 ΞΌg/m3 increase in PM2.5 and PM10 concentrations using different lag days in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.t003

Table 4 presents the results of stratified analyses. The estimates were similar in males (PM2.5: percentage increase, 0.31%; 95% CI 0.12%–0.49%; PM10: percentage increase, 0.20%; 95% CI 0.13%–0.28%) and females (PM2.5: percentage increase, 0.33%; 95% CI 0.13%–0.53%; PM10: percentage increase, 0.15%; 95% CI 0.05%–0.25%) (the P values for the differences in the estimates were 0.887 for PM2.5 and 0.422 for PM10). The associations between PM and pneumonia were stronger in people aged β‰₯75 years, but the difference between the strata was only significant for PM10 (the P values for the differences in the estimates were 0.201 for PM2.5 and P < 0.001 for PM10). The effect estimates were higher in the southern region (PM2.5: percentage increase, 0.40%; 95% CI 0.16%–0.65%; PM10: percentage increase, 0.35%; 95% CI 0.16%–0.54%) than in the northern region (PM2.5: percentage increase, 0.23%; 95% CI 0.04%–0.43%; PM10: percentage increase, 0.12%; 95% CI 0%–0.24%), but the between-region difference was not significant for PM2.5 (the P values for the regional differences in the estimates were 0.276 for PM2.5 and 0.045 for PM10). We grouped the cities further into six geographical areas; the average estimates for the six areas are presented in S3 Table. There was a significant heterogeneity in the PM–hospitalization associations across different regions. The associations were more evident in the Middle-south, East, and North regions.

thumbnail
Table 4. National-average percentage increase (PI) with 95% CI in daily hospital admissions for pneumonia associated with a 10 ΞΌg/m3 increase in PM2.5 and PM10 concentrations (lag 0–2), stratified by sex, age, and geographical region.

https://doi.org/10.1371/journal.pmed.1003010.t004

We noted a slightly nonlinear concentration-response curve between PM2.5 and pneumonia admissions, where there appears to be a plateau at higher concentrations. This evidence of linearity is consistent with studies that examined short-term PM exposure and other outcomes [13–15, 19, 27]. The curve of the association between PM10 and pneumonia increased sharply at concentrations below 100 ΞΌg/m3 and then climbed relatively moderately as concentrations increased (Fig 2). The estimated risk of PM10 concentrations for pneumonia was 0.38% (0.18% to 0.57%, P < 0.001) at PM10 < 100 ΞΌg/m3 and 0.10% (βˆ’0.10% to 0.30%, P = 0.754) at PM10 β‰₯ 100 ΞΌg/m3.

thumbnail
Fig 2. National-average concentration-response curves of 3-day moving average (lag 0–2) concentrations of PM2.5 and PM10 and hospital admissions for pneumonia in 184 cities in China, 2014–2017.

Note, y-axes of the two graphs are scaled differently. The solid line represents relative changes and the dashed lines represent the 95% CIs. The vertical lines represent the air quality guidelines or standards for 24-hour average concentrations of PM2.5 and PM10. PM2.5: 25 ΞΌg/m3 is the World Health Organization (WHO) air quality guideline for daily PM2.5 concentrations, 35 ΞΌg/m3 is the Chinese Grade I standard for daily PM2.5 concentrations, and 75 ΞΌg/m3 is the Chinese Grade II standard; PM10: 50 ΞΌg/m3 is both the WHO air quality guideline and Chinese Grade I standard for daily PM10 concentrations, and 150 ΞΌg/m3 is the Chinese Grade II standard for daily PM10 concentrations. PM2.5, particulate matter ≀2.5 ΞΌm in aerodynamic diameter; PM10, particulate matter ≀10 ΞΌm in aerodynamic diameter.

https://doi.org/10.1371/journal.pmed.1003010.g002

Table 5 summarizes the results of meta-regression analyses of effect modification on the associations between PM and pneumonia by city-level characteristics. We observed stronger associations between PM10 and pneumonia in cities with higher annual-average temperatures or relative humidity. For each 10-ΞΌg/m3 increase in PM10 concentrations, a city with 1Β°C higher annual-mean temperature with respect to another city would see an additional 0.031% (0.003% to 0.058%, P = 0.026) increase in pneumonia admissions; a city with 1% higher annual-mean relative humidity with respect to another city would see an additional 0.011% (0% to 0.022%, P = 0.045) increase in pneumonia admissions. We observed the same direction of effect modification by temperature (P = 0.129) and relative humidity (P = 0.378) on the association between PM2.5 and pneumonia, but the effects were not significant. No evidence was found for effect modification by annual-average PM2.5 (percentage increase, βˆ’0.019%; 95% CI βˆ’0.113% to 0.075%, P = 0.688) or PM10 (percentage increase, 0; 95% CI βˆ’0.064% to 0.064%, P = 0.997). Annual-average NO2 (the P values were 0.737 for PM2.5 and P = 0.861 for PM10) or CO (the P values were 0.466 for PM2.5 and P = 0.367 for PM10) concentrations, GDP per capita (the P values were 0.864 for PM2.5 and P = 0.655 for PM10), or the coverage rate of the population (the P values were 0.835 for PM2.5 and P = 0.665 for PM10) did not significantly modify the associations.

thumbnail
Table 5. Meta-regression results of the modification effects of city-level characteristics on the associations between PM2.5 and PM10 and hospital admissions for pneumonia in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.t005

The associations of PM2.5 and PM10 with pneumonia remained after adjusting for SO2 (PM2.5: percentage increase, 0.22%; 95% CI 0.07%–0.37%; P = 0.006; PM10: percentage increase, 0.12%; 95% CI 0.02%–0.22%; P = 0.015) and O3 (PM2.5: percentage increase, 0.31%; 95% CI 0.16%–0.46%; P < 0.001; PM10: percentage increase, 0.20%; 95% CI 0.09%–0.30%; P < 0.001); the effects remained positive but not significant after controlling for CO (PM2.5: percentage increase, 0.12%; 95% CI βˆ’0.03%–0.27%; P = 0.127; PM10: percentage increase, 0.06%; 95% CI βˆ’0.04%–0.16%; P = 0.256) and NO2 (PM2.5: percentage increase, 0.10%; 95% CI βˆ’0.05%–0.25%; P = 0.271; PM10: percentage increase, 0.04%; 95% CI βˆ’0.06%–0.14%; P = 0.431) (S4 Table).

S5 Table lists the results of the sensitivity analyses. Consistent associations were observed in datasets with different lengths of data-years. Specifically, for PM2.5, the increases in pneumonia admissions were 0.28% (0% to 0.56%, P = 0.044) in 78 cities with 3-year data and 0.37% (0.20% to 0.56%, P < 0.001) in 106 cities with 4-year data; for PM10, the corresponding values were 0.13% (βˆ’0.07% to 0.34%, P = 0.213) in cities with 3-year data and 0.24% (0.11% to 0.37%, P < 0.001) in cities with 4-year data. The estimates were similar in cities with different levels of population coverage (<20% and β‰₯20%). The estimates were slightly attenuated after excluding cities with ≀2 monitoring stations (PM2.5: percentage increase, 0.30%; 95% CI 0.15%–0.46%, P < 0.001; PM10: percentage increase, 0.21%; 95% CI 0.10%–0.32%, P < 0.001) or after adjustment of hospitalization of influenza (PM2.5: percentage increase, 0.29%; 95% CI 0.13%–0.45%, P < 0.001; PM10: percentage increase, 0.18%; 95% CI 0.07%–0.29%, P < 0.001). The national-average estimates generated from models with alternative df for time (6–10 per year, all P < 0.05) or with penalized spline functions for time and meteorological variables (PM2.5: percentage increase, 0.31%; 95% CI 0.17%–0.45%, P < 0.001; PM10: percentage increase, 0.22%; 95% CI 0.12%–0.32%, P < 0.001) were comparable to those generated from the base model.

Discussion

Using comprehensive data, we examined the associations between ambient PM pollution and pneumonia hospitalizations in 184 Chinese cities. We also investigated potential effect modifications by a variety of demographic, geographical, and meteorological characteristics. We found that short-term elevations in PM2.5 and PM10 were associated with higher pneumonia-related hospital admissions in Chinese adults. To our knowledge, this is the first study in China, or even in the low- or middle-income countries, to investigate whether short-term changes in ambient PM concentrations are related to increases pneumonia hospital admissions on a national scale.

We observed increased risk of hospital admissions for pneumonia in association with both PM2.5 and PM10 at lag 0 and lag 1, in line with previous findings from studies examining short-term PM exposures [14, 21]. The highest estimates for single-day lag models were observed for lag 0, indicating immediate PM effects on pneumonia. In addition, the estimate at lag 0–2 was slightly higher than at the estimate for lag 0, suggesting that there might be some cumulative effects of PM on pneumonia admissions as well. We also explored the concentration-response curve of the associations between PM and pneumonia hospitalizations. China is one of the most polluted countries worldwide [38]. The average city-specific annual-mean PM2.5 and PM10 concentrations during the study period were 50 ΞΌg/m3 and 89 ΞΌg/m3, respectively, allowing us to assess the concentration-response pattern at high concentrations. The curves plateau at high concentrations, indicating that the increase in risk of pneumonia admission is greater at lower compared with higher concentrations. A similar shape of the association was also observed in other studies [27]. The leveling off at high concentrations might be explained in that people vulnerable to PM exposure may have developed symptoms and sought treatment before PM concentrations reached high concentrations. Other reasons could also be that there may be different health risks associated with different PM sizes and that people avoid spending time outdoors or start wearing face masks when air pollution is severe.

The PM2.5 estimates were higher than those of PM10 at all lag days in this study, in line with previous findings on pneumonia [22, 39]. PM10 is representing exposure to more particles in the coarse range (PM10–2.5) than PM2.5. The sources, composition, and lung deposition patterns of PM10–2.5 vary from those of PM2.5 [40, 41]. The chemical composition of PM differs by size, with more crustal materials in PM10–2.5 and more combustion-related constituents in PM2.5 [31, 41]. It was reported that PM2.5 can penetrate deep into the lungs, reaching the bronchioles and depositing inside the alveoli. In addition, PM2.5 has a larger surface area than PM10 and thus can absorb more toxic substances per unit mass [42].

There are several possible mechanisms linking PM pollution to pneumonia. It has been postulated that air pollution may act as an irritant and evoke defensive responses in the airways, including increased mucus secretion and bronchial hyperreactivity [9]. PM is a potent oxidant that can produce free radicals and cause oxidative stress in lung cells [10]. An animal study demonstrated that short-term exposure to concentrated ambient particles could increase the levels of reactive oxygen species in the lung [11]. PM-induced oxidative stress might impair the cellular defense and immune system, increasing the susceptibility to infection [19]. In addition, PM might enhance the susceptibility to infection through the suppression of the immune response [43].These toxicological evidences were consistent with our findings on the short-term association between ambient PM concentrations and pneumonia hospitalizations. However, it should be noted that our study does not reveal the mechanisms related to PM that may be active in these processes, and thus further study is still required to verify the causality.

In the two-pollutant models, the associations between PM2.5 and PM10 and pneumonia remained positive but were not significant after adjusting for CO and NO2 concentrations. Similarly, a case-crossover study by Tsai and Yang demonstrated a positive association between PM2.5 and hospital admissions for pneumonia in a single-pollutant model, but the estimate decreased and became nonsignificant when controlling for CO and NO2 [44]. We previously reported positive effects of PM2.5 and PM10 on hospital admissions for respiratory disease, but the estimates decreased dramatically and even became negative after adjusting for NO2 [45]. NO2 and CO are closely associated with vehicle exhaust emissions [46, 47]. In China, a substantial fraction of PM2.5 originates in traffic emissions [48]. The high correlation between pollutants in the regression model may have weakened the observed effects, although the CIs of the estimates were not remarkably inflated in the two-pollutant model. Another possible explanation is that the effects of PM may be confounded by those of other air pollutants. Additional studies are warranted to investigate the independent effects of PM on pneumonia.

There were some effect modifiers in the PM-pneumonia associations. First, the elderly had a higher risk of PM-associated pneumonia. Biological functions decline in the elderly, and this population has a higher prevalence of chronic health conditions [9, 49] that may contribute to the increased vulnerability to air pollution. Second, the effects were stronger in cities with higher temperatures (or higher relative humidity, as they are highly correlated). There are several possible explanations. High temperature itself is associated with increased risk of pneumonia [50, 51]. It was reported that high temperature impacts the emission, transportation, dilution, chemical transformation, and deposition of pollutants [52]. Thus, there may exist the interaction between high temperatures and PM. Previous studies also show moderate evidence of the modifying effect of temperature on PM10. A meta-analysis reported that the relative risks for respiratory mortality per 10 ΞΌg/m3 increase in PM10 in the low, middle, and high temperature level were 1.005 (95% CI 1.000–1.010), 1.008 (1.006–1.010), and 1.019 (1.010–1.028), respectively [53]. This is consistent with our findings. In addition, the larger estimate might be associated with exposure patterns. People generally spend more time outdoors in warmer cities, resulting in smaller measurement errors [52, 54].

This study provided detailed estimates of the risk of PM-associated hospital admission for pneumonia through the use of data from a national health insurance program that serves approximately one fifth of the population in China, and offered a unique advantage in identifying potential modifiers of the associations. Our study also has some limitations. First, the use of citywide monitoring measurements as a proxy for population exposure could cause exposure measurement errors, which tend to bias the estimates downward [55]. Second, children (aged <18 years) were not included in this study because UEBMI is for urban working and retired employees, although young children may be more sensitive to air pollution, similarly to the elderly. In addition, only urban employed and retired individuals were included in this study based on the type of insurance enrolled. Due to the differences in sociodemographic characteristics and PM concentrations between rural and urban areas, the generalizability of our findings should be interpreted with caution. Third, both one-stage and two-stage methods could be applied in multicity analysis. In this study, the two-stage method was applied, considering the issue of the computational complexity [56], in line with previous studies [13, 27, 32]. The variances estimated in the first step in the two-stage analysis may not be exactly those in the one-stage method. However, it was reported that one-stage and two-stage methods often generate very similar results [57]. Fourth, in this study, we examined the associations between PM and hospital admissions for pneumonia, not including emergency room visits. As hospital admissions generally include the more severe cases, future investigations are needed to assess the effects of PM on other morbidity outcomes (e.g., emergency room visits). Fifth, as in other environmental health studies using a large administrative health database [13, 27, 30, 58], data on several patient-level variables, such as medical history and smoking status, were not available, limiting the ability to identify potentially susceptible populations. More detailed individual information would be needed to assess the modifying effects of individual characteristics on the associations between PM pollution and pneumonia.

In summary, we found that short-term elevations in PM were associated with increased pneumonia-related hospital admissions in Chinese adults. Our findings support the rationale to further limit PM concentrations in China.

Supporting information

S1 Table. Summary statistics on city-level characteristics in 184 Chinese cities.

https://doi.org/10.1371/journal.pmed.1003010.s001

(DOCX)

S2 Table. City-specific percentage increase with 95% CI in daily hospital admissions for pneumonia associated with a 10 ΞΌg/m3 increase in PM2.5 and PM10 concentrations (lag 0–2) in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.s002

(DOCX)

S3 Table. Regional-average percentage increase with 95% CI in daily hospital admissions for pneumonia associated with a 10 ΞΌg/m3 increase in PM2.5 and PM10 concentrations (lag 0–2) in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.s003

(DOCX)

S4 Table. National-average percentage increase with 95% CI in daily hospital admissions for pneumonia associated with a 10 ΞΌg/m3 increase in PM2.5 and PM10 concentrations (lag 0–2) in two-pollutant models in 184 Chinese cities, 2014–2017.

https://doi.org/10.1371/journal.pmed.1003010.s004

(DOCX)

References

  1. 1. Mandell LA, Wunderink RG, Anzueto A, Bartlett JG, Campbell GD, Dean NC, et al. Infectious Diseases Society of America/American Thoracic Society consensus guidelines on the management of community-acquired pneumonia in adults. Clin Infect Dis. 2007;1(44):511159.
  2. 2. Morimoto K, Suzuki M, Ishifuji T, Yaegashi M, Asoh N, Hamashige N, et al. The burden and etiology of community-onset pneumonia in the aging Japanese population: a multicenter prospective study. PLoS ONE. 2015;10(3):e0122247. pmid:25822890
  3. 3. Guan X, Silk BJ, Li W, Fleischauer AT, Xing X, Jiang X, et al. Pneumonia incidence and mortality in Mainland China: systematic review of Chinese and English literature, 1985–2008. PLoS ONE. 2010;5(7):e11721. pmid:20668535
  4. 4. Fine MJ, Smith MA, Carson CA, Mutha SS, Sankey SS, Weissfeld LA, et al. Prognosis and outcomes of patients with community-acquired pneumonia. A meta-analysis. JAMA. 1996;275(2):134–41. pmid:8531309
  5. 5. Corrales-Medina VF, Alvarez KN, Weissfeld LA, Angus DC, Chirinos JA, Chang C-CH, et al. Association between hospitalization for pneumonia and subsequent risk of cardiovascular disease. JAMA. 2015;313(3):264–74. pmid:25602997
  6. 6. Welte T, Torres A, Nathwani D. Clinical and economic burden of community-acquired pneumonia among adults in Europe. Thorax. 2012;67(1):71–9. pmid:20729232
  7. 7. Wroe PC, Finkelstein JA, Ray GT, Linder JA, Johnson KM, Rifas-Shiman S, et al. Aging population and future burden of pneumococcal pneumonia in the United States. J Infect Dis. 2012;205(10):1589–92. pmid:22448012
  8. 8. Brunekreef B, Holgate ST. Air pollution and health. Lancet. 2002;360(9341):1233–42. pmid:12401268
  9. 9. Anderson HR, Atkinson RW, Bremner SA, Marston L. Particulate air pollution and hospital admissions for cardiorespiratory diseases: are the elderly at greater risk? Eur Respir J Suppl. 2003;40:39s–46s. pmid:12762573
  10. 10. Gilmour PS, Rahman I, Donaldson K, MacNee W. Histone acetylation regulates epithelial IL-8 release mediated by oxidative stress from environmental particles. Am J Physiol Lung Cell Mol Physiol. 2003;284(3):22.
  11. 11. Gurgueira SA, Lawrence J, Coull B, Murthy GG, Gonzalez-Flecha B. Rapid increases in the steady-state concentration of reactive oxygen species in the lungs and heart after particulate air pollution inhalation. Environ Health Perspect. 2002;110(8):749–55. pmid:12153754
  12. 12. Hou W, Xu X, Lei Y, Cao J, Zhang Y, Chen L, et al. The role of the PM2.5-associated metals in pathogenesis of child Mycoplasma Pneumoniae infections: a systematic review. Environ Sci Pollut Res Int. 2016;23(11):10604–14. pmid:27040534
  13. 13. Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, et al. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA. 2006;295(10):1127–34. pmid:16522832
  14. 14. Medina-Ramon M, Zanobetti A, Schwartz J. The effect of ozone and PM10 on hospital admissions for pneumonia and chronic obstructive pulmonary disease: a national multicity study. Am J Epidemiol. 2006;163(6):579–88. pmid:16443803
  15. 15. Strosnider HM, Chang HH, Darrow LA, Liu Y, Vaidyanathan A, Strickland MJ. Age-Specific Associations of Ozone and Fine Particulate Matter with Respiratory Emergency Department Visits in the United States. Am J Respir Crit Care Med. 2019;199(7):882–90. pmid:30277796
  16. 16. Horne BD, Joy EA, Hofmann MG, Gesteland PH, Cannon JB, Lefler JS, et al. Short-Term Elevation of Fine Particulate Matter Air Pollution and Acute Lower Respiratory Infection. Am J Respir Crit Care Med. 2018;198(6):759–66. pmid:29652174
  17. 17. Host S, Larrieu S, Pascal L, Blanchard M, Declercq C, Fabre P, et al. Short-term associations between fine and coarse particles and hospital admissions for cardiorespiratory diseases in six French cities. Occup Environ Med. 2008;65(8):544–51. pmid:18056749
  18. 18. Tian L, Qiu H, Pun VC, Lin H, Ge E, Chan JC, et al. Ambient carbon monoxide associated with reduced risk of hospital admissions for respiratory tract infections. Am J Respir Crit Care Med. 2013;188(10):1240–5. pmid:23944864
  19. 19. Qiu H, Tian LW, Pun VC, Ho KF, Wong TW, Yu IT. Coarse particulate matter associated with increased risk of emergency hospital admissions for pneumonia in Hong Kong. Thorax. 2014;69(11):1027–33. pmid:25164925
  20. 20. Brauer M, Freedman G, Frostad J, van Donkelaar A, Martin RV, Dentener F, et al. Ambient Air Pollution Exposure Estimation for the Global Burden of Disease 2013. Environ Sci Technol. 2016;50(1):79–88. pmid:26595236
  21. 21. Xu Q, Li X, Wang S, Wang C, Huang F, Gao Q, et al. Fine Particulate Air Pollution and Hospital Emergency Room Visits for Respiratory Disease in Urban Areas in Beijing, China, in 2013. PLoS ONE. 2016;11(4):e0153099. pmid:27054582
  22. 22. Zhang Z, Hong Y, Liu N. Association of ambient Particulate matter 2.5 with intensive care unit admission due to pneumonia: a distributed lag non-linear model. Sci Rep. 2017;7(1):8679. pmid:28819316
  23. 23. Yip WC, Hsiao WC, Chen W, Hu S, Ma J, Maynard A. Early appraisal of China's huge and complex health-care reforms. Lancet. 2012;379(9818):833–42. pmid:22386036
  24. 24. Tian Y, Liu H, Zhao Z, Xiang X, Li M, Juan J, et al. Association between ambient air pollution and daily hospital admissions for ischemic stroke: A nationwide time series analysis. PLoS Med. 2018;15(10):e1002668. pmid:30286080
  25. 25. Tian Y, Liu H, Si Y, Cao Y, Song J, Li M, et al. Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: A national time-series study. PLoS Med. 2019;16(1):e1002738. pmid:30689640
  26. 26. Zhao B, Su Y, He S, Zhong M, Cui G. Evolution and comparative assessment of ambient air quality standards in China. J Integr Environ Sci. 2016;13(2–4):85–102.
  27. 27. Chen R, Yin P, Meng X, Liu C, Wang L, Xu X, et al. Fine Particulate Air Pollution and Daily Mortality. A Nationwide Analysis in 272 Chinese Cities. Am J Respir Crit Care Med. 2017;196(1):73–81. pmid:28248546
  28. 28. Yin P, Chen R, Wang L, Meng X, Liu C, Niu Y, et al. Ambient Ozone Pollution and Daily Mortality: A Nationwide Study in 272 Chinese Cities. Environ Health Perspect. 2017;125(11):117006. pmid:29212061
  29. 29. Chen R, Kan H, Chen B, Huang W, Bai Z, Song G, et al. Association of particulate air pollution with daily mortality: the China Air Pollution and Health Effects Study. Am J Epidemiol. 2012;175(11):1173–81. pmid:22510278
  30. 30. Yin P, He G, Fan M, Chiu KY, Liu C, Xue A, et al. Particulate air pollution and mortality in 38 of China's largest cities: time series analysis. BMJ. 2017;356:j667. pmid:28292780
  31. 31. Peng RD, Chang HH, Bell ML, McDermott A, Zeger SL, Samet JM, et al. Coarse particulate matter air pollution and hospital admissions for cardiovascular and respiratory diseases among Medicare patients. JAMA. 2008;299(18):2172–9. pmid:18477784
  32. 32. Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. Ambient Particulate Air Pollution and Daily Mortality in 652 Cities. N Engl J Med. 2019;381(8):705–15. pmid:31433918
  33. 33. Chen R, Yin P, Meng X, Wang L, Liu C, Niu Y, et al. Associations between Coarse Particulate Matter Air Pollution and Cause-Specific Mortality: A Nationwide Analysis in 272 Chinese Cities. Environ Health Perspect. 2019;127(1):17008. pmid:30702928
  34. 34. Samoli E, Analitis A, Touloumi G, Schwartz J, Anderson HR, Sunyer J, et al. Estimating the exposure-response relationships between particulate matter and mortality within the APHEA multicity project. Environ Health Perspect. 2005;113(1):88–95. pmid:15626653
  35. 35. Altman DG, Bland JM. Interaction revisited: the difference between two estimates. BMJ. 2003;326(7382):219. pmid:12543843
  36. 36. Knapp G, Hartung J. Improved tests for a random effects meta-regression with a single covariate. Stat Med. 2003;22(17):2693–710. pmid:12939780
  37. 37. Thompson SG, Sharp SJ. Explaining heterogeneity in meta-analysis: a comparison of methods. Stat Med. 1999;18(20):2693–708. pmid:10521860
  38. 38. Huang Y, Shen H, Chen H, Wang R, Zhang Y, Su S, et al. Quantification of global primary emissions of PM2.5, PM10, and TSP from combustion and industrial process sources. Environ Sci Technol. 2014;48(23):13834–43. pmid:25347079
  39. 39. Duan Z, Han X, Bai Z, Yuan Y. Fine particulate air pollution and hospitalization for pneumonia: a case-crossover study in Shijiazhuang, China. Air Qual Atmos Hlth. 2016;9(7):723–33.
  40. 40. Lippmann M, Schlesinger RB. Toxicological bases for the setting of health-related air pollution standards. Annu Rev Public Health. 2000;21:309–33. pmid:10884956
  41. 41. Wilson WE, Suh HH. Fine particles and coarse particles: concentration relationships relevant to epidemiologic studies. J Air Waste Manag Assoc. 1997;47(12):1238–49. pmid:9448515
  42. 42. Pope CA 3rd, Dockery DW. Health effects of fine particulate air pollution: lines that connect. J Air Waste Manag Assoc. 2006;56(6):709–42. pmid:16805397
  43. 43. Becker S, Soukup JM. Exposure to urban air particulates alters the macrophage-mediated inflammatory response to respiratory viral infection. J Toxicol Environ Health A. 1999;57(7):445–57. pmid:10494914
  44. 44. Tsai SS, Yang CY. Fine particulate air pollution and hospital admissions for pneumonia in a subtropical city: Taipei, Taiwan. J Toxicol Environ Health A. 2014;77(4):192–201. pmid:24555678
  45. 45. Liu H, Tian Y, Xiang X, Juan J, Song J, Cao Y, et al. Ambient Particulate Matter Concentrations and Hospital Admissions in 26 of China's Largest Cities: A Case-Crossover Study. Epidemiology. 2018;29(5):649–57. pmid:29870428
  46. 46. Seaton A, Dennekamp M. Hypothesis: ill health associated with low concentrations of nitrogen dioxideβ€”an effect of ultrafine particles?: Thorax. 2003;58(12):1012–5. pmid:14645960
  47. 47. Requia WJ, Koutrakis P, Roig HL, Adams MD. Spatiotemporal analysis of traffic emissions in over 5000 municipal districts in Brazil. J Air Waste Manag Assoc. 2016;66(12):1284–93. pmid:27623986
  48. 48. Gao J, Wang K, Wang Y, Liu S, Zhu C, Hao J, et al. Temporal-spatial characteristics and source apportionment of PM2.5 as well as its associated chemical species in the Beijing-Tianjin-Hebei region of China. Environ Pollut. 2018;233:714–24. pmid:29126093
  49. 49. Janssens JP. Aging of the respiratory system: impact on pulmonary function tests and adaptation to exertion. Clin Chest Med. 2005;26(3):469–84. pmid:16140139
  50. 50. Basu R, Pearson D, Malig B, Broadwin R, Green R. The effect of high ambient temperature on emergency room visits. Epidemiology. 2012;23(6):813–20. pmid:23007039
  51. 51. Onozuka D, Hashizume M, Hagihara A. Impact of weather factors on Mycoplasma pneumoniae pneumonia. Thorax. 2009;64(6):507–11. pmid:19318345
  52. 52. Macdonald RW, Harner T, Fyfe J. Recent climate change in the Arctic and its impact on contaminant pathways and interpretation of temporal trend data. Sci Total Environ. 2005;342(1–3):5–86. pmid:15866268
  53. 53. Chen F, Fan Z, Qiao Z, Cui Y, Zhang M, Zhao X, et al. Does temperature modify the effect of PM10 on mortality? A systematic review and meta-analysis. Environ Pollut. 2017;224:326–35. pmid:28215581
  54. 54. Chen R, Peng RD, 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;451:259–65.
  55. 55. Goldman GT, Mulholland JA, Russell AG, Strickland MJ, Klein M, Waller LA, et al. Impact of exposure measurement error in air pollution epidemiology: effect of error type in time-series studies. Environ Health. 2011;10(61):10–61.
  56. 56. Thompson S, Kaptoge S, White I, Wood A, Perry P, Danesh J. Statistical methods for the time-to-event analysis of individual participant data from multiple epidemiological studies. Int J Epidemiol. 2010;39(5):1345–59. pmid:20439481
  57. 57. Turner RM, Omar RZ, Yang M, Goldstein H, Thompson SG. A multilevel model framework for meta-analysis of clinical trials with binary outcomes. Stat Med. 2000;19(24):3417–32. pmid:11122505
  58. 58. Di Q, Dai L, Wang Y, Zanobetti A, Choirat C, Schwartz JD, et al. Association of Short-term Exposure to Air Pollution With Mortality in Older Adults. JAMA. 2017;318(24):2446–56. pmid:29279932