Next Article in Journal
Influence of Particulate Matter during Seasonal Smog on Quality of Life and Lung Function in Patients with Chronic Obstructive Pulmonary Disease
Previous Article in Journal
Does Europe Need an EU Product Safety Directive for Access Scaffolding?
Previous Article in Special Issue
Association Between PM2.5 and Daily Hospital Admissions for Heart Failure: A Time-Series Analysis in Beijing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Pollutant Air Quality Health Index (AQHI) Based on Short-Term Respiratory Effects in Stockholm, Sweden

1
Atmospheric Science Unit, Department of Environmental Science and Analytical Chemistry, Stockholm University, 11418 Stockholm, Sweden
2
Environment and Health Administration, SLB, Box 8136, 104 20 Stockholm, Sweden
3
Division of Occupational and Environmental Medicine, Department of Public Health and Clinical Medicine, Umeå University, 90187 Umeå, Sweden
4
Swedish Museum of Natural History, 114 18 Stockholm, Sweden
5
Umeå School of Business, Economics and Statistics (USBE), Umeå University, 90187 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(1), 105; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010105
Submission received: 27 September 2018 / Revised: 19 December 2018 / Accepted: 27 December 2018 / Published: 2 January 2019

Abstract

:
In this study, an Air Quality Health Index (AQHI) for Stockholm is introduced as a tool to capture the combined effects associated with multi-pollutant exposure. Public information regarding the expected health risks associated with current or forecasted concentrations of pollutants and pollen can be very useful for sensitive persons when planning their outdoor activities. For interventions, it can also be important to know the contribution from pollen and the specific air pollutants, judged to cause the risk. The AQHI is based on an epidemiological analysis of asthma emergency department visits (AEDV) and urban background concentrations of NOx, O3, PM10 and birch pollen in Stockholm during 2001–2005. This analysis showed per 10 µg·m–3 increase in the mean of same day and yesterday an increase in AEDV of 0.5% (95% CI: −1.2–2.2), 0.3% (95% CI: −1.4–2.0) and 2.5% (95% CI: 0.3–4.8) for NOx, O3 and PM10, respectively. For birch pollen, the AEDV increased with 0.26% (95% CI: 0.18–0.34) for 10 pollen grains·m–3. In comparison with the coefficients in a meta-analysis, the mean values of the coefficients obtained in Stockholm are smaller. The mean value of the risk increase associated with PM10 is somewhat smaller than the mean value of the meta-coefficient, while for O3, it is less than one fifth of the meta-coefficient. We have not found any meta-coefficient using NOx as an indicator of AEDV, but compared to the mean value associated with NO2, our value of NOx is less than half as large. The AQHI is expressed as the predicted percentage increase in AEDV without any threshold level. When comparing the relative contribution of each pollutant to the total AQHI, based on monthly averages concentrations during the period 2015–2017, there is a tangible pattern. The AQHI increase associated with NOx exhibits a relatively even distribution throughout the year, but with a clear decrease during the summer months due to less traffic. O3 contributes to an increase in AQHI during the spring. For PM10, there is a significant increase during early spring associated with increased suspension of road dust. For birch pollen, there is a remarkable peak during the late spring and early summer during the flowering period. Based on monthly averages, the total AQHI during 2015–2017 varies between 4 and 9%, but with a peak value of almost 16% during the birch pollen season in the spring 2016. Based on daily mean values, the most important risk contribution during the study period is from PM10 with 3.1%, followed by O3 with 2.0%.

1. Introduction

The air quality index (AQI) is a numerical scale that has been used in order to describe the air quality situation in relation to a reference value at a given time and place. An AQI is usually based on several pollutants, where each pollutant is compared with its air quality limit value, and the pollutant with the highest value in relation to its limit value determines the AQI-value [1]. In this case, one individual pollutant determines the index, but without considering the other pollutants with smaller values in relation to their limit values, which means that the cumulative effect is not captured. Another limitation is that the AQI does not capture the no-threshold concentration-response relationship between air pollution exposure and health [2]. The AQI is usually based on pollutant concentrations at a few monitoring stations, often located at hot-spots, but without reflecting the spatial variations in the city, and consequently not the health risks for the general populations associated with air pollution exposure.
The air quality health index (AQHI) is a tool that has been developed in order to reflect the health risks associated with simultaneous exposure to several different air pollutants. In a number of cities in Canada, the AQHI is based on the sum of the mortality risks associated with short-term exposure to different air pollutants. Different combinations of CO, SO2, NO2, O3, PM2.5 and PM10 are included in the AQHI, which is presented on an arbitrary scale of 0–10 [2]. A similar approach is used in Guangzhou in China, based on time-series studies of the associations between exposure to SO2, NO2, O3 and PM2.5 and mortality during the period 2012–2015. The excess mortalities associated with these pollutants are then added in order to construct an AQHI on an arbitrary scale [3]. The AQHI can be applied in different contexts. One application is to present the current air quality situation on a map as an index based on the coexistence of several different pollutants. Another application is for forecasting purposes to provide warnings regarding the air quality situation tomorrow or the next few days considering the risk as a combined effect of several pollutants.
Cairncross et al. [4] defined an Aggregate Risk Index (ARI; similar to AQHI in this paper) to describe the health effects associated with the coexistence of several different air pollutants in the city of Cape Town, South Africa. Daily mortality associated with combined short-term exposure to PM10, PM2.5, SO2, O3, NO2 and CO is calculated, and the total incremental daily mortality risk associated with exposure to these pollutants is scaled to an index value ranging from 0 to 10 [4].
The AQHI concept has also been evaluated by studying the associations between the AQHI value and some health outcomes on a daily basis. In Chen et al. [5], the associations between AQHI and emergency department visits for acute ischemic stroke, based on logistic regression, were investigated during the period 1998–2002 in the city of Edmonton in Alberta, Canada. The AQHI was based on mortality risks associated with a combination of NO2, O3 and PM2.5. A time-stratified case-crossover analysis, adjusted for temperature and humidity, was performed. The AQHI was significantly associated with emergency department visits for acute ischemic stroke during April to September, and with the strongest association in the age group 75 years and above.
In Kousha et al. [6], a time-stratified case-crossover analysis regarding emergency department visits for urticaria was performed during the period 2004–2010 in the city of Windsor in Ontario, Canada. The AQHI was developed based on the mortality risks associated with exposure to NO2, O3 and PM2.5. The correlation between this AQHI and emergency department visits for urticaria was calculated by using logistic regression. Odds ratios for the association between one unit AQHI increase and emergency department visits for urticaria were presented. Significant odds ratios (95% CI) were observed for lags of 2 and 3 days.
Szyszkowicz et al. [7] performed an analysis with the same methodology as the above mentioned studies (including NO2, O3 and PM2.5 in the AQHI), but where the association between one unit increase in AQHI and emergency department visits for asthma was calculated in the city of Windsor in Ontario, Canada during the period 2004–2010. Significant associations between the AQHI levels and emergency department visits for asthma were observed, with the largest value lagged by 9 days in the warm season during April to September.
In To et al. [8], asthma attributed hospitalizations, emergency department visits and outpatient visits, associated with one unit increase in AQHI based on NOx, O3 and PM2.5, were analyzed in Ontario, Canada during the period 2003–2006. Poisson regression was used to estimate the risk ratio for the AQHI, and where adjustments for a number of possible confounders were performed. Statistically significant associations were found for one unit increase in AQHI and asthma-related hospitalizations on the same day (lag 0), and for emergency department visits after two days (lag 2).
In another paper by To et al. [9], the associations between air pollutant concentrations and emergency department visits and hospitalizations due to 11 major chronic diseases were quantified in Ontario, Canada during the period 2003–2010. The AQHI was based on NOx, O3, and PM2.5. The short-term effects associated with an increase in AQHI were estimated with Poisson regression, adjusting for possible confounders. One unit increase in the 10-point AQHI-scale corresponded to increases in the outpatient visits ranging from 1 to 5%.
In summary, an AQHI can be based on the cumulative effects of several different air pollutants and include different types of health outcomes like daily mortality and emergency department visits for acute ischemic stroke, asthma, urticaria or other different chronic diseases. Since the health impact is different depending on pollutant and health outcome, the AQHI will be different.
Normally, synergistic effects of air pollutants are not considered, even though such effects may well exist [10]. In addition, subpopulations may be more sensitive due to other factors such as pollen. Some studies have shown that short-term exposure to high levels of air pollutants may enhance the asthmatic response to allergens [11], but this is seldom described in the epidemiological studies. However, Guilbert et al. [12] studied the Brussels Capital Region over 6 years and reported that the daily number of hospitalizations for asthma increased by 3.2% (95% CI: 1.1–5.3) for an interquartile range increase (40 pollen grains·m–3) in the birch pollen concentration. There was no confounding from pollutants, but an interaction with stronger effects of birch pollen at ozone concentrations above the 85th percentile. Ito et al. [13] found an association between birch pollen concentrations and emergency department visits in New York City. Statistically significant rate ratios for no exposure vs the 98th percentile (680 pollen grains·m–3) were found for lag 1 and lag 2.
Which air pollutants and health outcomes to consider may therefore be critical for the AQHI that is reported to the citizens. The purpose of this study is to construct an AQHI for Stockholm based on the effects on asthma emergency department visits (AEDV) associated with combined exposure to several different air pollutants, and also to evaluate the relative importance of each pollutant during different parts of the year. We also discuss the implication for the use of an AQHI as a tool to inform the general public about the air quality situation in a city.

2. Materials and Methods

The air quality health index in Stockholm has been calculated based on the city-specific beta-coefficients for AEDV associated with short-term exposure to NOx, O3, PM10 and birch pollen according to the following equation:
AQHI = i = 1 p 100   ( e β i X i 1 )
The beta-coefficient (βi) represents the increase in AEDV per 1 µg·m–3 increase of each individual air pollutant (Xi), except for birch pollen, which is the increase of one pollen grain per cubic meter. AQHI is the total percent increase in visits due to pollutants and birch pollen (p). The beta-coefficients are based on all age groups, and the exposure window is lag01, which means the concentration mean of the same day and yesterday. The city-specific beta-coefficients from Stockholm are also compared with the corresponding results from a meta-analysis [14], presented in Table 2. The reason for the comparisons with the meta-coefficients is to see to what extent they differ from our coefficients obtained in Stockholm.
For the assessment, we use concentrations of NOx, O3 and PM10, based on daily mean values from a central site in Stockholm. The monitoring station is a part of the city’s regulatory air pollution control network and equipped with reference instruments according to the EU air quality directive for NOx, O3 and PM10 (Table 1). It is located on a roof 20 m above the street level and represents urban background concentrations. Air quality dispersion modelling has shown that the urban background concentrations recorded are good indicators of the exposure to the population in the Greater Stockholm area, even though the NOx and PM10 exposures are slightly overestimated, while the opposite applies to O3 [15]. Birch pollen were sampled in a Burkard trap located on a roof at Stockholm University, about 4 km north of the city center. Daily average number of pollen grains per cubic meter were recorded manually by the Palynological Laboratory at the Swedish Natural History Museum. For birch pollen, the levels are based on measurements during 2015–2017, and the levels are reported as the daily average number of pollen grains per cubic meter.
Local epidemiological data in terms of daily number of AEDV (ICD10: J45-J46) in Greater Stockholm (population 1.2 million) from January 2001 until December 2005 were collected from the National Patient Register at the Swedish National Board of Health and Welfare. The average daily mean of AEDV was 22.2, with a range of 0–71, and with an interquartile range value of 25. Data on daily pollen counts in Stockholm and daily average urban background (roof level) pollutant levels of NOx and PM10 were obtained from the same stations as described above. For O3, the daily maximum 8-h mean concentration is used. Meteorological observations for Stockholm in terms of temperature and relative humidity were collected from the Swedish Meteorological and Hydrological Institute (SMHI). Data were analyzed by using additive Poisson regression models for over-dispersed counts. The three air pollutants and birch pollen were included as mean of the same day and the day before (lag01), and the meteorological variables as same day (lag 0) and the mean of two previous days (lag12). Time trend and seasonality were adjusted for by using a smooth function (with 3 degrees of freedom per year) using penalized regression splines. Day of week, month and holidays were included with dummy variables, and influenza (according to hospital admissions) was included as a smooth function.
As a validation of the AQHI index, a comparison was made between different regression models determination coefficients (R2adjusted), generated from inclusion of either AQHI or inclusion of the three air pollutants and birch pollen separately. For this purpose, a new set of data was achieved for the same region: daily counts of AEDV for the time period of 1 January 2005 to 30 November 2013 in Stockholm, Sweden. To clarify, in a first step, a base model was generated with daily counts of AEDV and adjusted for calendar variables, time trend, meteorological data and influenza admissions (as described above but without air pollutants and birch pollen). In a second step, we compared the change of the base models R2adjusted values generated when adding the AQHI to the base model, or adding the three air pollutants and birch pollen separately into the base model, for which a higher R2adjusted can be interpreted as a better model.

3. Results

3.1. Validation of the AQHI Index

In our validation analysis of how the AQHI index describes daily observations of AEDV between the years 2005–2013, a higher R2adjusted was obtained when including AQHI compared to the inclusion of NOx, O3, PM10 and birch pollen separately. For all age groups (which the AQHI was derived on), an inclusion of AQHI into the regression model (adjusting for calendar variables, meteorological data and influenza admissions) increased the R2adjusted from 0.743 to 0.765, which is an increase of 13% higher than when the three air pollutants and birch pollen were included separately.

3.2. The Risk-Coefficients that Are Used for the AQHI

In Table 2, the beta-coefficients for AEDV obtained in Stockholm are presented together with the beta-coefficients obtained in a meta-analysis [14]. In the absence of AEDV coefficients for NO2 and O3 in the meta-analysis, the risk increase of hospital admissions for asthma has been used as a substitute, while AEDV has been used for PM10. For Stockholm, the beta-coefficients associated with exposure to PM10 and birch pollen are statistically significant, while the coefficients associated with exposure to NOx and O3 are not. In the meta-analysis, the coefficients for NO2 and PM10 are significant, while the coefficient for O3 is not. However, all the point estimates representing the average values are positive, and these have been used for further calculations. When comparing the point estimates in Stockholm with the point estimates in the meta-analysis, the coefficients for PM10 are reasonable similar, while they are very different for O3. The meta-coefficient for NO2 is also considerably larger than the NOx coefficient in Stockholm. The studies in the meta-analysis are mainly based on studies from the 90’s, and with many studies from limited parts of the world. Possible reasons for the differences between the meta-coefficients and the coefficients obtained in Stockholm are further analyzed in the discussion section.

3.3. The Air Quality Health Index Resulting from the Levels of Air Pollutants and Birch Pollen

The monthly average AQHI during 2015–2017, based on NOx, O3, PM10 and birch pollen, has been calculated based on the observed concentrations and risk coefficients from Stockholm (Figure 1). The increase in AQHI resulting from the levels of birch pollen is most prominent during a short period associated with the flowering that occurs during April and May. But the importance of pollen is also quite different for the three years, as shown in Figure 1. The average total monthly risk increase is 6.3% in 2016, while 5.7% and 5.6% are found for 2015 and 2017.
In Table 3, the average AQHI on a daily basis during the period 2015–2017 is calculated for each pollutant. The calculations based on the coefficients in Stockholm are compared with the coefficients based on the meta-analysis [14]. The average daily AQHI based on the meta-coefficients are, especially for O3, considerably larger in comparison with the coefficients in Stockholm.
In Figure 2, the number of days exceeding the AQHI 90th percentile have been calculated on a monthly basis during the period 2015–2017. The number of exceedances are very different during the year, where both March, April and May stand out as there are relatively large numbers of exceedances associated with both O3 and PM10. April and May stand out with a large number of exceedances for birch pollen, in conjunction with the birch pollen season.

4. Discussion

4.1. Comparing AQHI and AQI

An AQI based on air quality goals or limit values provides information on the air quality situation in relation to limit values and goals, but it may not be representative of the public health risks associated with air pollution exposure, at least not in a direct way. In Stockholm, the AQI reported by the city’s Environment and Health Administration is based on the daily average concentrations of NO2 and PM10 obtained from street-level monitoring stations [17]. For PM10, the AQI is based on the average concentration from four curbside monitoring sites in central Stockholm. Here should be noted that the measurement data for the calculated AQI differ from the calculated AQHI, as for both NOx, O3 and PM10 are based on measurements from an urban background station. The AQI is “high” when the concentrations exceed the daily limit values, “moderately high” above the upper thresholds but below limit values, “moderate” between the lower threshold and upper threshold and “low” when concentrations are below the lower threshold value. Comparisons between AQI and AQHI in Stockholm show that the information to the public may be very different depending on if it is based on limit values or health risk coefficients (Table 4). Even though the mean AQHI increases as AQI increases, there is a large variation of the AQHI. E.g., for AQI = 3, AQHI varies from 1.7 to 9.5, i.e., from quite small increase in AEDV to a large increase in AEDV.
The AQHI based on the sum of NOx, O3, PM10 and birch pollen shows much larger variability than an AQI based on PM10 or NO2, and can become high even when AQI’s indicate low concentrations. This variability is especially prominent when the effect of birch pollen is considered.
This is illustrated in Figure A1, Figure A2 and Figure A3, which show 3-day running averages of AQHI (sum of all pollutants and birch pollen) as compared to the AQI using only PM10 for the years 2015–2017.

4.2. Interaction Effects

An AQHI based on the addition of the individual effect of each pollutant may be wrong if there are significant interactions between the pollutants, or if they are indicators of pollutant mixtures coming from the same source (are correlated), and relative risks used for the AQHI come from single pollutant models. There are human experimental studies indicating both synergistic and antagonistic effects associated with multi-pollutant exposures, but the extrapolation of these findings to environmental exposures should be done with care [10]. Most of these studies have used pollutant concentrations that have been significantly higher compared to the corresponding environmental exposures. The evidence of synergism at normally occurring air pollution exposure levels is limited, and in the absence of evidence of synergism or antagonism, the assumption of additive effects may be considered as the default for risk assessments [10]. In other words, despite possible interaction effects, an AQHI based on the sum of the individual effect of each pollutant in multi-pollutant models is the method that has the greatest scientific support. However, since the AQHI presented here also includes birch pollen, it is also a factor to take into account. Simultaneous air pollution exposure has been shown to be able to increase the allergic response associated with pollen exposure [11,12,18,19]. The potential interaction effects are not possible to quantify in the AQHI; this must be done already in the epidemiological studies, but it is important to be aware of such limitations when it comes to adopting additive effects.

4.3. The Contribution to the AQHI from the Different Pollutants

We have chosen to include pollutants that are indicators of health effects due to emissions from different sources. NOx is an indicator of vehicle exhaust emissions, and in Stockholm, both diesel and gasoline vehicles contribute to the NOx concentrations [20]. The annual variation of NOx is not as large as for PM10 and O3 (Figure 1), but it contributes to an increased AQHI during wintertime with more intensive traffic, more stable atmosphere and lower wind speeds. In contrary, during the summer months, the NOx emissions are lower as the traffic intensity decreases, and the meteorological conditions contribute to greater dilution effects. There is also some anti-correlation between NOx and O3, which is mainly caused by the net destruction of O3 due to the reaction between NO and O3.
For O3, the highest levels and the highest predicted risk increase during the year arise during April and May, which is shown in Figure 1. Tropospheric ozone concentrations in the Northern Europe maximize during spring, and that is caused by both long-distance transport of global tropospheric ozone, and by a large-scale peak in the local photochemical production [21]. In Sweden, the O3 concentrations are lower inside the cities compared to outside due to the titration by NOx [15,22].
The calculated increase in AQHI associated with PM10 is most important during late winter and early spring. This is related to the suspension of road dust particles (>1 µm in diameter) that occurs during late winter and spring in Stockholm, when the road surfaces become drier [23,24]. The road dust is due to the use of studded winter tires [25]. As much as 90% of the mass of PM10 in Stockholm during the spring originates from road abrasion [23]. Exhaust particles represent only a very small proportion of the total mass of PM10 because of their small size. As shown in Figure 1, PM10 contributes a large proportion of the total AQHI. Considering the fact that PM10 to a large extent is locally generated, local action strategies aimed at reducing the PM10 levels may potentially decrease the AQHI.
Birch pollen exhibits a pattern with very high levels during the pollen flowering season that takes place during April to June. As shown in Figure 1, the birch pollen concentrations, and the risk increases associated with those, exhibit a biennial pattern with high concentrations during 2016, but with considerable lower concentrations during 2015 and 2017. This phenomenon is mainly caused by a resource competition between vegetative and reproductive organs of the trees. When the leaves act as a sink of growth substances, this process will also inhibit the development of male inflorescences [26]. An abundant seed production during one year can also reduce the pollen production during the following year [27,28]. However, this biennial pattern is not always distinct, and deviations may sometimes occur [29].

4.4. The Difference between the Local Risk Coefficients and the Meta-Coefficients

The difference between the meta-coefficients and the coefficients in Stockholm (Table 2) indicates that there may be real differences between the locations and periods. The meta-coefficients in Anderson et al. [14], where the associations between AEDV and exposure to NO2, O3 and PM10 are presented, are based on studies performed during the 1990’s, and no studies from the Nordic countries are represented. The meta-coefficient for NO2 is based on five studies, representing Spain, Hong Kong, Italy, Canada and the UK For O3, the meta-coefficient is based on seven studies, representing Spain, Australia, Hong Kong, Italy, the Netherlands, the UK and France. For the meta-coefficient associated with PM10 exposure, two out of the four studies have been performed in the U.S., one study in Canada and one in the UK This causes the prevailing conditions in the represented countries to have relatively large impacts on the meta-coefficients. The coefficients from the meta-analysis have been used in order to see if the coefficients obtained in Stockholm are similar.
Despite that NOx is used in Stockholm, while NO2 is used in the meta-analysis, the mean value of the coefficient in the meta-analysis is almost three times larger compared to the mean value of the coefficient in Stockholm. NOx and NO2 are good indicators of traffic-related exhaust emissions in general, and the use of catalysts and particle filters, which gradually have begun to be used during the 1990’s and until now, may have some significance for this large difference between the coefficients.
For O3, the meta-coefficient is considerably larger than the coefficient found in Stockholm. One possible explanation may be that the O3 levels tend to by high during the birch pollen season. Considering the fact that only one out of the seven studies in the meta-analysis has made adjustments for pollen exposure may explain the relatively large value of the meta-coefficient, where the pollen exposure can contribute to the risk increase.
For PM10, the difference between the coefficients is less compared to the difference for NOx/NO2 and O3. In Stockholm, PM10 consists largely on mechanically generated particle from studded tires, which is not the case for most of the studies included in the meta-coefficient. The relatively small difference between these in both cases significant coefficients is therefore surprising. However, this relatively small difference can be a random occurrence.
For birch pollen, the coefficient found for Brussels [12] is three times larger compared to the coefficient found for Stockholm. There are several factors that can affect this difference, which makes it motivated to conduct further studies related to this phenomenon.
All the coefficients in Stockholm are obtained from a multi-pollutant model by adjusting the effect of each individual pollutant from the effects of the other pollutants, which is not the case for the meta-coefficients, and this can of course also contribute to differences between the coefficients. However, when comparing the coefficients from Stockholm with the meta-coefficients, it turns out that the point estimates of the meta-coefficients are within the confidence intervals of the coefficients in Stockholm.

4.5. Strengths and Limitations of this Study

An advantage of this study is that the coefficients for the different pollutants are adjusted for each other, which means that there is no risk for double calculations when the risk effects of the different pollutants are added. Another advantage is the inclusion of birch pollen, which have major effects during parts of the year, and where there are possible interaction effects with both exhaust emissions and ozone. The validation of the AQHI, presented in the beginning of in Section 3.1, indicate that the index is a better predictor of AEDV related to air pollution exposure in comparison with the inclusion of each pollutant separately, but where the effect of each pollutant is not adjusted for the effects of the others. As far as we can find, none of the cited studies using an AQHI have tried to validate this index with AEDV by comparing the levels of determination with air pollution data. However, a disadvantage of this study is the less precise local risk coefficients for NOx and O3. Follow-up studies, based on coefficients with larger amounts of data, may therefore be beneficial. The air pollution concentrations are based on the measurements from an urban background station, but without considering differences in population exposure, which is also a disadvantage. Another limitation is the use of AEDV, but without considering other more serious outcomes like e.g., mortality.

5. Conclusions

By shifting from a threshold-based index (AQI) to a health-risk based index (AQHI), there are several advantages that can be achieved. The cumulative effects arising from combined exposure to several different air pollutants are illustrated by the use of the AQHI. The relative importance of each pollutant during different parts of the year is also illustrated through this multi-pollutant approach. In Stockholm, the risk increase for AEDV associated with combined exposure varies throughout the year, with the largest risk increase during the period from March to June, in conjunction with relatively high levels of PM10, O3 and birch pollen.
From a policy point of view, there are two particularly important factors that can be achieved by incorporating the AQHI concept in the legislation. The first one is to geographically present the current air quality situation with an index based on the coexistence of several different pollutants. The second one is for forecasting purposes to provide warnings regarding the air quality situation tomorrow or the next few days considering the combined effect of several pollutants. Moreover, public information regarding the health risks associated with different pollutants and pollen can be very useful for sensitive persons when planning their outdoor activities. It may also be important to be able to know if it is pollen or some specific air pollutant that causes the adverse health risks.

Author Contributions

The review of former studies using the AQHI concept has been performed by H.O. The manuscript has mainly been written by H.O. together with C.J. and B.F. The tables and figures have been created by H.O. and C.J.; A.E. has provided the study with birch pollen data, and she has also been a counselor on all pollen-related issues. K.M. has performed the statistical calculations on which the risk coefficients in Stockholm are based. A.T. has performed the statistical analysis that is the basis for the validation of the AQHI. Conceptualization, H.O., C.J. and B.F.; Data curation, C.J., B.F., A.E., K.M. and A.T.; Formal analysis, H.O., C.J., B.F., A.E., K.M. and A.T.; Methodology, C.J., B.F., K.M. and A.T.; Supervision, C.J. and B.F.; Writing—original draft, H.O.; Writing—review & editing, C.J., B.F. and A.T. All authors have contributed to the writing of the paper.

Funding

This project was partly funded by ERA-PLANET (www.era-planet.eu), trans-national project SMURBS (www.smurbs.eu) (Grant Agreement n. 689443) under the EU Horizon 2020 Framework Programme, and partly by the city of Stockholm, as part of the strategy for developing the city’s digitization and information services.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2015. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Figure A1. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2015. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Ijerph 16 00105 g0a1
Figure A2. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2016. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Figure A2. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2016. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Ijerph 16 00105 g0a2
Figure A3. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2017. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Figure A3. Comparison of 3-day running mean values of AQHI total and AQI PM10 in Stockholm during 2017. AQHI is the percentage increase in AEDV associated with combined exposure to NOx, O3, PM10 and birch pollen. AQI is based on a four-point index scale, representing low (PM10 ≤ 25 µg·m–3), moderate (PM10: 25–35 µg·m–3), moderately high (PM10: 35–50 µg·m–3) and high (PM10 ≥ 50 µg·m–3).
Ijerph 16 00105 g0a3

References

  1. Hewings, J. Air Quality Indices: A Review; Pollution Probe: Toronto, ON, Canada, 2001; pp. 13–30. [Google Scholar]
  2. Stieb, D.M.; Burnett, R.T.; Smith-Doiron, M.; Brion, O.; Shin, H.H.; Economou, V. A new multipollutant, no-threshold air quality health index based on short-term associations observed in daily time-series analyses. J. Air Waste Manag. Assoc. 2008, 58, 435–450. [Google Scholar] [CrossRef] [PubMed]
  3. Li, X.; Xiao, J.; Lin, H.; Liu, T.; Qian, Z.; Zeng, W.; Guo, L.; Ma, W. The construction and validity analysis of AQHI based on mortality risk: A case study in Guangzhou, China. Environ. Pollut. 2017, 220A, 487–494. [Google Scholar] [CrossRef] [PubMed]
  4. Cairncross, E.K.; John, J.; Zunckel, M. A novel air pollution index based on the relative risk of daily mortality associated with short-term exposure to common air pollutants. Atmos. Environ. 2007, 41, 8442–8454. [Google Scholar] [CrossRef]
  5. Chen, L.; Villeneuve, P.J.; Rowe, B.H.; Liu, L.; Stieb, D.M. The Air Quality Health Index as a predictor of emergency department visits for ischemic stroke in Edmonton, Canada. J. Expos. Sci. Environ. Epidemiol. 2014, 24, 358–364. [Google Scholar] [CrossRef] [PubMed]
  6. Kousha, T.; Valacchi, G. The air quality health index and emergency department visits for urticaria in Windsor, Canada. J. Toxicol. Environ. Health A 2015, 78, 524–533. [Google Scholar] [CrossRef] [PubMed]
  7. Szyszkowicz, M.; Kousha, T. Emergency department visits for asthma in relation to the Air Quality Health Index: A case-crossover study in Windsor, Canada. Can. J. Public Health 2014, 105, 336–341. [Google Scholar] [CrossRef]
  8. To, T.; Shen, S.; Atenafu, E.G.; Guan, J.; McLimont, S.; Stocks, B.; Licskai, C. The air quality health index and asthma morbidity: A population-based study. Environ. Health Perspect. 2013, 121, 46–52. [Google Scholar] [CrossRef] [PubMed]
  9. To, T.; Feldman, L.; Simatovic, J.; Gershon, A.S.; Dell, S.; Su, J.; Foty, R.; Licskai, C. Health risk of air pollution on people living with major chronic diseases: A Canadian population-based study. BMJ Open 2015, 5, e009075. [Google Scholar] [CrossRef] [PubMed]
  10. Mauderly, J.L.; Samet, J.M. Is There Evidence for Synergy Among Air Pollutants in Causing Health Effects? Environ. Health Perspect. 2009, 117, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Svartengren, M.; Strand, V.; Bylin, G.; Jarup, L.; Pershagen, G. Short-term exposure to air pollution in a road tunnel enhances the asthmatic response to allergen. Eur. Respir. J. 2000, 15, 716–724. [Google Scholar] [CrossRef] [Green Version]
  12. Guilbert, A.; Cox, B.; Bruffaerts, N.; Hoebeke, L.; Packeu, A.; Hendrickx, M.; De Cremer, K.; Bladt, S.; Brasseur, O.; Van Nieuwenhuyse, A. Relationships between aeroallergen levels and hospital admissions for asthma in the Brussels-Capital Region: A daily time series analysis. Environ. Health 2018, 17, 35. [Google Scholar] [CrossRef]
  13. Ito, K.; Weinberger, K.R.; Robinson, G.S.; Sheffield, P.E.; Lall, R.; Mathes, R.; Ross, Z.; Kinney, P.L.; Matte, T.D. The associations between daily spring pollen counts, over-the-counter allergy medication sales, and asthma syndrome emergency department visits in New York City, 2002–2012. Environ. Health 2015, 14, 71. [Google Scholar] [CrossRef] [PubMed]
  14. Anderson, H.R.; Atkinson, R.W.; Bremner, S.A.; Carrinton, 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; Contract Number: 0020017, Report to Department of Health Revised Following First Review 2007; Division of Community Health Sciences St George’s, University of London: London, UK, 2007. [Google Scholar]
  15. Olstrup, H.; Forsberg, B.; Orru, H.; Spanne, M.; Nguyen, H.; Molnár, P.; Johansson, C. Trends in air pollutants and health impacts in three Swedish cities over the past three decades. Atmos. Chem. Phys. 2018. under review. [Google Scholar] [CrossRef]
  16. Soldevilla, C.G.; Gonzales, P.C.; Teno, P.A.; Vilches, D. Spanish Aerology Network (REA): Management and Quality Manual; Servicio De Publicaciones De La Universidad De Córdoba: Córdoba, Spain, 2007; ISBN 978-84-690-6354-5. [Google Scholar]
  17. SLB. Luften i Stockholm. Årsrapport 2011. SLB 1, 2012, SLB Analys. 2011. Available online: http://slb.nu/slb/rapporter/pdf8/slb2012_001.pdf (accessed on 13 August 2018).
  18. Rogerieux, F.; Godfrin, D.; Senechal, H.; Motta, A.C.; Marliere, M.; Peltre, G.; Lacroix, G. Modifications of Phleum pratense grass pollen allergens following artificial exposure to gaseous air pollutants (O3, NO3, SO2). Int. Arch. Allergy Immunol. 2007, 143, 127–134. [Google Scholar] [CrossRef]
  19. Strand, V.; Svartengren, M.; Rak, S.; Barck, C.; Bylin, G. Repeated exposure to an ambient level of NO2 enhances asthmatic response to a non-symptomatic allergen dose. Eur. Respir. J. 1998, 12, 6–12. [Google Scholar] [CrossRef] [PubMed]
  20. Krecl, P.; Johansson, C.; Targino, A.C.; Ström, J.; Burman, L. Trends in black carbon and size-resolved particle number concentrations and vehicle emission factors under real-world conditions. Atmos. Environ. 2017, 165, 155–168. [Google Scholar] [CrossRef]
  21. Monks, P.S. A review of the observations and origins of the spring ozone maximum. Atmos. Environ. 2000, 34, 3545–3561. [Google Scholar] [CrossRef]
  22. Hagenbjörk, A.; Malmqvist, E.; Mattisson Nilsson, K.; Sommar, J.; Modig, L. The spatial variation of O3, NO, NO2 and NOx and the relation between them in two Swedish cities. Environ. Monit. Assess. 2017, 189, 161. [Google Scholar] [CrossRef] [PubMed]
  23. Johansson, C.; Norman, M.; Gidhagen, L. Spatial & temporal variations of PM10 and particle number concentrations in urban air. Environ. Monit. Assess. 2007, 127, 477–487. [Google Scholar] [CrossRef] [PubMed]
  24. Denby, B.R.; Sundvor, I.; Johansson, C.; Pirjola, L.; Ketzel, M.; Norman, M.; Kupiainen, K.; Gustafsson, M.; Blomqvist, G.; Omstedt, G. A coupled road dust surface moisture model to predict non-exhaust road traffic induced particle emissions (NORTRIP). Part 1: Road dust loading and suspension modelling. Atmos. Environ. 2013, 77, 283–300. [Google Scholar] [CrossRef]
  25. Norman, M.; Johansson, C. Studies of some measures to reduce road dust emissions from paved roads in Scandinavia. Atmos. Environ. 2006, 40, 6154–6164. [Google Scholar] [CrossRef]
  26. Dahl, A.; Strandhede, S.O. Predicting the intensity of the birch pollen season. Aerobiologia 1996, 12, 97–106. [Google Scholar] [CrossRef]
  27. Latałowa, M.; Mietus, M.; Uruska, A. Seasonal variations in the atmospheric Betula pollen count in Gdansk (southern Baltic coast) in relation to meteorological parameters. Aerobiologia 2002, 18, 33–43. [Google Scholar] [CrossRef]
  28. Ranta, H.; Oksanen, A.; Hokkanen, T.; Bondestam, K.; Heino, S. Masting by Betula-species, applying the resource budget model to north European data sets. Int. J. Biometeorol. 2005, 49, 146–151. [Google Scholar] [CrossRef]
  29. Grewling, L.; Jackowiak, B.; Nowak, M.; Uruska, A.; Smith, M. Variations and trends of birch pollen seasons during 15 years (1996–2010) in relation to weather conditions in Poznań (western Poland). Grana 2012, 51, 280–292. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The average monthly AQHI resulting from levels of contributing air pollutants in Stockholm during 2015–2017. The risk coefficients in Stockholm have been used in the calculations.
Figure 1. The average monthly AQHI resulting from levels of contributing air pollutants in Stockholm during 2015–2017. The risk coefficients in Stockholm have been used in the calculations.
Ijerph 16 00105 g001
Figure 2. The number of days which exceed the 90th percentile of the risk increase. The calculations are based on the number of monthly exceedances for the whole period from 2015 to 2017.
Figure 2. The number of days which exceed the 90th percentile of the risk increase. The calculations are based on the number of monthly exceedances for the whole period from 2015 to 2017.
Ijerph 16 00105 g002
Table 1. Description of the measurement methods and the instruments that have been used to measure the pollutants included in the AQHI (Air Quality Health Index).
Table 1. Description of the measurement methods and the instruments that have been used to measure the pollutants included in the AQHI (Air Quality Health Index).
PollutantMethodInstrument
NOxChemiluminescenceAC32M, Environnement S.A., Envea, Poissy, France
O3UV absorptionO342M, Environnement S.A. Envea, Poissy, France
PM10GravimetricTEOM 1400A, Thermo Fisher Scientific, Waltham, MA, USA
Birch pollenManual counting in microscope [16]Burkard Seven Day Recording Volumetric Spore Trap (Burkard Manufacturing Co. Ltd., Rickmansworth UK)
Table 2. The percentage increase in asthma emergency department visits (AEDV) associated with 10 µg·m–3 increase of NOx, O3, and PM10 in Stockholm during the period 2001–2005, and from a meta-analysis [14], based on NO2, O3 and PM10. For birch pollen, the percentage increase associated with an increase of ten pollen grains per cubic meter has been used.
Table 2. The percentage increase in asthma emergency department visits (AEDV) associated with 10 µg·m–3 increase of NOx, O3, and PM10 in Stockholm during the period 2001–2005, and from a meta-analysis [14], based on NO2, O3 and PM10. For birch pollen, the percentage increase associated with an increase of ten pollen grains per cubic meter has been used.
PollutantStockholm (Our Analysis)Meta-Analysis [14]
NOx0.5 (95% CI: −1.2–2.2)Not included
NO2Not included1.4 (95% CI: 0.6–2.2)
O30.3 (95% CI: −1.4–2.0)1.6 (95% CI: −0.1–3.3)
PM102.5 (95% CI: 0.3–4.8)3.3 (95% CI: 1.0–5.5)
Birch pollen0.26 (95% CI: 0.18–0.34)Not included
Table 3. The average daily percentage risk increase of each pollutant during the period 2015–2017. The results based on the coefficients in Stockholm are compared with the results based on the meta-coefficients (Table 2). For birch pollen, the risk increase based on the coefficient in Stockholm is compared with the risk increase based on the coefficient obtained in the study from Brussels.
Table 3. The average daily percentage risk increase of each pollutant during the period 2015–2017. The results based on the coefficients in Stockholm are compared with the results based on the meta-coefficients (Table 2). For birch pollen, the risk increase based on the coefficient in Stockholm is compared with the risk increase based on the coefficient obtained in the study from Brussels.
PollutantStockholmMeta-Analysis [14]Brussels [12]
NOx/NO20.6 (NOx)2.0 (NO2)*
O32.010.4*
PM103.14.0*
Birch pollen0.4*1.2
Sum6.116.41.2
* Not included in the analysis.
Table 4. Comparison of AQHI (based on the sum of the risks associated with NOx, O3, PM10 and birch pollen) with AQI for PM10 in Stockholm 2015–2017. There are 1021 days with data for all parameters out of the total number of 1096 days.
Table 4. Comparison of AQHI (based on the sum of the risks associated with NOx, O3, PM10 and birch pollen) with AQI for PM10 in Stockholm 2015–2017. There are 1021 days with data for all parameters out of the total number of 1096 days.
PM10 Index (AQI)PM10 Concentration Intervals (µg·m−3) for AQI 1AQHINumber of Days
MeanMinMax
1 (low)<252.570.529839
2 (moderate)25–354.81.910.195
3 (moderaltely high)35–505.21.79.561
4 (high)≥506.72.711.626
1 Index = 1 corresponds to PM10 concentrations below the lower assessment threshold, index = 2 corresponds to PM10 between the lower assessment threshold and upper assessment threshold, index = 3 corresponds to PM10 above the upper assessment threshold, and index = 4 corresponds to PM10 above or equal to the EU air quality limit value.

Share and Cite

MDPI and ACS Style

Olstrup, H.; Johansson, C.; Forsberg, B.; Tornevi, A.; Ekebom, A.; Meister, K. A Multi-Pollutant Air Quality Health Index (AQHI) Based on Short-Term Respiratory Effects in Stockholm, Sweden. Int. J. Environ. Res. Public Health 2019, 16, 105. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010105

AMA Style

Olstrup H, Johansson C, Forsberg B, Tornevi A, Ekebom A, Meister K. A Multi-Pollutant Air Quality Health Index (AQHI) Based on Short-Term Respiratory Effects in Stockholm, Sweden. International Journal of Environmental Research and Public Health. 2019; 16(1):105. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010105

Chicago/Turabian Style

Olstrup, Henrik, Christer Johansson, Bertil Forsberg, Andreas Tornevi, Agneta Ekebom, and Kadri Meister. 2019. "A Multi-Pollutant Air Quality Health Index (AQHI) Based on Short-Term Respiratory Effects in Stockholm, Sweden" International Journal of Environmental Research and Public Health 16, no. 1: 105. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop