Health impact of the 2008 cold spell on mortality in subtropical China: the climate and health impact national assessment study (CHINAs)
Environmental Health volume 13, Article number: 60 (2014)
Many studies have investigated heat wave related mortality, but less attention has been given to the health effects of cold spells in the context of global warming. The 2008 cold spell in China provided a unique opportunity to estimate the effects of the 2008 cold spell on mortality in subtropical regions, spatial heterogeneity of the effects, stratification effect and added effects caused by sustained cold days.
Thirty-six study communities were selected from 15 provinces in subtropical China. Daily mortality and meteorological data were collected for each community from 2006 to 2010. A distributed lag linear non-linear model (DLNM) with a lag structure of up to 27 days was used to analyze the association between the 2008 cold spell and mortality. Multivariate meta-analyses were used to combine the cold effects across each community.
The 2008 cold spell increased mortality by 43.8% (95% CI: 34.8% ~ 53.4%) compared to non-cold spell days with the highest effects in southern and central China. The effects were more pronounced for respiratory mortality (RESP) than for cardiovascular (CVD) or cerebrovascular mortality (CBD), for females more than for males, and for the elderly aged ≥75 years old more than for younger people. Overall, 148,279 excess deaths were attributable to the 2008 cold spell. The cold effect was mainly from extreme low temperatures rather than sustained cold days during this 2008 cold spell.
The 2008 cold spell increased mortality in subtropical China, which was mainly attributable to the low temperature rather than the sustained duration of the cold spell. The cold effects were spatially heterogeneous and modified by individual-specific characteristics such as gender and age.
The Intergovernmental Panel on Climate Change (IPCC) has projected that in the coming decades, extreme weather events will become more frequent and more intense in some parts of the world and such events will impact on health . Generally, the health effects of extreme heat events are acute and some harvesting is observed, but the effects of extreme cold temperatures are generally more prolonged than heat without mortality displacement [2–4]. Due to the projections associated with climate change many more studies have been conducted on health effects of heat waves compared to cold-related health impacts [2, 5–11]. Moreover, most previous studies on health effects of cold spells were conducted in temperate climate developed countries with very few in tropical or subtropical regions [2, 8, 9]. However, the health effects of extreme cold spells may be larger in these warm regions because populations are not acclimatized to cold spells and are unprepared for such events .
Traditional quantitative approaches to investigate the health effects of temperature are of 2 types-episode analysis and continuous-temperature time-series analysis [10, 11]. Recently, Gasparrini and his colleagues combined the two approaches to investigate heat wave-related mortality, which divided heat effects into two parts: independent contributions of daily temperature occurrences and the effects of exposure to hot temperatures protracted for several days . This approach is helpful to better understand the mechanism of temperature on mortality, and provides more significant information on public health intervention and cold-related burden estimation under projected climate change scenarios. However, few studies used this new approach to investigate the cold spell effects .
In January and February 2008, the majority of China, especially fifteen provinces of subtropical southern China experienced a severe continuous cold spell of a long duration, with lower than normal temperatures, heavy precipitation and thick snow deposition. Average air temperature during the cold spell days was 2–4°C lower than that during the same period of neighboring years. This event is considered a once in 50–100 years event. The estimated direct economic losses were more than US $22.3 billion [13, 14]. This extreme weather event provides a unique opportunity to assess the health impact of extreme cold spells on populations in subtropical regions (between the Tropic of Cancer and Tropic of Capricorn) . Three previous Chinese studies, two from Shanghai and one from Guangdong Province have reported that the 2008 cold spell significantly increased mortality or morbidity risk [3, 16, 17]. However, some health issues associated with this event remain unclear. Firstly, these published studies only focused on a single city or province, which prevents an estimate of the total effect of this extreme cold spell on mortality in subtropical China . Secondly, the spatial distribution of effects of this event on mortality across southern China are unknown, and such information is helpful to identify vulnerable regions and populations . Thirdly, whether the duration of extreme low temperatures sustained for several consecutive days imposed an added effect to the independent effects of daily temperature levels is unclear, and as this has been observed in relation to heat wave related effects . This type of analysis is helpful to better understand how it impacts on health. Thus, it is necessary to conduct multi-community studies to comprehensively assess health impacts of the 2008 cold spell in China, which can expand our understanding of the health impacts of cold spell events in subtropical regions. This study aimed to estimate the effects of the 2008 cold spell on mortality in subtropical regions of China, spatial heterogeneity of the effects, stratification effect and added effects caused by sustained cold days.
Thirty-six communities in subtropical China affected by the 2008 cold spell were selected for this study. These communities were all selected from the China’s Disease Surveillance Points system (DSPs) according to mortality data quality (mortality > 0.6%) and population size (>200,000) in order to assure enough daily death counts for community-specific time series analysis. Of the 36 communities, 12 communities were in urban areas, and the others in rural areas. They are distributed across 4 geographical regions: eastern China (Jiangsu Province, Zhejiang Province, Anhui Province and Shanghai Municipality), central China (Henan Province, Hubei Province, Hunan Province and Jiangxi Province), southern China (Fujian Province, Guangdong Province and Guangxi Province) and southwestern China (Sichuan Province, Guizhou Province, Yunnan Province and Chongqing Municipality). The 4 regions include 15 provinces/municipalities, home to 784 million inhabitants. The number of study communities in each region was 12, 8, 7 and 9, respectively (Figure 1).
Daily non-accidental mortality data from each study community were obtained from the Chinese Center for Disease Control and Prevention for the period December to March of 2006–2010. The original data was from death certificates, which included the causes and dates of death, gender, age, and place of death (hospital wards, emergency rooms and home). Non-accidental causes of deaths were categorized using codes A00–R99 from the International Classification of Diseases 10th Revision (ICD–10) . The codes J00–J99, I00-I99 and I60-I69 represent respiratory diseases (RESP), cardiovascular diseases (CVD), and cerebrovascular diseases (CBD), respectively. Demographic data in each community were obtained from the sixth national population census which was conducted in 2010 .
The corresponding daily meteorological data were collected from the China Meteorological Administration and included latitude, daily average temperature (Tm), maximum temperature (Tmax), minimum temperature (Tmin), average wind speed (WS) and relative humidity (RH).
Study duration and definition of cold spell
In order to exclude the influence of high temperature on mortality, only winter seasons (December to February) in 2006–2010 were included in this study. The 2008 cold spell started in the middle of January, and ended in the middle of February. In order to estimate the expected long lag effects [3, 21], March was also included in the study.
Definitions of cold spell vary greatly across different studies [4, 11, 22, 23]. For instance, Kysely et al. defined a cold spell as “a period of days on which air temperature did not exceed -3.5°C” , and Hickey suggested that “a cold spell could consist of a period of 10 consecutive days when the minimum air temperature was 5°C or more below normal” . In our previous study, we defined “a weather fluctuation as a cold spell if the minimum daily temperature fell below the 5th percentile of temperatures recoded at that community from January 2006 through December 2009 for at least 5 consecutive days” . However, some researchers have argued that it may be more appropriate to use daily average temperature to define a cold spell because it reflects the exposure throughout the whole day, while minimum or maximum temperature only represent a short period. Furthermore, daily average temperature can be easily interpreted for decision making purposes [16, 24]. Therefore, in the current study, a cold spell was defined as five or more consecutive days with daily average temperature falling below the 5th percentile of daily mean temperatures recorded at each community from December to March of 2006–2010. According to this definition, study days were divided into three groups: 2008 cold spell days, cold spell days in 2006, 2007, 2009 and 2010, and non-cold spell days.
Estimation of community-specific excess risk of the 2008 cold-spell
Where t is the day of observation; E(Y t ) is the expected number of deaths on day t; Z represents cold spell exposure, which was defined as a categorical variable (2 = 2008 cold spell days, 1 = cold spell days in 2006, 2007, 2009 and 2010, and 0 = non-cold spell days). “Non-cold spell days” were defined as the reference group for calculating relative risks (RR). T t,l is a matrix obtained by applying the DLNM to cold spell; β is the vector of coefficients for T t,l, and l is the number of lag days. We employed a B-spline function and a natural cubic spline function to estimate the non-linear and lagged effect of cold spell, respectively. To completely capture the overall effects of the 2008 cold spell exposure, a lag structure of up to 27 days was fitted, which is consistent with previous studies [14, 21, 27]. Degrees of freedom (df) for the lag structure were chosen based on Akaike information criterion (AIC) . It was found that 2 df s for non-linear and 5 df s for the lag produced the best model fitting. ns() is a natural spline. Both df s for RH and WS were set to 3, consistent with previous studies [13, 21, 29]. Another 3 df s were used to smooth year, calendar month and day to control for secular and seasonal trends. DOW is a dummy variable representing day of the week, and η is vector of coefficients. COVs represents all other covariates in the model. In this model, excess risks (ER) were reported, which indicated the percentage of increased death risk due to exposure to the 2008 cold spell compared to the non-cold spell. The cumulative effect of the 2008 cold spell on mortality during the lag 0–27 days was defined as cumulative excess risk (CER) .
Estimation of community-specific main and added effects of the 2008 cold-spell
As demonstrated by Gasparrini et al., the effect of heat waves can be divided into the main effect induced by the independent effects of daily temperature levels, and an added effect related to the duration of high temperatures sustained for several consecutive days . Huang et al. also employed this method to estimate the added effect of cold spell on CVD mortality in Brisbane . In this study, we also estimated the main and added effects of the 2008 cold spell on mortality. Another DLNM was employed, which is expressed as:
Zm and Za represent the variables for the main and added effects of the 2008 cold spell; T mt,l and T at,l are two matrixes obtained by applying the DLNMs to the main and added effects of the 2008 cold spell, while β m and β a are vectors of coefficients for T mt,l and T at,l , respectively. l is the number of lag days. We set both the maximum lag duration of the main and added effects of the 2008 cold spell as 27 days. We employed a B-spline function to estimate the non-linear main effects, and natural cubic spline functions to estimate the lagged main and added effects of the 2008 cold spell. According to the AIC principles , df of 5 was employed for the B-spline function of the main effect. Both df s of 5 were set for natural cubic spline functions of the main as well as added effects. The main effect was estimated by predicting the relative risk for median temperature among 2008 cold spell days versus the 50th percentile of temperature among non-cold spell days. This reference was chosen as a temperature at which little if any adverse effect of temperature on mortality is expected. Before estimating the added effect of the 2008 cold spell, the variable was coded as a dichotomous variable, 0 indicating non-cold spell days and the first day of the 2008 cold spell, and 1 indicating the remaining days of the 2008 cold spell. Therefore, the added effect was computed by estimating the relative risk of 1 compared to 0 . Other parameter settings were the same as in equation (1).
Summary effects of the 2008 cold spell on mortality
A series of Bayesian meta-analysis models were used to estimate the summary effects of the 2008 cold spell on mortality in all study communities and different regions [31, 32]. The combined CER of the 2008 cold effect among different communities was defined as summary CER. Multivariate meta-analysis models were performed to estimate the summary lag structure of cold effects.
Calculation of excess deaths related to the 2008 cold spell
According to the sixth national population census  and method recommended by World Health Organization (WHO) , the number of excess deaths in each geographical region was calculated by the following formula:
where ED is the number of excess deaths attributable to the 2008 cold spell; ER is the cumulative ER of the 2008 cold spell on total mortality; P represents average daily mortality during non-cold spell days; n is the number of 2008 cold spell days, and N is the total population in each geographical region. The total number of excess deaths induced by this cold spell was calculated as the sum of excess deaths from 4 geographical regions.
A series of sensitivity analyses were performed to test the robustness of our results. We changed degrees of freedom for year (2–4) and calendar month (2–3). When estimating the summary CER of 2008 cold spell on mortality, the largest CER, the smallest CER and both in all communities were respectively removed from the meta-analyses. In addition, we also changed meta-analysis methods from a random effects model to a fixed effects model.
All statistical tests were two-sided, and p <0.05 was considered statistically significant. We used R software (version 2.15.2; R Development Core Team 2012, http://www.R-project.org/) to analyze the data. The “dlnm” package was used to fit Poisson regression . The “Metafor” package was used to fit Meta-analysis . According to the aims of the present study, we only report the excess mortality related to the 2008 cold spell. The effect of cold spell in 2006, 2007, 2009 and 2010 was shown in Additional file 1.
Table 1 presents the distributions of mortality, latitude and weather statistics during the 2008 cold spell, the cold spells in 2006, 2007, 2009 and 2010, and non-cold spell days in 36 communities. All the communities are located between 27.3 and 36.0 degrees north. The 2008 cold spell persisted from 14 to 36 days in the 36 communities, with an average duration of 26.5 days. Compared to cold spells in other years as well as non-cold spell days, the 2008 cold spell had a higher mortality and relative humidity (RH), and lower temperatures. Similar results were also observed in each geographical region (Figure 2).
Figure 3 illustrates the lag structure for the effect of the 2008 cold spell on non-accidental mortality for all 36 communities. The mortality risk increased to a maximum after three days exposure to the cold spell, approached 1.01 after 7 days, and then leveled off for the next 3 weeks. Similar patterns were also observed in southern, eastern and central China. However, the lag structure during the first week was different in Southwestern China. It rapidly declined in the first 2 days, and then increased to around 1.02 of RR at lag 7 days (Additional file 2). Figures 4 and 5 show the cumulative excess mortality (CER) of the 2008 cold spell on mortality at lag 0–27 days. There was a total of 43.8% (95% CI: 34.8% ~ 53.4%) excess deaths caused by the 2008 cold spell compared to non-cold spell days in all 36 communities, and higher death risks were observed in southern China (CER = 53.2%, 95% CI: 33.1% ~ 76.4%) and central China (CER = 53.8%, 95% CI: 26.1% ~ 87.6%) compared to southwestern China (CER = 36.0%, 95% CI: 17.6% ~ 57.4%) and eastern China (CER = 39.0%, 95% CI: 29.5% ~ 49.2%). The attributable numbers of deaths were 19,628, 56,936, 37,454, and 34,262 in eastern, central, southern and south western China, respectively. In total, 148,279 excess deaths were associated with the 2008 cold spell across 15 subtropical provinces of China.
Estimated CERs after stratification by demographic characteristics, location, causes of death and place of death are provided in Table 2. They were higher for RESP (CER = 61.9%, 95% CI: 41.5% ~ 85.3%) than for CVD (CER = 52.9%, 95% CI: 42.1% ~ 64.5%) or CBD (CER = 54.3%, 95% CI: 36.5% ~ 74.4%), and for females (CER = 48.9%, 95% CI: 38.6% ~ 60.1%) than for males (CER = 38.5%, 95% CI: 29.4% ~ 48.3%). The CERs were statistically significant for all age groups, and much higher for the elderly (CER = 65.4% for population ≥ 85 years compared to CER = 17.9% for those aged 0–64 years). The risk of people dying in emergency rooms was significantly higher (CER = 80.3%, 95% CI: 27.9% ~ 154.3%) than the risk of dying in hospital wards (CER = 22.9%, 95% CI: 5.7% ~ 42.9%) or at home (CER = 44.8%, 95% CI: 35.5% ~ 54.7%). CERs in rural areas were higher than in urban areas for RESP (66.0% versus 53.7%) and CBD (61.1% versus 39.0%), but the reverse for CVD (46.8% versus 65.7%).
Table 3 shows the main effect and added effect of the 2008 cold spell on non-accidental mortality at lag 0–27 days. The main effect (CER = 22.5%, 95% CI: 5.2% ~ 42.6%) was much larger than the added effect (CER = 9.4%, 95% CI: -7.8% ~ 29.9%). This phenomenon was also observed in eastern China (CER = 30.5% versus 4.6%), central China (CER = 45.1% versus 11.2%) and southern China (CER = 20.6% versus 13.7%), but not in southwestern China (CER = 8.1% versus 15.2%).
Summary CERs of the 2008 cold spell on mortality varied from 37.8% to 43.8% for changes in the df for smoothness of year and calendar month (Additional file 3), from 42.2% to 45.1% for the changes of fitness methods in the meta-analyses (Additional file 4), and from 42.2% to 44.6% for different lag days (Additional file 5). These findings indicated that our results were robust.
To understand how extreme low temperatures impact human health is crucial not only for policymakers who develop prevention and intervention strategies for extreme weather events, but also for the community, especially for those where extreme cold spells are uncommon. The 2008 cold spell that occurred in subtropical China provided a unique research opportunity to assess the health impacts of extreme cold spells in subtropical regions. In the present study, we assessed the effects of the 2008 cold spell on mortality, spatial heterogeneity of the effects and demographic modifying factors using data from 36 communities of subtropical China. To our knowledge, this is the first multi-community study to investigate effects of the 2008 cold spell on mortality in subtropical China.
As expected, we found the 2008 cold spell significantly increased mortality risk compared to non-cold spell periods, which was consistent with previous studies [2, 4, 8, 11, 34–37]. For example, an average excess mortality of 12.8% was found during the cold spells in 1979–1997 of Netherlands . The cumulative excess mortality of 75+ age people in Moscow were 9.9% and 8.9% during the two cold spells of 2006 . However, we found that the excess mortalities in previous studies were considerably lower than our findings [4, 11, 36, 37], which was also larger than the effect of cold spell in 2006, 2007, 2009 and 2010 (Additional file 1). There are several possible reasons for this difference. Firstly, the intensity and duration of cold spells varied across different studies, which largely modified cold-related effects [11, 38]. The 2008 cold spell appears to be more intensive and longer than cold spells in previous studies. For instance, the cold spells in 23 Italian cities in 2012 lasted 6–18 days, and the temperature dropped 3.5-7.9°C compared to the reference days . In our study, the average daily mean temperature during the 2008 cold spell dropped 7.3-8.7°C compared to the non-cold spell days, and the mean duration of the 2008 cold spell was 26.5 days, with a range of 14–36 days across all study communities. Secondly, methods and cold spell definition employed to estimate the excess mortality varied in different studies, which can lead to the difference in findings between studies [4, 11, 22, 23]. For instance, Kysely et al. defined a cold spell as a period of days when air temperature did not exceed -3.5°C, and calculated the expected (baseline) number of deaths using the mean annual cycle smoothed by 15-day running means . Thirdly, populations in subtropical regions may be more sensitive to cold spells and lack necessary preparedness because cold spells are not common in these regions.
Our study confirmed previous findings that health impacts of low temperatures appeared to be larger in warmer regions than in colder regions [2, 10, 34, 35]. Populations in southern and central China were found to be at higher risk than those in the eastern and southwestern China. This regional variation could be explained by intensity and duration of the cold spell, socio-demographic characteristics and adaptive capacity [9, 11, 35, 39, 40]. In eastern China, the average temperature in winter was lower than that in the other three regions and people may acclimatize to low temperatures and take more effective adaptive measures against cold temperatures during a cold spell . In contrast, populations in southern China are less well acclimatized to extreme low temperatures and few buildings in this region have heating systems to respond to low temperatures. The reasons for the lower effect of the 2008 cold spell in southwestern China is not clear, and require further study in the future.
We further found that morality risks of the 2008 cold spell were higher for people who died of respiratory illness compared with those where cardiovascular or cerebrovascular illness were the reported causes of death , which is consistent with some previous studies [2, 10, 11, 40, 41]. This may be due to the more rapid spread of infectious diseases, reduced response mechanisms of the upper respiratory tract, suppressed immune responses, exacerbated chronic respiratory diseases, and increment of fibrinogen concentration related to the respiratory infections [11, 38, 42]. Further stratification analyses revealed that higher respiratory and cerebrovascular mortality risks were found in rural areas, and higher cardiovascular mortality risks were observed in urban communities. The possible reason for this difference is that people in rural areas are used to staying indoors during winter, and hence are likely to have higher exposure to a complex pollutant mixture of particulate matter (PM) and other toxic compounds  associated with cigarette and biomass fuels or coals use for heating and cooking . These pollutants have been proved to be associated with some respiratory and cerebrovascular diseases, such as infections, asthma, chronic obstructive pulmonary disease (COPD) , hypertension , and strokes . In addition, the less air condition use in rural area could increase people’s vulnerability to extremely cold weather . The reasons for higher CVD mortality risks in urban areas may be partially attributable to the higher prevalence of coronary heart diseases in urban areas of China . Extreme low temperature exposure could increase the cardiovascular stress through changes in blood pressure, vasoconstriction, increase in blood viscosity and levels of red blood cell count, plasma cholesterol and plasma fibrinogen [4, 38, 42], hence increasing the risk of death for people with chronic cardiovascular diseases during the cold spell.
In the current study, we confirmed that age was another important modifier of the temperature-mortality relationship [2, 10]. The elderly were more susceptible to low temperatures than younger people, and this can be attributed to factors such as limited thermoregulation and concomitant chronic disease . In addition, we found that place of death with the highest risk was emergency rooms in urban cities, while in rural areas, the place of death with the highest risk was home. The possible reason is that in urban areas, people are sent to an emergency room quickly and with more severe diseases during an extreme cold spell, while people in rural areas are unable to access emergency services when they are ill. In addition, influenced by the Chinese tradition, people who are almost dying in rural areas tend to stay at home rather than go to hospital. These results indicate that increasing access to emergency services and early preparedness for health facilities during extreme cold days are crucial for preventing cold-related health effects in rural areas.
We also identified that females had higher mortality risk than males, a finding similar to some previous studies [35, 39, 49]. This phenomenon may be mainly due to a higher proportion of women in the elderly because women have a higher life expectancy than men in China, leading to a much larger proportion of vulnerable women . In addition, males usually have better thermoregulation ability than females due to the physiological differences , such as sex hormones, body water regulation and exercise capacity . However, some other studies found inverse effects between males and females [40, 41], indicating that there may be some other specific reasons behind sex, such as cultural practices, daily activities and exposure level and hence more studies of this differential vulnerability may be required in the future.
Few previous studies assessed the number of excess deaths of an extreme weather event [12, 25, 52]. This study observed that 148,279 deaths were attributable to the 2008 cold spell in subtropical China. However, the health effect of the 2008 cold spell was underestimated because we did not include morbidity (including injury) due to this event [13, 14]. Our findings found that health impacts of the 2008 cold spell were very large, suggesting that adapting to health effects of cold spells in subtropical regions is necessary in the context of future climate change.
In order to deepen our understanding of the effects of the 2008 cold spell, we further estimated the main and added effects of the 2008 cold spell. We found that cold effects were mostly attributable to extreme low temperatures rather than any added effect related to sustained duration of the low temperatures, which is similar to Huang et al.’s finding that added effects may contribute small proportions in the total effects of cold spells on CVD mortality in Brisbane . This finding suggests that it is the initial cold ‘snap’ that is the main effect and that as the event continued over time people took protective action. This implies that it is necessary to evaluate potential additional effects from cold spells by decomposing the cold effect into a temperature term and an added effect term so we can better understand the mechanism of the cold effects and better plan public health interventions or better estimate the future burden of temperature-related deaths under predicted climate change scenarios.
Some limitations of this study should be mentioned. Firstly, the meteorological variables used were a simple mean of data collected from various monitoring stations, and the variance of measurement may differ from station to station, which may induce larger heterogeneity of the results. Secondly, air pollutants were not adjusted for in this study, which may cause overestimation of the cold effects. Finally, mortality fluctuations in the winter season can partially be attributed to seasonal patterns of illness such as influenza . Unfortunately, influenza data were not available in this study. However, our previous study conducted in Guangdong Province of China indicated that the association between the 2008 cold spell and mortality was not significantly changed after adjusting for influenza death .
The 2008 cold spell significantly increased mortality in subtropical China, which was mainly explained by independent effects of low temperature rather than consecutive days of cold periods. The cold effect was spatially heterogeneous and modified by individual-specific factors. Considering increased weather variability and changes in frequency and intensity of extreme weather events in light of climate change in the future, the findings of the current study imply that it is necessary to develop adaptive plans for cold spells even in subtropical regions, where populations are generally most acclimatized to hot weather both biologically and behaviorally but more sensitive to cold weather due to lower adaptive capacity.
Climate and health impact national assessment study
Distributed lag linear non-linear model
The intergovernmental panel on climate change
China’s disease surveillance points system
International classification of diseases 10th revision
Daily average temperature
Average wind speed
Akaike information criterion
Day of the week
Cumulative excess risk
World health organization
Chronic obstructive pulmonary disease.
Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM: IPCC, 2013: climate change 2013: the physical science basis: contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. 2013, Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press
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.
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: Short-term effects of the 2008 cold spell on mortality in three subtropical cities in Guangdong Province, China. Environ Health Perspect. 2013, 121: 210-216.
Kysely J, Pokorna L, Kyncl J, Kriz B: Excess cardiovascular mortality associated with cold spells in the Czech Republic. BMC Public Health. 2009, 9: 19-
Whitman S, Good G, Donoghue ER, Benbow N, Shou W, Mou S: Mortality in Chicago attributed to the July 1995 heat wave. Am J Public Health. 1997, 87: 1515-1518.
Fouillet A, Rey G, Laurent F, Pavillon G, Bellec S, Guihenneuc-Jouyaux C, Clavel J, Jougla E, Hemon D: Excess mortality related to the August 2003 heat wave in France. Int Arch Occup Environ Health. 2006, 80: 16-24.
Huang W, Kan H, Kovats S: The impact of the 2003 heat wave on mortality in Shanghai, China. Sci Total Environ. 2010, 408: 2418-2420.
Cagle A, Hubbard R: Cold-related cardiac mortality in King County, Washington, USA 1980–2001. Ann Hum Biol. 2005, 32: 525-537.
Healy JD: Excess winter mortality in Europe: a cross country analysis identifying key risk factors. J Epidemiol Community Health. 2003, 57: 784-789.
Anderson BG, Bell ML: Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009, 20: 205-213.
De’donato FK, Leone M, Noce D, Davoli M, Michelozzi P: The impact of the February 2012 cold spell on health in Italy using surveillance data. PloS one. 2013, 8: e61720-
Gasparrini A, Armstrong B: The impact of heat waves on mortality. Epidemiology. 2011, 22: 68-73.
Zhou B, Gu L, Ding Y, Shao L, Wu Z, Yang X, Li C, Li Z, Wang X, Cao Y: The great 2008 Chinese ice storm: its socioeconomic-ecological impact and sustainability lessons learned. Bull Am Meteorol Soc. 2011, 92: 47-60.
Shi P, Liu L, Wang J, Xu W, Fang W, Wang M: Experiences and lessons of large-scale disaster governance in China—perspective to the response of wenchuan earthquake disaster. International Human Dimensions Program (IHDP) 2009. 2009, Bonn: Open Meeting, 26-30.
Sitnikov IG: Principal weather systems in subtropical and tropical zones. Environmental Structure and Function: Climate System. 2014,http://www.eolss.net/Sample-Chapters/C01/E4-03-04-03.pdf,
Guo Y, Jiang F, Peng L, Zhang J, Geng F, Xu J, Zhen C, Shen X, Tong S: The association between cold spells and pediatric outpatient visits for asthma in Shanghai, China. PloS one. 2012, 7: e42232-
Ma W, Yang C, Chu C, Li T, Tan J, Kan H: The impact of the 2008 cold spell on mortality in Shanghai, China. Int J Biometeorol. 2013, 57: 179-184.
McCarthy JJ, Canziani OF, Leary NA, Dokken DJ, White KS: Climate change 2001: impacts, adaptation, and vulnerability: contribution of working group II to the third assessment report of the intergovernmental panel on climate change. 2001, Cambridge: Cambridge University Press
Van Donkelaar A, Martin RV, Brauer M, Kahn R, Levy R, Verduzco C, Villeneuve PJ: Global estimates of ambient fine particulate matter concentrations from satellite-based aerosol optical depth: development and application. Environ Health Perspect. 2010, 118: 847-855.
Population Census Office under the State Council, Department of Population and Employment Statistics, National Bureau of Statistics: Tabulation on the 2010 population census of people’s Republic of China. 2012, Beijing: China Statistics Press
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.
Menzel A, Seifert H, Estrella N: Effects of recent warm and cold spells on European plant phenology. Int J Biometeorol. 2011, 55: 921-932.
Hickey K: The historic record of cold spells in Ireland. Ir Geogr. 2011, 44: 303-321.
Guo Y, Punnasiri K, Tong S: Effects of temperature on mortality in Chiang Mai city. Thailand: a time series study Environmental Health. 2012, 11: 36-
Gasparrini A: Distributed lag linear and non-linear models in R: the package dlnm. J Stat Softw. 2011, 43: 1-
Zeng W, Lao X, Shannon R, Xu Y, Xu X, Lin H, Liu T, Luo Y, Xiao J, Hu M, Chu C, Ma W: The effect of heat waves on mortality and effect modifiers in four communities of Guangdong Province, China. Sci Total Environ. 2014, 482–483: 214-221.
Armstrong B: Models for the relationship between ambient temperature and daily mortality. Epidemiology. 2006, 17: 624-631.
Peng RD, Dominici F, Louis TA: Model choice in time series studies of air pollution and mortality. J R Stat Soc. 2006, 169: 179-203.
Liu T, Li TT, Zhang YH, Xu YJ, Lao XQ, Rutherford S, Chu C, Luo Y, Zhu Q, Xu XJ, Xie HY, Liu ZR, Ma WJ: The short-term effect of ambient ozone on mortality is modified by temperature in Guangzhou, China. Atmos Environ. 2013, 76: 59-67.
Huang C, Barnett AG, Wang X, Tong S: Effects of extreme temperatures on years of life lost for cardiovascular deaths: a time series study in Brisbane, Australia. Circulation Cardiovascular quality and outcomes. 2012, 5: 609-614.
Viechtbauer W: Conducting meta-analyses in R with the metafor package. J Stat Softw. 2010, 36: 6-14.
Wu W, Xiao Y, Li G, Zeng W, Lin H, Rutherford S, Xu Y, Luo Y, Xu X, Chu C, Ma W: 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-362.
Campbell-Lendrum D, Woodruff R, Prüss-Üstün A, Corvalán C: Climate change: quantifying the health impact at national and local levels. (WHO Environmental Burden of Disease Series No. 14). Edited by: Prüss-üstün A, Corvalán C. 2007, Geneva: World Health Organization
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.
Barnett AG, Dobson AJ, McElduff P, Salomaa V, Kuulasmaa K, Sans S: Cold periods and coronary events: an analysis of populations worldwide. J Epidemiol Community Health. 2005, 59: 551-557.
Huynen MM, Martens P, Schram D, Weijenberg MP, Kunst AE: The impact of heat waves and cold spells on mortality rates in the Dutch population. Environ Health Perspect. 2001, 109: 463-470.
Revich B, Shaposhnikov D: Excess mortality during heat waves and cold spells in Moscow, Russia. Occup Environ Med. 2008, 65: 691-696.
Iñiguez C, Ballester F, Ferrandiz J, Pérez-Hoyos S, Sáez M, López A: Relation between temperature and mortality in thirteen Spanish cities. Int J Environ Res Public Health. 2010, 7: 3196-3210.
Schwartz J: Who is sensitive to extremes of temperature?: a case-only analysis. Epidemiology. 2005, 16: 67-72.
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.
The Eurowinter Group: 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.
Keatinge WR: Winter mortality and its causes. Int J Circumpolar Health. 2002, 61: 292-299.
Smith KR, Mehta S, Maeusezahl-Feuz M: Indoor air pollution from household use of solid fuels. Comparative quantification of health risks: global and regional burden of disease attributable to selected major risk factors. Edited by: Ezzati M, Lopez AD, Rodgers A, Murray CJL. 2004, Geneva: World Health Organization
Mestl HE, Aunan K, Seip HM, Wang S, Zhao Y, Zhang D: Urban and rural exposure to indoor air pollution from domestic biomass and coal burning across China. Sci Total Environ. 2007, 377: 12-26.
Baumgartner J, Schauer JJ, Ezzati M, Lu L, Cheng C, Patz JA, Bautista LE: Indoor air pollution and blood pressure in adult women living in rural China. Environ Health Perspect. 2011, 119: 1390-1395.
Zhang Z, Yu S-Z, Zhou G: Indoor air pollution of coal fumes as a risk factor of stroke, Shanghai. Am J Public Health. 1988, 78: 975-977.
Liu T, Xu YJ, Zhang YH, Yan QH, Song XL, Xie HY, Luo Y, Rutherford S, Chu C, Lin HL: Associations between risk perception, spontaneous adaptation behavior to heat waves and heatstroke in Guangdong province, China. BMC Public Health. 2013, 13: 913-
Zhang XH, Lu ZL, Liu L: Coronary heart disease in China. Heart. 2008, 94: 1126-1131.
Wilkinson P, Pattenden S, Armstrong B, Fletcher A, Kovats RS, Mangtani P, McMichael AJ: Vulnerability to winter mortality in elderly people in Britain: population based study. BMJ. 2004, 329: 647-
Marchand I, Johnson D, Montgomery D, Brisson GR, Perrault H: Gender differences in temperature and vascular characteristics during exercise recovery. Can J Appl Physiol. 2001, 26: 425-441.
Kaciuba-Uscilko H, Grucza R: Gender differences in thermoregulation. Curr Opin Clin Nutr Metab Care. 2001, 4: 533-536.
Huang DS, Zhang SQ: Valuation of the health effects of particulate matter pollution in the Pearl River Delta, China. Book valuation of the health effects of particulate matter pollution in the Pearl River Delta, China: City. 2013,http://www.cerdi.org/uploads/sfCmsContent/html/323/Huang.pdf,
Phu Pin S, Golmard JL, Cotto E, Rothan-Tondeur M, Chami K, Piette F: Excess winter mortality in France: influence of temperature, influenza like illness, and residential care status. J Am Med Dir Assoc. 2012, 13: 309-e301-307
This work was supported in part by the Adapting to Climate Change in China project (ACCC), funded by the UK Department for International Development, the Swiss Agency for Development and Cooperation, and the UK Department for Energy and Climate Change, and Natural Science Foundation of Guangdong Province (S2013010014670).
The authors declare they have no actual or potential competing financial interests.
MGZ participated in the design of this study and drafted the manuscript. LJW participated in the design of this study, collected the data and analyzed the data. TL participated in the design of this study, analyzed the data and drafted the manuscript. YHZ participated in the design of this study. HLL participated in the design of the study and draft revision. YL collected the data and analyzed the data. JPX collected the data and analyzed the data. WLZ collected the data and analyzed the data. XFW provided some data. YWZ collected the data. XG collected the data. SR participated in the design of the study and draft revision. CC participated in the design of the study and draft revision. WJM conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.
Mai Geng Zhou, Li Jun Wang contributed equally to this work.
Electronic supplementary material
Additional file 1:Summary of cumulative excess risks (CER, %) of the cold spells in 2006, 2007, 2009 and 2010 on non-accidental mortality at lag 0–27 in 36 communities of subtropical China, by cause of death, age, gender and place of death.(DOCX 14 KB)
Additional file 2:Summary single day RRs (95% CI) of the 2008 cold spell on non-accidental mortality along lag 0–27 days in four geographical regions of subtropical China. Note: All results were adjusted for secular trend, wind speed, day of week and relative humidity. (DOCX 281 KB)
Additional file 3:Sensitivity analyses of df per year/month on the associations between 2008 cold spell and mortality along lag 0–27 days in 36 communities of subtropical China. Note: All results were adjusted for secular trend, wind speed, day of week and relative humidity. (DOCX 143 KB)
Additional file 4:Sensitivity analyses of meta analysis methods on the association between 2008 cold spell and all mortality along lag 0–27 days in 36 communities of subtropical China. Note: All results were adjusted for secular trend, wind speed, day of week and relative humidity. A: Summary CER was estimated by random effect model. B: Summary CER was estimated by fixed effect model. C: Summary CER was estimated by random effect model with removing the largest ER in the total 36 selected communities. D: Summary CER was estimated by random effect model with removing the smallest ER in the total 36 selected communities. E: Summary CER was estimated by random effect model with simultaneously removing the largest and smallest ERs in the total 36 selected communities. (DOCX 43 KB)
Additional file 5:Sensitivity analyses on the association between 2008 cold spell and all mortality along different lag days in 36 communities of subtropical China. Note: All results were adjusted for secular trend, wind speed, day of week and relative humidity. (DOCX 144 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zhou, M.G., Wang, L.J., Liu, T. et al. Health impact of the 2008 cold spell on mortality in subtropical China: the climate and health impact national assessment study (CHINAs). Environ Health 13, 60 (2014). https://doi.org/10.1186/1476-069X-13-60