|
|
||||||||
RESEARCH |
Kristie L. Ebi is with EPRI, Palo Alto, Calif. K. Alex Exuzides, Edmund Lau, and Michael Kelsh are with the Exponent Health Group, Menlo Park, Calif. Anthony Barnston is with the International Research Institute for Climate Prediction, Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY.
Correspondence: Requests for reprints should be sent to Kristie L. Ebi, PhD, MPH, EPRI, 3412 Hillview Ave, Palo Alto, CA 94304 (e-mail: krisebi{at}epri.com).
| ABSTRACT |
|---|
|
|
|---|
Objectives. This study examined associations between weather and hospitalizations of females for viral pneumonia during normal weather periods and El Niño events in the California counties of Sacramento and Yolo, San Francisco and San Mateo, and Los Angeles and Orange.
Methods. Associations between weather and hospitalizations (lagged 7 days) for January 1983 through June 1998 were evaluated with Poisson regression models. Generalized estimating equations were used to adjust for autocorrelation and overdispersion. Data were summed over 4 days.
Results. Associations varied by region. Hospitalizations in San Francisco and Los Angeles increased significantly (30%50%) with a 5°F decrease in minimum temperature. Hospitalizations in Sacramento increased significantly (25%40%) with a 5°F decrease in maximum temperature difference. The associations were independent of season. El Niño events were associated with hospitalizations only in Sacramento, with significant decreases for girls and increases for women.
Conclusions. The results suggest that viral pneumonia could continue to be a major public health issue, with a significant association between weather and hospitalizations, even as the global mean temperature continues to rise. An understanding of population sensitivity under different weather conditions could lead to an improved understanding of virus transmission.
| INTRODUCTION |
|---|
|
|
|---|
It is projected that by 2100, global mean surface air temperatures will increase by 1°C to 3.5°C relative to 1900.3 The average rate of warming under the various scenarios considered by the Intergovernmental Panel on Climate Change will probably be greater than any seen in the last 10 000 years. The actual annual to decadal weather patterns will likely include considerable natural variability. Regional temperature changes could differ substantially from the global mean values. During the 1900s, the contiguous United States warmed by approximately 1°F and precipitation increased, with much of the increase due to increases in heavy precipitation events ( >5 cm/day) and decreases in light precipitation events.4,5 There is considerable interest, but a great deal of uncertainty, regarding whether increasing global mean temperatures will be associated with changes in the occurrence and intensity of El Niño Southern Oscillation events.6,7
The El Niño Southern Oscillation is a climate phenomenon that describes the periodic changes in oceanic and atmospheric circulation patterns in the tropical Pacific Ocean. El Niño and La Niña events are the 2 extremes of the El Niño Southern Oscillation. El Niño events, which are associated with a warming of the tropical Pacific sea surface temperature, occur approximately every 2 to 7 years. El Niño events can have a strong effect on local weather patterns. For example, in southern and central California, El Niño winters tend to get more rain than usual. During the years 1983 through 1998, there were 5 El Niño events (19821983, 19861987, 19911993, 19941995, and 19971998).
Although there is uncertainty as to whether the frequency and magnitude of El Niño events might change with increasing global mean temperatures, studying what happened during recent El Niño events provides information on the sensitivity of diseases to specific changes in weather patterns. El Niño events are "natural experiments" with worldwide impacts. Until recently, almost no research looked at how the distribution and rate of health outcomes varied, if at all, during El Niño events.8 Understanding these relationships is useful for developing current public health responses, for evaluating population vulnerability to these events, and for designing future adaptation measures.
This project was designed to describe and compare the existing associations between certain weather variables and hospitalizations for viral pneumonia during normal weather periods and El Niño events. Changing weather variables do not cause pneumonia, but they may set up conditions that facilitate increased or decreased viral transmission.
| METHODS |
|---|
|
|
|---|
We obtained weather data from the National Climatic Data Center of the National Oceanic and Atmospheric Administration. Data included daily minimum temperature, the difference between daily maximum and minimum temperatures, daily precipitation (in hundredths of an inch), monthly sea surface temperatures for 2 regions off the California coast, and monthly Southern Oscillation Index values. We focused on minimum rather than average temperatures because minimum temperatures are more informative for defining the cold season in California. We also included the maximum daily temperature difference, defined as the largest difference between the daily maximum and minimum temperatures.
The sea surface temperature data are from 10° quadrants of the world; one was centered on 40°N and 130°W (northern California), and the other was centered on 30°N and 120°W (southern California). Local sea surface temperatures are important predictors of general weather patterns in California and, as such, were used as markers of other weather variables not explicitly studied, such as wind and cloud cover. Southern Oscillation Index values (Niño Region 3.4 normalized air pressure differentials) were used to define normal weather periods (values of 0.5 to +0.5) and El Niño events (values of less than 0.5). We further categorized El Niño events as mild (less than 0.5 to 1.5) or strong (less than 1.5). By these definitions, during the study period there were 48 months with normal weather, 36 months with mild El Niño conditions, and 60 months with strong El Niño conditions.
The size of the data set required that we aggregate the weather data into larger units. We chose 4-day averages because large weather patterns typically take several days to move through the state and as a smoothing process. We defined 1 unit of temperature difference as a 5°F decrease, 1 unit of precipitation change as a 0.2-inch increase, and 1 unit of sea surface temperature change as a 1°C decrease. We chose these categories because they were about 1 SD from the mean. The maximum temperature difference was the largest difference between the maximum and minimum temperature for 1 day during the 4-day period.
We evaluated the associations between weather and hospitalizations, considering a 7-day lag period. The data were lagged before they were aggregated (i.e., hospitalizations on day 8 were associated with weather data on day 1; this was done for the entire data set before cutting the data into 4-day analysis units). The lag can be thought of as a temporal shift between 2 data streams. Because the 4-day grouping was not excessively broad relative to the lag, the 2 were largely independent of each other.
Hospitalization Data
Hospital admission data were provided by California's Office of Statewide Health Planning and Development for January 1983 through June 1998. Data are included from all nonfederal hospitals in the state. Records include data on patient demographics (age, sex, race, and zip code), date of admission, and disease diagnoses. We analyzed the broader category of viral pneumonia (International Classification of Diseases [ICD-9] 480, 487) instead of just influenza because of the possibility of misclassification from one category to another.
We analyzed data separately for males and females and considered 3 broad age groups: 0 to 17 years, 18 to 55 years, and 56 years and older. Here we show only the results for females. We analyzed total numbers of hospitalizations, not population rates, because using a rate would require adding the population as an offset parameter to the model described below. We felt this was unnecessary, because we included calendar year in the analysis to account for changes in hospitalizations over time due to changes in population age and density.
Statistical Methods
We evaluated the associations between weather variables and viral pneumonia hospitalizations in 3 geographic regions in California (the counties of Sacramento and Yolo, San Francisco and San Mateo, and Los Angeles and Orange, referred to as "Sacramento," "San Francisco," and "Los Angeles," respectively). These counties include portions of 3 of the major metropolitan areas in California. San Francisco and Los Angeles are coastal communities, about 400 miles apart, and Sacramento is inland, about 90 miles from San Francisco. Because weather patterns vary across the large geographic regions encompassed by the San Francisco and Los Angeles metropolitan areas, we decided to analyze only 2 counties in each area with similar weather.
We followed the statistical approach applied to air pollution and health outcomes by Morgan et al.9 and a collaborative European project.10 We used Poisson regression models that included a linear time factor to account for trends in hospital admissions over the 16-year period. We applied the generalized estimating equations (GEE) of Liang and Zeger to adjust for autocorrelation and overdispersion.11 Autocorrelation is an important factor that needs to be controlled for in time series data (weather measurements). Overdispersion may be present in count data (hospital admissions) that are assumed to follow a Poisson distribution.
The GEE approach is similar to generalized additive models in the use of a Poisson regression model to estimate daily hospital admission counts in relation to weather data. However, no a priori smoothing was performed for the weather time series. Instead, the GEE model allows for the removal of long-term patterns in the data by adjusting for overdispersion and autocorrelation. Furthermore, the analyst does not have to choose the optimal degree of smoothing by running a number of smoothing applications on the data. The main disadvantage of the GEE approach is that it can be computationally prohibitive when dealing with large databases. We used 4-day groupings for hospital admissions and weather variables both as a smoothing technique and to improve computational efficiency. All analyses were performed with SAS software (version 8.0; SAS Institute, Inc, Cary, NC).
The temperature and precipitation variables used in the statistical model may not fully account for all the important climatological factors during normal weather periods and El Niño events. To evaluate the impact of those latent factors on hospital admissions, we performed a sensitivity analysis by including a season indicator in our main model. Season was defined as winter (December, January, February), spring (March, April, May), summer (June, July, August), or fall (September, October, November). Our goal was to compare the association between weather variables and the El Niño indicator with hospital admissions, after adjustment for season, with the association without the season indicator included in the model. By the Fisher exact test, there was no association between season and El Niño events (P = .78). Therefore, inclusion of both season and El Niño indicators in our model should not introduce multicollinearity.
The sensitivity analysis was performed for Los Angeles and San Francisco. Although the season indicator was significant, the percentage change in hospital admissions due to weather variables, after adjustment for season, did not change substantially. The changes in hospital admissions were in the same direction and of about the same magnitude for both models (with and without the season indicator). The only exception was for temperature variables for the youngest group (<18 years) in San Francisco. There was a larger but not statistically significant reduction in hospital admissions in the model without the season indicator. Furthermore, the predictive power of our model did not increase significantly with the inclusion of the season indicator. As a result of this analysis, we present results without adjustment for season.
We assessed the goodness of fit of our models via deviance and Pearson
2 statistics. We estimated the percentage change in hospitalization per unit increase for each of the variables in the model and computed the corresponding 95% confidence intervals.
We performed 2 sets of analyses. The first set were strata-specific for normal weather periods and for El Niño events, to estimate the association of weather variables with hospitalizations for each of the periods separately. The second set included an El Niño indicator that defined mild and strong El Niño events and normal weather conditions. The purpose of this approach was to estimate the association of El Niño events with hospitalizations while controlling for other weather effects.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
In the comparison between El Niño events and normal weather periods in Sacramento, hospitalizations decreased for children by 21% to 28% and increased for adults by 10% to 48%. The associations were more pronounced during mild El Niño events.
The patterns of association were different in San Francisco (Table 4
). In this region, minimum temperature was strongly associated with hospitalizations; decreasing the minimum temperature by 5°F increased viral pneumonia hospitalizations by more than 30% in all age groups during both weather periods. As was found in Sacramento, sea surface temperature had a stronger association with hospitalizations in the younger age group (a 25% increase in hospitalizations during normal weather periods for each 1°C decrease in sea surface temperature); there was no significant association in adults. Maximum temperature difference was more important for adults (slightly decreased hospitalizations during normal weather periods, with increased hospitalizations during El Niño events). Hospitalizations did not change much between normal weather periods and El Niño events.
Similar patterns were found for Los Angeles (Table 5
), with 2 exceptions. First, weather variables had a weaker association with hospitalizations in Los Angeles than in San Francisco, although the association was still strong for decreasing minimum temperature. Second, the patterns for sea surface temperature and maximum temperature difference were weaker and inconsistent in adults. As was found for San Francisco, El Niño events had little additional association with hospitalizations in Los Angeles.
| DISCUSSION |
|---|
|
|
|---|
One temperature variable could not describe the hospitalization patterns found in the 3 geographic regions. In San Francisco and Los Angeles, a 5°F decrease in the minimum temperature was associated with a large increase in hospitalizations (about 30%50%). In Sacramento, a 5°F decrease in the maximum temperature difference was associated with a large increase in hospitalizations (about 25%40%).
There appeared to be no association of El Niño events with hospitalizations in San Francisco and Los Angeles, possibly because temperatures did not change much in these regions during El Niño events. In Sacramento, where seasonal temperatures did change during El Niño events, hospitalizations decreased among children and increased among adults. There is no obvious explanation for these patterns. Further investigation is required to understand the differences among the age groups in terms of changes in behavior, potential exposures, and other factors during El Niño events and normal weather periods. The recently developed capability to predict El Niño events means that public health interventions could be devised once these differences are understood.
More analyses are required to elucidate the relative importance of weather variables and specific weather patterns in the context of other factors associated with virus transmission.13 We did not have data on viral isolates and so were unable to determine whether patterns varied by virus type. Also, we did not have data on air pollution. Air pollution exposures could confound the associations present here, but it is not clear in which direction. For example, the increased precipitation during El Niño events may result in decreased exposure to some air pollutants. Much of the literature on air pollution and respiratory diseases has treated weather as a confounding variable to be controlled for in the analysis. Further analysis is warranted to discern the joint and separate associations of weather variables and air pollution with hospitalizations for viral pneumonia. Such analyses should take El Niño events into account where appropriate.
These data underscore the difficulties inherent in trying to understand the potential health effects of climate variability and change. The associations between viral pneumonia hospitalizations and specific weather factors varied across geographic regions. A model based on either the inland region or one of the coastal regions would not be predictive for the other regions.
There is interest in predicting how winter influenza mortality could change with warmer winters.1,2 There are, of course, significant problems with speculating on what might happen in a future world. Apart from weather, many social and economic changes are likely that could influence disease patterns. However, to inform decision makers and to prioritize research areas, it can be useful to consider what would happen if only the weather changed. In that unlikely event, if the current weatherhospitalization pattern found in Sacramento is similar to that for regions experiencing colder climates, then the results from El Niño events suggest that warmer winters may decrease hospitalizations for viral pneumonia in children and increase such hospitalizations in adults. El Niño events are short-term weather fluctuations, and the response of a population to these events may be different from the acclimatization response to long-term climate change. Different patterns could lead to different conclusions.
Although El Niño events were not associated with hospitalizations for viral pneumonia in San Francisco and Los Angeles, regions known for their moderate climates, there were associations with minimum temperature changes during both normal weather periods and El Niño events. These results do not directly address the question of whether warmer global temperatures could reduce winter influenza mortality. However, they do suggest that viral pneumonia would continue to have a seasonal pattern and a significant association with hospitalizations even as cities in colder climates begin to experience warmer temperatures.
Although it is not possible to directly intervene to alter the weather, understanding why specific weather factors are associated with increased or decreased hospitalizations can lead to an improved understanding of the dynamics of virus transmission, which can lead to better efforts to control disease outbreaks. One concern is how to reduce population vulnerability to the potential adverse consequences of climate variability and change. Vulnerability is a function of the sensitivity of a population to climate change and the ability of the population to adapt in anticipation of or in response to adverse impacts.14 A knowledge of ways to minimize risk or maximize adaptive capacity can be used to develop anticipatory adaptation measures. These results are a first step in understanding population sensitivity under different weather conditions.
| Acknowledgments |
|---|
We appreciate the scientific input into the design of these analyses provided by Guido Franco (California Energy Commission), Raymond Neutra (California Department of Health Services), and Bart Ostro (California Office of Environmental Health Hazard Assessment). We also appreciate the helpful comments of Susan Bernard (The Johns Hopkins University) and the reviewers.
| Footnotes |
|---|
K. L. Ebi planned the study and took the lead in writing the article. K. A. Exuzides developed the analysis approach, which was implemented by K. A. Exuzides and E. Lau. M. Kelsh contributed to the study design and analysis. A. Barnston provided climatological expertise and data. All authors contributed to evaluating the results and writing the paper.
Accepted for publication January 1, 2001.
| References |
|---|
|
|
|---|
2. McMichael AJ, Ando M, Carcavallo R, et al. Human population health. In: Watson RT, Zinyowera MC, Moss RH, eds. Climate Change 1995: Impacts, Adaptations, and Mitigation of Climate Change: Scientific-Technical Analyses. Contribution of Working Group II to the Second Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, England: Cambridge University Press; 1996:561584.
3. Houghton JT, Meira Filho LG, Callander BA, et al. Summary for policymakers. In: Houghton JT, Meira Filho LG, Callander BA, et al., eds. Climate Change 1995: The Science of Climate Change. Contribution of Working Group I to the Second Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, England: Cambridge University Press; 1996:37.
4. Karl TR, Knight RW, Easterling DR, Quayle RG. Indices of climate change for the United States. Bull Am Meteorol Soc. 1996;77:279303.
5. Karl TR, Knight RW. Secular trends of precipitation amount, frequency, and intensity in the USA. Bull Am Meteorol Soc. 1998;79:231241.
6. Karl TR, Knight RW, Plummer N. Trends in high-frequency climate variability in the twentieth century. Nature.1995;377:217220.
7. Meehl GA, Zwiers F, Evans J, et al. Trends in extreme weather and climate events: issues related to modeling extremes in projections of future climate change. Bull Am Meteorol Soc. In press.
8. Kovats S, Bouma M, Haines A. El Niño and Health. Geneva, Switzerland: World Health Organization Task Force on Climate and Health; 1999.
9.
Morgan G, Corbett S, Wlodarczyk J. Air pollution and hospital admissions in Sydney, Australia, 1990 to 1994. Am J Public Health.1988;88:17611766.
10. Katsouyanni K, Schwartz J, Spix C, et al. Short term effects of air pollution on health: a European approach using epidemiologic time series data: the APHEA protocol. J Epidemiol Community Health.1995;50(suppl 1):S12S18.
11.
Liang K, Zeger S. Longitudinal data analysis using generalised linear models. Biometrika.1986;73:1322.
12. Kilbourne ED. Influenza. New York, NY: Plenum Medical Book Co; 1987.
13. Hope-Simpson RE. The Transmission of Epidemic Influenza. New York, NY: Plenum Press; 1992.
14. Summary for policymakers. In: Watson RT, Zinyowera MC, Moss RH, eds. Climate Change 1995: Climate Change 1995: Impacts, Adaptations, and Mitigation of Climate Change: Scientific-Technical Analyses. Contribution of Working Group II to the Second Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, England: Cambridge University Press; 1996:118.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |