Article Text

Download PDFPDF

Linking environmental effects to health impacts: a computer modelling approach for air pollution
Free
  1. Jennifer Mindell1,
  2. Roger Barrowcliffe2
  1. 1Imperial College London, UK
  2. 2Environmental Resources Management, London, UK
  1. Correspondence to:
 Dr J Mindell
 Department of Epidemiology and Public Health, University College London, 1–19 Torrington Place, London WC1E 6BT, UK; j.mindellucl.ac.uk

Abstract

Study objective and Setting: To develop a computer model, using a geographical information system (GIS), to quantify potential health effects of air pollution from a new energy from waste facility on the surrounding urban population.

Design: Health impacts were included where evidence of causality is sufficiently convincing. The evidence for no threshold means that annual average increases in concentration can be used to model changes in outcome. The study combined the “contours” of additional pollutant concentrations for the new source generated by a dispersion model with a population database within a GIS, which is set up to calculate the product of the concentration increase with numbers of people exposed within each enumeration district exposure response coefficients, and the background rates of mortality and hospital admissions for several causes.

Main results: The magnitude of health effects might result from the increased PM10 exposure is small—about 0.03 deaths each year in a population of 3 500 000, with 0.04 extra hospital admissions for respiratory disease. Long term exposure might bring forward 1.8–7.8 deaths in 30 years.

Conclusions: This computer model is a feasible approach to estimating impacts on human health from environmental effects but sensitivity analyses are recommended.

Relevance to clinical or professional practice: The availability of GIS and dispersion models on personal computers enables quantification of health effects resulting from the additional air pollution new industrial development might cause. This approach could also be used in environmental impact assessment. Care must be taken in presenting results to emphasise methodological limitations and uncertainties in the numbers.

  • GIS, geographical information system
  • ED, enumeration district
  • air pollution
  • computer model
  • environmental exposure
  • geographical information system
  • health impacts

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Environmental impact assessment is a legal requirement in the European Union, as in many parts of the world, for a range of structural and industrial projects. Consideration of impacts on human health is required1–4 but is seldom done.2,5–7 Even when considered, health effects are seldom quantified.8

Many epidemiological studies show remarkably consistent associations between several health outcomes and exposure to increased concentrations of certain airborne pollutants, such as particulate matter (PM10, particles of aerodynamic diameter less than 10 microns).9–12 The availability of geographical information systems (GIS) and dispersion models on personal computers makes applying this knowledge a realistic possibility for quantifying health effects of the additional air pollution from a new industrial development. Although spreadsheets and other models have been produced to quantify health effects of specific scenarios,13,14 these have not been widely applied linked to GIS. By using a GIS, we are able to perform a large number of calculations for variables that display a large spatial variation. Without this tool, the task of aggregating the effects across such a wide area at this level of detail would be impossible.

In this case study, we have quantified some of the health effects of a new energy from waste facility on the surrounding population, in the context of environmental assessment for planning and licensing of the proposal.

METHOD

Our method was a development based on an existing computer spreadsheet model.13,15 That model was based on changes in ambient pollution at monitoring sites resulting from policies to improve air quality in a local authority, assuming that levels changed in parallel across the whole area. This new method links the spreadsheet to a GIS, so that the health impact calculations for the modelled additional exposure experienced by the resident population of each enumeration district can be made. (The enumeration district (ED), containing around 100 residents, was the smallest output level in the UK’s 1991 census. It has been superseded by output areas from the 2001 census.)

The quantification relies on the simple equation:

ΔE  =  β×ΔC×P×E,

where: (Δ)E  =  (change in) background rate of events;

β  =  exposure-response coefficient;

ΔC  =  change in concentration of pollutant;

P  =  population exposed.

The coefficients are derived from epidemiological studies, which show a clear association between increased exposure and increased effects. We have included only those pollution-health outcomes associations for which the evidence of a causal relation is at least moderately strong, based on systematic literature reviews13 (table 1).

Table 1

 Strength of evidence for causal relations for acute exposure to ambient levels

When assessing potential health impacts of exposure to raised levels of criterion air pollutants from the proposed plant, we examined PM10, sulphur dioxide, and nitrogen dioxide. As this paper aims to describe the quantification method of a computer model combining geographical information with calculations, we have restricted description of the methods and findings to PM10 as an exemplar.

There is strong evidence that exposure to PM10 brings forward non-traumatic deaths9,13 and moderate evidence13 that PM10 increases emergency hospital admissions for circulatory17 and respiratory17 diseases and accident and emergency attendances for asthma.19 There is also moderate evidence that long term exposure to particulate matter increases deaths.13,20 Within the range of ambient pollution encountered in the region, the health effects of particulate matter is linear and there is no good evidence of a threshold below which they have no effect on any person.15,21,22 This assumption of no threshold is made based on the epidemiological evidence but results in being able to use the modelled increase in annual average concentration rather than all the individual 24 hour average concentrations that make up a year.

The anticipated emissions from the proposed incinerator were entered into the ADMS air pollution dispersion model23,24 to obtain “contours” of additional concentrations for the new source. This output file was easily transferred into the GIS. A population database at ED level for a 20 km radius by age group (0–14, 15–64, 65+) was added to the GIS as another data layer.

The pollutant concentrations are represented by the dispersion model as lines. These were used to create small area polygons in the GIS, each defining an area with a single value for concentration. The polygons were formed by reference to the ED boundaries, so that each polygon had one total population value, with the assumption that the population was uniform across the whole area. Using an overlay function in GIS, the population data were combined with the PM10 concentration data, resulting in a single layer containing details of ED name, spatial area, population, and a PM10 concentration for many hundreds of non-uniform areas.

Using the GIS it was possible to see, at a glance, the population density variations across the study area and in conjunction with the distribution of pollutant concentrations. It was a straightforward exercise to add new data into the GIS and run the model again, for a larger area for example. Similarly, the model could be run with different pollutant concentration fields.

The results from the GIS were then exported into MSExcel so that further calculations could be made.

Aggregated data on current regional rates of mortality and hospital admissions for several causes from the Office for National Statistics (ONS) and Hospital Episode Statistics databases respectively were obtained from the Compendium of Clinical and Health Indicators 2002/Clinical and Health Outcomes Knowledge Base and the London Health Observatory. We used a three year average to reduce the effects of random variation.

The GIS is linked to a spreadsheet containing the effect estimates and the background rates for the various health outcomes.

A number of uncertainties exist in the information underlying this assessment. We therefore included three types of sensitivity analyses. Precise effect estimates from relevant, well conducted primary studies and from meta-analyses are of the same order of magnitude but vary (table 2), so the analyses were repeated using a number of coefficients. The analyses were repeated using a 15 km radius for impacts, to see the effect of changing the geographically determined population for which health impacts are assessed. Thirdly, for a “worst case scenario” approach, we used the event rates from the local government area within the 20 km circle that had the highest background rates of deaths and hospital admissions.

Table 2

 Estimates used of percentage change in mortality and hospital admissions for the effect of an acute change of 1 μg/m3 PM10

RESULTS

The modelled changes in annual mean concentrations of PM10 were small, being 0.08 µg/m3 at the most affected location (fig 1; a colour version of the figure is available on line at http://www.jech.com/supplemental). The average additional exposure over the area within 20 km would be 0.002 μg/m3. Even at the most affected location, concentration increases of this magnitude would be impossible to detect through the use of monitoring instruments, given current background concentrations of about 25 µg/m3 and the precision of current instruments (±2 μg/m3).

Figure 1

 Additional mean annual PM10 levels from the proposed incinerator (20 km radius).

At baseline, the 3 689 826 population living within 20 km of the proposed plant experienced an average of 29 984 non-traumatic deaths, 28 541 cardiovascular and 27 345 respiratory admissions per year. Table 3 shows the estimated number of deaths brought forward by operation of the facility.

Table 3

 Acute effects on premature deaths of extra PM10 within 20 km of the proposed facility

Using a 20 km radius, the additional PM10 emitted from the proposed plant was estimated to bring forward about 0.03 deaths per year—that is, 0.9 earlier death, including 0.39 earlier respiratory death, every 30 years, corresponding to 0.0001% of all deaths.

Table 4 shows extra or earlier hospital admissions. There would be about 1.3 extra or earlier emergency respiratory admissions, and 0.9 extra or earlier emergency cardiovascular admissions in 30 years, <0.0002% of all respiratory and 0.0001% of all circulatory admissions.

Table 4

 Acute effects on admissions and consultations of extra PM10 within 20 km of the proposed facility

Over a 30 year period, there would be 0.08 extra or earlier hospital admission among children aged less than 15 years and 0.17 in adults aged 16 to 64 for asthma. Visits to accident and emergency departments for asthma could not be calculated, as neither baseline numbers nor rates were available.

Long term effects

There is less evidence for the long term effects of particulate exposure. Three American33–35 and one Dutch36 cohort studies have found associations with increased mortality. Re-analysis adjusting for a wide range of sociodemographic and other potential confounder confirmed a probable effect, although with smaller effect estimates.37

The effect is probably causal but the magnitude of effect has not been measured in England. The English Committee on the Medical Effects of Air Pollution (COMEAP), stated in 200138 that the most probable effect estimate was the lower limit reported by the Health Effects Institute reanalysis –0.1% per 1 μg/m3 PM2.5.37

The long term effects of 30 year exposure to the estimated additional PM10 over the area within 20 km of the proposed plant would cause between an extra 1.8 deaths in 30 years (using as our lower limit the COMEAP “most likely” effect estimate of 0.1% per 1 μg/m3 PM2.5 and applying this effect estimate to the PM10 levels) and 7.8 deaths in 30 years (using the weighted average estimate from the WHO three country study39).

Sensitivity analyses

Table 5 summarises the results of sensitivity analyses for which the effect estimate, baseline mortality rate, or area affected were varied.

Table 5

 Examples of sensitivity analyses for increases in annual health outcomes

The most important effect on health is the increase in non-traumatic deaths. Different exposure effect estimates gave the range for the estimated acute impact on mortality as 0.55–2.91 earlier deaths every 30 years.

The primary calculation used the region’s average death rate, but there is pronounced geographical variation in this rate. Using the highest all cause death rate experienced in a borough within 20 km of the proposed plant as a worse case scenario had a similar effect to increasing the exposure effect estimate: it increased the estimate of the number of deaths brought forward to 1.26 in 30 years.

Considering the population living within 15 km (2 214 151) rather than 20 km of the proposed plant reduced the number of deaths brought forward by the additional PM10 emitted over a 30 year period from 0.91 to 0.72.

DISCUSSION

Exposure response coefficients from American and European epidemiological studies apply best to those populations and pollution types for which the studies were done. These have tended to be urban populations where pollution is influenced primarily by road traffic pollution. The influence of industrial pollution in urban areas is much less than it was several decades ago. For particles, there is be more reason to consider the source of pollution, as particulate matter is usually a composite of many chemicals, although the consistency of association between outcomes and mass of PM10 from studies in different parts of the world is striking.

Some factors that may increase the uncertainty are:

  • The temporal nature of the pollution. Industrial plumes expose a downwind population to intermittent, short periods of higher concentrations. Road traffic pollution in urban areas is more uniform, with occasional episodes of increased concentrations caused by unfavourable weather conditions.

  • The socioeconomic circumstances of the two population groups. In our example, the exposed population is mainly urban and reflects well the populations used in epidemiological studies of pollution mainly originating from transport.

The English Department of Health’s expert committee on health effects of air pollution (COMEAP) has confirmed that it is appropriate to use the results of time-series studies from transport pollution for estimating effects of industrial pollution but has reported on the additional uncertainties.40 In line with their recommendations, we have been explicit about assumptions made.

The more widely practised methods of quantifying health outcomes from exposure to air pollutants in environmental assessment have usually entailed a source-pathway-receptor approach, in which each various exposure route is treated as a series of compartments thorough which the pollutants are transferred. A large number of pollutants are amenable to this treatment and their effects on humans are calculated by reference to toxicological data. This approach is popular in North America, where the concept of risk assessment is embedded in regulatory practice.41,42 The approach we outline here is fundamentally different, not least with respect to the origin of the data relating to harm. The North American approaches to the evaluation of health risk rely on the use of large toxicology databases for an array of chemicals. Information on individual chemicals is typically based on animal and occupational exposure data. For carcinogens, a non-threshold approach is adopted and so a theoretical risk extends to all non-zero exposures, however small. For non-carcinogens, a reference dose is used, which signifies the point at which harm is assumed to occur. Often, the technique is applied to defined human receptors that are taken to represent specific exposure scenarios, such as a farmer or child. Significant parts of the modelling elements of the risk assessment approach are concerned with the estimation of exposure, either to inhalation or ingestion, perhaps through the food chain.

What this paper adds

  • Combining a GIS system, dispersion modelling, and a spreadsheet model to quantify potential health impacts of new or existing sources of pollution is feasible.

  • Using such an approach, sensitivity analyses are easily performed, to examine the effects of changing the underlying assumptions.

In contrast, the approach we describe here is limited to a small number of pollutants, but the basis of the estimation of harm is the reported effect from epidemiological studies on the wider human population. It is also not a measure of the risk of harm occurring, but rather an attempt to quantify an outcome assumed, from an assessment of causality from a systematic literature review, to occur. We assume that the exposure of the affected population is well represented by the dispersion model and we include only inhalation of airborne pollutants.

The California Environmental Protection Agency approach, and its HARP (hotspots analysis reporting program) software, can be linked to a GIS but assumes a threshold relation for non-carcinogenic air pollution effects.41 The impacts of a proposed development are therefore based on calculating the health impacts of air pollution with and without the proposal, the impacts of the development being the difference between the two. Because of the assumption of a threshold, annual averages cannot be used but each day’s particulate levels (or hourly levels for some pollutants) are required to be computed separately and the effects summed.

Our computer model uses a similar theoretical approach to quantifying health impacts as Health Canada’s revised risk determination framework.43 We have produced a tool to estimate these impacts.

The mathematical model developed by Biddulph and colleagues14,44 uses more detailed population dynamics than our model but in the absence of more detailed age specific exposure-effect estimates, we decided not to use population data broken down into more than the three age bands. If the population in the potentially affected areas within each of the adult age groups were older than average, a higher baseline number of events and therefore greater impact would have been found by using narrower age specific baseline event rates. Their model also incorporates baseline health status and models transition between these for a “no intervention” and a “proposed project” scenario, the health impact of the proposal being the difference between these. However, such detailed population and health status data are seldom available to those conducting environmental health impact assessments, nor have the time series studies on which the exposure-effect estimates are based used such data. Their model also does not seem to incorporate geographical information.

Policy implications

Where environmental effects are of concern because of their potential impacts on health, these should be identified and, where possible, quantified.

Like Health Canada’s approach to estimating health benefits,45 but unlike the tools mentioned above, our approach uses the change model46 rather than a burden of disease approach.46,47 This has the advantage of ignoring factors that do not change, thus simplifying the calculations. This is particularly important for aspects such as baseline air pollution levels, which are generally available only at a few specified locations. A model that requires calculation of the effects both of baseline pollution exposure and of exposure to the “baseline plus effects of proposal” would therefore require spatial modelling of pollutant levels at a small area resolution if modelled exposure to pollution from a proposed development is to be estimated accurately.

This approach is also applicable to other impacts, provided the change in “exposure”, the magnitude of effects of a unit change, and the baseline rates are known, and to other sources, if changes in exposure can be modelled. It could also be used to model differential effects on susceptible populations if effect estimates for different subgroups became available, together with the relevant population data. This model does already allow for the tendency for more deprived populations to be exposed to greater concentrations of pollutants because of its fine geographical resolution in modelling air pollution and its consequent health impacts. The additional refinement of using area specific baseline mortality and morbidity rates would add to the assessment of equity of effects.

It is possible to quantify the impact on some health outcomes of exposure to some additional levels of pollution, caused by a new point source. Using a GIS in this way can result in a more precise estimate, as the modelled exposure for each small area can be related to the resident population of that area, rather than applying a uniform average or maximum exposure to all. Our approach might be particularly suitable for the case of air quality changes introduced by a linear pollutant source, such as a road, or changes across a road network on the scale of a city, provided that such air quality changes can be captured within a GIS and matched with a population database at the appropriate spatial resolution.

The main advantage of this computer model approach is the ease with which the input data can be varied to reflect changing knowledge, such as more recent effect estimates48,49; other effect estimates to enable comparisons with those used elsewhere9; other underlying assumptions (for example, normal operating conditions or worse case scenarios); or the level of uncertainty in the data (and therefore the need to conduct sensitivity analyses). When more than one pollutant is considered, effect estimates obtained from analyses using two-pollutant models can be substituted.16

A similar model was used in the past to ascertain whether calculations using annual changes in sulphur dioxide differed from those using monthly changes.50 No significant difference in conclusions were found, so that empirical study confirms our theoretical basis for using annual changes for a pollutant with no threshold of effect.

This method can also show the effects of varying the geographical parameters within which exposure is defined. Eighty per cent of the potential health effects seen when a 20 km radius was used occur within 15 km of the proposed plant, although only 60% of the larger population live within this smaller area. Extending the area further would be unlikely to have more than a minimal effect.

Simple calculations of health effects of modelled changes in air pollution have been undertaken in the past, examining the benefits of improvements in air quality13 or potential hazards from increased pollution from a development.51 Our computerised approach permits more sophisticated modelling. Whereas those studies assumed that all the residents of the selected area would be exposed to the same estimated concentrations, the basis of our method is to model air pollution concentrations at a fine geographical resolution and use GIS to ascertain the effects of the modelled concentration of the pollutant for each enumeration district on the population living within that small area.

Additionally, our computer model enables different background rates to be applied to each area for more precise estimates. In this case, we used localised population and modelled air pollution data but did not incorporate specific background mortality or morbidity rates for each local authority area as the numbers were so small that using the regional average and the highest local authority rate within the affected area were deemed sufficient analysis.

The estimates are a useful order of magnitude guide to some health impacts. Application of the method in this case was able to show that the magnitude of quantifiable health effects that might result from the increased exposure to PM10 is small—about 0.03 deaths each year in a population of 3 500 000, with 0.04 extra hospital admissions for respiratory disease.

In comparison with the more deterministic modelling approaches to calculating health impact, our approach is restricted to those pollutants for which we have strong evidence of causality through epidemiological studies. However, these tend to be the pollutants associated with combustion processes and so are those of primary interest. In addition, particulate matter is not a pollutant that is easily treated by reference to toxicological data derived from occupational exposure or animal studies and its effects on the human population could be missed by such an approach.

This is a useful technique but is only a part of identifying potential health impacts. A similar model can be used to quantify any causal impact if the exposure can be modelled and the baseline event rate, susceptible population, and the magnitude of a linear exposure-effect relation are known.13 It can also be used to compare potential health impacts from different proposals as part of options appraisal or within an economic assessment. However, it cannot be used to compare the magnitude of health impacts from different proposals or caused by different aspects of a proposal (for example, from transport of waste, effects on visual amenity, noise) unless all the elements to be compared are quantifiable.13

The extent to which deaths are brought forward is debated.52–55 The long lasting effects of long term exposure to lower levels of pollution are generally held to be considerably larger than the acute effects.56 However, the effects of long term exposure to pollution are still an order of magnitude less than the premature deaths caused by smoking, where each smoker whose death is caused by tobacco loses an average of 10 years of life.57

Some care must be taken in presenting these results to emphasise the uncertainties in the numbers and the limitations of the method. Nevertheless, for decision makers who might be unclear on the scale of possible impacts, the methodology is a useful tool. For the public, context is crucial: for them even one premature death might be seen as an unnecessary burden from a new source appearing in their community.

Acknowledgments

We thank Kate Murray for her GIS work for this project; Dr Mike Joffe for help with developing the model and commenting on earlier drafts of this paper; and Dr Bob Humphries for helpful comments.

REFERENCES

Supplementary materials

  • Figure 1 in Colour

    Files in this Data Supplement:

Footnotes

  • Funding/competing interests: JM and RB were funded by the company that was proposing to build the waste facility.