If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Natural Resources Defense Council and Columbia Climate and Health Program, New York, New YorkDepartment of Environmental Health Sciences, Columbia University Mailman School of Public Health, New York, New York
The adverse respiratory effects of ground-level ozone are well established. Ozone is the air pollutant most consistently projected to increase under future climate change.
Purpose
To project future pediatric asthma emergency department visits associated with ground-level ozone changes, comparing 1990s to 2020s.
Methods
This study assessed future numbers of asthma emergency department visits for children aged 0–17 years using (1) baseline New York City metropolitan area emergency department rates; (2) a dose–response relationship between ozone levels and pediatric asthma emergency department visits; and (3) projected daily 8-hour maximum ozone concentrations for the 2020s as simulated by a global-to-regional climate change and atmospheric chemistry model. Sensitivity analyses included population projections and ozone precursor changes. This analysis occurred in 2010.
Results
In this model, climate change could cause an increase in regional summer ozone-related asthma emergency department visits for children aged 0–17 years of 7.3% across the New York City metropolitan region by the 2020s. This effect diminished with inclusion of ozone precursor changes. When population growth is included, the projections of morbidity related to ozone are even larger.
Conclusions
The results of this analysis demonstrate that the use of regional climate and atmospheric chemistry models make possible the projection of local climate change health effects for specific age groups and specific disease outcomes, such as emergency department visits for asthma. Efforts should be made to improve on this type of modeling to inform local and wider-scale climate change mitigation and adaptation policy.
Introduction
Asthma, the most common chronic disease seen in pediatrics, is a genetic and environmental disorder.
The lungs are exposed to air pollution more than other organ systems particularly for children because of their increased minute ventilation compared to adults. As such, asthma exacerbations may be viewed as the “canary in the coal mine” for children's environmental health. In recent years, public health literature has highlighted respiratory illness as a potential future health impact related to climate change.
in: Ezzati M. Lopez A. Rodgers A. Murray C. Comparative quantification of health risks: global and regional burden of disease due to selected major risk factors. WHO,
Geneva2004: 1543-1649
no models to date have produced regional projections of future pediatric asthma related to climate-driven changes in the U.S.
Climate change through both temperature and wind-pattern changes is projected to affect multiple air-pollutant levels and specifically ground-level ozone (O3).
Ozone is a secondary pollutant formed through photochemical reactions involving other pollutants such as nitrogen oxides (NOx) and volatile organic chemicals (VOC).
For example, during the Atlanta Olympics in 1996 when peak daily O3 dropped 28%, there was an 11% reduction in pediatric emergency department visits for asthma and an over-40% reduction in acute care asthma events not seen with non-asthma acute care events.
Asthma emergency department visits are one manifestation of uncontrolled asthma that is associated with both direct (services and medication) and indirect (such as missed work by parents) costs.
The objective of the present study was to demonstrate one method of projecting climate-driven, ozone-related pediatric asthma effects for the 2020s in an urban area of the U.S.
Methods
Climate Change and O3 Modeling
Projections of ground-level O3 were developed on a 36 × 36-km grid over the New York City metropolitan area by linking models for global climate, regional climate, and regional air quality.
Briefly, global climate was modeled using the Goddard Institute for Space Studies general circulation model. The regional climate model used the Penn State/National Center for Atmospheric Research Mesoscale Model 5. The atmospheric chemistry and O3 simulations came from the Community Multiscale Air Quality model. The specific greenhouse gas emissions projection came from the A2 scenario of the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES). The IPCC SRES involve future assumptions about energy use; population growth; and political, environmental, and social development and have been used to help standardize climate modeling. The A2 scenario assumes relatively rapid growth in greenhouse gas emissions and population.
previously published a description of these models, the greenhouse gas emission scenarios, the O3 simulations, the model evaluation compared to historic values, and uncertainties within models.
The present study focused on the 14 New York State counties that are considered part of the New York City (NYC) metropolitan area (Appendix A, available online at www.ajpmonline.org). The regional model simulated O3 for the summer season (June–August) for five consecutive mid-decadal years (e.g., 1993–1997) in the 1990s and 2020s.
Although the entire period of May–September can have high O3 levels, this study was restricted to model outputs for June–August because of computational constraints. The health outcome model used daily 8-hour maximum O3 concentrations, the same metric specified in the National Ambient Air Quality Standards 2008 revision.
This study used a health-impact assessment framework to assess changes in O3-related asthma emergency department visits in the population aged 0–17 years in the 2020s compared with the 1990s.
For a typical summer of each decade, county-level impacts were computed as
where M is the mean number of daily asthma emergency department visits among children aged 0–17 years attributable to daily 8-hour maximum O3 concentrations; P is the county population in that age group during time period of interest; B is the baseline county-level daily asthma emergency department visit rate among children aged 0–17 years in June–August (per 10,000 population); ERC is the exposure–risk coefficient of asthma morbidity for a given change in the mean O3 metric, as follows:
where b is the parameter estimate that reflects a 4% change in asthma emergency department visits per 20 parts per billion (ppb) change in ground-level O3 (daily 8-hour maximum) derived from Tolbert et al.,
examined associations between air pollutants and emergency department visits for children aged 0–16 years in Atlanta. Ideally, the ERC comes from pooled results that are geographically, diagnostically, and age matched to the study population. However, no such pooled study exists; thus, after considering the relevant epidemiologic literature, this study's ERC was derived from a single study.
to capture asthma events (asthma, ICD-9 code 493); wheezing (ICD-9 code 786.09); and reactive airway disease (ICD-9 code 519.1, relevant in only 1993 prevalence data) for pediatric emergency department visits among children aged 0–16 years. That study calculated the association of ozone and asthma events with 1-day lag. This study used their relative risk estimate for asthma of 1.04 per 20 ppb increase in daily maximum 8-hour ozone (95% CI=1.008, 1.073) to calculate the ERC.
Once the mean daily O3-related morbidity was derived, the summer morbidity was calculated by multiplying by 92 for the number of days in June–August. County populations in the mid-1990s were obtained from 2000 U.S. Census data.
Population was held constant in the 2020s for base-case calculations with the intent of isolating the climate influence and reserving exploration of this expected underestimation for the sensitivity analysis described below.
Baseline Morbidity Rate Estimation
This study used publicly available New York State Department of Health asthma emergency department visit data, which define an asthma emergency department visit as a primary emergency department diagnosis or an admitting diagnosis from the emergency department of ICD-9 code 493 for any age.
Use of this diagnosis in children aged <3 years is problematic because of the difficulty of asthma diagnosis in this age group, but more-refined population-level data were not available.
Despite these limitations of the state data set, these data were the most appropriate for this model. Average daily summer pediatric asthma emergency department visits for each of the 14 counties were calculated by adjusting the annual age-specific rates with a summer scaling factor. The scaling factor of 0.137 was derived from the proportion of annual asthma emergency department visits that occurred in June–August in New York State.
The summer seasonal rate was then converted to a daily rate. For comparison, the average summer scaling value for asthma hospitalizations for children aged 0–4 years in New York State from 1990 to 2004 was calculated to be 0.138.
The comparable summer asthma hospitalizations for all ages, including adults, from 2000 to 2005 in New York State was 0.173; these comparisons suggest that the scaling factor used was reasonable.
These scaling calculations were made because the historical daily pediatric emergency department data for these counties (i.e., at a comparable temporal and geographic scale) were not available to these researchers. Baseline morbidity rates were held constant in all analyses. Although morbidity rates will undoubtedly change in the future in response to changes in disease management, access to preventive and chronic disease care, and changing demographics (age, race/ethnicity, SES), projection of these shifts was beyond the scope of this study.
Impact Assessments
One primary health outcome assessment and two sensitivity analyses were performed. The health outcome assessment (HOA) used models of future O3 concentrations to project asthma emergency department visits resulting only from climate change with altered greenhouse gas emissions under the A2 SRES scenario. Sensitivity analyses examined alterations in two of the assumptions underlying the primary assessment.
Health Outcomes Assessment: 2020s Climate Change Only
The objective here was to assess how climate change alone might contribute to changes in summer O3 concentrations and associated pediatric asthma emergency department visits in the New York region over the next 10–15 years. County population totals were held constant at Census 2000 levels.
Similarly, anthropogenic O3 precursor emissions were held constant consistent with the 1996 county-level U.S. Environmental Protection Agency (EPA) National Emissions Trends inventory. Thus, no projected changes in anthropogenic precursor emissions were applied in the Community Multiscale Air Quality model projections of 2020s summer O3. The base case did allow for temperature-dependent changes in biogenic and mobile source emissions. No threshold for O3 impacts was assumed.
Two sensitivity analyses were conducted to better understand the effects on projections of changing individual modeling assumptions. The first analysis, S1, included the same assumptions of the HOA model but added population growth projections. For S1 in which population was allowed to grow, 2020s populations used the county and age-specific population projections from Cornell University Program on Applied Demographics.
This method projected a regional population increase in 2025 compared to 2000 of 3.9% among those aged 0–17 years with a range of negative 14.2% in suburban counties to positive 20% among the urban core counties.
The second analysis, S2, explored the effects of increasing anthropogenic O3-precursor emissions in addition to climate change. The first step used projections from a separate run of the climate model that included climate change and emissions changes that would be consistent with the A2 SRES story line (increases of 29.5% for all anthropogenic NOx emission sources and 9% for all anthropogenic VOC emissions). Of note, NOx and VOC have documented respiratory health effects but this model captures only O3-related impacts.
With this second set of O3 projections, isolation of the emissions effect was then possible by adding the difference of the two projections from the 2020s to the estimates of O3 from the 1990s and using those projections to model the health projections. These analyses occurred in 2010.
Results
Health Outcomes Assessment
County-specific O3 concentrations and the associated pediatric morbidity projections for emergency department visits for the primary HOA appear in Table 1. Percentage changes in emergency department visits between the 1990s and 2020s for the HOA model appear in Figure 1. The calculations of projections for O3-related asthma morbidity include the 95% CI from the epidemiologic study of Tolbert et al.
The O3 projections for average summer daily 8-hour maximum concentrations under the HOA model—which included only climate change—increased their median by 2.7–5.3 ppb across the 14 counties. The distribution of the O3 changes shows greater increases in those counties outside the urban core toward the coast and along the predominant air mass trajectory from the southwest. As a direct result, the distribution of percentage changes for emergency department visits follows a similar pattern (Figure 1). The underlying population of each county strongly influences the absolute numbers. Median regional O3-attributable asthma pediatric emergency department visits increased 7.3% in the 2020s compared to the 1990s. This analysis suggests that projected increases in emergency department visits occur both in the surrounding metropolitan counties as well as the urban core.
Table 1O3-related asthma emergency department visit change for children aged 0–17 years in New York City area
1990
2020s under A2 IPCC SRES climate change scenario (HOA)
HOA, Health Outcomes Assessment; IPCC SRES, Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios; O3, ozone; ppb, parts per billion
a Mean summer daily 8-hour maximum O3 concentration in ppb
Note: Percentage change is shown within each of 14 New York state counties of metropolitan NYC modeled under the HOA with A2 SRES.
A2, A2 scenario of IPCC SRES; HOA, health outcomes assessment; IPCC SRES, Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios
The results of the HOA model and the sensitivity analysis outcomes appear in Appendix B (available online at www.ajpmonline.org), which plots the median and 95th percentiles of the percentage change for emergency department visits. Sensitivity analysis S1 shows that population change accounts for the largest contribution to changes in the O3-attributable emergency department visits. This change is bidirectional as the population of those aged 0–17 years is projected to increase in urban counties—where asthma morbidity is already greatest—but decrease in some suburban counties. Including population growth projections (S1) resulted in a median of 10.6% increase in morbidity.
The second sensitivity analysis S2 isolated the effect of increased anthropogenic O3-precursor emissions. In the first step of this sensitivity analysis, the O3 projections from the model run that included both climate change and increased emissions showed a greater range by county, but the median increase in asthma emergency department visits was only 6.5%. The isolated emissions effect when climate change effect was removed also showed decreased emergency department visit projections compared to the HOA projections. This result occurred because of decreases in projected O3 concentrations for most of the urban metropolitan counties in this model, a result of air pollutant interactions that can cause O3 to decrease in areas where there are higher concentrations of other pollutants, such as NOx.
The median O3-related asthma emergency department visits in the S2 analysis was decreased by 0.4%.
Discussion
The results of this assessment suggest that, compared to the 1990s, by the 2020s climate change could cause a median increase of 7.3% in regional summer O3-related asthma emergency department visits for children aged 0–17 years across the New York City metropolitan region. However, when examining individual counties, O3-related emergency department increases ranged from 5.2% to 10.2%. Actual O3 precursor emissions have been decreasing as opposed to increasing—as detailed by U.S. EPA
—and updated precursor emissions assumptions would be useful in future analyses. A sensitivity analysis demonstrated the sometimes counterintuitive effects of air pollution dynamics. As with previous health-impact assessments using these climate models, these O3 simulations did not account for O3 precursor emission effects from outside of the modeled region (Eastern U.S.) on future air quality dynamics within the study area.
The assumption of uniform exposure to ozone within each county and the application of a single ERC to children of all ages are two other significant limitations of this study. Other U.S. studies have shown variable associations of O3 and asthma emergency department visits,
with a range of RR from 0.98 to 1.06 per 20 ppb increase in daily maximum 8-hour O3. These studies included a mixture of pediatric and adult patients and used lags ranging from 0 to 3 days. The study by Peel et al.
included the same data as Tolbert et al. but included all ages and the association was not significant. Confounding by other air pollutants and choice of lag days are discussed more fully within the original article and subsequent articles published by those same researchers.
In addition, the ERC came from a study in Atlanta, not New York, because of its availability as an age-matched study for asthma emergency department visits. To the extent that air conditioning is more prevalent in Atlanta than in New York (98% vs 85% according to 2003 and 2004 American Housing Survey data),
In this regard, the health-impacts modeling presented here is conservative. However, populations can adapt to climate through behavior changes such as increased time indoors or purchase of air conditioning units. Thus, levels of New York air conditioning use could match those of Atlanta in the future as the climate warms in which case such an ERC is arguably appropriate for projection of future health effects.
In a sensitivity analysis, an examination was made of the additional impact on future morbidity resulting from plausible increases in the exposed population. However, no examination was made of any scenarios of changing baseline morbidity rates, because of the high level of uncertainty regarding such estimates. It should be noted that projected cases scale directly to the baseline rate; thus a reduction by 50% in the baseline for asthma emergency department use would result in a 50% reduction in projected O3-related cases.
Prior modeling exercises similar to this one have mainly used mortality as an outcome; few have addressed morbidity or the unique vulnerabilities of children.
Another strength of this particular assessment is the use of county-specific inputs. Other pediatric respiratory effects have been well documented to be associated with O3 changes such as increased asthma hospitalization,
Effects of ambient air pollutants on asthma medication use and wheezing among fourth-grade school children from 12 Southern California communities enrolled in the Children's Health Study.
These additional documented health effects from O3 as well as a narrow definition for the baseline rates (i.e., excluding diagnoses of wheezing in young children) make these projections inherently conservative regarding the overall burden of disease from O3. Because the development of the chronic asthma phenotype is a function of both genotype and environmental insults,
the future incidence of asthma exacerbations and likely also emergency department visits will change as a function of the changing prevalence of chronic asthma (due in part to increased O3) as well as the increased rate of O3-induced exacerbations.
Clearly, this model simplifies complex relationships among change in climate, O3 precursors, population growth, asthma prevalence, healthcare utilization, and the exposure–risk coefficient of O3 and this specific outcome of asthma. Ethnicity and SES are important risk factors of uncontrolled asthma such as signaled by emergency department visits and are only indirectly included as part of the baseline morbidity estimates.
Improvements of these projections would include more geospatially refined and updated data for all of the above as well as different regional population and urbanization growth patterns, more current climate models—as will be used in the 2013 release of the IPCC's 5th Assessment Report, inclusion of bioaerosols such as pollen and other allergens, and additional scenarios beyond the single A2 SRES one used in this study.
Given that the projections of this single outcome of asthma emergency department visits are small in absolute number, the costs would not be staggering. However, given that a single asthma emergency department visit often represents far more uncontrolled asthma cases with associated increase in outpatient visits, medication use, and indirect costs such as missed school days for children and work days for parents, the combined costs would be much greater.
These projected effects add an important contribution to current research regarding climate-related disease in children. As the first model of climate-related, regional, pediatric morbidity, this study not only demonstrates an important modeling approach but also provides some quantitative projections to which future work can add and compare. Adaptation measures to climate change that work to reduce ozone levels should be coupled with ongoing efforts for better disease management of asthma. The projections from this study can inform the discussion of local, regional, and national policy.
PES received support from the NIH Research Training Program in Environmental Pediatrics (NIH 5T32 HD049311). PLK is supported by NIEHS Center grant ES09089 . This work was also partially supported by EPA STAR grant R-82873301 . The authors extend their gratitude to Philip J. Landrigan for his assistance in enabling this research to happen, and Christian Hogrefe of the Atmospheric Sciences Research Center at the University at Albany in New York for his guidance on interpreting the climate models.
No financial disclosures were reported by the authors of this paper.
in: Ezzati M. Lopez A. Rodgers A. Murray C. Comparative quantification of health risks: global and regional burden of disease due to selected major risk factors. WHO,
Geneva2004: 1543-1649
Effects of ambient air pollutants on asthma medication use and wheezing among fourth-grade school children from 12 Southern California communities enrolled in the Children's Health Study.