Original Research

A Global Analysis of Associations between Fine Particle Air Pollution and Cardiovascular Risk Factors: Feasibility Study on Data Linkage

Authors: {'first_name': 'Min', 'last_name': 'Zhao'},{'first_name': 'Gerard', 'last_name': 'Hoek'},{'first_name': 'Maciej', 'last_name': 'Strak'},{'first_name': 'Diederick E.', 'last_name': 'Grobbee'},{'first_name': 'Ian', 'last_name': 'Graham'},{'first_name': 'Kerstin', 'last_name': 'Klipstein-Grobusch'},{'first_name': 'Ilonca', 'last_name': 'Vaartjes'}


Background: This paper presents a feasibility study of data linkage between global air pollution data and clinical medical data to assess the associations of PM2.5 with cardiovascular risk factors.

Methods: Cardiovascular risk factor data were obtained from the SUrvey of Risk Factors (SURF) for coronary heart disease (CHD) patients from 10 countries in Europe, Asia, and the Middle-East. Annual average PM2.5 concentrations were estimated using recent global WHO PM2.5 maps combining satellite and surface monitoring data for the location of the 71 participating centers. Associations of PM2.5 with risk factors were assessed by mixed-effect generalized estimation equation models adjusted by sex, age, exercise, body mass index, and smoking. In the final model there was further adjustment for country.

Results: Linkage between cardiovascular risk factor data and PM2.5 via the postal address of participating hospitals was shown to be feasible, however with several limitations noted.

Eight thousand three hundred and ninety two patients (30% women) were included. Globally, an increase of 10 μg/m3 in PM2.5 was significantly associated with decreased BP and increased glucose. After controlling for country, an increase of 10 μg/m3 in PM2.5 was associated with decreased BP and increased LDL (SBP: –0.45 mmHg [95% CI: –0.85, –0.06]; DBP: –0.47 mmHg [–0.73, –0.20]; LDL: 0.04 mmol/L [0.01, 0.08]). The association with glucose attenuated (0.08 mmol/L [–0.23, 0.16]).

Conclusion: It is feasible to link PM2.5 and cardiovascular risk factors but it is still challenging to interpret these observed associations due to unavailability of potential confounders. After country adjustment, PM2.5 was associated with small increases in LDL and small decreases in BP.


Keywords: air pollutionenvironmental healthcardiovascular diseaserisk factorsfeasibilitydata linkage 
DOI: http://doi.org/10.5334/gh.877
 Accepted on 16 Jul 2020            Submitted on 16 Jul 2020


Cardiovascular disease (CVD) remains one of the leading causes of death worldwide with 18 million deaths in 2016 [1]. Traditionally, evidence based guidelines and daily practice on secondary prevention of CVD have focused on modifiable risk factor management [2, 3]. Several recent epidemiological studies have suggested air pollution could also be associated with CVD risks [4, 5, 6, 7]. The number of studies investigating the association between PM2.5 and modifiable cardiovascular risk factors is scarce [8, 9, 10, 11, 12, 13, 14]. These studies have predominantly been conducted in Western countries with rather low levels of PM2.5 concentrations [8, 9, 10]. In contrast, low- and middle-income countries, for which have limited data on the association of PM2.5 and risk factors, show much higher PM2.5 concentrations [15]. Existing evidence on the role of environmental exposure on cardiovascular risk factors may however not be generalizable to these settings since the chemical composition and characteristics of PM2.5 may differ significantly from those in Western countries [15]. This, together with a rapid increase of CVD prevalence in low- and middle-income countries, stresses the importance of a better understanding of global associations of PM2.5 with cardiovascular risk factors.

Conducting targeted studies on the association between PM2.5 and cardiovascular risk factors on a global scale is challenging. Data linkage is an efficient and cost-effective method to maximize the use of existing data for more health related research questions [16, 17]. Current study aims to assess the feasibility of linking global air pollution data with the cardiovascular risk factors data collected from an international audit to establish the technical and scientific possibilities on data linkage. This study also aims to investigate the potential association between PM2.5 and cardiovascular risk factors (Blood pressure <BP>, total cholesterol <TC>, low-density lipoprotein cholesterol <LDL>, high-density lipoprotein cholesterol <HDL>, and glucose) among patients with established coronary heart disease (CHD) in Europe, Asia, and the Middle East.


Study population and outcomes

We used cardiovascular risk factors from the SUrvey of Risk Factors (SURF). Details have been reported previously [18, 19]. Briefly, SURF was a clinical audit carried during routine cardiology visit in ten countries among three regions, including Europe (Croatia, Denmark, Ireland, Italy, Northern Ireland, Romania, Russia), Asia (Mainland of China and Taiwan), and Middle East (Saudi Arabia). Within each center, patients aged ≥18 years with a clinical diagnosis of CHD (coronary artery bypass surgery <CABG>, percutaneous coronary intervention <PCI>, acute coronary syndromes <ACS> or stable angina) were recruited between 2012 and 2013. Data on patient demographics (age, sex, and center location), lifestyles (smoking status and physical activity), physical and laboratory measurements (body anthropometry, BP, TC, LDL, HDL, and glucose), and medications were collected by trained research staffs using one-page data collection. BP, lipids, and glucose were measured according to local national guidelines and retrieved directly from medical records.

Air pollution data

We extracted annual average PM2.5 concentrations from the World Health Organization (WHO) database (http://www.who.int/phe/health_topics/outdoorair/databases/modelled-estimates/en/). The database provides estimates of annual average concentration of PM2.5 at a spatial resolution of 0.1° × 0.1°, which is approximately 11 × 11 km at the equator globally. Due to data availability, we used annual average of the year 2014. The estimates are based on the recently developed Data Integration Model for Air Quality [20]. The model estimates PM2.5 using satellite retrievals of aerosol optical depth, chemical transport models, population estimates, topography and ground measurements from 6003 stations worldwide. A Bayesian hierarchical model is used to integrate these information sources [20]. The major advantage of the model is that estimates are available from a consistent method globally, as opposed to ground measurements, which are concentrated in limited regions of the world.

We additionally collected data in 2013 for European centers from countries that report measurements data to the European Environment Agency using the Airbase database (https://www.eea.europa.eu/data-and-maps/data/airbase-the-european-air-quality-database-7). For the 17 districts in the city of Beijing we also obtained online PM2.5 data from the Beijing Municipal Environmental Protection Bureau for the year 2013.

Linkage of the data sources

The postal address of each clinic was transformed into geographical coordinates-the latitude, longitude coordinate system (5 digits)-using Google Earth. We first linked PM2.5 data from the background monitoring stations in the town itself. If no station was available, we estimated PM2.5 from the more frequently measured pollutant PM10 if available or used the average of the nearest two background stations if PM10 was also not available. We used country-specific ratios from EEA database to convert PM10 into PM2.5 fractions if available. If not available, we used PM2.5/PM10 = 0.60 from a large European project or a generic PM2.5/PM10 ratio of 0.60 from a large European project if no country-specific estimates were available [21]. For a small town, we used regional stations and for a large city urban stations.

Statistical analyses

The associations of cardiovascular risk factors with an increase of 10 μg/m3 in PM2.5 were assessed by adjusted mixed-effect generalized estimation equation models. Patient’s characteristics and lifestyles varied country by country and thus, we included all available potential confounding factors related to both cardiovascular risk factors and PM2.5, including sex, age, and individual risk factors (physical activity <low, moderate, vigorous>, smoking status <current smoker, ex-smoker, never>, and body mass index <BMI>) [22]. All patients with established CHD were expected to be on cardiovascular medications to prevent the recurrence of cardiac event irrespective of geographical areas. Thus, cardiovascular medications were not included as a confounder. We further adjusted for country as a fixed covariate as a proxy for potential unknown and known confounders for which we did not have individual information. All outcomes were also nested within center (the random effect) to allow for clustering within centers.

Imputed data were analyzed in the primary analysis. There were about less than 4% missing data for all variables (Appendix Table A). Ten datasets were imputed for missing data with multivariate imputation by chained equations (MICE package in R) [23]. Briefly, MICE predicts missing data by iteratively optimizing a series of regression models using other potentially predictive variables such as basic demographics and geographic area. The continuous variables including height, weight, BP, TC, LDL, HDL, and glucose were imputed by predictive mean matching and the categorical data including smoking status and physical activity were imputed with logistic regression.

Because of uncertainty of the shape of the concentration response function at high concentrations, we performed sensitivity analyses excluding the two countries with the highest PM2.5 levels (China and Saudi Arabia) (Appendix Figure A and B). We further analyzed associations of PM2.5 retrieved from the Airbase for European countries and the database from the Beijing Municipal Environmental Protection Bureau for China with the same statistical strategy.

Statistical analyses were performed by using ‘MICE’ and ‘GEEPACK’ packages in R [23, 24]. All tests were two tailed with statistical significance assumed at the 0.05 level.


We first describe the collected data and associations between air pollution and cardiovascular risk factors and then summarize the potential limitations of using existing audit data.

Baseline characteristics

A total of 8392 SURF patients were included. The mean age of all patients was 64.9 years; 29.6% were women; 16% reported current smoker (Table 1). The average overall systolic blood pressure (SBP), diastolic blood pressure (DBP), TC, LDL, HDL, and glucose were 131.1 mmHg, 75.8 mmHg, 4.2 mmol/L, 2.4 mmol/L, 1.1 mmol/L, and 7.5 mmol/L, respectively. The average PM2.5 level from WHO database was 38.1 μg/m3, ranging from 10.1 μg/m3 in Ireland to 92.7 μg/m3 in Saudi Arabia. Appendix Figure B illustrates the large variation of individual outcome variables, especially within countries.

Table 1

Description of patient characteristics, cardiovascular risk factors, and air pollutants.

Overall Europe Croatia Denmark Ireland Italy NI Russia Romania KSA Taiwan Beijing

No. 8392 5001 1223 300 1716 771 159 463 369 1509 732 1150
No. of Centre 71 51 8 1 11 14 2 8 7 4 4 12
Age, years 64.9 ± 11.2 64.9 ± 10.7 65.2 ± 10.8 65.5 ± 11.6 63.6 ± 10.4 68.2 ± 10.1 63.9 ± 11.0 65.1 ± 10.0 62.6 ± 11.8 62.2 ± 12.1 67.0 ± 13.0 66.7 ± 9.5
Women, % 29.6 26.8 30.8 28.0 23.0 21.5 20.8 38.4 28.5 26.4 30.7 45.3
PM2.5 WHO 38.1 ± 34.5 15.8 ± 5.7 20.7 ± 2.6 11.2* 10.1 ± 0.27 22.6 ± 7.4 10.6 ± 0.6 16.7 ± 4.4 20.5 ± 3.1 92.7 ± 31.6 34.1 ± 3.9 67.4 ± 14.7
PM2.5 Local NA 15.0 ± 5.9 19.6 ± 4.1 10.6* 10.5 ± 0.5 20.3 ± 7.2 10.2 ± 1.0 NA 16.4 ± 3.4 NA NA 86.3 ± 12.6
Smoker, % 16.1 17.9 20.5 23.2 14.6 14.6 20.0 25.5 17.3 10.3 20.0 12.3
Exercise, % 54.2 57.2 55.0 25.3 69.2 46.0 44.1 62.4 53.4 45 48.8 56.4
BMI, kg/m2 28.1 ± 34.5 28.6 ± 4.7 28.7 ± 4.2 28.3 ± 4.8 29.0 ± 5.0 26.9 ± 4.0 30.7 ± 5.5 29.5 ± 5.1 28.5 ± 4.5 30.3 ± 6.0 28.0 ± 4.0 25.0 ± 2.9
HTN, % 74.5 71.7 83.5 63.6 58.6 72.8 63.0 88.9 75.5 83.9 70.8 81.2
Diabetes, % 34.2 25.2 3.0 NA 17.9 30.2 14.8 24.0 26.7 76.1 35.8 45.8
Hyperlipidaemia, % 67.6 68.7 31.9 70.6 68.9 68.1 64.6 62.2 56.4 88.9 49.4 42.7
SBP, mmHg 131.1 ± 18.4 132.1 ± 19.3 133 ± 20.2 131 ± 19.4 134.2 ± 18.5 128.1 ± 18.9 121.8 ± 15.5 140.0 ± 14.5 136.5 ± 23.4 128.1 ± 17.0 131.5 ± 19.4 129.9 ± 13.2
DBP, mmHg 75.8 ± 10.9 76.7 ± 10.7 79.9 ± 11.6 76.0 ± 10.4 74.3 ± 10.0 75.4 ± 9.3 71.3 ± 8.1 79.9 ± 9.1 79.9 ± 12.7 71.3 ± 10.3 75.8 ± 12.7 77.5 ± 8.8
TC, mmol/L 4.2 ± 1.5 4.3 ± 1.7 4.7 ± 1.7 4.1 ± 0.9 4.0 ± 1.0 4.1 ± 2.2 3.9 ± 0.8 3.9 ± 1.8 4.9 ± 2.7 3.8 ± 1.0 4.3 ± 1.0 4.4 ± 1.0
LDL, mmol/L 2.4 ± 1.1 2.4 ± 1.2 2.8 ± 1.3 2.1 ± 0.8 2.1 ± 0.8 2.0 ± 1.0 1.9 ± 0.6 2.7 ± 1.3 3.5 ± 2.0 2.1 ± 0.8 2.7 ± 0.9 2.8 ± 0.9
HDL, mmol/L 1.1 ± 0.4 1.2 ± 0.4 1.1 ± 0.5 1.2 ± 0.4 1.2 ± 0.4 1.1 ± 0.4 1.2 ± 0.4 1.1 ± 0.4 1.1 ± 0.6 0.9 ± 0.3 1.2 ± 0.4 1.2 ± 0.4
Glucose, mmol/L 7.5 ± 1.5 6.1 ± 2.1 6.7 ± 2.6 6.4 ± 2.4 5.8 ± 1.4 6.1 ± 1.8 6.4 ± 2.6 5.3 ± 1.4 6.3 ± 2.7 7.8 ± 3.6 6.5 ± 2.5 6.2 ± 1.6

EU: Europe; NI: Northern Ireland; KSA: Saudi Arabia; SBP: systolic blood pressure; DBP: diastolic blood pressure; TC: total cholesterol; NA: not applicable. ‘Smoker’ was recorded as current smoker; ‘Exercise’ was recorded as adequate physical activities.

Numeric variables are mean ± standard deviation and categorical variables are percentage. Units are years for age, kg/m2 for BMI, mmHg for SBP and DBP, mmol/L for TC, LDL, HDL, and glucose, and μg/m3 for PM2.5.

* Only one centre from Denmark participated in SURF so standard deviation could not be provided.

Associations between PM2.5 and cardiovascular risk factors

Appendix Figure C shows the crude association between PM2.5 and cardiovascular risk factors, indicating weak associations if any.

Globally, a 0.26 mmHg decrease in SBP per 10 μg/m3 increase in PM2.5 was observed (Figure 1). After controlling for country, the observed inverse association with SBP was slightly stronger but with wider confidence intervals (–0.45 mmHg; 95% CI: –0.85, –0.06). There were no statistically significant associations with SBP when the analysis was restricted to the European centers (1.32 mmHg; 95% CI: –6.73, 4.08).

Figure 1 

Changes (95% CI) in blood pressure increase in PM2.5 derived from World Health Organization.

SBP: systolic blood pressure; DBP: diastolic blood pressure.

All analyses were applied with generalized estimating equation model with centre clustered. ‘Full’ adjustment was sex, age, and risk factors (exercise, smoking status, and body mass index). ‘Full+country’ was sex, age, risk factors (exercise, smoking status, and body mass index), and country. Results are presented as changes in mmHg (95% CI).

‘Global*’ presented results are based on all participating countries except China; ‘Global**’ presented results are based on all participating countries except Saudi Arabia.

Similar results were found for DBP: an increase of 10 μg/m3 in PM2.5 was associated with lower DBP (–0.36 mmHg; 95% CI: –0.10, 0.61) and the association tended to be stronger (–0.47 mmHg; –0.73, –0.20) after country adjustment on a global scale. On European level, a similar association between PM2.5 and DBP was observed which became non-significant after country adjustment.

Figure 2 shows the association between PM2.5 and lipid levels. Associations of PM2.5 with lipid levels were not statistically significant on a global scale. After controlling for country non-significant associations remained for TC and HDL; while, an increase of 10 μg/m3 in PM2.5 was associated with an increased LDL level (0.04 mmol/L, 95% CI: 0.01, 0.08). Weak positive associations of TC and LDL were observed among European participants (TC: 0.32 mmol/L; 95% CI: 0.01, 0.62; LDL: 0.30 mmol/L, 95% CI: 0.03, 0.58), which disappeared after adjustment for country. There was no significant association for HDL among European patients with or without adjustment for country.

Figure 2 

Changes (95% CI) in lipids (Total cholesterol, LDL-cholesterol, and HDL-cholesterol) increase in PM2.5 derived from World Health Organization.

TC: total cholesterol; LDL: low-density lipoprotein cholesterol; HDL: high-density lipoprotein cholesterol.

All analyses were applied with generalized estimating equation model with centre clustered. ‘Full’ adjustment was sex, age, and risk factors (exercise, smoking status, and body mass index). ‘Full+country’ was sex, age, risk factors (exercise, smoking status, and body mass index), and country. Results presented as changes in mmol/L (95% CI).

‘Global*’ presented results are based on all participating countries except China; ‘Global**’ presented results are based on all participating countries except Saudi Arabia.

Globally, an increase of 10 μg/m3 PM2.5 was associated with an increased glucose level by 0.10 mmol/L (95% CI: 0.03 to 0.16). For Europe the increase in glucose was 0.30 mmol/L (95% CI: 0.06 to 0.53) (Figure 3). These associations, however, disappeared after adjustment for country.

Figure 3 

Changes (95% CI) in glucose increase in PM2.5 derived from World Health Organization.

All analyses were applied with generalized estimating equation model with centre clustered. ‘Full’ adjustment was sex, age, and risk factors (exercise, smoking status, and body mass index). ‘Full+country’ was sex, age, risk factors (exercise, smoking status, and body mass index), and country. Results are presented as changes in mmol/L (95% CI).

‘Global*’ presented results are based on all participating countries except China; ‘Global**’ presented results are based on all participating countries except Saudi Arabia.

Sensitivity analyses

Separate analyses with exclusion of China and Saudi Arabia (called as ‘global*’ and ‘global**’ in Figures 1, 2, 3) and with local PM2.5 exposure data (Appendix Table B) did not alter the main findings.

Feasibility of data linkage

It is feasible to link existing cardiovascular risk factor data with PM2.5. During the linkage process, some limitations were identified in the various data sources. The air pollution data sources did not always contain data of the exact year of interest and thus we used the nearest by year. The SURF database did not contain individual addresses and hence we used the postal code of the hospital address of the patient as a proxy for the location of exposure to air pollution.


The analyses establish the technical feasibility of developing future data linkage studies but also point at challenges in their interpretation. In the current analysis, the long-term PM2.5 exposure from a consistent global exposure model was linked to individual data on routinely measured CVD risk factors from a large audit of 8,392 CHD patients from 71 centers in Europe, Asia, and the Middle East to explore potential association between air pollution and cardiovascular risk factors. Notably, taking country into account in the analyses materially affected the observed associations. While this adjustment may account for unmeasured confounding and lead to over adjustment.

Associations between air pollution and cardiovascular risk factors

We observed an inverse association of PM2.5 with BP globally and among European participants after adjustment for country, which is in contrast with several previous studies that found positive associations between long-term exposures to PM2.5 and elevated BP [4, 9, 25, 26]. Other studies found non-significant association [10]. For instance, findings from a national population-based study among 1024 elderly Taiwanese participants suggested that an interquartile increase in PM2.5 (48 μg/m3) is associated with 32.1 mmHg (95% CI 21.6–42.6) and 31.3 mmHg (95% CI 25.4–37.1) increases in SBP and DBP, respectively, after controlling age, sex, BMI, smoking, and drinking habitats [27]. A comprehensive meta-analysis among 113,926 patients from 15 European population-based cohort studies, ESCAPE, demonstrated inconsistent relationships between long-term exposure to modeled air pollutants including PM2.5 and BP in each cohort and the pooled results remained non-significant [10]. Studies on mechanisms have suggested that exposure to PM2.5 could instigate acute autonomic imbalance and then lead BP increases [4, 5, 25, 28, 29]. However, our study was conducted in CHD patients who all received cardiovascular medications to control potential risk factors. Consequently, we measured the potential impact of air pollution beyond medical treatment. Future linkage studies would need to include both treated and untreated patients to better investigate the association between air pollution and CHD risk factors.

Some previous evidence suggested that PM2.5 may affect lipid levels but the quantity and quality of these studies is still limited and results are not fully consistent [27, 30, 31]. A large cross-sectional study with 39,863 healthy participants in Denmark demonstrated that the interquartile range (11.3 μg/m3) of PM2.5 was associated with a higher level of TC (0.78 mg/dl; 95% CI: 0.22–1.34) [31]. An animal study also indicated that mice exposed to PM2.5 had significantly higher levels of TC and LDL than those exposed to filtered air [30]. However, effect estimates are typically small and may have little clinical implications.

We observed direct associations of PM2.5 with glucose in both global and European analyses, although these associations attenuated after country adjustment. These findings are in line with previous studies [32, 33]. A cross-sectional study based on Chinese populations reported that both elevated glucose levels and increased type II diabetes prevalence are significantly associated with increased PM2.5 [34]. A review from 21 published studies reported concentrations of PM2.5 to be associated with increased insulin resistance and higher rates of type II diabetes [32]. Mechanisms suggested to link glucose metabolism to PM2.5 with endothelial dysfunction, endoplasmic reticulum stress, insulin signaling abnormalities, and systematic inflammation [5, 12, 33, 34]. Differences in the study characteristics, population characteristics, and exposure duration in different geographic research areas may contribute to the discrepancies in these findings.

Feasibility and challenges

The current study has piloted feasibility to add air pollution exposure using a coherent methodology to the rich database of clinical observations on cardiovascular risk factors from SURF. Data linkage is a robust, valuable and cost-effective research tool for combining individual level data from different sources for maximizing use of these existing database and increasing amounts of data that are being produced in order to: 1) address clinical research questions that require large sample sizes, detailed data on hard-to-reach population, or specified measurements by using a single dataset; 2) generate evidence with a high level of external validity and applicability; 3) reduce participant burden and avoid duplication of effort [16, 35, 36, 37]. This study facilitates data-linkage possibilities to investigate the impact of air pollution on CHD on a global scale, which is important clinical practice. A recent study demonstrated that the contribution of air pollution to CVD is comparable to that of smoking [38]. Such efficient and cost-effective methods enable all healthcare providers to enrich clinical data to investigate novel health-related research questions.


There are several limitations in this study. SURF records CHD management in daily practice. Unlike other epidemiological studies, physical and laboratory measurements are not standardized. Some potential confounders were not available and thus could not be adjusted for. In addition, SURF collected anonymous data and thus only participating center’s locations at aggregated level were linked to air pollution data instead of individual level, which may not reflect actual exposure at the individual level. However, most routine cardiology visits were conducted in local hospitals with the distance between home and clinic generally being less than 10 km as confirmed by SURF national collaborators for 80%–90% of their patients had their residence near hospitals. The lack of individual addresses resulted in that only PM2.5 concentrations were assigned to each center, as PM2.5 is a regionally varying pollutant with limited small-scale spatial variation [21]. Finally, data on other spatially-correlated air pollution factors, such as traffic noise, greenness, and urbanity, were not taken into account and thus not adjusted for, potentially under- or overestimating results.

Furthermore, the air pollution data from WHO database was not available for 2013, which was our year of interest because this coincides with the year of observation for the SURF study. In sensitivity analysis, we further analyzed association with PM2.5 exposure data provided by local resource from 2013 and found that the results are broadly similar. Annual average concentrations may vary from year to year due to variations in weather, but the spatial contrasts in air pollution are typically stable over years and as such it may only have limitedly impacted our findings [22, 23, 24].

While most epidemiological studies of air pollution are based upon more individual exposure assessment, our approach does not invalidate the epidemiological study. First, the selected pollutant PM2.5 mostly varies on a regional scale with limited local variability. In a large monitoring study across Europe, we observed that 81% of the variance was due to between study area variability [21]. Second, people do not spend only time at their residence but in a wider neighborhood, arguing for exposure assessment at a larger scale. Third, the error made by assigning an area-level estimate may lead to Berkson rather than classical error which would not bias air pollution effect estimates but only increase imprecision [39]. Fourth, if the contrast in exposure is large between study areas, assigning an area-level value may be acceptable. Recent studies have applied this approach in settings with large exposure contrasts [40]. Therefore, we would like to clarify that current study was an attempt to use clinical audit data to investigate more health related research questions beyond cardiovascular risk factor management. Further research is needed to validate current findings due to these methodological limitations.

Further direction

We hope that our findings may stimulate linkage studies on cardiovascular risk and disease in primary prevention settings in which relationships may be stronger and the findings less likely to be confounded by medication. Even though clinicians may not be able to change patients’ living environment they should become more aware of the hazards of air pollution and take it into account in their risk assessment and recommendations, such as promoting exercise in less polluted areas.


The current study has demonstrated the feasibility of data. The approach exemplifies the opportunity to assess the impact of the environment on cardiovascular risk factors across large geographic areas. We noted that estimates were highly sensitive to adjustment for country. After country adjustment, PM2.5 levels are marginally associated with increases in LDL cholesterol and decreases in BP. The implication is that similar global studies should aim at multiple centers per country with sufficient within country exposure contrast to balance any effects of over adjustment.

Data Accessibility Statement

The data that support the findings of this study are available from the SURF project, which have been published previously. The references have been included in the current study.

Additional Files

The additional files for this article can be found as follows:

Appendix Table A

Missing data in SURF. DOI: https://doi.org/10.5334/gh.877.s1

Appendix Table B

Associations between CVD risk factors and PM2.5 retrieved from WHO, airbase for European countries and local government database for China. DOI: https://doi.org/10.5334/gh.877.s2

Appendix Table C

Risk factors, stratified by centres. DOI: https://doi.org/10.5334/gh.877.s3

Appendix Figure A

PM2.5 distribution. DOI: https://doi.org/10.5334/gh.877.s4

Appendix Figure B

Variations of individual outcome variable by country. DOI: https://doi.org/10.5334/gh.877.s5

Appendix Figure C

Crude association between PM2.5 and cardiovascular risk factors. DOI: https://doi.org/10.5334/gh.877.s6

Ethics and Consent

SURF is a clinical audit without intervention and follow-up involved. We are suggested that ethics approval was not needed.


We gratefully acknowledge the support from Jie Chen from Utrecht University for her support with regard to air pollutant data collection from the Chinese government.

Funding Information

MZ is supported by a grant from the Netherlands Organization for Scientific Research (NWO; grant number: NWO was not involved in the design, execution, analysis, or writing up of the study. IV is supported by a grant from the Dutch Heart Foundation (grant DHF project ‘Facts and Figures’).

Competing Interests

The authors have no competing interests to declare.

Author Contributions

MZ, GH, IG, DEG, KK, and IV conceived and designed the study. MZ, GH, and IV analysed and interpreted the data. MZ drafted the manuscript and all authors contributed to critical revision of the manuscript.


  1. Vos T, Abajobir AA, Abbafati C, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017; 390: 1211–1259. DOI: https://doi.org/10.1016/S0140-6736(17)32154-2 

  2. Perk J, De Backer G, Gohlke H, et al. European Guidelines on Cardiovascular Disease Prevention in Clinical Practice (Version 2012). Int J Behav Med. 2012; 19: 403–488. DOI: https://doi.org/10.1007/s12529-012-9242-5 

  3. Piepoli MF, Hoes AW, Agewall S, et al. 2016 European Guidelines on cardiovascular disease prevention in clinical practice: The Sixth Joint Task Force of the European Society of Cardiology and Other Societies on Cardiovascular Disease Prevention in Clinical Practice (constituted by rrepresentatives of 10 societies and by invited experts). Atherosclerosis. 2016; 252: 207–274. DOI: https://doi.org/10.1016/j.atherosclerosis.2016.05.037 

  4. Cosselman KE, Navas-acien A, Kaufman JD. Environmental factors in cardiovascular disease. Nat Rev Cardiol. 2015; 12: 627–642. DOI: https://doi.org/10.1038/nrcardio.2015.152 

  5. Newby DE, Mannucci PM, Tell GS, et al. Expert position paper on air pollution and cardiovascular disease. Eur Heart J. 2015; 36: 83–93. DOI: https://doi.org/10.1093/eurheartj/ehu458 

  6. Arden Pope C, Burnett RT, Turner MC, et al. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: Shape of the exposure-response relationships. Environ Health Perspect. 2011; 119: 1616–1621. DOI: https://doi.org/10.1289/ehp.1103639 

  7. Cesaroni G, Forastiere F, Stafoggia M, et al. Long term exposure to ambient air pollution and incidence of acute coronary events: Prospective cohort study and meta-analysis in 11 European cohorts from the ESCAPE project. BMJ. 2014; 348: f7412. 

  8. Shanley RP, Hayes RB, Cromar KR, et al. Particulate Air Pollution and Clinical Cardiovascular Disease Risk Factors. Epidemiology. 2015; 27: 291–298. DOI: https://doi.org/10.1097/EDE.0000000000000426 

  9. Brook RD, Rajagopalan S, Pope A, et al. Particulate matter air pollution and cardiovascular disease. Circulation. 2010; 121: 2331–2378. DOI: https://doi.org/10.1161/CIR.0b013e3181dbece1 

  10. Fuks K, Weinmayr G, Foraster M, et al. Arterial Blood Pressure and Long-Term Exposure to Traffic-Related Air Pollution: An Analysis in the European Study of Cohorts for Air Pollution Effects (ESCAPE). Environ Health Perspect. 2014; 896: 896–906. DOI: https://doi.org/10.1289/ehp.1307725 

  11. Gakidou E, Afshin A, Abajobir AA, et al. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990-2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017; 390: 1345–1422. DOI: https://doi.org/10.1016/S0140-6736(17)32366-8 

  12. Strak M, Janssen N, Beelen R, et al. Long-term exposure to particulate matter, NO2 and the oxidative potential of particulates and diabetes prevalence in a large national health survey. Environ Int. 2017; 108: 228–236. DOI: https://doi.org/10.1016/j.envint.2017.08.017 

  13. Eze IC, Hemkens LG, Bucher HC, et al. Association between ambient air pollution and diabetes mellitus in Europe and North America: systematic review and meta-analysis. Env Heal Perspect. 2015; 123: 381–389. DOI: https://doi.org/10.1289/ehp.1307823 

  14. Yang B, Qian Z, Howard SW, et al. Global association between ambient air pollution and blood pressure: A systematic review and meta-analysis. Enviromental Pollut. 2018; 235: 576–588. DOI: https://doi.org/10.1016/j.envpol.2018.01.001 

  15. The World Bank. Indicators WD. The World Bank. http://data.worldbank.org/data-catalog/world-development-indicators. 

  16. Andrew NE, Sundararajan V, Thrift AG, et al. Addressing the challenges of cross-jurisdictional data linkage between a national clinical quality registry and government-held health data. Aust N Z J Public Health. 2016; 40: 436–442. DOI: https://doi.org/10.1111/1753-6405.12576 

  17. Martin-Sanchez F, Verspoor K. Big data in medicine is driving big changes. Yearb Med Inform. 2014; 9: 14–20. DOI: https://doi.org/10.15265/IY-2014-0020 

  18. Cooney MT, Reiner Z, Sheu W, et al. SURF – SUrvey of Risk Factor management: First report of an international audit. Eur J Prev Cardiol. 2014; 21: 813–822. DOI: https://doi.org/10.1177/2047487312467870 

  19. Zhao M, Cooney MT, Klipstein-Grobusch K, et al. Simplifying the audit of risk factor recording and control: A report from an international study in 11 countries. Eur J Prev Cardiol. 2016; 23: 1202–1210. DOI: https://doi.org/10.1177/2047487316647827 

  20. Shaddick G, Thomas M, Green A, et al. Data integration model for air quality: A hierarchical approach to the global estimation of exposures to ambient air pollution. J R Stat Soc Ser C Appl Stat. 2017; 67: 231–253. DOI: https://doi.org/10.1111/rssc.12227 

  21. Eefrens M, Tsai M, Ampe C, et al. Sptial variation of PM2.5, PM10, PM2.5 absorbance and PM coarse concentrations between and within 20 European study areas and the relationship with NO2: Results from the ESCAPE project. Atmos Env. 2012; 62: 303–317. DOI: https://doi.org/10.1016/j.atmosenv.2012.08.038 

  22. Danaei G. Metabolic mediators of the effects of body-mass index, overweight, and obesity on coronary heart disease and stroke: A pooled analysis of 97 prospective cohorts with 1·8 million participants. Lancet. 2014; 383: 970–983. DOI: https://doi.org/10.1016/S0140-6736(13)61836-X 

  23. Groothuis-oudshoorn K. MICE: Multivariate Imputation by Chained. J Stat Softw. 2011; 45: 3. DOI: https://doi.org/10.18637/jss.v045.i03 

  24. Halekoh U, Hojsgaard S. The R package geepack for generalized estimating equations. Jounral Stat Softw. 2006; 15. DOI: https://doi.org/10.18637/jss.v015.i02 

  25. Chen H, Burnett R, Kwong J, et al. Spatial associations between ambient fine particulate matter and incident hypertension. Circulation. 2014; 129: 562–569. DOI: https://doi.org/10.1161/CIRCULATIONAHA.113.003532 

  26. Brook RD, Urch B, Dvonch JT, et al. Insights into the mechanisms and mediators of the effects of air pollution exposure on blood pressure and vascular function in healthy humans. Hypertension. 2009; 54: 659–667. DOI: https://doi.org/10.1161/HYPERTENSIONAHA.109.130237 

  27. Chuang K, Yan Y, Chiu S, et al. Long-term air pollution exposure and risk factors for cardiovascular diseases among the elderly in Taiwan. Occup Env Med. 2011; 68: 64–68. DOI: https://doi.org/10.1136/oem.2009.052704 

  28. Chahine T, Baccarelli A, Litonjua A, et al. Particulate air pollution, oxidative stress genes, and heart rate variability in an elderly cohorts. Env Heal Perspect. 2007; 115: 1617–1622. DOI: https://doi.org/10.1289/ehp.10318 

  29. Widdicombe J, Lee L-Y. Airway reflexes, autonomic function, and cardiovascular responses. Env Heal Perspect. 2001; 109: 579–584. DOI: https://doi.org/10.1289/ehp.01109s4579 

  30. Chen T, Jia G, Wei Y, et al. Beijing ambient particle exposure accelerates atherosclerosis in ApoE knockout mice. Toxicol Lett. 2013; 223: 146–153. DOI: https://doi.org/10.1016/j.toxlet.2013.09.004 

  31. Sørensen M, Hjortebjerg D, Eriksen KT, et al. Exposure to long-term air pollution and road traffic noise in relation to cholesterol: A cross-sectional study. Environ Int. 2015; 85: 238–243. DOI: https://doi.org/10.1016/j.envint.2015.09.021 

  32. Meo SA, Memon AN, Sheikh SA, et al. Effect of environmental air pollution on type 2 diabetes mellitus. Eur Rev Med parmacology Sci. 2015; 19: 123–128. 

  33. Rajagopalan S, Brook RD. Air Pollution and Type 2 Diabetes Mechanistic Insights. Diabetes. 2012; 61: 3037–3045. DOI: https://doi.org/10.2337/db12-0190 

  34. Liu C, Yang C, Zhao Y, et al. Associations between long-term exposure to ambient particulate air pollution and type 2 diabetes prevalence, blood glucose and glycosylated hemoglobin levels in China. Environ Int. 2016; 93: 416–421. DOI: https://doi.org/10.1016/j.envint.2016.03.028 

  35. Harron K, Dibben C, Boyd J, et al. Challenges in administrative data linkage for research. Big Data Soc. 2017; 1–12. DOI: https://doi.org/10.1002/9781119072454.ch1 

  36. Boyd JH, Randall SM, Ferrante AM, et al. Technical challenges of providing record linkage services for research. BMC Med Inform Decis Mak. 2014; 14: 23. DOI: https://doi.org/10.1186/1472-6947-14-23 

  37. Bohensky MA, Jolley D, Sundararajan V, et al. Data Linkage: A powerful research tool with potential problems. BMC Health Serv Res. 2010; 10: 346. DOI: https://doi.org/10.1186/1472-6963-10-346 

  38. Lelieveld J, Klingmu K, Pozzer A, et al. Cardiovascular disease burden from ambient air pollution in Europe reassessed using novel hazard ratio functions. Eur Heart J. 2019; 40: 1590–1596. DOI: https://doi.org/10.1093/eurheartj/ehz135 

  39. Hoek G. Methods for Assessing Long-Term Exposures to Outdoor Air Pollutants. Curr Environ Heal reports. 2017; 4: 450–462. DOI: https://doi.org/10.1007/s40572-017-0169-5 

  40. Yin P, Brauer M, Cohen A, et al. Long-term fine particulate matter exposure and nonaccidental and cause-specific mortality in a large national cohort of Chinese men. Environ Health Perspect. 2017; 125: 117002. DOI: https://doi.org/10.1289/EHP1673