Using Geographic Information Systems To Assess Individual Historical Exposure to Air Pollution from Traffic and House Heating in Stockholm.

A specific aim of a population-based case-control study of lung cancer in Stockholm, Sweden, was to use emission data, dispersion models, and geographic information systems (GIS) to assess historical exposure to several components of ambient air pollution. Data collected for 1,042 lung cancer cases and 2,364 population controls included information on residence from 1955 to the end of follow-up for each individual, 1990-1995. We assessed ambient air concentrations of pollutants from road traffic and heating throughout the study area for three points in time (1960, 1970, and 1980) using reconstructed emission data for the index pollutants nitrogen oxides ([NO.sub.x]/[NO.sub.2]) and sulfur dioxide together with dispersion modeling. [NO.sub.2] estimates for 1980 compared well with actual measurements, but no independently measured (study-external) data were available for [SO.sub.2], precluding similar validation. Subsequently, we used linear intra- and extrapolation to obtain estimates for all other years 1955-1990. Eleven thousand individual addresses were transformed into geographic coordinates through automatic and manual procedures, with an estimated error of [is less than] 100 m for 90% of the addresses. Finally, we linked annual air pollution estimates to annual residence coordinates, yielding long-term residential exposure indices for each individual. There was a wide range of individual long-term average exposure, with an 11-fold interindividual difference in [NO.sub.2] and an 18-fold difference in [SO.sub.2]. The 30-year average for all study subjects was 20 [micro]g/[m.sup.3] [NO.sub.2] from traffic and 53 [micro]g/[m.sup.3] [SO.sub.2] from heating. The results indicate that GIS can be useful for exposure assessment in environmental epidemiology studies, provided that detailed geographically related exposure data are available for relevant time periods. Key words: air pollution, combustion, epidemiology, exposure assessment, GIS, lung cancer, traffic. Environ Health Perspect 109:633-639 (2001). [Online 15 June 2001]

http://ehpnet1.niehs.nih.gov/docs/2001/ 109p633-639bellander/abstract.html

Recent developments of geographic information systems (GIS) have provided environmental epidemiologists with new tools to study associations between environmental exposures and disease. A GIS is a powerful computer mapping and analysis tool, which permits spatial linking of different types of data (e.g., residential addresses, environmental exposure levels, demographic information) (1,2). GIS can also be used for data management, such as automated address matching to geographic coordinates (geocoding).

One important GIS application is the mapping of environmental exposure (3,4). Modeling of air pollution data has been used to assess ambient concentration fields, usually considered proxies of exposure fields for residents in the area (5,6). Dispersion models are physical (deterministic) models that use existing data on emissions and meteorologic and topologic conditions to create maps of pollutant concentrations (7). Such models are generally based on Gaussian plume dispersion equations. Several models have been constructed, such as the BREEZE model developed for the U.S. Environmental Protection Agency, the CAR model developed by the Dutch Environmental ministry (4), and the AIRVIRO models developed by the Swedish Meteorological and Hydrological Institute (8-10).

It has long been argued that exposure to urban air pollution contributes to the development of lung cancer, although the risk increase supposedly is limited (11-13). Studies from many countries indicate a smoking-adjusted relative risk in urban areas over countryside areas of up to 1.5 (14,15). The extent to which urban air pollution contributes to this excess is uncertain, but studies on diesel-exposed occupational groups support a causative role by traffic-related pollution (14,16).

Lung cancer is the most common cause of cancer death among males in Sweden. Stockholm County has a higher incidence of lung cancer than Sweden as a whole, and in the central parts of Stockholm the incidence increases even further (17). A large case-control study of lung cancer in Stockholm males--the LUCAS (Lung Cancer in Stockholm) study--aimed to determine the role of known and potential risk factors, such as smoking, occupational exposures, and exposure to ambient air pollution. Our methodologic study, based on the LUCAS study, specifically aimed to develop methods to assess individual historical exposure to ambient pollution using reconstructed emission data, dispersion models, and GIS techniques. Nitrogen oxides ([NO.sub.x] and [NO.sub.2]) and sulfur dioxide were chosen as indicators of air pollution from motor vehicles and residential heating, respectively.

Materials and Methods

The source population consists of all men 40 to 75 years of age who lived in Stockholm County at any time between 1985 and 1990, and who had lived in the County since 1950, with a maximum of five years of residence outside the County. From this population we recruited three groups, which form the study population for this study: 1,042 cancer cases reported to the Stockholm County cancer registry 1985-1990, and two population-based control groups comprising 2,364 individuals, age-stratified to the case group (18,19).

We collected information on individual exposures using a postal questionnaire sent to living respondents or to next of kin when the study subject had died. The questionnaire included questions on occupations, smoking habits, and food consumption. To assess exposure to air pollution, we also collected information on all dwellings where the case or control had lived for more than one year from 1955 onward. In case of incomplete data in form of gaps, overlaps, or other inconsistencies in the residential histories, we sought supplemental information from the local offices for demography (parishes) and/or the tax authority. Respondents who had lived in apartment houses were also asked whether they had lived on the ground floor.

We transformed the addresses into geographic coordinates using the MAPINFO computer software (20) in conjunction with a regional geographic address database (21). The address database covers all urban areas in Stockholm County, defined as areas with more than 200 inhabitants and no more than 200 meters between houses (22). The database contains the names of all streets in the urban area. Each street is divided into segments (generally the length of a block), with geographic coordinates on both ends of the street segment and with street numbers on each side of the street at both endpoints. When an address with street and street number is matched to the address database, coordinates for street numbers between endpoints of the street segment are calculated through interpolation.

We geocoded the 10,800 addresses in the study in a three-step procedure. The first step used the automatic geocoding function in MAPINFO. Addresses that could not be coded directly by the program were coded interactively in a second step, with an operator adjusting spelling or street number. In the third step, addresses not matched interactively were coded manually.

We tested the reliability of the automatic coding function using two external suppliers of geocoding services [Telemedia (Stockholm, Sweden), a private company making the maps for the Swedish telephone directory, and CFD (Centrala fastighetsdata; Gavle, Sweden), central authority for real estate data]. Telemedia uses a technique similar to that of MAPINFO, whereas CFD sets the coordinates for an address to the center of the property.

We assessed exposure to traffic-related air pollution ([NO.sub.x] and [NO.sub.2]) for 1955-1990 for the whole territory of Stockholm County, approximately 6,500 [km.sup.3]. The exposure assessment was based on geographic modeling of three points in time: 1960, 1970, and 1980. We used the AIRVIRO dispersion modeling system in conjunction with retrospectively constructed emission databases for these years. We used a detailed regional emission database for 1993, collected by the Stockholm-Uppsala Air Quality Management Association (Stockholm), as a basis (23). It contains approximately 4,300 line sources related to road traffic, covering all roads with more than 1,000 vehicles/24 hr, which corresponds to 90% of the estimated total emission from road traffic. It also contains information on over 500 point sources, including major industries and energy plants as well as small industry and ferries in ports. Limited diffuse emission sources, e.g., air traffic and merchant vessels in commercial routes, are treated as area sources, and several population-density related sources, such as local heating, are mapped as grid-sources (250 m or 1,000 m grids). Work machine emissions are treated as grid sources, distributed with one-quarter over industrial areas, and the remainder according to population density. The total emission of [NO.sub.x] in the region in 1993 is estimated at 18,900 tons, of which 79% originated from road traffic.

To reconstruct comparable historical emission databases, we collected retrospectively data from 1955 to 1993 on the growth of the urban areas in the Stockholm county region, the development of the district heating system, and the growth and redistribution of the road traffic (Table 1). We constructed separate emission databases for regional emissions of [NO.sub.x] for 1960, 1970, and 1980. The two former databases contain only contributions from road traffic, whereas the 1980 database was constructed in two versions: with road traffic only and with all sources. The total emissions of [NO.sub.x] from road traffic in these databases are 13,400, 18,100, and 20,800 tons/year.

Table 1. Examples of input data for reconstruction of [NO.sub.x]
emission data in Stockholm County for 1960, 1970, and 1980, and
local interpolation of [NO.sub.2] levels 1955-1959, 1961-1969,
1971-1979, and 1981-1990.

Item, year                       Data/event

Traffic counts, Stockholm city
  and suburbs
  1955-1970                      Linear increase from 210,000 to
                                   490,000 vehicles/day passing city
                                   limits (+19,000/year)
  1970-1980                      Small changes
  1981-1989                      Linear increase from 470,000 to
                                   620,000 vehicles/day passing
                                   inner-city limits (+19,000/year)

Traffic counts, rest of county
  1955-1970                      Estimated at a linear increase
  1971-1980                      Small changes
  1981-1989                      Linear increase from 530,000 to
                                   720,000 vehicles/day passing
                                   outer city limits (+27,000/year)

Total traffic work,              608, 901, 882, and 1,008 vehicle
  Stockholm city                   Gm/year(a)
  1960, 1970, 1980, and 1993

Total traffic work, city         1,868, 2,832, 2,876, and 3,370
  and suburbs                      vehicle Gm/year
  1960, 1970, 1980, and 1993

Total traffic work,              2,678, 3,305, 4,016, and 5,103
  rest of county                   vehicle Gm/year
  1960, 1970, 1980, and 1993

Major through roads              Two parts of the western link open
  1966 and 1970

Public transport                 Unified taxes in the whole region
  1971

(a) 1 vehicle Gm/year = [10.sup.6] vehicle km/year.

Emissions from heating were assessed as [SO.sub.2] in a similar way as those from [NO.sub.x] were from road traffic. We reconstructed three historical regional emission databases for the same years as for the [NO.sub.x] emissions, with information again based on the 1993 database. Detailed data on the development of district heating (point sources) and other energy plants, and local restrictions of sulfur content in oil were used for the reconstruction (Table 2), as well as [SO.sub.2] measurement data. We used data from five monitoring stations (annual and decade mean levels) to calibrate the model for average emission levels from grid-type sources (mainly local oil-fueled house heating) for 1970 and 1980, through iterative model runs. Measured levels from the two central monitoring stations were allowed to influence the relative level of emissions from grid sources in the 1960 and 1970 databases, respectively.

Table 2. Examples of input data for reconstruction of [SO.sub.2]
emission data for 1960, 1970, and 1980, and local interpolation of
[SO.sub.2] levels 1955-1959, 1961-1969, 1971-1979, and 1981-1990.

Item, year                 Data/event

House heating fuel
  1955-1960                Gradual shift from carbon to oil

House heating methods
  1955-1969                Mainly independent for each building
                             or property
  1970-1979                Gradual shift to district heating
  1980                     50% of heating energy in Stockholm city
                             from district heating
  1990                     70% of heating energy in Stockholm city
                             from district heating, remaining areas
                             converted to electric heating

Energy conservation
  1973-1979                Oil crisis started a series of actions on
                             energy conservation in house heating

Sulfur limit in oil,
  Stockholm city
  Until 1968               Not regulated
  1969-1990                1% in the city center
  October 1969             1% in the rest of the city

Sulfur limit, Stockholm
  county except the city
  Until 1969               Not regulated
  1969-1970                3.5%
  1971-1983                1% gradually implemented
  1984-1990                1%

Estimated average sulfur
  content in fuel
  1955-1968                3% in all the county
  1969-1970                1% in the city; 3% in rest of the county
  1971-1990                1% in all the county

We performed dispersion calculations from the retrospectively constructed emission databases using the Gaussian model in the AIRVIRO system (8). We used the 1990s average distribution of climate in 180 different conditions over a year for all periods, and we used no correction factors. We calculated [SO.sub.2] and [NO.sub.x] annual mean concentrations for the relevant period, using only the contributions from road traffic and heating sources, respectively. The [NO.sub.x] concentrations were transformed to [NO.sub.2] data using the equation

[[NO.sub.2]] = [[NO.sub.x]]**{A+B/(C+[[NO.sub.x]])},

in which the parameters A, B, and C were estimated from measurements in Stockholm County in the early 1980s (hourly means from 10 fixed and five mobile stations). The parameters A and B were allowed to vary according to location, and the estimates were in the ranges 0.60-0.70 and 30-40, respectively. The parameter C was estimated at 100 throughout the county. The performance of these conversions was tested in 220,418 paired 1-hr measurements of [NO.sub.x] and [NO.sub.2] from the 1990s. The correlation between [NO.sub.2] levels estimated from [NO.sub.x] and measured [NO.sub.2] levels was 0.91.

We performed the dispersion calculations in four different resolutions: 2,000 m x 2,000 m, 500 m x 500 m, 200 m x 200 m, and 100 m x 100 m grids. We used the highest resolutions in central city areas and the lowest resolutions in the countryside. For each grid we estimated the average levels of the index pollutants. All illustrations are based on 500 m x 500 m grids. The isolines in the illustrations were not used for the calculations.

For main streets in the city center, we added a street canyon contribution to the assessed roof concentrations of [NO.sub.x]. We assessed this contribution by dispersion calculations with the AIRVIRO street canyon model, and summarized it in an 80% addition at the ground floor and a 40% addition for all other dwellings in the streets concerned. We gave 68 addresses (2% of 3,406) an addition to the 1960 dispersion estimate; 58 addresses an addition to the 1970 dispersion estimate, and 44 to the 1980 dispersion estimate. Additions to the assessed levels of [NO.sub.2] were made in the same way, but with a 50% addition at the ground floor and a 20% addition for other dwellings in these streets.

For 1980, we compared the dispersion model calculations including all sources of [NO.sub.x]--after transformation to [NO.sub.2] as indicated above--with measurements of [NO.sub.2] in five rooftop or background sites (one stationary and four mobile stations) taken in the early 1980s (Figure 1). In the first model run, all five annual values predicted by the model at the highest resolution (not shown in Figure 1) were within [+ or -] 20% of the value measured at their respective locations. The model was successively calibrated to minimize this deviation. All available measurements for [SO.sub.2] had been used to calibrate the model, so no formal validation of the [SO.sub.2] modeling was possible.

[ILLUSTRATION OMITTED]

To create annual levels of exposure to [SO.sub.2] and [NO.sub.2] from the modeled concentration fields, we computed extra- and intrapolated values for each address for each year between 1955 and 1990. We modeled the [NO.sub.2] concentration in a geographic point to increase linearly from 1955 to 1970, using the estimated values for 1960 and 1970 to obtain the slope. The rationale for linearity is that the traffic counts increased approximately linearly in this period. Between 1970 and 1980 the changes were small and were also modeled as linear. We modeled the increases for after 1980 as linear with the same slope as 1955-1970, motivated by a similar linear increase in traffic counts. In the vicinity of new major roads and the corresponding relieved roads, the slopes before and after the change were extracted from surrounding areas.

We modeled [SO.sub.2] concentrations as constant during the period 1955-1967 (1960 estimate), after which there was a sharp linear decrease until 1969, and we used the 1970 estimate for the 1969-1979 period. Another rapid linear decrease occurred for 1979-1980, and we obtained the estimates for 1981-1990 from a linear decrease to low levels in 1995. Measurements show that from that year the levels are about 5 [micro]g/[m.sup.3] in central Stockholm, 7 [micro]g/[m.sup.3] in the suburbs, and 2.5 [micro]g/[m.sup.3] in the rest of the county.

We then linked the annual air pollution data to 10,800 individual address coordinates for the relevant time intervals, yielding exposure indices for each of the three air pollution indicators for each year from 1955 to 1990. For 565 individuals (17%) there were gaps in the information on residency, on average for 2.3 years and never exceeding 5 years per individual. Individual average exposure was calculated for the three decades preceding the year in which follow-up ended (year of lung cancer diagnosis for cases and year of selection for controls).

Besides the individual exposure estimates, we calculated the average annual level of all addresses in the study, which may be interpreted as an annual area average, approximately weighted for average population.

The results from epidemiologic analyses based on these exposure estimates are presented elsewhere (19).

Results

The dispersion calculations for [NO.sub.2] from road traffic show a substantial increase of the average level with time (Figures 2 and 3), mainly from 1960 to 1970. The increase affects most parts of greater Stockholm and is caused mainly by the increase in traffic volume (Table 1). The dispersion calculations for [SO.sub.2] from heating show a substantial decrease of the levels with time, from above 130 [micro]g/[m.sup.3] in central Stockholm in the end of the 1960s to less than 50 [micro]g/[m.sup.3] in the whole area at the beginning of the 1980s (Figures 4 and 5). The main reasons for this decrease are restrictions in the sulfur contents of heating oil, massive investment in district heating, and investment in emission control at the energy plants (Table 2). The decrease is more pronounced in densely populated areas because of earlier implementation of more strict sulfur standards and a more rapid development of district heating.

[ILLUSTRATIONS OMITTED]

Of the 10,800 addresses in the study, 56% were geocoded by the automatic geocoding function in MAPINFO. A further 27% were geocoded by the computer program with minor operator assistance. The remaining 17% of the addresses had to be geocoded manually; these addresses were either outside urban areas or too inexact to be geocoded interactively. Three percent of all the addresses could not be found and were given approximate coordinates, e.g., of a village center. We checked the result of the geocoding visually by examining maps where the geocoded addressed were plotted.

For a subset (n = 100) of the addresses, we performed a reliability test of the automatic coding of geographic coordinates using two external suppliers of geocoding services. In this subset, 68 addresses could be geocoded automatically with all three procedures. Both external codings agreed reasonably well with our automatic coding. The average deviance from our coding was about 50 m in both; 80% of the addresses agreed within 100 m for all three codings; and our coding agreed within 100 m for almost 90% in comparison with any one of the other two codings. The maximal deviance was 168 m and 237 m, respectively, when our coordinates were compared with those from the external suppliers.

The estimated annual mean concentration levels for [NO.sub.2] from road traffic and [SO.sub.2] from heating, at all of the 10,800 addresses in the study (i.e., considering the whole period for every address), shows clearly the increase of [NO.sub.2] and decline of [SO.sub.2] in Stockholm County (Figures 6 and 7). The corresponding annual means for the study population individuals are similar to those of all the addresses. There is a slight tendency of temporal decrease in the individual means compared with the average address means, caused by a net migration from more polluted to less polluted areas (Figures 6 and 7).

[GRAPHS OMITTED]

All individual exposure indicators show approximate normal distribution (Table 3). Despite the increase in [NO.sub.2] levels during the period, the average estimated individual traffic-related [NO.sub.2] levels during the three decades before the end of follow-up are well below the current World Health Organization (WHO) ambient concentration guideline of 40 [micro]g/[m.sup.3] for [NO.sub.2], which refers to contributions from all sources. This guideline defines an annual limit value for the protection of human health, and is numerically identical to that of the European Union recommendation proposed for 2010 (24). The maximum of the individual period average is about 11 times higher than the minimum. The total range extends just over the WHO limit.

Despite the decline of emissions since 1960, the average estimated individual heating-related [SO.sub.2] level during the 30-year study period exceeds the current WHO recommended guideline of 50 [micro]g/[m.sup.3] for the annual mean ambient [SO.sub.2] concentration from all sources (25). About 58% of the individuals have study-period mean exposures from heating-related sources that exceed this level. The total range of individual average levels is quite large, with a maximum that is about 20 times higher than the minimum. That the levels decrease with time is reflected in differences in individual average levels in the three decades preceding individual end of follow-up. The mean and most other parameters decrease by about two-thirds from the first to the last decade (Table 3). The interquartile ratio is about 1.7 for all exposure metrics, and the ratio of the means of the upper and lower quartiles is about 3 (Table 3).

Table 3. Descriptive statistics for average estimated individual
exposure to [NO.sub.2] from road traffic and [SO.sub.2] from house
heating, used in a case-control study of lung cancer in Stockholm.
Individual exposure is reported for three decades before the year
in which follow-up ended. Lung cancer incidence was followed up in
1985-1990 and the individual decades fall within the calendar years
1955-1969, 1965-1979, and 1975-1989.

                          Valid n   Mean   Median   Minimum   Maximum

Average exposure to
    road traffic fumes
    ([micro]g/[m.sup.3]
    [NO.sub.2])
  Decade 1-3(a)            3,406     20      20        4         44
  Decade 1                 3,406     18      17        4         51
  Decade 2                 3,406     21      20        4         56
  Decade 3                 3,406     22      21        4         50
Average exposure to
    house heating fumes
    ([micro]g/[m.sup.3]
    [SO.sub.2])
  Decade 1-3(a)            3,406     53      53        6        110
  Decade 1                 3,406     87      88        8        185
  Decade 2                 3,406     46      44        6        108
  Decade 3                 3,406     27      25        4         65

                             Lower      Upper     Ratio of
                            quartile   quartile   quartile
                            cutpoint   cutpoint   cutpoints

Average exposure to
    road traffic fumes
    ([micro]g/[m.sup.3]
    [NO.sub.2])
  Decade 1-3(a)                15         25         1.6
  Decade 1                     13         23         1.8
  Decade 2                     15         26         1.7
  Decade 3                     16         28         1.8
Average exposure to
    house heating fumes
    ([micro]g/[m.sup.3]
    [SO.sub.2])
  Decade 1-3(a)                41         67         1.6
  Decade 1                     66        110         1.7
  Decade 2                     33         59         1.8
  Decade 3                     19         34         1.8

                             Lower      Upper     Ratio of
                            quartile   quartile   quartile
                            average    average    averages

Average exposure to
    road traffic fumes
    ([micro]g/[m.sup.3]
    [NO.sub.2])
  Decade 1-3(a)                11         30         2.8
  Decade 1                      9         29         3.1
  Decade 2                     11         32         3.0
  Decade 3                     11         34         3.1
Average exposure to
    house heating fumes
    ([micro]g/[m.sup.3]
    [SO.sub.2])
  Decade 1-3(a)                29         78         2.7
  Decade 1                     45        129         2.8
  Decade 2                     23         71         3.1
  Decade 3                     13         43         3.2

(a) Decades before the individual year of the end of follow-up.
Example: For an individual with end of follow-up in 1988,
Decades 1, 2, and 3 correspond to the calendar years 1958-1967,
1968-1977, and 1978-1987.

The individual 30-year average exposures to [NO.sub.2] and [SO.sub.2] are correlated (r = 0.65). Despite this correlation, there are obvious differences in individual exposure to the two agents. For example, at an individual average [SO.sub.2] exposure of 80 [micro]g/[m.sup.3], there is a 3-fold total variation (15-40 [micro]g/[m.sup.3]) of the individual average exposure to [NO.sub.2].

Discussion

Recently, atmospheric dispersion models have been applied in a retrospective study of health effects of a point source of arsenic dust (26). Some other epidemiologic studies have used regression models of measured air pollution levels in GIS-based exposure assessment--for example, the Small Area Variation in Air Pollution and Health (SAVIAH) project (27). To our knowledge, our study is the first attempt at retrospective application of atmospheric dispersion models to assess exposure to a wide range of emission sources that are followed over several decades. We assessed the exposure at the individual level, linked to residence over a 30-year period, and based our assessment on indices of complex mixtures of air pollution from road traffic and house heating in an entire region.

The geocoding of addresses is often central to the use of a GIS in environmental epidemiology research. The success of the geocoding depends on the quality of the databases available. It became clear in the present study that the regional address database contained errors such as spelling mistakes and local differences in the completeness of data. This is likely to occur in most detailed databases of this kind. Also, the interpolation technique, used in the present study, yields a set of address coordinates close to the road, which should describe road traffic exposure more accurately than the method using the property centroid. A more formal validation of the address coordinates (e.g., by using the satellite Global Positioning System technique) was not performed, but we compared our method to geocoding using two other address databases, for a subset of the addresses. There was a difference of less than 100 m for 90% of the addresses. Since the highest resolution of the air dispersion model was 100 m x 100 m, the geocoding accuracy seems satisfactory.

We assessed both [NO.sub.x] and [NO.sub.2] as indices of pollution from road traffic. [NO.sub.x] represents the direct emissions and is thus, in principle, a better index of traffic exhaust exposure. [NO.sub.2] is also important because of its toxicity and widespread use in exposure assessment. The relation between the two indices was derived from observations in Stockholm in the 1980s. The assessed relation is relevant for the period before catalytic converters were implemented, and showed good agreement with the measured relation in the beginning of the 1990s. However, it may not be used outside our study period and location, because it is dependent on composition of automobile fleet and climatologic conditions.

Gaussian dispersion models running on personal computers can now be implemented in several places (7). The main challenge is not the model itself, but rather the quality of the emission and meteorologic data. We used the same average yearly distribution of meteorologic conditions for all study periods, which should be sufficient for calculation of long-term average levels. The alternative is using time-series weather data, which would allow spatial assessment of the pollutant level day to day. However, this places enormous demands on computational capacity and was judged unnecessary. The ultimate answer on how well a model run performs is gained by comparison with independently measured pollutant levels. Unfortunately, we could perform this type of validation only for [NO.sub.2].

The 1993 emission database for Stockholm County contains very detailed descriptions of road traffic and house heating. We constructed historic databases for 1980, 1970, and 1960 mainly by successive additions and deletions to the 1993 database. Modeled ambient [NO.sub.2] data based on the 1980 database compared well with actual measurements for that period. The current database has been somewhat improved, and we used the 1995 edition to compare measured [NO.sub.2] data in 16 locations throughout Stockholm for 1994-1997. Modeled annual average levels correlated very well (r = 0.96) with the measurements [site averages in the range 4-26 [micro]g/[m.sup.3] (28)].

Assumptions regarding emission factors, proportion of heavy traffic, and different heating fuel qualities are likely to be less valid for the beginning than for the end of the study period. The average pollutant levels assessed by modeling based on the early databases may therefore be subject to greater measurement error. For traffic emissions, however, the main source of geographic exposure contrast at any given time is the difference in traffic work between different neighborhoods, which is well documented in traffic flow charts from the 1960s onward. The interindividual [NO.sub.x] or [NO.sub.2] exposure contrasts within limited time periods, e.g., decades, are therefore expected to be reasonably valid in the beginning of the study period.

In contrast, the differences in [SO.sub.2] levels depend on the distribution of different heating methods. The estimated emissions of [SO.sub.2] from heating in 1970 were lower per capita than those estimated for the whole country based on fuel consumption in the beginning of the 1970s (29,30). Despite this discrepancy, we judge the temporal development of average house heating to be quite accurate. However, the geographic variations are less well documented in the reconstructed emission databases. The interindividual [SO.sub.2] exposure contrasts in the early study period may therefore be less valid than those for traffic exhausts.

Despite the mobility of study subjects (on average about three residences in 30 years), the results of the present study show ranges of assessed individual exposure to be quite wide. Considering the average individual exposure for three decades before the end of follow-up, the ratio between maximum and minimum was approximately 20 for [SO.sub.2] and 10 for [NO.sub.2]. When the subjects are divided in quartiles, the ratio between the average exposure in the upper group and that in the lower group is almost 3 for both indices. Within the different decades, this contrast is somewhat higher. It thus seems possible to categorize this population in different levels of exposure, which permits the detection of exposure response relations. Furthermore, despite the overall correlation between the two agents, there were substantial differences in the relation between exposures to [NO.sub.2] and [SO.sub.2] for study individuals, which should make it possible to differ between these two pollutants in epidemiologic analyses.

However, the outdoor levels assessed in residential areas to a varying extent reflect true individual total exposure. Exposure contributions at work and in traveling are not included. In addition, perfectly measured individual levels of these pollutants have proxy properties, because they are single components in very complex mixtures of pollutants from road traffic and heating. The relation of [NO.sub.2] and [SO.sub.2] with the other components of these mixtures is likely to have changed over time, and may also vary across the territory. Any positive finding of exposure-response based on this exposure information is therefore likely to be an underestimate of the total effects of air pollution from these sources.

Conclusion

To our knowledge, this is the first attempt to use reconstructed emission data together with dispersion modeling and GIS for retrospective individual exposure assessment of multi-source air pollution in an epidemiologic study. The results show that the technique may be useful for exposure assessment in environmental epidemiologic studies, given that detailed emission data or environmental exposure data are available for study-relevant time periods and with good spatial resolution. Another prerequisite is that the biologically relevant true individual exposure contrast must be high enough for individual contrast to remain in a geographically based assessment of proxy exposures.

REFERENCES AND NOTES

(1.) Croner CM, Sperling J, Broome FR. Geographic information systems (GIS): new perspectives in understanding human health and environmental relationships. Stat Med 15:1961-1977 (1996).

(2.) Vine MF, Degnan D, Hanchette C. Geographic information systems: their use in environmental epidemiologic research. Environ Health Perspect 105:598-605 (1997).

(3.) Briggs DJ. Mapping environmental exposure. In: Geographical and Environmental Epidemiology: Methods for Small-area Studies (Elliott P, Cuzick J, English D, Stern R, eds). Oxford:Oxford Univesity Press, 1992;158-176.

(4.) Briggs DJ, Elliott P. The use of geographical information systems in studies on environment and health. World Health Stat Q 48:85-94 (1995).

(5.) Stockwell JR, Sorensen JW, Eckert JW Jr, Carreras EM. The U.S. EPA Geographic Information System for mapping environmental releases of Toxic Chemical Release Inventory (TRI) chemicals. Risk Anal 13(2):155-164 (1993).

(6.) Valjus J, Hongisto M, Verkasalo P, Jarvinen P, Heikkila K, Koskenvuo M. Residential exposure to magnetic fields generated by 110-400 kV power lines in Finland. Bioelectromagnetics 16(6):365-376 (1995).

(7.) Moussiopoulos N, Berge E, Bohler T, de Leeuw F, Gronskei K-E, Mylona S, Tombrou M. Ambient Air Quality, Pollutant Dispersion and Transport Models. Topic Report 19/96. Copenhagen, Denmark:European Environment Agency, 1996.

(8.) SMHI: Indic Airviro. Swedish Meteorological and Hydrological Institute. Technical Description of the Dispersion Models. Norrkoping, Sweden: SMHI: Indic. Swedish Meteorological and Hydrological Institute, 1993.

(9.) Swedish Meteorological and Hydrological Institute. Airviro User Documentation. Norrkoping, Sweden: Swedish Meteorological and Hydrological Institute, 1997.

(10.) SMHI: Indic Airviro. Swedish Meteorological and Hydrological Institute. Available: http://www.indic-airviro.smhi.se [cited 22 December 1999].

(11.) Higgins IT. Epidemiological evidence on the carcinogenic risk of air pollution [Review]. IARC Sci Publ 13:41-52 (1976).

(12.) Doll R. Atmospheric pollution and lung cancer. Environ Health Perspect 22:23-31 (1978).

(13.) Friberg L, Cederlof R. Late effects of air pollution with special reference to lung cancer. Environ Health Perspect 22:45-66 (1978).

(14.) Hemminki K, Pershagen G. Cancer risk of air pollution: epidemiological evidence. Environ Health Perspect 102(suppl 4):187-192 (1994).

(15.) Katsouyanni K, Pershagen G. Ambient air pollution exposure and cancer. Cancer Causes Control 8(3):284-291 (1997).

(16.) Cohen AJ, Pope CA III. Lung cancer and air pollution. Environ Health Perspect 103(suppl 8):219-224 (1995).

(17.) Nyberg F. Environmental and Genetic Factors in Lung Cancer: Epidemiological and Biomolecular Studies Focusing on Nonsmokers [PhD Thesis]. Stockholm, Sweden:Karolinska Institutet, 1998.

(18.) Gustavsson P, Jakobsson R, Nyberg F, Pershagen G, Jarup L, Scheele P. Occupational exposure and lung cancer risk: a population-based case-referent study in Sweden. Am J Epidemiol 152:32-40 (2000).

(19.) Nyberg F, Gustavsson P, Jarup L, Bellander T, Berglind N, Jakobsson R, Pershagen G. Urban air pollution and lung cancer: a population-based case-control study from Stockholm. Epidemiology 11(5):487-495 (1999).

(20.) MapInfo Professional (TM) User's Guide, Version 4.1. New York: MapInfo Corporation, 1995

(21.) Tatort 2000 [in Swedish]. Stockholm, Sweden: Kartcentrum, Lantmateriverket (National Land Survey of Sweden), 1996.

(22.) Statistics Sweden. Description of Previous Census of the Population and Housing [in Swedish]. Available: http://www.scb.se/statinfo/pbesk/be0205.asp [cited 22 December 1999].

(23.) Ekstrom A, Hadenius A, Jonson T. Emission database 93: a documentation [in Swedish]. Stockholms Ians luftvardsforbund 2:95. Stockholm luft- och bulleranalys. Stockholm:Stockholms Miljoforvaltning, 1995.

(24.) Commission Proposes New Air Quality Limit Values Available:http://europa.eu.int/comm/environment/press/ ip97853.htm [cited 22 December 1999].

(25.) WHO. Air Quality Guidelines. World Health Organization, Geneva. Available: http://www.who.int/peh/air/ Airqualitygd.htm [cited 22 December 1999].

(26.) Ihrig MM, Shalat SL, Baynes C. A hospital-based case-control study of stillbirths and environmental exposure to arsenic using an atmospheric dispersion model linked to a geographical information system. Epidemiology 9(3):290-294 (1998).

(27.) Briggs DJ, Collins S, Elliott P, Fischer P, Kingham S, Lebret E, Pryl K, van Reeuwijk H, Smallbone KL, van der Veen. Mapping urban air pollution using GIS: a regression-based approach. Int J Geogr Inf Sci 11:699-718 (1997).

(28.) Johansson C, Hadenius A, Johansson P-A, Jonson T. The Stockholm Study on Health Effects of Air Pollution and their Economic Consequences. I: [NO.sub.2] and Particulate Matter in Stockholm, Concentrations and Population Exposure. Stockholm:Stockholm Environment and Health Protection Administration, Air quality and Noise Analysis, 1999.

(29.) Rodesjo B. Sulphur Dioxide Emission Inventory for Sweden 1972 [In Swedish: Emissionsinventering av svaveldioxid: Sverige 1972]. Solna, Sweden:Swedish Environmental Protection Agency, 1973.

(30.) Bostrom C-E, Levander T, Persson B. Air Pollution in Sweden 1970-1980. [In Swedish: Luftfororeningar i Sverige 1970-1980]. (SNV PM 1521). Solna, Sweden:Swedish Environmental Protection Agency, 1982.

Tom Bellander,(1,2) Niklas Berglind,(1) Per Gustavsson,(3) Tage Jonson,(4) Fredrik Nyberg,(2,5) Goran Pershagen,(1,2) and Lars Jarup(1,2,6)

(1) Department of Environmental Health, Stockholm County Council, Stockholm, Sweden; (2) Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden; (3) Department of Occupational Health, Stockholm County Council, Stockholm, Sweden; (4) Environment Health Protection Administration, Slb-analys, Stockholm, Sweden; (5) International Agency for Research on Cancer, Lyon, France; (6) Department of Epidemiology and Public Health, Imperial College, London, UK

Address correspondence to T. Bellander, Department of Environmental Health, Norrbacka 3rd floor, Karolinska Hospital, S-176 76 Stockholm, Sweden. Telephone: [phone]. Fax [phone]. E-mail: [e-mail]

The work was performed at Department of Environmental Health, Karolinska Hospital, Norrbacka, S-171 76 Stockholm, Sweden.

The study was supported by the Swedish Environmental Protection Agency, the Swedish Council for Work Life Research, and the Stockholm County Council.

We are indebted to C. Bengtsson, A. Boberg, E-B. Gustavsson, M. Hugosson, L-G. Horte, R. Jakobsson, E. Junghahn, U. Klinga, A. Lauber, S. Norberg, M. Pettersson, C. Rudengren, P. Scheele, E. Skarke, A. Wahlbom, and K. Astrom for valuable contributions to various aspects of the study.

Received 8 February 2000; accepted 11 December 2000.