Estimated Effect of Temperature on Years of Life Lost: A Retrospective Time-Series Study of Low-, Middle-, and High-Income Regions.

BACKGROUND: Numerous studies have reported a strong association between temperature and mortality. Additional insights can be gained from investigating the effects of temperature on years of life lost (YLL), considering the life expectancy at the time of death.
OBJECTIVES: The goal of this work was to assess the association between temperature and YLL at seven low-, middle-, and high-income sites.
METHODS: We obtained meteorological and population data for at least nine years from four Health and Demographic Surveillance Sites in Kenya (western Kenya, Nairobi), Burkina Faso (Nouna), and India (Vadu), as well as data from cities in the United States (Philadelphia, Phoenix) and Sweden (Stockholm). A distributed lag nonlinear model was used to estimate the association of daily maximum temperature and daily YLL, lagged 0-14 d. The reference value was set for each site at the temperature with the lowest YLL.
RESULTS: Generally, YLL increased with higher temperature, starting day 0. In Nouna, the hottest location, with a minimum YLL temperature at the first percentile, YLL increased consistently with higher temperatures. In Vadu, YLL increased in association with heat, whereas in Nairobi, YLL increased in association with both low and high temperatures. Associations with cold and heat were evident for Phoenix (stronger for heat), Stockholm, and Philadelphia (both stronger for cold). Patterns of associations with mortality were generally similar to those with YLL.
CONCLUSIONS: Both high and low temperatures are associated with YLL in high-, middle-, and low-income countries. Policy guidance and health adaptation measures might be improved with more comprehensive indicators of the health burden of high and low temperatures such as YLL. https://doi.org/10.1289/EHP1745.

Introduction

Ambient temperature is an important determinant of health. Many studies have reported that mortality and morbidity increase with very high and/or very low ambient temperatures (e.g., Analitis et al. 2008; Baccini et al. 2008; Basagana et al. 2015; Basu and Samet 2002; Gasparrini et al. 2015; Kenny et al. 2010). Although evidence on the temperature–health association in wealthier countries is well established, research is still emerging on this topic from low- and middle-income countries (LMICs) (Azhar et al. 2014; Burkart et al. 2014; Diboulo et al. 2012; McMichael et al. 2008; Mrema et al. 2012). Because weather-related health outcomes are dependent on local contexts, the large volume of research generated on the temperature–health association in wealthier countries may not necessarily reflect the burden from temperature exposure in LMICs. In contrast to high-income countries, the burden of disease in LMICs is often characterized by higher youth mortality rates and high prevalence of communicable diseases. Health impacts from heat and cold in LMICs may therefore occur through causal pathways other than those identified in high-income countries. For example, lower life expectancy in low-income countries might limit the pool of vulnerable elderly individuals, who suffer disproportionately from heat–health risks (Kenny et al. 2010). In LMICs, a large proportion of the population is working in physically demanding jobs on the streets, in agriculture, or on construction sites. Greater occupational exposure to environmental risk factors in these countries might lead to increased adverse health outcomes in comparison with the outcomes in wealthier countries.

A second observation of the literature regarding temperature–health associations is that most research uses the total number of daily deaths as the primary outcome (Bunker et al. 2016; Hajat and Kosatky 2010). Daily deaths are attractive to use for a number of reasons, including the accessibility of records (in some, but not all, jurisdictions), as well as high interpretability of study findings. Death counts, however, do not necessarily provide an ideal representation of the total mortality burden attributable to high and low temperatures. A major shortcoming of using daily death counts as the outcome is that this measure does not indicate the extent to which lives are shortened by exposure to heat or cold. At short time scales, such as days to weeks, many studies suggest that a fraction of the attributable mortality to heat and cold is from individuals who would have died anyway within days or weeks, a phenomenon known as mortality displacement (e.g., Hajat et al. 2005). For these decedents, it may be more appropriate instead to largely attribute their death to pre-existing ill-health condition.

Beyond short-term displacement, an additional drawback of using total death counts is that the loss of a young person’s life is equally summed in the total burden as the loss of an old or elderly person’s life. Although all lives have value, the loss of a young person’s life leads to a greater potential loss of societal contributions (Rocklöv et al. 2009).

One approach to avoid these drawbacks is to consider years of life lost (YLL) as the outcome measure instead of death counts. YLL is an indicator of premature mortality used as a global burden of disease death metric (Lopez et al. 2006). The YLL approach accounts for the age at which death occurred by giving greater weight to deaths at a younger age. YLL is a compound measure combining the number of daily deaths with age at death; it is considered a more informative and differentiated measurement for assessing premature mortality than total or age-specific mortality rates alone. A comparison of health burdens by different exposures can facilitate risk ranking for preventive interventions, public health planning, and resource allocation (Aragón et al. 2008; Gardner and Sanborn 1990; Huang et al. 2012a, 2012b; Yang et al. 2015). Although common in other risk-assessment fields, to our knowledge YLL has rarely been applied in studies investigating the temperature–health association (Baccini et al. 2013; Egondi et al. 2015; Huang et al. 2012a, 2012b; Xu et al. 2014; Yang et al. 2015).

To expand our understanding of the impact of temperature on health and to address drawbacks of existing literature, we estimated immediate and lagged effects of high and low daily temperatures on YLL across seven high-, middle- and low-income study sites. These sites were selected to cover a range of socioeconomic and climatic settings in Africa, Asia, Europe, and North America, based on available population and mortality data. Research questions motivating the analyses were: a) What is the association of heat and cold with YLL? b) How do temperature–YLL associations vary among climatically and socioeconomically diverse locations?

Methods

Study Sites

The seven study sites are a mix of urban and rural locations, including Nouna (Burkina Faso), Nairobi and Kisumu (Kenya), Vadu (India), Phoenix and Philadelphia (USA) and Stockholm (Sweden) (Figure 1). Study locations in Africa and India are Health and Demographic Surveillance Sites (HDSS), which are part of the International Network for the Demographic Evaluation of Populations and their Health (INDEPTH) (Sankoh and Byass 2012). The sites in Sweden and the United States included the total municipal population. The locations cover a range of different climates, population characteristics, degree of economic development, and urbanization. The locations in the United States and Sweden were intended to be illustrative of cities in different climate zones.

Figure 1.

Map depicting the location of study sites. The map was generated using ArcMap (version 10.5; ESRI, Inc.) with topographical basemap content from ESRI, the United States Geological Survey, and the United States National Oceanic and Atmospheric Administration.

Study periods varied from site to site, but all covered at least nine consecutive years between 1990 and 2012 (Table 1). Because daily mortality data are difficult to obtain in low-income countries, and even in high-income countries, we chose study sites for which such data were available. Sites were selected based on existing collaboration of site representatives with the Umeå Centre of Global Health Research (Rocklöv et al. 2012) and on the availability of health and weather data for the relevant study period.

Table 1

Descriptive statistics of exposure and outcome measures across the seven study sites.

Population characteristics	Study site (study period, location)
	Kisumu (2003–2011, Villages Asembo, Gem, Karemo, rural Kenya)	Nairobi (2003–2012, Slums Korogocho, Viwandani, urban Kenya)	Nouna (2000–2010, Kossi province, rural Burkina Faso)	Vadu (2003–2012, 22 villages east of Pune, rural India)	Phoenix (1997–2007, Urban Arizona, USA)	Philadelphia (1997–2007, Urban Pennsylvania, USA)	Stockholm (1990–2002, Stockholm county, urban Sweden)
Life expectancy at birth (y)a	Men 54.8 Women 64.0	Men 72.0 Women 65.9	Men 58.0 Women 61.5	Men 68.1 Women 71.5	Men 75.3 Women 81.2	Men 74.1 Women 79.9	Men 77.3 Women 82.0
Population sizeb	240,633	63,639	ca. 78,000	110,085	3,251,876	1,517,550	1,823,210
Daily max. temperature (°C)
Minimum	15.0	15.0	22.8	21.1	7.2	−8.9	−13.4
5th percentile	26.0	21.1	30.0	26.6	16.0	1.1	−2.0
Median	29.0	26.0	36.4	31.3	30.6	18.3	9.6
Mean	29.2	25.8	36.1	32.0	29.9	17.9	10.5
95th percentile	33.0	29.6	41.7	38.8	42.8	32.2	24.6
Maximum	37.0	38.2	43.9	42.4	47.2	37.8	33.5
Standard deviation	2.2	2.6	3.6	3.8	9.0	10.1	8.6
Daily deaths
Minimum	0	0	0	0	24	16	19
5th percentile	2.0	0.0	0.0	0.0	46.0	31.0	29.0
Median	7.0	1.0	1.0	1.0	62.0	44.0	40.0
Mean	8.7	1.3	1.7	0.9	61.9	44.5	40.0
95th percentile	15.7	4.0	4.0	3.0	80.0	59.0	52.0
Maximum	124	57	32	24	117	79	74
Standard deviation	9.6	1.9	2.8	1.2	10.5	8.4	7.2
Daily YLL
Minimum	0.0	0.1	0.0	0.0	387.9	173.7	119.4
5th percentile	57.7	0.1	0.0	0.0	605.0	452.2	284.2
Median	212.0	40.3	41.5	9.3	893.6	695.1	427.2
Mean	272.1	56.6	66.2	24.4	904.0	703.9	431.7
95th percentile	1,357.7	179.2	187.3	91.8	1,227.7	991.1	598.2
Maximum	4,326.3	2,324.4	1,607.3	568.3	1,703.7	1,329.5	818.1
Standard deviation	314.6	82.0	125.5	38.0	189.9	164.7	94.7
Minimum YLL temperature
Temperature (°C)	30.0	26.6	27.8	29.8	29.4	22.2	27.0
Percentile	62	60	1	31	47	60	98

Note: YLL, years of life lost at the age at death, based on conditional life expectancies. YLL were calculated separately for each as the difference between years of age at the time of death and estimated life expectancy given the decedent’s age and gender. Daily YLL are the total YLL of all persons who died on the same day. Minimum YLL temperature is the calculated daily maximum temperature in °C at which average YLL were lowest. HDSS, Health and Demographic Surveillance Site.

Life expectancy at birth: Site-specific data for HDSS sites in last study year. State-level or national data from official registries: Phoenix data of Arizona (1999–2001); Philadelphia data of Pennsylvania (1999–2001); Stockholm county data of Sweden (2002). Not used for YLL analyses.

Population size: HDSS sites in last study year; Phoenix, Philadelphia, and Stockholm in 2000.

Population Data

Daily mortality data were collected from HDSS sites in Kisumu, Nairobi, Nouna, and Vadu. The HDSS data contain individual statistics including gender, date of birth, date of death, and dates of in- and out-migration to and from the HDSS site. Details of the data collection procedures at these INDEPTH sites are described in Sankoh and Byass (2012). Daily mortality data for Stockholm, Phoenix, and Philadelphia were obtained from state or national population registers and vital statistics offices (Hondula et al. 2015b; Rocklöv et al. 2009). For all seven study sites, information about date and age at death of all deceased individuals was available. All causes of death were included in the analyses.

As an indicator of general population health and infant mortality, life expectancy at birth was calculated for each HDSS site, and respective data for U.S. and Swedish sites were obtained from registers for State of Arizona (Phoenix site), the State of Pennsylvania (Philadelphia site) and Sweden (Stockholm site).

Calculation of Daily Years of Life Lost

A study protocol was developed and shared among coauthors to obtain standardized daily YLL for each study site. No city-specific data of conditional life expectancies were available for the three sites in the United States and Sweden. Therefore, national abridged conditional life expectancies for five-year age groups, based on age-specific mortality rates, were used for these sites. We used age-specific mortality rates for each HDSS site to generate site-specific five-year abridged life tables, stratified by gender, for estimating life expectancy. A separate age band was created for infants under 1 y of age to account for high infant mortality at the HDSS sites. For each death occurring at the different sites, we calculated individual YLL. These YLL were estimated as the site-specific remaining life expectancy on a population level at the individual age of death. For example, if the remaining life expectancy of women age 60 to 64 y was 17.9 y, the individual YLL for a woman who died at that age would be 17.9 y. All YLL values were greater than or equal to zero. Individual YLL were calculated separately for men and women due to differences in life expectancy between the sexes. Daily YLL were the sum of individual YLL occurring on each study day.

Weather Data

Daily maximum temperature data were obtained from the nearest weather station for each study site, including Kisumu: Kisumu airport; Nairobi: Moi Airbase Eastleigh weather station (if missing, data from Jomo Kenyatta International Airport); Nouna: Dedougou weather station; Vadu: Pune Airport; Phoenix: Sky Harbor International Airport; Philadelphia: Philadelphia International Airport; Stockholm: Bromma Airport. Maximum temperature was used as the exposure variable as it was the most consistently available commonly examined temperature metric available across all sites.

Statistical Analyses

The effect of daily maximum temperature on YLL was estimated using a distributed lag nonlinear modelling (DLNM) framework (Gasparrini 2014b). The DLNM is a flexible strategy that incorporates both the lag as well as the exposure dimension simultaneously in the model. This is performed through the specification of a cross-basis function defining the maximum lag of the exposure variable with a basis function for both the exposure variable and the lag dimension. The basis function can be defined to model the effect of the exposure variable as a nonlinear or linear function.

In our analyses, we constructed a cross-basis function of daily maximum temperature with a 2-degree b-spline basis with knot placements at the second and 98th percentile. We considered lags 0 to 14 d with knot placements for the lags at three equally spaced positions. The lag dimension was modeled with a natural cubic spline basis.

To control for trend and seasonality of YLL, we included a natural cubic spline function of time with 6 degrees of freedom per year. Models also included indicator variables for day of week. In addition, models for the HDSS sites included an indicator for “heaping days,” which are dates that are customarily assigned to deaths for which the actual date of death is unknown, specifically, the 15th of each month for Kisumu and Nairobi, the ninth of each month for Nouna, and the first of January for Vadu. The sites in Sweden and the United States did not require control for heaping.

The complete model included a cross-basis function of maximum temperature, spline function of time trend to control for seasonality, the day of week and heaping day in the HDSS sites. We modeled daily YLL using a censored Gaussian distribution with Tobit (McDonald and Moffitt 1980; Tobin 1958) to keep days with (which were common at the four HDSS sites; see Figures S1 and S2) in the model.

The degree of freedom for the trend function and the knot placement for maximum temperature were chosen after testing different combinations of values. We evaluated alternative options for placement of the knots for the b-spline of maximum daily temperature (second and 98th; fifth, 25th, and 95th; 10th, 70th, 90th; second, 70th, 98th; second, 25th, 70th, 90th; and fifth, 50th, 70th, and 95th percentiles, respectively) and for the degrees of freedom (df) for the seasonal trend function (4, 8, 6, or 10 df per year). To facilitate the interpretation of effect estimates across sites, we chose the parameters for the final models (b-spline knots at the second and 98th percentiles and 6 df for the seasonal trend) that produced the lowest total AIC summed over the seven sites (see Table S1). Although this metric does not necessarily yield the optimal model at each specific site, holding model parameters consistent across the sites facilitates interpretation of effect estimates.

We estimated differences in YLL with lower or higher daily maximum temperatures relative to the YLL at the daily maximum temperature with the smallest YLL value. To determine the reference temperature for each site, we predicted the YLL for each percentile of the daily maximum temperature distribution, and identified the percentile with the lowest YLL value.

Mortality comparison.

Associations between daily maximum temperature and mortality were estimated using a quasi-Poisson model with the same cross-basis function for maximum temperature as the Gaussian model for YLL. As for the YLL models, we included a cubic spline trend function for seasonality, indicator variables for day of week, and an indicator for heaping days at the HDSS sites. Associations with higher and lower temperatures were estimated relative to the reference temperature used for the YLL analyses at each site.

Sensitivity analysis

We performed several different sensitivity analyses by varying the degrees of freedom for the trend function and the percentile for the knot placement for the different sites and compared how the YLL maximum temperature relationship changed. In addition, we reran the analysis using a natural cubic spline basis for maximum temperature instead of a b-spline, applied similar knot percentile positions, and compared the exposure response relationship.

We used STATA software (versions 11 to 13; StataCorp) for data management and descriptive statistics. The R statistical software (R Development Core Team 2015) was used for all analyses, the DLNM package (Gasparrini 2011, 2014a) for DLNM modeling and the VGAM package for the Tobit regression models.

Results

Summary statistics of the seven study sites are presented in Table 1. Among the seven study sites, Nouna had the highest average daily maximum temperature (36.1°C), and Stockholm was the coldest location (average daily maximum temperature of 10.5°C). At the HDSS sites, which are located in equatorial regions, variation of temperature was considerably smaller than in Sweden and the United States. The standard deviation of temperature was lowest in Kisumu (2.2°C) and highest in Philadelphia (10.1°C). Table 1 shows also the site-specific temperature and temperature percentile at which the lowest YLL was observed (“minimum YLL temperature”).

The numbers of deaths per day varied among the sites, with HDSS sites frequently having days with zero deaths, consistent with their smaller population sizes (Table 1). Daily YLL as a consequence of daily deaths also varied among the study sites and within each site.

Seasonal Variations in Temperature and Daily YLL

Average daily YLL varied over the year in all locations, with colder months tending to have higher average daily YLL (see Figure 2 for monthly averages at each site, and Figure S1 for time-series data). However, seasonal patterns were less pronounced for Nouna and Vadu than for the other sites. Seasonal variations in daily death counts were similar to those for YLL at each site (data not shown).

Figure 2.

Average daily maximum temperature and YLL per month at seven study sites. The YLL range vary between the sites in the HDSS and in high income sites. YLL, years of life lost at the age at death, based on conditional life expectancies. YLL were calculated separately for each site as the difference between years of age at the time of death and estimated life expectancy given the decedent’s age and gender. Complete time-series for each site are shown in Figure S1.

Association between Daily Maximum Temperature and YLL

Low- and middle-income sites.

At the Kenyan sites (Kisumu and Nairobi), associations between YLL and daily maximum temperature were U-shaped (Figures 3A, B). The point of minimum YLL temperatures (and thus the lowest estimated YLL) was near the middle of the temperature range for each site (30°C and 27°C, the 62nd and 60th percentiles, respectively); estimated YLL for maximum daily temperatures were higher above and below these values. Estimates were imprecise, but differences in YLL relative to the minimum YLL were greater for daily maximum temperatures at the fifth percentile than temperatures at the 95th percentile (Table 2).

Figure 3.

Association between daily maximum temperature and years of life lost at low- and middle-income study sites. The overall estimates for each location (first column of panels A, B, C, and D, respectively) show YLL by temperature, cumulative over 14 days. The second and third columns show estimated differences in the YLL at the 5th or 95th percentile of maximum daily temperature, respectively, relative to the YLL at the reference temperature at each site, on the day of death () and up to 14 days prior to death. Tmax, daily maximum temperature; YLL, years of life lost. YLL estimates indicate the difference in years of life lost (the difference between life expectancy at the age of death and the age of death) relative to the YLL at the reference temperature for each site (the daily maximum temperature associated with the lowest average daily YLL). Gray shadings represent 95% confidence interval bands of estimates. Note: (A) Kisumu, (B) Nairobi, (C) Nouna, and (D) Vadu.

Table 2

Estimated cumulative difference in YLL (14-day lag) at different percentiles of maximum daily temperature relative to the YLL at the reference temperature.

Study site	Kisumu	Nairobi	Nouna	Vadu	Phoenix	Philadelphia	Stockholm
Minimum YLL temperature
Temperature (°C)	30.0	26.6	27.8	29.8	29.4	22.2	27.0
Percentile	62	60	1	31	47	60	98
5th percentile
Temp. (°C)	26.0	21.1	30.0	26.6	16.0	1.1	−2.0
YLL (CI)	19.47 (−26.73, 65.67)	34.36 (−1.03, 69.75)	17.67 (−0.87, 36.22)	3.38 (−6.87, 13.62)	46.46 (−5.07, 97.98)	48.25 (−10.19, 106.69)	46.2 (1.4, 91.01)
25th percentile
Temp. (°C)	28.0	24.1	33.3	29.2	22.2	10.0	3.6
YLL (CI)	4.82 (−12.58, 22.22)	7.17 (−4.17, 18.52)	43.33 (5.44, 81.21)	0.13 (−1.41, 1.67)	13.39 (−9.76, 36.55)	15.97 (−14.13, 46.06)	45.88 (6.01, 85.76)
50th percentile
Temp. (°C)	29.0	26.0	36.4	31.3	30.6	18.3	9.6
YLL (CI)	1.18 (−7.23, 9.6)	0.44 (−2.54, 3.42)	63.76 (17.89, 109.63)	0.7 (−2.63, 4.03)	0.38 (−3.39, 4.14)	1.55 (−10.03, 13.12)	40.98 (5.67, 76.29)
75th percentile
Temp. (°C)	30.4	27.6	38.9	35.1	37.8	27.2	17.5
YLL (CI)	0.21 (−3.59, 4.02)	1.07 (−4.92, 7.06)	77.66 (30.49, 124.83)	8.98 (−2.76, 20.73)	18.33 (−15.95, 52.61)	2.91 (−16.99, 22.8)	27.31 (2.6, 52.03)
95th percentile
Temp. (°C)	33.0	29.6	41.7	38.8	42.8	32.2	24.6
YLL (CI)	11.19 (−32.35, 54.73)	9.97 (−14.23, 34.16)	90.48 (41.92, 139.04)	26.03 (−0.36, 52.42)	46.6 (−19.31, 112.52)	11.31 (−35.56, 58.18)	8.36 (0.39, 16.33)

Note: The reference temperature is the maximum daily temperature with the lowest predicted YLL at each site (“Minimum YLL temperature,” first row). YLL, years of life lost at the age at death, based on conditional life expectancies. YLL were calculated separately for each site as the difference between years of age at the time of death and estimated life expectancy given the decedent’s age and gender. YLL, years of life lost; CI, 95% confidence interval; temp, daily maximum temperature.

(, 65.67)

(, 69.75)

(, 36.22)

(, 13.62)

(, 97.98)

(, 106.69)

(, 22.22)

(, 18.52)

(, 1.67)

(, 36.55)

(, 46.06)

(, 9.6)

(, 3.42)

(, 4.03)

(, 4.14)

(, 13.12)

(, 4.02)

(, 7.06)

(, 20.73)

(, 52.61)

(, 22.8)

(, 54.73)

(, 34.16)

(, 52.42)

(, 112.52)

(, 58.18)

The minimum YLL temperature for Nouna, the hottest site, was set at the first percentile of the daily maximum temperature distribution (28°C) (Table 1). Estimated YLL at temperatures above this value increased continuously (Figure 3C). Of all sites, Nouna showed the strongest association between YLL and temperature at the 95th percentile (42°C: 90 YLL; 95% CI: 42, 139) (Table 2). At Vadu, the second-hottest site, the minimum YLL temperature was at the 31st percentile of the daily maximum temperature distribution (30°C). In comparison with this reference value, YLL were notably higher at higher temperatures only (26 YLL; 95% CI: , 52 at the 95th percentile, 39°C) (Table 2, Figure 3D).

Although lag-specific estimates at the four HDSS sites were imprecise due to small population sizes, they hinted generally at immediate effects of high temperatures, and at the Kenyan sites for delayed responses (Figure 3A and 3B and Figure S3). YLL increased with heat in Nouna on the day of exposure and remained elevated up to day 5, but rose again after lag 10. At Kisumu, cold-related YLL increased at lag 6 to 10; in Nairobi, a small increase of YLL at low temperatures was observed with a two-week delay.

Patterns of associations between daily maximum temperatures and mortality (daily deaths) were generally similar to the patterns of associations with YLL at each HDSS site (Figure S4).

High-income sites.

At all three high-income sites, YLL increased at both high and low temperatures (Figure 4), relative to the temperature of minimum YLL. Philadelphia, the cooler of the two sites in the United States, had the widest range of temperatures ( to 38°C), 22°C (60th percentile) being the point of minimum YLL. At this site, the increase in YLL at the fifth percentile of temperature (1.1°C) and below was more pronounced than that at the 95th percentile and above (Figure 4B, Table 2). Estimated cold effects lasted for about one week (Figure S3). In Phoenix, the minimum YLL temperature was at the 47th percentile of the temperature range (29°C). Despite its hotter climate, estimated cold effects at the fifth percentile (16.0°C) were similar to heat effects at the 95th percentile (42.8°C) (Table 2). Estimated cold effects were most discernible at lag 12–14 d, whereas heat effects appeared to occur within the first four days (Figure S3).

Figure 4.

Association between daily maximum temperature and YLL at high-income study sites. The overall estimates for each location (first column of panels A, B, and C, respectively) show YLL by temperature, cumulative over 14 days. The second and third columns show estimated differences in the YLL at the 5th or 95th percentile of daily maximum temperature, respectively, relative to the YLL at the reference temperature at each site, on the day of death () and up to 14 days prior to death. Tmax, daily maximum temperature; YLL, years of life lost. YLL estimates indicate the difference in years of life lost (the difference between life expectancy at the age of death and the age of death) relative to the YLL at the reference temperature for each site (the daily maximum temperature associated with the lowest average daily YLL). Gray shadings represent 95% confidence interval bands of estimates. Note: (A) Phoenix, (B) Philadelphia, and (C) Stockholm.

Stockholm, the coolest site in our study, with a minimum YLL temperature of 27°C (98th percentile), had a sharp increase in YLL at high temperatures, but the association was stronger for temperatures below the fifth percentile () (Figure 4C). Both high and low temperatures were associated with increased YLL during the first two days (Figure S3).

The 14-d cumulative temperature-mortality associations (Figure S5) were also U-shaped, although for Philadelphia, estimated increases of mortality at low temperatures appeared weaker than those of YLL (Figure S5B). In Phoenix, YLL appeared more sensitive to heat than mortality (Figures 4A and S5A), and patterns for YLL and mortality were rather similar in Stockholm (Figures 4C and S5C).

The relationships described above did not substantially change when we varied the degrees of freedom used for the seasonal term, adjusted the knot placements to different temperature percentiles, or used a natural cubic spline basis instead of a b-spline (Figures S6-15). As expected, there were some small differences in the estimated associations on a site-by-site basis with different model parameters, especially the increase in the complexity of the modeled associations when we used more the three knots for the exposure term. However, the general nature of the differences between concerning the association between temperature and YLL remained consistent across the sensitivity analysis.

Discussion

This study investigated the association of daily maximum temperature with YLL at seven low- to high-income sites. We observed an association between high and low temperatures and YLL at many of the seven study sites, contributing to a very limited evidence base regarding the effect of heat and cold on the burden of heat beyond daily death counts.

Increasing knowledge of the role of temperature as a determinant of health, combined with projections of a warmer future (particularly in cities), has elevated public health concerns about the impacts of heat and cold (Hondula et al. 2015a; Huang et al. 2012a; Jones et al. 2015; Sheridan and Allen 2015). Enhancement of the evidence base regarding temperature–health effects that most effectively informs planning and policy in all geographic locations where heat and cold are of concern is an important goal for environmental health researchers.

The seven locations that we examined featured diverse climatic and demographic profiles, and patterns of associations between daily maximum temperature and YLL varied substantially among the sites. In HDSS sites with a hot climate—Nouna and Vadu—heat was associated with YLL, but there was no evidence of an association with cold. For Nouna, the hottest study site, YLL increased steadily beyond the reference temperature (first percentile of the temperature distribution). Patterns of associations were similar for daily temperature and mortality. The observed association was generally consistent with the study by Diboulo et al. (2012), who showed a linear association between mortality and temperature in Nouna. Those authors also reported the largest relative risks of high temperature exposure on young children, which is consistent with our observation of large temperature effects on YLL, given that children have higher life expectancy.

Of all sites, Nairobi showed the strongest association between YLL and temperatures at the fifth percentile. This HDSS site, an African city with relatively mild, subtropical climate, had the lowest estimated YLL at 27°C. In Guangzhou, China, another large, subtropical and rapidly developing city, Yang et al. (2015) estimated an immediate increase in YLL with high temperature and cold effects that lasted for one week. Overall, associations with YLL were larger for low than for high temperatures, particularly for those less than 75 years of age. We found similar results for temperature extremes in Nairobi, although unlike findings in Guangzhou, cold effects appeared to last longer than one week.

Phoenix experienced summer temperatures as high as those observed in the HDSS sites Vadu and Nouna, but much colder winters in which temperatures can fall below 10°C. In this high-income city, low temperature was associated with increased YLL. Estimated cold effects, however, were weaker than heat effects—unlike associations found in Nairobi, where YLL appeared to increase more at the lower end of the temperature range. Heat was also associated with higher YLL in Phoenix. In Brisbane, Australia, Huang et al. (2012a) investigated the association between daily mean temperature and YLL and found evidence consistent with immediate heat effects on YLL, whereas associations with cold were lagged by two days or longer—a pattern usually seen also in weather–mortality studies (Armstrong 2006; Braga et al. 2001; Guo et al. 2014; Pattenden et al. 2003). YLL from cold exposure outnumbered those from heat in Brisbane (Huang et al. 2012b). In Phoenix (the study location most like Brisbane in terms of socioeconomic and climatic features), we observed a U-shaped association, although estimated cold effects occurred over longer lags (12–14 d) than they did in Brisbane.

At the two coldest locations, Stockholm and Philadelphia, heat posed a health risk. High temperature was associated with higher YLL in Philadelphia at shorter lags, as well as in Stockholm, although to a lesser extent. One possible explanation is that summer temperatures rarely exceed a critical level for public health risks in Sweden; in Stockholm, the lowest YLL was observed at the 98th percentile of temperature, making it difficult to detect adverse effects of extreme heat on YLL. Another possible explanation is that the majority of temperature-related deaths in Stockholm occur among individuals age 65 or above, individuals who are more likely to suffer from chronic health conditions, and therefore are more vulnerable to heat effects (Oudin Åström et al. 2015; Åström et al. 2013). This hypothesis is supported by Rocklöv et al. (2011) who found short-term heat-related mortality risks in the very old, but not in younger persons. Many elderly people amounting to fewer YLL in Stockholm resulted in a low number of heat-related YLL, consistent with our results.

Interestingly, in Nouna (the hottest site) and Stockholm (the coldest site), the minimum YLL temperatures were at opposite ends of the spectrum (first and 98th percentile, respectively). However, effects of low temperature could not be estimated for Nouna due to its small population size and low death counts, whereas for Stockholm, heat above the 98th percentile was associated with a substantial increase in YLL.

In Stockholm and Philadelphia, where winter temperatures can reach as low as , YLL increased with cold, but to a lesser degree than in Nairobi. Experimental and observational studies have reported cardiovascular risk factors and diseases to be associated with cold (Mercer 2003; Mercer et al. 1999; Wolf et al. 2009). Wolf et al. (2009) report that the risk of fatal cardiac events in Germany increased with low-temperature exposure even in summer, and associations were stronger during milder years than colder ones. Heightened sensitivity in populations not adapted to low ambient temperature might explain the J- or U-shaped associations we found in warmer locations, such as Phoenix and Nairobi. Furthermore, influenza contributes to excess mortality during winter in temperate regions. In tropical and subtropical countries, however, high humidity during the rainy season rather than low temperature is the main driver of influenza incidence and mortality (Reichert et al. 2004; Tamerius et al. 2013). Although we lack cause-specific data, we assume that cardiovascular diseases are the main causes of cold-related mortality in high-income sites, whereas patterns for influenza and other respiratory diseases might differ between hot and cold regions.

The observed contrasts in estimated health effects across sites might be explained by differences in climate, socioeconomic status, demographic characteristics, and health profiles (although data-related bias cannot be excluded). Nairobi, which has mild winters with the coldest climate of all HDSS sites, was the only HDSS site with a pronounced association of low temperatures with YLL. A similar association was observed at Kisumu, but estimated effects lacked statistical precision. In a previous study, Egondi et al. (2015) described that in Nairobi YLL increased during cold waves, but not during heat waves. In our study, a modest increase in YLL was found at higher temperatures at this site, which hints at heat as the cause of relatively few deaths among people at young ages, resulting in low YLL. Physiological and behavioral adaptation as well as poverty-related factors might explain high vulnerability to cold in Nairobi, in comparison with U.S. and Swedish sites. Nairobi’s urban slum setting is characterized by high infant mortality, and infectious diseases are the dominant cause of death (Streatfield et al. 2014), both of which are related to poverty and adverse living conditions. Our study site in Vadu, on the other hand, is a setting with higher incomes in comparison with the other HDSS sites, and infant mortality is low (Streatfield et al. 2014). At Vadu, we did not find strong evidence of cold effects. Socioeconomic deprivation might therefore reinforce vulnerability to cold, even with modest reduction in temperature (as in Nairobi), making cold an underestimated public health concern in low-income settings. This speculation is supported by Yang’s study in Guangzhou, China, where larger risks were associated with heat and cold among people with lower education and lower income levels, which may be serving as a proxy for behavioral factors, exposure differences (both temperature and air quality), and/or housing quality (Yang et al. 2015).

Xu et al. (2014) described that in addition to heat, short-term variation in temperature was associated with YLL in Brisbane, Australia. This finding hinted at the importance of physiological adaptation to heat and cold and might explain observed differences in the magnitude of effect estimate among our seven study sites, where the local range and temporal fluctuation of temperature varied considerably. Building materials and indoor facilities may also play an important role in these settings to protect occupants indoor from extreme ambient temperatures.

It is likely that varying profiles of cause of death and preexisting morbidity are underlying mechanisms contributing to different temperature–health associations across sites. It is known that adaptation to extreme temperature is impaired in individuals with chronic illnesses, such as coronary heart disease (Kenny et al. 2010; Zanobetti et al. 2012), and death occurs often with short delay (Gasparrini et al. 2015). In Vadu, in contrast with the other HDSS sites of our study, noncommunicable diseases account for more than half of all deaths in the population age 12 or older, and heat-related mortality risk is largest for these causes of death (Ingole et al. 2015). A high number of fragile elderly people dying from noncommunicable diseases might therefore explain the larger heat-related YLL in Vadu, in comparison with such deaths in Nouna and Kisumu. Further inspection of lagged heat effects in Vadu (reduction and stabilization of YLL following the increase in YLL at lag 0) suggests the advancement of deaths in frail individuals.

We observed some variations in the lag structure of estimated effects at locations with a clear cold signal (Nairobi, Phoenix, Philadelphia, and Stockholm). In Philadelphia and Stockholm, associations of low temperature with YLL were larger at short lags and declined afterwards, whereas in Phoenix and Nairobi, YLL appeared to increase only in the second week. Considering that acclimatization and adaptive behavior contribute to vulnerability, subgroups affected by cold may have different characteristics in terms of behavior, physiology, and housing conditions at different locations. A closer look at cause-specific mortality related to low temperature could provide insights into the mechanisms behind differing lag structures.

Strengths and Limitations

To our knowledge, this is the first study to estimate the effects of temperature on YLL across low-, middle-, and high-income countries. Quantifying health outcomes as YLL instead of death counts has the advantage of accounting for individual life expectancy at the time of death, and as such, the magnitude of the dependent variable is disproportionately driven by death among younger people. The seven study sites cover a wide range of climates and populations with specific demographic and socioeconomic features. The INDEPTH HDSS sites provide unique data from countries where collection of regular nationwide health and demographic information is rare. These sites, however, comprise small areas and populations, reducing the occurrence of daily deaths and YLL, which could limit the statistical power and precision of results. The small sample sizes also led to imprecision of estimated life expectancies and individual YLL, adding uncertainty to the analyses. Furthermore, in particular at Vadu HDSS, the heaping day indicator might not fully control for missing dates of birth.

Our analyses accounted for nonlinear effects of temperature, such as the U-shaped association observed in Phoenix. Simultaneously, delayed effects on YLL were modeled (Armstrong 2006; Gasparrini et al. 2010). By applying a distributed lag nonlinear model, we were able to estimate the effects of both low and high temperature and differences in their lag structure to compare across sites.

Kinney (2012) noted that heat victims are likely to be more frail than others in the same age group, which implies that their remaining life expectancy might be much lower. Baccini et al. (2013) showed in a European multisite study that the magnitude of harvesting (advancement of deaths that were likely to occur within days or weeks due to frailty) is considerable, giving support for the hypothesis of differential life expectancies. Our study did not adjust for harvesting. In traditional temperature–mortality studies, failure to adjust for harvesting can result in an overestimate of the public health burden of heat or cold. We expect that the results of our study based on YLL are less sensitive to the displacement phenomena because the deaths of elderly people (who are more likely to have had deaths advanced forward in time by heat or cold by days to weeks due to higher frailty) are given less weight in calculating the total burden than the deaths of younger people. However, we recommend additional exploration of how the controlling for mortality displacement and calculations based on YLL might overlap.

For all analyses, daily maximum temperature was used as the exposure variable. The extent to which exposure-variable choice affects temperature–mortality relationships has been the subject of investigation of several recent studies in high-income countries (Barnett et al. 2010; Rodopoulou et al. 2015); sensitivity to exposure variable appears to vary on a city-by-city basis. It is possible in some of our locations that the true cold effect, especially in winter and at locations with large intradiurnal temperature variability, might have been underestimated by the use of daily maximum temperature. We expect that the general pattern of our results would hold with the choice of a different exposure variable but could not test their sensitivity because of the lack of data availability at all sites. We also examined the sensitivity of our conclusions to different modeling assumptions concerning degrees of freedom and the nature of the spline function used to represent the exposure–response relationship. By holding model parameters consistent across sites, we likely used site-specific models that were not individually optimized; this choice facilitated interpretation for this first cross-regional comparison of temperature–YLL effects. However, we encourage more research refining the appropriate statistical and theoretical frameworks to understand better how and why the exposure–response relationship varies from place to place.

Policy implications/recommendations.

Individual and environmental factors such as adaptive behavior, gender, socioeconomic status, access to health care, and building structure have been shown to contribute to the adverse health effects of heat (Dutta et al. 2015; Eisenman et al. 2016; Hondula et al. 2012; Laverdière et al. 2016; Näyhä et al. 2016; Smith et al. 2014; Son et al. 2016; Vandentorren et al. 2004). We still know little, however, about the role of these factors in exacerbating vulnerability to cold in low- and middle-income countries, where low temperature usually is not given priority as a health risk.

We found evidence of an effect of low and high temperature on YLL, an indicator of premature mortality, in climatically and socioeconomically diverse locations. By providing policy makers with an absolute number of YLL as a measure of health burden, the results are readily useful to policy makers for setting priorities to reduce heat and cold exposure of individuals in their populations. Many individual and environmental factors contributing to increased vulnerability from heat and cold can certainly be affected by simple measures (Delgado Cortez 2009; Dodman et al. 2010; WHO 2012). Apart from poverty reduction, interventions targeting improved housing conditions and awareness campaigns may help to reduce temperature-related health risks.

In high-income locations (United States and Sweden), it is important to increase public health awareness about the risks associated with temperature extremes. Even for the hot desert city of Phoenix, we found evidence of an adverse effect of low temperature, a risk factor that might be underrated by residents and health care providers. We also observed associations between low temperature and YLL in Philadelphia and Stockholm, where cold winters are common. In these locations, efforts to prevent cold-related deaths are likely already underway and should continue to be supported.

Conclusions

This study investigated the effect of temperature on daily YLL in seven locations with different socio-economic and climatic characteristics. Both heat and cold were associated with YLL in most locations, although effect sizes and lag structures varied. Our findings, using YLL as the outcome, are in line with studies showing weather effects on mortality in which the dependent variable is the number of daily deaths. Insights from this study can help to identify places and populations where burden of heat and cold is highest from a YLL perspective, which we hope can guide the improvement of local, national, and global-scale adaptation and mitigation policies and programs.

Click here for additional data file.