Skip to main content

Time-series analysis of daily ambient temperature and emergency department visits in five US cities with a comparison of exposure metrics derived from 1-km meteorology products



Ambient temperature observations from single monitoring stations (usually located at the major international airport serving a city) are routinely used to estimate heat exposures in epidemiologic studies. This method of exposure assessment does not account for potential spatial variability in ambient temperature. In environmental health research, there is increasing interest in utilizing spatially-resolved exposure estimates to minimize exposure measurement error.


We conducted time-series analyses to investigate short-term associations between daily temperature metrics and emergency department (ED) visits for well-established heat-related morbidities in five US cities that represent different climatic regions: Atlanta, Los Angeles, Phoenix, Salt Lake City, and San Francisco. In addition to airport monitoring stations, we derived several exposure estimates for each city using a national meteorology data product (Daymet) available at 1 km spatial resolution.


Across cities, we found positive associations between same-day temperature (maximum or minimum) and ED visits for heat-sensitive outcomes, including acute renal injury and fluid and electrolyte imbalance. We also found that exposure assessment methods accounting for spatial variability in temperature and at-risk population size often resulted in stronger relative risk estimates compared to the use of observations at airports. This pattern was most apparent when examining daily minimum temperature and in cities where the major airport is located further away from the urban center.


Epidemiologic studies based on single monitoring stations may underestimate the effect of temperature on morbidity when the station is less representative of the exposure of the at-risk population.

Peer Review reports


Studies of the relationship between ambient temperature and adverse health outcomes routinely utilize observations at single airport weather stations to derive exposure estimates. While airport stations have long time series and high-quality data, their measurements may not accurately reflect average exposure experienced by the surrounding study population. Specifically, airport monitors cannot capture intra-urban spatial variability driven by features of the built and natural environment [1, 2], especially when there exist urban-rural gradients due to urban heat island effects [3, 4].

Exposure measurement error that arises from insufficient characterization of exposure spatial variability may lead to attenuated health effect estimates [5, 6]. Previous investigations in air pollution epidemiology have demonstrated that using more spatially-resolved exposure estimates can result in stronger associations compared to using a single air quality monitor in the study region [7,8,9]. Accurate exposure assessment is an important component of population-based health studies, in part because the health association estimates are frequently used in subsequent risk assessment and health impact analyses [10, 11]. While the use of spatially-resolved exposures has become common in air pollution epidemiology, only a few studies on temperature and mortality have considered heat exposure’s spatial variation in short-term health effect studies [12,13,14]. Though spatially-resolved estimates have been utilized in studies that focus on spatial exposure contrasts [15,16,17,18].

Our study has two objectives. The first is to evaluate the potential benefits of using fine-scale meteorology data products when estimating acute health effects of high temperature, as opposed to a single monitoring station. We are particularly interested in comparing health effect estimates obtained from exposure assessment methods that account for spatial variability and those that assume spatially homogenous exposures. The second objective of this paper is to fill an important knowledge gap on the association between high temperature and morbidity as measured by emergency department (ED) visits. Previous work has predominantly focused on mortality [19,20,21], and hospitalizations [22,23,24,25], and many of the multi-city US morbidity studies have been restricted to the Medicare population ages 65 or above.

To achieve these objectives, we conducted time-series analyses of warm-season daily temperature (maximum and minimum) and all-age ED visits in five US cities from different climatic regions: Atlanta, Los Angeles, Phoenix, Salt Lake City, and San Francisco. To estimate exposures, we used Daymet [26], a publicly available product that provides gridded daily meteorology estimates at 1 km spatial resolution starting in 1980. Daymet data have been found to accurately describe ambient temperature and mean heat index at weather stations across the contiguous United States [27]. Daymet also provides complete spatial coverage over the United States such that spatial variability in at-risk population size can also be incorporated when analyzing aggregated health outcomes [28]. To our knowledge, this is the first study to explore the impact of spatially-resolved exposure variables in estimating short-term temperature-morbidity associations.

Materials and methods

Emergency department visit and meteorological data

Multi-year, patient-level ED records were obtained from individual hospitals in the Atlanta metropolitan area (1993–2004) and the Georgia Hospital Association (2005–2012), the Arizona Department of Health Services, Bureau of Public Health Statistics (2008–2016), the California Health and Human Services Agency, Office of Statewide Planning and Development; 2005–2016), and the Utah Department of Health, Office of Health Care Statistics (2005–2016). Records included admission date and International Classification of Diseases (ICD) diagnosis codes. ICD 9th revision (ICD-9) codes were used for ED visits prior to October 1, 2015, followed by the use of ICD 10th revision (ICD-10) codes. Our definition of an ED visit included patients seen in the ED as outpatients and discharged, as well as patients admitted as inpatients from the ED. The analysis was restricted to the warmest 6 months of the year, May to October.

We used both the primary and secondary diagnosis codes to identify ED visits for specific health outcomes that have been found to be associated with high temperature in previous studies [29]. The health outcomes of interest for this study were fluid and electrolyte imbalance (ICD-9: 276, ICD-10: E86-E87), acute renal injury (ICD-9: 584, ICD-10: N17), circulatory disease (ICD-9: 390–459, ICD-10: I00-I99), respiratory disease (ICD-9: 460–519, ICD-10: J00-J99), gastrointestinal infections (ICD-9: 001–009, ICD-10: A00-A09), and heat-related illnesses (ICD-9: 992, ICD-10: T67). For each of these outcomes, ED visits were aggregated over the metropolitan statistical area (MSA) by admission date.

Analytic datasets of daily ED visit counts for each city were obtained by aggregating ED visits by day based on patient residential location in each corresponding MSA. Specifically, an ED visit was included in the study if the patient’s residential ZIP code corresponded to a ZIP code tabulation area (ZCTA) that overlapped with the MSA definition for a city, defined by the U.S. Office of Management and Budget based on contiguous counties. The five MSAs varied in the number of counties (and area size): 20 counties (15,013 km2) for Atlanta, 2 counties (11,809 km2) for Los Angeles, 2 counties (35,406 km2) for Phoenix, 3 counties (23,577 km2) for Salt Lake City, and 5 counties (6136 km2) for San Francisco. Based on the Köppen-Geiger climate classification [30], the five cities are classified as: Atlanta (humid subtropical), Los Angeles (Mediterranean hot summer), Phoenix (acrid climate), San Francisco (Mediterranean warm/cool summer), Salt Lake City (hot summer continental).

Daily maximum temperature, minimum temperature, and average dewpoint temperature observations were obtained from monitors at the major international airport serving each MSA. We also obtained Daymet ( gridded surfaces and identified all 1 km × 1 km grid cells that are within the MSA boundary. The Daymet temperature product ingests all surface observations from the Global Historical Climatology Network. The algorithm also performs spatial-temporal interpolation using surrounding observations, as well as information on other weather variables, elevation, and land/water masks. A Daymet grid cell is linked to an MSA if the centroid of the 1 km grid cell falls within any of the MSA’s county boundary. Daymet temperature data were used to develop three different daily exposure temperature metrics: a daily simple spatial average of all 1 km grid cells over each MSA, a daily weighted average based on county population, and a daily weighted average based on ZCTA population. For time-series analysis of aggregated health outcomes, the optimal exposure metric should correspond to the average of exposures of all at-risk individuals. Data on population size were obtained from the US Census at the county-level (1990, 2000, and 2010) and at the ZCTA-level (2000 and 2010). Annual populations were based on linear interpolation of decadal Census data. Daily minimum and maximum temperature from the different exposure assessment methods were also averaged to obtain daily mean temperature.

Statistical analysis

To estimate the association between temperature and ED visits during warm seasons in each city, we used an over-dispersed Poisson log-linear model. The primary analysis focuses on the associations between same-day (lag 0) temperature and ED visit outcome counts, modeled using natural cubic splines with 4 degrees of freedom (3 equidistant interior knots). While mortality studies have usually utilized longer lags (e.g., 7 days), previous analyses have found that acute effects of heat on ED visits operate on shorter lags [29, 31]. We fitted separate models for each city, each ED visit outcome, each exposure variable (maximum, average or minimum temperature), and each exposure assessment method (airport observations, Daymet simple average, Daymet county population weighted average (PWA), and Daymet ZCTA PWA).

All time-series models included the following variables for confounder adjustments. We controlled for daily mean dew-point temperature using natural cubic splines with 4 degrees of freedom. Long-term and seasonal temporal trends were modeled as a smooth function of day-of-year with 6 degrees of freedom. We added interaction terms between time spline coefficients and indicators for year, resulting in year-specific seasonal trends. For heat-related illnesses, we only included interactions between linear day-of-year and indicators of year due to the smaller number of events. We also included indicators for day-of-the-week, federal holidays, and hospital-specific indicators to account for hospitals’ contributions to the total ED visits in the city.

To better compare the strength of non-linear associations between ED counts and temperature across different exposure assessment methods, we report relative risk based on two different contrasts along the exposure range: 75th percentile versus 25th percentile, and 95th percentile versus 50th percentile. The percentile values were specific to each city and exposure metric. Some studies have defined a common reference temperature by identifying the temperature value that corresponds to the minimum risk on the non-linear function [32]. We opted for a percentile-based method because of the restriction to the warm season and because the range of exposure varied across exposure methods.

It is difficult to statistically assess differences in relative risks estimated from to the use of different exposure assessment methods because these models were fit to the same health outcome and were non-nested. In addition to qualitatively examining the magnitude of associations, we considered two model comparison tools: (1) the Akaike information criterion (AIC) assuming the outcome is Poisson and (2) the estimated overdispersion parameter. Lower values of AIC and overdispersion are preferred, and may indicate better model fit and out-of-sample prediction performance. We note that these are general model comparison tools that do not necessarily reflect the accuracy and precision in inference for the relative risk parameters of interest. We conducted three sensitivity analyses to examine the robustness of the estimated associations and comparison across exposure assessment methods. First, we varied the natural cubic spline’s degrees of freedom for the exposure of interest to 3, 5, or 6. Second, we evaluated an additional exposure assessment method by considering only the Daymet grid cell linked to the airport monitor. Third, we modified the temperature exposure metric to a 3-day moving average (same-day, lag-1, lag-2).


Descriptive statistics (mean and standard deviation) of the four temperature metrics are given in Table 1. Study locations showed clear geographical differences in temperature. For all exposure metrics, San Francisco and Los Angeles had lower temperatures, while Phoenix had considerably higher temperatures. Supplementary Figure S1 shows the spatial variability in average Daymet maximum and minimum temperature (May to October) at 1 km resolution and locations of airport monitors. We assessed spatial variability by computing the standard deviation of Daymet values within each MSA on each day. Across the study period, the average standard deviation for maximum (minimum) temperature were: 0.33 (0.27) for Atlanta, 1.94 (0.99) for Los Angeles, 2.19 (1.76) for Phoenix, 3.29 (2.13) for Salt Lake City, and 0.79 (0.45) for San Francisco. Greater within-MSA variability was observed for cities with larger area and with more heterogeneous topography.

Table 1 Mean (standard deviation) in °C of daily maximum (TMX), mean (AVG) and minimum (TMN) temperature during May to October by study city for four exposure metrics: airport observations, simple average of Daymet 1 km data, and Daymet population-weighted averages (PWA) using county or ZIP code population

Pairwise correlations between the four temperature metrics are given in Supplementary Table S1. The three Daymet metrics were highly correlated (≥ 0.90) with airport observations, except for Los Angeles and San Francisco. For example, correlations between Daymet average and airport observations for minimum temperature were 0.93 for Atlanta, 0.87 for Los Angeles, 0.95 for Phoenix, 0.97 for Salt Lake City, and 0.79 for San Francisco. Airport monitors in Los Angeles and San Francisco are located by the water (Figure S1) and this likely contributed to their weaker correlations with other exposure metrics [33]. We also calculated correlations between on the subset of days when airport monitor observations were above its the 75th percentile. The correlations were consistently lower when restricted to the highest quartiles of temperature, especially for daily maximum temperature in Los Angeles.

Outcome-specific total and mean daily ED visits are shown in Table 2. Relative risks (RRs) of daily ED visits associated with same-day temperature between the 95th and the 50th percentiles are shown in Fig. 1 for minimum temperature and in Fig. 2 for maximum temperature. Overall, consistent positive associations were identified between ambient temperature and ED visits. We also found evidence that exposure assessment methods accounting for spatial variability in temperature and population are associated with stronger RR estimates compared to the use of airport monitor observations. This pattern is most apparent for Los Angeles and San Francisco, for outcomes that are particularly heat-sensitive (i.e., fluid and electrolyte imbalance, heat-related illnesses, and acute kidney injury), and for minimum temperature. Similar patterns were found for RRs between the 75th and 25th percentile of the exposure (Supplementary Figures S2 and S3). RRs for daily mean temperature (Supplemental Figures S4 and S5) are similar to RRs or daily maximum temperature. We continue to see stronger associations with the use of Daymet data projects in Los Angeles and San Francisco, and for fluid and electrolyte imbalance and heat-related illnesses.

Table 2 Total counts of emergency department visits during May to October in five US cities. Daily average counts are provided in parentheses
Fig. 1
figure 1

Relative risks (RR) of daily emergency department visits associated with same-day minimum temperature between the 95th and the 50th percentile, comparing four different exposure assessment methods: airport observations (), Daymet average (), Daymet county-level population-weighted average (), and Daymet ZCTA population-weighted average (▲). The y-axis ranges are different across outcomes

Fig. 2
figure 2

Relative risks (RR) of daily emergency department visits associated with same-day maximum temperature (Max) between the 95th and the 50th percentile, comparing four different exposure assessment methods: airport observation (), average of Daymet data (), county-level population-weighted average (), and ZIP code-level population-weighted average (▲). The y-axis ranges are different across outcomes

To better visualize differences in the entire exposure-response function due to exposure assessment methods, Fig. 3 shows the relative risk for two outcomes (fluid and electrolyte imbalance, and acute kidney injury) in Los Angeles and San Francisco, comparing the airport observations and Daymet ZCTA PWA exposures. The exposure-response function is centered at the 50th percentile of each exposure metric. We see evidence that accounting spatial variability resulted in steeper changes in relative risks across the entire exposure distribution in the cities with the largest differences in temperature distribution when comparing their Daymet exposures and airport observations. Differences in AIC and ratios of overdispersion for each city, outcome, and temperature metric are given in Supplementary Table S3. In general, we see evidence that the use of Daymet-derived exposures resulted in smaller AIC and overdispersion in Los Angeles and San Francisco. This suggests that the differences observed in RRs and exposure-response relationships lead to better model fit.

Fig. 3
figure 3

Exposure-response functions and 95% pointwise confidence intervals for warm-season daily minimum temperature and same-day emergency department visits for fluid and electrolyte imbalance and acute kidney injury, comparing two different exposure assessment methods: airport observations and Daymet ZCTA population-weighted average (PWA). Reference level for the relative risk (RR) is the median observed temperature value for each exposure metric

In sensitivity analyses, when considering using Daymet values at the airport monitors, we found that estimated RRs were either similar to those using airport observations (e.g., Los Angeles and San Francisco) or similar to those using population-weighted Daymet data (Supplementary Figures S4 and S5). The estimated RRs across exposure metrics were also robust against changing the degrees of freedom for the exposure-response function from 4 to 3, 5, or 6 (results for Daymet ZCTA PWA exposures given in Supplementary Figures S6 and S7). RRs of daily ED visits associated with 3-day moving-average temperature between the 95th and the 50th are given in Supplementary Figures S8 and S9. Generally, we found that the use of moving average temperature data reduced differences between using airport observations and Daymet ZCTA PWA exposures. One explanation may be that additional temporal smoothing increased the correlation between monitor measurements and Daymet-derived exposures. However, differences in relative risk estimates across exposure metrics were still more pronounced for minimum temperature and in San Francisco and Los Angeles. Numerical values of all estimated relative risks and 95% confidence intervals for all primary analyses are provided in the Supplementary Materials.


We conducted time-series analyses to examine associations between same-day temperature and ED visits in five US cities: Atlanta, Los Angeles, Phoenix, Salt Lake City and San Francisco. In addition to airport monitoring stations, we utilized the Daymet maximum and minimum temperature products at 1 km spatial resolution. Overall, we found consistent positive associations between ambient temperature and ED visits for various outcomes across cities, though associations were often stronger when using Daymet-based exposures as compared to those from the airport stations. Furthermore, incorporating spatial distributions of the at-risk population may provide a more accurate measure of the temperature experienced by the populations of each of the cities.

A few recent temperature and mortality studies have also considered spatially-varying temperature estimates developed by various modeling approaches. For example, a study in Brisbane, Australia showed that spatially-resolved exposure estimates obtained from geostatistical models gave a slightly better model fit and stronger exposure-response functions compared to observations at a central monitor or averages of multiple monitors [12]. Similar findings are reported by Lee et al. [13], in a study in Southeastern US with 1 km exposure estimated using a statistical model with satellite parameters. Finally, in a study of 113 US counties, Weinberger et al. [14], found that associations were largely similar with the use of weather station observations or population-weighted temperature estimates based on a 4 km temperature product. To our knowledge, our study is the first to compare impacts of different temperature exposure metrics using spatially-resolved temperature data products in analyzing ED visits instead of mortality.

Another unique aspect of our study is the evaluation of daily minimum temperature. Prolonged exposure to heat can overwork the body’s natural compensatory mechanisms. Studies have shown that hot nights measured by daily minimum temperature preceded or followed by hot days are associated with adverse health outcomes and excess mortality [34,35,36]. This has been attributed to the body not having time to recover from heat exposure. In order to maintain thermal homeostasis during heat events, cardiac output increases and the body redirects blood flow from vital organs to the skin, thus cooling the body [37]. Sweating also increases to cool the body, but excessive sweating can lead to reduced blood volume and dehydration. These mechanisms are often attenuated in the elderly or populations with comorbidities and could be a contributing factor to the increases seen in mortality risk when using a minimum temperature threshold [36, 38]. Murage et al. [36], found mortality risks were increased on hot nights followed by a hot day more than on cool nights followed by hot days. Overall, the above studies examined minimum temperature and mortality risks but could help explain the strong associations seen with minimum temperatures and morbidity for this analysis. We also note that maximum daily temperature typically occurs after half of the day has elapsed. Hence the weaker observed associations may reflect the temporal misalignment between outcome and short-term exposure. In particular, Davis et al. [39], examined associations between hourly temperature and mortality, and found that yesterday’s afternoon temperature had strong associations in Boston, Philadelphia and Seattle. Hence, future work may consider exploring the use of Daymet products in estimating exposure lag structures for different outcomes and temperature characteristics.

We found that differences in estimated RR’s between airport observations and Daymet-based exposures to be larger for minimum temperature than maximum temperature. One possible explanation is that minimum temperature corresponds to night-time exposure when individuals are most likely to be at home. Hence, using population-weighted average may better reflect average exposures. In contrast, population-weighted maximum temperature exposure may still be subject to considerable error due to variation in between-individual time-activity patterns. Another possible explanation is that urbanization has stronger impacts on night time temperature compared to daytime hours, contributing to larger spatial variation in minimum temperature that may not be captured by airport monitors. We also note that distributions of Daymet-based minimum temperatures exposures were more different than airport observations (Table S2). For example, in Los Angeles, the differences between the 95th and 50th percentile were 4.25 °C for Daymet ZCTA population-weighted average versus 3.33 °C for airport observations. Hence, the larger absolute change in exposure may also contribute to higher risk estimates.

Although Daymet is a well-established meteorology product, one limitation is its incomplete representation of urban heat islands (UHIs), which comes from several factors. First, Daymet uses stations compiled within the Global Historical Climatology Network-Daily (GHCN-D) dataset, from the US National Centers for Environmental Information. For the US, the GHCN-D dataset compiles US automated surface observing stations (ASOS) which include airports, cooperative network stations (e.g. US National Weather Service cooperative network stations), and some regional mesonets (e.g. Remote Automatic Weather Station [RAWS],), providing roughly 10–15,000 temperature observations a day across the contiguous US through the study period [40], which may have insufficient density in to fully observe UHIs. Additionally, Daymet uses a truncated Gaussian filter to perform the spatial interpolation of the station observations to the grid (on average 20 stations per grid point), effectively smoothing small scale observed variations in regions of high observation density. Finally, Daymet uses elevation as a predictor within their model (as elevation strongly influences temperature), but does not use any other geophysical attributes such as urban fraction [26]. These three factors limit the representation of UHIs in Daymet to resemble smooth bullseyes with possibly reduced magnitudes. We note that the generalized concepts of sparse stations, interpolation functions, and limited geophysical predictors are common to essentially all currently available in situ observation based gridded meteorology products.

Possible extensions to alleviate these concerns include use of supplemental station networks, use of other geophysical attributes as predictors in statistical models, or incorporation of urban models within the general interpolation framework to explicitly represent UHIs. For example, there are many thousands more stations available in the past 10–15 years through local mesonets that may improve representation of UHIs. There are also high-resolution land use/land cover maps that allow for inclusion of other geophysical predictors in a statistical model, or use of physically based urban models [41], that explicitly model the full energy balance of urban grid cells to better represent the spatial variability of UHIs [42]. These extensions are the subject of current research by this team [43].

Our study also has several limitations common to those that utilize administrative health databases. First, our exposure estimates do not account for individual activity patterns and indoor temperature. However, by using only temporal contrast in exposure, long-term trends that impact personal exposures should not confound estimated relative risks. Moreover, quantifying risks associated with ambient temperature is more relevant for designing local heat warning systems. Second, we did not control for ambient air pollution, such as ozone and fine particulate matter, which have been associated with various cardiorespiratory outcomes. Higher temperature often increases pollutant concentrations due to increased emission (e.g., from electricity generation) and favorable conditions for pollutant formation and transport. Hence, controlling for ambient pollution may reduce the total association of temperature on ED visits [44]. However, some multi-city mortality studies have also demonstrated that adjusting for ambient air pollution often lead to similar health effect estimates for temperature [20, 45, 46]. Finally, our empirical analyses do not allow for a rigorous quantification of attenuation in estimated relative risks due to spatial variation in heat exposures. A simulation study will help characterize how different degrees of spatial variation and the location of the single observation monitor can impact health effect estimation. The 1 km spatial resolution of Daymet can also be used to investigate whether this is an optimal spatial aggregation most practical for short-term exposure and health studies.


In summary, this study found positive health associations between high temperature and emergency department visits in five US cities located in different climate regions. The comparison of different exposure metrics suggests that epidemiological studies based on single monitoring stations may underestimate the effect of temperature on morbidity when monitoring station is less representative of the exposure of the at-risk population. Our results further demonstrate the potential advantages of using spatially-resolved, population-weighted exposure estimates in estimating health effects of short-term heat exposure, and can benefit from future studies that examine additional locations, health outcomes, and exposure lag structures.

Availability of data and materials

Temperature exposure data are publicly available from online databases. Health data can be requested from state-specific data custodians by establishing appropriate data use agreements.


  1. Zhang K, Oswald EM, Brown DG, Brines SJ, Gronlund CJ, White-Newsome JL, et al. Geostatistical exploration of spatial variation of summertime temperatures in the Detroit metropolitan region. Environ Res. 2011;111(8):1046–53. Epub 2011 Sept 14.

    Article  CAS  Google Scholar 

  2. Johnson S, Ross Z, Kheirbek I, Ito K. Characterization of intra-urban spatial variation in observed summer ambient temperature from the New York City Community air survey. Urban Clim. 2020;31:100583.

    Article  Google Scholar 

  3. Li D, Bou-Zeid E. Synergistic interactions between urban Heat Islands and heat waves: the impact in cities is larger than the sum of its parts. J Appl Meteorol Climatol. 2013;52(9):2051–6.

    Article  Google Scholar 

  4. Diem JE, Stauber CE, Rothenberg R. Heat in the southeastern United States: characteristics, trends, and potential health impact. PLoS One. 2017;12(5):e0177937.

    Article  CAS  Google Scholar 

  5. Goldman GT, Mulholland JA, Russell AG, Strickland MJ, Klein M, Waller LA, et al. Impact of exposure measurement error in air pollution epidemiology: effect of error type in time-series studies. Environ Health. 2011;10(1):61.

    Article  Google Scholar 

  6. Dionisio KL, Baxter LK, Chang HH. An empirical assessment of exposure measurement error and effect attenuation in bipollutant epidemiologic models. Environ Health Perspect. 2014;122(11):1216–24. Epub 2014 July 8.

    Article  CAS  Google Scholar 

  7. Sarnat SE, Sarnat JA, Mulholland J, Isakov V, Özkaynak H, Chang HH, et al. Application of alternative spatiotemporal metrics of ambient air pollution exposure in a time-series epidemiological study in Atlanta. J Expo Sci Environ Epidemiol. 2013;23(6):593–605. Epub 2013 Aug 21.

    Article  Google Scholar 

  8. Xiao Q, Chen H, Strickland MJ, Kan H, Chang HH, Klein M, et al. Associations between birth outcomes and maternal PM2.5 exposure in Shanghai: a comparison of three exposure assessment approaches. Environ Int. 2018;117:226–36. Epub 2018 May 12.

    Article  CAS  Google Scholar 

  9. Mohegh A, Goldberg D, Achakulwisut P, Anenberg SC. Sensitivity of estimated NO2-attributable pediatric asthma incidence to grid resolution and urbanicity. Environ Res Lett. 2020;16(1):014019.

    Article  CAS  Google Scholar 

  10. Field CB, Barros VR, Dokken DJ, et al. Climate change 2014 impacts, adaptation and vulnerability: part a: global and Sectoral aspects: working group II contribution to the fifth assessment report of the intergovernmental panel on. Climate Change. 2014.

  11. Hsiang S, Kopp R, Jina A, Rising J, Delgado M, Mohan S, et al. Estimating economic damage from climate change in the United States. Science. 2017;356(6345):1362–9.

    Article  CAS  Google Scholar 

  12. Guo Y, Barnett AG, Tong S. Spatiotemporal model or time series model for assessing city-wide temperature effects on mortality? Environ Res. 2013;120:55–62. Epub 2012 Sept 29.

    Article  CAS  Google Scholar 

  13. Lee M, Shi L, Zanobetti A, Schwartz JD. Study on the association between ambient temperature and mortality using spatially resolved exposure data. Environ Res. 2016;151:610–7. Epub 2016 Sept 7.

    Article  CAS  Google Scholar 

  14. Weinberger KR, Spangler KR, Zanobetti A, Schwartz JD, Wellenius GA. Comparison of temperature-mortality associations estimated with different exposure metrics. Environ Epidemiol. 2019;3(5):e072.

    Article  Google Scholar 

  15. Laaidi K, Zeghnoun A, Dousset B, Bretin P, Vandentorren S, Giraudet E, et al. The impact of heat islands on mortality in Paris during the august 2003 heat wave. Environ Health Perspect. 2012;120(2):254–9. Epub 2011 Sept 1.

    Article  Google Scholar 

  16. Harlan SL, Declet-Barreto JH, Stefanov WL, Petitti DB. Neighborhood effects on heat deaths: social and environmental predictors of vulnerability in Maricopa County, Arizona. Environ Health Perspect. 2013;121(2):197–204. Epub 2012 Nov 16.

    Article  Google Scholar 

  17. Klein Rosenthal J, Kinney PL, Metzger KB. Intra-urban vulnerability to heat-related mortality in New York City, 1997-2006. Health Place. 2014;30:45–60. Epub 2014 Sept 6.

    Article  Google Scholar 

  18. Jenerette GD, Harlan SL, Buyantuev A, Stefanov WL, Declet-Barreto J, Ruddell BL, et al. Micro-scale urban surface temperatures are related to land-cover features and residential heat related health impacts in Phoenix, AZ USA. Landscape Ecol. 2016;31(4):745–60.

    Article  Google Scholar 

  19. Anderson GB, Bell ML. Heat waves in the United States: mortality risk during heat waves and effect modification by heat wave characteristics in 43 U.S. communities. Environ Health Perspect. 2011;119(2):210–8. Epub 2010 Oct 7.

    Article  Google Scholar 

  20. Gasparrini A, Guo Y, Hashizume M, Lavigne E, Zanobetti A, Schwartz J, et al. Mortality risk attributable to high and low ambient temperature: a multicountry observational study. Lancet. 2015;386(9991):369–75. Epub 2015 May 20.

    Article  Google Scholar 

  21. Luo Q, Li S, Guo Y, Han X, Jaakkola JJK. A systematic review and meta-analysis of the association between daily mean temperature and mortality in China. Environ Res. 2019;173:281–99. Epub 2019 Mar 22.

    Article  CAS  Google Scholar 

  22. Ye X, Wolff R, Yu W, Vaneckova P, Pan X, Tong S. Ambient temperature and morbidity: a review of epidemiological evidence. Environ Health Perspect. 2012;120(1):19–28. Epub 2011 Aug 8.

    Article  Google Scholar 

  23. Bobb JF, Obermeyer Z, Wang Y, Dominici F. Cause-specific risk of hospital admission related to extreme heat in older adults. JAMA. 2014;312(24):2659–67.

    Article  CAS  Google Scholar 

  24. Phung D, Thai PK, Guo Y, Morawska L, Rutherford S, Chu C. Ambient temperature and risk of cardiovascular hospitalization: an updated systematic review and meta-analysis. Sci Total Environ. 2016;550:1084–102. Epub 2016 Feb 9.

    Article  CAS  Google Scholar 

  25. Vaidyanathan A, Saha S, Vicedo-Cabrera AM, Gasparrini A, Abdurehman N, Jordan R, et al. Assessment of extreme heat and hospitalizations to inform early warning systems. Proc Natl Acad Sci U S A. 2019;116(12):5420–7. Epub 2019 Mar 4.

    Article  CAS  Google Scholar 

  26. Thornton PE, Running SW, White MA. Generating surfaces of daily meteorological variables over large regions of complex terrain. J Hydrol. 1997;190(3):214–51.

    Article  Google Scholar 

  27. Spangler KR, Weinberger KR, Wellenius GA. Suitability of gridded climate datasets for use in environmental epidemiology. J Expo Sci Environ Epidemiol. 2019;29(6):777–89. Epub 2018 Dec 11.

    Article  CAS  Google Scholar 

  28. Ivy D, Mulholland JA, Russell AG. Development of ambient air quality population-weighted metrics for use in time-series health studies. J Air Waste Manag Assoc. 2008;58(5):711–20.

    Article  CAS  Google Scholar 

  29. Winquist A, Grundstein A, Chang HH, Hess J, Sarnat SE. Warm season temperatures and emergency department visits in Atlanta, Georgia. Environ Res. 2016;147:314–23. Epub 2016 Feb 27.

    Article  CAS  Google Scholar 

  30. Beck HE, Zimmermann NE, McVicar TR, Vergopolan N, Berg A, Wood EF. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci Data 2018;5:180214. doi: Erratum in: Sci Data. 2020 Aug 17;7(1):274.

  31. Lippmann SJ, Fuhrmann CM, Waller AE, Richardson DB. Ambient temperature and emergency department visits for heat-related illness in North Carolina, 2007-2008. Environ Res. 2013;124:35–42. Epub 2013 Apr 30.

    Article  CAS  Google Scholar 

  32. Tobías A, Armstrong B, Gasparrini A. Brief report: investigating uncertainty in the minimum mortality temperature: methods and application to 52 Spanish cities. Epidemiology. 2017;28(1):72–6.

    Article  Google Scholar 

  33. Patzert W, LaDochy S, Ramirez P, Willis J. Los Angeles Weather Station’s relocation impacts climate and weather records. The California Geographer, vol. 55; 2016.

    Google Scholar 

  34. Napoli C, Pappenberger F, Cloke H. Verification of heat stress thresholds for a health-based heat-wave definition. J Appl Meteorol Climatol. 2019;58(6):1177–94.

    Article  Google Scholar 

  35. Loughnan M, Nicholls N, Tapper N. Mortality-temperature thresholds for ten major population centres in rural Victoria, Australia. Health Place. 2010;16(6):1287–90. Epub 2010 Aug 14.

    Article  Google Scholar 

  36. Murage P, Hajat S, Kovats RS. Effect of night-time temperatures on cause and age-specific mortality in London. Environ Epidemiol. 2017;1(2):e005. Epub 2017 Dec 13.

    Article  Google Scholar 

  37. Osilla EV, Marsidi JL, Sharma S. Physiology, Temperature Regulation. In: StatPearls. Treasure Island (FL): StatPearls Publishing; 2021.

    Google Scholar 

  38. Kenney WL, Craighead DH, Alexander LM. Heat waves, aging, and human cardiovascular health. Med Sci Sports Exerc. 2014;46(10):1891–9.

    Article  Google Scholar 

  39. Davis RE, Hondula DM, Patel AP. Temperature observation time and type influence estimates of heat-related mortality in seven U.S. cities. Environ Health Perspect. 2016;124(6):795–804.

    Article  CAS  Google Scholar 

  40. Menne MJ, Durre I, Vose RS, Gleason BE, Houston TG. An overview of the global historical climatology network-daily database. J Atmos Ocean Technol. 2012;29(7):897–910.

    Article  Google Scholar 

  41. Kusaka H, Kondo H, Kikegawa Y, Kimura F. A simple single-layer urban canopy model for atmospheric models: comparison with multi-layer and slab models. Bound-Layer Meteorol. 2001;101(3):329–58.

    Article  Google Scholar 

  42. Monaghan AJ, Hu L, Brunsell NA, Barlage M, Wilhelmi OV. Evaluating the impact of urban morphology configurations on the accuracy of urban canopy model temperature simulations with MODIS. J Geophys Res. 2014;119(11):6376–92.

    Article  Google Scholar 

  43. Newman AJ, Monaghan AM, Holmes H, Chang H, Darrow L, Warren J, et al. Development of a novel high-resolution long-term gridded meteorology dataset including urban heat islands over the contiguous United States for health studies. San Francisco: AGU Fall Meeting; 2019.

    Google Scholar 

  44. Buckley JP, Samet JM, Richardson DB. Commentary: does air pollution confound studies of temperature? Epidemiology. 2014;25(2):242–5.

    Article  Google Scholar 

  45. Anderson BG, Bell ML. Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009;20(2):205–13.

    Article  Google Scholar 

  46. Chen R, Yin P, Wang L, Liu C, Niu Y, Wang W, et al. Association between ambient temperature and mortality risk and burden: time series study in 272 main Chinese cities. BMJ. 2018;363:k4306.

    Article  Google Scholar 

Download references


Not applicable.


This research was supported by the National Institute of Environmental Health Sciences (NIEHS) of the National Institutes of Health (NIH) under Award R01 ES027892, R01 ES028346, and R21 ES023763. The content of this manuscript is solely the responsibility of the authors and does not necessarily represent the official views of the NIH and does not endorse the purchase of any commercial products or services mentioned in this manuscript.

Author information

Authors and Affiliations



NT contributed in methodology, analysis and writing. STE contributed in funding acquisition, data acquisition and writing. AJN contributed to conceptualization and writing. RRD contributed to data analysis. NS, SM, JLW, and MJS contributed to manuscript review. LAD contributed to funding acquisition and manuscript review. HHC contributed in conceptualization, methodology, analysis, writing, supervision, and funding acquisition. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Howard H. Chang.

Ethics declarations

Ethics approval and consent to participate

The study used de-identified hospital billing records and participant consent is not required. The study has been approved by the Institutional Review Board at Emory University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Descriptive statistics of exposures by city, temperature metric and exposure assessment methods. Table S2. Quantile values of exposures by city and exposure metrics. Table S3. Difference in AIC and ratio of overdispersion by city, outcome and temperature exposure. Figure S1. Maps of average May-October daily Daymet 1km temperature estimates by city. Figure S2. Associations with minimum temperature between the 75th and 25th percentile of exposures. Figure S3. Associations with maximum temperature between the 75th and 25th percentile of exposures. Figure S4. Associations with average temperature between the 95th and 50h percentile of exposures. Figure S5. Associations with average temperature between the 75th and 25h percentile of exposures. Figure S6. Associations with minimum temperature with Daymet grid cell link to the monitor. Figure S7. Associations with maximum temperature with Daymet grid cell link to the monitor. Figure S8. Sensitivity analyses for different natural cubic spline degrees of freedom for minimum temperature. Figure S9. Sensitivity analyses for different natural cubic spline degrees of freedom for maximum temperature. Figure S10. Associations for 3-day moving average of minimum temperature between the 95th and 50th percentile of exposures. Figure S11. Associations for 3-day moving average of maximum temperature between the 95th and 50th percentile of exposures. Primary.csv Model results for primary analysis for same-day exposure (by city, outcome, temperature metric, exposure assessment methods, AIC, overdispersion, and exposure contrast referenced at the 50th exposure percentile). Primary.csv Model results for primary analysis for same-day exposure (by city, outcome, temperature metric, exposure assessment methods, AIC, overdispersion, and exposure contrast referenced at the 50th exposure percentile). Primary_Ref25.csv Model results for primary analysis for same-day exposure (by city, outcome, temperature metric, exposure assessment methods, AIC, overdispersion, and exposure contrast referenced at the 25th exposure percentile). PrimaryMV.csv Model results for primary analysis for 3-day moving average exposure (by city, outcome, temperature metric, exposure assessment methods, AIC, overdispersion, and exposure contrast referenced at the 50th exposure percentile).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Thomas, N., Ebelt, S.T., Newman, A.J. et al. Time-series analysis of daily ambient temperature and emergency department visits in five US cities with a comparison of exposure metrics derived from 1-km meteorology products. Environ Health 20, 55 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: