- Open Access
- Open Peer Review
Daily temperature and mortality: a study of distributed lag non-linear effect and effect modification in Guangzhou
Environmental Healthvolume 11, Article number: 63 (2012)
Although many studies have documented health effects of ambient temperature, little evidence is available in subtropical or tropical regions, and effect modifiers remain uncertain. We examined the effects of daily mean temperature on mortality and effect modification in the subtropical city of Guangzhou, China.
A Poisson regression model combined with distributed lag non-linear model was applied to assess the non-linear and lag patterns of the association between daily mean temperature and mortality from 2003 to 2007 in Guangzhou. The case-only approach was used to determine whether the effect of temperature was modified by individual characteristics, including sex, age, educational attainment and occupation class.
Hot effect was immediate and limited to the first 5 days, with an overall increase of 15.46% (95% confidence interval: 10.05% to 20.87%) in mortality risk comparing the 99th and the 90th percentile temperature. Cold effect persisted for approximately 12 days, with a 20.39% (11.78% to 29.01%) increase in risk comparing the first and the 10th percentile temperature. The effects were especially remarkable for cardiovascular and respiratory mortality. The effects of both hot and cold temperatures were greater among the elderly. Females suffered more from hot-associated mortality than males. We also found significant effect modification by educational attainment and occupation class.
There are significant mortality effects of hot and cold temperatures in Guangzhou. The elderly, females and subjects with low socioeconomic status have been identified as especially vulnerable to the effect of ambient temperatures.
The impact of weather on mortality has become public health significance, especially in light of climate change and rising frequency of adverse weather events (e.g., heat waves and floods). The relationship between temperature and mortality are extensively documented, with excess morality during days with extremely low or high ambient temperatures [1–6]. Substantive evidence has shown a delay between changes in daily temperature and changes in mortality, while the lag period considered are not consistent in different studies. Previous studies usually examined the effects for a single lag [4, 7], and some studies chose a priori a lag period and assessed a cumulative effect [3, 5]. To date, there is lack of a criterion for selecting the optimal lag. Exploring lag distribution of effects may provide some information for selecting an appropriate time frame when assessing temperature effects.
It has been shown that the association between temperature and mortality varied greatly by climate, geographic regions and populations [6, 8, 9], which indicates that it is necessary to assess the impacts of temperature in various regions. The relationship of temperature and mortality had been extensively studied in Europe [1, 10, 11] and the United States [3, 4, 6, 12]. Only a limited number of studies were conducted in developing countries or tropical/subtropical regions [7, 13–16]. More studies on health effects of weather in these regions are required to fully assess global impact of climate change and guide local public health policy.
To develop public health policies that protect those persons most vulnerable to extreme temperatures, researchers have identified factors that confer susceptibility. Previous studies indicated that black, the elderly and females were especially susceptible populations [8, 17, 18], while vulnerability by socioeconomic factors remains unclear. A few studies suggested that those with less education were at higher risk of temperature-related mortality [7, 19, 20], but some other studies showed little or no evidence for effect modification by area-based [21, 22] or individual-based measures of education level [23, 24]. Researchers have called for further research on the role of education level and other socio-economic measures to understand inequalities in health impact of ambient temperatures.
Guangzhou is the largest metropolis in Southern China with the latitude of 23° 7′ N. It is in a typical subtropical climate with mild winter and hot summer. In the present study, we sought to examine the association between daily ambient temperature and mortality in Guangzhou and identify the effect modification of temperature by individual characteristics, including age, sex, educational attainment and occupation class.
Health and environmental data
The Guangzhou Bureau of Health provided individual information for 112,280 deaths from 1 January 2003 to 31 December 2007, including date of birth, date of death, cause of death, sex, educational attainment and occupation. The underlying cause of death was classified by the Tenth Revision of the International Classification of Diseases (ICD-10). We considered non-accidental mortality (ICD-10: A00-R99), mortality due to cardiovascular diseases (I00-I99) and respiratory diseases (J00-J99), and three subcategories, including chronic obstructive pulmonary diseases (J40-J47), stroke (I60-I69) and ischemic heart diseases (I20-I25). Daily number of deaths was summarized by the underlying cause of death, sex, age group (0–64, 65–74, 75–84, 85 years old or above), educational attainment and occupation class, respectively. Educational attainment was defined as the highest degree of schooling completed before death. In the present study, it was classified into three groups: none (i.e. illiterate or semiliterate), primary education and secondary or higher education. Occupation was classified into unemployed (including housewife), blue-collar workers and white-collar workers.
We obtained the following meteorological data from China Meteorological Data Sharing Service System : daily values of minimum, mean and maximum temperature, relative humidity and barometric pressure in Guangzhou. Weather data were collected from the only basic weather station in Guangzhou, Wushan Station. The Guangzhou Bureau of Environmental Protection provided air pollution data for 24-hour average concentrations of particulate matter with diameters less than 10 μm (PM10), nitrogen dioxide (NO2) and sulphur dioxide (SO2). These criteria pollutants have been associated with cardio-respiratory mortality . The Ethics Committee of Southern Medical University where this study was conducted has approved the study proposal.
We fitted a distributed lag non-linear model (DLNM) to assess the association between daily number of deaths and ambient temperature. DLNM, proposed recently by Gasparrini et al.  is a flexible model to describe simultaneously a non-linear exposure-response relationship and delayed effect. It is appropriate to evaluate the characteristic of temperature-associated mortality in that temperature and mortality has shown a J-, W-, V- or U-shaped association [27–29]. A Poisson regression with quasi-Poisson function for daily counts of deaths was constructed, which was specified as
Log[E(Y t )] = α + NS(Time,5*7) + NS(Hum t ,3) + NS(Press t ,3) + NS(PM 10t ,3) + NS(SO 2t ,3) + NS(NO 2t ,3) + γDow t + νHoliday t + βTemp t , l
where Y t is the observed daily deaths at day t (t = 1,2,3…1826); α is the intercept; NS(.) means a natural cubic spline; 7 degrees of freedom (df) per year for time and 3 df for relative humidity (Hum), barometric pressure (Press), PM 10 , NO 2 and SO 2 at the current day were recommended by several previous studies [13, 16, 27]. Dichotomous variables indicating day of the week (Dow) and public holidays (Holiday) are also included in the model. Temp t l is a matrix produced by DLNM to model non-linear and distributed lag effects of ambient temperature over the current day (lag 0) to lag l days, and β is vector of coefficients for Temp t,l ; The maximum lag l was set to 25 days to explore the lag structure of temperature effect . The median value of temperature was used as the reference value to calculate the relative risks. Akaike’s Information Criterion for quasi-Poisson (Q-AIC) was used to choose the df for temperature and lag [13, 30]. The final composition of the function was a natural cubic spline of temperature with 5 df and a natural cubic spline with 5 df for lag days.
We fitted a model for each combination of ambient temperature measures (daily minimum, mean and maximum temperature) and mortality types using the above steps. The Q-AIC was used to choose the temperature measure that best predicted mortality. Mean temperature generally gave the lowest Q-AIC values based on our data Additional file (see 1), and mean temperature provides more easily interpreted results in a policy context and is more familiar to the public. Therefore, we chose mean temperature for subsequent analyses. We calculated percentage change in mortality risk comparing the first to the 10th percentile (cold effect) and the 99th to the 90th percentile (hot effect) of daily mean temperatures .
In order to identify subpopulations that are more susceptible to the effects of cold and hot temperatures, we conducted stratified analyses by sex, age group, educational attainment and occupation class. Furthermore, we used the case-only approach to determine the statistical significance of effect modification by testing the interaction between temperature and individual characteristics of interest. The case-only approach only examines cases. A logistic regression model was constructed with binary variable of individual characteristic of deaths as the dependent variable and mean temperature as the independent variable. The case-only approach has some apparent advantages over conventional time-series methods in assessing the interaction between a time-dependent variable and a time-fixed factor, such as practical simplification of modeling and reduction of potential time-invariant confounders related to death (e.g., smoking) [31, 32]. In previous case-only studies, a dummy variable (e.g., an indicator of extreme temperature) or a linear term of temperature was included in the model [32, 33]. However, the interaction between temperature and individual characteristics may be non-linear since the main effect of temperature on mortality is typically non-linear [4, 8, 27, 34]. In the present study, we applied the distributed lag non-linear model to the case-only approach, which can be more flexible to determine whether individual characteristics modify the non-linear effect of temperature. The seasonal pattern of mortality was captured by including a sin and cosine term with a 365.24-day period in the model [32, 33].
All statistical analyses and modeling were completed in R version 2.13.1 .
During the study period from January 1, 2003 to December 31, 2007, average daily maximum temperature was 27.6°C, mean temperature 23.0°C, minimum temperature 19.8°C in Guangzhou. The 1st, 10th, 90th and 99th percentile of daily mean temperatures was 8.2°C, 13.8°C, 29.9°C and 32.0°C, respectively. The mean concentrations of PM10, NO2 and SO2 were 88.2 μg/m3, 73.2 μg/m3 and 59.3 μg/m3, respectively. Among the total of 112,280 non-accidental deaths, 43,393 persons (38.7%) died from cardiovascular diseases and 21,071 (18.8%) from respiratory diseases. The descriptive statistics of daily number of deaths can be found in Table 1.
The distributed non-linear lag surface revealed a non-linear relationship between temperature and mortality, with higher mortality risk at hot and cold temperatures (Figure 1). Figure 2 shows lag structures of temperature effect over the same day to 25 days previous. We observed immediate effects of cold and hot temperatures on mortality with the strongest and statistically significant effects on the current day (lag 0). Hot effects usually disappeared after approximately 4–5 days. Cold effects could persist for 10–12 days. The initial increase in mortality risk related to hot temperatures was followed by a decrease, consistent with short-term mortality displacement (i.e., harvesting), while no apparent harvesting was found for cold effects (Figure 2).
Based on the lag structures of hot and cold effects, we presented cumulative effects of hot temperatures at lag 0–5 days and cumulative cold effects at lag 0–12 days. The dose–response curve showed that both cold and hot effects seemed to be non-linear for some mortality categories (Figure 3). Therefore, we estimated relative effects of temperatures based on specific portions of the curve. Hot temperatures were associated with a 15.46% (95% confidence interval: 10.05% to 20.87%) increase in non-accidental mortality risk comparing the 99th to the 90th percentile of daily mean temperatures. Cold temperatures were associated with a 20.39% (11.78% to 29.01%) increase in non-accidental mortality comparing the first to the 10th percentile of temperatures. Stronger association was found with mortality due to cardiovascular and respiratory diseases, particularly with mortality due to ischemic heart diseases. Positive but non-significant effects of cold and hot temperatures were found for all other mortality (Table 2).
To explore the potential confounding of air pollution, we performed a sensitivity analysis with adding one air pollution variable at a time or excluding all pollution variables from the model. There was only slight change in the estimate of temperature effect compared to the model controlling for three air pollution variables together Additional file (see 2), indicating that the confounding by air pollution was very small if there is any.
We evaluated the effects of temperatures on non-accidental mortality for subpopulations (Table 3) and further examined the statistical significance of effect modification by individual characteristics (Table 4). Hot effects were significantly larger for females than males. By contrary, males were at higher risk of cold effects compared to females but the difference was non-significant. The point estimates of temperature effects generally increased with age. There was statistically significant difference between the oldest age group (≥85 years old) and the youngest age group (<65 years old). A trend of decreased temperature-related mortality risk with increased education level was observed. No education conferred significant susceptibility. As to occupation class, blue-collar workers suffered more from hot and cold effects compared to white-collar worker, and the difference were statistically significant.
To quantify mortality effects of ambient temperature, linear threshold regression models were often used previously, which assumed linear or log-linear increase below and above the threshold. In fact, the dose–response curve of temperature and mortality is hardly linear below or above the threshold, so it is inappropriate to use linear threshold models to estimate the effects directly. Recently Gasparrini and Armstrong first rigorously developed DLNMs [27, 28]. With unifying many of the previous methods in one unique framework, DLNMs are flexible enough to describe non-linear dependencies and delayed effects of exposure at the same time. Gasparrini et al.  speculated in their article that the DLNMs could be easily translated in other study design and regression models. In the present study, we fitted DLMNs to fully understand the dose–response function and lag effects of temperature. Furthermore, we applied DLNMs to the case-only approach. This allows sophisticated non-linear and delayed effect modification to be estimated and tested statistically.
We found significant impacts of hot and cold temperatures on mortality in Guangzhou, China. Hot temperatures had an acute but short-term effect, whereas the effect of cold temperatures lasted 10–12 days. Similar lag structures with characteristics of short term for hot effects and long duration for cold effects were observed in England , US [17, 27] and Canada . This finding suggests that a longer timeframe are required to capture the cold impact, and that it may be inappropriate to specify a priori identical timeframe for cold and hot exposure. Guo et al.  stated that use of short lags may underestimate cold effects, but it may overestimate hot effects. In the present study, we estimated the effects of hot and cold temperatures for two different lag periods.
Some evidence in the literature shows that the magnitude of temperature effects varied by climate and population. Notably, we found a significant cold effect with a 20.39% increase in mortality risk comparing 8.2°C to 13.8°C. An analysis in the subtropical region of Chiang Mai, Thailand reported a similar cold effect with a 19% increase of mortality risk over lag 0–13 days comparing 19.35°C to 24.7°C . Multi-city studies reported a higher risk of mortality associated with cold exposure in regions having milder winter climates [4, 34]. Thus, the public perception that hot regions do not suffer from cold weather is completely mistaken. Our findings highlight a need to strengthen the awareness of combating cold exposure in the public.
By analyzing in subgroups, many studies have found that the elderly were at higher risk of mortality associated with cold [4, 10] or hot temperatures [4, 10, 19]. It may be due to their poorer physiological adaptation to changes in ambient temperatures. Our findings confirmed vulnerability of the elderly. An interesting finding is that females were more susceptible to hot but not cold temperatures compared to males, which is in agreement with the results in US , Europe [18, 36] and Korea . Gender difference may be dependent on location and population. For example, the impact of hot temperature on mortality was higher for women in Mexico, but higher for men in Sao Paulo .
Greater effects of ambient temperatures were observed for those with less education in the United States , Latin American , Korea [20, 37] as well as in Shanghai, China . However, no effect of education level on temperature-related mortality risk was found in California  and Australia . Moreover, all previous findings were based on stratified analysis and it remains uncertain whether effect modification by education was statistically significant or not. In the present study, educational attainment was classified into three groups. Subjects with no education were significantly more susceptible to the effects of both hot and cold temperatures compared to other education groups. Additionally, to our knowledge, this is the first study to examine the potential effect modification by occupation class. We found that blue-collar workers were at significantly higher risk of temperature-related mortality than white-collar workers.
When analyzing by cause of death, effect estimates were markedly higher for cardiovascular and respiratory deaths compared to all non-accidental deaths, consistent with previous studies [5, 13, 17]. From a public health point of views, this finding is important since cardio-respiratory diseases are the leading cause of death in Guangzhou, accounting for 58% of all registered deaths during the study period. This increase in mortality is likely related to the failure in thermoregulation and the physiological changes in circulatory system .
Some limitations should be mentioned. Firstly, the data are only from one city, so it should be cautious to generalize the findings to other geographic areas and other climates. Secondly, we considered several major mortality categories but not very fine categories. Notably, as in all time-series studies on temperature effects, we assigned each individual the same exposure level measured by ambient daily mean temperature, which would bring about measurement errors because indoor temperature may be not closely correlated with outdoor temperature due to the use of air condition. Lastly, temperature-associated mortality was calculated by comparing the first to the 10th percentile and the 99th to the 90th percentile temperatures. This accounted for the effects of extreme temperatures in some way. However, further research aimed at the definition of warning systems and prevention programs should assess the potential added effect of heat waves and cold spells especially in terms of persistency of extreme conditions, frequency of events and adaptation of the population throughout the summer season.
Our results confirmed that cold and hot temperatures were associated with increased risk of mortality in the subtropical city of Guangzhou, China. The elderly, females and subjects with low socioeconomic status have been identified as especially susceptible to temperature-associated mortality. The findings can contribute to focus community and individual prevention programs targeting on mitigating weather-related mortality.
Distributed lag non-linear models
- ICD 10:
International classification of diseases tenth revision
- NO2 :
- PM10 :
Particulate matter with aerodynamic diameters less than 10 μm
Quasi-likelihood Akaike information criteria
- SO2 :
Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, Analitis A, Anderson HR, Bisanti L, D'Ippoliti D, Danova J, Forsberg B, Medina S, Paldy A, Rabczenko D, Schindler C, Michelozzi P: Heat effects on mortality in 15 European cities. Epidemiology. 2008, 19: 711-719. 10.1097/EDE.0b013e318176bfcd.
Basu R, Samet JM: Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev. 2002, 24: 190-202. 10.1093/epirev/mxf007.
Braga AL, Zanobetti A, Schwartz J: The effect of weather on respiratory and cardiovascular deaths in 12 U.S. cities. Environ Health Perspect. 2002, 110: 859-863.
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-87. 10.1093/aje/155.1.80.
Hajat S, Armstrong BG, Gouveia N, Wilkinson P: Mortality Displacement of Heat-Related Deaths. Epidemiology. 2005, 16: 613-620. 10.1097/01.ede.0000164559.41092.2a.
Medina-Ramon M, Schwartz J: Temperature, temperature extremes, and mortality: a study of acclimatisation and effect modification in 50 US cities. Occup Environ Med. 2007, 64: 827-833. 10.1136/oem.2007.033175.
Bell ML, O'Neill MS, Ranjit N, Borja-Aburto VH, Cifuentes LA, Gouveia NC: Vulnerability to heat-related mortality in Latin America: a case-crossover study in Sao Paulo, Brazil, Santiago, Chile and Mexico City, Mexico. Int J Epidemiol. 2008, 37: 796-804.
Basu R: High ambient temperature and mortality: a review of epidemiologic studies from 2001 to 2008. Environ Health. 2009, 8: 40-10.1186/1476-069X-8-40.
Yu W, Mengersen K, Wang X, Ye X, Guo Y, Pan X, Tong S: Daily average temperature and mortality among the elderly: a meta-analysis and systematic review of epidemiological evidence. Int J Biometeorol. 2012, 56: 569-581. 10.1007/s00484-011-0497-3.
Almeida SP, Casimiro E, Calheiros J: Effects of apparent temperature on daily mortality in Lisbon and Oporto, Portugal. Environ Health. 2010, 9: 12-17. 10.1186/1476-069X-9-12.
Analitis A, Katsouyanni K, Biggeri A, Baccini M, Forsberg B, Bisanti L, Kirchmayer U, Ballester F, Cadum E, Goodman PG, Hojs A, Sunyer J, Tiittanen P, Michelozzi P: Effects of cold weather on mortality: results from 15 European cities within the PHEWE project. Am J Epidemiol. 2008, 168: 1397-1408. 10.1093/aje/kwn266.
Basu R, Dominici F, Samet JM: Temperature and Mortality Among the Elderly in the United States. Epidemiology. 2005, 16: 58-66. 10.1097/01.ede.0000147117.88386.fe.
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 Perspect. 2011, 119: 1719-1725. 10.1289/ehp.1103598.
Guo Y, Punnasiri K, Tong S: Effects of temperature on mortality in Chiang Mai city. Thailand: a time series study. Environ Health. 2012, 11: 36-10.1186/1476-069X-11-36.
McMichael AJ, Wilkinson P, Kovats RS, Pattenden S, Hajat S, Armstrong B, Vajanapoom N, Niciu EM, Mahomed H, Kingkeow C, Kosnik M, O’Neill MS, Romieu I, Ramirez-Aguilar M, Barreto ML, Gouveia N, Nikiforov B: International study of temperature, heat and urban mortality: the 'ISOTHURM' project. Int J Epidemiol. 2008, 37: 1121-1131. 10.1093/ije/dyn086.
Muggeo VM, Hajat S: Modelling the non-linear multiple-lag effects of ambient temperature on mortality in Santiago and Palermo: a constrained segmented distributed lag approach. Occup Environ Med. 2009, 66: 584-591. 10.1136/oem.2007.038653.
Anderson BG, Bell ML: Weather-Related Mortality: How Heat, Cold, and Heat Waves Affect Mortality in the United States. Epidemiology. 2009, 20: 205-213. 10.1097/EDE.0b013e318190ee08.
Stafoggia M, Forastiere F, Agostini D, Biggeri A, Bisanti L, Cadum E, Caranci N, De’Donato F, De Lisio S, De Maria M, Michelozzi P, Miglio R, Pandolfi P, Picciotto S, Rognoni M, Russo A, Scarnato C, Perucci CA: Vulnerability to Heat-Related Mortality: A Multicity, Population-Based, Case-Crossover Analysis. Epidemiology. 2006, 17: 315-323. 10.1097/01.ede.0000208477.36665.34.
O'Neill MS, Zanobetti A, Schwartz J: Modifiers of the temperature and mortality association in seven US cities. Am J Epidemiol. 2003, 157: 1074-1082. 10.1093/aje/kwg096.
Son JY, Lee JT, Anderson GB, Bell ML: The Impact of Heat Waves on Mortality in 7 Major Cities in Korea. Environ Health Perspect. 2012, 120: 566-571. 10.1289/ehp.1103759.
Gouveia N: Socioeconomic differentials in the temperature-mortality relationship in Sao Paulo, Brazil. Int J Epidemiol. 2003, 32: 390-397. 10.1093/ije/dyg077.
Ishigami A, Hajat S, Kovats RS, Bisanti L, Rognoni M, Russo A, Paldy A: An ecological time-series study of heat-related mortality in three European cities. Environ Health. 2008, 7: 5-10.1186/1476-069X-7-5.
Basu R, Ostro BD: A Multicounty Analysis Identifying the Populations Vulnerable to Mortality Associated with High Ambient Temperature in California. Am J Epidemiol. 2008, 168: 632-637. 10.1093/aje/kwn170.
Yu W, Vaneckova P, Mengersen K, Pan X, Tong S: Is the association between temperature and mortality modified by age, gender and socio-economic status?. Sci Total Environ. 2010, 408: 3513-3518. 10.1016/j.scitotenv.2010.04.058.
China Metrorological Adminstration: China Meteorological Data Sharing Service System. Available at: http://cdc.cma.gov.cn/(accessed on 12 July 2012)
Wong CM, Vichit-Vadakan N, Kan H, Qian Z: Public Health and Air Pollution in Asia (PAPA): a multicity study of short-term effects of air pollution on mortality. Environ Health Perspect. 2008, 116: 1195-1202. 10.1289/ehp.11257.
Gasparrini A, Armstrong B, Kenward MG: Distributed lag non-linear models. Stat Med. 2010, 29: 2224-2234. 10.1002/sim.3940.
Armstrong B: Models for the relationship between ambient temperature and daily mortality. Epidemiology. 2006, 17: 624-631. 10.1097/01.ede.0000239732.50999.8f.
Martin SL, Cakmak S, Hebbern CA, Avramescu M-L, Tremblay N: Climate change and future temperature-related mortality in 15 Canadian cities. Int J Biometeorol. 2012, 56: 605-619. 10.1007/s00484-011-0449-y.
Peng RD, Dominici F, Louis TA: Model choice in time series studies of air pollution and mortality. J R Stat Soc [Ser A]. 2006, 169: 179-203. 10.1111/j.1467-985X.2006.00410.x.
Armstrong BG: Fixed factors that modify the effects of time-varying factors: applying the case-only approach. Epidemiology. 2003, 14: 467-472.
Schwartz J: Who is sensitive to extremes of temperature? A case-only analysis. Epidemiology. 2005, 16: 67-72. 10.1097/01.ede.0000147114.25957.71.
Medina-Ramón M, Zanobetti A, Cavanagh DP, Schwartz J: Extreme temperatures and mortality: assessing effect modification by personal characteristics and specific cause of death in a multi-city case-only analysis. Environ Health Perspect. 2006, 114: 1331-1336. 10.1289/ehp.9074.
Group TE: Cold exposure and winter mortality from ischaemic heart disease, cerebrovascular disease, respiratory disease, and all causes in warm and cold regions of Europe. Lancet. 1997, 349: 1341-1346.
R Development Core Team R: A Language and Environment for Statistical Computing. Version 2.13.1.http://www.r-project
Hajat S, Kovats RS, Lachowycz K: Heat-related and cold-related deaths in England and Wales: who is at risk?. Occup Environ Med. 2006, 64: 93-100. 10.1136/oem.2006.029017.
Son JY, Lee JT, Anderson GB, Bell ML: vulnerability to temperature-related mortality in Seoul. Korea: Environ Res, in press
Ma W, Yang C, Tan J, Song W, Chen B, Kan H: Modifiers of the temperature–mortality association in Shanghai, China. Int J Biometeorol. 2011, 56: 205-207.
Nixdorf-Miller A, Hunsaker DM, Hunsaker JC: Hypothermia and hyperthermia medicolegal investigation of morbidity and mortality from exposure to environmental temperature extremes. Arch Pathol Lab Med. 2006, 130: 1297-1304.
This study was supported by National Nature Science Foundation of China  and Guangdong Nature Science Foundation [S2011040005355]. We thank Yuming Guo and Francesca de’Donato for helpful comments on earlier version of the manuscript.
The authors declare that they have no competing interests.
CQ and PY initiated the study and collected the data. JY, YD and YX cleaned the data and performed statistical analysis. JY and CQ drafted the manuscript. All authors read and approved the final manuscript.