Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Frequency of Extreme Heat Event as a Surrogate Exposure Metric for Examining the Human Health Effects of Climate Change

  • Crystal Romeo Upperman,

    Affiliations Maryland Institute for Applied Environmental Health, University of Maryland School of Public Health, College Park, Maryland, United States of America, Marine Estuarine Environmental Science Program, University of Maryland, College Park, Maryland, United States of America

  • Jennifer Parker,

    Affiliation National Center for Health Statistics, Centers for Disease Control and Prevention, Hyattsville, Maryland, United States of America

  • Chengsheng Jiang,

    Affiliation Marine Estuarine Environmental Science Program, University of Maryland, College Park, Maryland, United States of America

  • Xin He,

    Affiliation Department of Epidemiology and Biostatistics, University of Maryland School of Public Health, College Park, Maryland, United States of America

  • Raghuram Murtugudde,

    Affiliation Department of Atmospheric and Oceanic Science, University of Maryland, College Park, Maryland, United States of America

  • Amir Sapkota

    amirsap@umd.edu

    Affiliation Maryland Institute for Applied Environmental Health, University of Maryland School of Public Health, College Park, Maryland, United States of America

Abstract

Epidemiological investigation of the impact of climate change on human health, particularly chronic diseases, is hindered by the lack of exposure metrics that can be used as a marker of climate change that are compatible with health data. Here, we present a surrogate exposure metric created using a 30-year baseline (1960–1989) that allows users to quantify long-term changes in exposure to frequency of extreme heat events with near unabridged spatial coverage in a scale that is compatible with national/state health outcome data. We evaluate the exposure metric by decade, seasonality, area of the country, and its ability to capture long-term changes in weather (climate), including natural climate modes. Our findings show that this generic exposure metric is potentially useful to monitor trends in the frequency of extreme heat events across varying regions because it captures long-term changes; is sensitive to the natural climate modes (ENSO events); responds well to spatial variability, and; is amenable to spatial/temporal aggregation, making it useful for epidemiological studies.

Introduction

Climate change is expected to cause approximately 250,000 deaths per year between 2030 and 2050 with direct damage costs totaling $ 2–4 billion USD per year by 2030 [1]. Chronic diseases, that may be exacerbated by climate change, disproportionately affect more vulnerable populations—including children, older adults, the socially isolated, and those with mental health issues [2]. Epidemiological investigation of the impact of climate change on human health is hindered by the differing temporal scale of the primary exposure of interest (climate change: decadal scale) and health outcomes that have varying sensitive time windows (days to years) in epidemiological studies that are based on a few years of data [3,4]. There is a need for a set of suitable exposure metrics that can capture the subtle attributes of changing climate (e.g. frequency, duration, and intensity of extreme events that are expected to rise), and would allow for comparisons across different geographical locations and time periods. Such an exposure metric should have enough flexibility for temporal aggregation to meet the needs of different types of epidemiological studies. Furthermore, for national health studies, there is a need for spatial compatibility, as meteorological data are often available at station level or a specific grid, whereas, national health data (e.g., behavioral risk factor surveillance system (BRFSS), Centers for Medicare and Medicaid Services (CMS), healthcare cost and utilization project (HCUP), CDC’s Public Health Tracking data) are commonly available at zip code, county or state levels, often with a non-uniform spatial resolution.

Public health researchers are increasingly using temperature measures (maximum, minimum, heat index, and apparent temperature), and heat wave episodes to identify the acute health outcomes associated with weather. While the linkage of frequency and intensity of heat waves with acute health outcomes provide important information about long-term climate trends and health, heat waves are still relatively uncommon in most locations in the US [5,6], limiting the generalizability of study results across locations and over time. Moreover, heat wave measurements are designed to capture physical phenomenon: most common definitions include a certain number of consecutive days exceeding a location-specific threshold [7]. Hence, by definition, heat wave does not capture isolated days where temperatures may have been high and, possibly, affecting health. Extreme heat events, on the contrary, will capture such isolated event.

For chronic health outcomes, the meaningful window of exposure may vary from months to several years. Not surprisingly, the relationships between climate change and chronic health outcomes are less understood than acute health outcomes, such as mortality or an emergency department (ED) visit, owing to the difficulty in defining appropriate exposure metrics that characterize underlying and long-term climate change at varying temporal and spatial resolutions and that are appropriate for chronic health outcomes [8,9].

In this paper we describe an indicator designed to capture exposure to climate variability and change at a spatial scale that is consistent with publically available county-level health outcome data. This indicator, henceforth referred to as an exposure metric or “extreme heat event”, captures positive anomaly derived using distributions of county- and month-specific climatology using a 30-year reference period (1960–1989). Although both extreme heat and extreme cold are important, this manuscript focuses exclusively on extreme heat because warmer temperature is associated with etiology of many infectious (higher rates of pathogen replication) as well as chronic diseases (increases in concentration of pollutant such as ozone). The exposure metric enables users to look at spatial and temporal changes over time using location specific baselines and serves as an additional resource to investigate the potential relationships between climate change and human health. We tabulate this exposure metric by time period, season, Census division, and 2006 urban-rural classification, documenting how the exposure metric is amenable to spatial and temporal aggregation across factors that are known to be associated with variability in temperature. Finally, we evaluate this exposure metric by assessing its correspondence to the different phases of El-Niño-Southern Oscillation (ENSO), a natural oscillation patterns that affect the weather phenomenon in the continental US and other parts of the world.

Materials and Methods

Meteorological data were acquired from the National Climatic Data Center (NCDC) branch of the National Oceanic and Atmospheric Association (NOAA) that maintains the world’s largest archive of meteorological data from the past 150 years. The data used to develop the metric are archived in two broad categories: DSI-3200 and DSI-3210. The DSI-3210 network is a smaller subset of DSI-3200 stations that collect several additional weather variables besides temperature and precipitation (e.g., barometric pressure, wind speed, wind direction), but has a poor spatial coverage. Therefore we chose the DSI-3200 database that contains approximately 8,000 active stations, with up to 23,000 stations for various years. The stations cover all 50 states plus Puerto Rico, US Virgin Islands and Pacific Island territories. Each dataset underwent quality control measures through both automated and manual edits by the NCDC, which consisted of internal consistency checks and evaluation against adjacent stations [10]. To develop and evaluate the metric, we used climate data for the 48 contiguous states and the District of Columbia. The county boundaries used for all years were defined by the 2000 Federal Processing Standards (FIPS) codes.

For Urban-Rural status, we used the 2006 county level National Center for Health Statistics (NCHS) Urban-Rural Classification Scheme. The 2006 NCHS urban-rural classification scheme was developed for use in studying and monitoring health disparities across the urban-rural continuum. The 2006 scheme consists of four levels of metropolitan counties (large central, large fringe, medium and small metro) and two levels of nonmetropolitan counties (micropolitan and non-core). This scheme is based on the December 2005 Office of Management and Budget delineations of county classification, Metropolitan Statistical Area (MSA) and principal city, MSA population-size cut points, and classification rules formulated by NCHS [11]. In congruence with other studies [1217], we opted to use this relatively recent classification scheme for the complete 51 years of data (1960–2010). However, because the 2006 scheme is not as accurate for earlier time periods as it is for more recent time periods [11], we also tabulated the extreme heat events using a 1990 urban rural classification [18] to conduct a sensitivity analysis.

The census division classification used in this study consists of groups of contiguous states as defined by the US Bureau of the Census (as: New England–CT, ME, MA, NH, RI, VT; Middle Atlantic–NJ, NY, PA; South Atlantic–DE, DC, FL, GA, MD, NC, SC, VA, WV; East South Central–AL, KY, MS, TN; West South Central–AR, LA, OK, TX; East North Central- Il, IN, MI, OH, WI; West North Central–IA, KS, MN, MO, NE, ND, SD; Mountain–AZ, CO, ID, MT, NV, NM, UT, WY; and Pacific–CA, WA, OR). Seasons were defined as: Winter–December, January, and February; Spring–March, April, and May; Summer–June, July, and August; and Autumn–September, October, and November.

El-Niño-Southern Oscillation (ENSO) indicator data—Oceanic Niño Index (ONI)—were obtained from the National Oceanic and Atmospheric Administration (NOAA), National Weather Service Climate Prediction Center. The Climate Prediction Center is a coordinated program that monitors, assesses and predicts climate phenomena and their linkage to weather events. Warm and cold episodes are based on a threshold of +/- 0.5°C for the Oceanic Niño Index (ONI)—a 3-month running median anomalies in the sea surface temperature in the Niño 3.4 region (5°N-5°S, 120°-170°W). The threshold values are based on centered 30-year base periods and are updated every 5 years. El Niño and La Niña episodes are defined when the threshold is met for a minimum of 5 consecutive over-lapping seasons. This ENSO indices data were used to categorize the months as La Niña, El Niño and Neutral months [19].

Extreme Heat Events

We assigned daily maximum temperature for all counties using the following rules: 1) average of daily maximum temperatures from all stations within the county, 2) if no station data were available from the county, the daily maximum temperature used was from the closest available station within a 20 km radius of the county boundary, and 3) a missing value was assigned if the previous two criteria were not met. In the complete dataset of observations, 99% of all counties had less than 1.5% missing data and there was no spatial pattern to the location of missing data. To compute extreme heat events, we used 1960–1989 as a reference period. This time period was chosen because the weather data were recorded consistently with current methods of NCDC measurement and the 30-year time period is generally accepted as the epoch (per the IPCC report) to represent the standardization of a climate regime [20]. For each county within the continental US, we compiled daily maximum temperatures (Tmax) by calendar months (e.g., Jan 1st to Jan 31st). For the 30-year reference period with no missing data, the total number of values would be approximately 900 observations (30 years by ~30 days in a month) for each county and calendar month. Using this distribution of daily Tmax, we calculated the month specific 95th percentile thresholds for each county. Using this cutoff value, we computed the calendar month and year specific extreme heat events for each county as: (1) where Ejk is the total number of extreme heat events for county j in calendar month k; Tijk-max is the daily maximum temperature (Tmax) in county j for day i of calendar month k; Tjk-95 is the 95th percentile Tmax value for county j for calendar month k for the 1960–1989 period; and Iijk represents the indicator of whether or not Tijk-max is greater than Tjk−95.

Evaluation

The units of analysis for our evaluation of the indicator were the annual and monthly total number of events; these are the metrics that are referenced throughout the paper. All covariates of interest were defined at the county, year and month levels. We computed descriptive statistics of the spatial (2006 urban-rural classification, Census division) and temporal (seasonal, decadal) characteristics. Additional descriptive statistics were calculated for ENSO periods. After checking the normality assumption, comparisons of means were performed using one-way analysis of variance (ANOVA) and post-hoc Tukey’s honest significant difference (HSD) tests [21]. We further investigated the temporal and spatial dependency of the exposure metric using negative binomial generalized estimation equation (GEE) models [22,23]. The year and monthly total extreme heat event anomalies in each county were modeled as a function of seasonality, ENSO, 2006 urban-rural classification, and Census division. We identified findings as statistically significant with a p-value <0.05. Most statistical analyses were performed using SAS 9.3 (SAS Institute, Cary, NC). In particular, PROC GENMOD was used to fit the negative binomial GEE models using a first-order autoregressive covariance structure. The exponent of the estimated regression coefficients was calculated to estimate the percent change in the mean response (number of extreme heat events) associated with changes in the covariates. Regression maps were created using ArcGIS 10 (esri, Redlands, CA) to display the county level regression parameter estimate for the impact of ENSO on the number of extreme heat events after adjusting for seasonal and 2006 urban-rural classification.

Results

The final extreme heat event dataset consisted of 3,109 counties over 51 years (1960 to 2010) located in the continental US (Table 1). In general, we observed significantly higher frequency of extreme heat events during the 1990s and 2000s compared to the reference period (1960–1989). This trend was consistent across season, 2006 urban-rural classifications and most Census divisions, with few exceptions. Within the two time periods (1990s and 2000s) the large central metro areas observed higher number extreme heat events compared to small metro and micropolitan areas. We also found an increasing trend in extreme heat events that varied considerably by area of the country, with the most pronounced trend observed for the New England, Middle Atlantic and Mountain divisions with lesser increases in the East and West North Central divisions (Fig 1). Interactions between time periods and Census divisions, ENSO, and Seasons were found to be highly significant and justified the stratification of the analysis by time period.

thumbnail
Fig 1. Temporal trend in extreme heat events across census division for the 1960–2010 periods.

https://doi.org/10.1371/journal.pone.0144202.g001

thumbnail
Table 1. County-level annual frequency of extreme heat events (mean (standard deviation, SD)), excluding Alaska and Hawaii.

https://doi.org/10.1371/journal.pone.0144202.t001

Table 2 provides the frequency of extreme heat events stratified by phases of ENSO for the 3 time periods (1960–1989, 1990–1999, and 2000–2010), across season, 2006 urban rural classification, and Census Division. In general, La Niña periods were characterized by significantly higher frequency of extreme heat events while El Niño periods showed significantly lower frequency of extreme heat events for all seasons, 2006 urban-rural classifications, and most Census divisions when compared to the ENSO Neutral years. A noted exception to this pattern appeared in the 1990–1999 when the frequency of extreme heat events during El Niño were not lower than those observed during the ENSO neutral period, with some census divisions such as New England, Middle Atlantic, East/West North Central, Mountain and Pacific regions observing higher frequency of extreme heat events. Exceptions were also noted for winter of 1990s and 2000s, as well as summer of 1990s.

thumbnail
Table 2. County-level annual frequency of extreme heat events (mean (SD)) overall and by season, urbanization and Census Division, across decades and ENSO periods.

https://doi.org/10.1371/journal.pone.0144202.t002

Table 3 presents the results from three negative binomial GEE models of monthly frequency of extreme heat events stratified by three time periods: 1960–1989, 1990–1999, and 2000–2010. Compared to ENSO neutral periods, El Niño periods were associated with significantly fewer events at the national scale, ranging from 9% fewer (estimated eβ = 0.91, p<0.001) during the 1990s to 24% fewer (estimated eβ = 0.76, p<0.001) during the 2000s, after adjusting for season, 2006 urban-rural classification and Census division (Table 3). By comparison, La Niña periods were associated with as much as 29% higher frequency of extreme heat events at the national level (estimated eβ = 1.29, p<0.001 for 1990–1999 & 2000–2010). For the 1990s and 2000s, counties that were large metropolitan areas based on the 2006 urban-rural classification tended to have a higher frequency of extreme heat events (estimated eβ >1.0) compared to non-core counties; although, this urban-rural difference was statistically significant only during the 1990s (p<0.05). Compared to New England, the other Census divisions of the country had significantly fewer differences in the extreme heat events for each of the three time periods (estimated eβ for all Census divisions <1.0), with the exception of the Middle Atlantic division.

thumbnail
Table 3. Relative percent change in extreme heat events, by time period, for the continental United States, excluding Alaska and Hawaii.

https://doi.org/10.1371/journal.pone.0144202.t003

The analysis for continental US were further broken down by Census division (Table 4). Overall, the Census division results agreed with the country level analysis presented in Table 3 with few noted exceptions. For example, compared to ENSO neutral periods, El Niño years were associated with significantly lower frequency of extreme heat events across Census divisions during 1960–1989 and 2000–2010 period. However, during 1990–1999, El Niño years were associated with increased frequency of extreme heat events compared to ENSO neutral years in several Census divisions (New England, Mid Atlantic and the Pacific divisions). La Niña periods were associated with a higher frequency of extreme heat events than ENSO neutral periods across most Census divisions; with the largest effect (75%) observed for the West South Central division during the 1990–1999 time periods. An exception to this pattern was in the Pacific division, where the La Niña period was associated with 15% lower frequency of extreme heat events compared to the ENSO neutral period during 2000–2010 (estimated eβ = 0.85, p<0.001).

thumbnail
Table 4. Relative percent change in count of extreme heat events, by Census Division (excluding Alaska and Hawaii).

https://doi.org/10.1371/journal.pone.0144202.t004

Spatial heterogeneity in these findings was further investigated using the county level regression coefficients (Fig 2). The findings presented in Fig 2 are in agreement with the results of the divisional model (Table 4), but the finer county level resolution allows for the identification of additional counties whose results were masked in the divisional level analysis (e.g., selected counties in TX and ME experienced larger percent changes in frequency of extreme heat events during El Niño periods). We conducted sensitivity analysis using 1990 NCHS Urban-Rural Classification Scheme for Counties instead of the 2006 schemes; however this did not change our conclusions (results not shown).

thumbnail
Fig 2. Relative percent change in monthly total extreme heat events for La Niña and El Niño months in 1960–2010 compared to ENSO Neutral months, adjusted for seasonal and 2006 land-use classification type.

https://doi.org/10.1371/journal.pone.0144202.g002

Discussion

We developed a generic surrogate exposure metric (extreme heat events) based on climatology that has a broad spatial coverage for the contiguous US with a county level geographic resolution. This exposure metric can correspond to county-level end-points, including many health outcome data such as behavioral risk factor surveillance system (BRFSS), Centers for Medicare and Medicaid Services (CMS), healthcare cost and utilization project (HCUP), CDC’s Public Health Tracking data and others. We assessed the exposure metric for its ability to capture the ENSO events, while controlling for other temporal, seasonal, divisional, and urban-rural classification influences. The results showed the ability of the exposure metric to capture salient features of climate variability and change (long term change in the frequency of extreme heat events) including the effect of natural variability such as El Niño-Southern Oscillation (ENSO) patterns that have distinct heterogeneous effects across geographical regions. We also demonstrated how the exposure metric could provide flexibility in spatial and temporal aggregation of exposure—an ideal attribute for epidemiological studies. Our county level approach enables a straightforward linkage of the exposure metric to many publically available national health outcome data collected at the county level, facilitating investigations of the possible impacts of climate change on chronic health outcomes [2426].

The threshold method we used has been used for defining extremes in studies looking at temperature and precipitation extremes; however, exceedences have not been quantified on a county level for the entire US [2730]. The metric we developed captured the local impacts of ENSO. This oscillation between warm (El Niño) and cold (La Niña) conditions in the equatorial Pacific Ocean can alter weather patterns and latent heat release into the atmosphere. Such changes lead to widespread remodeling in atmospheric circulation patterns far removed from the Pacific Ocean [31]. ENSO events have been linked to droughts, rainfall and the alteration of temperature and sunlight availability across the globe [3234]. In North America, the statistically significant relationships between ENSO and seasonal temperature extremes have occurred mostly in winter [35]. In some divisions and times of year, El Niño and La Niña conditions modify the probabilities of very warm or very cold seasons [35]. The effects of climate change can manifest through natural forcing systems such as ENSO [36]. Globally, ENSO impacts are largely symmetric. The warm state (El Niño) is generally associated with increased precipitation and cooler temperature anomalies and the cold state (La Niña) changes the sign of the anomalies, to a decrease in precipitation and increase in temperature [19]. The known temperature-related impacts of ENSO were expressed in the results of our analysis. Though it should be noted that the surface expressions of the ENSO anomalies in the tropical Pacific are alleged to have changed since 2000 (albeit with similar onsets [37]) and the interconnectedness of these events (teleconnections) over the US appear to have changed [38]. On a continental scale, we found that La Niña periods were consistently associated with an increase in extreme heat events and El Niño periods led to a decrease in extreme heat events. In certain region of the country, the magnitude as well as direction of the associations between ENSO periods and extreme heat events differed from this trend. In this context, the effects of climate change may be local for health endpoints that may manifest via local weather changes [39,40].

Land use factors are an important contributor to divisional climate. Urbanization affects divisional climates through changes in surface energy and water balance. The change in land use can alter the effects of net radiation through the division of energy into sensible and latent heat, and the partitioning of precipitation into soil water, evapotranspiration and runoff [37]. The urban “heat island” effect is an extreme case of how land use modifies divisional climate [41,42]. Previous studies have suggested that a major portion of the reduction in diurnal temperature range observed during the last several decades to urbanization and other land use changes [43,44]. In congruence with available literature, using the 2006 urban-rural classification, we found that more urbanized areas generally experienced relatively high proportional change in extreme heat events compared to the less urbanized areas. However, this pattern was not consistent across Census divisions and was only present during the latter 2 decades. These results may be due, in part, to the classification scheme used in this analysis. This scheme was developed based on the 2006 census statistics and applied in our study for time periods that span more than 4 decades prior. However, sensitivity analyses using the 1990 census scheme produced similar results.

More attention has been paid to the effects of hot temperature anomalies in the summer and spring particularly because these changes can have an impact on biotic factors (e.g., pollen) and industrial air pollution along with heat waves [4549]. Using the metric, we identified larger differences in extreme heat events occurring during the winter, spring and summer months on a continental scale. Yet, at the divisional level, the patterns differ considerably, with the New England and Middle Atlantic divisions experiencing the largest differences in extreme heat events during winter and lowest level during autumn. In the Mountain and Pacific divisions, the largest differences in extreme heat events were observed during spring and lowest level observed during winter season.

Overall, the exposure metric captured subtle variability across geographic division, season, and urban-rural categorization. More importantly, the exposure metric was sensitive to large scale phenomenon such as ENSO that are known to govern local weather patterns. As stated previously, the flexibility of this exposure metric lends itself to epidemiological studies of both infectious and chronic diseases. For example, in a recent study investigating the link between changing climate and Salmonellosis, Jiang et al. (2015) showed that the frequency of extreme heat and precipitation event was directly related to increased risk of Salmonellosis in Maryland, and that the risk was more pronounced among the coastal communities compared to inland communities [50]. Since the precise date of disease onset in the Jiang et al. (2015) was not known, the authors linked monthly count of Salmonellosis with number of extreme heat and precipitation event on the same month and employed negative binomial regression for the statistical analysis. In the instances where the precise date of onset is known (e.g., hospitalization for asthma, or stroke), investigators can use case-crossover approach looking at presence/absence of extreme events in the case period compared to control period with adequate lag structure that are determined based on current knowledge about the disease etiology. In addition, the frequency of extreme heat events can also be used to investigate the spatio-temporal pattern of vector borne diseases (e.g., Lyme disease) that are sensitive to temperature changes. Previous studies have shown that the frequency as well as intensity of extreme events will continue to rise in the near future [51,52]. The exposure metric we have presented in this manuscript allows investigators to document how increases in the frequency of extreme heat event impacts human health.

Conclusion

We report on the development of a novel temperature-related exposure metric and quantify its ability to capture small and large changes in climatic variability across the US and over time. Findings from this study suggest that natural modes of forcing, seasonality, urban-rural classification, and division of country have an impact on the number extreme heat events recorded. We observed that the increases in frequency of extreme heat events differ across the geographical region and time periods. Likewise, we observed higher frequency of extreme heat events during La Niña period and lower frequencies during the El Niño. At regional level, exceptions to this trend were noted for El Niño years in selected geographical areas. This county level exposure metric generated based on location specific climatology data is versatile and can be easily extended to developing metrics for different time periods and county based geographic aggregations. To facilitate research in this area, we will make this exposure metric freely available to potential users through a web portal.

Acknowledgments

The findings and conclusions in this paper are those of the authors and do not necessarily represent the views of the National Center for Health Statistics or the Centers for Disease Control and Prevention. The funding for this work was provided by the Centers for Disease Control and Prevention and the National Institute of Environmental Health Sciences (NIEHS: 1R21ES021422-01A1)

Author Contributions

Conceived and designed the experiments: AS JP RM. Analyzed the data: CRU CJ XH. Wrote the paper: CRU JP AS RM CJ XH.

References

  1. 1. Hales S, Kovats S, Lloyd S, Campbell-Lendrum D. Quantitative risk assessment of the effects of climate change on selected causes of death, 2030s and 2050s. World Health Organization. 2014.
  2. 2. Luber G, Knowlton K, Balbus J, Frumkin H, Hayden M, Hess J, et al. Climate Change Impacts in the United States: The Third National Climate Assessment, Melillo J. M., Terese (T.C.) Richmond, and Yohe G. W., Eds. US Glob Change Res Program. 2014; 220–256.
  3. 3. Haines A, McMichael AJ. Climate change and health: implications for research, monitoring, and policy. BMJ. 1997;315: 870–874. pmid:9353511
  4. 4. McMichael AJ. Global Environmental Change as “Risk Factor”: Can Epidemiology Cope? Am J Public Health. 2001;91: 1172–1174. pmid:11499093
  5. 5. Gutowski WJ, Hegerl GC, Holland GJ, Knutson TR, Mearns LO, Stouffer RJ, et al. Causes of Observed Changes in Extremes and Projections of Future Changes in Weather and Climate Extremes in a Changing Climate. Regions of Focus: North America, Hawaii, Caribbean, and U.S. Pacific Islands. T.R. Karl, G.A. Meehl, C.D. Miller, S.J. Hassol, A.M. Waple, and W.L. Murray (eds.). A Report by the U.S. Climate Change Science Program and the Subcommittee on Global Change Research, Washington, DC. 2008.
  6. 6. Climate Communication | Heat Waves: The Details [Internet]. [cited 20 Aug 2014]. Available: http://www.climatecommunication.org/new/articles/heat-waves-and-climate-change/heat-waves-the-details/.
  7. 7. Perkins SE, Alexander LV. On the Measurement of Heat Waves. J Clim. 2013;26: 4500–4517.
  8. 8. Frumkin H, Hess J, Luber G, Malilay J, McGeehin M. Climate Change: The Public Health Response. Am J Public Health. 2008;98: 435–445. pmid:18235058
  9. 9. Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, et al. Managing the health effects of climate change: Lancet and University College London Institute for Global Health Commission. Lancet. 2009;373: 1693–1733. pmid:19447250
  10. 10. National Climatic Data Center—td3200.pdf [Internet]. Available: http://www1.ncdc.noaa.gov/pub/data/documentlibrary/tddoc/td3200.pdf.
  11. 11. Ingram DD, Franco SJ. NCHS urban-rural classification scheme for counties. Vital Health Stat 2. 2012; 1–65.
  12. 12. Barnett E, Strogatz D, Armstrong D, Wing S. Urbanisation and coronary heart disease mortality among African Americans in the US South. J Epidemiol Community Health. 1996;50: 252–257. pmid:8935454
  13. 13. Ingram DD, Gillum RF. Regional and urbanization differentials in coronary heart disease mortality in the United States, 1968–1985. J Clin Epidemiol. 1989;42: 857–868. pmid:2778467
  14. 14. Knobeloch L, Imm P. Hypertensive Heart Disease Mortality in Wisconsin, 1979–2004. WMJ-MADISON-. 2007;106: 137.
  15. 15. Savitz DA, Stein CR, Ye F, Kellerman L, Silverman M. The Epidemiology of Hospitalized Postpartum Depression in New York State, 1995–2004. Ann Epidemiol. 2011;21: 399–406. pmid:21549277
  16. 16. Savitz DA, Danilack VA, Engel SM, Elston B, Lipkind HS. Descriptive Epidemiology of Chronic Hypertension, Gestational Hypertension, and Preeclampsia in New York State, 1995–2004. Matern Child Health J. 2014;18: 829–838. pmid:23793484
  17. 17. Rossen LM, Khan D, Warner M. Trends and Geographic Patterns in Drug-Poisoning Death Rates in the U.S., 1999–2009. Am J Prev Med. 2013;45: e19–e25. pmid:24237925
  18. 18. Data Access—Urban Rural Classification Scheme for Counties [Internet]. [cited 25 Jun 2014]. Available: http://www.cdc.gov/nchs/data_access/urban_rural.htm.
  19. 19. Climate Prediction Center—El Nino Southern Oscillation: Cold and Warm Episodes by Seasons [Internet]. [cited 29 Oct 2013]. Available: http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml.
  20. 20. Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, et al. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, 2007 [Internet]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA.; 2007. Available: http://www.ipcc.ch/publications_and_data/ar4/wg1/en/contents.html.
  21. 21. Abdi H, Williams LJ. Tukey’s honestly significant difference (HSD) test. Encycl Res Des Thousand Oaks CA Sage. 2010; 1–5.
  22. 22. Greene WH. Accounting for Excess Zeros and Sample Selection in Poisson and Negative Binomial Regression Models [Internet]. Rochester, NY: Social Science Research Network; 1994 Mar. Report No.: ID 1293115. Available: http://papers.ssrn.com/abstract=1293115.
  23. 23. Byers AL, Allore H, Gill TM, Peduzzi PN. Application of negative binomial modeling for discrete outcomes: A case study in aging research. J Clin Epidemiol. 2003;56: 559–564. pmid:12873651
  24. 24. Parker JD, Kravets N, Woodruff TJ. Linkage of the National Health Interview Survey to air quality data. Vital Health Stat 2. 2008; 1.
  25. 25. Parker JD, Akinbami LJ, Woodruff TJ. Air pollution and childhood respiratory allergies in the United States. Environ Health Perspect. 2009;117: 140. pmid:19165401
  26. 26. Akinbami LJ, Lynch CD, Parker JD, Woodruff TJ. The association between childhood asthma prevalence and monitored air pollutants in metropolitan areas, United States, 2001–2004. Environ Res. 2010;110: 294–301. pmid:20117766
  27. 27. Zhai P, Pan X. Change in Extreme Temperature and Precipitation over Northern China During the Second Half of the20th Century. Acta Geogr Sin. 2003;58.
  28. 28. Alexander LV, Zhang X, Peterson TC, Caesar J, Gleason B, Klein Tank AMG, et al. Global observed changes in daily climate extremes of temperature and precipitation. J Geophys Res Atmospheres. 2006;111: D05109.
  29. 29. Beniston M. Trends in joint quantiles of temperature and precipitation in Europe since 1901 and projected for 2100. Geophys Res Lett. 2009;36: L07707.
  30. 30. Hao Z, AghaKouchak A, Phillips TJ. Changes in concurrent monthly precipitation and temperature extremes. Environ Res Lett. 2013;8: 034014.
  31. 31. McPhaden MJ, Zebiak SE, Glantz MH. ENSO as an Integrating Concept in Earth Science. Science. 2006;314: 1740–1745. pmid:17170296
  32. 32. Stenseth NC, Mysterud A, Ottersen G, Hurrell JW, Chan K-S, Lima M. Ecological Effects of Climate Fluctuations. Science. 2002;297: 1292–1296. pmid:12193777
  33. 33. Stenseth NC, Ottersen G, Hurrell JW, Mysterud A, Lima M, Chan K-S, et al. Studying climate effects on ecology through the use of climate indices: the North Atlantic Oscillation, El Nino Southern Oscillation and beyond. Proc R Soc Lond B Biol Sci. 2003;270: 2087–2096.
  34. 34. Wolff C, Haug GH, Timmermann A, Damsté JSS, Brauer A, Sigman DM, et al. Reduced Interannual Rainfall Variability in East Africa During the Last Ice Age. Science. 2011;333: 743–747. pmid:21817050
  35. 35. Wolter K, Dole RM, Smith CA. Short-Term Climate Extremes over the Continental United States and ENSO. Part I: Seasonal Temperatures. J Clim. 1999;12: 3255–3272.
  36. 36. Corti S, Molteni F, Palmer TN. Signature of recent climate change in frequencies of natural atmospheric circulation regimes. Nature. 1999;398: 799–802.
  37. 37. Ramesh N, Murtugudde R. All flavours of El Niño have similar early subsurface origins. Nat Clim Change. 2012;3: 42–46.
  38. 38. Larkin NK, Harrison DE. On the definition of El Niño and associated seasonal average US weather anomalies. Geophys Res Lett. 2005;32: L13705.
  39. 39. Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. Impact of regional climate change on human health. Nature. 2005;438: 310–317. pmid:16292302
  40. 40. Murtugudde R. Regional Earth System prediction: a decision-making tool for sustainability? Curr Opin Environ Sustain. 2009;1: 37–45.
  41. 41. Luber G, McGeehin M. Climate Change and Extreme Heat Events. Am J Prev Med. 2008;35: 429–435. pmid:18929969
  42. 42. Ganeshan M, Murtugudde R, Imhoff ML. A multi-city analysis of the UHI-influence on warm season rainfall. Urban Clim.
  43. 43. Kalnay E, Cai M. Impact of urbanization and land-use change on climate. Nature. 2003;423: 528–531. pmid:12774119
  44. 44. Zhou L, Dickinson RE, Tian Y, Fang J, Li Q, Kaufmann RK, et al. Evidence for a significant urbanization effect on climate in China. Proc Natl Acad Sci U S A. 2004;101: 9540–9544. pmid:15205480
  45. 45. Westerling AL, Hidalgo HG, Cayan DR, Swetnam TW. Warming and earlier spring increase western US forest wildfire activity. Science. 2006;313: 940–943. pmid:16825536
  46. 46. Ziska LH, Beggs PJ. Anthropogenic climate change and allergen exposure: the role of plant biology. J Allergy Clin Immunol. 2012;129: 27–32. pmid:22104602
  47. 47. Menzel A, Fabian P. Growing season extended in Europe. Nature. 1999;397: 659–659.
  48. 48. Cayan DR, Dettinger MD, Kammerdiener SA, Caprio JM, Peterson DH. Changes in the onset of spring in the western United States. Bull Am Meteorol Soc. 2001;82: 399–415.
  49. 49. Leung LR, Gustafson WI. Potential regional climate change and implications to US air quality. Geophys Res Lett. 2005;32.
  50. 50. Jiang C, Shaw KS, Upperman CR, Blythe D, Mitchell C, Murtugudde R, et al. Climate change, extreme events and increased risk of salmonellosis in Maryland, USA: Evidence for coastal vulnerability. Environ Int. 2015;83: 58–62. pmid:26093493
  51. 51. Easterling DR, Meehl GA, Parmesan C, Changnon SA, Karl TR, Mearns LO. Climate extremes: observations, modeling, and impacts. science. 2000;289: 2068–2074. pmid:11000103
  52. 52. McCarthy JJ. Climate change 2001: impacts, adaptation, and vulnerability: contribution of Working Group II to the third assessment report of the Intergovernmental Panel on Climate Change. Cambridge University Press; 2001.