Black, dark gray, and light gray horizontal bars represent follow-up periods for the 1993-2001, 1996-2004, and 2006-2014 cohorts, respectively. Follow-up for the 2006-2014 cohort is truncated on the graph at 2011, the last year with air pollution data. ppb indicates parts per billion.
Symbol colors indicate community and the size—small, medium, and large—indicates the data are from the 1993-2001, 1996-2004, and 2006-2014 cohorts, respectively. Simple linear regression models based on asthma incidence and air pollution concentration used to generate regression lines separately for each community. ppb indicates parts per billion.
eMethods. Data Analysis
eTable 1. Distribution of Selected Characteristics Among Participants Included and Excluded From the Current Study of the Children's Health Study, 1993-2014
eTable 2. P Values for Evidence of Effect Modification of the Asthma Incidence and Air Pollution Association Based on Partial Likelihood Ratio Tests for Interaction Terms With Pollutant and Potential Effect Modifier in the Children's Health Study, 1993-2014
eTable 3. Results of Sensitivity Analyses for Incidence Rate Ratios (IRR) Associated With Reduction in Nitrogen Dioxide (NO2) and Particulate Matter <2.5 μm (PM2.5) in the Children's Health Study, 1993-2014
eTable 4. Results of Sensitivity Analyses Using Cox Proportional Hazards Models for Ozone (O3), Nitrogen Dioxide (NO2), Particulate Matter <10 μm (PM10) and <2.5 μm (PM2.5) in the Children's Health Study, 1993-2014
Customize your JAMA Network experience by selecting one or more topics from the list below.
Garcia E, Berhane KT, Islam T, et al. Association of Changes in Air Quality With Incident Asthma in Children in California, 1993-2014. JAMA. 2019;321(19):1906–1915. doi:10.1001/jama.2019.5357
What is the association between reductions in regional air pollutant concentration and incidence of childhood asthma?
In this longitudinal study that included 4140 children, each 4.3–parts-per-billion decrease in nitrogen dioxide was associated with a reduction of 0.83 cases per 100 person-years in asthma incidence; each 8.1-μg/m3 decrease in particulate matter less than 2.5 μm was associated with a reduction of 1.53 cases per 100 person-years in asthma incidence. There were no statistically significant associations with change in ozone and particulate matter less than 10 μm.
Declines in nitrogen dioxide and particulate matter less than 2.5 μm may be associated with decreased childhood asthma incidence.
Exposure to air pollutants is a well-established cause of asthma exacerbation in children; whether air pollutants play a role in the development of childhood asthma, however, remains uncertain.
To examine whether decreasing regional air pollutants were associated with reduced incidence of childhood asthma.
Design, Setting, and Participants
A multilevel longitudinal cohort drawn from 3 waves of the Southern California Children’s Health Study over a period of air pollution decline. Each cohort was followed up from 4th to 12th grade (8 years): 1993-2001, 1996-2004, and 2006-2014. Final follow-up for these data was June 2014. Population-based recruitment was from public elementary schools. A total of 4140 children with no history of asthma and residing in 1 of 9 Children’s Health Study communities at baseline were included.
Annual mean community-level ozone, nitrogen dioxide, and particulate matter less than 10 μm (PM10) and less than 2.5 μm (PM2.5) in the baseline year for each of 3 cohorts.
Main Outcomes and Measures
Prospectively identified incident asthma, collected via questionnaires during follow-up.
Among the 4140 children included in this study (mean [SD] age at baseline, 9.5 [0.6] years; 52.6% female [n = 2 179]; 58.6% white [n = 2273]; and 42.2% Hispanic [n = 1686]), 525 incident asthma cases were identified. For nitrogen dioxide, the incidence rate ratio (IRR) for asthma was 0.80 (95% CI, 0.71-0.90) for a median reduction of 4.3 parts per billion, with an absolute incidence rate decrease of 0.83 cases per 100 person-years. For PM2.5, the IRR was 0.81 (95% CI, 0.67-0.98) for a median reduction of 8.1 μg/m3, with an absolute incidence rate decrease of 1.53 cases per 100 person-years. For ozone, the IRR for asthma was 0.85 (95% CI, 0.71-1.02) for a median reduction of 8.9 parts per billion, with an absolute incidence rate decrease of 0.78 cases per 100 person-years. For PM10, the IRR was 0.93 (95% CI, 0.82-1.07) for a median reduction of 4.0 μg/m3, with an absolute incidence rate decrease of 0.46 cases per 100 person-years.
Conclusions and Relevance
Among children in Southern California, decreases in ambient nitrogen dioxide and PM2.5 between 1993 and 2014 were significantly associated with lower asthma incidence. There were no statistically significant associations for ozone or PM10.
Asthma is the most common pediatric chronic disease, estimated to have affected 14% of children globally in 2002-2003.1 There has been much interest in the effect of outdoor air pollution on asthma risk given its ubiquity and high levels in urban areas, which, compared with rural areas, have higher rates of asthma.2 Globally, outdoor air pollution has been recognized as a major public health concern and was estimated to have contributed 6.8% of the annual disability-adjusted life-years lost in 2016.3 Although ambient air pollution exposure has been causally linked to asthma exacerbations in children,4 evidence has been limited for a role in asthma development.2,4 There has been support for a link with close proximity to busy roads,4-6 but studies of regional pollutants, such as ozone, nitrogen dioxide, and particulate matter (PM), have provided less robust evidence.
This study was designed to take advantage of secular trends in air pollution to examine asthma incidence among children recruited and followed up longitudinally within the same set of communities over a period of air pollution decline. Since the early 1990s, air pollutant concentrations have decreased in Southern California.7 During this time, several successive cohorts of schoolchildren were enrolled from the same set of communities as part of the Southern California Children’s Health Study (CHS), a long-term study of cardiopulmonary pediatric health outcomes.8Quiz Ref ID Leveraging this unique data resource, the study examined whether observed reductions in regional air pollutants, specifically ozone, nitrogen dioxide, and PM less than 10 μm (PM10) or less than 2.5 μm (PM2.5), were associated with asthma incidence rates within these CHS communities.
All parents or guardians of participating children provided written informed consent. The study protocol was approved by the institutional review board of the University of Southern California. The study population was drawn from 3 successive cohorts of the CHS. The CHS has been described in detail elsewhere.8,9 Briefly, 12 communities in Southern California were selected in 1993 based on historical air pollutant levels.8 Children were recruited from participating communities through public schools and followed up prospectively until 12th grade. Fourth graders aged 9 to 10 years were recruited in 1993 (n = 1798) and in 1996 (n = 2061) from the 12 communities. In 2003, kindergarteners and first graders (n = 5736) were recruited from 13 communities, resulting in a total of 16 communities contributing data to the 3 cohorts.9 These 3 cohorts will hereafter be referred to as the 1993-2001, 1996-2004, and 2006-2014 cohorts. Nine communities participated in all 3 cohorts (Alpine, Lake Elsinore, Lake Gregory, Long Beach, Mira Loma, Riverside, San Dimas, Santa Maria, and Upland; n = 6858).
Baseline questionnaires regarding the children’s health and exposures, as well as demographic information, were completed by the parents or guardians. Annual follow-up questionnaires assessing changes in the children’s health, among other factors, were completed initially by the parents or guardians and later by the participating children starting at approximately age 11 years. To better facilitate the comparison of these 3 cohorts in this analysis, follow-up for the 2006-2014 cohort was realigned to begin in the fourth year (2006) of that cohort’s original timeline, when most 2006-2014 cohort participants (>46%) were in the fourth grade. Final follow-up for these data was June 2014 and the data reported here are the most recent available at the time of this study from these cohorts.
Quiz Ref IDAmbient air pollutant monitoring stations were established in each of the study communities and have been continuously measuring regional air pollution since the inception of the CHS. Data on concentrations of ozone, nitrogen dioxide, PM10, and PM2.5 were routinely collected, as previously described.8,10 Community-specific annual mean concentrations in the baseline year for each cohort (ie, 1993, 1996, and 2006 for the 1993-2001, 1996-2004, and 2006-2014 cohorts, respectively) were calculated based on 24-hour means for nitrogen dioxide, PM10, and PM2.5, and on the 10 am to 6 pm mean for ozone due to its marked diurnal variation.
Data were obtained from the baseline questionnaires on children’s date of birth, sex, race and ethnicity, history of asthma, participation in team sports, presence of a gas stove in the home, exposure to smoking in utero, exposure to secondhand smoke, parental education, parental history of asthma, and residential address. Determination of children’s race/ethnicity was made by the parents or guardians who completed the baseline questionnaire based on a question with fixed categories. Race/ethnicity was included as a covariate due to its role as a potential confounder of the air pollution–asthma relationship. Exposure to secondhand smoke was classified based on a positive response to either of the following questions: “Does anyone living in this child’s home currently smoke cigarettes, cigars or pipes on a daily basis inside the home?” or “In the past, has anyone living in this child’s home ever smoked cigarettes on a daily basis inside the home while the child was living there?” Use of a Spanish-language questionnaire by the parent or guardian at baseline was also recorded. Genetic ancestry, including Native American ancestry, was estimated using 233 ancestral informative markers and the STRUCTURE program.11-13
Baseline residential address was used to estimate exposure to local near-roadway pollution, based on a line source dispersion model as previously described.6,9 Inputs to this model included distance to roadways, vehicle counts, vehicle nitrogen oxide emission rates, wind speed and direction, and height of the mixing layer in each community. Mean temperature data were also collected during follow-up from measurements at monitoring stations. For each cohort baseline year, community-specific mean temperature was computed from monthly means in that year. Temperature data were not available for 1993, consequently 1994 levels were used instead for all 1993-2001 cohort participants. Additionally, 2006 temperature data were not available for 2 communities (Lake Gregory and San Dimas); therefore, data from 2001, the closest previous year with complete data, were used for 2006-2014 cohort participants in these 2 communities. Community-level annual mean temperatures did not vary greatly from year to year; the mean coefficient of variation from 1994 to 2001 and 2006, for those communities with data, was 4% (range, 2%-6%).
Quiz Ref IDIncident asthma was defined as a newly reported physician-diagnosed case of asthma on an annual questionnaire during follow-up (ie, first time answered “yes” to the question “Has a doctor ever diagnosed this child as having asthma?” when the parent or guardian was asked or “Has a doctor ever said you have asthma?” when the child was asked). Because incident asthma cases were defined with these annual questionnaires, specific dates of diagnoses were unknown. We imputed the date of diagnosis using the midpoint of the interval between the date of the questionnaire on which asthma diagnosis was first reported and the date of the questionnaire prior to reporting asthma status. This date was used for calculating follow-up time for all statistical analyses. Children with missing questionnaires during follow-up continued to contribute person-time until they reported an asthma diagnosis or were lost to follow-up.
To assess association between changes in regional air quality and asthma incidence in children within community over the course of follow-up for these 3 cohorts, we fitted multilevel Poisson regression models to estimate asthma incidence rate ratios (IRRs) and 95% CIs associated with exposure to regional air pollution.10,14,15 Models included an offset term for person-time (natural log-transformed) and a fixed effect for community. Additionally, to account for clustering effects of children by cohort and community, a random effect for cohort nested within community with an unstructured covariance matrix was included in the model. Follow-up time was calculated as the number of days between joining the cohort (ie, baseline questionnaire date) and either imputed date of asthma diagnosis or date of last completed questionnaire (either 12th grade or earlier if lost to follow-up), whichever came first.
Regional air pollution exposures were defined as the community-level annual mean concentrations in the baseline year for each cohort (ie, 1993, 1996, and 2006). Data were not available for 1993 on PM10 in 4 communities (Alpine, Lake Gregory, Riverside, and Upland) and PM2.5 in any community, therefore, 1994 concentrations were used. Point estimates were scaled to the median change in community-level annual mean concentration among the 9 communities from 1993 to 2006. These models were designed to make inferences regarding within-community changes in regional air pollution and asthma incidence rates. Incidence rate differences were calculated to provide context on absolute change; these models used sampling weights for communities to make results interpretable for the entire sample. Additional details on modeling approach are reported in eMethods in the Supplement.
Potential confounders were identified a priori based on a directed acyclic graph.16 These were baseline age (continuous), sex (female, male), ethnicity (Hispanic, non-Hispanic), race (Asian/Pacific Islander, black, Native American/other, white, mixed), presence of gas stove in home (yes, no), physical activity defined here as team sports participation (yes, no), temperature defined here as community-level mean temperature for cohort baseline year (continuous), and exposure to local near-roadway pollution (continuous). To avoid loss of sample size, missing indicators were included as needed for any categorical adjustment variable. Three sets of models were fitted for each pollutant: adjusted only for community (fixed effect), additionally adjusted for all potential confounders except local near-roadway pollution, and additionally adjusted for local near-roadway pollution. Adjustment for local pollutants was conducted separately in the third model because 198 children whose residential addresses could not be geocoded were missing these data, decreasing the sample size.
We assessed heterogeneity of the regional air pollution associations by comparing nested models using a partial likelihood ratio test with and without interaction terms for the following potential effect modifiers: sex, ethnicity, race, exposure to smoking in utero, secondhand smoke exposure, parental education, parental history of asthma, Native American ancestry (only among Hispanic children), and designation of high vs low air pollution community—based on whether a community was above or below corresponding median annual mean concentration in 1993. For evaluation of effect modification, Native American ancestry among Hispanic participants was categorized into 2 groups based on having less than or greater than 50% Native American ancestry.
Robustness of the main study findings were tested with the following sensitivity analyses: (1) excluding 1 community at a time, (2) excluding participants who reported wheeze and/or 3 or more months of cough in the prior 12 months at baseline, (3) excluding data from the first year of follow-up, (4) excluding 2006-2014 cohort participants whose baseline asthma status was defined based on the year 3 rather than year 4 questionnaire, (5) reincluding participants with missing baseline asthma status, (6) imputing asthma diagnosis date to 6 months after completion date of prior questionnaire, (7) restricting to participants with longer follow-up (followed to year 5 or later, or to year 7 or later), (8) including additional potential covariates, (9) omitting the random effect for cohort nested within community and instead bootstrapping17 at the community level to assess modeling assumptions, and (10) including a fixed effect for cohort to adjust for potential temporal confounding.
Two pollutant models were fitted whenever the correlations between covariates were found to be sufficiently low to avoid multicollinearity. In addition, a set of sensitivity analyses were conducted using Cox proportional hazards regression, using the same modeling approach as the main fully adjusted model but with no random effect. These models were used to evaluate (1) the inclusion of time-varying calendar year to adjust for potential temporal confounding and (2) the use of a time-varying air pollution exposure variable. No apparent violation of the underlying assumption of proportional hazards was detected based on inclusion of a time-dependent covariate for air pollution. Due to missing air pollution data in earlier years (as noted here) and no air pollution data after 2011 as well as missing PM2.5 data for 1 community in 2005, air pollution for these years was imputed by extending the closest years’ air pollution data (ie, 1994 for 1993, 2006 for 2005, and 2011 for 2012 and later years).
All hypotheses were tested assuming a .05 significance level and a 2-sided alternative hypothesis. P values were not adjusted for multiple comparisons because the tests were hypothesis driven. All analyses were conducted using SAS software version 9.4 (SAS Institute).
The characteristics of the 4140 children included in this study are described overall and by cohort in Table 1. Because the outcome was asthma incidence during follow-up, we excluded participants (in a hierarchical manner) who had no follow-up questionnaire (n = 1503), had physician-diagnosed asthma at baseline (n = 804), or were missing baseline asthma status (n = 143). For 2006-2014 cohort participants, data from the fourth year, according to the original 2006-2014 cohort timeline, were used to define baseline asthma status. If no questionnaire was completed that year, data from the prior year were used (n = 467). If data from neither the fourth-year nor third-year questionnaire were available, those participants were considered as having missing baseline asthma status and excluded from the analysis (n = 268). The final study population comprised 4140 children, including 1093, 1170, and 1877 from the 1993-2001, 1996-2004, and 2006-2014 cohorts, respectively. Descriptive statistics on participants excluded from the analysis are given in eTable 1 in the Supplement. Mean person-years observed per child were similar across the cohorts: 5.7, 5.8, and 6.0 for the 1993-2001, 1996-2004, and 2006-2014 cohorts, respectively. The crude incidence rate for asthma was the highest for the 1996-2004 cohort (2.69 cases per 100 person-years) and lowest for the 2006-2014 cohort (1.80 cases per 100 person-years).
Regional air pollution concentrations generally decreased among the 9 communities over the course of the study period (Figure 1). The median changes in community-level annual mean concentration among the 9 communities from 1993 to 2006 were −8.9 parts per billion (ppb) (range, −21.4 to 4.8) for ozone, −4.3 ppb (range, −14.1 to −0.8) for nitrogen dioxide, −4.0 μg/m3 (range, −10.9 to 4.3) for PM10, and −8.1 μg/m3 (range, −15.2 to 0.7) for PM2.5. Reductions in air pollution were larger in communities with higher 1993 concentrations.
Plots of the unadjusted community-level data, along with community-specific regression lines, comparing asthma incidence rates with regional air pollution concentrations across the 3 cohorts are shown in Figure 2. Greater reductions in asthma incidence rates were observed in communities with larger decreases in either nitrogen dioxide or PM2.5 concentrations. Results were less consistent for ozone and PM10.
Reductions in regional nitrogen dioxide and PM2.5, but not ozone or PM10, levels were statistically significantly associated with reductions in asthma incidence rate among children (Table 2). For community-level nitrogen dioxide, the IRR for asthma was 0.83 (95% CI, 0.74-0.92) when only adjusted for community (for a reduction of 4.3 ppb). When adjusted for additional potential confounders, including near-roadway pollution, the IRR was 0.81 (95% CI, 0.72-0.91). For community-level PM2.5, the IRR was 0.82 (95% CI, 0.69-0.98) (for a reduction of 8.1 μg/m3). In the adjusted model with near-roadway pollution, the IRR was 0.82 (95% CI, 0.67-0.99). Reduced risks associated with decreasing ozone and PM10 did not reach statistical significance. Results for incidence rate differences showed absolute decreases of 0.83, 1.53, 0.78, and 0.46 cases per 100 person-years for nitrogen dioxide, PM2.5, ozone, and PM10, respectively.
Associations for nitrogen dioxide and PM2.5 did not substantially vary by sex, ethnicity, race, exposure to smoking in utero, exposure to secondhand smoke, parental education, parental history of asthma, Native American ancestry (among Hispanic children), or high or low 1993 air pollution level (eTable 2 in the Supplement).
Sensitivity analyses demonstrated results for nitrogen dioxide were robust to a variety of analytical decisions as reported in eTables 3 and 4 in the Supplement. Results for analyses excluding participants who reported wheeze and/or 3 or more months of cough in the prior 12 months at baseline (IRR, 0.81 [95% CI, 0.71-0.92]), targeting both cough- and wheeze-variant potential asthmatic cases, as well as in analyses excluding the first year of follow-up (IRR, 0.77 [95% CI, 0.68-0.88]), targeting prevalent cases, remained similar to those of the main analysis. Two pollutant models were fitted for both nitrogen dioxide and PM2.5 with PM10 (Pearson correlation coefficients were 0.50 and 0.55, respectively) and for nitrogen dioxide with ozone (correlation = 0.54). Results for nitrogen dioxide remained robust. Results for PM2.5 were not statistically significant. Two pollutant models were not fitted for nitrogen dioxide with PM2.5 (correlation = 0.60), nor PM2.5 with ozone (correlation = 0.62). Results for nitrogen dioxide and PM2.5 based on models with time-varying exposure remained statistically significant and point estimates were similar, although attenuated, compared with the results for the main baseline exposure models (eTable 4 in the Supplement). Findings for nitrogen dioxide based on models adjusted for calendar time remained robust, while those for PM2.5 did not and results for this model were null (eTable 4 in the Supplement). Overall in sensitivity analyses, associations for nitrogen dioxide remained statistically significant and generally similar in magnitude to those in the primary analysis. These analyses revealed, however, that the findings for PM2.5 were not consistently robust.
Quiz Ref IDReductions in levels of regional nitrogen dioxide from 1993 to 2014 were statistically significantly associated with improvements in asthma incidence rates in Southern Californian children. These results were independent of changes in exposure to near-roadway pollution. Findings from this study also suggested a potential association with regional PM2.5; these results, however, were less robust to sensitivity analyses. Nitrogen dioxide results remained robust in sensitivity analyses. Associations did not appear to be substantially influenced by a single community. The inclusion of prevalent asthma cases in the study population could have resulted in exacerbation in addition to incidence being captured in the definition of the outcome. Sensitivity analyses excluding potential prevalent cases, including cough- and wheeze-variant potential asthmatic cases, did not show marked differences compared with main results, suggesting this source of bias was not of major concern. Findings for PM2.5 should be interpreted with caution because these appeared more sensitive to analytical choices, as demonstrated in sensitivity analyses. Point estimates for PM2.5 from sensitivity analyses were generally similar to those of the main models, but 95% CIs were wider, and several included the null.
This study provides evidence of a robust association between children’s exposure to community-level nitrogen dioxide and development of asthma in childhood. Previous review of the development of childhood asthma and environmental exposures concluded that although several ambient pollutants were associated with increased asthma, none were consistently identified.18 Several meta-analysis studies, however, have concluded that nitrogen dioxide exposure was associated with asthma incidence among children.5,19-22 A 2017 systematic review and meta-analysis reported an overall risk estimate of 1.05 (95% CI, 1.02-1.07) per 4 μg/m3 (2.1 ppb) nitrogen dioxide. The IRR from the current study scaled to a 2.1-ppb change in nitrogen dioxide was 1.11 (95% CI, 1.05-1.18; fully adjusted model including near-roadway pollution).
It is unclear whether nitrogen dioxide is the causal agent or rather is serving as a marker for the traffic-related air pollution mixture. One CHS study, using data from the 2006-2014 cohort only, reported an association between exposure to nitrogen dioxide and asthma incidence.6,23 Ambient nitrogen dioxide measured at community central sites was associated with increased asthma risk (hazard ratio, 2.18 [95% CI, 1.18-4.01]), although this association was attenuated after accounting for traffic-related pollutions (hazard ratio, 1.37 [95% CI, 0.69-2.71]). Results for nitrogen dioxide in the present study did not change when adjusted for near-roadway pollution. While positive associations were observed in both studies, the current study benefitted from a within-community, across-time design, which provided increased statistical power and allowed for the control of community-level unmeasured confounders. These results suggest that nitrogen dioxide is capturing effects of air pollution exposure beyond local near-roadway pollution. PM2.5 has also been found to be associated with asthma incidence,5,20,24 and while PM2.5 was positively associated with incident asthma, the results were not consistently robust to sensitivity analyses. PM2.5 mass in the current study comprises traffic-related PM2.5 and that from other sources, such as dust and ocean spray. Incorporating data on PM2.5 mass sources may help specify the relation between PM2.5 and incident asthma. From 1990 to 2012, a period corresponding to the period of study in the present analyses, levels of diesel PM have decreased by 68% in California25 and this downward trend was observed in Los Angeles.26 Furthermore, measures of diesel PM, such as elemental carbon and black carbon, are correlated with nitrogen dioxide particularly near major roadways.27,28 Studies have reported positive associations between measures of diesel PM and incident childhood asthma,5 but a lack of elemental carbon or black carbon data across all 3 cohorts precluded the evaluation of this relation in the present study.
There is evidence for the plausibility of a biological mechanism specifically for nitrogen dioxide. Studies indicate that at concentrations typical of high-income countries, exposure to nitrogen dioxide induces airway inflammation,29,30 airway hyperresponsiveness,31 and oxidative stress.32 In healthy adult humans, controlled exposure to nitrogen dioxide produced enhanced pulmonary neutrophilic inflammation and the promotion of a Th2 phenotype.33 The UK Committee on the Medical Effects of Air Pollutants identified 4 mechanisms for how air pollution might contribute to the pathogenesis of asthma: (1) oxidative stress and damage, including the depletion of antioxidants; (2) airway wall remodeling, leading to structural changes in the airways; (3) inflammatory pathways and immunological effects, including effects on the expression of inflammatory mediators; and (4) enhancement of respiratory sensitization to allergens.34
A benefit of the modeling framework used here was that communities were compared with themselves at 3 points in time, thus reducing the potential for confounding by spatial factors, under the assumption that contextual variables in the community did not change. A concern remains, however, for temporal confounding. Controlling for factors with trends across the 3 study cohorts (eg, health insurance, tobacco exposure, ethnicity) did not change results. Adjusting for cohort or time-varying calendar year showed the nitrogen dioxide findings to be robust while the PM2.5 findings were not. Furthermore, secular trends in asthma rates have been increasing over the study period,35 which would bias the results toward the null. In addition, reductions in asthma incidence rates were larger in communities with greater reductions in nitrogen dioxide concentrations, further indicating results were not likely to be simply an artifact of secular trends in asthma diagnosis or other potential temporal confounders.
The study has several limitations. First, baseline, rather than time-varying, community-level annual average pollutant concentration was used as the exposure in the main models. Although a model with time-varying exposure was implemented using Cox proportional hazards regression in sensitivity analyses, this was not used as the main model because it was not possible to obtain estimates with the multilevel modeling approach (eg, when a fixed effect for community and a random effect of cohort nested within community was included). The sensitivity analyses using a Cox model with no random effect and time-varying air pollutant exposure generated results that were similar, although attenuated, compared with the Poisson models using baseline air pollutant exposure. These time-varying exposure data had more missing exposure data that necessitated imputation and, as such, may have resulted in greater exposure misclassification compared with the use of baseline exposure data. Additionally, these models may not capture the best time window of exposure because only a 1-year lag could be used given that no exposure data were available prior to the start of the study.
Second, the modeling framework, which controls for spatial confounding based on a fixed effect for community, remained susceptible to temporal confounding. This source of confounding was evaluated in sensitivity analyses, in which the findings for nitrogen dioxide were robust. Third, the definition of asthma incidence depended on a questionnaire-based assessment of physician-diagnosed asthma, rather than a clinical evaluation of asthma (eg, methacholine challenge test). Studies examining the validity of questionnaire-based asthma diagnosis in children, using questions similar to those used in the current study, have reported a specificity of 96% compared with health claims as the reference standard36 and a specificity of 87% compared with a clinical assessment as the standard.37 Fourth, a lack of data on measures of diesel PM (eg, elemental carbon or black carbon) and PM2.5 mass sources precluded the investigation of temporal trends in concentrations of these air pollutants and the incidence of childhood asthma.
Among children in Southern California, decreases in ambient nitrogen dioxide and PM2.5 between 1993 and 2014 were significantly associated with lower asthma incidence. There were no statistically significant associations for ozone or PM10.
Corresponding Author: Erika Garcia, PhD, Department of Preventive Medicine, Keck School of Medicine, University of Southern California, 2001 N Soto St, SSB1 225F MC 9237, Los Angeles, CA 90089-9237 (firstname.lastname@example.org).
Accepted for Publication: April 17, 2019.
Author Contributions: Drs Garcia and Gilliland had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Garcia, Berhane, McConnell, Urman, Gilliland.
Acquisition, analysis, or interpretation of data: Garcia, Berhane, Islam, Urman, Chen, Gilliland.
Drafting of the manuscript: Garcia, Berhane, Islam.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Garcia, Berhane, Islam, Urman, Chen.
Obtained funding: Gilliland.
Administrative, technical, or material support: Gilliland.
Supervision: Berhane, Islam, Gilliland.
Conflict of Interest Disclosures: Dr McConnell reported receiving grants from the National Institutes of Health during the conduct of the study and outside the submitted work. No other disclosures were reported.
Funding/Support: This work was supported by the National Institute of Environmental Health Sciences (grants P30ES007048, P01ES009581, R01ES021801, and R01ES025786), the National Heart, Lung, and Blood Institute (grant R01HL118455), the US Environmental Protection Agency (grants R826708 and RD831861), and the Hastings Foundation.
Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Additional Contributions: We thank Ed Avol, MS, from the Department of Preventive Medicine at the University of Southern California, and Fred Lurmann, MS, from Sonoma Technology Inc, for their review of the manuscript. They did not receive financial compensation for their contributions. We also thank the participating students and their families, the school staff and administrators, and the members of the study field team for their efforts.
Create a personal account or sign in to: