The influence of temperature on mortality and its Lag effect: a study in four Chinese cities with different latitudes
BMC Public Health volume 16, Article number: 375 (2016)
Global climate change is one of the most serious environmental issues faced by humanity, and the resultant change in frequency and intensity of heat waves and cold spells could increase mortality. The influence of temperature on human health could be immediate or delayed. Latitude, relative humidity, and air pollution may influence the temperature–mortality relationship. We studied the influence of temperature on mortality and its lag effect in four Chinese cities with a range of latitudes over 2008–2011, adjusting for relative humidity and air pollution.
We recorded the city-specific distributions of temperature and mortality by month and adopted a Poisson regression model combined with a distributed lag nonlinear model to investigate the lag effect of temperature on mortality.
We found that the coldest months in the study area are December through March and the hottest months are June through September. The ratios of deaths during cold months to hot months were 1.43, 1.54, 1.37 and 1.12 for the cities of Wuhan, Changsha, Guilin and Haikou, respectively. The effects of extremely high temperatures generally persisted for 3 days, whereas the risk of extremely low temperatures could persist for 21 days. Compared with the optimum temperature of each city, at a lag of 21 days, the relative risks (95 % confidence interval) of extreme cold temperatures were 4.78 (3.63, 6.29), 2.38 (1.35, 4.19), 2.62 (1.15, 5.95) and 2.62 (1.44, 4.79) for Wuhan, Changsha, Guilin and Haikou, respectively. The respective risks were 1.35 (1.18, 1.55), 1.19 (0.96, 1.48), 1.22 (0.82, 1.82) and 2.47 (1.61, 3.78) for extreme hot temperatures, at a lag of 3 days.
Temperature–mortality relationships vary among cities at different latitudes. Local governments should establish regional prevention and protection measures to more effectively confront and adapt to local climate change. The effects of hot temperatures predominantly occur over the short term, whereas those of cold temperatures can persist for an extended number of days.
Global climate change is one of the most serious environmental issues confronting humanity. The Intergovernmental Panel on Climate Change (IPCC) reported that average temperature data showed a global warming trend of 0.85 (0.65 to 1.06) °C from 1880 to 2012. The IPCC projected that there would be an increase in the duration, intensity and spatial extent of heat waves and warm spells during the next decade . Many studies have found that heat waves and cold spells could lead to an increase in the number of temperature-related deaths [2, 3]. Temperature–morality curves have been described as U-, V-, or W-shaped, with U- or V-shaped representing the peaks of deaths occurring at very high and low temperatures; the U-shaped curve has a more pronounced decrease and increase in mortality with a rise in temperature than the V-shaped curve. The W-shaped curve represents the fact that increased mortality levels are distributed at high and low temperatures. There is a small peak of deaths between high and low temperatures, which may result from the response to a rise in temperature. People may ignore measures necessary to adapt to temperature change, resulting in a slight increase in mortality at the beginning, which then decreases after some adaptation measures are taken [4–8]. Previous studies have found that the adverse effects of temperatures on human health not only occur in the short term but also after a delay [9, 10]. A distributed lag nonlinear model (DLNM) has been developed to study these delayed effects and the nonlinear exposure–response relationship simultaneously .
The influence of temperature on human health varies among countries and regions . For more effective adaptation to climate change, governments should establish regional prevention and control systems based on the corresponding influences of temperature on human health. Many studies have shown that latitude exerts an important influence on the relationship between temperature and mortality. Generally, people living at low latitudes have a relatively weak adaptive capacity to cold and a strong adaptive capacity to heat, and the reverse has also been shown [5, 13]. Additionally, humidity [14, 15] and air pollution [16, 17] have been reported to have an effect on mortality rates.
We studied the relationship between temperature and mortality in four Chinese cities with different latitudes, and adjusted for the possible influence of relative humidity and air pollution. We used a quasi-Poisson analysis to investigate the relationship between temperature and mortality and a DLNM to study the delayed effect of temperature on mortality. Our findings will provide useful information for understanding the health effects of temperature at different latitudes in China.
We studied four cities with different latitudes (Wuhan, Changsha, Guilin, and Haikou, as shown in Fig. 1). These cities are the provincial capitals of their respective provinces in China, except for Guilin, which is a famous tourist destination in Guangxi. Wuhan has a subtropical monsoon humid climate; the city is located where the world’s third longest river, the Yangtze, joins its largest tributary, the Han River. Wuhan is an important transportation hub and a scientific research and education center. The city has four distinct seasons, with cold and wet winters and hot, humid summers. Changsha has a subtropical monsoon climate, with four distinct seasons including cold winters and hot summers; it is a major and highly industrialized city in China. The climate in Guilin also has subtropical monsoons but is relatively mild, with neither extreme cold in winter nor excessive heat in summers. Haikou has a tropical marine climate, with hot and rainy summers and mild winters.
Owing to data availability, we were only able to use data from two districts in Wuhan (Qiaokou and Jiangan) and one district each in Changsha (Tianxin), Guilin (Xiufeng), and Haikou (Meilan); according to the Sixth National Population Census of China, in 2010 the populations in these four districts were 1.73 million, 0.47 million, 0.16 million and 0.62 million, respectively. The regional gross domestic product was 551.58, 454.70, 110.86 and 59.06 billion CNY, respectively, and per capita annual disposable income was 20.81, 22.81, 17.95 and 16.72 thousand CNY, respectively .
Meteorological data of the four cities for 2008–2011, including daily mean temperature, minimum temperature, maximum temperature, and relative humidity, were obtained from the China Meteorological Data Sharing Service System .
Daily mortality data of the four cities were collected from the Provincial Centers for Disease Control and Prevention. Causes of death were encoded according to the 10th International Classification of Diseases (ICD-10). The term “non-accidental death” referred to all causes of death, excluding external causes such as injury and poisoning (ICD-10 codes: A00-R99). The Air Pollution Index (API) is an indicator of the ambient air quality, which is calculated using the measured concentrations of sulfur dioxide (SO2), nitrogen dioxide (NO2) and respirable particulate matter (PM10), according to the methodology of the World Health Organization and European Centre for Environment and Health, at all monitoring stations throughout each city-specific territory. A lower API indicates better air quality . API data were collected from the Ministry of Environmental Protection of the People’s Republic of China.
Relationships of mortality and certain risk factors are typically investigated using Poisson regression [21, 22]. However, some researchers have hypothesized that quasi-Poisson analysis should be used if overdispersion is found with Poisson regression, and that the R qcc package provides a test for this phenomenon [23, 24]. We used a DLNM to assess the effects of daily temperature on mortality with different lag days [10, 11]. To explore the effects of temperature on daily mortality, factors such as long-term and seasonal trends, day of the week (DOW), relative humidity (RH), and API were introduced into the model as covariates.
Describes the model we selected, where E[Y t ] denotes the expected number of nonaccidental deaths on day t; α is the intercept; T t,l is a matrix obtained by applying the temperatures, with l referring to the lag days and β the coefficient. NS represents the natural spline function; RH t is the daily relative humidity on day t; the degrees of freedom (df) is 3, and the df of each variant is confirmed by the Akaike information criterion for the quasi-Poisson models . Time represents the long-term temporal trend; Season represents the seasonal trend; API t is the air pollution index on day t; DOW t is the day of the week on day t.
We used lag-stratified natural cubic spline models to explore the nonlinear and delayed temperature–mortality association. We adopted city-specific optimum temperature as the reference value, which was determined by comparing the relative risks (RR) of temperatures at 3-day lag, and the temperature with the lowest RR was set as the optimum temperature. A natural cubic spline with 4 df was applied to the daily mean temperature. Knots of the mean temperatures were placed at equally spaced quantiles. Lag stratification was defined as 5 df for the mean temperature. The knots of the lag calculation were set at equally spaced values on the log scale of the lags. We plotted RR against the temperature and lags to show the entire relationship between mean temperature and mortality. Extreme cold temperatures were set as the first percentile of the daily mean temperature, and the lag days were set to 21. Extreme hot temperatures were set as the 99th percentile of daily mean temperature, and the lag days were set to 3.
The sensitivity analyses were performed by changing the df of long-term trends and API, with df (time) = 5, 6, 7 or 8/year, and df (API) = 2, 3, 4, or 5. We used R 3.0.2, software for the analyses. The distributed lag nonlinear models were specified with the “dlnm” package, and tests for over dispersion were performed using the “qcc” package.
Table 1 shows the city-specific distribution of temperature, relative humidity, API, and daily mortality counts. The temperature of the four cities showed an increasing trend from north to south, whereas the API had a decreasing trend. The relative humidity of Wuhan and Haikou were higher than those of the other two cities.
Figure 2 shows the city-specific distribution of daily mean temperature by month. We found that the primary cold months were December, January, February and March, and the hot months were June through September. Figure 3 shows that most deaths occurred in the cold months, and ratios of deaths in the cold months to hot months were 1.43, 1.54, 1.37 and 1.12 for Wuhan, Changsha, Guilin and Haikou, respectively. The risk ratios were higher in Wuhan and Changsha than in Guilin and Haikou.
Contour plots (Fig. 4) show the relationships between mean temperature and mortality at a lag of 21 days. For extremely high temperatures, the effect was the strongest at lag0 and then declined rapidly. Its effect was generally retained for three days at Wuhan, Changsha and Guilin, whereas for Haikou, the effect of hot temperatures could be maintained longer. For the extremely low temperatures, the risk was the smallest at lag0, followed by an incremental increase, except at Guilin, where the low temperature remained at a relatively low risk level within a lag of 21 days.
We drew city-specific temperature–mortality curves for different lag times (Fig. 5). The curves were similarly V-shaped at a lag of 3 days, which may indicate the short-term effects of temperature on mortality. Then, with increasing lag time, the risk of cold temperatures increased gradually, whereas the increase in risk of hot temperatures was much smaller. Table 2 shows the RR of extreme cold and hot temperatures compared with the optimum temperatures, which were 22.9, 26.9, 24.5 and 15.0 °C for Wuhan, Changsha, Guilin and Haikou. RRs (95 % confidence interval) of the extreme cold temperatures were 4.78 (3.63, 6.29), 2.38 (1.35, 4.19), 2.62 (1.15, 5.95) and 2.62 (1.44, 4.79), respectively, at a lag of 21 days. RRs of the extreme hot temperatures were 1.35 (1.18, 1.55), 1.19 (0.96, 1.48), 1.22 (0.82, 1.82) and 2.47 (1.61, 3.78), respectively, at a lag of 3 days.
Sensitivity analyses showed that with changing the df of time trends and API, the extreme cold and hot effects were robust to alternative models (Additional file 1: Table S1).
We found that mortality was higher in the cold months than in the hot months, and that the seasonal pattern was lower for mortality in the southern cities in China than in the central cities. This finding concurs with a number of previous studies [26–29]. The seasonality of mortality rates might be attributed to the following reasons: exposure to cold increases energy expenditure, peripheral vasoconstriction and cardiac afterload, which can trigger myocardial ischemia and stroke [26, 30]; influenza is more prevalent in winter ; and differences in lifestyle factors such as dietary habits  and physical activity  between cold and hot months. Considering global climate change, some researchers have suggested that an increase of heat-related mortality would be somewhat offset by reductions in cold-related mortality, and that even overall mortality might decrease for a certain period [34, 35], but other researchers have disagreed [36, 37]. It is possible that in the future, mortality might be reduced with winter-dominant warming and increased with pronounced summer warming .
We found that high and low temperatures were associated with increased mortality in the four studied cities, and the effect of high temperatures on mortality was predominantly within 3 days. The effect of low temperatures could persist for 21 days. This finding is consistent with some studies [10, 38]. For Haikou, the effect of high temperatures could last for an extended period, possibly because the temperature is persistently high in summer. We only examined a small sample, and this phenomenon requires further research.
We found that at a lag of 3 days, which reflects short-term effects of the temperature on mortality rate, the curve of their relationship had a V shape. As lag time increased, the shape of the curves changed, and RRs associated with low temperatures increased gradually, whereas the risks associated with high temperatures did not substantially increase but actually decreased. At a lag of 21 days, the temperature–mortality curve appeared to be U-shaped for Wuhan and Changsha. For Haikou, the curve appeared to have a W shape, which is consistent with a previous study in China by Jinan . For Guilin, the mortality risk decreased with temperature after a lag of 21 days, whereas the curve was V-shaped at a lag of 3 days. This might be because summer temperatures are not excessively high in Guilin; the effect of high temperatures typically persists for 3 days and the mortality effect of those temperatures after that is weak, so the effect estimated by the model may not have been exact. A similar phenomenon was found in Charlotte, a city in the southeastern United States of America . The study sample in Guilin was relatively small, so the generalizability of this result is limited. Because of the long lag effect of cold temperatures and relatively short lag effect of hot temperatures, we believed that the effect of extreme cold temperatures could be calculated at a lag of 21 days, whereas the effect of extreme hot temperatures could be calculated at a lag of 3 days. This method was also used in a previous study .
The estimated effects of extreme cold and hot temperatures on mortality varied by city. This might be because the four cities have different temperature distributions (Fig. 2). Generally, the cumulative effects of extreme cold were stronger than extreme heat, in agreement with the seasonal pattern discussed above. This phenomenon was also found in a previous study in China . The greatest effect of extreme cold was in Wuhan, where the mean temperature was the coldest of the four cities. The weakest effect was in Haikou, where both the mean and minimum temperatures were the highest of the four cities. The greatest effect of extreme heat was also in Haikou.
Our study has the following limitations. First, because of the limited availability of data, we only included cities in central and southern China, ignoring northern cities. We will try to correct this limitation in a future study. Second, samples were confined to one or two districts in the studied cities, restricting the generalizability of our results. Finally, socio-demographic factors, including age, sex, and income, were not considered. Those factors may influence the temperature–mortality relationship.
More deaths occurred during cold months in the studied cities, and we found that cities at different latitudes have diverse temperature–mortality relationships. The effects of hot temperatures are predominantly sudden, whereas those of cold temperatures can persist for an extended period. Local governments should establish regional prevention and protection measures to more effectively confront and adapt to climate change.
Ethics approval and consent to participate
Ethical approval was given by the medical ethics committee of Wuhan University School of Medicine with the following reference number: YXBGW0016.
Consent for publication
The mortality data used in this study is the number of deaths per day, and our paper does not contain any individual personal’s data in any form.
Availability of data and materials
Because that we have signed nondisclosure agreements with Chinese Center for Disease Control and Prevention, so the data used in this paper will not be shared.
air pollution index
degrees of freedom
distributed lag non-linear model
day of the week
10th International Classification of Diseases
natural cubic spline
Intergovernmental Panel on Climate Change. Climate Change 2013: the Physical Science Basis. Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 2013; 37–89.
Anderson BG, Bell ML. Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009;20:205–13.
Yu WW, Mengersen K, Wang XY, Ye XF, Guo YM, Pan XC, et al. Daily average temperature and mortality among the elderly: a meta-analysis and systematic review of epidemiological evidence. Int J Biometeorol. 2012;56:569–81.
Yang J, Ou CQ, Ding Y, Zhou YX, Chen PY. Daily temperature and mortality: a study of distributed lag non-linear effect and effect modification in Guangzhou. Environ Health-Glob. 2012;11:63.
Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA. Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002;155:80–7.
Burkart K, Schneider A, Breitner S, Khan MH, Kramer A, Endlicher W. The effect of atmospheric thermal conditions and urban thermal pollution on all-cause and cardiovascular mortality in Bangladesh. Environ Pollut. 2011;159:2035–43.
Li M, Zhou M, Zhang X, Huang J, Bai L, Sang S, et al. Impact of temperature on non-accidental deaths and cause-specific mortality in four districts of Jinan. Zhonghua Liu Xing Bing Xue Za Zhi. 2014;35:684–8.
Guo Y, Gasparrini A, Armstrong B, Li S, Tawatsupa B, Tobias A, et al. Global variation in the effects of ambient temperature on mortality: a systematic evaluation. Epidemiology. 2014;25:781–9.
Ha J, Shin Y, Kim H. Distributed lag effects in the relationship between temperature and mortality in three major cities in South Korea. Sci Total Environ. 2011;409:3274–80.
Wu W, Xiao Y, Li G, Zeng W, Lin H, Rutherford S, et al. Temperature-mortality relationship in four subtropical chinese cities: a time-series study using a distributed lag non-linear model. Sci Total Environ. 2013;449:355–62.
Gasparrini A. Distributed lag linear and non-linear models in R: the package dlnm. J Stat Softw. 2011;43:1–20.
Hajat S, Kosatky T. Heat-related mortality: a review and exploration of heterogeneity. J Epidemiol Commun H. 2010;64:753–60.
Hu MJ, Ma WJ, Zhang YH, Liu T, Lin HL, Luo Y, et al. Relationship between temperature and the risks of mortality in China: a meta-analysis. Zhonghua Liu Xing Bing Xue Za Zhi. 2013;34:922–26.
Fernandez-Raga M, Tomas C, Fraile R. Human mortality seasonality in Castile-Leon, Spain, between 1980 and 1998: the influence of temperature, pressure and humidity. Int J Biometeorol. 2010;54:379–92.
Barreca AI. Climate change, humidity, and mortality in the United States. J Environ Econ Manage. 2012;63:19–34.
Dholakia HH, Bhadra D, Garg A. Short term association between ambient air pollution and mortality and modification by temperature in five indian cities. Atmos Environ. 2014;99:168–74.
Guo YM, Li SS, Tawatsupa B, Punnasiri K, Jaakkola JJK, Williams G. The association between air pollution and mortality in Thailand. Sci Rep-Uk. 2014;4:5509.
National Bureau of statistics of China. China statistic year book 2010. Beijing: China Statistics Press; 2010.
China Meteorological Data Sharing Service System. 2016. http://data.cma.cn/data/index/6d1b5efbdcbf9a58.html. Accessed 30 April 2016.
Luo YP, Chen H, Zhu QA, Peng CH, Yang G, Yang YZ, et al. Relationship between Air pollutants and economic development of the provincial capital cities in China during the past decade. PLoS One. 2014;9:104013.
Ramis R, Vidal E, Garcia-Perez J, Lope V, Aragones N, Perez-Gomez B, et al. Study of non-Hodgkin’s lymphoma mortality associated with industrial pollution in Spain, using poisson models. BMC Public Health. 2009;9:26.
Frost G, Harding AH, Darnton A, McElvenny D, Morgan D. Occupational exposure to asbestos and mortality among asbestos removal workers: a poisson regression analysis. Br J Cancer. 2008;99:822–29.
Kabacoff RI. R in action: data analysis and graphics with R. Shelter Island: Manning Publications; 2011. p. 327–8.
Ver Hoef JM, Boveng PL. Quasi-poisson vs. negative binomial regression: how should we model overdispersed count data? Ecology. 2007;88:2766–72.
Guo Y, Barnett AG, Pan X, Yu W, Tong S. The impact of temperature on mortality in Tianjin, China: a case–crossover design with a distributed lag nonlinear model. Environ Health Persp. 2011;119:1719–25.
Marti–Soler H, Gonseth S, Gubelmann C, Stringhini S, Bovet P, Chen PC, et al. Seasonal variation of overall and cardiovascular mortality: a study in 19 countries from different geographic locations. PLoS One. 2014;9:113500.
Douglas S, Rawles J. Latitude–related changes in the amplitude of annual mortality rhythm. The biological equator in man. Chronobiol Int. 1999;16:199–212.
Huang CR, Chu C, Wang XM, Barnett AG. Unusually cold and dry winters increase mortality in Australia. Environ Res. 2015;136:1–7.
Kendrovski VT. The impact of ambient temperature on mortality among the urban population in Skopje, Macedonia during the period 1996–2000. BMC Public Health. 2006;6:44.
Wang Y, Levi CR, Attia JR, D’Este CA, Spratt N, Fisher J. Seasonal variation in stroke in the hunter region, Australia - a 5-year hospital–based study, 1995–2000. Stroke. 2003;34:1144–50.
Madjid M, Miller CC, Zarubaev VV, Marinich IG, Kiselev OI, Lobzin YV, et al. Influenza epidemics and acute respiratory disease activity are associated with a surge in autopsy-confirmed coronary heart disease death: results from 8 years of autopsies in 34 892 subjects. Eur Heart J. 2007;28:1205–10.
Faure H, Preziosi P, Roussel AM, Bertrais S, Galan P, Hercberg S, et al. Factors influencing blood concentration of retinol, alpha-tocopherol, vitamin C, and beta-carotene in the participants of the SU.VI.MAX trial. Eur J Clin Nutr. 2006;60:706–17.
Matthews CE, Freedson PS, Hebert JR, Stanek EJ, Merriam PA, Rosal MC, et al. Seasonal variation in household, occupational, and leisure time physical activity: longitudinal analyses from the seasonal variation of blood cholesterol study. Am J Epidemiol. 2001;153:172–83.
Guest CS, Willson K, Woodward AJ, Hennessy K, Kalkstein LS, Skinner C, et al. Climate and mortality in Australia: retrospective study, 1979–1990, and predicted impacts in five major cities in 2030. Climate Res. 1999;13:1–15.
Davis RE, Knappenberger PC, Michaels PJ, Novicoff WM. Seasonality of climate-human mortality relationships in US cities and impacts of climate change. Climate Res. 2004;26:61–76.
Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, et al. Managing the health effects of climate change. Lancet. 2009;373:1693–733.
Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. Impact of regional climate change on human health. Nature. 2005;438:310–17.
Lin YK, Ho TJ, Wang YC. Mortality risk associated with temperature and prolonged temperature extremes in elderly populations in Taiwan. Environ Res. 2011;111:1156–63.
Wang CC, Chen RJ, Kuang XY, Duan XL, Kan HD. Temperature and daily mortality in Suzhou, China: a time series analysis. Sci Total Environ. 2014;466:985–90.
The authors would like to thank Dr. Songbo Hu for his assistance with the data analysis in this work and to the participants who contributed to this research.
This study was supported by the National Basic Research Program of China (973 Program) (Grant No. 2012CB955504), the Fundamental Research Funds for the Central Universities (Grant No. 2015305020201) and the Modeling and forewarning the spatial-temporally distribution of the spread of Oncomelania hupensis in Jianghan plain, Central China (Grant No. 41571487).
The authors declare that they have no competing interests.
CHY, XDL and JZB conceived of the study and participated in its design. JZB and ZKW performed the experiments, analyzed the data, and wrote the manuscript. CHY and XDL reviewed the manuscript for major intellectual content. All authors read and approved the final manuscript.
About this article
Cite this article
Bao, J., Wang, Z., Yu, C. et al. The influence of temperature on mortality and its Lag effect: a study in four Chinese cities with different latitudes. BMC Public Health 16, 375 (2016). https://doi.org/10.1186/s12889-016-3031-z