ResearchOpen Access

Traffic-Related Air Pollution and Asthma Onset in Children: A Prospective Cohort Study with Individual Exposure Measurement

    Published:https://doi.org/10.1289/ehp.10968Cited by:236

    Abstract

    Background

    The question of whether air pollution contributes to asthma onset remains unresolved.

    Objectives

    In this study, we assessed the association between asthma onset in children and traffic-related air pollution.

    Methods

    We selected a sample of 217 children from participants in the Southern California Children’s Health Study, a prospective cohort designed to investigate associations between air pollution and respiratory health in children 10–18 years of age. Individual covariates and new asthma incidence (30 cases) were reported annually through questionnaires during 8 years of follow-up. Children had nitrogen dioxide monitors placed outside their home for 2 weeks in the summer and 2 weeks in the fall–winter season as a marker of traffic-related air pollution. We used multilevel Cox models to test the associations between asthma and air pollution.

    Results

    In models controlling for confounders, incident asthma was positively associated with traffic pollution, with a hazard ratio (HR) of 1.29 [95% confidence interval (CI), 1.07–1.56] across the average within-community interquartile range of 6.2 ppb in annual residential NO2. Using the total interquartile range for all measurements of 28.9 ppb increased the HR to 3.25 (95% CI, 1.35–7.85).

    Conclusions

    In this cohort, markers of traffic-related air pollution were associated with the onset of asthma. The risks observed suggest that air pollution exposure contributes to new-onset asthma.

    Toxicologic and epidemiologic research suggests that air pollution exacerbates asthma symptoms (Delfino 2002; Eder et al. 2006), but few prospective studies have addressed the question of whether this pervasive exposure contributes to disease onset in children. Recent research has focused on the contribution of traffic-related air pollution, partly due to the toxicologic effects of the pollution mixture from mobile sources (Nel 2005). Studies based on prevalent asthma have been inconsistent, with some reporting positive associations between traffic-related exposure and asthma or associated symptoms (Brunekreef et al. 1997; Ciccone et al. 1998; Duhme et al. 1996; Edwards et al. 1994; Gauderman et al. 2005; McConnell et al. 2006; Migliaretti and Cavallo 2004; Nitta et al. 1993; Oosterlee et al. 1996; van Vliet et al. 1997; Weiland et al. 1994; Wjst et al. 1993) and others reporting no significant associations (Lewis et al. 2005; Livingstone et al. 1996; Wilkinson et al. 1999).

    A few prospective studies have assessed the relation between traffic pollution and asthma onset (Brunekreef and Sunyer 2003; Eder et al. 2006; Gold and Wright 2005). A case–control study from France that selected subjects by medical clinic found that early-life exposure to traffic pollutants may contribute to asthma incidence (Zmirou et al. 2004). In a prospective birth cohort study from the Netherlands, traffic-related pollutants were associated with incident wheezing and doctor-diagnosed asthma in young children up to 4 years of age (Brauer et al. 2002, 2007). An earlier study from Japan showed positive associations between community nitrogen dioxide levels measured at a central monitor and asthma incidence in school children (Shima et al. 2002). None of these studies used individually based measures of exposure, and modeled exposures may have large errors that attenuate effects (Molitor et al. 2006). In this context, we investigated the relationship between childhood asthma onset and measured markers of traffic-related air pollution in Southern California.

    Materials and Methods

    Sample selection

    We selected a sample of 217 children from 917 eligible subjects living in 11 communities in the Southern California Children’s Health Study (CHS), a prospective cohort study of air pollution and respiratory health. Subjects were selected randomly from the larger cohort based on different levels of traffic-related exposure. Specifically, we divided subjects into two strata (above and below the median traffic exposures within each community) and then selected an equal number randomly from each of the strata. Although our study is similar to an earlier cross-sectional study from this population that focused on the effect of freeway-related traffic in 10 CHS communities (Gauderman et al. 2005), in this incident asthma analysis we excluded children with a lifetime history of asthma at study entry, and we included an 11th community, Lompoc, which had local traffic influence but no major freeways. A 12th community, Lake Arrowhead, had little local traffic variability and was not included in either analysis.

    Case ascertainment was as described previously (Gauderman et al. 2005; McConnell et al. 2002a). We excluded children with prevalent asthma from this analysis if a parent answered yes to the question “Has a doctor ever diagnosed this child as having asthma?” on a baseline survey sent home with children at the time of study enrollment (at 10 years of age) in 1993 or 1996. We also excluded children who answered yes to the question “Has a doctor ever said you had asthma?” on a questionnaire administered during the first year of the study by a trained interviewer. We then classified as incident cases children who answered yes to this question on any annual interview during up to 8 years of follow-up until high school graduation. Information on demographic and other known risk factors for asthma development, such as pets, molds in the home, and environmental tobacco smoke, was obtained at baseline and updated annually. We obtained informed consent from parents, and the study was approved by the University of Southern California Institutional Review Board.

    Exposure assignment

    As described previously (Gauderman et al. 2005), children had NO2 monitors (Palmes tubes; Gradko, Winchester, UK) placed in the front or back yard of their home in the year 2000 for 2 weeks in the summer (mid-August) and 2 weeks in the fall–winter season (mid-November). We deployed all monitors within 4 days to minimize temporal variation of the exposure estimates by change of season or meteorologic events. Within each community, monitors were deployed within a 4-hr period. We obtained valid readings for 209 subjects for the ambient household level of NO2 measurements during fall–winter months and 204 for summer months. We assessed exposures by season and as an average of the two measurement periods.

    Contextual confounding variables

    We also derived a series of potential “contextual” confounders at the community level. Contextual characteristics of the community, particularly socioeconomic position (SEP), have been associated with respiratory health effects (Haan et al. 1987; Jin et al. 1994; Juhn et al. 2005; Wissow et al. 1988). For this study, we tested the SEP measures of median household income, proportion of respondents with low education (i.e., no high school diploma), percent of males unemployed (as a marker for full-time income instability), and percent living in poverty as confounders. We used male unemployment (as opposed to general or female unemployment) because the relationship between employment and health is stronger for males than for females (Jin et al. 1994; Sorensen and Verbrugge 1987). Sensitivity analyses were also conducted with total unemployment. Definitions of unemployment and poverty followed those developed by the U.S. Census Bureau (1992). We used these variables to determine whether variation in asthma could be explained by the social conditions in the neighborhood because such contextual effects have been suggested as another risk factor for asthma incidence (Gold and Wright 2005).

    Characteristics of SEP were aggregated from census blocks where study subjects lived. We used data from the U.S. 1990 census to estimate community conditions at the time of study enrollment, which began in 1993. Because we recruited subjects through schools, residential locations clustered together. To compensate for clustering of study subjects, we weighted these community-level data by the proportion of the census blocks included in a community-specific bounding rectangle that contained 95% of local study subjects (Tatalovich et al. 2006). This spatial bounding captures the social characteristics around the residence of the children but excludes a few spatial outliers.

    Meteorologic variables

    Meteorologic conditions, especially temperature and humidity, may also be potential risk factors for asthma (Guo et al. 1999; Weiland et al. 2004). Based on data collected daily in 1995 (close to the midpoint of baseline enrollment) from local meteorologic stations in each community, we derived 10 markers of humidity and temperature, including minimum and average temperature in winter and summer months, maximum and average humidity in winter and summer months, and annual average temperature and humidity. We tested these as possible confounders of the asthma–air pollution relationship.

    Statistical model

    We used random-effects Cox proportional hazards models to assess the risk of asthma onset in relation to pollution exposures, while stratifying for age in years and sex in the baseline hazard. Community-level random effects allowed for clustering and assessment of residual community variation in time to asthma onset (Jerrett et al. 2005; Ma et al. 2003). The model took the following form:

    where hij(t) is the hazard function for the ith subject in jth community, h0s(t) is the baseline hazard function for stratum s (i.e., age and sex), ηj represents the positive random effects for community j with expectation 1 and variance σ2, Xij is the air pollution exposure (i.e., NO2) for individual i in community j, and Zij represents covariates (e.g., health insurance) for individual i in community j.

    We evaluated pollution exposures with other covariates such as relative humidity and medical insurance, which we treated as confounders if they had an association with asthma incidence and changed the pollution coefficient by ≥ 5% for the annual NO2. We report results across the average within-community interquartile range of annual average NO2 exposure.

    To determine whether exposures could be assessed at the individual level across the entire range, we tested for significant differences in effects from community aggregated means from the individually assigned exposure deviated from those community means. The model that allows for separate between-community means and the within-community deviations had the following form:

    where Xij is the individually assigned exposure (i.e., NO2) and Xj is the community average of the individual measurements. We then assessed the equivalence of these two pollution effects by testing the hypothesis H0: β1 = β2. If the two exposure β values were not statistically different based on the log-likelihood ratio test, we concluded that there was no need to distinguish between the within-community and between-community effects of pollution.

    We examined this alternative model to determine whether the individual exposures specified in Equation 1 were appropriate or whether we needed a model specifying different within- and between-community effects for correct interpretation of results.

    Results

    Table 1 lists the descriptive statistics and bivariate associations with asthma for the sample and all the variables used in the analysis. Hispanic Americans and those of white race were well represented in this sample, but there were few Asian Americans, African Americans, or others. Males represented just about 40% (86) of this sample. Many of the children had potential risk factors for asthma in the home, with 95% (183) having carpet in the child’s bedroom, 83% (181) having a pet, and 79% (168) having a gas stove. About 17% (35) of the children had a parental history of asthma. The median age at baseline was 9.6 years. Except for a higher proportion of females, there were no significant differences between the frequency of these variables among the subjects of this study and the remainder of eligible subjects from which we drew the study sample. The study sample and other eligible subjects also had similar traffic exposures (results not shown).

    Table 1 Subject characteristics and bivariate associations with asthma onset.

    Risk factorNo. (%)aHazard ratio (95% CI)b
    Subject characteristics
    Race/ethnicity
     Hispanic ethnicity66 (32)2.15 (1.00–4.58)
     African-American3 (1.5)NAc
     Asian14 (6.8)1.05 (0.20–5.61)
     White non-Hispanicd120 (59)1.00
     Mixed/other race14 (6.5)0.42 (0.04–4.57)
    BMIe
     Below 10th percentile27 (12)0.54 (0.14–2.08)
     10th–90th percentiled177 (82)1.00
     90th percentile13 (6.0)2.46 (0.77–7.83)
    Male sex86 (40)0.89 (0.42–1.88)
    Hay feverf27 (13)1.16 (0.40–3.35)
    Ever wheeze52 (26)0.99 (0.40–2.46)
    Medical care and SEP
    Medical insurance coverage187 (87)2.44 (0.57–10.49)
    Parental educationg
     High school or less69 (33)1.19 (0.51–2.73)
     Some colleged87 (41)1.00
     College or greater54 (26)0.93 (0.35–2.45)
    Home characteristics
    Mildew in home59 (29)0.59 (0.23–1.51)
    Carpet in child’s bedroom183 (95)0.76 (0.30–1.91)
    Plants in home72 (34)1.12 (0.52–2.41)
    Pets at home181 (83)0.80 (0.32–2.02)
    Gas stove in home168 (79)1.52 (0.54–4.27)
    Current daily smoker in home29 (13)0.24 (0.03–1.65)
    Family history
    Maternal smoking during pregnancy24 (13)0.95 (0.29–3.07)
    Parental history of asthmah35 (17)1.31 (0.51–3.35)

    Abbreviations: CI, confidence interval; NA, not applicable.

    aTotal subject number = 217; subject numbers in the table vary because of missing covariate values.

    bCrude estimates with baseline strata for age and sex and community random effects, except model for male sex.

    cCould not derive hazard ratio because of small counts in this stratum.

    dReference group, total subject number = 217; subject numbers in the table vary because of missing covariate values.

    eAge- and sex-adjusted body mass index.

    fSymptoms in the preceding 12 months.

    gCategories for parental educational attainment were less than 12th grade or completion of high school, some college or technical school, and completed 4 years of college or some graduate training after college.

    hMedical diagnosis of asthma in either biologic parent.

    The bivariate associations of individual and household variables as predictors of asthma were largely not statistically significant (Table 1). Associations with body mass index (BMI) and other covariates were generally consistent with associations observed in previous analyses of the entire cohort from which we drew this study sample (Gilliland et al. 2003; McConnell et al. 2002a). In this sample, Hispanics were more likely to report incident asthma during the follow-up than were whites. Because of small counts, we could not derive stable estimates for African Americans. Less than a third of the homes had mildew or a smoker, which may have led to unstable estimates.

    Table 2 shows the contextual confounders tested in the modeling process. In general, lower community SEP was associated with higher risk for asthma onset. Higher winter temperature was associated with higher risk for asthma. The effects of temperature were less consistent than those for humidity, which was a risk factor in all seasons. Higher relative humidity was generally associated with higher risk for incident asthma.

    Table 2 Community SEP and meteorologic characteristics.

    Risk factoraMean ± SDHazard ratio (95% CI)b
    SEP characteristics
    Percent low education20.07 ± 5.901.07 (0.97–1.18)
    Percent male unemployment4.45 ± 1.051.21 (0.71–2.07)
    Median household income39,498 ± 7,0030.44 (0.17–1.12)
    Percent poverty9.15 ± 3.201.11 (0.96–1.29)
    Climatic characteristics
    Mean daily average relative humidity (%)
     Fall–winter months63.08 ± 14.131.02 (0.99–1.06)
     Summer months56.89 ± 15.401.03 (1.00–1.06)
     Annual59.59 ± 13.761.03 (1.00–1.07)
    Mean daily maximum relative humidity (%)
     Fall–winter months78.07 ± 13.461.03 (1.00–1.07)
     Summer months75.23 ± 15.671.03 (1.00–1.07)
    Mean daily average temperature (°C)
     Fall–winter months13.67 ± 2.151.25 (0.82–1.93)
     Summer months18.77 ± 2.390.90 (0.72–1.13)
     Annual16.23 ± 1.930.96 (0.67–1.38)
    Mean daily minimum temperature (°C)
     Fall–winter months9.12 ± 2.791.29 (0.97–1.72)
     Summer months13.35 ± 1.791.09 (0.78–1.54)

    CI, confidence interval.

    aWe derived SEP characteristics using data from the 1990 U.S. Census, and meteorologic characteristics using daily data in 1995 from local meteorologic stations in each community.

    bAll models feature baseline hazard stratified for age and sex, and random effects for community of residence. Estimates are per one unit change in all cases (i.e., % or °C), except for median household income, which is per $10,000.

    In Figure 1, the bars illustrate individual residential NO2 levels within each study community. Variation exists within communities, and a broader regional pattern is also present. Communities in the Los Angeles metropolitan area have generally higher levels than do communities to the north and south, showing a spatial pattern consistent with earlier measurements based on central monitors (Gauderman et al. 2004). We examined the incidence rate by tertile of exposure. Dividing the sample into tertiles of annual NO2, the crude asthma incidence rates, from lowest to highest tertile of exposure, were 14.2, 19.1, and 20.7 per 1,000 person-years, respectively [with wide confidence intervals (CIs)].

    Figure 1 Measured NO2 at subjects’ homes in the 11 study communities: summer/winter average.

    Table 3 presents the descriptive statistics for the exposure variables within each community. There was a wide range of average exposure across communities but also substantial variation within communities (mean interquartile range of 6.2 ppb). Communities in the Los Angeles metropolitan area (Lake Elsinore, Long Beach, Mira Loma, Riverside, San Dimas, and Upland) generally had higher NO2.

    Table 3 Descriptive statistics for exposure to measured NO2 pollution (ppb) by community (mean ± SD).

    CommunityNo.aFall–WinterSummerAnnual
    Alpine1917.5 ± 4.119.9 ± 5.118.7 ± 3.8
    Atascadero1215.3 ± 8.811.0 ± 4.013.3 ± 5.7
    Lake Elsinore2227.3 ± 3.317.6 ± 3.122.5 ± 3.1
    Lancaster1321.0 ± 4.115.8 ± 3.918.4 ± 3.9
    Lompoc2713.8 ± 4.05.4 ± 1.69.6 ± 2.5
    Long Beach2250.0 ± 10.134.1 ± 5.641.5 ± 6.5
    Mira Loma1548.4 ± 3.837.7 ± 4.243.1 ± 2.9
    Riverside2642.7 ± 7.738.3 ± 7.540.7 ± 6.3
    San Dimas3050.6 ± 6.251.3 ± 5.751.3 ± 4.4
    Santa Maria1717.4 ± 3.111.9 ± 3.614.8 ± 2.8
    Upland1435.5 ± 6.146.3 ± 7.140.8 ± 6.4

    aNumber of subjects varied because of invalid measurements during fall–winter and summer months.

    We tested each of the individual variables shown in Tables 1 and 2 with an association to asthma to determine whether confounding reduced the pollution coefficient by ≥ 5%. Table 4 contains the results of the final models, in which we retained Hispanic ethnicity, medical insurance, enrollment group, and relative humidity as confounding adjustment variables. All pollution metrics were positively associated with incident asthma. For average residential NO2 the hazard ratio (HR) was 1.29 (95% CI, 1.07–1.56) across the average within-community interquartile range of 6.2 ppb. The fall–winter HR of 1.29 (95% CI, 1.11–1.49) was nearly the same as the annual average, whereas the summer HR of 1.27 (95% CI, 1.03–1.59) was slightly smaller. The sample size was 196 with 26 cases (for the annual NO2), so it is possible that influential outliers might have affected the results. However, there were no cases with exceptionally high exposure to NO2 and no noncases with notably low exposure to NO2.

    Table 4 Association between incident asthma and measured NO2 pollution with varying levels of control for confounding.

    Model 3: pollution model without relative humidity
    Model 4: pollution model with relative humidity
    ExposureNo.aNull model (σ2b)Model 1: individual covariates (σ2b)Model 2: individual covariates,relative humidity (σ2b)HR (95% CI)cσ2bHR (95% CI)cσ2b
    Measured NO2 pollution
     Fall–Winter2090.571610.523190.166681.17 (0.98–1.41)0.269241.29 (1.11–1.49)0.00016
     Summer2040.337050.264210.077701.03 (0.86–1.23)0.251221.27 (1.03–1.57)0.00020
     Annual1960.339280.272780.073031.10 (0.91–1.33)0.217051.29 (1.07–1.56)0.00018

    Null model contains only community random effects. Model 1 contains Hispanic ethnicity, enrollment group, and medical insurance coverage, with baseline strata for age and sex. Model 2 we further adjusted for community annual mean relative humidity. Model 3 we further adjusted for pollution but did not adjust for community annual mean relative humidity. Model 4 we further adjusted for community annual mean relative humidity.

    aNumber of subjects varied because of invalid measurements during winter and summer months.

    bValues of σ2 are the random-effects variance estimates, which represent the residual variation in asthma incidence attributable to the community cluster after we include all other fixed predictors in the model.

    cMeasured NO2 pollution estimates are over a 6.2-ppb exposure contrast, which is the average within-community interquartile range for average annual measured NO2.

    The σ2 values are the random-effects variance estimates, which represent the residual variation in asthma incidence attributable to the community cluster after we included all other fixed predictors in the model (Table 4). Addition of the individual covariates in model 1 reduced the σ2 only slightly from a null model containing only baseline age and sex strata (results not shown). Models 2 and 3 show the large decline in σ2 when we included community relative humidity and pollution, respectively, indicating that both variables account for a large proportion of the residual between-community variation in asthma incidence.

    We assessed exposures both as the deviation of the individual specific exposure metric from the corresponding community mean and as the community average (i.e., average of household NO2 per community). We used these variables to compare within- and between-community pollution effects. Table 5 lists results of the tests shown in Equation 2 for difference in the between-community exposures and deviation within community exposures. The between-community summer NO2 point estimate was larger than the within-community estimate, but these estimates had large CI values and were not significantly different based on the likelihood ratio test. For yearly average and winter NO2, the within- and between-community effects were similar.

    Table 5 Association between incident asthma and measured NO2 pollution: within-community versus between-community models.

    Measured NO2 pollutionNo.aHR (95% CI)bσ2cLikelihood ratio test p-valued
    Fall–winter2091.29 (1.11–1.49)0.00016
     Within-community deviation2091.32 (0.91–1.92)0.00016> 0.1
     Community mean1.28 (1.09–1.51)
    Summer2041.27 (1.03–1.57)0.00020
     Within-community deviation2040.99 (0.57–1.72)0.00019> 0.1
     Community mean1.32 (1.05–1.66)
    Annual1961.29 (1.07–1.56)0.00018
     Within-community deviation1961.31 (0.76–2.26)0.00019> 0.1
     Community mean1.28 (1.05–1.57)

    aNumber of subjects varied because of invalid measurements during fall–winter and summer months.

    bWe adjusted models for Hispanic ethnicity, enrollment group, medical insurance coverage, and community annual mean relative humidity, with baseline strata for age and sex. Measured NO2 pollution estimates are over a 6.2-ppb exposure contrast, which is the average within-community interquartile range for average annual measured NO2.

    cValues of σ2 are the random-effects variance estimates, which represent the residual variation in asthma incidence attributable to the community cluster after we included all other fixed predictors in the model.

    dLikelihood ratio tests for the null hypothesis that there is no difference between the random-effects models and between/within-community models were not rejected.

    Interactions between humidity and air pollution were tested, but found not to be significant, although pollution did appear to exert larger effects in more humid areas (t-values were between 0.9 and 1.3 depending on the humidity interaction term with various measures of pollution).

    We conducted further sensitivity analyses to assess whether excluding early wheeze at baseline affected the results (Table 6). Excluding baseline wheeze had a small impact on the asthma effect of air pollution. We implemented another sensitivity test by removing data for all children with a history of early chest illness, including croup, bronchitis, bronchiolitis, or pneumonia (Table 6). The HR values were mildly attenuated but still elevated and similar to the HR values without these exclusions.

    Table 6 Association between incident asthma and measured NO2 pollution with various exclusions to test for confounding by undiagnosed asthma and early chest illness at study baseline.

    All cases (n = 217)
    Excluding subjects with wheeze at study baseline (n = 152)
    Excluding subjects with early childhood chest illness (n = 180)
    Measured NO2 pollutionNo.aHR (95% CI)bNo.a(95% CI)bNo.a(95% CI)b
    Fall–winter2091.29 (1.11–1.49)1461.29 (1.06–1.51)1631.22 (1.04–1.45)
    Summer2041.27 (1.03–1.57)1441.21 (0.89–1.65)1731.26 (1.00–1.58)
    Annual1961.29 (1.07–1.56)1381.28 (1.00–1.64)1701.26 (1.02–1.55)

    aNumber of subjects varied because of invalid measurements during fall–winter and summer months.

    bMeasured NO2 pollution estimates are over a 6.2-ppb exposure contrast, which is the average within-community interquartile range for average annual measured NO2.

    Discussion

    In this prospective cohort of children from 11 Southern California communities, we found significant associations between incident asthma and measured outdoor residential exposure to NO2. A strength of the study was the individual assignment of exposure, based on measurements of NO2 outdoors at the homes of the children. Previous research on incident asthma and air pollution has largely relied on modeled exposure (Brauer et al. 2002; Zmirou et al. 2004) or on exposure measured at central site locations (McConnell et al. 2002a; McDonnell et al. 1999; Shima et al. 2002). To our knowledge, this is the first study to employ individually measured pollutant exposures with incident disease assessment. Although the results must be interpreted in light of the limited sample size and small number of cases, they highlight a need for further examination of relationships between incident disease and individually measured exposures.

    We tested numerous individual and contextual confounding variables. Some of the variables that had an association with asthma were not confounders because they exerted little or no effect on the pollution coefficient. Parsimonious models controlling for age and sex in the baseline hazard indicated that Hispanic ethnicity, medical insurance, and relative humidity were confounders in models with NO2. The air pollution exposures generally accounted for larger portions of the residual variance in the community random effects (i.e., σ2) than did the individual-level variables included as confounders. The results for residential NO2 were sensitive to inclusion of the relative humidity variable. This variable met our criteria for confounding and increased the pollution effect (Table 4). Because adjustment for humidity elevated the risks associated with air pollution exposure, we further investigated an interaction between pollution exposures and humidity. Although the interaction was not significant, we did observe generally higher risks for pollution in more humid communities.

    In the absence of interaction, the influence of humidity may represent a case of variance suppression (Tabachnick and Fidell 1985), whereby confounding variables may increase the effect size of the primary risk variable in the model. In this case, relative humidity and the pollutants share some variance, and this shared component also partially overlaps with the asthma outcome. The part of the variance shared by relative humidity and pollution is the component of variance with asthma that has a relatively less precise association with pollution. With the inclusion of relative humidity, this less precise component of the variance overlap between pollution and asthma is removed, and the remaining relationship becomes more pronounced.

    We further investigated the possibility of positive confounding by examining the correlation between humidity and NO2 and the impact of the humidity risk with NO2 excluded from the model. There was a moderately high, significant negative correlation between NO2 and relative humidity. Evaluated at the household level for the annual estimates, the Pearson correlation coefficient was −0.57, and at the community average level of NO2 it was −0.61. Except for Long Beach, most of the coastal communities had lower air pollution and higher humidity. Because humidity was positively associated with asthma but negatively associated with NO2, we would expect humidity to be a negative confounder (i.e., excluding humidity from the model should result in a negative bias on the NO2 estimate). In this case, we did observe a negative bias when humidity was not included in the model. We would also expect a negative bias in the humidity effect when we exclude NO2 from the model. Table 7 shows the humidity risks with and without NO2 in the model. Without NO2, the risks were attenuated in all models, and the negative bias was strongest in the summer season. Thus, the impact of humidity was negatively biased when we excluded NO2 from the model, which is consistent with the positive confounding effect observed in the NO2 models. These results imply that simultaneous consideration of meteorology and air pollution may be important in models assessing the effects of chronic air pollution exposure on respiratory disease.

    Table 7 Association between incident asthma and mean daily average relative humidity (annual) with and without control for measured NO2 pollution.

    HR (95% CI)a
    Mean daily average relative humidity (annual)No.bModel 1: not adjusted for measured NO2 pollutionModel 2: adjusted for measured NO2 pollution
    Fall–winter2092.45 (1.11–5.43)3.75 (1.85–7.62)
    Summer2041.98 (0.96–4.08)4.22 (1.58–11.28)
    Annual1962.09 (1.00–4.35)3.90 (1.70–8.97)

    aWe scaled HR values and 95% CIs across the interquartile range in relative humidity across all communities (23.3%). Models contain adjustment for Hispanic ethnicity, enrollment group, and medical insurance coverage, with baseline strata for age and sex. bNumber of subjects varied because of invalid measurements during winter and summer months.

    Interpreting the household NO2 findings is complicated because it represents a mixture of local and regionally transported pollutants in Southern California. In an earlier analysis, Gauderman et al. (2005) reported moderate to high correlations between the measured residential NO2 and various measures of traffic proximity or modeled concentrations. Although some toxicologic evidence supports the potential for NO2 as a direct causative agent for respiratory health effects at high concentrations (Kelly and Sandstrom 2004), the more likely interpretation is that NO2 was a proxy for other components of the local air pollution mixture with high toxic potential. Because NO2 pollutants may be transported from neighboring areas, some of the effects observed could represent local and regional contributions. If NO2 represents other toxic traffic-related pollutants, the transported component would have different toxicologic properties than purely local contributions.

    Rather than the specific causal agent, NO2 may represent a mixture of local traffic-related pollution and regional pollutants, which in Southern California include transported primary traffic-related pollutants and regional products of photochemistry. The design of the study allowed us to examine the independent effects of each of these contributions to the pollutant mixture by modeling the effects of the within- and between-community variation in NO2, as shown in Table 5. The within-community effects indicative of long-term local traffic sources were similar in magnitude to effects of community average NO2 across communities, suggesting that both regional and local pollution contributed to the associations with asthma. However, the range of variation of NO2 within communities was smaller than between communities, and the HR values were therefore smaller when constrained to the range within communities. For example, using the average interquartile range across all measurements of 28.9 ppb for annual NO2 increased the HR to 3.25. For the average within-community range of 16.4 ppb, the HR was 1.95.

    Compared with earlier studies on prevalent asthma, with a similar population and the same exposure measurements for NO2, the odds ratio per interquartile range increment of 5.7 ppb was 1.83 for prevalent asthma (Gauderman et al. 2005). Thus, the incident effects were smaller than those observed for prevalent asthma over a similar exposure contrast. The overall effects we observed here nonetheless support results from our earlier prevalence study. In both studies, we observed effects on asthma from NO2.

    Notable strengths of this study are the prospective design, individually measured pollutant exposures, incident disease assessment, and the large effects. Limitations include an inability to identify the specific constituents of pollution that were responsible for the observed health effects, the relatively small sample size, and questionnaire case ascertainment. Although a potential limitation, self-reported physician diagnosis is widely used to define asthma in epidemiologic studies (Burr 1992), is reproducible (Ehrlich et al. 1995; Peat et al. 1992), and reflects what physicians actually reported to patients (Burney et al. 1989; Greer et al. 1993). We conducted sensitivity analyses with alternative case definitions, excluding children with baseline wheeze and early childhood chest illness. The summer estimates were attenuated mildly (HR values range from 1.21 to 1.26), and after excluding baseline wheeze became insignificant (Table 6). Many of the communities had lower NO2 levels in the summer (Table 3), and the effects during the summer season were smaller than in the other seasons. In contrast, the HR values for NO2 were only slightly smaller for the fall–winter and annual estimates and remained significant (HR values range from 1.2 to 1.9). The general conclusions, therefore, are similar even after excluding children with wheeze at baseline or with a history of early childhood illness.

    The overall asthma incidence rate per thousand person-years of 18.8 cases resulting from parent- or child-reported physician diagnosis was high but is similar to rates reported previously in the CHS cohort (McConnell et al. 2002b) and in other studies of asthma or wheeze in school children (Lombardi et al. 1997; Ronmark et al. 2001). Lower rates have been found in other studies (Broder et al. 1974; Dodge and Burrows 1980; Yunginger et al. 1992). Comparison across studies, however, is complicated by varying ascertainment methods and temporal trends in incidence rates. The rate in this sample from the larger CHS cohort was similar to those in comparable studies.

    Further field studies measuring plausibly causal ultrafine particles, metals, and poly-cyclic aromatic hydrocarbons simultaneously (Nel 2005) might identify the constituents responsible for the observed association with NO2. Although the specific pollutant or related mixture of constituents responsible for the effects remains the focus of future research, the large risks observed suggest that traffic-related air pollution exposure contributes to new-onset asthma.

    This work was funded by the National Institute of Environmental Health Sciences (grants 5P30ES007048, 5P01ES011627, 5P01ES009581, and 1R03ES014046), the U.S. Environmental Protection Agency (grants RD831861, R826708, and R831845), California Air Resources Board contract 94-331, and the Hastings Foundation.

    References

    • Brauer M, Hoek G, Smit HA, de Jongste JC, Gerritsen J, Postma DSet al.. 2007. Air pollution and development of asthma, allergy and infections in a birth cohort. Eur Respir J 29(5):879-88817251230. Crossref, MedlineGoogle Scholar
    • Brauer M, Hoek G, Van Vliet P, Meliefste K, Fischer PH, Wijga Aet al.. 2002. Air pollution from traffic and the development of respiratory infections and asthmatic and allergic symptoms in children. Am J Respir Crit Care Med 166(8):1092-109812379553. Crossref, MedlineGoogle Scholar
    • Broder I, Higgins MW, Mathews KP, Keller JP. 1974. Epidemiology of asthma and allergic rhinitis in a total community, Tecumseh, Michigan. IV. Natural history. J Allergy Clin Immunol 54(2):100-1104851886. Crossref, MedlineGoogle Scholar
    • Brunekreef B, Janssen NA, de Hartog J, Harssema H, Knape M, van Vliet P. 1997. Air pollution from truck traffic and lung function in children living near motorways. Epidemiology 8(3):298-3039115026. Crossref, MedlineGoogle Scholar
    • Brunekreef B, Sunyer J. 2003. Asthma, rhinitis and air pollution: is traffic to blame? [Comment]. Eur Respir J 21(6):913-91512797480. Crossref, MedlineGoogle Scholar
    • Burney PG, Laitinen LA, Perdrizet S, Huckauf H, Tattersfield AE, Chinn Set al.. 1989. Validity and repeatability of the IUATLD (1984) bronchial symptoms questionnaire: an international comparison. Eur Respir J 2(10):940-9452606194. Crossref, MedlineGoogle Scholar
    • Burr ML. 1992. Diagnosing asthma by questionnaire in epidemiological surveys. Clin Exp Allergy 22(5):509-5101628247. Crossref, MedlineGoogle Scholar
    • Ciccone G, Forastiere F, Agabiti N, Biggeri A, Bisanti L, Chellini Eet al.. 1998. Road traffic and adverse respiratory effects in children. SIDRIA Collaborative Group. Occup Environ Med 55(11):771-7789924455. Crossref, MedlineGoogle Scholar
    • Delfino RJ. 2002. Epidemiologic evidence for asthma and exposure to air toxics: Linkages between occupational, indoor, and community air pollution research. Environ Health Perspect 110(suppl 4):573-58912194890. LinkGoogle Scholar
    • Dodge RR, Burrows B. 1980. The prevalence and incidence of asthma and asthma-like symptoms in a general population sample. Am Rev Respir Dis 122(4):567-5757436122. MedlineGoogle Scholar
    • Duhme H, Weiland S, Keil U, Kraemer B, Schmid M, Stender Met al.. 1996. The association between self-reported symptoms of asthma and allergic rhinitis and self-reported traffic density on street residence in adolescents. Epidemiology 7(6):578-5828899382. Crossref, MedlineGoogle Scholar
    • Eder W, Ege MJ, von Mutius E. 2006. The asthma epidemic. N Engl J Med 355(21):2226-223517124020. Crossref, MedlineGoogle Scholar
    • Edwards J, Walters S, Griffiths RK. 1994. Hospital admissions for asthma in preschool children: relationship to major roads in Birmingham, United Kingdom. Arch Environ Health 49(4):223-2277518223. Crossref, MedlineGoogle Scholar
    • Ehrlich RI, Du Toit D, Jordaan E, Volmink JA, Weinberg EG, Zwarenstein M. 1995. Prevalence and reliability of asthma symptoms in primary school children in Cape Town. Int J Epidemiol 24(6):1138-11458824855. Crossref, MedlineGoogle Scholar
    • Gauderman WJ, Avol E, Gilliland F, Vora H, Thomas D, Berhane Ket al.. 2004. The effect of air pollution on lung development from 10 to 18 years of age. N Engl J Med 351(11):1057-106715356303. Crossref, MedlineGoogle Scholar
    • Gauderman WJ, Avol E, Lurmann F, Kuenzli N, Gilliland F, Peters Jet al.. 2005. Childhood asthma and exposure to traffic and nitrogen dioxide. Epidemiology 16(6):737-74316222162. Crossref, MedlineGoogle Scholar
    • Gilliland FD, Berhane K, Islam T, McConnell R, Gauderman WJ, Gilliland SSet al.. 2003. Obesity and the risk of newly diagnosed asthma in school-age children. Am J Epidemiol 158(5):406-41512936895. Crossref, MedlineGoogle Scholar
    • Gold DR, Wright R. 2005. Population disparities in asthma. Annu Rev Public Health 26:89-11315760282. Crossref, MedlineGoogle Scholar
    • Greer JR, Abbey DE, Burchette RJ. 1993. Asthma related to occupational and ambient air pollutants in nonsmokers. J Occup Med 35(9):909-9158229343. Crossref, MedlineGoogle Scholar
    • Guo YL, Lin YC, Sung FC, Huang SL, Ko YC, Lai JSet al.. 1999. Climate, traffic-related air pollutants, and asthma prevalence in middle-school children in Taiwan. Environ Health Perspect 107:1001-100610585904. LinkGoogle Scholar
    • Haan M, Kaplan GA, Camacho T. 1987. Poverty and health. Prospective evidence from the Alameda County Study. Am J Epidemiol 125(6):989-9983578257. Crossref, MedlineGoogle Scholar
    • Jerrett M, Burnett RT, Ma R, Pope CA, Krewski D, Newbold KBet al.. 2005. Spatial analysis of air pollution and mortality in Los Angeles. Epidemiology 16(6):727-73616222161. Crossref, MedlineGoogle Scholar
    • Jin RL, Shah CP, Svaboda TJ. 1994. The Health Impact of Unemployment: A Review and Application of Research Evidence. Working paper for the Population Health CommitteeToronto, ONOntario Medical Association. Google Scholar
    • Juhn YJ, Sauver JS, Katusic S, Vargas D, Weaver A, Yunginger J. 2005. The influence of neighborhood environment on the incidence of childhood asthma: a multilevel approach. Soc Sci Med 60(11):2453-246415814171. Crossref, MedlineGoogle Scholar
    • Kelly FJ, Sandstrom T. 2004. Air pollution, oxidative stress, and allergic response. Lancet 363(9403):95-9614726159. Crossref, MedlineGoogle Scholar
    • Lewis SA, Antoniak M, Venn AJ, Davies L, Goodwin A, Salfield Net al.. 2005. Secondhand smoke, dietary fruit intake, road traffic exposures, and the prevalence of asthma: a cross-sectional study in young children. Am J Epidemiol 161(5):406-41115718476. Crossref, MedlineGoogle Scholar
    • Livingstone AE, Shaddick G, Grundy C, Elliott P. 1996. Do people living near inner city main roads have more asthma needing treatment? Case control study. BMJ 312(7032):676-6778597735. Crossref, MedlineGoogle Scholar
    • Lombardi E, Morgan WJ, Wright AL, Stein RT, Holberg CJ, Martinez FD. 1997. Cold air challenge at age 6 and subsequent incidence of asthma. A longitudinal study. Am J Respir Crit Care Med 156(6):1863-18699412567. Crossref, MedlineGoogle Scholar
    • Ma R, Krewski D, Burnett RT. 2003. Random effects Cox models: a Poisson modelling approach. Biometrika 90(1):157-169. CrossrefGoogle Scholar
    • McConnell R, Berhane K, Gilliland F, Islam T, Gauderman WJ, London SJet al.. 2002a. Indoor risk factors for asthma in a prospective study of adolescents. Epidemiology 13(3):288-29511964930. Crossref, MedlineGoogle Scholar
    • McConnell R, Berhane K, Gilliland F, London SJ, Islam T, Gauderman WJet al.. 2002b. Asthma in exercising children exposed to ozone: a cohort study. Lancet 359(9304):386-39111844508. Crossref, MedlineGoogle Scholar
    • McConnell R, Berhane K, Yao L, Jerrett M, Lurmann F, Gilliland Fet al.. 2006. Traffic, susceptibility, and childhood asthma. Environ Health Perspect 114:766-77216675435. LinkGoogle Scholar
    • McDonnell WF, Abbey DE, Nishino N, Lebowitz MD. 1999. Long-term ambient ozone concentration and the incidence of asthma in nonsmoking adults: the AHSMOG study. Environ Res 80(2 pt 1):110-12110092402. Crossref, MedlineGoogle Scholar
    • Migliaretti G, Cavallo F. 2004. Urban air pollution and asthma in children. Pediatr Pulmonol 38(3):198-20315274097. Crossref, MedlineGoogle Scholar
    • Molitor J, Molitor NT, Jerrett M, McConnell R, Gauderman J, Berhane Ket al.. 2006. Bayesian modeling of air pollution health effects with missing exposure data. Am J Epidemiol 164(1):69-7616624966. Crossref, MedlineGoogle Scholar
    • Nel A. 2005. Atmosphere. Air pollution-related illness: effects of particles. Science 308(5723):804-80615879201. Crossref, MedlineGoogle Scholar
    • Nitta H, Sato T, Nakai S, Maeda K, Aoki S, Ono M. 1993. Respiratory health associated with exposure to automobile exhaust. I. Results of cross-sectional studies in 1979, 1982, and 1983. Arch Environ Health 48(1):53-587680850. Crossref, MedlineGoogle Scholar
    • Oosterlee A, Drijver M, Lebret E, Brunekreef B. 1996. Chronic respiratory symptoms in children and adults living along streets with high traffic density. Occup Environ Med 53:241-2478664961. Crossref, MedlineGoogle Scholar
    • Peat JK, Salome CM, Toelle BG, Bauman A, Woolcock AJ. 1992. Reliability of a respiratory history questionnaire and effect of mode of administration on classification of asthma in children. Chest 102(1):153-1571623744. Crossref, MedlineGoogle Scholar
    • Ronmark E, Jonsson E, Platts-Mills T, Lundback B. 2001. Incidence and remission of asthma in schoolchildren: report from the obstructive lung disease in northern Sweden studies. Pediatrics 107(3):e3711230618. Crossref, MedlineGoogle Scholar
    • Shima M, Nitta Y, Ando M, Adachi M. 2002. Effects of air pollution on the prevalence and incidence of asthma in children. Arch Environ Health 57(6):529-53512696649. Crossref, MedlineGoogle Scholar
    • Sorensen G, Verbrugge LM. 1987. Women, work, and health. Annu Rev Public Health 8:235-2513555524. Crossref, MedlineGoogle Scholar
    • Tabachnick BG, Fidell LS. 1985. Using Multivariate StatisticsSan FranciscoHarper and Row. Google Scholar
    • Tatalovich T, Wilson J, Jerrett M, Milam J, McConnell R. 2006. Competing definitions of contextual environments. Int J Health Geogr 7:55-7017156433. Crossref, MedlineGoogle Scholar
    • U.S. Census Bureau. 1992. 1990 Summary Tape File 3 Technical DocumentationAvailable: http://factfinder.census.gov/metadoc/1990stf3td.pdf[accessed 19 March 2005]. Google Scholar
    • van Vliet P, Knape M, de Hartog J, Janssen N, Harssema H, Brunekreef B. 1997. Motor vehicle exhaust and chronic respiratory symptoms in children living near freeways. Environ Res 74(2):122-1329339225. Crossref, MedlineGoogle Scholar
    • Weiland SK, Husing A, Strachan DP, Rzehak P, Pearce N. 2004. Climate and the prevalence of symptoms of asthma, allergic rhinitis, and atopic eczema in children. Occup Environ Med 61(7):609-61515208377. Crossref, MedlineGoogle Scholar
    • Weiland SK, Mundt KA, Ruckmann A, Keil U. 1994. Self-reported wheezing and allergic rhinitis in children and traffic density on street of residence. Ann Epidemiol 4:243-2477519948. Crossref, MedlineGoogle Scholar
    • Wilkinson P, Elliott P, Grundy C, Shaddick G, Thakrar B, Walls Pet al.. 1999. Case-control study of hospital admission with asthma in children aged 5–14 years: relation with road traffic in north west London. Thorax 54:1070-107410567625. Crossref, MedlineGoogle Scholar
    • Wissow LS, Gittelsohn AM, Szklo M, Starfield B, Mussman M. 1988. Poverty, race, and hospitalization for childhood asthma. Am J Public Health 78(7):777-7823381951. Crossref, MedlineGoogle Scholar
    • Wjst M, Reitmeir P, Dold S, Wulff A, Nicolai T, von Loeffelholz-Colberg EFet al.. 1993. Road traffic and adverse effects on respiratory health in children. BMJ 307(6904):596-6007691304. Crossref, MedlineGoogle Scholar
    • Yunginger JW, Reed CE, O’Connell EJ, Melton LJ, O’Fallon WM, Silverstein MD. 1992. A community-based study of the epidemiology of asthma. Incidence rates, 1964–1983. Am Rev Respir Dis 146(4):888-8941416415. Crossref, MedlineGoogle Scholar
    • Zmirou D, Gauvin S, Pin I, Momas I, Sahraoui F, Just Jet al.. 2004. Traffic related air pollution and incidence of childhood asthma: results of the vesta case-control study. J Epidemiol Community Health 58(1):18-2314684722. Crossref, MedlineGoogle Scholar