Environ HealthEnviron HealthEnvironmental Health1476-069XBioMed CentralLondon4352547510.1186/s12940-015-0005-zResearchSpatial variation in hyperthermia emergency department visits among those with employer-based insurance in the United States – a case-crossover analysisSahaShubhayu
ssaha@cdc.gov
BrockJohn W
jbrock@warren-wilson.edu
VaidyanathanAmbarish
dvq3@cdc.gov
EasterlingDavid R
david.easterling@noaa.gov
LuberGeorge
gcl4@cdc.gov
4320154320152015142026620141122015© Saha et al.; licensee BioMed Central. 2015This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.Background

Predictions of intense heat waves across the United States will lead to localized health impacts, most of which are preventable. There is a need to better understand the spatial variation in the morbidity impacts associated with extreme heat across the country to prevent such adverse health outcomes.

Methods

Hyperthermia-related emergency department (ED) visits were obtained from the Truven Health MarketScan® Research dataset for 2000-2010. Three measures of daily ambient heat were constructed using meteorological observations from the National Climatic Data Center (maximum temperature, heat index) and the Spatial Synoptic Classification. Using a time-stratified case crossover approach, odds ratio of hyperthermia-related ED visit were estimated for the three different heat measures. Random effects meta-analysis was used to combine the odds ratios for 94 Metropolitan Statistical Areas (MSA) to examine the spatial variation by eight latitude categories and nine U.S. climate regions.

Results

Examination of lags for all three temperature measures showed that the odds ratio of ED visit was statistically significant and highest on the day of the ED visit. For heat waves lasting two or more days, additional statistically significant association was observed when heat index and synoptic classification was used as the temperature measure. These results were insensitive to the inclusion of air pollution measures. On average, the maximum temperature on the day of an ED visit was 93.4oF in ‘South’ and 81.9oF in the ‘Northwest’ climatic regions of United States. The meta-analysis showed higher odds ratios of hyperthermia ED visit in the central and the northern parts of the country compared to the south and southwest.

Conclusion

The results showed spatial variation in average temperature on days of ED visit and odds ratio for hyperthermia ED visits associated with extreme heat across United States. This suggests that heat response plans need to be customized for different regions and the potential role of hyperthermia ED visits in syndromic surveillance for extreme heat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12940-015-0005-z) contains supplementary material, which is available to authorized users.

KeywordsCase-crossoverExtreme heatHyperthermiaMeta-analysisSpatial riskissue-copyright-statement© The Author(s) 2015
Introduction

Recent climate assessments indicate more frequent, more intense, and longer-lasting heat waves for most of the United States [1,2]. Though there is evidence suggesting a declining trend in heat-related mortality in the country [3,4], an aging population and inadequate use of residential air conditioning imply that specific sections of the population will continue to remain vulnerable to the projected increase in extreme heat resulting in a large burden of associated adverse health outcomes [5,6]. While the spatial variation in mortality risk from heat waves is shown to vary across different parts of United States [7], little is known about the spatial variation in risk of morbidity outcomes from extreme heat. Since adverse health impacts from extreme heat are preventable [8], public health agencies may need to consider a range of localized health outcomes for syndromic surveillance as part of designing heat response plans [9-12]. Hyperthermia is a direct physiologic impact from exposure to external heat that impair the thermoregulation mechanisms in the body leading to severe consequences [13]. Syndromic systems based on hyperthermia-related emergency department (ED) visits have been found to be effective in early detection of health impacts during heat waves [14]. During the 2006 heat wave in California, 13% (2134 out of 16166) of the estimated excess ED visits were reported for hyperthermia [15]. A report found that 80% of all hyperthermia-related hospital admissions began in an ED setting, supporting the need to further examine impacts of ambient heat on ED visits [16]. Yet, in the search of the literature on heat and morbidity outcomes [17,18], only one study was found that specifically examined measures of ambient heat and the risk of a hyperthermia-related ED visit [19].

This study therefore examines the associations between measures of ambient heat and hyperthermia-related ED visits during April through September from 2000–2010. Novel patient-level health data was obtained from a health insurance database across multiple U.S. metropolitan statistical areas (MSA). The association between hyperthermia-related ED visit with various measures of ambient heat as found in previous epidemiologic studies was examined – maximum temperature, heat index and Spatial Synoptic Classification [7,19-22], since there is no unanimity on any one measure being superior to others [23]. The additional health impact associated with a heat wave [24] was estimated. A case-crossover design as used in recent studies of heat morbidity [19,21,25] was used to estimate the odds ratio of hyperthermia-related ED visit while controlling for the impact of ambient PM2.5 and ozone concentrations. Finally, the geographic variation in the odds ratio of ED visit was examined by U.S. national climatic regions [26] and latitude zones [27] using random-effects meta-analysis.

MethodsHealth data

For the years 2000–2010, health data from the Truven Health MarketScan® Research database including (i) Commercial Claims and Encounters (CCAE) and (ii) Medicare Supplemental and Coordination of Benefits were obtained. The database is a large convenience sample representative of the US population with employer-based health insurance. The large sample size provides a rare opportunity to examine geographic variation in health effects for low-prevalence conditions like hyperthermia. The database captures de-identified patient-level episode with age, gender, date of healthcare service, county of residence, county of emergency department visited, and a list of diagnoses codes based on the ninth revision of the International Classification of Disease (ICD-9-CM). From all available ED visits in the database, patients were included in the study with (i) an ED visit with a diagnosis of hyperthermia (ICD-9-CM = 992.0-992.9 or injury code = E900.0) anywhere in the record; (ii) episodes occurring between April 1 and September 30 in a calendar year; (iii) same county of residence and the county where the ED was located. Patients with injuries due to excessive heat of man-made origin (ICD-9-CM E900.1) were not included. All five-digit county Federal Information Processing Standard (FIPS) information was converted to a MSA for each hyperthermia-related ED visit. Each of the 141 MSAs for which health data was available had a weather-station located within the jurisdiction.

Measures of ambient heat

Hourly weather data was obtained from the land-based weather stations in these 141 MSAs for each day of the year for 2000–2010 (Integrated Surface Hourly, National Climatic Data Center). Climate Normal information was obtained for each weather station for each day of the year for the period 1980–2010 [28]. Records for weather stations located at major urban airports were used and all stations had less than 5% of days with missing data for the 30 year period. Within each day, missing data were interpolated for up to 4 sequential missing data points by the hour. Days with more than 4 sequential missing data points were excluded from further analyses. The highest recorded temperature among the 24 hourly records (midnight to midnight) was used for daily maximum temperature. The Steadman heat index [29] was used except when temperature or humidity was below 70°F and 40% respectively, the index was set to the dry-bulb air temperature [30]. Daily spatial synoptic classification (SSC) of weather was available for only 107 of these MSAs (http://sheridan.geog.kent.edu/ssc.html). SSC information was available for 90% of the cases. An indicator variable was created to designate a day as extremely hot if the SSC was characterized by the following air mass types - DT,MT+, MT++ for locations in the East of Rockies; DT, MT, MT+, MT++ for locations in the West of Rockies (based on personal communication with Dr. Scott Sheridan). Any day with missing SSC information was not included in the analysis. A heat wave was defined as a block of two or more consecutive days leading to the ED visit when daily maximum temperature or heat index was above the 95th percentile of the 30-year normal for the station [24]. For the SSC, any day that was part of two or more consecutive days with the SSC indicator of extreme heat was deemed to be part of a ‘heat wave’.

Air quality data

PM2.5 daily maximum concentrations (μg/m3) and daily maximum 8-hour average ozone concentrations (ppb) were obtained from the EPA’s Data Mart (http://www.epa.gov/ttn/airs/aqsdatamart/) for monitors designated as Federal Reference Methods (FRM). A daily MSA level dataset was created by retaining the maximum concentration among all monitors within the MSA for each monitored day. Although ozone monitors operate on a daily basis during the summer months, PM monitors typically operate on a once-every-third-day schedule [31]. Since ozone and PM2.5 were simultaneously factored in the analysis, information on both pollutants was available for about 35% of the cases.

The health, temperature and air pollution datasets were them merged by MSA and date. In the final dataset, days that were designated as U. S. federal holidays during the study period were identified.

Statistical analysis

A case-crossover design was used for the analysis [32]. Given the ED visit day for each patient, control days were identified using time-stratified approach [33]. In the time-stratified approach, each month was a priori divided into two halves, (i) days 1–15 and (ii) days 16-end of month, and control day(s) were selected on the same weekday as the case within the same half of the month.

Assessment of odds ratio of ED visit

Conditional logistic regressions (using SAS software ®) was used to assess the odds ratio of hyperthermia-related ED visit and 3 different temperature metrics using the entire dataset. Separate models were estimated for exposures on day0 to day−6 and the smallest Akaike Information Criteria (AIC) was obtained for day0 for all three temperature metrics. In subsequent analysis, day0 temperature metric was used. In a series of regressions, the odds ratio for the measure of ambient heat alone was estimated first (Model A). Then, the heat wave indicator was added to assess if there was any additional heat wave effect (Model B). To model B, measures of air pollution and an indicator for U.S. federal holidays were added (Model C) to see if the association between temperature metric and ED visit was sensitive to those factors. Since estimates of PM2.5 and ozone were missing for some days, the sample size in Models C and D were smaller than those in Models A and B. Finally, Model B was replicated using the smaller dataset with non-missing air pollution measures (Model D) to check for consistency in results between the pooled sample and the sub-sample with complete air pollution information.

Estimation of geographic variation in temperature profile and odds ratio for ED visits

Each MSA was classified into one of the (i) nine US climatic regions and (ii) eight latitude categories spanning the continental U.S. (<30°N, >42°N, and six categories in between in 2° intervals). The average temperature profile (using daily maximum temperature, maximum heat index and SSC) were calculated for the days of ED visit and compared with the ‘control’ days by the climate regions. Some areas of the country are sparsely represented in the health database. Since the estimate of the odds ratio for some MSAs with low number of cases of hyperthermia ED visits could be statistically unreliable, estimates were pooled together geographically to derive more statistically robust results. Random-effects meta-analysis was used to combine the MSA-specific odds ratios by the climate regions and latitude zones. An a priori decision was made to include MSAs which had at least 50 cases of ED visits in the dataset. Comprehensive Meta-analysis Software ® was used to perform this task.

Sensitivity analyses

For selecting control days for the case-crossover analysis, the symmetric bidirectional approach was also used. Following this approach, control days were chosen in a two-week window before and after the case day for the same weekday as the ED visit. The results were comparable to those obtained using the time-stratified approach. An alternative definition of a heat wave - 2 or more consecutive days above the 98th percentile of maximum temperature – yielded comparable results. To check if the results were sensitive to the case definition of hyperthermia, regression results using the sub-sample with only principal diagnosis of hyperthermia produced similar results compared to the entire sample.

ResultsCharacteristics of study population

The dataset comprised of 11,031 episodes of hyperthermia-related ED visits spanning 141 MSAs. The characteristics of the study population are presented in Table 1. A majority (61%) of the patients belonged to the 25–64 years age group, followed by the patients in the 6–18 years group (22%). The age distribution of all individuals included in the MarketScan database for the study years was similar. 68% of the patients were male. 80% of the ED visits occurred in June, July and August. Specifically, 52% of the ED visits within the 13–18 year age group of junior- and high school-aged children were in August and September (Additional file 1: Table S1). 82% of the patients lived either in the South, Southeast, Central or Northeast climate regions. Comparing the percentage of all individuals in the MarketScan database with the hyperthermia ED cases by the climate regions, hyperthermia ED cases were higher from the South and Southeast, and lower from the Northeast and West. A relatively higher percentage of ED visits occurred on Saturdays (19%) compared to the average on other days of the week (~14%).

Descriptive statistics of patients with hyperthermia-related emergency department visit 2000–2010 (N = 11031)

Age N % Month N % Latitude ( o N) N %
<= 5 years1491April3543<30110710
6-12 years5615May10099<=30 & <326666
13-18 years191817June250523<=32 & <34304228
18-24 years8167July308928<=34 & <36151314
25-50 years450241August314529<=36 & <388648
51-64 years222020September9298<=38 & <40157114
65+ years8658<=40 & <42151414
US Climate region*> = 427547
GenderCentral224920
Female350232East North Central5315Weekday
Male752968Northeast123911Sunday154914
Northwest831Monday151514
Hyperthermia DXSouth325329Tuesday158914
Primary926282Southeast245422Wednesday149814
Secondary172716Southwest3403Thursday141913
West8037Friday140213
West North Central791Saturday205919

*Climate regions are comprised of the following states:

Central: KY, IL, IN, MO, OH, TN, WV;

East North Central: IA, MI, MN, WI;

Northeast: CT, DE, ME, MD, MA, NH, NJ, NY, PA, RI, VT;

Northwest: ID, OR, WA;

South: AR, LA, KS, MS, OK, TX;

Southeast: AL, FL, GA, NC, SC, VA;

Southwest: AZ, CO, NM, UT;

West: CA, NV;

West North Central: MT, NE, ND, SD, WY.

Assessment of odds ratio of ED visit

For the different measures of ambient temperature, the temperature measures on day0 were most strongly correlated with the likelihood of ED visit (smallest Akaike Information Criteria across all the models) (Additional file 1: Table S2). In all subsequent analyses, the same day measure of ambient heat was used. These results concurred with findings in other studies where same-day temperature measures were most influential for heat-related morbidity [19].

The odds ratio of ED visit for maximum temperature and heat index remain unchanged with inclusion of other variables (Table 2). For example, a 1°F increase in maximum temperature was associated with an odds ratio of 1.15 (95% confidence interval: 1,14, 1,16) in model A. Inclusion of the heat wave indicator variable in model B or the air pollution variables in model C did not change the odds ratio estimates. In contrast, when the SSC is used, the association attenuates when other covariates are included in the model. Evidence for an additional heat wave impact (using the 95th percentile threshold) is observed when using heat index (odds ratio = 1.46, 95% confidence interval: 1.26, 1.70) and SSC indicator (odds ratio = 1.25, 95% confidence interval: 1.09, 1.44) as the temperature metric. The associations of the air pollution variables remain consistent across the models with a small positive magnitude. Since the estimates for maximum temperature were similar across the model specifications, model A with the complete sample of observations was used to calculate MSA-specific odds ratios.

Regression results for time-stratified case-crossover analyses (pooled data 2000–2010)

Model A Model B Model C Model D
OR 95% CI OR 95% CI OR 95% CI OR 95% CI
Maximum temperature °F1.151.14, 1.161.151.14, 1.161.151.14, 1.171.161.15, 1.17
Heat wave indicatora 1.130.96, 1.331.140.84, 1.551.150.84, 1.54
PM2.5 concentration (g/m3)1.021.01, 1.03
Ozone concentration (ppb)0.990.98, 1.01
Holiday indicatorb 1.230.82, 1.88
Model AIC148331483049724983
Cases110311103137563756
Control days137741377446554655
Maximum heat index1.121.11, 1.211.111.11, 1.121.111.10, 1.121.121.11, 1.13
Heat wave indicatora 1.461.26, 1.701.471.09, 1.981.581.17, 2.12
PM2.5 concentration (g/m3)1.011.01, 1.02
Ozone concentration (ppb)1.011.00, 1.01
Holiday indicator1.250.82, 1.90
Model AIC149021487949625007
Cases110311103137563756
Control days137741377446554655
Spatial synoptic classification2.562.39, 2.742.141.89, 2.431.681.36, 2.081.991.61, 2.45
Heat wave indicatorc 1.251.09, 1.441.321.04, 1.651.321.05, 1.67
PM2.5 concentration (g/m3)1.031.02, 1.04
Ozone concentration (ppb)1.011.01, 1.02
Holiday indicator1.270.85, 1.91
Model AIC149961498653495621
Cases9978997836723672
Control days124421244245374537

Maximum daily temperature, maximum daily heat index and SSC correspond to day of ED visit (day0). Model B indicates if there is an additional ‘heat wave’ duration effect. Since PM2.5 and Ozone concentrations are not available for all days, the sample sizes for Models C & D are smaller. Comparison of Models B & D indicate if the restricted sample used in Model D (observations with air pollution variables) produce different effect estimate than those obtained using the full sample in Model B.

aHeat wave indicator denotes any day that is part of 2 or more consecutive days when maximum temperature was above the 95th percentile for the MSA.

bHoliday indicator denotes a U.S. Federal holiday. cHeat wave indicator denotes any day that is part of 2 or more consecutive days when SSC indicator of extreme heat.

Geographic variation in temperature profile and odds ratio for ED visits

For all three temperature metrics, temperature profiles on the day of the ED visit varied across regions. For example, maximum temperature values on the day of an ED visit were higher in the southern climatic regions compared to the northern regions on average, with a range of 93.4°F in ‘South’ and 81.9°F in the ‘Northwest’ (Table 3). Temperature profiles on days of the ED visit were hotter compared to the ‘control’ days across all climate regions. For example, among the MSAs in the ‘Central’ region, the maximum temperature on the day of an ED visit was 4.9°F higher on average compared to the ‘control’ days. However, this pattern was not uniform across regions, as case days in the north were much hotter compared to the south on average. This spatial variation in temperature profiles on days of ED visit led to the examination of the geographic variation in the odds ratio of hyperthermia-related ED visits. The odds ratio for hyperthermia-related ED visit was estimated for 94 MSAs individually which had at least 50 cases of ED visits each in the dataset (Additional file 1: Table S3). Then random-effects meta-analysis was used to derive estimates for each of the U.S. climatic regions and latitude zones. Test of heterogeneity [34] showed that the odds ratios derived for the different geographic groups of the MSAs were significantly different from each other (I2 statistic = 31.1, P < 0.003). The odds ratios for each climatic region were statistically significant, and varied from 1.07 [95% confidence interval: 1.01, 1.15] in the Southwest region to 1.18 [95% confidence interval: 1.15, 1.21] in the East North Central region (Figure 1). The odds ratios for each latitude category were statistically significant, and varied from 1.11 [95% confidence interval: 1.07, 1.14] in the southernmost MSAs to 1.17 [95% confidence interval: 1.14, 1.19] in the northernmost MSAs (Figure 2). Odds ratio of ED visits were relatively higher in the northern MSAs (within 36°N and 42°N latitudes), compared to those below 36°N. Summary results from the meta-analyses using heat index and SSC produced very similar spatial patterns which are presented in the (Additional file 1: Table S4).

Average temperature profile on days with hyperthermia ED visit (case day) compared with control days by us climate regions

Case day Control day
Climate region Average maximum temperature Average maximum Heat Index % of days indicated as ‘extremely hot’ by SSC Average maximum temperature Average maximum Heat Index % of days indicated as ‘extremely hot’ by SSC
Central88.593.33783.688.717
East North Central85.489.83978.884.411
Northeast87.490.64380.685.412
Northwest81.981.24371.479.313
South93.497.03891.395.229
Southeast90.495.13287.792.419
Southwest90.684.67089.084.166
West87.281.84682.380.730
West North Central89.893.43983.686.117

Results from the random effect meta-analysis of odds ratios of hyperthermia-related ED visit associated with maximum temperature grouped by US climate regions. Any MSA with less than 50 observations were excluded. The number in the boxes on top show the number of MSAs included in the meta-analysis for each latitude category. Hyperthermia data on ED visit was obtained from MarketScan research database for 2000–2010.

Results from the random effect meta-analysis of odds ratios of hyperthermia-related ED visit associated with maximum temperature grouped by latitude categories. Any MSA with less than 50 observations were excluded. The number in the boxes on top show the number of MSAs included in the meta-analysis for each latitude category. Hyperthermia data on ED visit was obtained from MarketScan research database for 2000–2010.

Discussion

A unique patient-level insurance claims database was combined with different measures of daily ambient heat and analyzed using a case-crossover study design to estimate odds ratio of hyperthermia-related ED visit. Our analyses show a gradient in the impacts of heat across the continental United States - higher odds ratios of hyperthermia ED visits were estimated in the central and the northern parts of the country compared to the south and southwest. This finding aligns with a previous study on heat mortality in 48 cities in the United States [7] that found the risk of mortality to be higher in the Northeast and Midwest compared to the South. Patients in the northern parts of the country visit the ED with a hyperthermia diagnosis on days with lower daily maximum temperatures than in southern parts on average. People living in higher latitudes have previously been found to be more susceptible to extreme high temperatures [11]. Further, our results suggest that patients in the northern parts may be more sensitive to unusually hot days, as the temperature difference between ED visit days and control days was higher on average in northern parts compared to southern parts. The potential of greater tolerance and adaptability among people in the southern regions of US is also mentioned in analysis of heat-related 911 emergency dispatches in Phoenix and Chicago [35]. While different heat metrics produced different magnitudes of odds ratios for hyperthermia ED visits, the spatial pattern in regional estimates of these odds ratios were similar. The range of temperature values associated with mortality across different parts of the US [7] however were different compared to the temperature profile we find on days of hyperthermia ED visits in our data.

Our finding of the strongest association of ambient heat on the day of an ED visit and a significant association up to 3 days prior is similar to that found in previous studies [18]. While several studies have examined the impact of extreme temperature on morbidity, few specifically studied hyperthermia as an outcome. The risk estimates of hyperthermia ED visits from extreme heat found for California [19] (% excess risk per 10°F was 393.3) is similar to what was found for the Western states of California and Nevada in this study (% excess risk for 10°F was 210.5). Additionally, these results suggest that middle and high school aged youth have a higher frequency of hyperthermia-related ED visits in August and September. The literature supports these findings with evidence on increased risks among adolescents playing outdoor sports and thus active youth are deemed to be among the at-risk population during heat waves [36,37]. The association between heat and ED visit was not found to be sensitive to inclusion of air pollution factors in the analysis. While PM2.5 was found to have no effect on the association between apparent temperature and emergency department visit [19], the effect of heat wave on mortality was found to be higher on high ozone and PM10 days [38]. In a review article on morbidity impacts of heat [18], little agreement was found on potential confounding or modification role of air pollution on heat-related morbidity.

The study has several limitations. The health data used in this study include only individuals with employer-based health insurance or Medicaid living in urban areas, thus limiting the generalizability of the results to the broader US population. As the health risk from excess heat is considerable among the elderly, the uninsured [16] and those living in rural regions [39], the regional estimates derived in this analysis could represent a lower bound. Using the online query system provided by the Healthcare Cost and Utilization Project (http://hcupnet.ahrq.gov/), an average of 50% of patients with hyperthermia-related ED visit in the US between 2006 and 2010 were found to have had health insurance through private sources or Medicaid. MarketScan also slightly over represents the ‘Central’ and ‘Southern’ climatic regions and under represents the ‘Northeast’ as compared to the distribution of the overall population. The fact that 71% of the hyperthermia ED visits used in this study were from the ‘Central’, ‘South’ and ‘Southeast’ climatic regions could be a result of this artifact of the database. However, national rates of hyperthermia-related ED visits were found to be highest in the South and Midwest regions [40] in agreement with the finding in this study. While previous studies have examined the impact of extreme heat on a broader suite of diseases [41], this analysis only examines hyperthermia outcome. However, hyperthermia is a direct physiologic manifestation when exposed to high ambient heat and often a part of suite of health conditions included in a heat surveillance system. In the absence of information on individual level exposure to ambient heat for each patient on each day included in the analysis, one of the criteria to be included in the analysis was that the county of patient’s residence and the ED that was visited needed to be the same. This was the only alternative to ensure that the patient was in the same location on the days included in the analysis, and station-based temperature data for each MSA were most likely what they were exposed to. This study only models the linear effect of same day temperature on the health outcome. While this helped establish the utility of health insurance databases in studying health effects of exposure to heat by producing results comparable with results obtained using commonly available hospital records [19], future applications of this dataset will explore the potential non-linear effect of cumulative temperature exposure over a period of time [42].

Conclusions

Increase in heat-related illness associated with rising temperatures in the future remains a public health concern. A study estimated a two to six-fold increase in excess respiratory illness during warmer summers in 2080–99 compared to 1991–2004 in the state of New York [22]. Establishing syndromic surveillance systems for extreme heat using hyperthermia-related ED visit information has proved to be an effective response strategy [14]. This study shows hyperthermia ED visits to be sensitive to high temperatures across different parts of country, and availability of real-time ED visit information on hyperthermia could be considered an important indicator for heat surveillance systems. Since we find that hyperthermia ED visits were sensitive to extreme heat and the associated temperature profiles were different across regions of the US, public health agencies would need to assess these risks analyzing locally available data and use those findings in the design of comprehensive heat response plans. The Building Resilience Against Climate Effects (BRACE) framework developed by the Centers for Disease Control for public health agencies to prepare for climate change suggest derivation of exposure-outcome relationships to estimate the disease burden for climate-sensitive environmental exposures [43]. The regional exposure-outcome associations derived for hyperthermia add to the available suite of locally-specific estimates to assess heat-related morbidity.

Additional file

Distribution of hyperthermia-related cases by age group and month. Table S2. Results from separate regression models with different temperature metric and different temporal lag. Table S3. MSA-specific odds ratio estimates for hyperthermia-related ED visit and maximum temperature on day of ED visit. Table S4. Summary of odds ratios for three temperature metrics and hyperthermia ED visit by US climatic regions.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors have contributed meaningfully to the development of the manuscript, have read the manuscript, and accept the responsibility for the manuscript’s content. SS and JWB conceived the projected, performed data analysis and drafted the article. AV helped with acquiring the air pollution, weather data and drafting the article. DRE helped with the acquisition and analysis of weather data. GL assisted in drafting and analyzing the health data. All authors read and approved the final manuscript.

Acknowledgements

The authors would like to thank Dr. Dana Flanders, Ethel Taylor, Dr. Fuyuen Yip at Centers for Disease Control for helpful comments on the paper. John W. Brock acknowledges funding from the Climate and Health Program within the Centers for Disease Control and Prevention and sabbatical release from Warren Wilson College for research. We thank Dr. Jay Livermore at National Oceanic and Atmospheric Administration’s National Climatic Data Center (NCDC) for consultation, office space, use of computer time and access to archived weather data. We thank Anthony Arguez for assistance with the U.S. climate normals. We also thank Steven Anthony from STG Inc. for technical support at NCDC.

ReferencesKunkel, K.E, L.E. Stevens, S.E. Stevens, L. Sun, E. Janssen, D. Wuebbles, K.T. Redmond and J.G. Dobson, 2013. Regional Climate Trends and Scenarios for the U.S. National Climate Assessment. Part 5. Climate of the Southwest U.S., NOAA Technical Report NESDIS 142-5, 79 pp. http://scenarios.globalchange.gov/report/regional-climate-trends-and-scenarios-us-national-climate-assessment-part-9-climate.MeehlGATebaldiCWaltonGEasterlingDMcDanielLRelative increase of record high maximum temperatures compared to record low minimum temperatures in the U.SGeophys Res Lett20093623L2370110.1029/2009GL040736DavisREKnappenbergerPCMichaelsPJNovicoffWMChanging heat-related mortality in the United StatesEnviron Health Perspect2003111141712171810.1289/ehp.633614594620PetkovaEPGasparriniAKinneyPLHeat and mortality in New York City since the beginning of the 20th centuryEpidemiology (Cambridge, Mass)201425455456010.1097/EDE.0000000000000123BobbJFPengRDBellMLDominiciFHeat-related mortality and adaptation to heat in the United StatesEnviron Health Perspect201412281181624780880SheridanSKalksteinAKalksteinLTrends in heat-related mortality in the United States, 1975–2004Nat Hazards200950114516010.1007/s11069-008-9327-2AndersonBGBellMLHeat waves in the United States: mortality risk during heat waves and effect modification by heat wave characteristics in 43 U.S. communitiesEnviron Health Perspect2011119221021810.1289/ehp.100231321084239LuberGMcGeehinMClimate change and extreme heat eventsAm J Prev Med200835542943510.1016/j.amepre.2008.08.02118929969BernardSMMcGeehinMAMunicipal heat wave response plansAm J Public Health20049491520152210.2105/AJPH.94.9.152015333307HajatSSheridanSCAllenMJPascalMLaaidiKYagoutiAHeat-health warning systems: a comparison of the predictive capacity of different approaches to identifying dangerously hot daysAm J Public Health201010061137114410.2105/AJPH.2009.16974820395585MartielloMAGiacchiMVHigh temperatures and health outcomes: a review of the literatureScandinavian Journal of Public Health20103888263710.1177/140349481037768520688791SemenzaJCHallDEWilsonDJBontempoBDSailorDJGeorgeLAPublic perception of climate change voluntary mitigation and barriers to behavior changeAm J Prev Med200835547948710.1016/j.amepre.2008.08.02018929974BouchamaAKnochelJPHeat strokeN Engl J Med2002346251978198810.1056/NEJMra01108912075060JosseranLFouilletACaillèreNBrun-NeyDIlefDBruckerGAssessment of a syndromic surveillance system based on morbidity data: results from the Oscour® Network during a heat wavePLoS One201058e1198410.1371/journal.pone.001198420711252KnowltonKRotkin-EllmanMKingGMargolisHGSmithDSolomonGThe 2006 California heat wave: impacts on hospitalizations and emergency department visitsEnviron Health Perspect20091171616710.1289/ehp.1159419165388Merrill CT, Miller M, Steiner C. Hospital Stays Resulting from Excessive Heat and Cold Exposure Due to Weather Conditions in U.S. Community Hospitals, 2005. In: Statistical Brief #55. Healthcare Cost and Utilization Project; 2008. http://www.hcup-us.ahrq.gov/reports/statbriefs/sb55.pdf.Oudin ÅströmDBertilFJoacimRHeat wave impact on morbidity and mortality in the elderly population: a review of recent studiesMaturitas20116929910510.1016/j.maturitas.2011.03.00821477954YeXWolffRYuWVaneckovaPPanXTongSAmbient temperature and morbidity: a review of epidemiological evidenceEnviron Health Perspect20121201192810.1289/ehp.100319821824855BasuRPearsonDMaligBBroadwinRGreenRThe effect of high ambient temperature on emergency room visitsEpidemiology (Cambridge, Mass)201223681382010.1097/EDE.0b013e31826b7f97AndersonGBDominiciFWangYMcCormackMCBellMLPengRDHeat-related emergency hospitalizations for respiratory diseases in the Medicare populationAm J Respir Crit Care Med2013187101098110310.1164/rccm.201211-1969OC23491405FletcherBALinSFitzgeraldEFHwangSAAssociation of summer temperatures with hospital admissions for renal diseases in New York State: a case-crossover studyAm J Epidemiol2012175990791610.1093/aje/kwr41722455834LinSHsuWHVan ZutphenARSahaSLuberGHwangSAExcessive heat and respiratory hospitalizations in New York State: estimating current and future public health burden related to climate changeEnviron Health Perspect2012120111571157710.1289/ehp.110472822922791BarnettAGTongSClementsACWhat measure of temperature is the best predictor of mortality?Environ Res2010110660461110.1016/j.envres.2010.05.00620519131HajatSArmstrongBBacciniMBiggeriABisantiLRussoAImpact of high temperatures on mortality: is there an added heat wave effect?Epidemiology (Cambridge, Mass)200617663263810.1097/01.ede.0000239688.70829.63BhaskaranKArmstrongBHajatSHainesAWilkinsonPSmeethLHeat and risk of myocardial infarction: hourly level case-crossover analysis of MINAP databaseBMJ (Clinical research ed)2012345e8050KarlTRKnightRWSecular trends of precipitation amount, frequency, and intensity in the United StatesBull Am Meteorol Soc199879223124110.1175/1520-0477(1998)079<0231:STOPAF>2.0.CO;2Xiao J, Peng J, Zhang Y, Liu T, Rutherford S, Lin H et al. How much does latitude modify temperature-mortality relationship in 13 eastern US cities? International journal of biometeorology 2015;59(3):365–372. http://link.springer.com/article/10.1007%2Fs00484-014-0848-y.ArguezADurreIApplequistSVoseRSSquiresMFYinXNOAA's 1981–2010 U.S. Climate normals: an overviewBull Am Meteorol Soc201293111687169710.1175/BAMS-D-11-00197.1Rothfusz LP: The Heat Index “Equation”. Technical attachment, SR 90–23. In. Edited by Scientific Services Division NSRH. Fort Worth, TX; 1990. http://www.srh.noaa.gov/images/ffc/pdf/ta_htindx.PDF.RobinsonPJOn the definition of a heat waveJ Appl Meteorol20014076277510.1175/1520-0450(2001)040<0762:OTDOAH>2.0.CO;2VaidyanathanADimmickWKeglerSQualtersJStatistical air quality predictions for public health surveillance: evaluation and generation of county level metrics of PM2.5 for the environmental public health tracking networkInt J Health Geogr201312111310.1186/1476-072X-12-1223305074Carracedo-MartinezETaracidoMTobiasASaezMFigueirasACase-crossover analysis of air pollution health effects: a systematic review of methodology and applicationEnviron Health Perspect201011881173118210.1289/ehp.090148520356818JanesHSheppardLLumleyTCase-crossover analyses of air pollution exposure data: referent selection strategies and their implications for biasEpidemiology (Cambridge, Mass)200516671772610.1097/01.ede.0000181315.18836.9dHigginsJPTThompsonSGQuantifying heterogeneity in a metaanalysisStat Med2002211539155810.1002/sim.118612111919HartzDABrazelAJGoldenJSA comparative climate analysis of heat-related emergency 911 dispatches: Chicago, Illinois and Phoenix, Arizona USA 2003 to 2006Int J Biometeorol20125756697810.1007/s00484-012-0593-z23053064GilchristJHaileyesusTMurphyMYardENonfatal sports and recreation heat illness treated in hospital emergency departments–United States, 2001–2009MMWR Morb Mortal Wkly Rep2011602997798021796094KerrZYCasaDJMarshallSWComstockRDEpidemiology of exertional heat illness among U.S. high school athletesAm J Prev Med201344181410.1016/j.amepre.2012.09.05823253644AnalitisAMichelozziPD’IppolitiDDe’DonatoFMenneBMatthiesFEffects of heat waves on mortality: effect modification and confounding by air pollutantsEpidemiology (Cambridge, Mass)2014251152210.1097/EDE.0b013e31828ac01bSheridanSCDolneyTJHeat, mortality, and level of urbanization: measuring vulnerability across Ohio, USAClimate Res20032425526510.3354/cr024255HessJJSahaSLuberGSummertime acute heat illness in U.S. Emergency Departments from 2006 through 2010: analysis of a nationally representative sampleEnviron Health Perspect2014122111209121524937159LinSLuoMWalkerRJLiuXHwangSAChineryRExtreme high temperatures and hospital admissions for respiratory and cardiovascular diseasesEpidemiology (Cambridge, Mass)200920573874610.1097/EDE.0b013e3181ad5522GasparriniAModeling exposure-lag-response associations with distributed lag non-linear modelsStat Med201433588189910.1002/sim.596324027094MarinucciGLuberGUejioCSahaSHessJBuilding resilience against climate effects—a novel framework to facilitate climate readiness in Public Health AgenciesInt J Environ Res Public Health20141166433645810.3390/ijerph11060643324991665