Aslak Harbo Poulsen,Mette Sørensen,Ulla Arthur Hvidtfeldt,Jesper H. Christensen,Jørgen Brandt,Lise Marie Frohn,Matthias Ketzel,Christopher Andersen,Ole Raaschou-Nielsen
Source-Specific Air Pollution Including Ultrafine Particles and Risk of Myocardial Infarction: A Nationwide Cohort Study from Denmark
Background: Air pollution is negatively associated with cardiovascular health. Impediments to efficient regulation include lack of knowledge about which sources of air pollution contributes most to health burden and few studies on effects of the potentially more potent ultrafine particles (UFP). Objective: The authors aimed to investigate myocardial infarction (MI) morbidity and specific types and sources of air pollution. Methods: We identified all persons living in Denmark in the period 2005–2017, age y and never diagnosed with MI. We quantified 5-y running time-weighted mean concentrations of air pollution at residencies, both total and apportioned to traffic and nontraffic sources. We evaluated particulate matter (PM) with aerodynamic diameter (), (UFP), elemental carbon (EC), and nitrogen dioxide (). We used Cox proportional hazards models, with adjustment for time-varying exposures, and personal and area-level demographic and socioeconomic covariates from high-quality administrative registers. Results: In this nationwide cohort of 1,964,702 persons (with person-years of follow-up and 71,285 cases of MI), UFP and were associated with increased risk of MI with hazard ratios (HRs) per interquartile range (IQR) of 1.040 [95% confidence interval (CI): 1.025, 1.055] and 1.053 (95% CI: 1.035, 1.071), respectively. HRs per IQR of UFP and from nontraffic sources were similar to the total (1.034 and 1.051), whereas HRs for UFP and from traffic sources were smaller (1.011 and 1.011). The HR for EC from traffic sources was 1.013 (95% CI: 1.003, 1.023). from nontraffic sources was associated with MI (; 95% CI: 1.034, 1.062) but not from traffic sources. In general, nontraffic sources contributed more to total air pollution levels than national traffic sources. Conclusions: and UFP from traffic and nontraffic sources were associated with increased risk of MI, with nontraffic sources being the dominant source of exposure and morbidity. https://doi.org/10.1289/EHP10556
Source-Specific Air Pollution Including Ultrafine Particles and Risk of Myocardial Infarction: A Nationwide Cohort Study from Denmark

Air pollution is negatively associated with cardiovascular health. Impediments to efficient regulation include lack of knowledge about which sources of air pollution contributes most to health burden and few studies on effects of the potentially more potent ultrafine particles (UFP).

The authors aimed to investigate myocardial infarction (MI) morbidity and specific types and sources of air pollution.

We identified all persons living in Denmark in the period 2005–2017, age y and never diagnosed with MI. We quantified 5-y running time-weighted mean concentrations of air pollution at residencies, both total and apportioned to traffic and nontraffic sources. We evaluated particulate matter (PM) with aerodynamic diameter (), (UFP), elemental carbon (EC), and nitrogen dioxide (). We used Cox proportional hazards models, with adjustment for time-varying exposures, and personal and area-level demographic and socioeconomic covariates from high-quality administrative registers.

In this nationwide cohort of 1,964,702 persons (with person-years of follow-up and 71,285 cases of MI), UFP and were associated with increased risk of MI with hazard ratios (HRs) per interquartile range (IQR) of 1.040 [95% confidence interval (CI): 1.025, 1.055] and 1.053 (95% CI: 1.035, 1.071), respectively. HRs per IQR of UFP and from nontraffic sources were similar to the total (1.034 and 1.051), whereas HRs for UFP and from traffic sources were smaller (1.011 and 1.011). The HR for EC from traffic sources was 1.013 (95% CI: 1.003, 1.023). from nontraffic sources was associated with MI (; 95% CI: 1.034, 1.062) but not from traffic sources. In general, nontraffic sources contributed more to total air pollution levels than national traffic sources.

and UFP from traffic and nontraffic sources were associated with increased risk of MI, with nontraffic sources being the dominant source of exposure and morbidity. https://doi.org/10.1289/EHP10556

A recent review of experimental and clinical studies concluded that air pollution can cause oxidative stress as well as local and systemic inflammation and that it affects vascular tone, fibrinolysis, platelet activity, and plaque formation. Correspondingly, many epidemiological studies have found associations between air pollution and a range of cardiovascular end points with stronger evidence for particulate matter (PM) than gaseous components. With regard to incident myocardial infarction (MI), some studies found higher risk in association with exposure to nitrogen dioxide (), whereas other studies did not. For PM with aerodynamic diameter (), a recent meta-analysis found a hazard ratio (HR) of 1.08 per (95% confidence intervals (CI): 0.99, 1.18), with high heterogeneity between studies and concluded that there was suggestive evidence that long-term exposure to is associated with MI. Most subsequent studies have reported an elevated risk with a few exceptions. Ultrafine particles (UFP), PM , may have health effects disproportionate to their mass contribution to , because they have a large reactive surface area and may penetrate beyond the respiratory tract. Studies of short-term exposure to UFP suggest an association with stroke, blood pressure, systemic inflammation, and autonomic tone. Few epidemiological studies have investigated long-term cardiovascular effects of UFP. Two studies found associations between UFP and risk of ischemic heart disease and stroke, two studies found associations between UFP and risk of MI, and one study found an association between UFP and hypertension. These previous studies all indicate an association between UFP and cardiovascular diseases, but they each have their limitations, e.g., small size, exposure assessed only at single point in time or with at spatial resolution at the kilometer scale, or having only area-level covariates. Several studies have linked air pollution from elemental carbon (EC) or the closely related measures black smoke (BS) and absorbance, with higher CVD mortality. Two studies showed associations between absorbance and coronary heart disease, whereas another study found no association between black carbon and ischemic heart disease. The chemical composition and toxicity of air pollution depends on the emission source, as does the relative contribution to total air pollution exposure. Knowledge about the relative contribution to disease burden of sources such as traffic may therefore facilitate efficient political regulation to improve public health. Most long-term studies have, however, not made such subdivisions, and the few results relating to MI are inconsistent. In a nationwide cohort with extensive data on personal and area-level covariates, we assessed time-varying total and source-specific air pollution at individual residences with a state-of-the-art multiscale air pollution model system and evaluated the association of UFP, , elemental carbon (EC), and with risk of MI.

The study was conducted in Denmark, where a unique personal identification number has been allocated to all Danish residents since 1968. The identification number allows all residents to be traced in and across all health and administrative registers. We identified Danish residents in the Civil Registration System, which also provided complete residential histories (including dates of moving in and out) and information about emigration from Denmark. We identified all who were born after 1920 (information on education is not available for those born before 1921), who lived in Denmark from birth or 1979 onward, and who were present in Denmark on 1 January 2005, and were 50 y of age or older any time after that date. From this population (), we excluded 60,967 with prevalent MI at baseline (see below) and 22,613 with missing information about one or more covariates, leaving a final cohort of 1,964,702 persons. The cohort was followed up from 1 January 2005 to 31 December 2017. Thus, study entry (baseline) of cohort members was 1 January 2005 or the date of turning 50 y of age, whichever came last. In accordance with Danish law, entirely register-based studies do not require ethical approval.

We identified all cases of MI [International Classification of Diseases, 8th Revision (ICD8) (ICD8: 410 or ICD10: I21)] recorded as primary diagnosis in the Danish National Patient Register or with the same codes as primary cause of death in the Danish Register of Cause of Death. Individuals registered with an MI in one of those registries, from start of the registration (patient registry: 1977; mortality registry: 1970) and until study entry were excluded as prevalent cases. All other people with a first, incident diagnosis of MI after study entry (age 50 y or 1 January 2005, whichever came last) were counted as cases.

We established exact geocodes of all residential addresses in Denmark from the Building and Housing Registry and modeled outdoor air pollution for the period 1995–2017 for these dwellings. We used the fully deterministic Danish DEHM/UBM/AirGIS modeling system, which is based on detailed emission inventories and meteorological data from a weather forecast model and then physical (atmospheric transport and dispersion) and chemical processes in the atmosphere are modeled. The system adds contributions to air pollution concentrations at three scales: a) the long-range transported regional background, modeled with the Danish Eulerian Hemispheric Model (DEHM) assessing the contribution from emissions from the entire northern hemisphere; b) the local background contribution at a resolution calculated from Danish emissions of air pollution (more details below) modeled with the Urban Background Model (UBM); and c) the contribution from the traffic on the address street modeled with the Operational Street Pollution Model (OSPM), with detailed input parameters including traffic load, composition, and emission factors as well as street and building configuration, and meteorology. For each address one final estimate was achieved by adding DEHM and UBM, and for addresses with vehicles per day, contribution from local traffic calculated by OSPM was also added. In this coupling of models, it was ensured that emission sources were not counted twice. For all investigated pollutants hourly, address-level concentrations of air pollutants were then aggregated as monthly means. We modeled , EC, and mass concentration as well as particle number concentration (PNC) as an indicator for UFP (denoted in this article as UFP); further details about the modeling system can be found elsewhere. In brief, we modeled PNC based on an extension of the DEHM/UBM/AirGIS system with the aerosol dynamics module M7. The M7 module represents the aerosol in four size fractions: Nucleation mode (typical size below ), Aitken mode (), Accumulation Mode () and Coarse Mode (); and it handles conversion between size modes. The modes output from the DEHM feed into the UBM and OSPM models without including the particle dynamics. The size range was chosen to facilitate validation against measurement data. We scaled UFP particle number emissions to source-specific and EC mass emissions to derive a consistent emission inventory for PNC. As input data for the modeling system, we used high-resolution emission inventories for Denmark for the period 1990–2020 that are based on registers of high quality and detail, e.g., the National Road and Traffic Database. The emission inventories are subject to an annual review from international experts to ensure that the calculations are done in accordance with best practices and meet the internationally established quality criteria. The inventories are based on the standardized European Selected Nomenclature for Air Pollution (SNAP) code classification of emission categories. The Danish emission estimations are done at the most detailed SNAP level and can be aggregated to the 11 SNAP main categories (Table S1). For large point sources such as power plants or road traffic high-quality estimates are ensured by detailed measurements and detailed information on technology and abatement. This allowed us to model air pollution both with and without the emissions from Danish road traffic (SNAP 07 codes) to the local background (UBM) and the traffic on the address street. Subtracting the two modeling results from each other enabled apportioning of total pollutant levels into national road traffic (denoted as “traffic” in this article) and nontraffic sources, with the latter also including traffic and other sources from surrounding countries. The contribution from road traffic to secondary particles (e.g., nitrate) is therefore not included in the Danish road traffic scenario in this setup. Based on available emission factors, the contribution to from national traffic was further divided into contributions from tailpipe (from motorized road transport vehicles) and nontailpipe (resuspended dust, dust from brakes, clutches, tires, road abrasion, and evaporation from vehicles) by running air pollution models with and without these emissions. Using modeled monthly mean concentrations, we then calculated running mean exposure 1, 5, and 10 y back in time for each cohort member. These calculations were done as monthly updated, running averages of air pollution at the residential address(es) held by each cohort member. These averages took into account the time each person stayed at each address and provided time-weighted averages (TWA). The running 5-y TWA was our a priori main exposure metric. The final result of the complete modeling exercise was time-varying estimates of , UFP, EC, and in total and from traffic and nontraffic sources. For , traffic estimates were further subdivided into tailpipe and nontailpipe emissions. For all 14 estimates, we estimated TWA over 1, 5, and 10 y, with 5 y being the primary averaging period. Comparisons between modeled and measured air pollution concentrations show generally good performance of the Danish DEHM/UBM/AirGIS modeling system. For total and EC, the correlation coefficients were in the range of 0.67–0.85 and 0.77–0.79, respectively, across different locations and measurement periods, and the correlations were similarly high for . For annual means of UFP, correlations between predictions and measurements were 0.86, 0.87, and 0.95 for validations at the regional, urban, and street scale, respectively. The source apportioning cannot be validated directly against measurements, but as described previously, the input data are of highest standard, and the model reproduces concentrations in high traffic roads well (Figure S1).

We selected potential confounders a priori. Statistics Denmark provided annually updated information about country of origin (“Danish origin” defined as being a Danish citizen or having at least one Danish citizen parent; all other classified as “other country of origin”), civil status (married/cohabiting, other), highest attained education (mandatory, short: vocational training, secondary or y of postsecondary education, medium/long: y of postsecondary education), occupational status (high-level white-collar, low-level white-collar, blue-collar, unemployed, retired), personal income (quintiles), and household per capita income (quintiles). Sex- and calendar year–specific quintiles were used for personal and household income to allow for inflation and sex-specific income disparities. In 2017, the 20th, 40th, 60th, and 80th personal income percentiles in Danish kroner (kr) were 142,926, 198,363, 249,053, and 312,717 kr for women and 147,235, 219,142, 281,097, and 370,134 kr for men. The corresponding percentiles of household income per adult in household were 174,497, 229,476, 286,079, and 365,263 kr for women and 171,817, 232,821, 290,162, and 369,949 kr for men. To account for potential residual confounding and potential socioenvironmental disparities, we included a comprehensive array of parish-level covariates reflecting aspects of deprivation. In the year 2017, there were 2,160 parishes in Denmark with a mean area of and a median population of 1,032 persons. For socioeconomic area-level covariates, we obtained annually updated information at the parish level about the proportion of inhabitants with a criminal record, living in single-parent households, with only basic education, with manual labor, with income in lowest quartile, living in social housing, and with non-Western background. A Western background was defined as being a citizen of or having at least one parent who is a citizen of a Western country (Denmark and the rest of the EU, the UK, Norway, Iceland, Liechtenstein, Switzerland, Andorra, San Marino, the Vatican, Australia, New Zealand, the United States, and Canada). Persons with missing covariates were excluded.

We used Spearman rank correlation coefficient to evaluate correlations between air pollutants and covariates. Cox-proportional hazards models with age as time axis were used to calculate HRs and associated 95% CI for 5-y TWA air pollution and MI. Cohort members were followed from the age of 50 y or 1 January 2005, whichever came last, until death, MI, more than 14 consecutive days of unknown address, emigration, or 31 December 2017, whichever came first. We evaluated associations linearly per interquartile range (IQR) of exposure to facilitate comparison of risk estimates between pollutants (our primary analyses) and per fixed increment to facilitate comparison with results from other studies. We also performed categorical analysis based on percentiles of exposure (, , , , , , and 95–100) as a graphical illustration of the shape of the exposure–response function. We used Akaike information criteria (AIC) to compare the model fit of linear and categorical representation of pollutants. We evaluated three models with increasing levels of adjustment. Model 1 was the crude model only adjusted for age, sex, and calendar year (categorical, in 2-y categories). Model 2 was additionally adjusted for personal covariates (country of origin, civil status, educational level, occupational status, personal and household income). Model 3 (our main model) was additionally adjusted for all eight parish-level factors: proportion of inhabitants with non-Western background, living in single-parent households, with only basic education, with manual labor, with personal income in the lowest quartile, with household income in the lowest quartile, living in social housing, and with a criminal record. All variables (except sex and country of origin) were time-dependent, and subjects were allowed to change between levels of exposure and covariates over time. Therefore, exposures and covariates changed with time and, at any moment during follow-up, reflected exactly the conditions pertaining to that specific instance. In sensitivity analyses, we evaluated shorter (1-y) and longer (10-y) TWA periods, and we evaluated effect stratified by age (50–69, 70–80, y) and sex. Because both air pollution and stroke incidence decreased over the study period, we ran three models as sensitivity analyses: a) including exact year of entry as strata; b) including running calendar year of observation (2-y period) as strata rather than as a covariate; and c) analysis restricted to those entering the study in the year 2005. All statistical analyses were performed in SAS (version 9.4; SAS Institute Inc.).

The final cohort of 1,964,702 persons accrued a total of 18,309,318 person-years of follow-up and 71,285 MI events. Most (71%) cohort members entered the cohort in 2005 (Figure S2). Table 1 describes the cohort members at entry, where median 5-y average concentrations of air pollutants were: (10th–90th percentile: 9.1–12.2), (10th–90th percentile: 7,963–15,695), (10th–90th percentile: 0.5–1.0), and (10th–90th percentile: 10.5–23.7). High concentrations of UFP or were associated with higher age, being retired, being female, not living with a spouse, not being a blue-collar worker, and being of non-Danish origin. Further, higher concentrations at the residence were associated with living in a parish with more social housing and non-Western immigrants (Table S2). For the vast majority of addresses, national road traffic sources contributed less to air pollution than nontraffic sources, which was most pronounced for UFP and (Table 1; Figure S3). Spearman rank correlations between pollutants ranged from 0.26 to 0.99. Contributions to from traffic tailpipes and traffic nontailpipe sources were highly correlated with each other and with their sum ( from traffic) () but only moderately correlated with from all sources () (Table S3). The correlation between area-level covariates and air pollutants and individual socioeconomic covariates were generally very low to moderate (Table S4). Pollutants averaged over 1, 5, and 10 y were highly correlated () (Table S5). Results from our main model (full adjustment) showed HRs per IQR of 1.053 (95% CI: 1.035, 1.071) per , 1.040 (95% CI: 1.025, 1.055) per UFP, 1.027 (95% CI: 1.013, 1.040) per , and 1.009 (95% CI: 1.000, 1.019) per EC (Table 2). Results from the crude model showed HRs below 1. The HRs per IQR were stronger for the nontraffic contributions than for the traffic contributions to , UFP, and . The pattern was opposite for EC, where only the association with the contribution from traffic deviated from the null. The tailpipe and nontailpipe contributions to from traffic showed similar associations with MI. Table S6 shows HRs per fixed increase in air pollution. Results based on categories for , UFP, and , indicated increase in risk over the entire exposure range (Figure 1A–D), which was similar for the nontraffic contribution to UFP and . For from nontraffic sources, there was some indication of a leveling off at higher exposure levels (Figure S4A–D; Table S7). For traffic contributions, categorical analysis indicated a lower threshold for the associations with the traffic contributions. (Figure S5 and S6; Tables S7 and S8). Generally, the linear representation of exposure offered a better or only marginally worse fit of data than the categorical model (: 0.00–1.54). Exceptions included traffic, tailpipe and nontailpipe, and nontraffic, where the categorical model was superior (Table S9). Sensitivity analyses showed that the association between air pollution and MI was stronger in older age groups, with only little evidence of an association before age 70 y (Table 3). Except for , associations tended to be stronger in males, albeit with widely overlapping CIs. HRs changed little when averaging air pollution over 1 or 10 y (Table S10). Including year of entry or running calendar year as strata produced results identical to those of the main model (Table S11).

In this nationwide register-based cohort, with person-years of follow-up, 71,285 cases, and exposure levels well below EU-regulatory limits (annual mean: , ), air pollution was associated with a higher risk of MI. The associations were stronger among older people. and UFP showed the highest HRs per IQR. Contributions to and UFP from nontraffic sources showed higher HRs than contributions from national road traffic when expressed per IQR. In general, nontraffic sources contributed more to total air pollution levels than national road traffic sources. For ease of comparison, risk estimates from other studies have in the following discussion been rescaled to the IQR of the present study by means of this formula: . We found to be associated with MI with a HR of 1.053 (95% CI: 1.035, 1.071) per IQR (), which is statistically compatible with the majority of previous studies., There are, however, also studies suggesting both weaker and stronger associations. The studies with weaker associations generally modeled exposure with lower temporal or spatial resolution than the present study, and two of them found stronger risks after relatively crude adjustment for potential covariates. Exposure misclassification and residual confounding may thus have driven results toward the null. Furthermore, the study populations in these studies were much younger than those in the present study, which may have attenuated risk, as our and some other studies suggested a stronger association among older age groups. Stronger associations than those of our study were found in a Dutch cross-sectional study, oversampling elderly participants, and showing to be associated with self-reported heart attack (, 95% CI: 1.12, 1.24) and in a Korean cohort dominated by younger people, where was associated with higher risk of MI hospitalization (, 95% CI: 1.38, 2.09). The authors of the Korean study hypothesized that the high risk may be attributable to higher levels of or different composition of air pollution in comparison with air pollution levels or composition in Western cities. The small size of UFP allows systemic penetration beyond the respiratory tract and provides a large reactive surface area. UFPs may therefore be a more potent risk factor for cardiovascular disease than , of which they constitute a weight-wise small fraction. The present study showed an association between UFP and higher risk for MI. The few previous studies on long-term exposure to UFP and risk of MI have suggested a positive association with MI, which is corroborated by results on self-reported stroke/ischemic heart disease and ischemic heart disease mortality. These results support our UFP and MI findings, although they tended to show weaker associations than those of our study, which may reflect differences in age composition, available covariates, or exposure composition and assessment. The physical, chemical, and toxicological characteristics of PM depend on their sources, which may have implications for development of cost-efficient prevention strategies. Few other studies have investigated source-specific air pollution in relation to MI. A case–control study from Massachusetts, in the United States, found that the strongest risk factor per IQR for acute MI hospitalization was from area-level sources, whereas there was only a weak association with local traffic-dominated sources, which corroborates our observation for . In contrast, in the Californian Teachers Study, and UFP from gasoline and diesel combustion were associated with similar HRs per IQR for ischemic heart disease as total and UFP. These studies had smaller populations and more crude exposure assessment than the present study. In the present study, nontraffic sources of and UFP seemed most important when expressed per IQR. Altogether, and UFP from both traffic and nontraffic sources seem associated with risk of MI, with the highest HRs per IQR for nontraffic sources. This finding, combined with nontraffic sources contributing most to exposure, suggests that nontraffic sources may contribute more to the public health burden from MI than national traffic sources. The main nontraffic sources of emissions of and UFP air pollution in Denmark are nonindustrial combustion plants (such as wood stoves) and ships (Figure S7; Table S12). Further, the categorical analyses and to some extent the AIC analyses indicated that some of the linear HRs should be treated with caution due to a possible nonlinear shape of the function. The present study found an association between and higher risk of MI, which has also been found in several previous studies. We found associations for from other sources than traffic but virtually no association with hailing from traffic sources, indicating that correlated air pollutants rather than per se might have influenced these results for . We found indicative associations between total EC and traffic EC and risk of MI, but no association with EC from nontraffic sources. One previous study found no association between BC and ischemic heart disease, another study found BC to be associated with coronary heart disease, and a third study found a higher incidence of coronary events in association with absorbance although not statistically significant. Altogether, the literature is inconclusive concerning a possible association between EC and MI.

The prospective nationwide cohort design was a major strength of our study. Furthermore, we benefited from extensive and continuously updated high-quality public registers providing detailed information on both outcome, residential history, and a wide array of covariates. Similar to two recent state or nationwide studies with millions of participants, our study could only discern a positive relationship between air pollution and MI after adjusting for individual- and area-level factors. This phenomenon could be partly due to higher air pollution levels at the residences of individuals of higher socioeconomic position in Denmark, combined with higher risk of MI among individuals of low socioeconomic position. The association of socioeconomic position and risk of MI is likely through lifestyle. We could not adjust our analyses for individual lifestyle factors, such as smoking, diet, and physical activity, which could potentially have confounded our results. In a previous study, we analyzed the consequences of lack of information about such lifestyle factors among 246,766 randomly selected residents in Denmark for whom questionnaire data on lifestyle were also available. We found that when including all individual- and area-level covariates of the present study, information on body mass index, smoking, physical activity, and diet bestowed no additional impact on HRs for air pollution and MI. Even so, we cannot entirely rule out that residual confounding from lifestyle may have affected our results in either direction. Air pollution is a complex mixture of many chemical species that might correlate depending on the source pattern. Thus we cannot exclude the possibility that associations with MI, which was identified with one pollutant, might be (partly) due to exposure to correlated pollutants. Two other qualities of our study was the state-of-the-art, address-level exposure modeling covering an entire country and the ability of the modeling to assess source-specific contributions to air pollution. The validity of separation of the traffic and nontraffic contributions to air pollution depends to a high degree on the quality of the emission inventory for road traffic. The uncertainty of the Danish emission inventory depends on the emission sector and the pollutant considered. For some sectors, e.g., large point sources such as power plants and road transport, the data quality is high because a large share of the emission is covered by direct measurements or very detailed information on technology and abatement. For other sectors, such as small-scale combustion, the uncertainty is higher due to a large number of small installations with limited knowledge about technology and abatement measures. Regarding pollutants, the uncertainty for is lower than for PM and EC because of a higher data quality of main emitters of than of PM. Greater uncertainty may also have affected UFP estimates because there are no inventory standards and no mandatory national registration. However, correlations with measurement data were similar to those for other pollutants. Direct validation against measurements was not possible for source-specific concentrations. Circumstantial evidence for their quality is, however, provided by our model’s ability to predict concentrations also in conditions where traffic is the primary source of air pollution and by the quality of input data inventories. Our modeling system only allowed separate estimation of primarily emitted national road-traffic pollution, and we can thus not determine to what degree long-distance transported secondary species related to traffic contribute to MI risk. It was a particular strength that our modeling included also UFP, which has not previously been modeled in detail on such a large scale. A limitation of our study is exposure misclassification. First, we did not have information about occupational or other nonresidential exposures. Second, some uncertainty is inevitable when modeling exposure. We would expect the resulting error in the exposure assessment to be a mixture of Berkson and Classical error, which could affect risk estimate in either direction. We recommend caution when generalizing the results of the present nationwide study to other populations. Potential caveats include differences in the age composition, levels and composition of air pollution, population susceptibility (e.g., due to genetics and comorbidity), and possible differences in correlation between air pollution and sociodemographic characteristics at the individual and area level.

This study adds to the evidence that PM air pollution, including UFP, is associated with higher risk of MI. The association was apparent for PM from both traffic and nontraffic sources, with highest HRs per IQR for nontraffic sources, which also contributed most of the air pollution.

Click here for additional data file. Click here for additional data file.