Skip to main content


We're creating a new version of this page. See preview

  • Research
  • Open Access
  • Open Peer Review

Particulate air pollution on cardiovascular mortality in the tropics: impact on the elderly

Contributed equally
Environmental Health201918:34

  • Received: 8 November 2018
  • Accepted: 4 April 2019
  • Published:
Open Peer Review reports



Air pollution has a significant health impact. Most data originate from temperate regions. We aim to study the health impact of air pollution, particularly among the elderly, in a tropical region.


A daily time-series analysis was performed to estimate excess risk (ER) of various air pollutants on daily death counts amongst the general population in Singapore from 2001 to 2013. Air pollutants included particulate matters smaller than 10 μm, and 2.5 μm (PM10, PM2.5), carbon monoxide (CO), nitrogen dioxide (NO2), ozone (O3) and sulphur dioxide (SO2). The studied outcomes were non-accidental and cardiovascular mortality. Single-day lag and distributed lag models were studied and adjusted for confounders.


In single-day lag models, a 10 μg/m3 increase in particulate matter was associated with significant increases in non-accidental (PM10 ER: 0.627%; 95% confidence interval (CI): 0.260–0.995% and PM2.5 ER: 0.660%; 95% CI: 0.204–1.118%) and cardiovascular mortality (PM10 ER: 0.897; 95% CI: 0.283–1.516 and PM2.5 ER: 0.883%; 95% CI: 0.121–1.621%). This was significant in the elderly ≥ 65 years but not in those < 65 years and were seen in the acute phase of lag 0-5 days. Effects by other pollutants were minimal. For cardiovascular mortality, the effects turned protective at a cumulative lag of 30 days in the elderly and could due to “harvesting”.


These first contemporary population-based data from an equatorial country with tropical climate show that exposure to particulate air pollution was significantly associated with non-accidental mortality and cardiovascular mortality, especially in the elderly.


  • Air pollution
  • Cardiovascular mortality
  • Tropics
  • Elderly


Air pollution has a significant health impact on mortality and morbidity worldwide, resulting in an estimated 4.2 million deaths and 103.1 million disability–adjusted life-years lost [1]. Air pollution has wide-ranging systemic effects on the human body, impacting both the respiratory, and cardiovascular systems via multiple mechanisms including oxidative stress, inflammation and endothelial dysfunction [2].

Many prospective cohort and daily time-series studies published globally have consistently demonstrated the negative associations between long and short-term exposure to air pollution and human health [39]. These studies did not only establish associations between ambient particulate matters and respiratory health, but also on cardiovascular health, with the elderly being an especially susceptible group [1012]. However, the vast majority of these studies were conducted in temperate regions [39] rather than in the tropics. As seasonal variations and temperature changes have been shown to impact the relationship of air pollution on health outcomes [1315], we aim to study the impact of air pollution on non-accidental and cardiovascular mortality in the general population, as well as the elderly, in an equatorial country with tropical climate and no seasons.


Data sources

Singapore is an Asian city-state, situated near the equator with a tropical climate, comprising 5.40 million of people. Besides being exposed to daily ambient air pollution generated from domestic sources, Singapore is also exposed almost yearly to haze episodes of about a month long duration, whereby smoke from regional forest fires especially during the dry seasons is blown by winds from neighbouring countries [16].

Air pollutant and meteorological data are comprehensively collected in Singapore. Daily average of 24-h concentrations for particulate matters smaller than 10 μm and 2.5 μm (PM10 and PM2.5), 8-h carbon monoxide (CO), 24-h nitrogen dioxide (NO2), 8-h ozone (O3) and 24-h sulphur dioxide (SO2) from the years of 2001 to 2013 were obtained from the National Environment Agency (NEA) Singapore. During this 13-year study period, air pollution in Singapore was monitored in air monitoring stations located at various sites around Singapore. Four stations located at road-sides were excluded from the study because they did not reflect the daily exposure to air pollutants amongst the general population. Our final analysis included data from 18 studied stations and data completeness for each studied station was assessed by calculating the proportion of days on which data was collected out of the number of days in operation. Aggregated daily air pollutant concentrations were calculated following the Air Pollution on Health: European Approach (APHEA) protocol with added modifications [13]. In summary, imputed annual mean concentrations for each station were subtracted from the station-specific daily concentrations of the same year to generate a set of ‘centred’ values. Using the centred values, the arithmetic mean was calculated across all stations by day and the average of the imputed annual mean concentrations were then added back to derive the daily values used in the analysis. Average meteorological values were derived from daily means of dry bulb temperature and relative humidity from 5 stations which selected based on the shortest proximity to the air monitoring stations provided by NEA’s Meteorological Service division.

Mortality data from the Registry of Births & Deaths were extracted to calculate aggregated daily counts of all non-accidental deaths (International Classification of Diseases (ICD)-9000–799 from the period of 2001–2012 & ICD-10 A00-R99 for the year of 2013) and cardiovascular deaths (ICD-9390–459 & ICD-10 I00-I99) by age groups (all-age, < 65 and ≥ 65 years) over the 13 years study period. A separate analysis was also conducted looking at the subset of subjects ≥ 80 years. Ethics approval was obtained from the SingHealth Centralised Institutional Review Board.

Statistical analysis

A quasi-Poisson generalized additive model was used for the analysis [17]. Single-day lag models and distributed lag models (DLMs) were respectively built to analyse the lagged day effects of the different air pollutants first on all ages combined non-accidental and cardiovascular mortality [1820] respectively. Analyses of PM10 and PM2.5 effects stratified by two age groups (< 65 and ≥ 65 years) were then repeated using the DLMs to further study age-specific effects.

Separate core models were firstly built (for each age group and mortality type) without adding pollutant variable to explain the variations, as much as possible, due to long term trends and other potential time-varying confounding factors:
$$ Log\ E\left[(Y)\right]=\alpha +s(T)+{\beta}_1 DOW+{\beta}_2 SARS+{\beta}_3 FLU+{\beta}_4 PH+{\beta}_5 afterPH+s(DBT)+s(RH) $$

Where Log is the natural logarithm, Y is the daily counts of non-accidental or cardiovascular deaths, DOW refers to the day of week and s(T), s(DBT) and s(RH) refer to the penalised cubic regression smoothers for long term trends, dry bulb temperature (°C) and relative humidity (%) respectively. 1 to 3 days lag were tested for variables DBT and RH. Selection of lag-day variables and smoothing spline parameters were based on models giving the lowest quasi-Akaike’s Information Criterion [21]. Dummy variables of SARS and FLU are used to control for impact due to the periods of 2003 severe acute respiratory syndrome (SARS) and of the 2009 influenza A(H1N1) pandemics respectively while PH and afterPH are terms for public holidays and the day following public holidays.

After the core model was established and chosen for each mortality outcome, variables depicting the various air pollutants were added for the single-day lag models. Two-pollutant models were further constructed for pollutants with significant effects in the single pollutant model. For the DLMs, basis functions were applied with the lag-response relationship defined a priori using a 3rd degree polynomial. Polynomials had been used in previous air pollution studies to analyse the lag-response associations [2224]. A 30-day lag structure, signifying a medium-term period of about a month, was used to observe cumulative effects and mortality displacement if any. As a sensitivity analysis, cumulative lag-response charts generated with DLMs using 2nd and 4th degrees polynomials are shown in the Additional file 1: Figures S1 to S4. Percentage changes in mortality risk, or excess risks (ERs), associated with a 0.1 mg/m3 for CO and 10 μg/m3 for other pollutants were calculated using (RR-1) × 100%, where RR denotes relative risk estimated from the regression coefficient of air pollutant variable of the models. Statistical significance was assessed using 95% confidence intervals (CI). Residual autocorrelation (ACF) and partial autocorrelation (PACF) charts for the core models were performed to assess model fit and are shown in Additional file 1: Figure S5. Data collection process might be interrupted for certain stations over the study period due to site relocation or site closure, and hence annual mean pollutant concentrations used for centring the data were imputed using spatial interpolation with the “fields” package version 8.3–6 in the R statistical software [25]. All statistical analyses, including plotting of lag-response relationships and residual diagnostics to assess model fit, were carried out using the “mgcv” and “dlnm” packages in R version 3.3.1 [17, 19, 26].


Mean daily concentrations of PM10 & PM2.5 were 29.4 μg/m3 and 20.0 μg/m3 respectively, and reached as high as 336.6 μg/m3 and 275.8 μg/m3 during the haze episode in 2013. CO concentration also reached a high of 3.6 mg/m3 during the same episode. Given Singapore’s tropical climate, daily mean of dry bulb temperatures ranged narrowly from 23.4 °C to 30.9 °C while daily mean of relative humidity ranged from 63.6 to 98.1%. Mean daily counts of non-accidental and cardiovascular mortality for all-ages were 43.5 and 15.4 respectively while figures of 30.6 for non-accidental mortality and 11.2 for cardiovascular mortality were recorded in the elderly aged 65 years old and above (Table 1). The correlation between air pollutant and meteorological data is found in the Additional file 1: Table S1. There is a very high linear correlation between levels of PM10 and PM2.5 (r = 0.923).
Table 1

Summary statistics of daily air pollutant concentrations, meteorological data and mortality.2001–2013











Air Pollutants

 PM10 (μg/m3)









 PM2.5 (μg/m3)









 CO (mg/m3)









 NO2 (μg/m3)









 O3 (μg/m3)









 SO2 (μg/m3)









Meteorological variables

 DBT (°C)









 Relative humidity (%)









Daily non-accidental mortality

 All ages









 < 65 years









 > 65 years









Daily cardiovascular mortality

 All ages









 < 65 years (No.)









 > 65 years (No.)









aNo data for ozone on 11/2/2012

Completeness of the air pollutant data ranged from 58.2 to 100.0%. For each pollutant, there were 10 to 17 stations that collected data on more than 80% of the days that they were in operation during the study period. There was only one day in 2012 when data for zone was not available from any of the stations.

Air pollution and mortality

The low residual ACF and PACF values (Additional file 1: Figure S5) showed that the core models were able to explain the temporal trends of the data adequately. Significant effects were observed for PM10 & PM2.5 in both the single-day and distributed lag models. In the single-day lag models, the highest ER associated with a 10 μg/m3 increase occurred in the 3rd day after exposure for both non-accidental (PM10 ER: 0.627, 95% CI: 0.260–0.995% and PM2.5 ER: 0.660%; 95% CI: 0.204–1.118%) and cardiovascular mortality (PM10 ER: 0.897%; 95% CI: 0.283–1.516% and PM2.5 ER: 0.883%; 95% CI: 0.121–1.621%). The estimated ERs generally remained significant for PM10 and PM2.5 even after controlling for additional pollutants (Additional file 1: Table S2a,b).

O3 also exhibited a significant effect at day 1 lag for non-accidental mortality (ER: 0.354%; 95% CI: 0.011–0.698%). However, this effect was insignificant after adjusting for PM10 and PM2.5 (Additional file 1: Table S2c). Effects by other pollutants were observed to be minimal or insignificant (Table 2).
Table 2

Single-day lag models for association between air pollutant levels and mortalitya


Percentage change (95% confidence interval)

Lag 0

Lag 1

Lag 2

Lag 3

Non-accidental mortality


0.485 (0.109, 0.863)

0.487 (0.114, 0.862)

0.481 (0.111, 0.852)

0.627 (0.260, 0.995)


0.593 (0.129, 1.058)

0.518 (0.056, 0.981)

0.504 (0.047, 0.964)

0.660 (0.204, 1.118)


0.064 (−0.237, 0.366)

0.133 (− 0.166, 0.433)

0.240 (− 0.059, 0.539)

0.219 (− 0.077, 0.514)


− 0.124 (− 0.985, 0.744)

− 0.025 (− 0.867, 0.824)

− 0.082 (− 0.947, 0.791)

0.042 (− 0.800, 0.890)


0.015 (− 0.372, 0.403)

0.354 (0.011, 0.698)

0.037 (− 0.356, 0.432)

0.219 (− 0.117, 0.557)


− 0.218 (− 1.044, 0.616)

0.221 (− 0.603, 1.053)

0.679 (− 0.147, 1.512)

−0.390 (− 1.208, 0.435)

Cardiovascular mortality


0.553 (− 0.074, 1.185)

0.633 (0.007, 1.262)

0.610 (− 0.010, 1.234)

0.897 (0.283, 1.516)


0.632 (−0.142, 1.412)

0.603 (−0.172, 1.385)

0.541 (− 0.229, 1.318)

0.883 (0.121, 1.651)


0.167 (−0.328, 0.665)

0.071 (−0.418, 0.561)

0.135 (− 0.356, 0.629)

0.362 (− 0.132, 0.858)


− 0.344 (− 1.707, 1.038)

−0.151 (− 1.486, 1.203)

−0.897 (− 2.217, 0.441)

−0.352 (− 1.726, 1.041)


0.500 (− 0.134, 1.137)

0.292 (− 0.266, 0.854)

0.105 (− 0.454, 0.667)

0.489 (− 0.154, 1.137)


0.223 (− 1.114, 1.577)

−0.325 (− 1.647, 1.014)

0.857 (− 0.479, 2.211)

−0.197 (− 1.519, 1.143)

aPer 0.1 mg/m3 increase in CO and per 10 μg/m3 increase for other pollutants

Results highlighted in bold are statistically significant

Results from the distributed lag models are presented in Table 3 and Fig. 1. Significant cumulative effects were observed over days 0–5 but not over subsequent longer periods. Cumulative ERs over days 0–5 for PM10 & PM2.5 were estimated to be 0.700% (95% CI: 0.276–1.126%) & 0.829% (0.276–1.386%) respectively for non-accidental [30] mortality, and 0.921% (0.218–1.629%) & 1.073% (0.157–1.998%) respectively for cardiovascular mortality. Sensitivity analysis using 2nd or 4th degree polynomials showed similar results (Additional file 1: Figures S1-S4).
Table 3

Distributed lag models for association between air pollutant levels and mortalitya


Percentage change (95% confidence interval)

Lags 0–5

Lags 0–15

Lags 0–30

Non-accidental mortality


0.700 (0.276, 1.126)

0.316 (−0.277, 0.912)

0.220 (−0.700, 1.148)


0.829 (0.276, 1.386)

0.362 (−0.419, 1.150)

0.297 (−0.904, 1.512)


0.260 (−0.051, 0.572)

0.074 (−0.368, 0.518)

0.660 (− 0.015, 1.339)


−0.622 (−1.713, 0.480)

−1.560 (−3.119, 0.024)

−0.106 (−2.239, 2.073)


0.123 (−0.373, 0.622)

− 0.628 (− 1.305, 0.054)

− 1.184 (− 2.200, − 0.159)


−0.125 (− 1.380, 1.146)

0.120 (− 1.715, 1.989)

1.048 (− 1.514, 3.677)

Cardiovascular mortality


0.921 (0.218, 1.629)

0.011 (−0.948, 0.989)

− 1.078 (− 2.489, 0.354)


1.073 (0.157, 1.998)

−0.196 (− 1.451, 1.074)

− 1.450 (− 3.263, 0.396)


0.278 (− 0.232, 0.791)

− 0.151 (− 0.854, 0.557)

0.110 (− 0.922, 1.152)


−1.345 (− 3.034, 0.375)

−3.112 (−5.365, − 0.805)

−1.685 (−4.498, 1.211)


0.435 (− 0.350, 1.227)

−0.785 (− 1.809, 0.249)

−2.165 (− 3.586, − 0.722)


−0.095 (− 2.098, 1.948)

−0.963 (− 3.747, 1.901)

−1.032 (− 4.653, 2.727)

aPer 0.1 mg/m3 increase in CO and per 10 μg/m3 increase for other pollutants

Results highlighted in bold are statistically significant

Fig. 1
Fig. 1

a Cumulative percent change (%) in non-accidental mortality for 10 μg/m3 pollutant concentration increase in PM10 (top row) and PM2.5 (bottom row). b Cumulative percent change (%) in cardiovascular mortality for 10 μg/m3 pollutant concentration increase in PM10 (top row) and PM2.5 (bottom row)

Impact by age

Cumulative effects were significant in the elderly aged 65 years old and above, with every 10 μg/m3 increase in pollutant concentration resulting in a 0.771% (95% CI: 0.265–1.279%) and a 0.955% (0.297–1.618%) increase over days 0–5 in non-accidental mortality risk for PM10 & PM2.5 respectively. Cardiovascular mortality in turn showed a 1.236% (0.436–2.042%) and a 1.478% (0.437–2.530%) increase in risk of PM10 & PM2.5 over days 0–5 (Table 4).
Table 4

Distributed lag models for association between air pollutant levels and mortality by age groupa


Percentage change (95% confidence interval)

Lags 0–5

Lags 0–15

Lags 0–30

Non Accidental Mortality

 < 65 years


0.105 (− 0.641, 0.856)

− 0.726 (− 1.659, 0.216)

− 1.214 (− 2.442, 0.029)


−0.073 (− 1.039, 0.903)

− 1.045 (− 2.246, 0.171)

−1.694 (− 3.220, − 0.143)

 > 65 years


0.771 (0.265, 1.279)

0.376 (− 0.345, 1.102)

0.146 (− 0.980, 1.284)


0.955 (0.297, 1.618)

0.447 (− 0.501, 1.405)

0.241 (− 1.229, 1.733)

Cardiovascular Mortality

 < 65 years


− 0.706 (− 2.038, 0.645)

− 0.885 (− 2.466, 0.721)

−1.456 (− 3.411, 0.539)


−1.145 (−2.872, 0.613)

− 1.302 (− 3.332, 0.771)

−2.015 (− 4.417, 0.447)

 > 65 years


1.236 (0.436, 2.042)

− 0.288 (− 1.366, 0.801)

−2.085 (− 3.614, − 0.532)


1.478 (0.437, 2.530)

−0.655 (− 2.060, 0.771)

−2.721 (− 4.667, − 0.736)

aPer 10 μg/m3 increase in PM10 and PM2.5

Results highlighted in bold are statistically significant

No significant effects were seen over days 0–15, but interestingly significant negative cumulative effects were seen for cardiovascular mortality over days 0–30 in the elderly group (PM10 ER: -2.085%; 95% CI: -3.614% – -0.532% and PM2.5 ER: -2.721%; 95% CI: -4.667% – -0.736%). A similar negative effect was also seen for the younger age group for non-accidental mortality over days 0–30 for PM2.5 (ER:-1.694%; 95% CI: -3.220% – -0.143%). Sensitivity analysis using up to 2nd or 4th degree of polynomials showed similar results (Additional file 1: Figures S1-S4).

In a separate analysis looking at the very elderly aged 80 years and above, cumulative effects were significant, with every 10 μg/m3 increase in pollutant concentration resulting in a 0.749% (95% CI: 0.095–1.408%) and a 0.955% (0.105–1.812%) increase over days 0–5 in non-accidental mortality risk for PM10 & PM2.5 respectively (data not shown). There was a trend towards increased cardiovascular mortality at days 0–5, although this association was not significant. Significant negative cumulative effects were seen for cardiovascular mortality over days 0–30 in this very elderly group for PM10 (ER: -2.060%; 95% CI: -3.987% – -0.094%) and a similar non-signficant trend for PM2.5.

Lag-response association

Graphs showing non-cumulative percentage change in mortality over days 0–30 are shown in the Additional file 1: Figure S6, S7. The graphs for the all-age group and the ≥ 65 years age group showed that impact on mortality was mostly observed in the immediate days (days 0–5) after exposure. After this initial period, significant protective effects were actually observed in both non-accidental & cardiovascular mortality. However, while non-accidental mortality showed effects approaching the null around days 15–20, cardiovascular mortality only showed a rebound around days 20–25, thus indicating a slower approach.


In this daily time series analysis performed in an equatorial Asian city-state with tropical climate, we found significant associations between particulate air pollutants (PM2.5, PM10) and non-accidental mortality, as well as cardiovascular mortality, which corroborated with other previous studies conducted in temperate countries, but our estimated effect sizes are closer to the lower end of their estimates. Greater effects were especially found in the elderly. The impacts of other pollutants were generally insignificant however.

The strongest evidence for the impact of air pollutants on mortality lie in the particulate air pollutants. In a recent meta-analysis, a 10 μg/m3 increment in PM2.5 was associated with a 1.04% increase in mortality with substantial regional variations (ranging from 0.25 to 2.08%) [27]. In a study from 20 cities in the United States, a 10 μg/m3 increased in PM10 was associated with significant increase of 0.51% in overall mortality, with no impact from the other air pollutants [7]. Another study from 38 cities in China similarly showed a 10 μg/m3 change in PM10 concentrations was associated with a 0.44% increase in daily number of deaths [3]. In Europe, PM2.5 was significantly associated with overall mortality while NO2 was not [5]. The ESCALA study (Estudio de Salud y Contaminación del Aire en Latinoamérica) which included 9 Latin American cities in the tropics (albeit non-equatorial) found a similar significant association between PM10 and mortality [28]. In our study, we found an increase of 0.63 and 0.66% in overall mortality with every 10 μg/m3 increase in PM10 and PM2.5 respectively from the single-day lag model. No significant impact on overall mortality was found from the other air pollutants. The lack of impact from other air pollutants despite some findings from other studies may in part be explained by the lower concentrations of these air pollutants in our study. For example, the NO2 concentrations in our study were about half that of several published studies [22, 29] which found significant associations between NO2 and mortality. Similarly, significant lower concentrations of O3 was noted in our study compared to others [28].

Since the hypothesis was first postulated by Seaton et al. (1995), that ultra-fine particulates, less than 0.3 μm in aerodynamic diameter, may provoke more alveolar inflammation causing exacerbation of existing lung disease and increased blood coagulability than larger respired particles, the impact of air pollution on cardiovascular health has received great attention. In addition, air pollutants result in increase in risk of cardiovascular events via various pathways including oxidative stress, systemic inflammation, endothelial dysfunction, atherothrombosis, and arrhythmogenesis [2, 29, 30]. Increased levels of air pollutants have been shown to be associated with progression in coronary calcification on computed tomography scans [8]. The same meta-analysis mentioned in the previous paragraph found a 0.84% increase in cardiovascular mortality with a 10 μg/m3 in PM2.5 [27]. The ESCALA study from tropical Latin American also found a significant association between PM10 and cardiovascular events [28]. In study of 8 cities in China, a 10 μg/m3 in PM10 was associated with 0.36% increase in coronary heart disease mortality, with significant associations also found for NO2 and SO2 [31]. A study from Korea found that exposure to PM2.5, PM10, CO, SO2, and NO2, but not O3, increased the risks of cardiovascular events and mortality [6].

In our analysis of the single-day lag model, we found an increase of 0.90 and 0.88% in cardiovascular mortality with every 10 μg/m3 increase in PM10 and PM2.5 respectively. Interestingly, in the single-day lag model, the earliest significant effects on overall mortality were seen on the same day as exposure while the effects on cardiovascular mortality were primarily seen on lag day 3. Although both lags are fairly acute, this may indicate a relatively more delayed effect of air pollution on the cardiovascular system compared to the respiratory system (the respiratory system being in direct contact with the inhaled pollutants). Furthermore, although the pro-ischemic and pro-thrombotic effects of pollution may occur within hours of exposure [32], it may take longer (eg. days) for cardiovascular mortality to occur. A similar finding was demonstrated by Costa et al., whereby the impact of air pollution on cardiovascular mortality was evident only on lag day 3 while the impact on respiratory mortality was more immediate [22]. We did not find an impact of other pollutants on cardiovascular mortality.

The elderly is an increasingly vulnerable group and sensitive to the effects of particulate matter. The elderly tends to be frailer, have more co-morbidities and less physiological reserves. In the Chinese study, the effect of PM10 on overall mortality in those less than 60 years old was insignificant but became significant in those aged ≥ 60 years old (0.57% per 10μg/m3 increase) [3]. In the Netherlands, PM10 caused increased overall mortality in all age groups but only caused significant circulatory disease mortality in those aged > 65 years old [4]. Multiple other studies have also demonstrated the ill effects from air pollution in the elderly age group. In our study, we found similar significant associations with overall and cardiovascular mortality with PM10 and PM2.5 in the elderly aged ≥65 years but not in those < 65 years. The significant effects on mortality were seen primarily in the acute phase of lag 0-5 days. Interestingly, for cardiovascular mortality, the effects turned protective at a cumulative lag of 30 days in the elderly. This suggests that the number of cardiovascular deaths was lower as compared to another 30-day window period where no exposure to particulate matter had taken place. This could be attributed to the “mortality displacement” or “harvesting” effect [23, 33]. In the elderly population, high level of short term exposure to air pollution may cause greater mortality, advancing deaths and depleting the ‘at-risk’ group early on, resulting in a follow-on period with a mortality rate that is lower than expected [23, 33]. In a study on air pollution and myocardial infarction, higher levels of PM10 were associated with short-term increased risk of MI, but later reductions in risk suggest that air pollution may be associated with bringing events forward in time (“short-term displacement”) [34]. Another reason for the apparent lower number of deaths later on might be the subsequent protective measures (e.g. wearing masks, staying indoors) that alerted individuals or organizations put in place during periods of high air pollution levels. This may in part explain the lower non-accidental mortality in the younger age group at cumulative lag of 0–30 days.

The vast majority of studies on the impact of air pollution were performed in temperate regions. The meta-analysis used data mainly from European, American and Western Pacific regions, with only one study on hospital admissions that included data from south-east Asia [27]. The Global Burden of Diseases Study included data from tropical regions but air pollution burden was largely estimated from satellite data and also extrapolated and applied risk estimates of air pollution from other studies for places where such estimates were not available [1]. Moreover, studies have shown that seasonal changes, temperature and humidity variation may impact the relationships between air pollution effect and health outcomes [1315], emphasizing the importance of validating these findings in the tropics, where the climate is starkly different. In the APHEA study, there was higher effects of PM10 on daily mortality in cities with a higher temperatures [13] and similarly in Korea, hotter temperatures accentuated the effect of PM10 on mortality [14]. In our study, we found similar significant findings on the negative impact of particulate air pollution on overall and cardiovascular mortality. The estimates from our study, although not directly comparable to other studies, appear to be on the lower end of the spectrum.


There are several limitations. Firstly, due to the time-series analysis nature of the study, unmeasured clinical variables (eg. comorbidities), which could affect the actual estimates of air pollution on mortality, could not be fully accounted for despite our best efforts in controlling for potential confounding factors in the models. Our models also did not adjust for the transboundary haze episodes as we wanted to analyse the general exposure of the public to air pollutants, and not to distinguish between domestic and transboundary air pollution. Secondly, the impact of indoor air pollution was not assessed as this data was not readily available, but we believe such effect would be minimal as residents in Singapore do not burn coals for cooking in the houses. Thirdly, a 3rd degree of polynomials was used to define the lag-response relationship over a period of 30 days. A different parametric shape or a different degree parameter would result in different ER estimates. However, a sensitivity analysis using up to 2nd or 4th degree of polynomials showed that these differences were minute and did not affect our final interpretation. A fourth limitation is that we did not explore the inclusion of different concentration-response functions to describe the exposure-response relationship [35]. In our analysis, we assumed a log-linear relationship via a quasi-Poisson model. Lastly, correlation between the paired pollutants in the two-pollutant models potentially resulted in multicollinearity that produced some unstable estimates.


These first contemporary population-based data from an equatorial country with tropical climate and no seasons show that exposure to particulate air pollution (PM2.5, PM10) was significantly associated with non-accidental mortality and cardiovascular mortality, especially in the elderly.




Air Pollution on Health: European Approach


Confidence interval


Carbon monoxide


Distributed lag model


Excess risk


International Classification of Diseases


National Environment Agency


Nitrogen dioxide




Particulate matters < 10 μm


Particulate matters < 2.5 μm


Relative risk


Sulphur dioxide



Not applicable


No sources of funding required.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available as these are obtained with special permission from the relevant national agencies in Singapore but are available from the corresponding author on reasonable request.

Authors’ contributions

JY, YN, SM designed the study, conducted the literature review, performed the analysis and wrote the mansucript. YKK, VL designed the study, performed the analysis and made critical revisions to the manuscript. AS and CSPL designed the study and made critical revisions to the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Ethics approval was obtained from the Singhealth Centralised Institutional Review Board.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Department of Cardiology, National Heart Centre Singapore, 5 Hospital Drive, Singapore, 169609, Singapore
Public Health Group, Ministry of Health, Singapore, Singapore
Duke-NUS Graduate Medical School, Singapore, Singapore
Karolinska Institutet, Stockholm, Sweden
Saw Swee Hock School of Public Health, National University of Singapore, Singapore, Singapore


  1. Cohen AJ, Brauer M, Burnett R, et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015. Lancet. 2017;389:1907–18.View ArticleGoogle Scholar
  2. Mills NL, Donaldson K, Hadoke PW, et al. Adverse cardiovascular effects of air pollution. Nat Clin Pract Cardiovasc Med. 2009;6:36–44.View ArticleGoogle Scholar
  3. Yin P, He G, Fan M, Chiu KY, et al. Particulate air pollution and mortality in 38 of China's largest cities: time series analysis. BMJ. 2017;356:j667.View ArticleGoogle Scholar
  4. Fischer PH, Marra M, Ameling CB, et al. Air pollution and mortality in seven million adults: the Dutch environmental longitudinal study (DUELS). Environ Health Perspect. 2015;123:697–704.View ArticleGoogle Scholar
  5. Beelen R, Raaschou-Nielsen O, Stafoggia M, et al. 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–95.View ArticleGoogle Scholar
  6. Kim H, Kim J, Kim S, et al. Cardiovascular effects of long-term exposure to air pollution: a population-based study with 900 845 person-years of follow-up. J Am Heart Assoc. 2017;6:e007170.Google Scholar
  7. Samet JM, Dominici F, Curriero FC, et al. Fine particulate air pollution and mortality in 20 U.S. cities, 1987-1994. N Engl J Med. 2000;343:1742–9.View ArticleGoogle Scholar
  8. Kaufman JD, Adar SD, Barr RG, et al. Association between air pollution and coronary artery calcification within six metropolitan areas in the USA (the multi-ethnic study of atherosclerosis and air pollution): a longitudinal cohort study. Lancet. 2016;388:696–704.View ArticleGoogle Scholar
  9. Di Q, Wang Y, Zanobetti A, et al. Air pollution and mortality in the Medicare population. N Engl J Med. 2017;376:2513–22.View ArticleGoogle Scholar
  10. Zeng Y, Gu D, Purser J, et al. Associations of environmental factors with elderly health and mortality in China. Am J Public Health. 2010;100:298–305.View ArticleGoogle Scholar
  11. Cakmak S, Dales RE, Vidal CB. Air pollution and mortality in Chile: Susceptibility among the elderly. Environ Health Perspect. 2007;115:524–7.View ArticleGoogle Scholar
  12. Fischer P, Hoek G, Brunekreef B, et al. Air pollution and mortality in the Netherlands: are the elderly more at risk? Eur Respir J Suppl. 2003;40:34s–8s.View ArticleGoogle Scholar
  13. Katsouyanni K, Touloumi G, Spix C, et al. Short term effects of ambient Sulphur dioxide and particulate matter on mortality in 12 European cities: results from time series data from the APHEA project. BMJ. 1997;314:1658–63.View ArticleGoogle Scholar
  14. Kim SE, Lim YH, Kim H. Temperature modifies the association between particulate air pollution and mortality: a multi-city study in South Korea. Sci Total Environ. 2015;524-525:376–83.View ArticleGoogle Scholar
  15. Roberts S. Interactions between particulate air pollution and temperature in air pollution mortality time series studies. Environ Res. 2004;96:328–37.View ArticleGoogle Scholar
  16. National Environmental Agency, Singapore. Haze. Available at Accessed 1 Mar 2019.
  17. Wood SN. Generalized additive models: an introduction with R. Chapman and Hall/CRC; 2006.View ArticleGoogle Scholar
  18. Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology. 2000;11:320–6.View ArticleGoogle Scholar
  19. Gasparrini A. Distributed lag linear and non-linear models in R: the package dlnm. J Stat Softw. 2011;43:1–20.View ArticleGoogle Scholar
  20. Roberts S. An investigation of distributed lag models in the context of air pollution and mortality time series analysis. J Air Waste Manage Assoc. 2005;55:273–82.View ArticleGoogle Scholar
  21. Anderson DR, Burnhum KP, White GC. AIC model selection in overdispersed capture-recapture data. Ecology. 1994;75:1780–93.View ArticleGoogle Scholar
  22. Costa AF, Hoek G, Brunekreef B, et al. Air pollution and deaths among elderly residents of São Paulo, Brazil: an analysis of mortality displacement. Environ Health Perspect. 2017;125:349–54.View ArticleGoogle Scholar
  23. Zanobetti A, Wand MP, Schwartz J, Ryan LM. Generalized additive distributed lag models: quantifying mortality displacement. Biostatistics. 2000;1:279–92.View ArticleGoogle Scholar
  24. Schwartz J. Harvesting and long term exposure effects in the relation between air pollution and mortality. Am J Epidemiol. 2000;151:440–8.View ArticleGoogle Scholar
  25. Nychka D, Furrer R, Paige J, et al. Fields: tools for spatial data. R package version 8.3-6. 2016. Available at Scholar
  26. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. Available at Google Scholar
  27. Atkinson RW, Kang S, Anderson HR, et al. Epidemiological time series studies of PM2.5 and daily mortality and hospital admissions: a systematic review and meta-analysis. Thorax. 2014;69:660–5.View ArticleGoogle Scholar
  28. A. Romieu I, Gouveia N, Cifuentes LA, de Leon AP, Junger W, Vera J, Strappa V, Hurtado-Díaz M, Miranda-Soberanis V, Rojas-Bracho L, Carbajal-Arroyo L, Tzintzun-Cervantes G, HEI Health Review Committee. Multicity study of air pollution and mortality in Latin America (the ESCALA study). Res Rep Health Eff Inst. 2012;171:5–86.Google Scholar
  29. Goggins WB, Chan EY, Yang C, Chong M. Associations between mortality and meteorological and pollutant variables during the cool season in two Asian cities with sub-tropical climates: Hong Kong and Taipei. Environ Health. 2013 Jul;12:59.View ArticleGoogle Scholar
  30. Newby DE, Mannucci PM, Tell GS, et al. Expert position paper on air pollution and cardiovascular disease. Eur Heart J. 2015;36:83–93b.View ArticleGoogle Scholar
  31. Li H, Chen R, Meng X, et al. Short-term exposure to ambient air pollution and coronary heart disease mortality in 8 Chinese cities. Int J Cardiol. 2015;197:265–70.View ArticleGoogle Scholar
  32. Mills NL, Törnqvist H, Gonzalez MC, et al. Ischemic and thrombotic effects of dilute diesel-exhaust inhalation in men with coronary heart disease. N Engl J Med. 2007;357:1075–82.View ArticleGoogle Scholar
  33. Zeger SL, Dominici F, Samet J. Harvesting- resistant estimates of air pollution effects on mortality. Epidemiology. 1999;10:171–5.View ArticleGoogle Scholar
  34. Bhaskaran K, Hajat S, Armstrong B, et al. The effects of hourly differences in air pollution on the risk of myocardial infarction: case crossover analysis of the MINAP database. BMJ. 2011;343:d5531.View ArticleGoogle Scholar
  35. Apte JS, Marshall JD, Cohen AJ, et al. Addressing global mortality from ambient PM2.5. Environ Sci Technol. 2015;49:8057–66.View ArticleGoogle Scholar


© The Author(s). 2019


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.