- 1Jiading District Center for Disease Control and Prevention, Shanghai, China
- 2Fudan University School of Public Health, Shanghai, China
- 3Key Laboratory of Public Health Safety, Fudan University, Ministry of Education, Shanghai, China
- 4Fudan University Center for Tropical Disease Research, Shanghai, China
- 5Faculty of Medicine, School of Epidemiology and Public Health, University of Ottawa, Ottawa, ON, Canada
Background: Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of death in China. Although numerous studies have been conducted to determine the risk factors for COPD mortality such as ambient air pollution, the results are not fully consistent.
Methods: This study included mortality analysis and a case-control design by using the data extracted from the Mortality Registration System in Jiading District, Shanghai. Traditional logistic regression, geographically weighted logistic regression (GWLR), and spatial scan statistical analysis were performed to explore the geographic variation of COPD mortality and the possible influencing factors.
Results: Traditional logistic regression showed that extreme lower temperature in the month prior to death, shorter distance to highway, lower GDP level were associated with increased COPD mortality. GWRL model further demonstrated obvious geographical discrepancies for the above associations. We additionally identified a significant cluster of low COPD mortality (OR = 0.36, P = 0.002) in the southwest region of Jiading District with a radius of 3.55 km by using the Bernoulli model. The geographical variation in age-standardized mortality rate for COPD in Jiading District was explained to a certain degree by these factors.
Conclusion: The risk of COPD mortality in Jiading District showed obvious geographical variation, which were partially explained by the geographical variations in effects of the extreme low temperature in the month prior to death, residential proximity to highway, and GDP level.
Introduction
The 2017 Global Burden of Science Study (GBD) found that the total number of deaths due to chronic respiratory diseases increased by 18.0%, from 3.32 (95% CI: 3.01 to 3.43) million in 1990 to 3.91 (95% CI: 3.79 to 4.04) million in 2017 (1). Chronic obstructive pulmonary disease (COPD), the most common type of chronic respiratory disease, is characterized by airflow obstruction that can progress to pulmonary heart disease and respiratory failure. The China Pulmonary Health (CPH) Study, which was carried out among 50,991 adults aged 20 years or older during 2012-2015, reported an overall prevalence of COPD being 8.6% (11.9% in men and 5.4% in women) (2). COPD is currently the fourth leading cause of death in China (3), making it a one of most pressing health issues in the country.
Smoking is a major cause of COPD (4) and ambient air pollutants play an important part in COPD exacerbation (5, 6), and both are related to increased mortality and decreased quality of life. Previous studies have shown high heterogeneity in geographic variation in COPD mortality, which may be associated with personal characteristics, air pollution and other factors such as temperature (7) or may be due to study limitations. Many studies used ecological time-series method to determine the effect of air pollutants on COPD (7) and low spatiotemporal resolution of measurements of ambient air pollution obtained from fixed monitoring stations as a surrogate for individual exposure, and ecological fallacy is a concern. A study found traffic-related pollution was affected by many factors including dynamic emissions, measurement distance and height from the roadway, street canyon, obstructions along the way, and meteorological conditions such as wind direction and speed (8). When the exposure concentration was highly spatially and geographically heterogeneous, mis-estimation of exposures would be introduced into analysis (9). In recent years, multi-pollutant indicators have been adopted to investigate the association with health outcomes to better evaluate real-world exposure to pollutants such as residential distance to the closest roads (10), factory density (11), and traffic density (12). Residential proximity to road has been found to be a good proxy for traffic-related air pollutant concentrations (13).
The current study aimed to determine the geographical variation in COPD mortality risk, and to explore the geographical variations in the effects of ambient air pollution indicators and other factors on COPD mortality. We used geographically weighted regression model (GWR) to explore the spatial variation in the risk of COPD death and relevant driving factors at a local level. GWR is an extension of weighted regression by establishing the local regression equation at each point in the space scope, where the weights are defined by the geographical distance between objectives (14), and has been to study different diseases (15, 16).
Materials and Methods
Study Site
This study was conducted in Jiading District, located in the northwest of Shanghai. The district covers an area of 464.2 km2 and 12 towns with a registered population of 658,200. Jiading District is situated in the northern edge of northern subtropical region with abundant precipitation, warm and humid climate, and adequate sunshine. The annual mean temperature was 17.6°C in 2018, while the lowest temperature reached −6.4°C in the same year (17). Since 2010, the roads system in Jiading District has been gradually developing into a network, with a total length of 1,407 km. Meanwhile, the main traffic network of the district currently includes five vertical and six horizontal lines.
Study Design
We performed a population-based case-control study by using the data extracted from the Mortality Registration System in Jiading District, Shanghai. The underlying causes of death coded by the International Statistical Classification of Diseases, 10th revision (ICD-10), were reclassified following the method of GBD study (3). Cases were the deaths due to COPD between 2011 and 2018 with complete information on age, sex, education, occupation, marriage, date of death, residential address. Two groups of controls were selected from deaths due to cerebrovascular disease and due to tumor. Frequency matching was used to match the two control groups to the case group with age (5 years or less), sex and year of death. The study included a total of 3,438 deaths (1,146 COPD cases, 1,146 cerebrovascular controls, and 1,146 tumor controls).
Data Collection
Population Data
Population data from 2011 to 2018 were taken from the Public Security Bureau in Jiading District. Age-standarized mortality rates (ASMR) were calculated according to sex, town, and 5-year time period.
Vector Data
We collected the following vector data. The map of Jiading District and roads data in 2010 were purchased from the Shanghai Digital Bitmap Information Technology Co., LTD. The roads in Jiading District were divided into main roads (national, provincial, and some county roads with a traffic volume of 10,000 counts/day) and highways (Chinese Expressways). The geographical addresses of factories and public hospitals were obtained from the Jiading District Center for Disease Control and Prevention (CDC), and then we screened out factories that produce air pollutants according to the detection report of factory and technological process. Finally, the geographical addresses of factories and public hospitals were converted into latitude and longitude by using Geocoding software (https://www.esri.com/en-us/arcgis/products/geocoding) and then transformed into vector point data in ArcMap 10.5 (Environmental Systems Research Institute, Inc, Redlands, CA).
Raster Data
Gross domestic product (GDP), population density, and extreme low temperature in the month prior to death were stored as raster data. China's GDP and population per kilometer grid in 2010 were downloaded from the Global Change Research Data Publishing & Repository, with a spatial resolution of 1 km. Previous studies have found that extreme low temperature can cause more serious disease burden of COPD than high temperature (18), and the cold effect can be delayed for up to a month (19). Therefore, the average extreme low temperature in the month prior to death was used to explore the impact of temperature on COPD death. The extreme low temperature in China by day from December 1, 2010 to November 31, 2018 was downloaded from the Chinese Meteorological Data Network (http://data.cma.cn). Meanwhile, raster data of GDP, population and temperature in Jiading District were extracted by the mask extraction tool in ArcMap10.5.
Data Processing
Firstly, we used the nearest neighbor analysis tool in ArcMap10.5 to calculate distance. The distance between residential address and the closest road was calculated to evaluate individual traffic-related air pollution, the distance to the closest factory was used to represent individual industrial-related air pollution, and the distance to the closest public hospital was used as an indicator of individual health care accessibility. Then, we also calculated the factory density according to geographical distribution of factories, in which the deceased once lived, by using a kernel density analysis tool in ArcMap 10.5. Finally, we superimposed the residential address with raster data and evaluated individual exposure level by using the Extract Values to Points tool in ArcMap 10.5. In this study, GDP was used as an indicator of the economic status of the dead and population per kilometer was adopted to describe population distribution in Jiading District. Besides, temperature represented the average extreme low temperature in the month prior to death in this study.
Exposure Classification
We classified socioeconomic factors including marriage (married or others), education (senior high school or above, junior middle school, primary school or below) and occupation (industry-related, agriculture-related, others). Extreme low temperature in the month prior to death was classified into four levels (<5, 5–15, 16–20, >20°C) with >20°C as the reference group. We classified the following distances from the residential address to the closest highways (<200, 200–350, 351–650, >650 m), to the closest major road (<100, 100–300, 301–500, >500 m), to the closest public hospital (<200, 200–500, 501–1,000, 1,001–2,000, >2,000), and to the closest factory (<50, 50–200, 201–500,>500 m), with the most distant group as the reference group for each. We classified residential factory density, GDP, and population density into quartiles (Q1, Q2, Q3, Q4), with Q4 served as the reference group.
Analytic Methods
Age-Standardized Mortality
The average ASMR for COPD between 2011 and 2018 was calculated using the Chinese population in 2010 as standard population, and visualized in ArcMap10.5.
Non-spatial Logistic Regression
For the case-control study, we firstly used simple logistical regression to analyze the association between COPD and each risk factor by using SPSS (SPSS Inc., USA), and significant factors were then included in multivariable logistical regression. Considering the seasonality in COPD mortality, we also divided the death data into cold season (from September to February) and warm season (from March to August), and then we conducted non-spatial logistical regressions in the same method as above-mentioned, respectively.
Geographically Weighted Logistic Regression (GWLR)
A GWLR model was used to analyze the associations between COPD mortality and all the risk factors using the data from the case-control study. The relationship may be described as:
where yi and xij, denote the cause of death and explanatory variables for individual i with location coordinates (ui, vi), and β0 (ui, vi) and βj (ui, vi) are the location-specific intercept and coefficients.
In this GWLR model, we combined two control groups. Adaptive bi-square was used as the kernel type, golden section search as the bandwidth selection method and AICc as criteria. Variables included in the model were extreme low temperature in the month prior to death, distance to the closest highway, and GDP. Data analysis was conducted by using GWR4.0 software. ArcMap10.5 software was used to map the results of the GWLR model.
Spatial Scan Statistical Analysis
The spatial scan statistical analysis was widely applied to detect clustering of cases by using cylindrical scanning windows with dynamic changes in size and position (20). The expected and observed number of cases were used to generate the log likelihood ratio (LLR). Then, the scanning window with the largest LLR value was selected as the high-clustering window. The spatial and temporal information involved in the window was collected and determined. Finally, the statistical significance of those clusters was tested by using the Monte Carlo testing method. In this study, data from 1,146 COPD cases and 2,292 non-COPD controls fit the Bernoulli model with a 6-month interval. The analysis was carried out by SaTScan9.6 software (https://www.satscan.org/) and then visualized the significant spatial cluster in ArcMap 10.5.
Results
The Average ASMR
Figure 1 showed that the average ASMR due to COPD was the lowest in Jiading town. We observed a lower risk of COPD mortality from the northeast corner of Jiading District across the southwest corner of this district. Besides, the ASMR in the north, northwest, and southeast of Jiading District appeared to be higher than other areas.
Non-spatial Logistic Regression
Case and control groups are frequently matched by sex (male: 69.1%, female: 30.9% and age (83.8 ± 7.7 for COPD deaths, 83.4 ± 7.4 for cerebrovascular disease deaths, and 82.8 ± 7.1 for tumor deaths, respectively), and there were no significant differences in sex and age distributions between cases and controls.
Table 1 showed the results from logistic regression models based on data from 1,146 COPD cases and 1,146 cerebrovascular controls. In the adjusted model, industry- related occupation (vs. others, aOR = 1.414, 95% CI: 1.079–1.853), the extreme lower temperature (<5°C vs. >20°C, aOR = 1.496, 95% CI: 1.169–1.915), shorter distance to the closest highway (<200 m vs. >650 m, aOR = 1.819, 95% CI: 1.134–2.919), and lower GDP level (Q1 vs. Q4, aOR = 1.505, 95% CI: 1.019–2.225) were significantly associated with increased risk of COPD death. In addition, population density Q3 was significantly and negatively associated with the risk of COPD death (aOR = 0.449, 95% CI: 0.322–0.627) compared to Q4.
Table 1. The association of COPD mortality and risk factors from non-spatial logistical regressions using cerebrovascular disease deaths as controls.
When we used tumor deaths as controls, the results were similar (Table 2). In the multivariable logistic regression model, the extreme lower temperature (<5°C vs. >20°C, aOR = 2.052, 95% CI: 1.611–2.613), shorter distance to the closest highway (200–350 m vs. >650 m, aOR = 2.119, 95% CI: 1.521–2.951), shorter distance to the closest factory (200–500 m vs. >500 m, aOR = 1.595, 95% CI: 1.218–2.090), and lower GDP level (Q1 vs. Q4, aOR = 1.795, 95% CI: 1.179–2.733) were significantly associated with increased risk of COPD death. Besides, married people (aOR = 0.754, 95% CI: 0.625–0.909) showed lower risk of COPD death (Table 2). Furthermore, although we divided the deaths into cold season (from September to February) and warm season (from March to August), we still found that extreme lower temperature, shorter distance to the closest highway, lower GDP level were associated with higher risk of COPD mortality (Supplementary Table 1).
Table 2. The association of COPD mortality and risk factors from non-spatial logistical regressions using tumor deaths as controls.
GWLR Model
In this study, the best bandwidth size was 654. The value of AICC in the global regression was 3172.3, while the value was 2927.6 for GWLR, which indicated that the GWLR model fitted the data better than the regular regression model. Geographical distribution of risk factors associated with COPD mortality from GWLR model was presented in Figure 2. The OR values from GWR indicated that extreme lower temperature, shorter distance to highway and lower GDP increased the risk of COPD death, which were consistent with results from non-spatial logistic regression. However, GWR added that the spatial distribution of significant clusters for each risk factors were obviously different. Extreme low temperature (per 1°C reduce) increased COPD mortality (P < 0.05) for 40.2% of study individuals and the significant cluster areas marked in red (Figure 2A). Residential proximity to highway (per 1 km reduce) was positively (P < 0.05) associated with COPD mortality for 44.2% individuals (Figure 2B). GDP (per 1 unit reduce) also increased COPD mortality for 35.2% study individuals, marked in red (Figure 2C).
Figure 2. Geographical distribution of risk factors associated with COPD death from GWLR model: (A) (temperature), (B) (distance to the closest highway), (C) (GDP).
Spatial Scan Statistical Analysis
Figure 3 presented the results of spatial scan statistical analysis by using Bernoulli model with 6-month interval, the triangles in red denoted the residential location of COPD cases while points in green denoted the residential location of controls. There was a significant cluster of lower risk of COPD mortality (OR = 0.36, P = 0.002) in the southwest of Jiading District, the cluster spanned the period from 2011/1/1 to 2014/12/31 with a radius of 3.55 km.
Figure 3. Distribution of COPD case and controls in Jiading District with the location of significant spatial cluster, 2011–2018.
Discussion
According to the geographic distribution of ASMR in Jiading district, Shanghai, lower ASMR appeared from the northeast to the southwest of Jiading District, while higher ASMR distributed in the north, northwest, and southeast of Jiading District, which highlighted the geographical variation of COPD mortality in Jiading District. The results from non-spatial logistic regression indicated extreme lower temperature, shorter distance to highway, lower GDP level might be associated with increased risk of COPD mortality. Besides, the results from GWLR model not only further verified the association between COPD mortality and the above risk factors but also added that the spatial distribution of significant clusters for each risk factors were obviously different. Finally, the Bernoulli model showed a significant cluster of reduced COPD mortality, which was consistent with the map of ASMR and could be interpreted as few traffic routes (according to road distribution in Jiading District), relatively higher temperature and higher GDP level in this area.
Previous study conducted in Ningbo found a reverse J-shape association between COPD mortality and average air temperature, both extreme low temperature (RR = 2.767, 95% CI: 1.950–3.928) and high temperature (RR = 1.197, 95% CI: 1.021–1.404) over lags of 0–14 days could increase COPD mortality (21). Another study also revealed that low temperature caused more serious effect on COPD mortality than high temperature by using the distributed lag non-linear model (18). A time-series panel study was carried out by measuring peak expiratory flow (PEF) and forced expiratory volume in one second (FEV1) in 28 male COPD patients repeatedly to estimate the cumulative effect of temperature on pulmonary function, and the results showed robust evidence that both low and high temperatures were significantly related to the decline of pulmonary function, PEF in particular (22).
Traffic becomes one of the major contributors to ambient air pollution in urban environments (23). A growing number of studies have used traffic as a measure of air pollution exposure to evaluate its negative effects on lung cancer (9), childhood asthma (24), and cardiovascular diseases (25). In this study, residential proximity to highways was a risk factor for COPD mortality while residential proximity to roads showed no significant association with COPD mortality. Previous studies showed consistent results for the relationship between COPD mortality and traffic-related pollution. A cross-sectional study conducted in Jiading District from 2001 to 2010 found that arterial road density was significantly associated with the mortality of COPD (26). Proximity to roads (0–200 m) and road density are predictors for COPD morbidity in Kandy (27). A cohort study performed by Cakmak et al. in Canada found that exposure to higher road density and proximity to major traffic roads were related to increased COPD mortality (28). However, they selected the length of local roads within a 200 m radius of postal code centroid rather than residential distance to the closest road as indicator, which might explain part of the geographical heterogeneity observed since different traffic parameters led to various sensitivity in estimates of exposure (11). The discrepancy between the results of residential proximity to highway and road might be due to a wide emission range of pollutants on expressways and a declined concentration with increasing distance (29).
In this study, we found people engaged in agriculture or industrial-related occupations, unmarried, and with lower education level showed higher risk of COPD death, which were consistent with previous studies. A national cross-sectional study of COPD showed that middle or high school education (OR = 0.76, 95% CI: 0.64–0.90) and college or higher education (OR = 0.47, 95% CI: 0.33–0.66) had lower risk of COPD (7). Previous studies have found that exposure to pollutants in the occupational environment would exacerbate airflow limitations in COPD patients and increase the risk of COPD death (30, 31). In addition, we used GDP as an indicator for social-economic status, and the results showed that lower GDP level was risk factor for COPD mortality.
The results in GWLR model further verified the relationship between COPD mortality and extreme low temperature in the month prior to death, residential proximity to highway, and GDP. We also found significant geographical variations in the effects of these factors. For example, the significant clusters of residential proximity to highway mainly located from the southeast across the northwest of Jiading District, adjacent to the highway. COPD mortality in northeast of Jiading District was less affected. In addition, a previous study showed that wind direction greatly affected the concentration of ambient pollutants (32). Due to the influence of southeast monsoon, the dominant wind direction of Jiading District is from the southeast in summer but northwest in winter, which could partially explain the distributions of the clusters. In addition, we selected a representative map of extreme low temperature in January 1, 2011 to describe the distribution of temperature in Jiading District. Relatively lower temperature and lower GDP might partially explain the higher ASMR in northern and northwest of Jiading District. Because of urban heat island effect (33), rapid social-economic development, and the policy that avoid building highways in downtown area, the center of Jiading District (Jiading town) had relatively higher temperature, higher GDP level, and longer distance to highway, which might attributed to the lowest ASMR of COPD. The Bernoulli model showed a significant cluster of reduced COPD mortality, which partially overlapped with the significant clusters of the extreme low temperature in the month prior to death, proximity to highway and GDP in the GWLR model, which indicated few traffic routes (according to road distribution in Jiading District), relatively higher temperature and higher GDP level in this area may explain the reduced COPD mortality cluster.
Our study had several strengths and limitations. In this study, conventional data and various land use variables at fine spatial resolutions were combined to explore the factors associated with COPD death and spatial distribution difference in a relatively novel method. Besides, comprehensive indicators were used to better reflect individual exposure to pollutants such as distance to the closest highway, factory density. However, our study also had several limitations, Firstly, we used China's GDP and population per kilometer grids in 2010 to represent the social-economic status and population density of the whole research period, and used roads data in 2010 to evaluate traffic-related pollution, which would lead to inaccurate assessment of present exposure. Secondly, due to the lack of individual information collection, the potential risk factors included in our study were incomplete, particularly the lack of data on smoking (34) and body mass index (BMI) (35), which might confound the associations in the current study. However, a study conducted in China found 2 years after being diagnosed, patients with COPD had smoking cessation rates of 62.8% (36). In addition, smoking is a predominant behavior for men, frequency matching by sex might reduce the confounding effects of smoking.
Conclusion
The risk of COPD mortality in Jiading District showed an obvious geographical variation, which might partially be explained by geographical variations in the effects of the extreme low temperature in the month prior to death, residential proximity to highway, and GDP level. These factors should be taken into consideration when targeted interventions are developed at the local level to reduce COPD mortality risk.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
Written informed consent was obtained from the individual(s), and minor(s)' legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author Contributions
YZho is the project leader and contributed to all aspects of this work. QP and NZ are joint first authors, who contributed equally to data assembly, analysis, and drafted the manuscript of the present paper. All authors contributed intellectually to this manuscript and have approved this final version.
Funding
This work was supported by the Jiading District, Shanghai, Epidemiology of Chronic Diseases Foundation (grant number GW01).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We gratefully acknowledge the participants for their continuous and enthusiastic participation in the investigation. Also, we appreciate doctors and nurses involved in this study for their technical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2021.627312/full#supplementary-material
References
1. Li X, Cao X, Guo M, Xie M, Liu X. Trends and risk factors of mortality and disability adjusted life years for chronic respiratory diseases from 1990 to 2017: systematic analysis for the Global Burden of Disease Study 2017. BMJ. (2020) 368:m234. doi: 10.1136/bmj.m234
2. Wang C, Xu J, Yang L, Xu Y, Zhang X, Bai C, et al. Prevalence and risk factors of chronic obstructive pulmonary disease in China (the China Pulmonary Health [CPH] study): a national cross-sectional study. Lancet. (2018) 391:1706–17. doi: 10.1016/S0140-6736(18)30841-9
3. Zhou M, Wang H, Zeng X, Yin P, Zhu J, Chen W, et al. Mortality, morbidity, and risk factors in China and its provinces, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. (2019) 394:1145–58. doi: 10.1016/S0140-6736(19)30427-1
4. Wheaton AG, Liu Y, Croft JB, VanFrank B, Croxton TL, Punturieri A, et al. Chronic obstructive pulmonary disease and smoking status - United States, 2017. MMWR Morb Mortal Wkly Rep. (2019) 68:533–8. doi: 10.15585/mmwr.mm6824a1
5. Liang L, Cai Y, Barratt B, Lyu B, Chan Q, Hansell AL, et al. Associations between daily air quality and hospitalisations for acute exacerbation of chronic obstructive pulmonary disease in Beijing, 2013-17: an ecological analysis. Lancet Planet Health. (2019) 3:e270–9. doi: 10.1016/S2542-5196(19)30085-3
6. Pfeffer PE, Donaldson GC, Mackay AJ, Wedzicha JA. Increased chronic obstructive pulmonary disease exacerbations of likely viral etiology follow elevated ambient nitrogen oxides. Am J Respir Crit Care Med. (2019) 199:581–91. doi: 10.1164/rccm.201712-2506OC
7. Moore E, Chatzidiakou L, Kuku MO, Jones RL, Smeeth L, Beevers S, et al. Global associations between air pollutants and chronic obstructive pulmonary disease hospitalizations. A systematic review. Ann Am Thorac Soc. (2016) 13:1814–27. doi: 10.1513/AnnalsATS.201601-064OC
8. Wang JM, Jeong CH, Hilker N, Shairsingh KK, Healy RM, Sofowote U, et al. Near-road air pollutant measurements: accounting for inter-site variability using emission factors. Environ Sci Technol. (2018) 52:9495–04. doi: 10.1021/acs.est.8b01914
9. Lee HJ, Chatfield RB, Bell ML. Spatial analysis of concentrations of multiple air pollutants using NASA DISCOVER-AQ aircraft measurements: implications for exposure assessment. Environ Res. (2018) 160:487–98. doi: 10.1016/j.envres.2017.10.017
10. Shao Y, Wang Y, Yu H, Zhang Y, Xiang F, Yang Y, et al. Geographical variation in lung cancer risk associated with road traffics in Jiading District, Shanghai. Sci Total Environ. (2019) 652:729–35. doi: 10.1016/j.scitotenv.2018.10.266
11. Li X, Tan H, Luo M, Wu X, Huang X, Zhou S, et al. Exposure to firework chemicals from production factories in pregnant women and risk of preterm birth occurrence in Liuyang, China. J Toxicol Environ Health. (2018) 81:154–9. doi: 10.1080/15287394.2017.1415579
12. Liu SV, Chen FL, Xue J. Evaluation of traffic density parameters as an indicator of vehicle emission-related near-road air pollution: a case study with NEXUS measurement data on black carbon. Int J Environ Res Public Health. (2017) 14:1581. doi: 10.3390/ijerph14121581
13. Jhun I, Kim J, Cho B, Gold DR, Schwartz J, Coull BA, et al. Synthesis of Harvard Environmental Protection Agency (EPA) Center studies on traffic-related particulate pollution and cardiovascular outcomes in the Greater Boston Area. J Air Waste Manage Assoc. (2019) 69:900–17. doi: 10.1080/10962247.2019.1596994
14. Fotheringham A, Brunsdon C. Geographically Weighted Regression-The Analysis of Spatially Varying Relationships. Chichester: John Wiley & Sons (2002).
15. Park SY, Kwak JM, Seo EW, Lee KS. Spatial analysis of the regional variation of hypertensive disease mortality and its socio-economic correlates in South Korea. Geospat Health. (2016) 11:420. doi: 10.4081/gh.2016.420
16. Mena C, Sepúlveda C, Fuentes E, Ormazábal Y, Palomo I. Spatial analysis for the epidemiological study of cardiovascular diseases: a systematic literature search. Geospat Health. (2018) 13:587. doi: 10.4081/gh.2018.587
17. Shanghai Bureau of Statistics. Shanghai Statistical. Year Book. (2018). Available online at: http://www.stats-sh.gov.cn/html/sjfb/201801/1000339.html (accessed January 31, 2019)
18. Luan G, Yin P, Wang L, Zhou M. Association between ambient temperature and chronic obstructive pulmonary disease: a population-based study of the years of life lost. Int J Environ Health Res. (2019) 29:246–54. doi: 10.1080/09603123.2018.1533533
19. Zhao Y, Huang Z, Wang S, Hu J, Xiao J, Li X, et al. Morbidity burden of respiratory diseases attributable to ambient temperature: a case study in a subtropical city in China. Environ Health. (2019) 18:89. doi: 10.1186/s12940-019-0529-8
21. Lu BB, Gu SH, Wang AH, Ge T, Wang Y, Li XH, et al. [Study on influence of air temperature on daily chronic obstructive pulmonary disease mortality in Ningbo]. Zhonghua Liu Xing Bing Xue Za Zhi. (2017) 38:1528–32. doi: 10.3760/cma.j.issn.0254-6450.2017.11.018
22. Lin Z, Gu Y, Liu C, Song Y, Bai C, Chen R, et al. Effects of ambient temperature on lung function in patients with chronic obstructive pulmonary disease: a time-series panel study. Sci Total Environ. (2018) 619–20:360–5. doi: 10.1016/j.scitotenv.2017.11.035
23. Sinharay R, Gong J, Barratt B, Ohman-Strickland P, Ernst S, Kelly FJ, et al. Respiratory and cardiovascular responses to walking down a traffic-polluted road compared with walking in a traffic-free area in participants aged 60 years and older with chronic lung or heart disease and age-matched healthy controls: a randomised, crossover study. Lancet. (2018) 391:339–49. doi: 10.1016/S0140-6736(17)32643-0
24. Khreis H, Kelly C, Tate J, Parslow R, Lucas K, Nieuwenhuijsen M. Exposure to traffic-related air pollution and risk of development of childhood asthma: a systematic review and meta-analysis. Environ Int. (2017) 100:1–31. doi: 10.1016/j.envint.2016.11.012
25. Cai Y, Hodgson S, Blangiardo M, Gulliver J, Morley D, Fecht D, et al. Road traffic noise, air pollution and incident cardiovascular disease: a joint analysis of the HUNT, EPIC-Oxford and UK Biobank cohorts. Environ Int. (2018) 114:191–201. doi: 10.1016/j.envint.2018.02.048
26. Wang L, Chen R, Sun W, Yang X, Li X. Impact of high-density urban built environment on chronic obstructive pulmonary disease: a case study of Jing'an District, Shanghai. Int J Environ Res Public Health. (2019) 17:252. doi: 10.3390/ijerph17010252
27. Kumarihamy RMK, Tripathi NK. Geostatistical predictive modeling for asthma and chronic obstructive pulmonary disease using socioeconomic and environmental determinants. Environ Monit Assess. (2019) 191(Suppl. 2):366. doi: 10.1007/s10661-019-7417-0
28. Cakmak S, Hebbern C, Vanos J, Crouse DL, Tjepkema M. Exposure to traffic and mortality risk in the 1991-2011 Canadian Census Health and Environment Cohort (CanCHEC). Environ Int. (2019) 124:16–24. doi: 10.1016/j.envint.2018.12.045
29. Wang ZL, Lu QC, He HD, Wang DS, Gao Y, Peng ZR. Investigation of the spatiotemporal variation and influencing factors on fine particulate matter and carbon monoxide concentrations near a road intersection. Front Earth Sci. (2017) 11:63–75. doi: 10.1007/s11707-016-0564-5
30. Dodd KE, Wood J, Mazurek JM. Mortality among persons with both asthma and chronic obstructive pulmonary disease aged ≥25 years, by industry and occupation - United States, 1999-2016. MMWR Morb Mortal Wkly Rep. (2020) 69:670–9. doi: 10.15585/mmwr.mm6922a3
31. Syamlal G, Doney B, Mazurek JM. Chronic obstructive pulmonary disease prevalence among adults who have never smoked, by industry and occupation - United States, 2013-2017. MMWR Morb Mortal Wkly Rep. (2019) 68:303–7. doi: 10.15585/mmwr.mm6813a2
32. Rönkkö T, Kuuluvainen H, Karjalainen P, Keskinen J, Hillamo R, Niemi JV, et al. Traffic is a major source of atmospheric nanocluster aerosol. Proc Natl Acad Sci USA. (2017) 114:7549–54. doi: 10.1073/pnas.1700830114
33. Mohajerani A, Bakaric J, Jeffrey-Bailey T. The urban heat island effect, its causes, and mitigation, with reference to the thermal properties of asphalt concrete. J Environ Manage. (2017) 197:522–38. doi: 10.1016/j.jenvman.2017.03.095
34. van Gemert F, Chavannes N, Kirenga B, Jones R, Williams S, Tsiligianni I, et al. Socio-economic factors, gender and smoking as determinants of COPD in a low-income country of sub-Saharan Africa: FRESH AIR Uganda. NPJ Prim Care Respir Med. (2016) 26:16050. doi: 10.1038/npjpcrm.2016.50
35. Masoud O, Ramsay M, Suh ES, Kaltsakas G, Srivastava S, Pattani H, et al. Long-term adherence to home mechanical ventilation: a 10-year retrospective, single-centre cohort study. J Thorac Dis. (2020) 12(Suppl. 2):S120–8. doi: 10.21037/jtd-cus-2020-003
Keywords: COPD, traffic-related pollutant, temperature, GWLR, mortality
Citation: Peng Q, Zhang N, Yu H, Shao Y, Ji Y, Jin Y, Zhong P, Zhang Y, Jiang H, Li C, Shi Y, Zheng Y, Xiong Y, Wang Z, Jiang F, Chen Y, Jiang Q and Zhou Y (2021) Geographical Variation of COPD Mortality and Related Risk Factors in Jiading District, Shanghai. Front. Public Health 9:627312. doi: 10.3389/fpubh.2021.627312
Received: 09 November 2020; Accepted: 06 January 2021;
Published: 03 February 2021.
Edited by:
Anke Huels, Emory University, United StatesReviewed by:
Howard Chang, Emory University, United StatesXiaoxu Xie, Fujian Medical University, China
Copyright © 2021 Peng, Zhang, Yu, Shao, Ji, Jin, Zhong, Zhang, Jiang, Li, Shi, Zheng, Xiong, Wang, Jiang, Chen, Jiang and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yibiao Zhou, z_yibiao@hotmail.com
†These authors have contributed equally to this work