Research Volume 121 | 2013
Short-Term Effects of the 2008 Cold Spell on Mortality in Three Subtropical Cities in Guangdong Province, China
Huiyan Xie,1 Zhibin Yao,2 Yonghui Zhang,3 Yanjun Xu,3 Xiaojun Xu,3 Tao Liu,1,3 Hualiang Lin,1,3 Xiangqian Lao,4 Shannon Rutherford,5 Cordia Chu,5 Cunrui Huang,1,6 Scott Baum,5 and Wenjun Ma1,3
Background: Few studies have been conducted to investigate the impact of extreme cold events on mortality in subtropical regions.
Objective: In the present study we aimed to investigate the effects of the 2008 cold spell on mortality and the possibility of mortality displacement in three subtropical cities in China.
Methods: Daily mortality, air pollution, and weather data were collected from 2006 to 2009 in Guangzhou, Nanxiong (no air pollutants), and Taishan. We used a polynomial distributed lag model (DLM) to analyze the relationship between the 2008 cold spell and mortality. To observe the mortality displacement of the cold spell, we estimated the cumulative effects at lag0, lag0–6, lag0–13, lag0–20, and lag0–27 separately.
Results: During the 2008 cold spell, the cumulative risk of nonaccidental mortality increased significantly in Guangzhou [relative risk (RR) = 1.60; 95% CI: 1.19, 2.14] and Taishan (RR = 1.60; 95% CI: 1.06, 2.40) when lagged up to 4 weeks after the cold spell ended. Estimated effects at lag0–27 were more pronounced for males than for females, for respiratory mortality than for cardiovascular mortality, and for the elderly (≥ 75 years of age) than for those 0–64 years of age. Most of the cumulative RRs increased with longer lag times in Guangzhou and Taishan. However, in Nanxiong, the trend with cumulative RRs was less consistent, and we observed no statistically significant associations at lag0–27.
Conclusion: We found associations between the 2008 cold spell and increased mortality in the three subtropical cities of China. The lag effect structure of the cold spell varied with location and the type of mortality, and evidence of short-term mortality displacement was inconsistent. These findings suggest that extreme cold is an important public health problem in subtropical regions.
Citation: Xie H, Yao Z, Zhang Y, Xu Y, Xu X, Liu T, Lin H, Lao X, Rutherford S, Chu C, Huang C, Baum S, Ma W. 2013. Short-Term Effects of the 2008 Cold Spell on Mortality in Three Subtropical Cities in Guangdong Province, China. Environ Health Perspect 121:210–216; http://dx.doi.org/10.1289/ehp.1104541
Address correspondence to W. Ma, Guangdong Provincial Institute of Public Health, Center for Disease Control and Prevention of Guangdong Province, 176, Xin Gang Xi Rd., Guangzhou, China. Telephone: (8620) 84458530. E-mail: firstname.lastname@example.org
This work was supported in part by the Adapting to Climate Change in China project, funded by the UK Department for International Development, the Swiss Agency for Development and Cooperation, and the UK Department for Energy and Climate Change.
The authors declare they have no actual or potential competing financial interests.
Received: 27 September 2012
Accepted: 31 October 2012
Advance Publication: 31 October 2012
Final Publication: 1 February 2013
Climate change is likely to cause increased occurrence of extreme weather events, including both heat waves and cold spells (Molloy et al. 2008). Many studies have examined the relationship between extreme temperature events and mortality (Gómez-Acebo et al. 2010; Hajat et al. 2005; Iniguez et al. 2010; Kaiser et al. 2007; Rooney et al. 1998; Sartor et al. 1995), mainly focusing on heat waves to demonstrate the effects of global warming (Gasparrini and Armstrong 2011; Knowlton et al. 2009; Le Tertre et al. 2006; Semenza et al. 1996; Tong et al. 2010), but fewer studies have examined the health effects of extreme cold spells (Kysely et al. 2009; Montero et al. 2010). As noted by Ballester et al. (2003), some studies have reported greater cold-related mortality than heat-related mortality; in addition, heat wave effects appear to last for a few days at most, whereas effects of cold spells may persist for up to 2 months. Most studies on the impact of extreme cold events have been conducted in temperate cities in developed countries (Analitis et al. 2008; Cagle and Hubbard 2005; Healy 2003; O’Neill et al. 2003). Estimated effects of temperature on mortality may be heterogeneous across areas with differing socioeconomic status and education level (Basu and Samet 2002; Bell et al. 2008). However, few studies have been conducted in tropical or subtropical cities in developing countries.
Guangdong, a subtropical province in China, experienced an unusually persistent and widespread severe cold spell in 2008. This event also affected 20 other provinces across southern China. The daily mean temperature during this extreme weather event was much lower than that for the same period in previous years. Although intensive public attention was focused on the adverse impact of this cold spell on ecological, social, and economic systems, health impacts on local residents have not been studied (Liangxun et al. 2009).
Many previous studies on associations of temperature with mortality have considered delayed effects (Bell et al. 2008; Hajat et al. 2005; Hertel et al. 2009; Huynen et al. 2001; Kysely 2004), including lagged effects of temperature on single days, and of moving average temperature on subsequent days. For example, Bell et al. (2008) estimated the association between high temperature and mortality using single day lags of 0, 1, 2, and 3 days and cumulative lags up to 1 week (lags 0–6) using a moving average. This approach could overestimate the effects of current-day exposure by ignoring effects of exposure on previous days (Gasparrini et al. 2010). However, it may also underestimate effects of exposure on mortality if effects persist longer than the observed lag period (Roberts and Martin 2007; Schwartz 2000). Distributed lag models (DLMs), which allow a detailed representation of the time course of the exposure–response relationship while avoiding problems related to colinearity among lagged exposure variables, have been proposed for analyses of delayed effects (Schwartz 2000). Numerous studies have applied DLMs to analyze lagged health effects of temperature, primarily for continuous-temperature time-series analysis (Analitis et al. 2008; Ha et al. 2011; Hajat et al. 2005; Liu et al. 2011).
In the present study we aimed to assess the health impacts of the 2008 cold spell in three subtropical cities of Guangdong by analyzing extended time-series data for daily mortality and modeling lagged effects using distributed lag models. The findings of this study will improve our understanding of relationships between extreme cold events and mortality in subtropical areas and provide evidence to support the need to develop adaptation strategies to mitigate the adverse effects of cold climate extremes in the context of climate change.
Materials and Methods
Study settings. Guangdong is one of China’s southernmost provinces. It has a typical subtropical climate with an average annual temperature of 22°C. Data were collected for three cities located in different parts of the province (Figure 1): Nanxiong, the northernmost city, with a population of > 400,000 by the end of 2009; Guangzhou, the centrally located capital of Guangdong Province, with a total population of > 7 million; and Taishan, a coastal city in southern Guangdong, with a population > 900,000 by the end of 2009. On the basis of data availability, we used data from two districts of Guangzhou (Yue Xiu and Li Wan, with an estimated population of 1.86 million in 2009) for this study.
Data sources. Daily nonaccidental mortality data from 1 January 2006 through 31 December 2009 were obtained from the Guangdong Provincial Center for Disease Control and Prevention. The original data source was death certificates, which included the age and sex of the deceased and the date and cause of death. Nonaccidental causes of deaths were categorized using codes A00–R99 from the International Classification of Diseases, 10th Revision (ICD–10; World Health Organization 2007). ICD-10 codes J00–J99 represent respiratory diseases, and codes I00–I99 represent cardiovascular diseases.
Daily meteorological data were collected from the local meteorological bureaus of each city from 1 January 2006 through 31 December 2009. We used maximum temperature, mean temperature, minimum temperature (Tmin), and relative humidity (RH) for the analysis.
Air pollution data obtained for the same period included daily average concentrations of particulate matter with aerodynamic diameter < 10 μm (PM10), nitrogen dioxide (NO2), and sulfur dioxide (SO2), all given in milligrams per cubic meter. Air pollution data were measured at single environmental monitoring sites located in the centers of Guangzhou (monitored continuously) and Taishan (monitored hourly). Air pollution data were not available for Nanxiong.
Definition of cold spell. A variety of approaches have been used to define a cold spell (Kysely et al. 2009; Lin et al. 2011; Montero et al. 2010), but there is no universally accepted definition based on specific temperatures. The Chinese National Bureau of Meteorology defines a cold spell as a period with a temperature decrease of at least 8°C over 48 hr that results in a Tmin < 4°C (China Meteorological Administration 2008). However, this definition was inappropriate for Guangdong because the province is a subtropical region. Because the Tmin of the three cities correlated more closely with their respective mortality than did maximum temperature and mean temperature (results not shown), the definition of cold spell that we adopted for this analysis was based on the daily Tmin. Therefore, in this study, we defined a weather fluctuation as a cold spell if the minimum daily temperature fell below the 5th percentile of temperatures recorded at that location from January 2006 through December 2009 for at least 5 consecutive days. This definition was very similar to that used in a previous meteorological study in China (Wu et al. 2008). According to this definition, the 2008 cold spell lasted between 18 and 21 days in the three sampled Guangdong cities. Table 1 shows detailed information on this cold spell in these three cities.
Calculation of excess mortality. To estimate excess mortality attributable to the 2008 cold spell, we calculated 31-day moving averages of daily mortality during the cold spell and during the same days in the 2 years before the cold spell and the year after the spell combined (Rooney et al. 1998). Excess mortality was assessed as the difference between the number of deaths observed on a given day during the 2008 cold spell and the corresponding moving average values for 2006, 2007, and 2009 combined. We calculated an approximate conﬁdence interval (CI) for the excess mortality by treating the total number of deaths during the cold spell as a Poisson distribution and comparing the upper and lower 95% conﬁdence bounds of this value with the expected number of deaths.
Estimation of city-specific relative risk. We evaluated the association between the 2008 cold spell and daily mortality using Poisson regression with a distributed lag model. For Poisson regression, the unconstrained distributed lag model may be written as
Log(μt) = α + COVs + β0Zt + β1Zt – 1 + … + βj Zt – q, 
where COVs represents all other covariates in the model, and Zt represents cold-spell exposure delayed over time for j = 0 … q days. In this study, we defined Zt as a binary variable that equals 1 for the 2008 cold spell days and 0 for other days.
To gain more precision in the estimate of the distributed lag curve, a polynomial distributed lag constrains the βj to follow a polynomial pattern in the lag number:
where j is the number of lag days and d is the degree of the polynomial. We chose to use a third-degree polynomial in this study to ensure enough degrees of freedom (df) to fit the pattern of response over time. We specified the lagged effect of the cold spell up to 27 days, consistent with previous studies (Armstrong 2006; Guo et al. 2011).
We estimated the cumulative mortality risk associated with the cold spell using the cross-basis functions for the spaces of the cold spell and the lag dimension [based on a bidimensional functional space expressed by the combination of two sets of basis functions, which specify the relationships in the dimensions of the cold spell and lags (Gasparrini 2011)] as a covariate in the Poisson regression model. Cumulative mortality risk and 95% CIs were estimated by comparing mortality during the cold spell with mortality during the non–cold spell periods. To observe mortality displacement, we estimated cumulative effects for lag 0, lag 0–6, lag 0–13, lag 0–20, and lag 0–27 days.
Relative humidity, PM10, NO2, and SO2 were modeled as natural cubic splines with 3 df in models for Guangzhou and Taishan, as described previously (Anderson and Bell 2009; Guo et al. 2011). However, we did not adjust for air pollutants in Nanxiong because air pollutant data were not available. We also modeled a binary variable assigned as 1 on days when any influenza deaths were reported and 0 otherwise (ICD-10 codes J10–J11) to account for influenza viral activity, similar to the approach used by Braga et al. (2000). To control for seasonality and long-term trends, we included a smooth function of time and also included day of the week as a covariate in the models. Therefore, the complete Poisson regression model was
Yt~Poisson (μt) Log(μt) = α + β0Zt + β1Zt – 1 + … + βjZt – q + S(RHt, 3) + S(PM10t, 3) + S(SO2t, 3) + S(NO2t, 3) + S(Time,8/year) + ηDOWt + υInfluenzat = α + β0Zt + β1Zt – 1 +…+ βjZt – q + COVs. 
Here, t is the date of the observation; Yt is the observed daily death count on day t; α is the intercept; Zt is the cold spell exposure on the same day (lag0), with Zt – 1 on the previous day (lag1), and so on. S() is a natural cubic spline. RHt, PM10t, SO2t, and NO2t represent the relative humidity and concentrations of PM10, NO2, and SO2, respectively, at time t with 3 df; S(time,8/year) is the natural cubic spline of time with 8 df per year, which was chosen by minimizing the Akaike information criterion (AIC; Akaike 1974). DOWt is the day of the week on day t, and η is the vector of coefficients. Influenza is a binary variable that is 1 if there are any influenza deaths on day t and 0 if there are not.
All statistical tests were two-sided, and values of p < 0.05 were considered statistically significant. We used R software (version 2.11.0; R Development Core Team, http://www.R-project.org/) and SAS software (version 9.1; SAS Institute Inc., Cary, NC, USA; Daly 1992) to analyze the data. The dlnm package in R software was used to construct the polynomial distributed lag basis.
Table 2 summarizes the weather, air pollutant, and mortality statistics during the 2008 cold spell and the corresponding days during 2006, 2007, and 2009 in Guangzhou, Nanxiong, and Taishan. Compared with the same periods in 2006, 2007, and 2009, the mean daily Tmin during the 2008 cold spell was > 7°C lower, with the lowest mean daily Tmin reaching 2.26°C in Nanxiong. Mean daily death counts during 2006, 2007, and 2009 of 37, 7, and 23 in Guangzhou, Nanxiong, and Taishan, respectively, increased to 55, 11, and 32 during the 2008 cold spell.
Figure 2 shows the increase of daily death counts observed in the three cities during the 2008 cold spell relative to the means for corresponding days in 2006, 2007, and 2009. The largest increase in mortality was observed in Nanxiong, with 52% more deaths than the average for the corresponding days in 2006, 2007 and 2009, and the smallest increase in deaths was observed in Taishan, with 35% more deaths than in 2006, 2007, and 2009 (Table 3). The excess mortality rate increased dramatically with age in all three cities, and was highest for residents > 75 years of age in Nanxiong (427.2 excess deaths per 100,000; 95% CI: 336.6, 543.7).
Table 4 presents the cumulative relative risks (RRs) for lag0–27 days by age group, sex, and cause of death. There was a significant increase in nonaccidental mortality during the 2008 cold spell for all ages combined in Guangzhou and Taishan, both before and after adjustment for air pollution, but the increase in mortality in Nanxiong was not statistically significant for any age group or according to sex or cause of death (p > 0.05). In Guangzhou and Taishan, the estimated mortality associated with the 2008 cold spell was higher for males than for females (RR = 1.56; 95% CI: 1.07, 2.28 vs. RR = 1.31; 95% CI: 0.85, 2.02 in Guangzhou; RR = 1.93; 95% CI: 1.14, 3.27 vs. RR = 1.27; 95% CI: 0.71, 2.27 in Taishan), for respiratory mortality than for cardiovascular mortality (RR = 2.15; 95% CI: 1.11, 4.16 vs. RR = 1.40; 95% CI: 0.87, 2.25 in Guangzhou; RR = 3.23; 95% CI: 1.38, 7.58 vs. RR = 1.67; 95% CI: 1.02, 2.73 in Taishan), and for those ≥ 75 years of age than for those 0–64 years of age (RR = 1.48; 95% CI: 1.01, 2.17 vs. RR = 1.42; 95% CI: 0.77, 2.60 in Guangzhou; RR = 2.19; 95% CI: 1.34, 2.60 vs. RR = 1.10; 95% CI: 0.50, 2.44 in Taishan).
To evaluate the lag structure of effects of the cold spell on mortality, including potential effects of mortality displacement, we estimated cumulative effects by age group, sex, and cause of death for different lags using the distributed lag model (Figure 3). The cumulative RRs based on these analyses can be interpreted as the net effects of the cold spell after accounting for mortality displacement, which is characterized by an increasing trend of cumulative RRs for exposures at lower lags (resulting in part from deaths that occurred earlier in time as a consequence of exposure) followed by a decline in cumulative RRs at higher lags (because of the relative deficit in deaths that have been displaced forward in time) (Hajat et al. 2005; Roberts and Switzer 2004). In general, RRs were lowest at lag0; in Guangzhou and Taishan, cumulative RRs increased with longer cumulative lags, with the highest RR at lag27. However, in Nanxiong, the highest cumulative RRs (except for those affecting residents < 75 years of age or females) were observed at lag0–13, after which they decreased slowly, suggesting a deficit offset for only part of the overall excess after 2 weeks of exposure.
In view of the global change in climate predicted for future decades, the frequency, intensity, and duration of extreme climate events are expected to change (Albritton et al. 2001). Understanding the relationship between extreme climate events, such as heat waves and cold spells, and their potential health impacts is the first step in managing and reducing the adverse impact of such events. To our knowledge, the present study is unique in estimating the short-term effects of a cold spell on mortality in multiple subtropical cities in China using distributed lag models.
Our estimates of increased mortality during the 2008 cold spell are much higher than those from studies in the Netherlands (Huynen et al. 2001), Russia (Revich and Shaposhnikov 2008), or the Czech Republic (Kysely et al. 2009). There are several possible reasons for this. First, techniques used to measure increased mortality varied across studies. For example, different reference baselines can lead to different estimated values for increased mortality. For example, Kysely et al. (2009) calculated the expected (baseline) number of deaths using the mean annual cycle smoothed by 15-day running means adjusting for the observed year-to-year changes in mortality. Second, populations who live in cold climates may be more accustomed to and prepared for extreme cold weather than subtropical residents. The 2008 cold spell that occurred in Guangdong was popularly considered to be the most extreme cold spell in five decades, and hence populations of these three subtropical cities may be more sensitive to extreme cold weather (Analitis et al. 2008). Third, in subtropical areas such as Guangdong Province, few buildings have heating systems equipped to provide enough heat in extremely cold weather; this may contribute to a greater risk for vulnerable populations, such as the elderly. Finally, the health care system, especially the emergency service, could not meet the sudden increase in the need for care during the extreme cold event in southern China, when the number of ambulance calls rose to such a high number that nearly one-fourth of the calls received no response (Xiong et al. 2010).
Our findings suggest that the elderly suffered the most during the 2008 cold spell, with those > 75 years of age being the most susceptible, consistent with some previous studies (Analitis et al. 2008; Huynen et al. 2001; Iniguez et al. 2010; Revich and Shaposhnikov 2008). However, our findings differed from a study conducted in the Czech Republic, in which cold spells had the greatest effect on middle-aged men who died from cardiovascular disease (Kysely et al. 2009). The authors attributed this finding to occupational exposure to the cold in men 25–59 years of age, whereas the elderly tended to stay indoors during the cold spell and thus avoided direct exposure to ambient temperatures.
We also found that the effects of the cold spell appeared to be more pronounced for respiratory disease patients than that for cardiovascular patients, consistent with previous findings from Russia (Revich and Shaposhnikov 2008) and Europe (Analitis et al. 2008). Some researchers have attributed the increased deaths from respiratory diseases to increased infection from indoor crowding, adverse effects of cold weather on the immune system, and the fact that low temperatures may facilitate the survival of bacteria in droplets (The Eurowinter Group 1997; Handley and Webster 1995).
In the present study, the effects of the 2008 cold spell on nonaccidental mortality appeared to be greater within 2 weeks of exposure in Nanxiong compared with the other two cities. One possible explanation for this is that Nanxiong is the northernmost city of the three cities and had the lowest average Tmin during the 2008 cold spell. Another possible reason is a lower adaptive capacity to extreme weather events in Nanxiong, as reflected by the lower socioeconomic status of the city, for which gross national product per capita was much lower than in the other two cities. To explore this further, it would be necessary to study the social determinants of adaptive capacity for extreme weather events.
We also found an apparent rise in respiratory disease mortality that lasted up to 4 weeks after the cold spell. This was consistent with two studies in which the apparent effects of low temperatures continued for a longer period of time than did estimated effects of heat waves (Healy 2003; Iniguez et al. 2010). In Nanxiong and Guangzhou, estimated effects decreased after 2 weeks of the cold spell for persons > 75 years of age, suggesting some compensatory risk reduction consistent with a harvesting phenomenon (Hajat et al. 2005). However, this phenomenon was not observed in Taishan. This discrepancy needs to be explored further.
We controlled for air pollution effects in Guangzhou and Taishan (pollution data were not available for Nanxiong), and the estimated effects of the cold spell showed a slight reduction after adjustment. This might be explained by the relative impacts of air pollution on mortality during the cold spell, because air pollutant concentrations are likely to increase during such episodes (O’Neill et al. 2005; Schwartz 2000; Tong et al. 2010). Previous studies have suggested that the effects of air pollution on mortality are much lower than the effects of temperature (Ren et al. 2011; Spickett et al. 2011). Thus, the relationship we observed between mortality and the cold spell was not likely to have been substantially confounded by the effects of air pollution. However, we did not consider the possible effects of indoor air pollution from smoking, cooking, and home-heating fuels on mortality.
In the present study, we used a binary indicator to describe both cold-spell and non–cold-spell days to attain unique DLM coefficients representing overall mortality effects of the cold spell period. This approach is different from modeling mortality as a continuous function of temperature, as has been done in previous studies (Hertel et al. 2009; Kaiser et al. 2007). Although we believe these approaches are comparable from a purely conceptual viewpoint, the validity of our approach should be evaluated in future studies.
Our findings suggest that further research is needed. First, studies should be conducted based on longer time-series data with multiple cold spells to estimate the impact of cold spells according to their duration or intensity. Second, the estimated effects of weather on mortality may have been influenced by the age structure, sociodemographic characteristics, and environmental conditions of each population. Further research on factors that determine vulnerability to cold would help inform the development and implementation of cold-weather emergency plans. Information on the effects of indoor environments, energy usage, and human thermal comfort thresholds on vulnerability would also help determine appropriate strategies for adapting to a changing climate.
The 2008 cold spell was associated with an increase in daily mortality in three subtropical cities of Guangdong Province, China. As a subtropical region, Guangdong is relatively ill equipped to adapt to extreme cold events. For example, most temperature control systems in buildings in Guangdong Province were designed for cooling, not heating. Climate models indicate that seasonal weather patterns and conditions will continue to vary from current climate conditions as average global temperatures increase (Albritton et al. 2001), and climate change is expected to contribute to an increase in the intensity of extreme cold events as well as heat waves (Lionello et al. 2008). It is both necessary and timely for governments and relevant sectors to develop adaptive plans for such extreme events. Similar to the heat-watch warning system adopted in the United States (Kalkstein 2000), subtropical cities need to develop cold weather emergency plans to improve the delivery of health emergency services, and also to issue timely weather alerts when extreme events are expected. On the basis of findings in this study, decision makers from subtropical regions not only should pay attention to heat waves but also must consider adaptive measures to protect vulnerable populations from extreme cold events.
Albritton DL, Allen MR, Baede APM, Church JA, Cubasch U, Xiaosu D, et al. 2001. Summary for policymakers: a report of Working Group I of the Intergovernmental Panel on Climate Change. In: Climate Change 2001: The Scientific Basis: Summary for Policymakers (Houghton JT, Ding Y, Griggs DJ, Noguer M, van der Linden PJ, Dai X, et al., eds). Cambridge, UK:Cambridge University Press, 1–20. Available: http://www.acrim.com/%5C/Reference%20Files/CLIMATECHANGE%202001%20-%20The%20Scientific%20Basis.pdf [accessed 21 December 2012].
Analitis A, Katsouyanni K, Biggeri A, Baccini M, Forsberg B, Bisanti L, et al. 2008. Effects of cold weather on mortality: results from 15 European cities within the PHEWE project. Am J Epidemiol 168:1397–1408.
Bell ML, O’Neill MS, Ranjit N, Borja-Aburto VH, Cifuentes LA, Gouveia NC. 2008. Vulnerability to heat-related mortality in Latin America: a case-crossover study in São Paulo, Brazil, Santiago, Chile and Mexico City, Mexico. Int J Epidemiol 37:796–804.
China Meteorological Administration. 2008. Cold Spell Prevention. Available: http://www.cma.gov.cn/fzjz/qxzhyjxhyfyzn/200805/t20080508_6578.html [accessed 21 October 2010].
Guo Y, Barnett AG, Pan X, Yu W, Tong S. 2011. The impact of temperature on mortality in Tianjin, China: a case-crossover design with a distributed lag nonlinear model. Environ Health Perspect 119:1719–1725.
Kaiser R, Le Tertre A, Schwartz J, Gotway CA, Daley WR, Rubin CH. 2007. The effect of the 1995 heat wave in Chicago on all-cause and cause-specific mortality. Am J Public Health 97(suppl 1):S158–S162.
Knowlton K, Rotkin-Ellman M, King G, Margolis HG, Smith D, Solomon G, et al. 2009. The 2006 California heat wave: impacts on hospitalizations and emergency department visits. Environ Health Perspect 117:61–67.
Kysely J, Pokorna L, Kyncl J, Kriz B. 2009. Excess cardiovascular mortality associated with cold spells in the Czech Republic. BMC Public Health 9:19; doi:10.1186/1471-2458-9-19 [Online 15 January 2009].
Liu L, Breitner S, Pan X, Franck U, Leitte AM, Wiedensohler A, et al. 2011. Associations between air temperature and cardio-respiratory mortality in the urban area of Beijing, China: a time-series analysis. Environ Health 10:51; doi:10.1186/1476-069X-10-51 [Online 25 May 2011].
O’Neill MS, Hajat S, Zanobetti A, Ramirez-Aguilar M, Schwartz J. 2005. Impact of control for air pollution and respiratory epidemics on the estimated associations of temperature and daily mortality. Int J Biometeorol 50:121–129.
The Eurowinter Group 1997. Cold exposure and winter mortality from ischaemic heart disease, cerebrovascular disease, respiratory disease, and all causes in warm and cold regions of Europe. Lancet 349:1341–1346.
World Health Organization. 2007. International Statistical Classification of Diseases and Related Health Problems, 10th Revision, Version for 2007. Geneva:World Health Organization. Available: http://apps.who.int/classifications/icd10/browse/2010/en [accessed 21 December 2012].