Abstract

The authors investigated the relation between ambient concentrations of five of the Environmental Protection Agency's criteria pollutants and asthma exacerbations (daily symptoms and use of rescue inhalers) among 990 children in eight North American cities during the 22-month prerandomization phase (November 1993–September 1995) of the Childhood Asthma Management Program. Short-term effects of carbon monoxide, nitrogen dioxide, particulate matter less than 10 μm in aerodynamic diameter (PM10), sulfur dioxide, and warm-season ozone were examined in both one-pollutant and two-pollutant models, using lags of up to 2 days. Lags in carbon monoxide and nitrogen dioxide were positively associated with both measures of asthma exacerbation, and the 3-day moving sum of sulfur dioxide levels was marginally related to asthma symptoms. PM10 and ozone were unrelated to exacerbations. The strongest effects tended to be seen with 2-day lags, where a 1-parts-per-million change in carbon monoxide and a 20-parts-per-billion change in nitrogen dioxide were associated with symptom odds ratios of 1.08 (95% confidence interval (CI): 1.02, 1.15) and 1.09 (95% CI: 1.03, 1.15), respectively, and with rate ratios for rescue inhaler use of 1.06 (95% CI: 1.01, 1.10) and 1.05 (95% CI: 1.01, 1.09), respectively. The authors believe that the observed carbon monoxide and nitrogen dioxide associations can probably be attributed to mobile-source emissions, though more research is required.

Many recent panel studies of the relation between air pollution concentrations and asthma aggravation in children have shown positive associations. Of the Environmental Protection Agency's six criteria pollutants, particulate matter and ozone have exhibited the most evidence of associations with asthma exacerbation (112). Still, carbon monoxide, nitrogen dioxide, and sulfur dioxide have also been linked to asthma aggravation (3, 5, 8, 10, 12). Because much of the evidence for a relation between air pollution and asthma aggravation has arisen from standard panel studies, which are carried out in a single geographic region and over a short time frame (e.g., a season), individual study results may not apply to different populations and longer periods of time. An analysis of 846 children participating in the National Cooperative Inner-City Asthma Study found that during the summer months, levels of ozone, nitrogen dioxide, and sulfur dioxide were associated with morning asthma symptoms in single-pollutant models (13), but particulate matter less than 10 μm in aerodynamic diameter (PM10) was not. In contrast to the smaller panel studies, the National Cooperative Inner-City Asthma Study examined these relations in children living in eight distinct cities, so its results may be more generally applicable.

We investigated the relation between daily ambient concentrations of carbon monoxide, nitrogen dioxide, ozone, PM10, and sulfur dioxide and daily asthma symptoms and rates of rescue inhaler use among 990 children in eight North American cities during the prerandomization phase of the Childhood Asthma Management Program (CAMP) (14). This report extends previous analyses of children who participated in CAMP at the Seattle, Washington, center (8, 12). In those analyses, Yu et al. (12) and Slaughter et al. (8) reported associations between carbon monoxide, PM10, and fine particulate matter (defined as <1 μm) and the risks of any asthma symptoms and more severe asthma symptoms, respectively. For example, a 1.11-fold increase (95 percent confidence interval (CI): 1.03, 1.20) in the odds of any symptoms (12) and a 1.10-fold increase (95 percent CI: 1.01, 1.18) in the odds of more severe symptoms (8) were observed per 10-μg/m3 difference in the 1-day lag of PM10.

The present report provides a unique contribution in that it can be considered a meta-analysis of eight large, within-city panel studies; however, it does not suffer from many of the challenges associated with meta-analyses in the published literature (e.g., between-study heterogeneity and obvious publication bias) (15). The large and geographically diverse panel of children participating in CAMP during the prerandomization phase was observed from November 1993 to September 1995, and individuals were followed for 1–6 months. This allowed us to examine the health effects of ambient concentrations of carbon monoxide, nitrogen dioxide, PM10, and sulfur dioxide across seasons and geographic regions. Because the health effects of ozone are thought to occur during the warm season, we limited the analysis of ozone to the months of May through September in 1994 and 1995. Like the findings from the National Cooperative Inner-City Asthma Study, results from this study may be applicable to a broad population, but with the exception of those pertaining to ozone, they are intended to represent year-round associations (e.g., averages of season-specific associations).

MATERIALS AND METHODS

CAMP is a National Heart, Lung, and Blood Institute-sponsored clinical trial that recruited children from eight North American cities: Albuquerque, New Mexico; Baltimore, Maryland; Boston, Massachusetts; Denver, Colorado; San Diego, California; Seattle, Washington; St. Louis, Missouri; and Toronto, Ontario, Canada. The primary purpose of CAMP is to evaluate the long-term effects of daily use of inhaled antiinflammatory medication on asthma status and lung growth in children with mild-to-moderate asthma (14). From the time when diary card completion began until randomization, children used albuterol as needed for symptoms and oral prednisone for exacerbations but were restricted from daily controller medication. At the end of each day, under the supervision of their families, the children completed daily diary cards to report their experience with asthma symptoms. The children and their families were carefully instructed on the importance of accurately reporting symptoms and medication use on the daily diary cards.

Inclusion criteria

All subjects studied here were randomized into the primary CAMP study, and inclusion criteria have been reported previously (14). Participants were excluded if we were able to verify that the centroid of the ZIP or postal code in which they lived was greater than 50 miles (80 km) from the nearest Aerometric Information Retrieval System pollutant monitor. Twelve children from Albuquerque, one from Baltimore, 16 from Denver, 10 from Seattle, four from St. Louis, and five from Toronto were excluded for this reason. Of the 1,041 children randomized into CAMP, 48 were excluded because of our 50-mile rule. Three additional children were excluded because information on a key measure of sensitivity, the methacholine concentration required to reduce forced expiratory volume in 1 second by 20 percent, was not available.

Ambient air pollution

Data on pollution concentrations were obtained from the Aerometric Information Retrieval System (www.epa.gov/air/data/index.html) for US cities and from Environment Canada (www.ec.gc.ca/envhome.html) for Toronto. Within each city, multiple monitors were used to measure ambient pollution concentrations. Monitor-specific concentrations were 24-hour averages for PM10, carbon monoxide, nitrogen dioxide, and sulfur dioxide and a 1-hour maximum for ozone. Daily, citywide concentrations were calculated by averaging monitor-specific concentrations, and an adjusted average was used on days on which only a fraction of monitors were in operation. Individual monitors that appeared to be recalibrated during the study (e.g., an obvious shift in pollution concentrations) were excluded from ambient pollution exposure calculations. For carbon monoxide, nitrogen dioxide, sulfur dioxide, and ozone, monitors that were in operation for at least half of the study days and, for PM10, monitors that were in operation for at least 90 of the study days were included in the ambient pollution exposure calculations.

The National Oceanic and Atmospheric Administration provided meteorologic data for the eight participating cities. Daily temperature and humidity measures were summarized using the maxima of hourly measurements. Wind speed was a 24-hour average, precipitation a daily sum, and stagnation a percentage of hourly measurements with wind speed measured below 6 miles/hour (9.6 km/hour).

Outcomes

Children recorded the severity of their asthma symptoms on daily diary cards completed at the end of each day. Asthma severity was coded as follows: 0 = no asthma symptoms; 1 = 1–3 mild asthma episodes, each lasting 2 hours or less; 2 = four or more mild asthma episodes or one or more episodes that temporarily interfered with activity, play, school, or sleep; and 3 = one or more asthma episodes lasting longer than 2 hours or resulting in shortening of normal activity, seeing a doctor for acute care, or going to a hospital for acute care. Daily diary cards on which symptom severity was not recorded (1.6 percent) were excluded from all analyses. Severity codes 1, 2, and 3 were collapsed into a single category to create a dichotomous response variable of any daily asthma symptoms versus none. On the diary cards, children also reported their daily number of rescue inhaler puffs. Since the standard dose of inhaled medication is two puffs, we defined one “use” of a rescue inhaler as half the standard number of puffs (i.e., one puff), and if an odd number of puffs was reported, we rounded up. The dichotomous variable of any asthma symptoms versus none and the number of rescue inhaler uses were the outcomes studied.

Statistical analysis

Models.

We sought to examine the relation between asthma aggravation and daily ambient air pollution concentrations for five criteria pollutants—carbon monoxide, nitrogen dioxide, ozone, PM10, and sulfur dioxide—across all cities. Towards this end, we constructed models within each city and then combined the city-specific parameter estimates to obtain the reported study-wide estimates using meta-analytic procedures (16). Since there is no general consensus regarding the appropriate day-lag of ambient concentrations to consider when studying their relations with asthma exacerbations, we examined the effects of 0-, 1-, and 2-day lags and the 3-day moving sum in separate models. The 3-day moving sum is the sum of the 0-, 1-, and 2-day lags in concentrations, and is therefore three times the 3-day moving average. One- and two-pollutant models were considered.

The generalized estimating equations procedure (17, 18) was used to estimate the marginal relation between asthma symptoms and rescue inhaler uses and ambient pollution concentrations in each city. Logistic and Poisson regression models were constructed for asthma symptoms and rescue inhaler uses, respectively. We report increases in the odds of asthma symptoms and in the rates of rescue inhaler use associated with large differences in exposure concentrations: carbon monoxide, 1.0 parts per million (ppm); nitrogen dioxide, 20 parts per billion (ppb); ozone, 30 ppb; PM10, 25 μg/m3; and sulfur dioxide, 10 ppb. For estimation, we used a working independence covariance structure, since multiple lags in exposure were potentially associated with both responses. In this setting, nonindependence covariance weighting can result in biased estimates (19, 20). Robust standard errors were used to account for the intraindividual correlations.

When constructing within-city regression models, the potential for confounding by a number of factors was considered carefully. We controlled for day of the week, ethnicity (White, Black, Hispanic/Latin, and other), and annual family income (<$15,000, $15,000–<$30,000, $30,000–<$50,000, ≥$50,000, and “no answer”), as well as flexible functions of age and log-transformed sensitivity to the methacholine challenge using natural splines with knots fixed at the 25th, 50th, and 75th percentiles (quartiles 1, 2, and 3) of the respective distributions. Confounding due to seasonal factors is a primary consideration in air pollution epidemiology, and we adjusted for “season” in three ways. First, we controlled for same-day maximum hourly temperature, humidity, and their interaction by using natural splines with knots placed at quartiles 1, 2, and 3 of each distribution. Second, we adjusted for calendar date with natural splines, fixing knots at 45-day separations, which allowed the effect of latent, seasonal factors to change smoothly over the study period. Finally, we decomposed daily ambient pollution concentrations into between- and within-subject exposures and modeled them separately; we report estimates of within-subject exposure effects. Specifically, subject i's exposure to pollutant x1 on day t was decomposed as follows:
\[x_{1,it}{=}{\bar{x}}_{1,i}{+}(x_{1,it}{-}{\bar{x}}_{1,i}).\]
The between-subject component,
\({\bar{x}}_{1,i},\)
is subject i's average exposure during the time in which he or she was observed, and the within-subject component,
\(x_{1,it}{-}{\bar{x}}_{1,i},\)
is the deviation from the average on day t. Within-subject effects are likely to be far less confounded by season than between-subject effects (21), and they are more comparable to the parameters estimated in a traditional panel study, where all subjects are observed over a single time period. For further discussion regarding the decomposition of time-varying exposures into between- and within-subject components, see the papers by Sheppard (21) and Neuhaus and Kalbfleisch (22).

Daily ambient air pollution records were nearly complete for most pollutants, although data on sulfur dioxide in Albuquerque and data on nitrogen dioxide in Seattle were not collected and therefore were not included in our analyses. PM10 was measured on a daily basis in Seattle (97 percent of study days) and Albuquerque (99 percent of study days). In all other cities, it was measured regularly but only on a subset of days: 50 percent in Baltimore, 23 percent in Boston, 37 percent in Denver, 24 percent in San Diego, 19 percent in St. Louis, and 47 percent in Toronto. To recover a fraction of the missing PM10 data, we constructed 10 imputation data sets per city for analyses involving PM10. Meteorologic values (wind speed, stagnation, temperature, humidity), other pollutant concentrations, including 0- to 4-day lags, and daily summaries of asthma aggravation (e.g., the fraction of children experiencing asthma symptoms and the mean number of rescue inhaler uses) were used to impute missing PM10 values. Depending upon the city, the prediction models accounted for 50–75 percent of the variation in PM10 on days on which PM10 data were available. Within each city, regression analysis was done on each imputation data set, and results were combined to ascertain city-specific estimates using standard multiple imputation techniques (23).

Analysis approach.

Within each city, one- and two-pollutant models were constructed with ambient concentrations decomposed into between- and within-subject components. One-pollutant models were given by
\[g(\mathrm{{\mu}}_{it}){=}\mathrm{{\beta}}_{0}{+}\mathrm{{\beta}}_{w}(x_{1,it}{-}{\bar{x}}_{1,i}){+}\mathrm{{\beta}}_{b}{\bar{x}}_{1,i}{+}{\ldots},\]
where μit is the expected response for subject i on day t, g is the link function, and interest is in the parameter βw. Similarly, two-pollutant models were given by
\begin{eqnarray*}&&g(\mathrm{{\mu}}_{it}){=}\mathrm{{\gamma}}_{0}{+}\mathrm{{\gamma}}_{1,w}(x_{1,it}{-}{\bar{x}}_{1,i}){+}\mathrm{{\gamma}}_{1,b}{\bar{x}}_{1,i}\\&&{+}\mathrm{{\gamma}}_{2,w}(x_{2,it}{-}{\bar{x}}_{2,i}){+}\mathrm{{\gamma}}_{2,b}{\bar{x}}_{2,i}{+}{\ldots},\end{eqnarray*}
and to capture the effect of simultaneous shifts of within-subject concentrations in two pollutants, interest is in the sum, γ1,w + γ2,w. When PM10 was modeled, multiple imputation was necessary to capture city-specific effects, and when ozone was modeled, we limited the analysis period to the warm months (May through September). Once city-specific estimates had been obtained, meta-analytic techniques were used to combine them.

We used SAS, version 8.02 (24), to manage all data and to build the multiple-imputation data sets. All regression modeling and the summarization of estimates was done with the R programming language (25), using the geepack (26), mitools (27), and rmeta (28) packages.

RESULTS

All subjects considered in this analysis were randomized into CAMP. A total of 990 children were studied (table 1). They were followed for a median of nearly 2 months and were 5–12 years old when diary card completion began (interquartile range, 7–10). Children reported asthma symptoms and rescue inhaler uses on 52 percent and 44 percent of all person-days, respectively, though in each city there were subjects who experienced symptoms or used the rescue inhaler on every day under study or on no days.

TABLE 1.

Characteristics of participants in the prerandomization phase of the Childhood Asthma Management Program, November 1993–September 1995




Albuquerque, New Mexico

Baltimore, Maryland

Boston, Massachusetts

Denver, Colorado

San Diego, California

Seattle, Washington

St. Louis, Missouri

Toronto, Ontario, Canada

All cities
No. of children109127123128122134128119990
No. of study days per child
    Median575651395754685955
    IQR*46–7448–7043–6433–4449–7747–6755–7948–7444–70
    Range31–20128–12529–9621–11530–11528–11228–11830–12121–201
Age (years)
    Median8.98.19.28.79.38.79.07.58.7
    IQR7.5–10.66.9–10.07.8–10.56.5–10.37.3–11.26.9–10.27.3–10.76.6–9.57.0–10.5
    Range5.1–12.75.0–12.85.0–12.85.0–12.95.0–<13.05.0–<13.05.0–<13.05.0–12.65.0–<13.0
Sensitivity to the methacholine challenge (mg/ml)
    Median1.0470.9051.0941.7510.9240.6911.5701.1391.095
    IQR0.326–2.8660.499–2.1400.518–3.0021.007–4.3770.444–1.8940.322–1.7270.657–4.1640.617–2.6350.501–2.703
    Range0.023–10.1530.023–11.5480.058–10.0130.176–11.6380.108–9.6630.057–11.9650.08–13.2650.133–10.5840.023–13.265
Symptom score,§
    00.510.580.520.540.450.400.390.460.48
    10.420.380.420.410.450.510.520.460.45
    20.060.040.060.040.100.080.080.070.07
    30.010.010.010.000.010.010.010.010.01
Subject-specific proportion of days with symptom score >0
    Median0.460.370.440.430.530.570.600.510.48
    IQR0.32–0.610.25–0.610.29–0.620.26–0.640.34–0.790.4–0.870.42–0.850.38–0.730.32–0.72
    Range0.07–1.000.00–1.000.00–1.000.00–1.000.04–1.000.03–1.000.15–1.000.04–1.000.00–1.00
No. of rescue inhaler puffs per day
    00.560.600.700.570.540.440.470.600.56
    10.370.360.250.380.370.450.410.330.37
    20.050.030.040.050.080.100.100.060.07
    >20.010.000.010.000.010.010.020.010.01
Subject-specific proportion of days with ≥1 rescue inhaler puff
    Median0.400.340.230.400.440.550.480.360.40
    IQR0.24–0.590.20–0.540.09–0.420.19–0.650.25–0.670.33–0.810.33–0.780.20–0.560.22–0.65
    Range
0.00–0.98
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.02–1.00
0.00–1.00
0.00–1.00



Albuquerque, New Mexico

Baltimore, Maryland

Boston, Massachusetts

Denver, Colorado

San Diego, California

Seattle, Washington

St. Louis, Missouri

Toronto, Ontario, Canada

All cities
No. of children109127123128122134128119990
No. of study days per child
    Median575651395754685955
    IQR*46–7448–7043–6433–4449–7747–6755–7948–7444–70
    Range31–20128–12529–9621–11530–11528–11228–11830–12121–201
Age (years)
    Median8.98.19.28.79.38.79.07.58.7
    IQR7.5–10.66.9–10.07.8–10.56.5–10.37.3–11.26.9–10.27.3–10.76.6–9.57.0–10.5
    Range5.1–12.75.0–12.85.0–12.85.0–12.95.0–<13.05.0–<13.05.0–<13.05.0–12.65.0–<13.0
Sensitivity to the methacholine challenge (mg/ml)
    Median1.0470.9051.0941.7510.9240.6911.5701.1391.095
    IQR0.326–2.8660.499–2.1400.518–3.0021.007–4.3770.444–1.8940.322–1.7270.657–4.1640.617–2.6350.501–2.703
    Range0.023–10.1530.023–11.5480.058–10.0130.176–11.6380.108–9.6630.057–11.9650.08–13.2650.133–10.5840.023–13.265
Symptom score,§
    00.510.580.520.540.450.400.390.460.48
    10.420.380.420.410.450.510.520.460.45
    20.060.040.060.040.100.080.080.070.07
    30.010.010.010.000.010.010.010.010.01
Subject-specific proportion of days with symptom score >0
    Median0.460.370.440.430.530.570.600.510.48
    IQR0.32–0.610.25–0.610.29–0.620.26–0.640.34–0.790.4–0.870.42–0.850.38–0.730.32–0.72
    Range0.07–1.000.00–1.000.00–1.000.00–1.000.04–1.000.03–1.000.15–1.000.04–1.000.00–1.00
No. of rescue inhaler puffs per day
    00.560.600.700.570.540.440.470.600.56
    10.370.360.250.380.370.450.410.330.37
    20.050.030.040.050.080.100.100.060.07
    >20.010.000.010.000.010.010.020.010.01
Subject-specific proportion of days with ≥1 rescue inhaler puff
    Median0.400.340.230.400.440.550.480.360.40
    IQR0.24–0.590.20–0.540.09–0.420.19–0.650.25–0.670.33–0.810.33–0.780.20–0.560.22–0.65
    Range
0.00–0.98
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.02–1.00
0.00–1.00
0.00–1.00
*

IQR, interquartile range.

The dose of methacholine required to decrease forced expiratory volume in 1 second by 20%.

Symptom score and rescue inhaler uses are given as the proportion of daily diary cards completed.

§

Symptom scores: 0 = no asthma symptoms; 1 = 1–3 mild asthma episodes, each lasting 2 hours or less; 2 = four or more mild asthma episodes or one or more episodes that temporarily interfered with activity, play, school, or sleep; and 3 = one or more asthma episodes lasting longer than 2 hours or resulting in shortening of normal activity, seeing a doctor for acute care, or going to a hospital for acute care.

In these calculations, a summary measure is represented for each subject (e.g., fraction of days with asthma symptom severity score X), and percentiles of the distribution of summary measures are reported.

TABLE 1.

Characteristics of participants in the prerandomization phase of the Childhood Asthma Management Program, November 1993–September 1995




Albuquerque, New Mexico

Baltimore, Maryland

Boston, Massachusetts

Denver, Colorado

San Diego, California

Seattle, Washington

St. Louis, Missouri

Toronto, Ontario, Canada

All cities
No. of children109127123128122134128119990
No. of study days per child
    Median575651395754685955
    IQR*46–7448–7043–6433–4449–7747–6755–7948–7444–70
    Range31–20128–12529–9621–11530–11528–11228–11830–12121–201
Age (years)
    Median8.98.19.28.79.38.79.07.58.7
    IQR7.5–10.66.9–10.07.8–10.56.5–10.37.3–11.26.9–10.27.3–10.76.6–9.57.0–10.5
    Range5.1–12.75.0–12.85.0–12.85.0–12.95.0–<13.05.0–<13.05.0–<13.05.0–12.65.0–<13.0
Sensitivity to the methacholine challenge (mg/ml)
    Median1.0470.9051.0941.7510.9240.6911.5701.1391.095
    IQR0.326–2.8660.499–2.1400.518–3.0021.007–4.3770.444–1.8940.322–1.7270.657–4.1640.617–2.6350.501–2.703
    Range0.023–10.1530.023–11.5480.058–10.0130.176–11.6380.108–9.6630.057–11.9650.08–13.2650.133–10.5840.023–13.265
Symptom score,§
    00.510.580.520.540.450.400.390.460.48
    10.420.380.420.410.450.510.520.460.45
    20.060.040.060.040.100.080.080.070.07
    30.010.010.010.000.010.010.010.010.01
Subject-specific proportion of days with symptom score >0
    Median0.460.370.440.430.530.570.600.510.48
    IQR0.32–0.610.25–0.610.29–0.620.26–0.640.34–0.790.4–0.870.42–0.850.38–0.730.32–0.72
    Range0.07–1.000.00–1.000.00–1.000.00–1.000.04–1.000.03–1.000.15–1.000.04–1.000.00–1.00
No. of rescue inhaler puffs per day
    00.560.600.700.570.540.440.470.600.56
    10.370.360.250.380.370.450.410.330.37
    20.050.030.040.050.080.100.100.060.07
    >20.010.000.010.000.010.010.020.010.01
Subject-specific proportion of days with ≥1 rescue inhaler puff
    Median0.400.340.230.400.440.550.480.360.40
    IQR0.24–0.590.20–0.540.09–0.420.19–0.650.25–0.670.33–0.810.33–0.780.20–0.560.22–0.65
    Range
0.00–0.98
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.02–1.00
0.00–1.00
0.00–1.00



Albuquerque, New Mexico

Baltimore, Maryland

Boston, Massachusetts

Denver, Colorado

San Diego, California

Seattle, Washington

St. Louis, Missouri

Toronto, Ontario, Canada

All cities
No. of children109127123128122134128119990
No. of study days per child
    Median575651395754685955
    IQR*46–7448–7043–6433–4449–7747–6755–7948–7444–70
    Range31–20128–12529–9621–11530–11528–11228–11830–12121–201
Age (years)
    Median8.98.19.28.79.38.79.07.58.7
    IQR7.5–10.66.9–10.07.8–10.56.5–10.37.3–11.26.9–10.27.3–10.76.6–9.57.0–10.5
    Range5.1–12.75.0–12.85.0–12.85.0–12.95.0–<13.05.0–<13.05.0–<13.05.0–12.65.0–<13.0
Sensitivity to the methacholine challenge (mg/ml)
    Median1.0470.9051.0941.7510.9240.6911.5701.1391.095
    IQR0.326–2.8660.499–2.1400.518–3.0021.007–4.3770.444–1.8940.322–1.7270.657–4.1640.617–2.6350.501–2.703
    Range0.023–10.1530.023–11.5480.058–10.0130.176–11.6380.108–9.6630.057–11.9650.08–13.2650.133–10.5840.023–13.265
Symptom score,§
    00.510.580.520.540.450.400.390.460.48
    10.420.380.420.410.450.510.520.460.45
    20.060.040.060.040.100.080.080.070.07
    30.010.010.010.000.010.010.010.010.01
Subject-specific proportion of days with symptom score >0
    Median0.460.370.440.430.530.570.600.510.48
    IQR0.32–0.610.25–0.610.29–0.620.26–0.640.34–0.790.4–0.870.42–0.850.38–0.730.32–0.72
    Range0.07–1.000.00–1.000.00–1.000.00–1.000.04–1.000.03–1.000.15–1.000.04–1.000.00–1.00
No. of rescue inhaler puffs per day
    00.560.600.700.570.540.440.470.600.56
    10.370.360.250.380.370.450.410.330.37
    20.050.030.040.050.080.100.100.060.07
    >20.010.000.010.000.010.010.020.010.01
Subject-specific proportion of days with ≥1 rescue inhaler puff
    Median0.400.340.230.400.440.550.480.360.40
    IQR0.24–0.590.20–0.540.09–0.420.19–0.650.25–0.670.33–0.810.33–0.780.20–0.560.22–0.65
    Range
0.00–0.98
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.00–1.00
0.02–1.00
0.00–1.00
0.00–1.00
*

IQR, interquartile range.

The dose of methacholine required to decrease forced expiratory volume in 1 second by 20%.

Symptom score and rescue inhaler uses are given as the proportion of daily diary cards completed.

§

Symptom scores: 0 = no asthma symptoms; 1 = 1–3 mild asthma episodes, each lasting 2 hours or less; 2 = four or more mild asthma episodes or one or more episodes that temporarily interfered with activity, play, school, or sleep; and 3 = one or more asthma episodes lasting longer than 2 hours or resulting in shortening of normal activity, seeing a doctor for acute care, or going to a hospital for acute care.

In these calculations, a summary measure is represented for each subject (e.g., fraction of days with asthma symptom severity score X), and percentiles of the distribution of summary measures are reported.

Pollution concentrations during the study period are summarized in table 2. We report the number of monitors used in the calculations of daily ambient concentrations, percentiles of daily concentrations, and season-adjusted (partial) correlations among pollutants. Since PM10 data were recorded on a subset of study days in all cities except Seattle and Albuquerque, we summarized the distribution of PM10 values in each city with a single, complete imputation sample. All pollutants tended to be positively correlated with one another (note, however, that correlations with ozone correspond only to the months of May through September in 1994 and 1995). Carbon monoxide and nitrogen dioxide tended to be most highly correlated, with values ranging from 0.63 in Toronto to 0.92 in San Diego. Ozone tended to be less positively correlated with the other gases, though in several cities it was highly correlated with PM10.

TABLE 2.

Ambient air pollution concentrations and season-adjusted correlations in pollutants during the prerandomization phase of the Childhood Asthma Management Program, November 1993–September 1995*


City

No. of monitors

Percentile

Season-adjusted correlation§
10th
25th
50th
75th
90th
Nitrogen dioxide
Ozone
PM10
Sulfur dioxide
Albuquerque, New Mexico (November 25, 1993–August 29, 1995)
    Carbon monoxide (ppm × 10)34.36.08.912.618.20.760.090.31
    Nitrogen dioxide (ppb)19.812.817.825.134.310.040.26
    Ozone (ppb)545.150.055.061.268.410.03
    PM10 (μg/m3)211.515.520.627.033.51
    Sulfur dioxide (ppb)0
Baltimore, Maryland (November 9, 1993–August 13, 1995)
    Carbon monoxide (ppm × 10)64.15.67.39.914.90.690.210.510.50
    Nitrogen dioxide (ppb)216.520.825.831.536.910.440.620.49
    Ozone (ppb)1341.452.165.882.794.710.730.30
    PM10 (μg/m3)1212.417.626.435.444.010.31
    Sulfur dioxide (ppb)33.24.76.79.814.2
Boston, Massachusetts (January 6, 1994–September 12, 1995)
    Carbon monoxide (ppm × 10)46.48.110.012.414.60.800.200.280.67
    Nitrogen dioxide (ppb)416.620.124.830.635.710.470.480.68
    Ozone (ppb)433.038.252.266.082.110.640.21
    PM10 (μg/m3)59.814.520.426.232.510.36
    Sulfur dioxide (ppb)52.73.75.89.114.1
Denver, Colorado (November 3, 1993–June 15, 1995)
    Carbon monoxide (ppm × 10)75.16.38.112.518.80.85−0.280.660.53
    Nitrogen dioxide (ppb)414.919.023.328.634.610.240.640.56
    Ozone (ppb)840.251.760.567.973.010.490.22
    PM10 (μg/m3)76.812.520.127.937.310.51
    Sulfur dioxide (ppb)21.22.54.46.79.5
San Diego, California (November 24, 1993–July 2, 1995)
    Carbon monoxide (ppm × 10)85.26.38.813.618.90.920.340.520.19
    Nitrogen dioxide (ppb)811.013.919.225.532.510.390.550.23
    Ozone (ppb)943.750.859.369.975.210.44−0.03
    PM10 (μg/m3)713.722.432.442.753.910.32
    Sulfur dioxide (ppb)31.21.72.23.14.4
Seattle, Washington (November 3, 1993–August 15, 1995)
    Carbon monoxide (ppm × 10)810.112.314.919.024.00.340.880.62
    Nitrogen dioxide (ppb)01
    Ozone (ppb)426.734.043.053.861.510.570.43
    PM10 (μg/m3)411.415.222.030.244.010.65
    Sulfur dioxide (ppb)43.74.76.07.59.5
St. Louis, Missouri (November 12, 1993–August 21, 1995)
    Carbon monoxide (ppm × 10)64.65.46.37.79.50.710.440.240.34
    Nitrogen dioxide (ppb)811.514.518.322.426.710.420.450.58
    Ozone (ppb)1338.548.359.373.486.410.690.33
    PM10 (μg/m3)614.021.328.536.545.010.38
    Sulfur dioxide (ppb)93.95.37.410.713.6
Toronto, Ontario, Canada (December 10, 1993–July 25, 1995)
    Carbon monoxide (ppm × 10)45.78.09.711.312.70.630.160.370.55
    Nitrogen dioxide (ppb)417.521.326.031.537.010.400.640.63
    Ozone (ppb)427.533.643.556.567.110.700.44
    PM10 (μg/m3)16.812.017.727.237.810.48
    Sulfur dioxide (ppb)
4
0.2
1.0
2.5
4.8
8.8





City

No. of monitors

Percentile

Season-adjusted correlation§
10th
25th
50th
75th
90th
Nitrogen dioxide
Ozone
PM10
Sulfur dioxide
Albuquerque, New Mexico (November 25, 1993–August 29, 1995)
    Carbon monoxide (ppm × 10)34.36.08.912.618.20.760.090.31
    Nitrogen dioxide (ppb)19.812.817.825.134.310.040.26
    Ozone (ppb)545.150.055.061.268.410.03
    PM10 (μg/m3)211.515.520.627.033.51
    Sulfur dioxide (ppb)0
Baltimore, Maryland (November 9, 1993–August 13, 1995)
    Carbon monoxide (ppm × 10)64.15.67.39.914.90.690.210.510.50
    Nitrogen dioxide (ppb)216.520.825.831.536.910.440.620.49
    Ozone (ppb)1341.452.165.882.794.710.730.30
    PM10 (μg/m3)1212.417.626.435.444.010.31
    Sulfur dioxide (ppb)33.24.76.79.814.2
Boston, Massachusetts (January 6, 1994–September 12, 1995)
    Carbon monoxide (ppm × 10)46.48.110.012.414.60.800.200.280.67
    Nitrogen dioxide (ppb)416.620.124.830.635.710.470.480.68
    Ozone (ppb)433.038.252.266.082.110.640.21
    PM10 (μg/m3)59.814.520.426.232.510.36
    Sulfur dioxide (ppb)52.73.75.89.114.1
Denver, Colorado (November 3, 1993–June 15, 1995)
    Carbon monoxide (ppm × 10)75.16.38.112.518.80.85−0.280.660.53
    Nitrogen dioxide (ppb)414.919.023.328.634.610.240.640.56
    Ozone (ppb)840.251.760.567.973.010.490.22
    PM10 (μg/m3)76.812.520.127.937.310.51
    Sulfur dioxide (ppb)21.22.54.46.79.5
San Diego, California (November 24, 1993–July 2, 1995)
    Carbon monoxide (ppm × 10)85.26.38.813.618.90.920.340.520.19
    Nitrogen dioxide (ppb)811.013.919.225.532.510.390.550.23
    Ozone (ppb)943.750.859.369.975.210.44−0.03
    PM10 (μg/m3)713.722.432.442.753.910.32
    Sulfur dioxide (ppb)31.21.72.23.14.4
Seattle, Washington (November 3, 1993–August 15, 1995)
    Carbon monoxide (ppm × 10)810.112.314.919.024.00.340.880.62
    Nitrogen dioxide (ppb)01
    Ozone (ppb)426.734.043.053.861.510.570.43
    PM10 (μg/m3)411.415.222.030.244.010.65
    Sulfur dioxide (ppb)43.74.76.07.59.5
St. Louis, Missouri (November 12, 1993–August 21, 1995)
    Carbon monoxide (ppm × 10)64.65.46.37.79.50.710.440.240.34
    Nitrogen dioxide (ppb)811.514.518.322.426.710.420.450.58
    Ozone (ppb)1338.548.359.373.486.410.690.33
    PM10 (μg/m3)614.021.328.536.545.010.38
    Sulfur dioxide (ppb)93.95.37.410.713.6
Toronto, Ontario, Canada (December 10, 1993–July 25, 1995)
    Carbon monoxide (ppm × 10)45.78.09.711.312.70.630.160.370.55
    Nitrogen dioxide (ppb)417.521.326.031.537.010.400.640.63
    Ozone (ppb)427.533.643.556.567.110.700.44
    PM10 (μg/m3)16.812.017.727.237.810.48
    Sulfur dioxide (ppb)
4
0.2
1.0
2.5
4.8
8.8




*

Data on concentrations of particulate matter less than 10 μm in aerodynamic diameter (PM10) were available for a subset of study days in Baltimore, Boston, Denver, San Diego, St. Louis, and Toronto. Percentiles and correlations with other pollutants were based on an imputation sample used in the primary analysis for PM10. Ozone values and correlations with other pollutants were limited to the warm season, which lasted from May through September in 1994 and 1995.

To be included in the calculation of daily ambient concentrations, monitors must have been functioning for at least half of the study days. Since PM10 was not measured daily in all cities, monitors were included in the calculation of daily ambient concentrations if they were functioning for at least 90 of the study days.

Daily, monitor-specific values were 24-hour averages for carbon monoxide, nitrogen dioxide, and sulfur dioxide and a 1-hour maximum for ozone. Daily, citywide ambient concentrations for all pollutants were calculated with an adjusted average over monitor-specific values.

§

Individual seasons used in the adjustments are given by consecutive 45-day time windows.

ppm, parts per million; ppb, parts per billion.

TABLE 2.

Ambient air pollution concentrations and season-adjusted correlations in pollutants during the prerandomization phase of the Childhood Asthma Management Program, November 1993–September 1995*


City

No. of monitors

Percentile

Season-adjusted correlation§
10th
25th
50th
75th
90th
Nitrogen dioxide
Ozone
PM10
Sulfur dioxide
Albuquerque, New Mexico (November 25, 1993–August 29, 1995)
    Carbon monoxide (ppm × 10)34.36.08.912.618.20.760.090.31
    Nitrogen dioxide (ppb)19.812.817.825.134.310.040.26
    Ozone (ppb)545.150.055.061.268.410.03
    PM10 (μg/m3)211.515.520.627.033.51
    Sulfur dioxide (ppb)0
Baltimore, Maryland (November 9, 1993–August 13, 1995)
    Carbon monoxide (ppm × 10)64.15.67.39.914.90.690.210.510.50
    Nitrogen dioxide (ppb)216.520.825.831.536.910.440.620.49
    Ozone (ppb)1341.452.165.882.794.710.730.30
    PM10 (μg/m3)1212.417.626.435.444.010.31
    Sulfur dioxide (ppb)33.24.76.79.814.2
Boston, Massachusetts (January 6, 1994–September 12, 1995)
    Carbon monoxide (ppm × 10)46.48.110.012.414.60.800.200.280.67
    Nitrogen dioxide (ppb)416.620.124.830.635.710.470.480.68
    Ozone (ppb)433.038.252.266.082.110.640.21
    PM10 (μg/m3)59.814.520.426.232.510.36
    Sulfur dioxide (ppb)52.73.75.89.114.1
Denver, Colorado (November 3, 1993–June 15, 1995)
    Carbon monoxide (ppm × 10)75.16.38.112.518.80.85−0.280.660.53
    Nitrogen dioxide (ppb)414.919.023.328.634.610.240.640.56
    Ozone (ppb)840.251.760.567.973.010.490.22
    PM10 (μg/m3)76.812.520.127.937.310.51
    Sulfur dioxide (ppb)21.22.54.46.79.5
San Diego, California (November 24, 1993–July 2, 1995)
    Carbon monoxide (ppm × 10)85.26.38.813.618.90.920.340.520.19
    Nitrogen dioxide (ppb)811.013.919.225.532.510.390.550.23
    Ozone (ppb)943.750.859.369.975.210.44−0.03
    PM10 (μg/m3)713.722.432.442.753.910.32
    Sulfur dioxide (ppb)31.21.72.23.14.4
Seattle, Washington (November 3, 1993–August 15, 1995)
    Carbon monoxide (ppm × 10)810.112.314.919.024.00.340.880.62
    Nitrogen dioxide (ppb)01
    Ozone (ppb)426.734.043.053.861.510.570.43
    PM10 (μg/m3)411.415.222.030.244.010.65
    Sulfur dioxide (ppb)43.74.76.07.59.5
St. Louis, Missouri (November 12, 1993–August 21, 1995)
    Carbon monoxide (ppm × 10)64.65.46.37.79.50.710.440.240.34
    Nitrogen dioxide (ppb)811.514.518.322.426.710.420.450.58
    Ozone (ppb)1338.548.359.373.486.410.690.33
    PM10 (μg/m3)614.021.328.536.545.010.38
    Sulfur dioxide (ppb)93.95.37.410.713.6
Toronto, Ontario, Canada (December 10, 1993–July 25, 1995)
    Carbon monoxide (ppm × 10)45.78.09.711.312.70.630.160.370.55
    Nitrogen dioxide (ppb)417.521.326.031.537.010.400.640.63
    Ozone (ppb)427.533.643.556.567.110.700.44
    PM10 (μg/m3)16.812.017.727.237.810.48
    Sulfur dioxide (ppb)
4
0.2
1.0
2.5
4.8
8.8





City

No. of monitors

Percentile

Season-adjusted correlation§
10th
25th
50th
75th
90th
Nitrogen dioxide
Ozone
PM10
Sulfur dioxide
Albuquerque, New Mexico (November 25, 1993–August 29, 1995)
    Carbon monoxide (ppm × 10)34.36.08.912.618.20.760.090.31
    Nitrogen dioxide (ppb)19.812.817.825.134.310.040.26
    Ozone (ppb)545.150.055.061.268.410.03
    PM10 (μg/m3)211.515.520.627.033.51
    Sulfur dioxide (ppb)0
Baltimore, Maryland (November 9, 1993–August 13, 1995)
    Carbon monoxide (ppm × 10)64.15.67.39.914.90.690.210.510.50
    Nitrogen dioxide (ppb)216.520.825.831.536.910.440.620.49
    Ozone (ppb)1341.452.165.882.794.710.730.30
    PM10 (μg/m3)1212.417.626.435.444.010.31
    Sulfur dioxide (ppb)33.24.76.79.814.2
Boston, Massachusetts (January 6, 1994–September 12, 1995)
    Carbon monoxide (ppm × 10)46.48.110.012.414.60.800.200.280.67
    Nitrogen dioxide (ppb)416.620.124.830.635.710.470.480.68
    Ozone (ppb)433.038.252.266.082.110.640.21
    PM10 (μg/m3)59.814.520.426.232.510.36
    Sulfur dioxide (ppb)52.73.75.89.114.1
Denver, Colorado (November 3, 1993–June 15, 1995)
    Carbon monoxide (ppm × 10)75.16.38.112.518.80.85−0.280.660.53
    Nitrogen dioxide (ppb)414.919.023.328.634.610.240.640.56
    Ozone (ppb)840.251.760.567.973.010.490.22
    PM10 (μg/m3)76.812.520.127.937.310.51
    Sulfur dioxide (ppb)21.22.54.46.79.5
San Diego, California (November 24, 1993–July 2, 1995)
    Carbon monoxide (ppm × 10)85.26.38.813.618.90.920.340.520.19
    Nitrogen dioxide (ppb)811.013.919.225.532.510.390.550.23
    Ozone (ppb)943.750.859.369.975.210.44−0.03
    PM10 (μg/m3)713.722.432.442.753.910.32
    Sulfur dioxide (ppb)31.21.72.23.14.4
Seattle, Washington (November 3, 1993–August 15, 1995)
    Carbon monoxide (ppm × 10)810.112.314.919.024.00.340.880.62
    Nitrogen dioxide (ppb)01
    Ozone (ppb)426.734.043.053.861.510.570.43
    PM10 (μg/m3)411.415.222.030.244.010.65
    Sulfur dioxide (ppb)43.74.76.07.59.5
St. Louis, Missouri (November 12, 1993–August 21, 1995)
    Carbon monoxide (ppm × 10)64.65.46.37.79.50.710.440.240.34
    Nitrogen dioxide (ppb)811.514.518.322.426.710.420.450.58
    Ozone (ppb)1338.548.359.373.486.410.690.33
    PM10 (μg/m3)614.021.328.536.545.010.38
    Sulfur dioxide (ppb)93.95.37.410.713.6
Toronto, Ontario, Canada (December 10, 1993–July 25, 1995)
    Carbon monoxide (ppm × 10)45.78.09.711.312.70.630.160.370.55
    Nitrogen dioxide (ppb)417.521.326.031.537.010.400.640.63
    Ozone (ppb)427.533.643.556.567.110.700.44
    PM10 (μg/m3)16.812.017.727.237.810.48
    Sulfur dioxide (ppb)
4
0.2
1.0
2.5
4.8
8.8




*

Data on concentrations of particulate matter less than 10 μm in aerodynamic diameter (PM10) were available for a subset of study days in Baltimore, Boston, Denver, San Diego, St. Louis, and Toronto. Percentiles and correlations with other pollutants were based on an imputation sample used in the primary analysis for PM10. Ozone values and correlations with other pollutants were limited to the warm season, which lasted from May through September in 1994 and 1995.

To be included in the calculation of daily ambient concentrations, monitors must have been functioning for at least half of the study days. Since PM10 was not measured daily in all cities, monitors were included in the calculation of daily ambient concentrations if they were functioning for at least 90 of the study days.

Daily, monitor-specific values were 24-hour averages for carbon monoxide, nitrogen dioxide, and sulfur dioxide and a 1-hour maximum for ozone. Daily, citywide ambient concentrations for all pollutants were calculated with an adjusted average over monitor-specific values.

§

Individual seasons used in the adjustments are given by consecutive 45-day time windows.

ppm, parts per million; ppb, parts per billion.

One-pollutant models

Study-wide estimates for asthma symptoms and rescue inhaler uses from one-pollutant models are displayed in figures 1 and 2, respectively. Lags in within-subject carbon monoxide and nitrogen dioxide concentrations showed the most evidence of a positive relation with both responses. The 3-day moving sum of sulfur dioxide was also positively associated with symptoms, but PM10 and warm-season ozone were not associated with exacerbations.

FIGURE 1.

Odds ratios for daily asthma symptoms associated with shifts in within-subject concentrations of one air pollutant, Childhood Asthma Management Program, November 1993–September 1995. Effect sizes reported correspond to the following shifts: carbon monoxide, 1 part per million; nitrogen dioxide, 20 parts per billion (ppb); ozone, 30 ppb; particulate matter less than 10 μm in aerodynamic diameter (PM10), 25 μg/m3; and sulfur dioxide, 10 ppb. All city-specific estimates of pollutant effects were included in calculations of study-wide effects except sulfur dioxide in Albuquerque, New Mexico, and nitrogen dioxide in Seattle, Washington. Models corresponding to ozone were limited to the warm season (May through September in 1994 and 1995). Horizontal lines, 95% confidence interval (with limits specified at ends).

FIGURE 2.

Rate ratios for number of rescue inhaler uses associated with shifts in within-subject concentrations of one air pollutant, Childhood Asthma Management Program, November 1993–September 1995. Effect sizes reported correspond to the following shifts: carbon monoxide, 1 part per million; nitrogen dioxide, 20 parts per billion (ppb); ozone, 30 ppb; particulate matter less than 10 μm in aerodynamic diameter (PM10), 25 μg/m3; and sulfur dioxide, 10 ppb. All city-specific estimates of pollutant effects were included in calculations of study-wide effects except sulfur dioxide in Albuquerque, New Mexico, and nitrogen dioxide in Seattle, Washington. Models corresponding to ozone were limited to the warm season (May through September in 1994 and 1995). Horizontal lines, 95% confidence interval (with limits specified at ends).

One-ppm differences in within-subject carbon monoxide concentrations were associated with symptom odds ratios of 1.08 (95 percent CI: 1.01, 1.14), 1.07 (95 percent CI: 0.99, 1.16), 1.08 (95 percent CI: 1.02, 1.15), and 1.05 (95 percent CI: 1.01, 1.09) for the 0-, 1-, and 2-day lags and the 3-day moving sum, respectively. Odds ratios associated with 20-ppb differences in nitrogen dioxide were 1.06 (95 percent CI: 1.00, 1.13), 1.04 (95 percent CI: 0.97, 1.10), 1.09 (95 percent CI: 1.03, 1.15), and 1.04 (95 percent CI: 1.01, 1.07), respectively. Although sulfur dioxide lags were positively related to risk of asthma symptoms, only the 3-day moving sum effect would be considered (marginally) statistically significant at the 0.05 level. A 1-ppb change in the 3-day moving sum of sulfur dioxide levels led to a 1.04-fold increase in the odds of symptoms (95 percent CI: 1.00, 1.08).

The results from single-pollutant models of rescue inhaler use were consistent with those of asthma symptoms. Relative rates of rescue inhaler use associated with a 1-ppm change in 0-, 1- and 2- day lags of carbon monoxide and the 3-day moving sum of carbon monoxide were 1.07 (95 percent CI: 1.01, 1.13), 1.05 (95 percent CI: 0.99, 1.10), 1.06 (95 percent CI: 1.01, 1.10), and 1.04 (95 percent CI: 1.01, 1.07), respectively, while 20-ppb shifts in the 0-, 1-, and 2-day lags of nitrogen dioxide and the 3-day moving sum of nitrogen dioxide led to relative rates of rescue inhaler use of 1.04 (95 percent CI: 1.00, 1.08), 1.03 (95 percent CI: 0.99, 1.08), 1.05 (95 percent CI: 1.01, 1.09), and 1.03 (95 percent CI: 1.01, 1.05), respectively. Sulfur dioxide was not associated with rescue inhaler use rates.

Two-pollutant models

Results from two-pollutant models were based on the sum of the two within-subject pollutant effects, as described in the “Analysis approach” section. Thus, they are intended to provide insight into the increased risk of asthma exacerbation associated with (simultaneous) shifts in two pollutants. As with the one-pollutant models, effects are characterized by increases in the odds of asthma symptoms and the rates of rescue inhaler use.

The effects of shifts in pairs of carbon monoxide, nitrogen dioxide, PM10, and sulfur dioxide levels on asthma symptoms and rescue inhaler uses are displayed in figures 3 and 4, respectively. Ozone was not included in two-pollutant models, since it is a warm-season pollutant, and our primary goal was to study year-round effects. While all lags in the pairs of pollutants exhibited positive relations with both responses, simultaneous shifts in 2-day lags of pairs of pollutants tended to have the largest impact on asthma aggravation. Among combinations of carbon monoxide, nitrogen dioxide, PM10, and sulfur dioxide, the 2-day lags in all models involving either carbon monoxide or nitrogen dioxide were associated with both responses, with 95 percent confidence intervals excluding 1. Simultaneous standardized shifts in the 2-day lags of these two pollutants led to a 1.09-fold increase in the odds of asthma symptoms (95 percent CI: 1.02, 1.16) and a 1.07-fold increase in the daily rate of rescue inhaler uses (95 percent CI: 1.02, 1.12). Simultaneous shifts in 3-day moving sums of two pollutants were also frequently (marginally) related to exacerbations.

FIGURE 3.

Odds ratios for daily asthma symptoms associated with shifts in within-subject concentrations of two air pollutants, Childhood Asthma Management Program, November 1993–September 1995. Effect sizes reported correspond to the following shifts: carbon monoxide, 1 part per million; nitrogen dioxide, 20 parts per billion (ppb); particulate matter less than 10 μm in aerodynamic diameter (PM10), 25 μg/m3; and sulfur dioxide, 10 ppb. All city-specific estimates of pollutant effects were included in calculations of study-wide effects except sulfur dioxide in Albuquerque, New Mexico, and nitrogen dioxide in Seattle, Washington. Horizontal lines, 95% confidence interval (with limits specified at ends).

FIGURE 4.

Rate ratios for number of rescue inhaler uses associated with shifts in within-subject concentrations of two air pollutants, Childhood Asthma Management Program, November 1993–September 1995. Effect sizes reported correspond to the following shifts: carbon monoxide, 1 part per million; nitrogen dioxide, 20 parts per billion (ppb); particulate matter less than 10 μm in aerodynamic diameter (PM10), 25 μg/m3; and sulfur dioxide, 10 ppb. All city-specific estimates of pollutant effects were included in calculations of study-wide effects except sulfur dioxide in Albuquerque, New Mexico, and nitrogen dioxide in Seattle, Washington. Horizontal lines, 95% confidence interval (with limits specified at ends).

Although the relation between the responses and the 0- and 1-day lags in ambient concentrations did not appear to be as strong as the 2-day lags, they were consistent with them (e.g., estimated odds ratios for symptoms and rate ratios for inhaler uses were greater than 1). The strongest among the 0- and 1-day lagged effects was observed with asthma symptoms: Shifts in the 0-day lag of carbon monoxide and PM10 resulted in an odds ratio of 1.08 (95 percent CI: 1.01, 1.15).

DISCUSSION

The goal of this panel-study analysis was to describe the relations between ambient concentrations of five of the Environmental Protection Agency's criteria pollutants and asthma exacerbation in children with mild-to-moderate asthma. Two important features of this study are that children lived in cities widely separated across North America and the entire observation period covered approximately 22 months. Thus, results reflect ambient air pollution effects averaged across seasons and geographic regions. Most previous panel studies of this subject examined these relations within a season and/or geographic region (112). Although the relations we estimated are implicitly different from those reported in previous research, acknowledging the large exposure differences used to characterize effect size, the observed associations were generally smaller than those from previous reports.

Among the pollutants studied, carbon monoxide and nitrogen dioxide appeared to capture the most relevant health information. Sulfur dioxide showed less evidence for a relation with asthma exacerbations unless it was considered in combination with carbon monoxide or nitrogen dioxide. There was no evidence of a warm-season effect of ozone or of a year-round effect of PM10. Mobile-source emissions contribute directly to concentrations of ambient carbon monoxide, nitrogen dioxide, and particulate matter (especially the fine fraction of PM10) and secondarily to ozone concentrations. In the cities we studied, carbon monoxide and nitrogen dioxide were highly correlated with one another, which is consistent with mobile emissions being a primary source of these pollutants. Carbon monoxide can interfere with oxygen transport through the formation of carboxyhemoglobin, and nitrogen dioxide has been shown to sensitize reactions to inhaled allergens (2931); however, the observed association between carbon monoxide and nitrogen dioxide and risk of symptoms and rescue inhaler use should perhaps be interpreted with caution. Correlates of carbon monoxide and nitrogen dioxide emanating from mobile sources, including the fine fraction of PM10, which was not measured explicitly in the majority of our study cities, have also been shown to be associated with asthma exacerbation in children (3, 7, 8, 12). Another report suggested that ambient carbon monoxide, and perhaps nitrogen dioxide in the winter, could serve as a surrogate measure for the fine fraction of ambient PM10 (32). We believe carbon monoxide and nitrogen dioxide effects are more likely the effects of other pollutants emanating from mobile sources. Air stagnation can, in certain circumstances, be a reasonable surrogate for fine particulate matter (33, 34). To evaluate the hypothesis that the carbon monoxide and nitrogen dioxide effects were driven by mobile-source emissions, we replicated the one-pollutant analyses using number of stagnant hours per day as the predictor of interest. Although all lags of stagnation showed positive associations with exacerbations, none were statistically significant at the 0.05 level.

In contrast to several past studies (112), ambient PM10 and warm-season ozone concentrations were unrelated to exacerbations. PM10 concentrations were measured on less than 50 percent of study days in all cities except Seattle and Albuquerque. Thus, referring to it as PM10 may be may be misleading, and perhaps “estimated PM10” would be more appropriate. Imputation techniques can be used to recover a fraction of missing data; however, we acknowledge that there was a significant amount of missing information that was not recovered. In single-pollutant models, the estimated relative increase in variance due to missing PM10 data in cities other than Seattle and Albuquerque ranged from 10 percent to 24 percent for asthma symptoms and from 5 percent to 33 percent for rescue inhaler uses. Attenuation of parameter estimates was also possible because of the error associated with imputing PM10 concentrations using other pollutants, meteorologic variables, and daily summaries of children's asthma experiences. In the presence of other pollutant concentrations and meteorologic variables, daily summaries of individual responses were not informative, and our choice to include them when constructing imputation data sets had no impact on our conclusions. While PM10 effects were not observed for the entire panel of children, they were observed in recent reports on the children participating at the Seattle center (8, 12).

While the lack of observed association between ozone and exacerbations during the warm season was not expected, the strength of this study was in the total number of children observed and in the length of the observation period. A total of 990 children were followed over the course of 22 months; however, on a given day, the average number of children observed was approximately 12 per city, making season-specific effects difficult to capture. More thorough city-specific analyses may also be appropriate for the analysis of ozone.

Although year-round follow-up of children made this study unique and allowed us to address novel questions, relative to many panel studies, it was more susceptible to seasonal confounding. We constructed regression models with this concern in mind and took care to adjust for seasonal factors. We also conducted sensitivity analyses, examining more flexible seasonal confounding control (e.g., more seasons), which led to similar qualitative conclusions. Nevertheless, we must acknowledge that seasonal confounding may persist.

The National Center for Health Statistics estimated that, in 2003, 12 percent of children under age 18 years had ever been diagnosed with asthma and 6 percent had had an asthma attack in the past 12 months (35). In another study, children aged 13–14 years from three US study sites (in two cities) reported wheezing in the past 12 months and ever having asthma at rates of 20 percent and 17 percent, respectively (36, 37). In 1996, an estimated 2.5 million US children sought health care for asthma (38) and 6.3 million school absence days were attributed to asthma (38), and in 1998, the cost of asthma in the United States was estimated to be $12.8 billion per year (39). While asthma prevalence had been on the rise over the past several decades (40), this trend appears to have plateaued in recent years (41). We observed year-round associations between ambient concentrations of carbon monoxide and nitrogen dioxide and asthma exacerbations which we believe may represent mobile-source emission (fine particulate matter) effects. Though these relations are modest, they are clinically relevant, considering the enormous burden asthma places on the population at large. Further investigation is required, however, to explore the potential surrogacy of carbon monoxide and nitrogen dioxide for mobile-source emissions.

This research was partially funded by Environmental Protection Agency grant R827355. The Childhood Asthma Management Program is supported by contracts N01-HR-16044, -16045, -16046, -16047, -16048, -16049, -16050, -16051, and -16052 with the National Heart, Lung, and Blood Institute.

A portion of Dr. Jonathan Schildcrout's salary is funded by Pfizer Inc. (New York, New York).

References

1.

Delfino RJ, Coate BD, Zeiger RS, et al. Daily asthma severity in relation to personal ozone exposure and outdoor fungal spores.

Am J Respir Crit Care Med
1996
;
154
:
633
–41.

2.

Delfino RJ, Zeiger RS, Seltzer JM, et al. Symptoms in pediatric asthmatics and air pollution: differences in effects by symptom severity, anti-inflammatory medication use and particulate averaging time.

Environ Health Perspect
1998
;
106
:
751
–61.

3.

Ostro B, Lipsett M, Mann J, et al. Air pollution and exacerbation of asthma in African-American children in Los Angeles.

Epidemiology
2001
;
12
:
200
–8.

4.

Ostro B, Lipsett MJ, Mann JK, et al. Air pollution and asthma exacerbation among African-American children in Los Angeles.

Inhal Toxicol
1995
;
7
:
711
–22.

5.

Peters A, Goldstein IF, Beyer U, et al. Acute health effects of exposure to high levels of air pollution in Eastern Europe.

Am J Epidemiol
1996
;
144
:
570
–81.

6.

Romieu I, Meneses F, Ruiz S, et al. Effects of intermittent ozone exposure on peak expiratory flow and respiratory symptoms among asthmatic children in Mexico City.

Arch Environ Health
1997
;
52
:
368
–76.

7.

Romieu I, Meneses F, Ruiz S, et al. Effects of air pollution on the respiratory health of asthmatic children living in Mexico City.

Am J Respir Crit Care Med
1996
;
154
:
300
–7.

8.

Slaughter JC, Lumley T, Sheppard L, et al. Effects of ambient air pollution on symptom severity and medication use in children with asthma.

Ann Allergy Asthma Immunol
2003
;
91
:
346
–53.

9.

Thurston GD, Lippmann M, Scott MB, et al. Summertime haze air pollution and children with asthma.

Am J Respir Crit Care Med
1997
;
155
:
654
–60.

10.

Timonen KL, Pekkanen J. Air pollution and respiratory health among children with asthmatic or cough symptoms.

Am J Respir Crit Care Med
1997
;
156
:
546
–52.

11.

Vedal S, Petkau J, White R, et al. Acute effects of ambient inhalable particles in asthmatic and nonasthmatic children.

Am J Respir Crit Care Med
1998
;
157
:
1034
–43.

12.

Yu O, Sheppard L, Lumley T, et al. Effects of ambient air pollution on symptoms of asthma in Seattle-area children enrolled in the CAMP study.

Environ Health Perspect
2000
;
108
:
1209
–14.

13.

Mortimer KM, Neas LM, Dockery DW, et al. The effect of air pollution on inner-city children with asthma.

Eur Respir J
2002
;
19
:
699
–705.

14.

Long-term effects of budesonide or nedocromil in children with asthma. The Childhood Asthma Management Program Research Group.

N Engl J Med
2000
;
343
:
1054
–63.

15.

Ward DJ, Ayres JG. Particulate air pollution and panel studies in children: a systematic review.

Occup Environ Med
2004
;
61
:
e13
. (Electronic article).

16.

DerSimonian R, Laird N. Meta-analysis in clinical trials.

Control Clin Trials
1986
;
7
:
177
–88.

17.

Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models.

Biometrika
1986
:
13
–22.

18.

Zeger SL, Liang KY. Longitudinal data analysis for discrete and continuous outcomes.

Biometrics
1986
;
42
:
121
–30.

19.

Pepe MS, Anderson GI. A cautionary note on inference for marginal regression models with longitudinal data and general correlated response data.

Commun Stat Simul Comput
1994
;
23
:
939
–51.

20.

Schildcrout JS, Heagerty PJ. Regression analysis of longitudinal binary data with time-dependent environmental covariates: bias and efficiency.

Biostatistics
2005
;
6
:
633
–52.

21.

Sheppard L. Insights on bias and information in group-level studies.

Biostatistics
2003
;
4
:
265
–78.

22.

Neuhaus JM, Kalbfleisch JD. Between- and within-cluster covariate effects in the analysis of clustered data.

Biometrics
1998
;
54
:
638
–45.

23.

Rubin DB. Multiple imputation for nonresponse in surveys. New York, NY: John Wiley & Sons, Inc,

1987
.

24.

SAS/Stat user's guide, release 8. Cary, NC: SAS Institute, Inc,

2001
.

25.

R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing,

2006
.

26.

Yan J, Fine J. Estimating equations for association structures.

Stat Med
2004
;
23
:
859
–74.

27.

Lumley T. mitools: tools for multiple imputation of missing data. (R package, version 1.0). Vienna, Austria: R Foundation for Statistical Computing,

2004
. (http://cran.r-project.org/doc/packages/mitools.pdf).

28.

Lumley T. rmeta: meta-analysis. (R package, version 2.12). Vienna, Austria: R Foundation for Statistical Computing,

2005
. (http://cran.r-project.org/doc/packages/rmeta.pdf).

29.

Bernstein JA, Alexis N, Barnes C, et al. Health effects of air pollution.

J Allergy Clin Immunol
2004
;
114
:
1116
–23.

30.

Strand V, Rak S, Svartengren M, et al. Nitrogen dioxide exposure enhances asthmatic reaction to inhaled allergen in subjects with asthma.

Am J Respir Crit Care Med
1997
;
155
:
881
–7.

31.

Tunnicliffe WS, Burge PS, Ayres JG. Effect of domestic concentrations of nitrogen dioxide on airway responses to inhaled allergen in asthmatic patients.

Lancet
1994
;
344
:
1733
–6.

32.

Sarnat JA, Brown KW, Schwartz J, et al. Ambient gas concentrations and personal particulate matter exposures: implications for studying the health effects of particles.

Epidemiology
2005
;
16
:
385
–95.

33.

Norris G, Larson T, Koenig J, et al. Asthma aggravation, combustion, and stagnant air.

Thorax
2000
;
55
:
466
–70.

34.

Slaughter JC, Kim E, Sheppard L, et al. Association between particulate matter and emergency room visits, hospital admissions and mortality in Spokane, Washington.

J Expo Anal Environ Epidemiol
2005
;
15
:
153
–9.

35.

Dey AN, Bloom B. Summary health statistics for U.S. children: National Health Interview Survey,

2003
. (Vital and health statistics, series 10). Hyattsville, MD: National Center for Health Statistics, 2005:1–78.

36.

Worldwide variations in the prevalence of asthma symptoms: The International Study of Asthma and Allergies in Childhood (ISAAC).

Eur Respir J
1998
;
12
:
315
–35.

37.

Worldwide variation in prevalence of symptoms of asthma, allergic rhinoconjunctivitis, and atopic eczema: ISAAC. The International Study of Asthma and Allergies in Childhood (ISAAC) Steering Committee.

Lancet
1998
;
351
:
1225
–32.

38.

Wang LY, Zhong Y, Wheeler L. Direct and indirect costs of asthma in school-age children.

Prev Chronic Dis
2005
;
2
:
A11
. (Electronic article).

39.

Weiss KB, Sullivan SD. The health economics of asthma and rhinitis. I. Assessing the economic impact.

J Allergy Clin Immunol
2001
;
107
:
3
–8.

40.

Beasley R. The burden of asthma with specific reference to the United States.

J Allergy Clin Immunol
2002
;
109
(suppl):
S482
–9.

41.

Akinbami LJ, Schoendorf KC. Trends in childhood asthma: prevalence, health care utilization, and mortality.

Pediatrics
2002
;
110
:
315
–22.