The effect of secondary inorganic aerosols, soot and the geographical origin of air mass on acute myocardial infarction hospitalisations in Gothenburg, Sweden during 1985–2010: a case-crossover study
Environmental Health volume 13, Article number: 61 (2014)
The relative importance of different sources of air pollution for cardiovascular disease is unclear. The aims were to compare the associations between acute myocardial infarction (AMI) hospitalisations in Gothenburg, Sweden and 1) the long-range transported (LRT) particle fraction, 2) the remaining particle fraction, 3) geographical air mass origin, and 4) influence of local dispersion during 1985–2010.
A case-crossover design was applied using lag0 (the exposure the same day as hospitalisation), lag1 (exposure one day prior hospitalisation) and 2-day cumulative average exposure (CA2) (mean of lag0 and lag1). The LRT fractions included PMion (sum of sulphate, nitrate and ammonium) and soot measured at a rural site. The difference between urban PM10 (particulate matter with an aerodynamic diameter smaller than 10 μm) and rural PMion was a proxy for locally generated PM10 (PMrest). The daily geographical origin of air mass was estimated as well as days with limited or effective local dispersion. The entire year was considered, as well as warm and cold periods, and different time periods.
In total 28 215 AMI hospitalisations occurred during 26 years. PM10, PMion, PMrest and soot did not influence AMI for the entire year. In the cold period, the association was somewhat stronger for PMrest than for urban PM10; the strongest associations were observed during 1990–2000 between AMI and CA2 of PMrest (6.6% per inter-quartile range (IQR), 95% confidence interval 2.1 to 11.4%) and PM10 (4.1%, 95% CI 0.2% − 8.2%). Regarding the geographical air mass origins there were few associations. Days with limited local dispersion showed an association with AMI in the cold period of 2001–2010 (6.7%, 95% CI 0.0% − 13.0%).
In the cold period, locally generated PM and days with limited local dispersion affected AMI hospitalisations, indicating importance of local emissions from e.g. traffic.
Cardiovascular disease (CVD) is the number one cause of death globally and also in European countries, such as Sweden . Ischemic heart disease (IHD) is responsible for 61% of the 2 million CVD deaths in Europe annually, and stroke accounts for most of the rest . Acute myocardial infarction (AMI) is the most important manifestation of ischemic heart disease.
The major risk factors for CVD are age, sex, smoking, low physical activity, increased waist circumference, diabetes and hypertension [2, 3]. Evidence is increasing on the effects of non-traditional risk factors, such as air pollution and weather (e.g. temperature) on CVD mortality and morbidity, specifically as short-term risk factors [4–6]. Nevertheless, many research gaps remain .
Numerous short-term epidemiological studies as well as in-vivo and in-vitro toxicological studies reported on the toxicity of the physical and chemical properties of particulate matter, but effects seem not to be limited to single properties or sources [4, 8, 9].
Polluted air may be local in origin or transported over long distances and the composition may vary . The composition of air masses from Eastern Europe, where coal heating continues to dominate, differs from that of Western Europe . Emissions have also dramatically changed over the past 20–30 years, as in the case of sulphur dioxide. The geographical origin of the air mass can be traced through a combination of measurements and calculations using meteorological models (trajectory analysis) . Large differences in concentrations may occur; e.g. with respect to trace elements and particulate matter with an aerodynamic diameter smaller than 2.5 μm (PM2.5) and black smoke [13–15]. Transition metals and polyaromatic hydrocarbons are possible causes of oxidative stress and inflammation . Trajectory analyses have been used only sporadically for studies on health effects of air pollutants and not previously in a case-crossover design .
Another research gap is that few epidemiological studies focused on subtypes of CVD, such as AMI hospitalisations [4–6]. It is assumed that the pathophysiology of subtypes of CVD is different and therefore different associations between air pollution and these subtypes can be anticipated.
The REVIHAAP report concluded that there is currently no threshold level (i.e. no safe level) for PM10, PM2.5, nitrogen dioxide (NO2) and ground-level ozone (O3) and that the concentration-response functions are mostly linear [7, 18]. Hence it is appropriate to conduct studies in cities with ‘relatively low levels’ like Gothenburg.
The objective of this study was to investigate whether there is an association between AMI hospitalisations in Gothenburg, Sweden and surrogates for long-range transported (LRT) air pollutants (ions and soot) and the more locally emitted parts of PM10. Another aim was to determine the association between AMI hospitalisations and geographical air mass origin (i.e. a surrogate for LRT air pollution and its composition). Gothenburg is a suitable area for such studies, since this region is characterised by a rather large meteorological variability resulting in air mass transport both from clean and polluted areas [14, 15, 19]. A case-crossover study design was applied to the study period of 26 years (1985–2010).
Population and register data
Sweden has a publicly financed health care system with hospital care available to all citizens at low cost. Swedish hospitals register principal and contributory discharge diagnoses for all patients in the national hospital discharge register. In the present study we used data from this register with coverage of all hospitalisations in Gothenburg since 1970.
The study population includes all cases of AMI events registered in Gothenburg from 1 January 1985 to 31 December 2010. The study period was determined by the availability of air pollution data. The International Classification of Diseases (ICD) versions 8 and 9 (ICD-8, ICD-9) were used from 1985 to 1996 and version 10 (ICD 10) from 1997 and onwards.
Hospitalisation for an AMI was defined as a discharge (dead or alive) with a principal diagnosis of ICD 8–9 410 (until 1996) or ICD 10 I21. Only emergency hospitalisations to any of the five hospitals within Gothenburg (1.7 to 5.8 km from the urban background monitoring station) were included. AMI hospitalisations that occurred within 28 days after a previous AMI hospitalisation were excluded (3 068 admissions) as readmissions following discharge for AMI are quite high . Due to the design of the case-crossover study, overlap between being a case and a control may occur for these 3 068 admissions. For information regarding the case-crossover study, see statistical analysis section.
Registry based health outcome data were applied in this study and all identifying variables were deleted. The study was approved by the institutional review board of the University of Gothenburg which waived the need for written informed consent.
Long-range transported air pollution data
PM ion , sulphate, nitrate and ammonium ions and soot
PMion (estimated sum of sulphate (SO4 2−), nitrate (NO3 −) and ammonium (NH4 +) ions) and the three individual ions were used as surrogates for LRT secondary PM .
where the subscripts S and N denote that the mass has been given as the equivalent mass of sulphur or nitrogen, respectively. In order to treat the particulate concentration variables in a comparable manner, these values were converted to equivalent masses of the SO4 2−, NO3 − and NH4 + ions, using the conversion factors 3.0, 4.4 and 1.3, respectively . The three ions were measured by IVL Swedish Environmental Research Institute on behalf of the Swedish Environmental Protection Agency, within the national air quality monitoring network at a rural background site located in Råö, 45 km south of the city centre [19, 22]. At the same place soot was monitored by a reflectance method .
Geographical origin of air masses
The geographical origins of the air masses that pass Gothenburg were applied as another surrogate for the LRT fraction of air pollution and its composition. For all 9 496 days two-dimensional trajectories were calculated and further analysed, using the HYSPLIT model to show the geographical origin of the air masses for these days.
The trajectories were divided into six classes which define the geographical origin of the air mass (Southern Scandinavia, Northern Scandinavia, Baltic Sea, Eastern Europe, United Kingdom (UK)/Denmark (DK)/North Sea and North Atlantic, Additional file 1), based on large scale daily synoptic weather data from NCEP/NCAR (National Centre for Environmental Prediction/National Centre for Atmospheric Research [19, 24, 25].
Days were also classified as either limited local dispersion days or effective local dispersion days. Limited local dispersion day means a day with inversion and/or stable air stratification suppressing the dispersion of locally emitted air pollutants. On the limited local dispersion days the concentration contribution from LRT is almost unchanged. However, the local contribution will be higher. We applied the Lambs weather type classification , along with the concentration of PMion as a surrogate for LRT  and nitrogen oxides (NOx) as a surrogate for limited urban dispersion  for the classification.
Meteorological and urban air pollution data
Temperature, relative humidity (HMP45a probe), PM10 and PM2.5 (tapered element oscillating microbalance instrument) were measured at the urban background monitoring station by the Environment Office in Gothenburg. NO2, NOx and O3 were also measured at the same station .
Daily averages were applied in the statistical analyses. Missing values were not imputed. For 437 days the PM10 levels were lower than 5 μg.m−3 and were set as 5 μg.m−3. The difference between the urban PM10 and the rural PMion levels was used as a surrogate for the more locally generated PM fraction (hereafter PMrest). PM10 and PM2.5 data are available only from 1990 and 2006, respectively.
The time-stratified case-crossover design was applied to investigate the association between AMI hospitalisations and PM10, PM2.5, surrogates for LRT air pollutants (PMion fraction, three individual ions, soot, geographical origin of the air masses), PMrest fraction and limited/effective local dispersion. The case-crossover design was developed as a variant of the case–control design to study the effects of transient exposures on emergency events, comparing each person’s exposure in a time period just prior to a case-defining event with the person’s exposure at other times . A time-stratified approach was applied to select the control days, defining the day of hospitalisation as the case day and the same day of the other weeks in the same month and year as control days. With this approach even very strong confounding of exposure by seasonal patterns is controlled by design [28–30]. The data were analysed using conditional logistic regression analysis (PROC PHREG in SAS 9.2, SAS Institute, Cary, NC).
Lag0 (the exposure the same day as hospitalisation) and lag1 (exposure one day prior hospitalisation) were investigated, as well as the 2-day cumulative average (CA2, mean of lag0-1). The values of CA2 were set as missing if any of the values of lag0 or lag1 were missing. Control days for lag1 and CA2 were defined as for lag0.
All models were adjusted for public holidays (as a binary variable), temperature and relative humidity (same lag as PM10, PMion, PMrest, the three individual ions, PM2.5, soot, geographical origin of air mass and limited/effective local dispersion). PM10, PMion, PMrest, the three individual ions, PM2.5, soot and meteorological variables were included as linear terms in the models. Previous studies reported a linear relationship between PM10 and AMI hospitalisations [31–34]. The linearity of the relationship between AMI hospitalisations and the meteorological variables (temperature and relative humidity) was confirmed in another study .
Models that included the geographical origin of the air mass were run for lag0 and lag1, but not for CA2, as the geographical origin of the air mass is a nominal variable. Southern Scandinavia was selected as the reference category as this is where Gothenburg is located, although such air masses could also come from adjacent countries (Additional file 1) . Models that included limited/effective local dispersion days were run in the same way.
Although individual factors (e.g. sex, age, season) cannot be examined due to the nature of the case-crossover design where each person is his/her own control, the susceptibility to a pollutant may differ between subgroups (sex, age, season). Warm and cold months have different temperatures, which many lead to different time-activity behaviour and also a difference in the composition of the air pollution mix due to fewer emissions from heating. Hence, association was investigated in stratified analyses, where the effect of the pollutant was compared between the spring/summer half-year (April − September) (hereafter warm period) and the autumn/winter half-year (October − March) (hereafter cold period). This approach was applied in previous studies [35–37].
A further stratification was made into time periods by 1985(1990) − 2000 and 2001–2010. Data were available from 1985 for soot, ions and geographical air mass origins as well as for days with limited or effective local dispersion. As mentioned before PM10 data were only available from 1990. A reason for this stratification is that air pollution emissions changed in Europe during the last two decades, with large reductions in SOx, NOx and NHx emissions during the late 1990s . The emission changes might have influenced the composition and toxicology of air pollution, such as PM, in the atmosphere. The earlier time period (1985/1990 − 2000) cannot be investigated for PM2.5 therefore we did not investigate this second stratification for PM2.5. As mentioned before PM2.5 data were only available from 2006. To test for different effects of the pollutants in different subgroups, models with interaction terms were applied.
For the PM10, PMion, PMrest, the three individual ions, PM2.5 and soot, odds ratios (OR) and the 95% confidence intervals (CI) were calculated per inter-quartile range (IQR) increase, which provide the magnitude-of-risk estimates that are comparable across the exposure variables. The results are presented as the per cent excess risk in AMI hospitalisations per IQR increase in the air pollution levels using the following calculation: (exp(βxIQR) – 1) × 100%, where β is the model estimate. For analysis of a given lagged exposure, a case was automatically dropped if exposure data were not available for the case and at least one control day.
Table 1 indicates the characteristics of the 28 215 AMI hospitalisations (of 22 475 people) during the study period. The mean age for AMI hospitalisations was 74 years. More hospitalisations occurred among men, people older than 75 years and during the cold period. Additional file 2 illustrates the time-series of the daily number of AMI hospitalisations during the study period. The number of AMI hospitalisations varied from 0 to 23 per day.
Table 2 provides an overview of the daily levels of PM10, PMion, PMrest, the three individual ions, PM2.5, soot, temperature and relative humidity. The daily WHO and EU air quality limits (50 μg.m−3) for the PM10 levels were exceeded on 60 days at the urban background site during the 21-year study period (Additional file 3) .
Additional files 4 and 5 illustrate the PM10, PMion, PMrest, the three individual ions, PM2.5, soot, temperature and relative humidity and additionally the NOx, NO2, O3 and wind speed levels by the six geographical air mass origins. Descriptive results of NOx, NO2 and O3 were reported in another publication . Compared to Southern Scandinavia, NO2 and NOx levels were higher when the air mass originated from Northern Scandinavia. In contrast, PM10, PMion, the three individual ions, PM2.5, soot, temperature and relative humidity were lower. Higher PM10, PMion, sulphate ion, ammonium ion, PM2.5, soot and relative humidity but lower NO2 and NOx levels were observed when the air mass origin was from Eastern Europe compared to Southern Scandinavia. The geographical air mass origin variable is thus not just a surrogate for LRT air pollution and its composition, but a combination of both the large-scale meteorological conditions as well as the large-scale air pollution transportation pattern.
Table 3 displays the Spearman correlations between PM10, PMion, PMrest, the three individual ions, PM2.5, soot, temperature and relative humidity during the entire year, and the warm and cold periods. Relative humidity was in general inversely correlated with the other exposure variables, except with the three individual ions and PMion. Temperature had a weak positive correlation with the other exposure variables, except with PMrest in the entire year and the cold period. PM10 had a positive and stronger correlation with PMrest and PM2.5 than with PMion. PMion had a stronger correlation with PM2.5 than with PM10. PMion had a weak inverse correlation with PMrest in the entire year and the cold period. The three individual ions were strongly correlated. Relatively high correlations were found between soot and sulphate and nitrate ions.
Figure 1 illustrates the% change in the AMI hospitalisations per IQR increase in the lags of PM10, PMion, PMrest, PM2.5 and soot during the entire year, warm and cold periods. The point estimates for PM10 (CA2) and PM2.5 (CA2) for the entire year were similar, 1.4% and 1.5% per IQR, but the 95% CI included zero. The 95% CI intervals of PM2.5 are wider due to the shorter study period (2006–2010). The point estimate for the increase in AMI in the cold period was somewhat larger for PMrest (3.2%, 95% CI −0.3% − 6.8%) than for PM10 (2.6%, 95% CI −0.4% − 5.7%). When the cold period was restricted to 1990–2000, stronger positive associations were observed between AMI hospitalisations and PMrest (6.6%, 95% CI 2.1% − 11.4%) and PM10 (4.1%, 95% CI 0.2% − 8.2%) (Figure 2).
Although there was no association between AMI and PMion, a positive association with sulphate ion was observed in the warm period (1.9%, 95% CI −0.4% − 4.3%) (Figure 2, Additional file 6), which was stronger when the warm period was restricted to 1985–2000 (Figure 2). In the latter restricted period, the sample sizes were however much smaller than those in Figure 1. Caution should however be applied to the stratified results as the interaction terms were not statistically significant.
Figure 3 indicates the% change in the AMI hospitalisations for the five geographical air mass origins relative to that of Southern Scandinavia during the entire year, warm and cold periods. There was no clear association between AMI hospitalisations and lag0 of the geographical air mass origin when the entire year was considered. In the unadjusted model, a 5.1% (95% CI 0.1% − 10.3%) increase in AMI was observed when the air mass originated the day before admission (lag1) in Northern Scandinavia (relative to Southern Scandinavia) (Additional file 7). Adjusting further for temperature, relative humidity and public holidays resulted in a weaker association (4.1%, 95% CI −1.4% − 10.0%, entire year, lag1), but with somewhat higher point estimates in the cold period (Figure 3). When the cold period was analysed separately for different time periods, stronger positive associations were observed for 2001–2010 when the air mass originated in Northern Scandinavia (12.6% increase in AMI for lag1), UK/DK/North Sea (16.0% for lag0) and North Atlantic (13.4% significant increase at lag0, 15.1% at lag1) relative to Southern Scandinavia (Additional file 8). The sample sizes were however much smaller than those in Figure 3. Again, caution should be applied to the stratified results as the interaction terms were not significant.
Additional file 9 illustrates the PM10 and PMrest levels by the six geographical air mass origins, after stratification by the entire year, warm and cold periods and also by the time periods 1990–2010, 1990–2000, and 2001–2010. The PM10 and PMrest levels were highest during the 2001–2010 period (Additional file 10).
Finally we examined the influence of air pollution dispersion. Days with limited local dispersion had a strong association at lag0 with AMI hospitalisations during 2001–2010 in the entire year (5.4%, 95% CI 0.7% − 9.8%), but not for lag1 (Figure 4). Caution should be applied to the stratified results as the interaction terms were not significant. For the entire year, warm and cold periods, days with local limited dispersion were much more common when the air masses originated from Eastern Europe, UK/DK/North sea and North Atlantic (<60% of the days) (Table 4).
We set out to address some of the research gaps on the association between CVD and air pollution [5, 7, 34]. First, we investigated whether there is an association between air pollution and AMI emergency room hospitalisations (a specific subtype of CVD). Secondly, we considered surrogates for LRT air pollutants (soot, ammonium, nitrate and sulphate ions), more locally emitted PM10 (PMrest), days with limited dispersion as well as the geographical origins of air masses that pass Gothenburg. Thirdly, the entire year as well as the warm (April − September) and cold periods (October − March) were investigated. Fourthly, we considered different time periods 1985(1990) − 2000, 2001–2010 and 1985(1990) − 2010.
We observed conflicting results with regards to seasonal and time periods, but found that locally generated PM was more important for AMI hospitalisations.
An IQR (10 μg.m−3) increase in lag0 of PM10 was associated with an insignificant increase in AMI hospitalisation (1%) when the entire 26-year period was considered. The strength of the association is similar to that of a meta-analysis: 0.6% significant increase in AMI hospitalisations per 10 μg.m−3 increase in lag0 of PM10. In contrast to the warm period, where no obvious association was observed, an IQR (10 μg.m−3) increase in lag0 of PM10 tended to increase AMI hospitalisations by 2% in the cold period, with borderline significance. The association between AMI hospitalisations and PM10 was strongest in the 1990–2000 period. This is supported by the association between the local fraction of PM10 (PMrest) and AMI hospitalisations in the cold period (2% increase). Reasons may include the possible differences in chemical composition and toxicity of PM10 in the warm and cold periods and during different time periods. In spring and summer dry road surfaces may increase the fraction of road dust in PM10, which may in turn be less toxic than combustion-related PM, which should be more common in winter-time. The mean PM10 level in our study (entire year average about 15 μg/m3) is much lower than PM10 in the most influential studies in the meta-analysis. It is, however, similar to recent population levels of PM10 in the Scandinavian countries, as shown in the ESCAPE study .
The strength of the association between AMI hospitalisations and PM2.5 (lag0) (1.7% per 10 μg.m−3 increase) is insignificant and somewhat weaker than that of a meta-analysis of 13 studies: 2.5% significant increase in AMI hospitalisations per 10 μg.m−3 increase in lag0 of PM2.5. No significant associations were observed when stratified in the warm and cold periods.
Regarding rural background soot levels, a surrogate for LRT air pollution, no clear association was observed with AMI hospitalisation during any seasonal or year periods. A study from Germany (1999–2003) reported a significant increase in AMI hospitalisations with increasing levels of estimated personal soot exposure (relative risk 1.3 per 1.1 m−1 × 10−5) , but not with ambient soot. There might be an association between urban soot and AMI hospitalisations, but data on urban soot levels are not available for Gothenburg.
A recent WHO report discussed the current toxicological and epidemiological evidence regarding secondary aerosols . The expert group proposed that inorganic secondary aerosols, such as ammonium, sulphates and nitrates, are valuable additional air quality metric to evaluate the health risks. We found, however, no associations for total nitrates or total ammonium ions. We observed an increase of 3% in AMI hospitalisations in the warm period during 1985–2000 per IQR increase in sulphate ions (lag0). One of the few epidemiological studies reported a borderline increase of 0.6% and 2.0% in CVD emergency room hospitalisations with an increase in sulphate (in PM2.5 fraction at urban sites, lag0) in the warm (April − September) and cold (October − March) periods, respectively . If sulphate ions are a surrogate for the LRT PM, then the association between sulphate ions and AMI hospitalisations in the warm period (1985–2000) is in contrast with those observed for PM10 and PMrest in the cold period (1990–2000). Hence, the local fraction of PM may be more important for AMI hospitalisations in the cold period, whilst the LTR fraction is of importance in the warm period.
The influence of geographical air mass origin on human health effects has been investigated sporadically, but never in a case-crossover design over a long study period . We compared the likelihood of AMI hospitalisations for six different trajectories using Southern Scandinavia as the reference. The geographical air mass origin variable is not just a surrogate for LRT air pollution and its composition, but a combination of both the large-scale meteorological conditions as well as the large-scale air pollution transportation pattern. PM levels, including also sulphates and soot, were higher when air masses originated from Eastern Europe. The fact that NO2 and NOx were higher when the air masses originated from Northern Europe could be due to lower wind speeds favouring the influence of locally emitted NOx.
We observed the largest increase in AMI hospitalisations at lag1 for air masses from Northern Scandinavia, UK/DK/North Sea and North Atlantic during the 2001–2010 in the cold period, although there was no significant difference between the seasons. This is in conflict with the association between PM10 and PMrest, and AMI hospitalisations, where the strongest associations were observed in the cold period in the time period (1990–2000) and not during 2001–2010, and we cannot provide any logical explanation for this.
In contrast no clear increase in AMI hospitalisations was observed when the air mass originated in Eastern Europe, where coal heating continues to dominate . The adverse association at lag1 when the air mass originated in Northern Scandinavia (more days with effective local dispersion) might be caused by LRT of wood smoke particles in the cold season from rural areas. Both PM10 and PMrest average levels were much higher in 2001–2010 than in 1990–2000 for air masses originated from Northern Scandinavia which can be caused by the increased use of wood burning for heating in Sweden .
The classification according to local meteorology (e.g. temperature inversions) indicated that days with limited local dispersion had a strong association with AMI hospitalisations at lag0 during 2001–2010 in the cold period. This is in accordance with the positive association in the same periods between air masses (UK/DK/North Sea and North Atlantic) having rather high percentage of days with limited local dispersion and thus more influence from local emissions. This indicates significant contribution of local emissions within Gothenburg.
Advantages of our study include accurate meteorological, air pollution and AMI hospitalisation data. Data from the Swedish Hospital Registry have high validity for a diagnosis of MI . Autopsy rates for persons dying outside hospital from IHD are high for younger, but not older people . Some disease misclassification is possible, but it is unlikely to be related to air pollution levels and the geographical origin of air masses. Our study period of 26 years is longer than the other studies that investigated the association between air pollution levels and AMI, which had study periods of 2–16 years . Unlike our study, most studies do not report the actual number of cases or days included in their case-crossover analyses, due to missing exposure data.
Another advantage of our study is that unlike most of the other studies, we investigated the effect of seasonal and time periods by interaction analyses. We also investigated a more refined PM estimate (PMrest), geographical air mass origins and limited/effective local dispersion days. To our knowledge no similar classifications have been used before in health-related studies.
As with all ecological epidemiology study designs, a disadvantage of our study is a possible exposure misclassification, i.e. the assumption that the ambient temperature, humidity and air pollution measured in the city are the same across Gothenburg. The exposure error resulting from using ambient temperature and air pollution as a surrogate for personal exposure can potentially lead to bias in the estimated association, and this can be more pronounced among the elderly and other frail groups who generally spend most of their time indoors. However, in a case-crossover study investigating the risk of AMI, more refined exposure PM2.5 surrogates that account for human activities did not result in larger health effects compared with analyses that used a central monitoring site .
A second limitation is the inability to consider the chemical composition of PM10 and PM2.5 and the relatively short measurement period of PM2.5 (only 5 years) .
A third limitation is that of the uncertainties of trajectory models, which have been discussed in many studies [46, 47]. The sources of uncertainty for a single trajectory can be the meteorological field, starting height, trajectory calculation or trajectory clustering. In this study, an ensemble method was carried out to diminish the uncertainty (as regards starting height, trajectory calculation) of a single trajectory . However, the uncertainties from other sources cannot be avoided.
A fourth limitation is that information on effect modifiers, e.g. the use of medications, pre-existing CVD or comorbidities, was not available in our study [48–50]. Such effect modifiers may bias the association between the air pollutants and AMI hospital admissions in either direction.
We performed a large number of tests which increases the probability of obtaining spurious significant associations. When multiple hypotheses are tested, the risk of false positive conclusions can increase and be higher than the intended level of e.g. 0.05. There are various ways of controlling the family-wise error rate (the risk of at least one false positive), by e.g. Bonferroni correction, the Holm step-down procedure or splitting and recycling the test mass [51, 52]. However, our study is very exploratory and can be considered as hypothesis-generating where one aim is to find patterns for how the composition of air masses from different geographical origins is associated with AMI. The p-values can be considered descriptive rather than decisive and the results can be used as a guide for a future study that can focus on the cold season − where this study showed the strongest results.
AMI hospitalisations may not accurately reflect AMI incidence. Fatal coronary heart disease events that occurred out-of-hospital could have been included in our analyses. Mortality in AMI has decreased much more (by about 50% in the past 15 years) than incidence of AMI (about 16%) . Much of the decrease in mortality is, however, caused by better in-hospital treatment, and nearly one in three patients with a major coronary heart disease event dies without reaching a hospital. Nevertheless, AMI hospitalisations are a reasonable surrogate for major coronary heart disease events.
Locally generated PM (PMrest) is associated with an increase in AMI hospitalisations during 1990–2000 in the cold period, whilst days with limited local dispersion within Gothenburg were important during 2001–2010, also in the cold period. Altogether this indicates the significance of local emissions from e.g. traffic within Gothenburg.
Acute myocardial infarction
International classification of diseases
Ischemic heart disease
- NOx :
- NO2 :
- O3 :
- PMion :
Sum of sulphate, nitrate and ammonium
- PMrest :
Difference between urban PM10 and rural PMion was a surrogate for locally generated PM10
- PM2.5 :
Particulate matter with an aerodynamic diameter smaller than 2.5 μm
- PM10 :
Particulate matter with an aerodynamic diameter smaller than 10 μm
- 95% CI:
95% confidence interval
World Health Organisation.
Rayner M, Allender S, Scarborough P: Cardiovascular disease in Europe. Eur J Cardiovasc Prev Rehabil. 2009, 16: S43-S47.
Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, Lanas F, McQueen M, Budaj A, Pais P, Varigos J, Lisheng L, INTERHEART Study Investigators: Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case–control study. Lancet. 2004, 364: 937-952.
Lloyd-Jones DM, Leip EP, Larson MG, D’Agostino RB, Beiser A, Wilson PW, Wolf PA, Levy D: Prediction of lifetime risk for cardiovascular disease by risk factor burden at 50 years of age. Circulation. 2006, 113: 791-798.
Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, Diez-Roux AV, Holguin F, Hong Y, Luepker RV, Mittleman MA, Peters A, Siscovick D, Smith SC, Whitsel L, Kaufman JD, American Heart Association Council on Epidemiology and Prevention, Council on the Kidney in Cardiovascular Disease, and Council on Nutrition, Physical Activity and Metabolism: Particulate matter air pollution and cardiovascular disease: An update to the scientific statement from the American Heart Association. Circulation. 2010, 121: 2331-2378.
Mustafic H, Jabre P, Caussin C, Murad MH, Escolano S, Tafflet M, Périer MC, Marijon E, Vernerey D, Empana JP, Jouven X: Main air pollutants and myocardial infarction: a systematic review and meta-analysis. JAMA. 2012, 307: 713-721.
Turner LR, Barnett AG, Connell D, Tong S: Ambient temperature and cardiorespiratory morbidity: a systematic review and meta-analysis. Epidemiology. 2012, 23: 594-606.
World Health Organisation: Review of evidance on health aspects of air pollution- REVIHAAP Project 2013. http://www.euro.who.int/en/health-topics/environment-and-health/air-quality/publications/2013/review-of-evidence-on-health-aspects-of-air-pollution-revihaap-project-final-technical-report. Accessed 14 January 2014
World Health Organisation: Air Quality Guideline Global Update 2005. 2006, http://www.who.int/phe/health_topics/outdoorair_aqg/en/. Accessed 14 January 2014
Rückerl R, Schneider A, Breitner S, Cyrys J, Peters A: Health effects of particulate air pollution: A review of epidemiological evidence. Inhal Toxicol. 2011, 23: 555-592.
Antilla P, Frolova M, Heidam N, Lövblag G, Sjöberg K, Pajuste K, Schaug J, Sopauskiene D: The Nordic-Baltic Regional Assessment of Long-range Transboundary Air Pollution 1980–2000. Nordic Council of Ministers. 2003, Report AIP, 418-http://www.norden.org/en/publications/publikationer/2003-418. Accessed 14 January 2014
Amann M, Bertok I, Borken-Kleefeld J, Cofala J, Heyes C, Höglund-Isaksson L, Klimont Z, Rafaj P, Schöpp W, Wagner F: Cost-effective Emission Reductions to Improve Air Quality in Europe in 2020. Scenarios for the Negotiations on the Revision of the Gothenburg Protocol under the Convention on Long-Range Transboundary Air Pollution. 2011, http://www.unece.org/fileadmin/DAM/env/documents/2011/eb/wg5/WGSR48/Informal%20docs/Info.doc.8_CIAM_report_on_Cost_effective_emission_reductions_to_improve_air_quality_in_Europe_in_2010.pdf. Accessed 14 January 2014
Gustafsson M, Rayner D, Chen D: Extreme rainfall events in southern Sweden: where does the moisture come from?. Tellus. 2010, 62A: 605-616.
Molnár P, Gustafson P, Johannesson S, Boman J, Barregård L, Sällsten G: Domestic wood burning and PM2.5 trace elements: Personal exposure, indoor and outdoor levels. Atmos Environ. 2005, 39: 2659-2669.
Molnár P, Johannesson S, Boman J, Barregard L, Sällsten G: Personal exposures, and indoor, residential outdoor and urban background levels of fine trace elements in the general population. J Environ Monit. 2006, 8: 543-551.
Johannesson S, Gustafson P, Molnár P, Barregard L, Sällsten G: Exposure to fine particles (PM2.5 and PM1) and black smoke in the general population: personal, indoor, and outdoor levels. J Expo Sci Environ Epidemiol. 2007, 17: 613-624.
Schlesinger RB, Kunzli N, Hidy GM, Gotschi T, Jerrett M: The health relevance of ambient particulate matter characteristics: coherens of toxicological and epidemiological inferences. Inhal Toxicol. 2006, 18: 95-125.
Park SK, O’Neill MS, Stunder BJB, Vokonas PS, Sparrow D, Koutrakis P, Schwartz J: Source location of air pollution and cardiac autonomic function: Trajectory cluster analysis for exposure assessment. JESEE. 2007, 17: 488-497.
Cesaroni G, Forastiere F, Stafoggia M, Andersen ZJ, Badaloni C, Beelen R, Caracciolo B, de Faire U, Erbel R, Eriksen KT, Fratiglioni L, Galassi C, Hampel R, Heier M, Hennig F, Hilding A, Hoffmann B, Houthuijs D, Jöckel KH, Korek M, Lanki T, Leander K, Magnusson PK, Migliore E, Ostenson CG, Overvad K, Pedersen NL, JP J, Penell J, Pershagen G: 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-doi:10.1136/bmj.f7412
Tang L, Haeger-Eugensson M, Sjöberg K, Wichmann J, Molnár P, Sallsten G: Estimation of the long-range transport contribution from secondary inorganic components to urban background PM10 concentrations in south-western Sweden during 1986–2010. Atmos Environ. 2014, 89: 93-101.
Kociol RD, Lopes RD, Clare R, Thomas L, Mehta RH, Kaul P, Pieper KS, Hochman JS, Weaver WD, Armstrong PW, Granger CB, Patel MR: International variation in and factors associated with hospital readmission after myocardial infarction. JAMA. 2012, 307: 66-74.
Kukkonen J, Sokhi R, Luhana L, Härkönen J, Salmi T, Sofiev M, Karppinen A: Evaluation and application of a statistical model for assessment of long-range transported proportion of PM2.5 in the United Kingdom and in Finland. Atmos Environ. 2008, 42: 3980-3991.
Ferm M, Hellsten S: Trends in atmospheric ammonia and particulate ammonium concentrations in Sweden and its causes. Atmos Environ. 2012, 61: 30-39.
Organisation for Economic Co-operation and Development: Methods of measuring air pollution. 1964, Paris: OECD
Chen D: A monthly circulation climatology for Sweden and its application to a winter temperature case study. Int J Climatol. 2000, 20: 1067-1076.
Haeger-Eugensson M: Vertical Interactions in a Nocturnal Multi-scale Wind System Influenced by Atmospheric Stability in a Coastal area. Theor Applied Climatol. 1999, 64: 69-82.
Wichmann J, Rosengren A, Sjöberg K, Barregard L, Sallsten G: Association between Ambient Temperature and Acute Myocardial Infarction Hospitalisations in Gothenburg, Sweden: 1985–2010. PLoS ONE. 2013, 8: e62059-doi:10.1371/journal.pone.0062059
Maclure M: The case-crossover design: a method for studying transient effects on the risk of acute events. Am J Epidemiol. 1991, 133: 144-153.
Bateson TF, Schwartz J: Control for seasonal variation and time trend in case-crossover studies of acute effects of environmental exposures. Epidemiology. 1999, 10: 539-544.
Bateson TF, Schwartz J: Selection bias and confounding in case-crossover analyses of environmental time-series data. Epidemiology. 2001, 12: 654-661.
Levy D, Lumley T, Sheppard L, Kaufman J, Checkoway H: Referent selection in casecrossover analyses of acute health effects of air pollution. Epidemiology. 2001, 12: 186-192.
Anderson HR, Atkinson RW, Bremner SA, Carrington J, Peacock J: Quantitative systematic review of short term associations between ambient air pollution (particulate matter, ozone, nitrogen dioxide, sulphur dioxide and carbon monoxide), and mortality and morbidity. Report to Department of Health revised following first review. 2007, http://www.dh.gov.uk/en/Publicationsandstatistics/Publications/PublicationsPolicyAndGuidance/DH_121200. Accessed 14 January 2014
Abrignani MG, Corrao S, Biondo GB, Renda N, Braschi A, Novo G, Di Girolamo A, Braschi GB, Novo S: Influence of climatic variables on acute myocardial infarction hospital admissions. Int J Cardiol. 2009, 137: 123-129.
Wolf K, Schneider A, Breitner S, von Klot S, Meisinger C, Cyrys J, Hymer H, Wichmann HE, Peters A, Cooperative Health Research in the Region of Augsburg Study Group: Air temperature and the occurrence of myocardial infarction in Augsburg, Germany. Circulation. 2009, 120: 735-742.
Bhaskaran K, Hajat S, Haines A, Herrett E, Wilkinson P, Smeeth L: Effects of air pollution on the incidence of myocardial infarction. Heart. 2009, 95: 1746-1759.
D’Ippoliti D, Forastiere F, Ancona C, Agabiti N, Fusco D, Michelozzi P, Perucci CA: Air pollution and myocardial infarction in Rome: a case-crossover analysis. Epidemiology. 2003, 14: 528-535.
Goldberg MS, Burnett RT, Yale JF, Valois MF, Brook JR: Associations between ambient air pollution and daily mortality among persons with diabetes and cardiovascular disease. Environ Res. 2006, 100: 255-267.
Nawrot TS, Torfs R, Fierens F, De Henauw S, Hoet PH, Van Kersschaever G, De Backer G, Nemery B: Stronger associations between daily mortality and fine particulate air pollution in summer than in winter: Evidence from a heavily polluted region in Western Europe. J Epidemiol Community Health. 2007, 61: 146-149.
Tørseth K, Aas W, Breivik K, Fjæraa AM, Fiebig M, Hjellbrekke AG, Lund Myhre C, Solberg S, Yttri KE: Introduction to the European Monitoring and Evaluation Programme (EMEP) and observed atmospheric composition change during 1972–2009. Atmos Chem Phys. 2012, 12: 5447-5481.
Beelen R, Raaschou-Nielsen O, Stafoggia M, Andersen ZJ, Weinmayr G, Hoffmann B, Wolf K, Samoli E, Fischer P, Nieuwenhuijsen M, Vineis P, Xun WW, Katsouyanni K, Dimakopoulou K, Oudin A, Forsberg B, Modig L, Havulinna AS, Lanki T, Turunen A, Oftedal B, Nystad W, Nafstad P, De Faire U, Pedersen NL, Östenson CG, Fratiglioni L, Penell J, Korek M, Pershagen G: Effects of long-term exposure to air pollution on natural-cause mortality: an analysis of 22 European cohorts within the multicentre ESCAPE project. Lancet. 2014, 383: 785-795.
von Klot S, Cyrys J, Hoek G, Kühnel B, Pitz M, Kuhn U, Kuch B, Meisinger C, Hörmann A, Wichmann HE, Peters A: Estimated personal soot exposure is associated with acute myocardial infarction onset in a case-crossover study. Prog Cardiovasc Dis. 2011, 53: 361-368.
Ito K, Mathes R, Ross Z, Nádas A, Thurston G, Matte T: Fine particulate matter constituents associated with cardiovascular hospitalizations and mortality in New York City. Environ Health Perspect. 2011, 119: 467-473.
Molnár P, Sallsten G: Contribution to PM from domestic wood burning in a small community in Sweden. Environ Sci Process Impacts. 2013, 15: 833-838.
Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim JL, Reuterwall C, Heurgren M, Olausson PO: External review and validation of the Swedish national inpatient register. BMC Public Health. 2011, 11: 450-doi:10.1186/1471-2458-11-450
Dudas K, Lappas G, Stewart S, Rosengren A: Trends in out-of-hospital deaths due to coronary heart disease in Sweden (1991 to 2006). Circulation. 2011, 123: 46-52.
Hodas N, Turpin BJ, Lunden MM, Baxter LK, Özkaynak H, Burke J, Ohman-Strickland P, Thevenet-Morrison K, Kostis JB, for the MIDAS 21 Study group and Rich D: Refined ambient PM2.5 exposure surrogates and the risk of myocardial infarction. JESEE. 2013, 23: 573-580.
Stohl A: Computation, accuracy and applications of trajectories - a review and bibliography. Atmos Environ. 1998, 32: 947-966.
Harris JM, Drexler RR, Oltmans SJ: Trajectory model sensitivity to differences in input data and vertical transport method. J Geophys Res. 2005, 110: D14109-doi:10.1029/2004JD005750
Miller KA, Siscovick DS, Sheppard L, Shepherd K, Sullivan JH, Anderson GL, Kaufman JD: Long-term exposure to air pollution and incidence of cardiovascular events in women. N Engl J Med. 2007, 356: 447-458.
O’Neill MS, Veves A, Sarnat JA, Zanobetti A, Gold DR, Economides PA, Horton ES, Schwartz J: Air pollution and inflammation in type 2 diabetes: a mechanism for susceptibility. Occup Environ Med. 2007, 64: 373-379.
Park SK, O’Neill MS, Wright RO, Hu H, Vokonas PS, Sparrow D, Suh H, Schwartz J: HFE genotype, particulate air pollution, and heart rate variability: a gene-environment interaction. Circulation. 2006, 114: 2798-2805.
Holm S: A simple sequentially rejective multiple test procedure. Scan J Stat. 1979, 6: 65-70.
Burman C-F, Sonesson C, Guilbaud O: A recycling framework for the construction of Bonferroni-based multiple tests. Stat Med. 2009, 28: 739-761.
The National Board of Health and Welfare: Sweden, http://www.socialstyrelsen.se. Accessed 15 April 2014
The authors would like to thank Georgios Lappas (Department of Emergency and Cardiovascular Medicine, Sahlgrenska University Hospital and University of Gothenburg) for data management along with Peter Molnár and Sara Gustavsson (Occupational and Environmental Medicine, Sahlgrenska University Hospital and University of Gothenburg) for discussions regarding the study. The authors gratefully thank the NOAA Air resources Laboratory (ARL) for the providing the HYSPLIT model; the NOAA/OAR/ESRL PSD, Boulder, Colorado, USA, for providing the NCEP Reanalysis datasets and the Environmental Office in Gothenburg for providing air quality data. This study was funded by the Swedish Research Council, FORMAS.
The authors declare that they have no competing interests. The Swedish Research Council Formas funded the study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
JW, EA, LB and GS designed the study. AR cleaned and contributed the acute myocardial infarction hospital admission data. KS, MHE and LT provided and ensured the quality of the air pollution and meteorological data. JW analysed the data. All authors contributed to writing and revising the manuscript and approved the final manuscript.
Janine Wichmann, Karin Sjöberg, Lin Tang, Marie Haeger-Eugensson, Annika Rosengren, Eva M Andersson, Lars Barregard and Gerd Sallsten contributed equally to this work.
Electronic supplementary material
Additional file 2:Time-series of acute myocardial infarction hospitalisation in Gothenburg, Sweden (1 January 1985 − 31 December 2010).(DOCX 555 KB)
Additional file 3:Time-series of available PM 10 urban levels in Gothenburg, Sweden during the study period 1 January 1990 − 31 December 2010.(DOCX 647 KB)
Additional file 4:(a) Average PM 10 *, PM ion , PM rest *, PM 2.5 *, soot, temperature and wind speed and (b) NO x , NO 2 , O 3 and relative humidity levels by origin of air mass in Gothenburg, Sweden (1 January 1985 − 31 December 2010).(DOCX 1 MB)
Additional file 5:(a) Average PM ion , sulphate, total nitrate and total ammonium by origin of air mass in Gothenburg, Sweden (1 January 1985 − 31 December 2010).(DOCX 407 KB)
Additional file 6:Association between the lag0, lag1 and the 2-day cumulative average of sulphate, total nitrate, total ammonium and acute myocardial infarction hospitalisation in Gothenburg (1985−2010) as percentage change in risk (%) and 95% confidence intervals during (a) the entire year, (b) warm period (April − September) and (c) cold period (October − March).(DOCX 2 MB)
Additional file 7:Unadjusted association between the lag0 and lag1 of the origin of the air masses and acute myocardial infarction hospitalisation in Gothenburg during the entire year (1985−2010) as percentage change in risk (%) and 95% confidence intervals.(DOCX 702 KB)
Additional file 8:Association between the lag0 and lag1 of the origin of the air masses and acute myocardial infarction hospitalisations in Gothenburg, Sweden as percentage change in risk (%) and 95% confidence intervals during the cold period (October−March) for (a) 1985−2010, (b) 1985−2000 and (c) 2001−2010.(DOCX 1 MB)
Additional file 9:Average PM 10 and PM rest levels by origin of air mass in Gothenburg, Sweden by (a,d) the entire year, (b,e) warm (April − September) and (c,f) cold periods (October − March). These periods were further stratified by 1990−2010, 1990−2000 and 2001−2010. (DOCX 3 MB)
Additional file 10:Descriptive statistics for daily PM 10 *, PM ion , PM rest *, PM 2.5 * and soot levels (lag0) in Gothenburg, Sweden by seasonal and year periods (1 January 1985−31 December 2010).(DOCX 18 KB)
About this article
Cite this article
Wichmann, J., Sjöberg, K., Tang, L. et al. The effect of secondary inorganic aerosols, soot and the geographical origin of air mass on acute myocardial infarction hospitalisations in Gothenburg, Sweden during 1985–2010: a case-crossover study. Environ Health 13, 61 (2014). https://doi.org/10.1186/1476-069X-13-61
- Air pollution
- Secondary inorganic aerosols
- Geographical air mass origin
- Acute myocardial infarction