The association between long-term exposure to low-level PM2.5 and mortality in the state of Queensland, Australia: A modelling study with the difference-in-differences approach
Authors:
Wenhua Yu aff001; Yuming Guo aff001; Liuhua Shi aff003; Shanshan Li aff002
Authors place of work:
Department of Epidemiology, School of Public Health and Management, Binzhou Medical University, Yantai, Shandong, China
aff001; Department of Epidemiology and Preventive Medicine, School of Public Health and Preventive Medicine, Monash University, Melbourne, Australia
aff002; Gangarosa Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, Georgia, United States of America
aff003
Published in the journal:
The association between long-term exposure to low-level PM2.5 and mortality in the state of Queensland, Australia: A modelling study with the difference-in-differences approach. PLoS Med 17(6): e32767. doi:10.1371/journal.pmed.1003141
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pmed.1003141
Summary
Background
To date, few studies have investigated the causal relationship between mortality and long-term exposure to a low level of fine particulate matter (PM2.5) concentrations.
Methods and findings
We studied 242,320 registered deaths in Queensland between January 1, 1998, and December 31, 2013, with satellite-retrieved annual average PM2.5 concentrations to each postcode. A variant of difference-in-differences (DID) approach was used to investigate the association of long-term PM2.5 exposure with total mortality and cause-specific (cardiovascular, respiratory, and non-accidental) mortality. We observed 217,510 non-accidental deaths, 133,661 cardiovascular deaths, and 30,748 respiratory deaths in Queensland during the study period. The annual average PM2.5 concentrations ranged from 1.6 to 9.0 μg/m3, which were well below the current World Health Organization (WHO) annual standard (10 μg/m3). Long-term exposure to PM2.5 was associated with increased total mortality and cause-specific mortality. For each 1 μg/m3 increase in annual PM2.5, we found a 2.02% (95% CI 1.41%–2.63%; p < 0.01) increase in total mortality. Higher effect estimates were observed in Brisbane than those in Queensland for all types of mortality. A major limitation of our study is that the DID design is under the assumption that no predictors other than seasonal temperature exhibit different spatial-temporal variations in relation to PM2.5 exposure. However, if this assumption is violated (e.g., socioeconomic status [SES] and outdoor physical activities), the DID design is still subject to confounding.
Conclusions
Long-term exposure to PM2.5 was associated with total, non-accidental, cardiovascular, and respiratory mortality in Queensland, Australia, where PM2.5 levels were measured well below the WHO air quality standard.
Keywords:
Death rates – Physical activity – Socioeconomic aspects of health – Age groups – Particulates – Census – Observational studies – Air quality
Introduction
Ambient particulate matter with diameters less than 2.5 micrometers (fine particulate matter; PM2.5) has been identified as the fifth leading mortality risk factor in 2015, contributing to 4.2 million deaths [1]. While a growing body of literature [2–6] has examined the causal effect of PM2.5 on mortality, scientific evidence remains weak. Causal model approaches seek to mimic randomized controlled trials, whereby exposure is measured independent of the other predictors of the health outcome. This effectively eliminates the possibility of confounding. Propensity score matching and inverse probability weighting are the most common approaches to reweight the study populations to ensure exposure is independent of all measured confounders [7]. Wu and colleagues applied a propensity score approach to estimate the causal effect of long-term PM2.5 exposures on mortality in New England [8]. Wang and colleagues employed a doubly robust causal modelling approach with inverse probability weights to estimate the hazards of long-term exposure to PM2.5 on survival in the southeast United States [9]. However, this method only accounts for measured confounders and unmeasured biases that are highly correlated with the measured confounders [10]. Another approach to estimating the causal effect is the regression discontinuity design. This approach compares observations lying closely on either side of a threshold to estimate the average treatment effect [11]. In one such study, Avraham Ebenstein and colleagues employed the regression discontinuity design to investigate the causal effect of particulate matter with a diameter of less than 10 micrometers (PM10) on Chinese life expectancy by assigning different coal subsidy policies for indoor heating on both sides of Huai River [12]. However, the estimated effects could still be misled by other potential confounders that occur at the same threshold [13].
To address this limitation, a difference-in-differences (DID) approach has been proposed [14]. DID estimates the effect of an exposure or treatment on an outcome by comparing the average change over time in the outcome variable for a treatment group, compared to the average change over time for a control group. A DID method mimics an experimental research design using observational study data to provide a causal estimate by adjusting for unmeasured confounders. It assumes that differences between outcomes that change over time are caused by the differences between the observed and counterfactual exposures, rather than by other factors such as socioeconomic status (SES), population, smoking, and obesity, because such factors among other unmeasured confounders are similar between locations across time. Several studies have utilized the DID method to explore the association between PM and mortality [2,15–18]. For example, Corrigan and colleagues used a DID approach to examine the association between changes in PM2.5 and changes in cardiovascular mortality rates before and after the implementation of a new PM2.5 National Ambient Air Quality Standards (NAAQS) in the US [17]. Wang and colleagues estimated the causal effect of long-term exposure to PM2.5 on mortality in New Jersey by developing a variant of the DID approach [2]. Similar models were used to assess the long-term PM2.5-mortality association in 207 cities across the US [18].
Most of the DID studies estimated the effect of long-term exposure to PM and mortality at relatively high concentrations [15,16,19]. For example, Matteo Renzi and colleagues used a DID approach to estimate the effect of annual PM10 exposure with a range of 21.9 ±4.9 μg/m3 in the Latium region, Italy [15]. Another study’s estimated mortality changes were specifically attributable to high exposure to PM2.5 in China between 2000 and 2010 [19]. However, growing studies have shown that there is no obvious PM2.5 threshold for the PM2.5-mortality association even at PM2.5 levels under the World Health Organization (WHO) air quality guideline (10 μg/m3 of annual average PM2.5). A meta-analysis including 14 studies conducted on participants with average exposure to PM2.5 below 10 μg/m3 supported the nonlinear PM2.5-mortality exposure-response association, where the effect increased rapidly at lower concentrations [20]. However, few studies [3,4] have explored the association of PM2.5 with mortality below a concentration of 10 μg/m3. Therefore, exploring the relationship between long-term low levels of PM2.5 and cause-specific mortality is warranted, particularly in areas with consistently low PM2.5 concentrations in Australia.
In addition, significant differences exist between rural and urban areas in terms of population characteristics, air pollution concentration, and chemical and physical composition. These may contribute to different health effects on cause-specific mortality, especially in Queensland, where approximately half of the population lives in the Brisbane metropolitan area. To address this gap, we assessed the association of long-term exposure to low-level PM2.5 with total mortality and cause-specific mortality from cardiovascular, respiratory, and non-accidental causes in Queensland and Brisbane during 1998–2013 using a DID approach [2].
Methods
This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 STROBE Checklist). We did not include a formal prospective analysis plan, but the study methodology and analysis were planned before conducting the DID approach, with the exception of the analysis of relative change rate and the effect modification by stratifying the age groups. This study was approved by the Monash University Human Research Ethics Committee.
Study area
Queensland is the second largest state in Australia, with an area of 1,852,642 square kilometers, and located in the northeast of the country [21] (Fig 1). Typically, Queensland experiences two weather seasons: a winter with mild temperatures and minimal rainfall, and a humid summer with both high temperatures and high levels of rainfall. As of 2016, it was composed of 4,689,134 inhabitants (449 postcode zones), mainly concentrated along the coast and in the state’s South East. Brisbane was the largest city in the state, with 2,109,466 residents in 2016. In this study, Brisbane was divided into a total of 119 postcode zones.
Data resources
Individual-level data, such as sex, age, and total and cause-specific mortality from 1 January, 1998, to 31 December 2013, were obtained from Queensland Health (https://www.health.qld.gov.au/public-health). Daily data on temperature were obtained from the Scientific Information for Land Owners (SILO) data set (https://www.longpaddock.qld.gov.au/silo/). Annual population and SES census data (including weekly income per person, Socio-Economic Indexes for Areas [SEIFA]) were provided by the Australian Bureau of Statistics (https://www.abs.gov.au/).
All mortality data were coded with the ICD-9 code before 1999 or ICD-10 code after that time. Specifically, the diseases of interest included the following: mental and behavioral disorders (F00–F99), diseases of the nervous system (G00–G99), diseases of the circulatory system (I00–I99), diseases of the respiratory system (J00–J99), diseases of the digestive system (K00–K93), diseases of the genitourinary system (N00–N99), and external causes of morbidity and mortality (V01–Y98). In this study, we focus on the annual counts of death from non-accidental causes (all above diseases except for V01–Y98), cardiovascular (ICD-9: 390–459, ICD-10: I00–I99), and respiratory causes (ICD-9: 460–519, ICD-10: J00–J99) during the study period (1998–2013).
Exposure data
The annual mean PM2.5 data were derived from the Atmospheric Composition Analysis Group, which estimated annual PM2.5 concentrations at 0.01° × 0.01° (approximately 1 km × 1 km) spatial resolution globally using a Geographically Weight Regression with an out-of-sample cross-validated R2 of 0.81. The regression incorporated satellite data, simulated aerosol composition, and land use information [22]. The population-weighted annual mean PM2.5 concentrations were assigned to each postcode area.
Statistical analyses
DID method
This study was based on a variant of DID design [2] to investigate the relationship of PM2.5 with total mortality and cause-specific mortality (non-accidental, cardiovascular, and respiratory causes). The substance of the DID design is that through comparing the same population to itself at different time points, some unmeasured individual and behavioral factors that remain constant over time have been controlled. In this study, the DID design was used to estimate the effect of PM2.5 on mortality by comparing the concordance between differences in counts of cause-specific deaths and differences in PM2.5 over time (from 1998 to 2013) in a given location (449 spatial units of postcode areas in this study). Specifically, a typical DID model is as follows:
In our study, we fit a model as follows:
Our model is based on the following assumptions: (1) The DID relies on a parallel trend assumption that in the absence of intervention (e.g., in the absence of the impact of PM2.5 or the concentration of PM2.5 remains constant in this study), the unobserved differences among locations are constant over time. Although there is no statistical test for this assumption, annual trends were evaluated visually for relative changes in PM2.5 and mortality rate over 16 years. We applied a relative change rate (RCc,t) to calculate the percentage change for each area unit in each year with the following equation: RCc,t = (Rc,t−Ec)/Ec×100; Ec=∑t=19982013Rc,t16, where RCc,t denotes the annual percent changes of PM2.5 or mortality rate in area c, year t; Rc,t is the annual concentration of PM2.5 or the cause-specific mortality rate in the same stratum; and Ec denotes the average of Rc,t from 1998 to 2013 in each area unit. (2) We assumed that, apart from seasonal temperature, no predictors exhibit differential spatial-temporal variations in relation to the exposure [2,15]. Under that assumption, overall spatial and temporal confounding is removed from the DID design. However, if this assumption is violated, the DID design is still subject to confounding.
Conditional Poisson regression
We applied a conditional Poisson regression model [24] to perform the DID design, to estimate the association of long-term exposure to PM2.5 with mortality using the “gnm” package in R software (version 3.2.5). We adjusted for summer and winter temperatures and SEIFA effect, conditional on strata of spatial units. Additionally, we also estimated whether the effect was modified by different proportions of age and sex based on Census 2016 data using subgroup analyses.
Sensitivity analysis
We performed sensitivity analyses to test the robustness of the results using a random-effects meta-analysis, estimating whether the effect estimate in a specific population was disparate from the pooled effect estimate in Queensland. The potential nonlinearity of the association between PM2.5 and total and cause-specific mortality was examined using cubic splines with 3 degrees of freedom. We also modelled the summer and winter temperatures using natural splines with 3 and 4 degrees of freedom, respectively, to assess for modification of health effects of PM2.5 by season. In order to control the impact of the population age structure, we tested the effect modification by stratifying the population into two age groups: <65 years and ≥65 years, using the offset term of age-specific person-years. All analyses were conducted using R software (version 3.2.5).
Results
We studied 242,320 deaths in 7 categories of diseases (ICD10: F00–F99, G00–G99, I00–I99, J00–J99, K00–K93, N00–N99, V01–Y98) from 1998 to 2013, which accounted for 60.5% of registered deaths during the study period. Specifically, 217,510 non-accidental deaths, 133,661 deaths from cardiovascular diseases, and 30,748 deaths from respiratory diseases in the Queensland region were investigated. Table 1 displays the distribution of deaths, PM2.5, and temperature in Brisbane and Queensland over the study period. In short, 81.9% of deaths were over 65 years old and 50.8% were male. The average PM2.5 concentrations were 6.0 μg/m3 (range: 2.13–8.00 μg/m3) in Brisbane (S1 Fig) and 3.63 μg/m3 (range: 1.63–9.00 μg/m3) in Queensland (Fig 1), with interquartile ranges (IQRs) of 2.00 μg/m3 and 2.23 μg/m3, respectively. There was a decrease of 9.33% in the standardized death rate in Queensland over the 16-year study period (S1 Table). There was no significant change in PM2.5 or temperature during the study period. The annual average PM2.5 ranged from 2.01 to 5.28 μg/m3, along with a downward trend from 1998 to 2008 and with slight increase between 2008 and 2013. The range of seasonal mean temperature varied from 24.5°C to 26.5°C in summer and from 15.0°C to 16.4°C in winter (S1 Table). In terms of the relative changes of PM2.5 and mortality rate, we calculated the percentage changes (RCc,t) of both the annual concentration of PM2.5 and the total mortality rates for all area units in the study period, as presented in Fig 2. As depicted, changes in mortality rates loosely follow changes in the PM2.5 concentrations, which can support our parallel trend assumption, although a reverse trend was observed in 2010.
In this study, we found a significant association between long-term exposure to PM2.5 and total mortality, with 2.02% (95% CI 1.41%–2.63%; p < 0.01) and 5.65% (95% CI 4.08–7.25%; p < 0.01) increases in total mortality per 1 μg/m3 increase in annual PM2.5 in Queensland and Brisbane, respectively. We also observed increases in cause-specific mortality associated with elevated PM2.5 levels (Table 2). Higher effect estimates were observed in Brisbane than those in Queensland for all types of mortality.
We additionally investigated the association between long-term PM2.5 and cause-specific mortality grouped by age and sex (Table 3). In Queensland, individuals under the age of 65 had a high risk of total, non-accidental, and cardiovascular death, with 5.76% (95% CI 4.29%–7.25%; p < 0.01), 6.18% (95% CI 4.29%–8.11%; p < 0.01), and 6.50% (95% CI 4.08%–9.20%; p < 0.01), respectively. In contrast, residents in Brisbane over the age of 65 were more likely to have elevated risk with 5.02% (95% CI 3.46%–6.50%; p < 0.01) and 4.08% (95% CI 2.22%–6.08%; p < 0.01) in non-accidental and cardiovascular mortality, respectively, for every 1 μg/m3 increase in annual average PM2.5. We also found larger effect estimates for the total mortality in females versus males in Queensland (2.84%, 95% CI 1.92%–3.67%; p < 0.01 versus 1.31%, 95% CI 0.50%–2.22%; p < 0.01).
Sensitivity analyses indicate our main findings are robust. With a random effect meta-analysis, we found that there were statistically significant pooled effects with a 2.76% (95% CI 0.61%–3.97%; p < 0.01) increase in total mortality, a 2.28% (95% CI 0.61%–3.97%; p < 0.01) increase in cardiovascular mortality, a 5.55% (95% CI 3.75%–7.39%; p < 0.01) increase in respiratory mortality, and a 2.88% (95% CI 1.37%–4.41%; p < 0.01) increase in non-accidental mortality (Fig 3 and S2 Table). The associations between PM2.5 and total/cause-specific mortality tend to be nonlinear, with a threshold around 4.5 μg/m3 for PM2.5 exposure (Fig 4). In addition, after modification by age group, we estimated a 5.21% (95% CI 3.20%–7.25%; p < 0.01) increase in total mortality among the <65 age group, whereas we estimated a 0.95% (95% CI −1.20% to 3.14%; p = 0.39) increase among subjects ≥65 years old, which were consistent with our model results (S3 Table). Moreover, the results were similar after adjusting for both the annual average and SD of temperature using natural splines with 3 and 4 degrees of freedom, respectively (S4 Table).
Discussion
This study examined the association of low levels of PM2.5 (<9.0 μg/m3) with cause-specific mortality using a DID approach, which controls for potential unmeasured and omitted confounders. We found that long-term exposure to PM2.5 was associated with increased risks of total mortality and cause-specific mortality, despite low-level concentrations falling consistently below the current WHO annual standard (10 μg/m3). Furthermore, the effect estimates were higher in the urban area (Brisbane) when compared to statewide estimates for the study period.
In the presence of unmeasured confounders, it is difficult to determine causality in observational studies, especially for relatively weak health risk factors such as particulate matter. The PM2.5-mortality association may be affected by factors that are not perfectly measured or routinely collected, such as SES and outdoor physical activities. Therefore, if these unmeasured factors are omitted or unavailable, it might result in biased effect estimates [25,26]. It should be noted that many factors affect PM2.5-attributable mortality rates [27]. In addition to the variations witnessed in nonlinear exposure–response relationship across different causes of deaths [20,27–29], regional variations in ambient PM2.5, population characteristics, and baseline disease incidence rates all contribute to variations of PM2.5-attributable mortality [27]. Even though controlled exposure studies [30] and randomized participants [31] could reduce certain sources of bias and control for the unmeasured confounders, ethical considerations could be given on account of the potential for toxic components in the chemical composition of PM2.5.
However, a well-designed observational study such as a DID method may overcome the limitations of a nonrandomized observational study [18,32]. A case in a DID study is compared to itself at different time points so that certain confounders (including the unmeasured ones) such as population structure and lifestyle factors that remain stable or rarely varied during the study period are cancelled out because the comparisons occur among the population in the same places [15,19,33,34]. In our case, we controlled for certain slow-changing spatial-temporal variations by design, such as population and SES, which play a role in confounding the PM2.5-mortality association [18]. We also included seasonal temperature, which varies in different periods across the study areas and correlated to PM2.5 exposure and mortality as well as other potential factors. Controlling for temperature may indirectly remove the impacts of other factors like influenza epidemics. Therefore, we assumed that most unmeasured potential and omitted confounders had been controlled.
Several studies have used a similar DID design to investigate the potential association between long-term PM exposure and mortality [2,16–19]. For comparison purposes, we have converted the percent change in our study to a 10 μg/m3 increase and found a 22.14% (95% CI 15.02%–29.69%) increase in total mortality for the entire population of Queensland. This finding was slightly higher than that of Wang and colleagues [2], who developed a variant of DID approach to study 1938 census tracts in New Jersey from 2004 to 2009. They reported that each 10 μg/m3 increase in annual PM2.5 was associated with a 15.5% (95% CI 0.8%–32.3%) increase in the natural-cause mortality. Different outcome metrics might explain differences in results, e.g., we used 7 causes of death in Queensland instead of the natural-cause mortality. Another study estimated the city-specific health effects of PM2.5 on mortality in 207 US cities between 2000 and 2010. They observed a link between long-term PM2.5 and mortality with a hazard ratio (HR) of 1.2 (95% CI 1.1–1.3) for each 10 μg/m3 increase in annual PM2.5. However, direct comparisons between these results and the current study should be cautious because they combined the city-specific exposure, which tends to have different mixture composition of ambient particles [35].
The yields in this study were higher compared to many previous cohort studies [36–40]. For example, a recent representative cohort of American adults [40] estimated the HR of 1.12 (95% CI 1.08–1.15) for all-cause mortality and 1.23 (95% CI 1.17–1.29) for cardiopulmonary mortality per 10 μg/m3 long-term exposure to PM2.5. Most of the cohort studies are based on the exposure levels at a fixed time or area and consider some potential confounders that have been measured or observed. Our study considers the variability in exposure and potential influences in time and space, and the unmeasured confounders were captured by design. While cohort studies primarily rely on survey design and methodologies to obtain individual-level data, our study focuses on the spatiotemporal effects of PM2.5 on mortality at a macro scale by comparing differences in mortality and exposure of PM2.5 by postal code. In a systematic review, Vodonos and colleagues observed higher estimates in the hybrid space-time model than other methods [20]. In addition, another possible explanation for the relative high results in our study might be that the PM2.5 concentrations were much lower than in other studies, with a mean exposure of 3.63 μg/m3 (<9.00 μg/m3) in Queensland. Previous studies [20,41] have provided evidence for a nonlinear PM2.5-mortality association, where mortality increased sharply with low exposure levels and leveled off at higher exposure. A recent systematic review [20] estimated a 1.29% (95% CI 1.09–1.50) increase in all-cause mortality per 1 μg/m3 increase in PM2.5 at a mean exposure of 10 μg/m3, which decreased to 1.03% (95% CI 0.97–1.11) at 15.7 μg/m3 and to 0.82% (95% CI 0.52–1.12) at 30 μg/m3. They also restricted concentrations of PM2.5 to <10 μg/m3 and found a 2.4% (95% CI 0.80–4.00) increase per 1 μg/m3, which is consistent with our results.
Our results show a PM2.5-mortality association at levels below the current WHO air quality standard. Throughout the study period, the annual PM2.5 concentrations across Queensland areas were well below the current the WHO standard (10 μg/m3 of annual average PM2.5) and the US Environmental Protection Agency (EPA) standard (12 μg/m3 of annual average PM2.5); however, the PM2.5-mortality association remained present in spite of this. A growing body of research has suggested that low-level PM2.5 exposure may increase mortality. Shi and colleagues estimated a 9.28% increase in mortality for every 10-μg/m3 increase in annual PM2.5 for populations ≥65 years in the US. [42]. Aligning with our study, Markar and colleagues [3] and Schwartz and colleagues [4] used different causal reference methods to estimate the effects on mortality at concentrations of PM2.5 below the standards.
Mounting toxicological literature has provided evidence for the causal effect of PM2.5 on mortality [43–46]. For example, one study on animals [45] reported that mice with 6 months exposure to a low concentration of PM2.5 compared with animals exposed to filtered air demonstrated marked increases in atherosclerotic plaque, macrophage infiltration, and vasoconstrictor responses in the aortic arch. Another study [43] found more severe lung dysfunction in mice exposed for 8 months to PM2.5 of 16.8 μg/m3 when compared to animals exposed to PM2.5 of 2.9 μg/m3.
Several underlying biological mechanisms have been investigated for the damaging effects of PM2.5 on organ systems, especially on the respiratory and cardiovascular system [47–50]. The first is the direct pathway. Ultrafine particles directly translocate into the bloodstream and into specific organs, which aggravate the local oxidative stress and inflammation, causing the atherosclerotic plaque instability, and ultimately induce the cardiotoxicity effects and increase the risk of congestive heart failure, arrhythmias, and cardiovascular mortality [47–49]. Another pathway may increase oxidative stress and activate inflammation. The free radicals and organic components of PM2.5 can generate a rich milieu of inflammatory mediators and induce free radical production to oxidize lung cells, which directly causes cell injury in the lungs [51] and may indirectly release into the blood or systemic circulation, leading to cardiovascular and pulmonary disease and even death [50,52].
We also found that PM2.5 had larger effects on mortality in the urban area (Brisbane) when compared to statewide. In our study, the elderly inhabitants in Brisbane displayed a higher increase in death relative to statewide rates (4.92%, 95% CI 3.46%–6.61%; p < 0.01, versus 1.21%, 95% CI 0.06%–1.92%; p < 0.01). One possible explanation is the consistently higher concentrations of PM2.5 found in Brisbane (with an annual average of 6.0 μg/m3) compared to Queensland (with an annual average of 3.63 μg/m3). Furthermore, chemical and physical particle compositional differences in rural and urban areas also show different health effects on mortality [35,53]. Increasing evidence suggests that observed rural-urban disparities in population density could result in differences in PM2.5-associated mortality rates [54]. In our study, compared with the rest of Queensland, approximately half of the overall population (44.9% based on Census 2016 data) resides in Brisbane, resulting in more widespread exposure to PM2.5 than in rural areas (S2 Fig and S3 Fig).
Growing evidence has shown that the elderly are at a higher risk of total mortality attributable to particulate matter [42,55]. Contrary to these findings, however, our study found that young people (below the age of 65) in Queensland experienced higher risks comparatively (5.76%, 95% CI 4.29%–7.25%; p < 0.01, versus 1.21%, 95% CI 0.06%–1.92%; p < 0.01). In Brisbane, those over 65 had a higher risk of mortality associated with PM2.5 in comparison with people below 65. This may be explained by differential exposure and physical activity between urban and rural residents, which may vary by age [56]. Compared with older adults, young people in rural areas may have higher rates of outdoor physical activity, which may lead to greater air pollution exposure [57]. With regards to sex sub-analysis, in our study, females in Brisbane experienced a greater risk than females in Queensland. We found that females in our study had higher PM2.5-attributable mortality than the male population in both Brisbane and statewide. Even though an inverse relationship between all-cause mortality and PM2.5 exposure among female farmers was found in an American study [58], it is noted that there is a significant difference in individual-level behaviors between urban and rural residents, such as physical activity and smoking [59].
Our study has several strengths. To the best of our knowledge, no other studies have comprehensively explored the link between long-term exposure to PM2.5 and cause-specific mortality in Queensland when controlling for unmeasured confounding by design. Our study is based on 7 categories of disease-specific mortality, with a relevantly long study period (16 years). Additionally, we employed a variant of DID approach to explore the association between long-term exposure and mortality. Our results have provided compelling evidence for an association between long-term exposure to PM2.5 and mortality at levels below the current WHO air quality standard (10 μg/m3 of annual average PM2.5). Finally, we investigated the PM2.5-attributable mortality in the highly populated metropolis and found more severe risks on cause-specific mortality in Brisbane.
There are also some limitations to this study. We assumed that no predictors other than seasonal temperature exhibit different spatial-temporal variations in relation to PM2.5 exposure. However, there are other potential spatial and temporal confounders that vary by periods across the study areas and correlate with the PM2.5 exposure, such as employment rate [2,60] and influenza epidemic [15]. Furthermore, even though the DID study design may eliminate most unmeasured confounders, these factors held as ideal assume that confounders such as population, SEIFA, and some other unmeasurable factors like behavior habits in one area maintained unchanged over the study period. Additionally, our study was unable to measure individual-level exposure and potentially introduced ecological bias because the study was conducted on a population-level scale. Because the assessment of the environmental exposures was based on individual zip code, errors in geocoding and invariable yearly exposure in every postal code for all residents may lead to potential exposure misclassification. Moreover, we cannot estimate the association between PM2.5 and natural-cause mortality because of limited data availability. Furthermore, the basic population size and social-economic data in our study were based on Census 2016, which is likely to change over time.
This study provides evidence that long-term exposure to PM2.5, even at low levels well below the current WHO air quality standard, is associated with non-accidental, cardiovascular, and respiratory mortality in Queensland and Brisbane, by using a variant of DID approach to control the unmeasured confounding. Even though the explanation of the DID approach relies on several assumptions that are theoretical and, to some extent, unprovable, the findings are important for scientific understanding of the health effects of air pollution and to inform policy makers.
Supporting information
S1 STROBE Checklist [doc]
STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.
S1 Table [sd]
Environmental and mortality data from 1998 to 2013.
S2 Table [docx]
The pooled effects of PM-mortality associations for different death types by using a random effect meta-analysis.
S3 Table [docx]
The effect modification of PM on mortality by age group in Queensland and Brisbane.
S4 Table [docx]
Assessing the robustness of natural splines with different degrees for the summer and winter temperature in Queensland.
S1 Fig [asgs]
The annual average PM concentrations in Brisbane during 1998–2013.
S2 Fig [asgs]
The population density in Queensland based on 2016 census.
S3 Fig [asgs]
The population density in Brisbane based on 2016 census.
Zdroje
1. Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, 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(10082): 1907–18. doi: 10.1016/S0140-6736(17)30505-6 28408086
2. Wang Y, Kloog I, Coull BA, Kosheleva A, Zanobetti A, Schwartz JD. Estimating Causal Effects of Long-Term PM2.5 Exposure on Mortality in New Jersey. Environ Health Perspect. 2016;124(8): 1182–8. doi: 10.1289/ehp.1409671 27082965
3. Makar M, Antonelli J, Di Q, Cutler D, Schwartz J, Dominici F. Estimating the Causal Effect of Low Levels of Fine Particulate Matter on Hospitalization. Epidemiology. 2017;28(5): 627–34. doi: 10.1097/EDE.0000000000000690 28768298
4. Schwartz J, Bind M-A, Koutrakis P. Estimating causal effects of local air pollution on daily deaths: effect of low levels. Environ Health Perspect. 2016;125(1): 23–9. doi: 10.1289/EHP232 27203595
5. HEI Health Review Committee. Causal inference methods for estimating long-term health effects of air quality regulations: Health Effects Institute; 2016 [cited 2019 Apr 17]. Available from: https://www.healtheffects.org/publication/causal-inference-methods-estimating-long-term-health-effects-air-quality-regulations
6. Schwartz JD, Wang Y, Kloog I, Yitshak-Sade Ma, Dominici F, Zanobetti A. Estimating the effects of PM 2.5 on life expectancy using causal modeling methods. Environ Health Perspect. 2018;126(12): 127002. doi: 10.1289/EHP3130 30675798
7. Rubin DB. The design versus the analysis of observational studies for causal effects: parallels with the design of randomized trials. Stat Med. 2007;26(1): 20–36. doi: 10.1002/sim.2739 17072897
8. Wu X, Mealli F, Kioumourtzoglou M-A, Dominici F, Braun D. Matching on Generalized Propensity Scores with Continuous Exposures. arXiv:1812.06575 [Preprint]. 2018 [cited 2019 Jul 15]. Available from: https://arxiv.org/abs/1812.06575
9. Wang Y, Lee M, Liu P, Shi L, Yu Z, Awad YA, et al. Doubly robust additive hazards models to estimate effects of a continuous exposure on survival. Epidemiology. 2017;28(6): 771. doi: 10.1097/EDE.0000000000000742 28832358
10. Garrido MM, Kelley AS, Paris J, Roza K, Meier DE, Morrison RS, et al. Methods for constructing and assessing propensity scores. Health Serv Res. 2014;49(5): 1701–20. doi: 10.1111/1475-6773.12182 24779867
11. Imbens GW, Lemieux T. Regression discontinuity designs: A guide to practice. J Econom. 2008;142(2): 615–35. doi: 10.1016/j.jeconom.2007.05.001
12. Ebenstein A, Fan M, Greenstone M, He G, Zhou M. New evidence on the impact of sustained exposure to air pollution on life expectancy from China’s Huai River Policy. Proc Natl Acad Sci U S A. 2017;114(39): 10384–9. doi: 10.1073/pnas.1616784114 28893980
13. Lee DS. Randomized experiments from non-random selection in US House elections. J Econom. 2008;142(2): 675–97. doi: 10.1016/j.jeconom.2007.05.004
14. Card D, Krueger AB. Minimum wages and employment: A case study of the fast food industry in New Jersey and Pennsylvania. National Bureau of Economic Research. 1993. doi: 10.3386/w4509
15. Renzi M, Forastiere F, Schwartz J, Davoli M, Michelozzi P, Stafoggia M. Long-Term PM10 Exposure and Cause-Specific Mortality in the Latium Region (Italy): A Difference-in-Differences Approach. Environ Health Perspect. 2019;127(6): 67004. doi: 10.1289/EHP3759 31166133
16. Leogrande S, Alessandrini ER, Stafoggia M, Morabito A, Nocioni A, Ancona C, et al. Industrial air pollution and mortality in the Taranto area, Southern Italy: A difference-in-differences approach. Environ Int. 2019;132: 105030. doi: 10.1016/j.envint.2019.105030 31398654
17. Corrigan AE, Becker MM, Neas LM, Cascio WE, Rappold AG. Fine particulate matters: The impact of air quality standards on cardiovascular mortality. Environ Res. 2018;161: 364–9. doi: 10.1016/j.envres.2017.11.025 29195185
18. Kioumourtzoglou M-A, Schwartz J, James P, Dominici F, Zanobetti A. PM2. 5 and mortality in 207 US cities: modification by temperature and city characteristics. Epidemiology. 2016;27(2): 221. doi: 10.1097/EDE.0000000000000422 26600257
19. Xue T, Zhu T, Zheng Y, Liu J, Li X, Zhang Q. Change in the number of PM2. 5-attributed deaths in China from 2000 to 2010: Comparison between estimations from census-based epidemiology and pre-established exposure-response functions. Environ Int. 2019;129: 430–7. doi: 10.1016/j.envint.2019.05.067 31154145
20. Vodonos A, Awad YA, Schwartz J. The concentration-response between long-term PM2. 5 exposure and mortality; A meta-regression approach. Environ Res. 2018;166: 677–89. doi: 10.1016/j.envres.2018.06.021 30077140
21. Xu Z, FitzGerald G, Guo Y, Jalaludin B, Tong S. Assessing heatwave impacts on cause-specific emergency department visits in urban and rural communities of Queensland, Australia. Environ Res. 2019;168: 414–9. doi: 10.1016/j.envres.2018.10.013 30388498
22. Van Donkelaar A, Martin RV, Brauer M, Hsu NC, Kahn RA, Levy RC, et al. Global estimates of fine particulate matter using a combined geophysical-statistical method with information from satellites, models, and monitors. Environ Sci Technol. 2016;50(7): 3762–72. doi: 10.1021/acs.est.5b05833 26953851
23. Shi L, Kloog I, Zanobetti A, Liu P, Schwartz JD. Impacts of temperature and its variability on mortality in New England. Nat Clim Chang 2015;5(11): 988–91. doi: 10.1038/nclimate2704 26640524
24. Armstrong BG, Gasparrini A, Tobias A. Conditional Poisson models: a flexible alternative to conditional logistic case cross-over analysis. BMC Med Res Methodol. 2014;14(1): 122. doi: 10.1186/1471-2288-14-122 25417555
25. Barrowman MA, Peek N, Lambie M, Martin GP, Sperrin M. How unmeasured confounding in a competing risks setting can affect treatment effect estimates in observational studies. BMC Med Res Methodol. 2019;19(1): 166. doi: 10.1186/s12874-019-0808-7 31366331
26. Fewell Z, Davey Smith G, Sterne JA. The impact of residual and unmeasured confounding in epidemiologic studies: a simulation study. Am J Epidemiol. 2007;166(6): 646–55. doi: 10.1093/aje/kwm165 17615092
27. Apte JS, Marshall JD, Cohen AJ, Brauer M. Addressing global mortality from ambient PM2. 5. Environ Sci Technol. 2015;49(13): 8057–66. doi: 10.1021/acs.est.5b01236 26077815
28. Burnett RT, Pope CA III, Ezzati M, Olives C, Lim SS, Mehta S, et al. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Perspect. 2014;122(4): 397–403. doi: 10.1289/ehp.1307049 24518036
29. Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. Ambient particulate air pollution and daily mortality in 652 cities. N Engl J Med. 2019;381(8): 705–15. doi: 10.1056/NEJMoa1817364 31433918
30. Mills NL, Törnqvist H, Gonzalez MC, Vink E, Robinson SD, Söderberg S, et al. Ischemic and thrombotic effects of dilute diesel-exhaust inhalation in men with coronary heart disease. N Engl J Med. 2007;357(11): 1075–82. doi: 10.1056/NEJMoa066314 17855668
31. Chuang H-C, Ho K-F, Lin L-Y, Chang T-Y, Hong G-B, Ma C-M, et al. Long-term indoor air conditioner filtration and cardiovascular health: A randomized crossover intervention study. Environ Int. 2017;106: 91–6. doi: 10.1016/j.envint.2017.06.008 28624750
32. Streeter AJ, Lin NX, Crathorne L, Haasova M, Hyde C, Melzer D, et al. Adjusting for unmeasured confounding in nonrandomized longitudinal studies: a methodological review. J Clin Epidemiol. 2017;87: 23–34. doi: 10.1016/j.jclinepi.2017.04.022 28460857
33. Callaway B, Sant'Anna PH. Difference-in-differences with multiple time periods. arXiv: 1803.09015[Preprint]. 2018 [cited 2019 Jul 19]. Available from: https://arxiv.org/abs/1803.09015.
34. French B, Heagerty PJ. Analysis of longitudinal data to evaluate a policy change. Stat Med. 2008;27(24): 5005–25. doi: 10.1002/sim.3340 18618416
35. Valavanidis A, Fiotakis K, Vlachogianni T. Airborne particulate matter and human health: toxicological assessment and importance of size and composition of particles for oxidative damage and carcinogenic mechanisms. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev 2008;26(4): 339–62. doi: 10.1080/10590500802494538 19034792
36. Lepeule J, Laden F, Dockery D, Schwartz J. Chronic exposure to fine particles and mortality: an extended follow-up of the Harvard Six Cities study from 1974 to 2009. Environ Health Perspect. 2012;120(7): 965–70. doi: 10.1289/ehp.1104660 22456598
37. Beelen R, Raaschou-Nielsen O, Stafoggia M, Andersen ZJ, Weinmayr G, Hoffmann B, 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(9919): 785–95. doi: 10.1016/S0140-6736(13)62158-3 24332274
38. Dehbi H-M, Blangiardo M, Gulliver J, Fecht D, de Hoogh K, Al-Kanaani Z, et al. Air pollution and cardiovascular mortality with over 25 years follow-up: A combined analysis of two British cohorts. Environ Int. 2017;99: 275–81. doi: 10.1016/j.envint.2016.12.004 27939045
39. Wang Y, Shi L, Lee M, Liu P, Di Q, Zanobetti A, et al. Long-term exposure to PM2. 5 and mortality among older adults in the southeastern US. Epidemiology. 2017;28(2): 207. doi: 10.1097/EDE.0000000000000614 28005571
40. Pope CA III, Lefler JS, Ezzati M, Higbee JD, Marshall JD, Kim S-Y, et al. Mortality risk and fine particulate air pollution in a large, representative cohort of US adults. Environ Health Perspect. 2019;127(7): 077007. doi: 10.1289/EHP4438 31339350
41. Pope CA III, Burnett RT, Turner MC, Cohen A, Krewski D, Jerrett M, et al. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: shape of the exposure-response relationships. Environ Health Perspect. 2011;119(11): 1616–21. doi: 10.1289/ehp.1103639 21768054
42. Shi L, Zanobetti A, Kloog I, Coull BA, Koutrakis P, Melly SJ, et al. Low-concentration PM2.5 and mortality: estimating acute and chronic effects in a population-based study. Environ Health Perspect. 2015;124(1): 46–52. doi: 10.1289/ehp.1409111 26038801
43. Mauad T, Rivero DHRF, de Oliveira RC, de Faria Coimbra Lichtenfels AJ, Guimarães ET, de Andre PA, et al. Chronic exposure to ambient levels of urban particles affects mouse lung development. Am J Respir Crit Care Med. 2008;178(7): 721–8. doi: 10.1164/rccm.200803-436OC 18596224
44. Soares SRC, Carvalho-Oliveira R, Ramos-Sanchez E, Catanozi S, da Silva LFF, Mauad T, et al. Air pollution and antibodies against modified lipoproteins are associated with atherosclerosis and vascular remodeling in hyperlipemic mice. Atherosclerosis. 2009;207(2): 368–73. doi: 10.1016/j.atherosclerosis.2009.04.041 19486979
45. Sun Q, Wang A, Jin X, Natanzon A, Duquaine D, Brook RD, et al. Long-term air pollution exposure and acceleration of atherosclerosis and vascular inflammation in an animal model. JAMA. 2005;294(23): 3003–10. doi: 10.1001/jama.294.23.3003 16414948
46. Sun Q, Yue P, Kirk RI, Wang A, Moatti D, Jin X, et al. Ambient air particulate matter exposure and tissue factor expression in atherosclerosis. Inhal Toxicol. 2008;20(2): 127–37. doi: 10.1080/08958370701821482 18236227
47. Nemmar A, Hoet PM, Vanquickenborne B, Dinsdale D, Thomeer M, Hoylaerts M, et al. Passage of inhaled particles into the blood circulation in humans. Circulation. 2002;105(4): 411–4. doi: 10.1161/hc0402.104118 11815420
48. Furuyama A, Kanno S, Kobayashi T, Hirano S. Extrapulmonary translocation of intratracheally instilled fine and ultrafine particles via direct and alveolar macrophage-associated routes. Arch Toxicol. 2009;83(5): 429–37. doi: 10.1007/s00204-008-0371-1 18953527
49. Du Y, Xu X, Chu M, Guo Y, Wang J. Air particulate matter and cardiovascular disease: the epidemiological, biomedical and clinical evidence. J Thorac Dis. 2016;8(1): E8. doi: 10.3978/j.issn.2072-1439.2015.11.37 26904258
50. Xing Y-F, Xu Y-H, Shi M-H, Lian Y-X. The impact of PM2. 5 on the human respiratory system. J Thorac Dis. 2016;8(1): E69. doi: 10.3978/j.issn.2072-1439.2016.01.19 26904255
51. Mills NL, Donaldson K, Hadoke PW, Boon NA, MacNee W, Cassee FR, et al. Adverse cardiovascular effects of air pollution. Nat Rev Cardiol. 2009;6(1): 36. doi: 10.1038/ncpcardio1399 19029991
52. Landrigan PJ, Fuller R, Acosta NJ, Adeyi O, Arnold R, Baldé AB, et al. The Lancet Commission on pollution and health. Lancet. 2018;391(10119): 462–512. doi: 10.1016/S0140-6736(17)32345-0 29056410
53. Viana M, Kuhlbusch T, Querol X, Alastuey A, Harrison R, Hopke P, et al. Source apportionment of particulate matter in Europe: a review of methods and results. J Aerosol Sci. 2008;39(10): 827–49. doi: 10.1016/j.jaerosci.2008.05.007
54. Garcia CA, Yap P-S, Park H-Y, Weller BL. Association of long-term PM2. 5 exposure with mortality using different air pollution exposure models: impacts in rural and urban California. Int J Environ Health Res. 2016;26(2): 145–57. doi: 10.1080/09603123.2015.1061113 26184093
55. Simoni M, Baldacci S, Maio S, Cerrai S, Sarno G, Viegi G. Adverse effects of outdoor pollution in the elderly. J Thorac Dis. 2015;7(1): 34. doi: 10.3978/j.issn.2072-1439.2014.12.10 25694816
56. Roberts JD, Voss JD, Knight B. The association of ambient air pollution and physical inactivity in the United States. PLoS ONE. 2014;9(3). doi: 10.1371/journal.pone.0090143 24598907
57. Yu W, Sulistyoningrum DC, Gasevic D, Xu R, Julia M, Murni IK, et al. Long-term exposure to PM2. 5 and fasting plasma glucose in non-diabetic adolescents in Yogyakarta, Indonesia. Environ Pollut. 2020;257: 113423. doi: 10.1016/j.envpol.2019.113423 31677868
58. Weichenthal S, Villeneuve PJ, Burnett RT, van Donkelaar A, Martin RV, Jones RR, et al. Long-term exposure to fine particulate matter: association with nonaccidental and cardiovascular mortality in the agricultural health study cohort. Environ Health Perspect. 2014;122(6): 609–15. doi: 10.1289/ehp.1307277 24633320
59. Trivedi T, Liu J, Probst JC, Martin AB. The metabolic syndrome: are rural residents at increased risk? J Rural Health. 2013;29(2): 188–97. doi: 10.1111/j.1748-0361.2012.00422.x 23551649
60. Branis M, Linhartova M. Association between unemployment, income, education level, population size and air pollution in Czech cities: evidence for environmental inequality? A pilot national scale analysis. Health Place. 2012;18(5): 1110–4. doi: 10.1016/j.healthplace.2012.04.011 22632903
Článek vyšel v časopise
PLOS Medicine
2020 Číslo 6
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
- Co obsahují nové pokyny WHO pro dozor nad antibiotickou rezistencí?
- AI místo lékaře? Výzkum zkoumal ochotu Čechů konzultovat své zdraví s umělou inteligencí
- Prof. Jan Škrha: Metformin je bezpečný, ale je třeba jej bezpečně užívat a léčbu kontrolovat
- Ibuprofen jako alternativa antibiotik při léčbě infekcí močových cest
Nejčtenější v tomto čísle
- Distinct subtypes of polycystic ovary syndrome with novel genetic associations: An unsupervised, phenotypic clustering analysis
- Fatty acids in the de novo lipogenesis pathway and incidence of type 2 diabetes: A pooled analysis of prospective cohort studies
- The potential impact of COVID-19 in refugee camps in Bangladesh and beyond: A modeling study
- Chronic pain diagnosis in refugee torture survivors: A prospective, blinded diagnostic accuracy study