Objectives First, we present a general analytical approach to estimating the association between medium-term changes in air pollution and health across small areas. As a specific illustration, we then applied the approach to data on London residents from a 4-year period to test whether reductions in traffic-related air pollution were associated with reductions in cardio-respiratory hospital admissions.
Methods A binomial distribution was used to model change in admissions over time in each small area, which was measured as the proportion of admissions in 2003–2004 out of admissions over all study years (2001–2004). Annual average concentrations of nitrogen oxides (NOx) were modelled using an emissions-dispersion model. The association between change in NOx and change in hospital admissions was estimated using logistic regression and an instrumental variable approach.
Results For some diagnostic groups, suggestive associations between reductions in NOx and reductions in admissions were observed, for example, OR=0.97 (95% CI 0.96 to 0.99) for an IQR decrease in NOx (3 μg/m3) and all respiratory admissions. Accounting for spatial dependence attenuated several of the associations; for respiratory admissions, the OR was 1.00 (95% CI 0.98 to 1.02), leaving only that for bronchiolitis significant (OR=0.91; 95% CI 0.84 to 0.99). In this particular illustration, the instrumental variable approach did not appear to add information.
Conclusions In this illustration, there was relatively limited power to detect an association between changes in air pollution and hospital admissions over time. However, the analytical approach could deliver more robust estimates of the health effects of changes in air pollution in settings with greater spatial contrast in changes in air pollution over time.
- Air pollution
- mathematical models
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.
Statistics from Altmetric.com
Ambient air pollution has decreased over the past several decades throughout much of the developed world,1 2 largely due to increasingly stringent air pollution controls. While there is substantial evidence of an association between air pollution and adverse health outcomes,3 4 direct evidence regarding the extent to which control measures have improved public health is limited.5 Whether past efforts to control air pollution have resulted in measurable improvements in health, and the likelihood that future efforts will continue to do so, is of great interest.5
Although several studies have observed an association between short-term variation in air pollution and hospitalisations for cardio-respiratory disease,6–9 few studies have investigated associations with long-term time trends in air pollution. Associations between long-term trends in air pollution and health are frequently confounded by slowly time-varying factors, which are often difficult to measure.10 Recent efforts to estimate the association between longer term trends in air pollution and health have therefore applied analytical approaches that specifically aim to minimise confounding.11 12 One such approach has been to estimate the association of local trends in air pollution and health outcomes, where time-varying factors would also have to vary across areas in order to be confounders.12–14 Such approaches should be less vulnerable to confounding by time trends that have little variability across areas as well as by confounders of purely spatial contrasts in exposure, such as smoking, which change little over time.
This paper has two objectives. First, we present a general analytical approach to estimating the association between changes in air pollution and changes in health outcomes over time across small areas. As one specific illustration, we apply the approach to data on traffic-related air pollution and cardio-respiratory hospital admissions from a 4-year period surrounding the introduction of the London Congestion Charge Scheme. We selected this period because the scheme may have led to changes in traffic-related pollutants in only some areas of the city. The scheme, introduced in February 2003, charges drivers to travel within a specified zone. The second objective is to test the hypothesis that areas with larger decreases in traffic-related air pollution had larger decreases in cardio-respiratory hospital admissions over time.
Analysis was based on Hospital Episodes Statistics data, which record all episodes of care provided by NHS hospitals in England.15 We included patients resident in Greater London who were admitted as an emergency with cardio-respiratory disease as the primary diagnosis during 2001–2004. As a patient may have more than one episode of care during an admission, we restricted analysis to first, finished episodes. To remove duplicates, we included only the first record among those with an identical patient identifier, health care provider code, date of admission and discharge. For records with the same identifier and date of admission, but different date of discharge, we used only the record with the later date of discharge. Admissions within 28 days of discharge from the previous admission for the same patient were excluded since these were likely to be readmissions due to complications rather than independent admissions. We considered admissions for the following primary diagnoses according to ICD-10 codes: all cardiovascular diseases (CVD) (I), ischaemic heart disease (IHD) (I20–I25), stroke (I60–I69), heart failure (I50), all respiratory (J), chronic obstructive pulmonary disease (COPD) (J43–J44), asthma (J45–J46) and bronchiolitis (J21).
Our outcome was the proportion of admissions occurring in 2003–2004 out of admissions in 2001–2004. This proportion can be considered as a measure of change in admissions over time: a proportion of 0.5 represents no change and greater than 0.5 represents an increase over time.
Air pollution exposure
We considered nitrogen oxides (NOx) as a surrogate for traffic-related air pollution. Annual average concentrations of NOx were modelled for each year from 2001 to 2004 using a hybrid emission-dispersion/regression modelling approach. The model is described in detail elsewhere.16 Briefly, ADMS Roads and OSPM were used to model the NOx contribution from roadways within a 500 m buffer around 31 monitoring locations. ADMS3 was used to model the NOx contribution from urban background pollution. The ADMS models are second generation Gaussian dispersion models17 18; OSPM is a street canyon model.19 Model input included detailed data on meteorology, traffic flow, speeds and vehicle type.20
As a calibration, we then fit a regression model to measurements of NOx from 31 monitoring sites. The model included two predictors: the contribution of NOx from roadways and the urban background as estimated from the emission-dispersion modelling. The regression model was then used to predict NOx concentrations on a uniform 20 m×20 m grid. The Spearman correlation between modelled difference in annual means between 2001–2002 and 2003–2004 and measured difference at 24 NOx monitors not included in the calibration step was 0.80.
Annual average NOx concentrations for each year and ward were calculated as the mean of NOx concentrations across model grid points falling within each ward's boundary. Change in ward-level NOx was defined as the mean concentration in 2001–2002 minus that in 2003–2004. Wards within or adjacent to the boundary of the congestion charging zone (CCZ) were characterised with a binary CCZ indicator.
Income deprivation was based on the Index of Multiple Deprivation (IMD) calculated at the level of lower layer super output areas (LSOAs), geographical units with a mean population of 1500.21 IMD is a composite score of deprivation made up of seven domains: income, employment, health and disability, education, skills and training, barriers to housing and services, living environment and crime. The income domain is measured as the proportion of residents living in income deprivation; therefore, calculating change over time is more appropriate for the income domain than for others which are based on relative scores.22 IMD and income deprivation at the ward level were calculated as population-weighted averages of LSOA deprivation using weights based on 2001 population estimates.23 We then calculated a measure of change in deprivation: the proportion of ward residents living in income deprivation in 200124 minus the proportion in 2005.25
We constructed demographic measures that reflected ward-level population at a single time as well as changes over time. We calculated percentages of ward residents by ethnicity in 2001 using the Census 2001 data.26 Changes in ward-level population were calculated as the ward mean of 2001–2002 minus that of 2003–2004 for residents aged 65 and older, younger than 20 and for the total population.23
To estimate the association between change in NOx and change in admissions, we used a binomial distribution to model the number of admissions in 2003–2004 as a proportion of admissions over all study years (2001–2004) by ward (analogous to modelling the number of events out of total trials). The association between the log odds of 2003–2004 admissions with change in NOx between 2001–2002 and 2003–2004 was estimated using Proc Logistic in the SAS software package. We adjusted for potential confounders that were independent predictors of change in admissions and correlated with changes in NOx: average IMD in 2001, change in proportion of residents living in income deprivation, change in number of residents aged 65 years and older, and percentages of residents who were non-white, black, Asian and Chinese. All covariates were continuous and modelled with a linear functional form.
The unit of analysis in the logistic model was the ward (n=628), a census based geographical unit that is nested within London boroughs (n=33). Wards within the same borough are likely to have more similar hospitalisation rates compared to wards in other boroughs, which may result in incorrect effect estimates and standard errors. We therefore conducted an analysis to account for this spatial dependence by including a random intercept for each borough using Proc Glimmix. To describe the extent of such clustering we also estimated the intraclass correlation, which indicates the proportion of the total variability in the outcome at the borough level.
Instrumental variable analysis
We further evaluated whether charging zone status could be used effectively as an instrumental variable to test the hypothesis stated in our second objective. An instrumental variable estimate may be less vulnerable to confounding than the association using modelled change in NOx as exposure if several assumptions are met.27 These assumptions are primarily that (1) the instrumental variable is associated with changes in NOx, (2) the instrumental variable affects changes in admissions only through changes in NOx and (3) the instrumental variable is not associated with other unmodelled predictors of change in admissions. We evaluated whether each assumption was reasonable using charging zone status as an instrumental variable. Specifically, we used a two-sample t test to test the association between the CCZ indicator and change in NOx. Assumptions 2 and 3 are not empirically verifiable; therefore, we evaluated them indirectly. We assessed the strength of association between the CCZ indicator and changes in admissions that was not mediated by changes in NOx. To evaluate the third assumption, we compared observable socio-demographic variables for wards within and outside the charging zone.
We then used a 'pseudo two-stage least squares' method to estimate the instrumental variable association between changes in NOx and admissions. In the first stage, we regressed change in NOx on the CCZ indicator (and potential confounders in the fully adjusted model) to estimate an adjusted mean difference in change in NOx inside minus outside the zone. In the second stage, we regressed the proportion of admissions in 2003–2004 on predicted values from the first stage (and potential confounders in the adjusted model). Although in linear models, instrumental variable estimates are consistent, with non-linear link functions they may be biased.28 Nonetheless, the bias is often small and two-stage logistic methods have been shown in practice to provide similar results as generalised method of moments estimation.29 In this application, the proportions predicted by the logistic models were near 0.5 (5th percentile: 0.48; 95th percentile: 0.55); therefore, the link function is very nearly linear and the bias is expected to be small.
Across London, NOx concentrations decreased by 3.5 μg/m3 between the two time periods (table 1). The average NOx concentration in 2001–2002 was 77 μg/m3 (SD 17 μg/m3). Across wards, there was relatively modest variation in decreases in NOx over time; the difference between the 90th and 10th percentile was 4.7 μg/m3. NOx concentrations increased in <5% of wards, particularly in the western half of the city (figure 1). For wards within or adjacent to the CCZ boundary, NOx decreased by 5.8 μg/m3. On average, there was an unadjusted increase in admissions for respiratory disease, COPD and asthma (table 1). However, there appeared to be little or no unadjusted change over time for other cardio-respiratory disease categories. The greatest variation across wards in the proportion of hospital admissions occurring in 2003–2004 was observed for bronchiolitis (90th–10th percentile=0.44) and heart failure (90th–10th percentile=0.26).
Table 2 presents correlations between decreases in NOx concentration, deprivation, population and ethnicity. Wards with the largest decrease in NOx had increasing income deprivation over time, higher deprivation in 2001, a decreasing number of residents aged 65 years or older and a higher percentage of black residents.
The relative odds of hospital admissions occurring in 2003–2004 per interquartile range (IQR; 2.66 μg/m3) decrease in NOx is presented according to diagnostic group in table 3. After adjusting for potential confounders, an IQR decrease in NOx was associated with a 9% (95% CI 2% to 15%) decrease in the odds of admission for bronchiolitis in the main analysis. Similarly, decreasing NOx was associated with decreasing admissions for respiratory disease and heart failure. Further adjustment for change in total population or in those aged 20 or younger between 2001–2002 and 2003–2004 did not alter the results of the main analysis. However, accounting for spatial dependence at the borough level substantially attenuated the associations. Of the fully adjusted estimates, only the association for bronchiolitis, which had the least spatial dependence, remained statistically significant. Much of the total variation in the proportion of admissions in 2003-04 was between boroughs, particularly for all cardiovascular and respiratory admissions.
The average change in NOx concentration differed significantly between the CCZ and non-CCZ wards: mean difference was 2.5 μg/m3 (SD 1.8), p<0.0001 from t test. However, the R2 was only 0.10 in a model regressing change in NOx on charging zone status, indicating that the instrumental variable was weakly associated with exposure. Furthermore, only for bronchiolitis was charging zone status significantly associated with the proportion of admissions in 2003–2004 after adjusting for confounding variables (table 4). This association was attenuated to non-significance after including change in NOx, providing little evidence against the assumption that the instrumental variable was associated with the outcome only through its association with exposure. Although we could not directly test whether there were unmeasured differences between the CCZ wards and the rest of London, we observed significant differences between wards within and outside the zone for several observable socio-demographic variables (table 5). For example, on average, wards within the zone were more deprived in 2001 and had a larger increase in total population over time. The instrumental variable estimates were broadly comparable to the estimates from the main analysis but had much wider confidence intervals (table 3).
We presented a general analytical approach to investigate whether spatial variation in medium-term changes in air pollution were associated with changes in the occurrence of a health outcome. This approach is particularly useful for outcomes which are counts (eg, hospital admissions or deaths), since the difference between counts may be negative and therefore not Poisson distributed. By comparison, changes over time in normally distributed outcomes can be modelled in a relatively straightforward fashion as the difference between the outcome at two points in time.
In the main analysis, wards with larger decreases in NOx concentration had, on average, larger decreases (or smaller increases) in admissions for heart failure, respiratory disease and bronchiolitis after adjusting for potential confounders. However, there was substantial spatial dependence in the data, and only the adjusted association for bronchiolitis remained after accounting for the dependence. That we observed a significant association consistent with our hypothesis for bronchiolitis, but not other disease groups, may reflect the greater variation in our outcome for bronchiolitis and thus better statistical power. Whether the association for bronchiolitis could be causal rather than chance requires further investigation.
While studies focussing on contrasts in local changes in exposure may be less vulnerable to confounding by time-invariant factors than purely spatial contrasts in exposure, confounding by time-invariant factors is still a concern. For example, in our data, the correlation between IMD for 2001 and annual average ward-level NOx in 2001 was 0.53, while the correlation with change in ward-level NOx was 0.27. Whether time-invariant factors are proxies for other variables that change over time or are important in their own right requires further investigation. We observed large differences in several socio-demographic variables between wards within versus outside the charging zone (table 5), indicating that wards with the largest decrease in NOx are different in many respects from wards with smaller decreases in NOx. Although we were able to adjust for several socio-demographic variables, and used an analysis specifically designed to reduce confounding, residual confounding of our estimated odds ratios by unmeasured socio-demographic variables cannot be ruled out.
In our illustration, unmeasured changes over time in individual hospitals' admission practices are also likely to be important determinants of the outcome since they would affect a given hospital's catchment area rather than London overall. Applying the proposed approach to other routinely collected health outcomes, such as mortality, would avoid possible confounding due to this source. Our choice to investigate hospital admissions rather than deaths was based on power considerations, as there are approximately 70 000 more hospital admissions per year in London compared to deaths from cardio-respiratory disease.
In the specific context of the congestion charge, the introduction of the scheme did not fully meet the assumptions necessary to be an effective instrumental variable. Our analysis indicated that charging zone status was associated with exposure, although the association was relatively weak. There were significant differences between several observable socio-demographic variables in wards within and outside the zone. Although not a formal test, this suggests that wards in the charging zone are different from remaining wards in several respects and might be so for factors for which we did not have data. The third assumption required for an instrumental variable to deliver less biased estimates would not be met if such unmeasured factors were associated with changes in admissions. A further limitation of the instrumental variable estimates is their relatively large standard errors, a common feature of the method and a reflection of the weakness of the association between the charging zone and change in NOx.30
Several limitations of the specific illustration we chose warrant consideration. We selected data from a period surrounding the introduction of a traffic management scheme that may have led to changes in pollutants in only some areas of London in order to maximise spatial contrast in changes in traffic-related air pollution. Nonetheless, the variations in modelled changes in NOx were small, limiting our power to detect associations. Modelled changes in exposure were also small relative to modelling error. The timing between a reduction in long-term exposure to air pollution and a possible reduction in the occurrence of cardio-respiratory hospitalisations remains uncertain; the 4-year period we considered may not be sufficiently long.
Relatively few studies have investigated the impact of decreasing medium- to long-term trends in air pollution on health because of the potential for these associations to be confounded by slowly time-varying factors. These few studies have developed or applied analytical approaches that specifically minimise confounding. To date, none have evaluated hospital admissions. Janes et al decomposed the association between trends in fine particles and mortality in US counties into the association between (1) national trends and (2) county-specific trends.12 The second, which indicates whether counties with steeper declines in particles also have steeper declines in mortality, was expected to be less confounded. Strong associations were observed between national trends in particles and mortality, but no association was observed for the county-specific trends. Illustrating another approach, Chay et al used non-attainment status designated under the 1970 Clean Air Act as an instrumental variable.11 Non-attainment status of the county was associated with reductions in total suspended particles, although no association with reductions in adult mortality was observed. Other studies investigated whether cities with the largest decrease in particles had the largest reduction in premature mortality,14 31 arguing that long-term trends in mortality were unlikely to explain their results since such trends would affect all cities, rather than have a differential effect in cities with larger reductions in pollution.
In conclusion, we developed an analytical approach to assess the association between changes in air pollution and health outcomes over time across small areas. In this illustration, which was restricted to a relatively short follow-up period and single metropolitan area, we had limited statistical power. We observed a significant association between reductions in NOx and reductions in admissions only for bronchiolitis. Nonetheless, the analytical approach could deliver more robust estimates of the health effects of changes in air pollution if data were available over a longer period or with greater contrasts in changes in air pollution.
What this paper adds
Few studies have evaluated whether reductions in air pollution on the scale of years are associated with improved population health.
This study provides a methodological approach for estimating the association between changes over time in air pollution and health outcomes using a binomial distribution to model changes over time in outcomes that are counts.
We provide a specific illustration of the approach, in which we test whether reductions in traffic pollutants over a 4-year period are associated with reductions in hospital admissions in London.
Funding This work was supported by the U.S. National Institute of Environmental Health Sciences (F32ES015432), and the Health Effects Institute.
Competing interests None.
Provenance and peer review Not commissioned; externally peer reviewed.