Spatial variability of climate effects on ischemic heart disease hospitalization rates for the period 1989-2006 in Quebec, Canada
© Bayentin et al; licensee BioMed Central Ltd. 2010
Received: 27 November 2009
Accepted: 8 February 2010
Published: 8 February 2010
Studies have suggested an association between climate variables and circulatory diseases. The short-term effect of climate conditions on the incidence of ischemic heart disease (IHD) over the 1989-2006 period was examined for Quebec's 18 health regions.
Analyses were carried out for two age groups. A GAM statistical model, that blends the properties of generalized linear models with additive models, was used to fit the standardized daily hospitalization rates for IHD and their relationship with climatic conditions up to two weeks prior to the day of admission, controlling for time trends, day of the season and gender.
Results show that, in most of Quebec's regions, cold temperatures during winter months and hot episodes during the summer months are associated with an increase of up to 12% in the daily hospital admission rate for IHD but also show decreased risks in some areas. The risk of hospitalization is higher for men and women of 45-64 years and varies spatially. In most regions, exposure to a continuous period of cold or hot temperature was more harmful than just one isolated day of extreme weather. Men aged 45-64 years showed higher risk levels of IHD than women of the same age group. In most regions, the annual maximum of daily IHD admissions for 65 years old was reached earlier in the season for both genders and both seasons compared to younger age groups. The effects of meteorological variables on the daily IHD admissions rate were more pronounced in regions with high smoking prevalence and high deprivation index.
This study highlights the differential effects of cold and hot periods on IHD in Quebec health regions depending on age, sex, and other factors such as smoking, behaviour and deprivation levels.
A decline in the rates of deaths and of hospital admissions is reported for cardiovascular diseases in Canada , while most risk factors, except smoking have been rising simultaneously throughout the country [2, 3]. Aside from known risk factors, climate variables, especially temperature and relative humidity, are documented to have an association with circulatory diseases [4–6]. The overall mean temperature as well as extreme weather events across the planet are projected to increase, according to future climate simulations . Therefore, the definition of meteorological variables effects on the population health is regarded as important in the governmental policy making process.
In Quebec, a study has confirmed a significant association between the historic temperature and mortality . Most studies have so far focused on mortality and reported a non-linear relation (U, J, or V shaped) between daily mortality and extreme temperatures. The few studies that focused on morbidity reported an increase in daily hospitalizations for cardiovascular and respiratory diseases during extreme weather conditions [5, 9–11]. It has also been reported that temperature effects are observed up to a few days after initial exposure .
In this study we investigate the impact of meteorological variables on the incidence of ischemic heart diseases (IHD, ICD9:410-414.9) in the province of Quebec (Canada) over the period 1989-2006. The main objective is to estimate the effects of climate variables on IHD hospital admissions while taking into account as much as possible confounding variables or risk factors such as gender, age group, material and social deprivation, smoking and the geographic health region. This represents a first step towards establishing future projections for hospital admissions under climate change scenarios.
Health and related data
Data on smoking prevalence by age group and by sex for the periods 1998, 2000-2001 and 2003 were obtained from the Institut national de santé publique du Québec (INSPQ); no comparable data was available for the period 1989-1997. Data on the deprivation index was also provided by INSPQ. This index has two dimensions: material and social [14–16] and is available at the dissemination area level  or its former equivalent, every five years, since 1996. Again no comparable data was available for the period 1989-1995. The indicators for other risk factors (e.g. obesity, diabetes, and hypertension) were unavailable at the spatial scale of interest for this study All preceding indicators were thus not included in the model and they were used for discussion purposes only.
Daily mean temperature
A generalized additive model (GAM)  was used to model climate effects on IDH hospital admissions. Parameter estimation is carried out with the approach combining smoothing parameter λ with the B-splines basis function . The variables to be included in the model are selected by a forward stepwise selection technique and an F-test is conducted in the comparison of two different models . Prior to variables selection, a lag search is performed using the generalized cross-validation (GCV) criterion. For each variable, the mean of past observations up to two weeks is computed and the number of days providing the smallest model GCV is added to the list of covariables. The parameter λ, is estimated by comparing the Akaïke information criterion (AIC) values corresponding to different values of λ (λ = 0.001, 5, 15, 30, 60, 90 and 180). The excess risk is the relative change in the daily rate of IHD admission due to a decrease of 1°C below a temperature threshold for winter and an increase of 1°C above a threshold in summer.
Models for winter
Excess risk (%) due to a 1°C drop below a threshold in winter
65 and up
65 and up
For men, the risk was higher (1.03% to 12.32%) in the 45-64 years age group in most regions, compared to older men (0.53% to 2.98%). There was a protective effect (decreasing risk) of cold temperature on IHD admissions rate in few regions.
The highest excess risks for women were in region Côte-Nord for both age groups. The protective effect of cold was more common for women across the province. Note that for both genders, the regions with the highest excess risk were in the group with the highest IHD annual rate. The low annual rate group also presented the lowest excess risk for men.
Maximum daily rates for men were observed in January or February. For most regions January was the high risk period whereas in a few, it was late February. Regions with February peaks were North-Eastern regions and North-Western regions. Peaks in hospital admissions were observed earlier in the season for men 65 years and over when compared to the 45-64 years age group. For women, maximum rates were observed at the end of December and early January. High rates were also observed earlier for 65 years old and up women compared to younger ones.
Models for summer
Excess risk (%) due to a 1°C rise above a threshold in summer
65 and up
65 and up
During summer, for both genders, higher values of the daily admission rate were recorded earlier in the season for the 65 years and over age group in comparison to the 45-64 years age group. Most regions, especially regions with large cities, presented a decrease in 65 years old men's daily admission rate around mid-July, while an increase is observed in the 45-64 years old men's rate.
Risk factors and health determinants
Men and women smoking prevalence mean (%) by age group and health region (1998, 2000-2001, 2003)*
65 and up
65 and up
Percentage of the region's population that scores in the highest quintile of material and social deprivation indexes* (2001)
ith quintile for material dimension
ith quintile for social dimension
As reflected in the literature, the relationship between hospital admission for IHD and climate is probably a very complex one, modulated by several factors. Nonetheless, a strong relationship was observed in this study between IHD hospital admission rates and climatic conditions of the days preceding the hospitalization date. Delayed variables corresponding to ground snow, precipitation and daily mean temperature were among the first to be selected in the GAM model, confirming a cumulative effect of cold weather on IHD.
Some counter-intuitive findings also showed up. Cold temperatures resulted in a protective effect for women except for most Northern regions, pointing to possible links with protective behaviour (such as staying put at home during a cold spell), or the lower participation of women in the building or similar work sectors with compulsory outdoor exposure even during cold spells. In summer, for men 45-64 years old, a significant effect of hot temperatures was observed only in a few regions, whereas for women of the same age group the excess risk was recorded in most regions. In general the 45-64 age group was more at risk of IHD admission than the group of 65 years and over for both genders and in most regions. However, the older group was affected earlier in the season. Again behavioural factors are probably involved here, such as younger people having to work even in hot weather compared to retirees and seniors with activity limitations, as well as differential access to air conditioning, both already documented for the province [21, 22].
A decline in the effects of meteorological variables on IHD daily admission rates was observed over the period of 1989-2006. This can partly be explained by the changes in surface air temperature over the period of our study. Studies have reported warming of the air temperature in Canada over the last few decades [23, 24]. In Quebec, a warming of the surface air temperature is reported in the southern part of the province over the period of 1960 and 2005, taking place mostly after 1995 . Also the increase trend during that period shows a gradient from South to North and West to East , which is consistent with our data. Also winters have been steadily warmer but summers have yet to become hotter for most regions .
The association between cold temperatures and IHD admission rates takes the same V or U shape as reported in the literature. However that shape was less common for summer models. The increase in daily hospitalizations with extreme temperatures was not observed in all cases (age groups, genders, health regions).
A decrease in hospital admissions was observed in all regions, for both genders and for both age groups, after an initial increase in the nineties. This is very consistent with the national trends of deaths and hospital admissions for heart diseases . Tu et al. (2009) suggested that better control of lipid levels and associated coronary plaque stabilization through increasing rates of statin use in Canada could be significant factors in the decrease, with other medical treatment of risk factors also increasing and potentially involved in the trend . Climate trends should be added as an explanatory variable of importance.
Regions with a significant excess risk also displayed high smoking prevalence and higher deprivation indexes (either material or social). A significant excess risk was also noted for the 65 years and over age group for some regions in which the smoking prevalence is moderate or low. This suggests that, in regions with high excess risk but moderate and low smoking prevalence, others risk factors not considered here are involved in the incidence of IHD. Diet, smoking and physical inactivity, considered modifiable risk factors, also contribute to the development of IHD. Only smoking could be addressed here, although the examination of the deprivation index allows for an assessment of the role of income, education, social support and similar important health determinants. Studies dealing with the impact of comorbidity on IHD, have reported hypertension, diabetes and obesity as known risk factors of cardiovascular diseases. Those risk factors increased in all age groups over the period of 1994-2004 and the highest increase was observed in the lower income groups . Lee et al. (2009) conclude that the rising burden of those risk factors could result in future increases in cardiovascular diseases. They also reported that the smoking prevalence decreased over that period in all Canadian provinces. This can partly explain the drop in hospital admissions in all regions by the beginning of 2000.
The physiological processes that explain the effect of climate of ischemic heart disease are complex and not completely understood. Nonetheless, vitamin D is reported to be associated to the incidence of cardiovascular disease  and the mechanisms of its effect well documented . Grimes, Hindle and Dyer  used the positive slope of the relationship between the concentration of blood cholesterol and latitude, the negative relationship between hours of sunshine per year and the death rate from coronary heart diseases (CHD) to suggest that the geographical variation of CHD is influenced by sunlight. Though we did not use sunlight exposure as an independent variable in our study, the separation of the seasons and the inclusion of the ith day of the season in the models account for the seasonality of sunlight. Also, because of the known relationship between latitude and sunlight and thus, the cutaneous synthesis of vitamin D , the separation of data by health region can be construed as a form of stratification based on geographical coordinates. Figure 4 clearly shows a latitude effect on IHD hospital admission rate, with the annual IHD admission rate decreasing with the latitude.
The proportion of the population that is active is the same in all regions and cannot be considered as an explanatory variable for IHD hospitalisations, but the proportion of workers involved outdoors likely varies from region to region; unfortunately no data were available on this topic.
Limitations of the study
Our study has some limitations. Because the health data were from existing administrative databases, patient history of heart diseases, personal habits such as smoking and comorbidity could not be taken into account. The limitations of data on smoking and the deprivation index allowed only for a qualitative analysis. Another limitation of our study is our inability to distinguish between new admissions and repeated ones, which could lead to increased rates of hospitalisation. Our study did not assess either the role of air pollution in the genesis of IHD episodes, although such a role seems significant [31–33].
While these results should be considered preliminary, they nonetheless point out important messages for physicians throughout the country. Weather extremes do play a significant role in the genesis of acute and severe episodes of IHD, and the longer the heat wave or the cold spell, the more likely the public health impact. Several already vulnerable subgroups, based on a variety of indicators, seem to present more sensitivity to weather extremes, but adaptive behaviours can modify this excess risk. Adapted treatment regimes by physicians could probably help improve the score, although much more research at the individual level is warranted before reaching the stage of formal recommendations. However, older people and lower income groups could already benefit from supplementary medical attention in periods of meteorological extremes.
We thank Christian Charron and Steve Toutant for their assistance in the extraction of data.
This project was funded by Ouranos, Santé Canada, the Institut National de Santé Publique du Québec and the Ministère de la santé et des services sociaux du Québec.
- Tu JV, Nardi L, Fang J, Liu J, Khalid L, Johansen H, for the Canadian Cardiovascular Outcomes Research Team: National trends in rates of death and hospital admissions related to acute myocardial infarction, heart failure and stroke, 1994-2004. CMAJ. 2009, 180: E118-125.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee DS, Chiu M, Manuel DG, Tu K, Wang X, Austin PC, Mattern MY, Mitiku TF, Svenson LW, Putnam W, Flanagan WM, Tu JV, Canadian Cardiovascular Outcomes Research Team: Trends in risk factors for cardiovascular disease in Canada: temporal, socio-demographic and geographic factors. CMAJ. 2009, 181: E55-66.PubMedPubMed CentralView ArticleGoogle Scholar
- Dai S, Bancej C, Bienek A, Walsh P, Stewart P, Weilgosz A: Tracking heart disease and stroke in Canada 2009. Chronic Dis Can. 2009, 29: 192-193.Google Scholar
- Lin S, Luo M, Walker RJ, Liu X, Hwang S-A, Chinery R: Extreme High Temperatures and Hospital Admissions for Respiratory and Cardiovascular Diseases. Epidemiology. 2009, 20: 738-746. 10.1097/EDE.0b013e3181ad5522.PubMedView ArticleGoogle Scholar
- Saez M, Sunyer J, Tobias A, Ballester S, Antó JM: Ischemic heart disease mortality and weather temperature in Barcelona, Spain. Eur J Public Health. 2000, 10: 58-63. 10.1093/eurpub/10.1.58.View ArticleGoogle Scholar
- Chang CL, Shipley M, Marmot M, Poulter N: Lower ambient temperature was associated with an increased risk of hospitalization for stroke and acute myocardial infarction in young women. J Clin Epidemiol. 2004, 57: 749-757. 10.1016/j.jclinepi.2003.10.016.PubMedView ArticleGoogle Scholar
- Schnoor JL: The IPCC fourth assessment. Environ Sci Technol. 2007, 41: 1503-10.1021/es072475x.PubMedView ArticleGoogle Scholar
- Doyon B, Belanger D, Gosselin P: The potential impact of climate change on annual and seasonal mortality for three cities in Quebec, Canada. Int J Health Geogr. 2008, 7: 23-10.1186/1476-072X-7-23.PubMedPubMed CentralView ArticleGoogle Scholar
- Shao L, Luo M, Walker R, Liu X, Hwang S, Robert C: Impact of Hot Weather Conditions on Respiratory and Cardiovascular Hospital Admissions in New York City, USA. Epidemiology. 2008, 19: S302-S303. 10.1097/EDE.0b013e3181632c3d.View ArticleGoogle Scholar
- Tiina MMk, Raija J, Jari J, Terttu HH, Ari P, Aini B, Sylvi S-K, Maija L, Juhani H: Cold temperature and low humidity are associated with increased occurrence of respiratory tract infections. Respir Med. 2009, 103: 456-462. 10.1016/j.rmed.2008.09.011.View ArticleGoogle Scholar
- Friger M, Yackerson N, Bolotin A, Kordysh E: Meteorological Factors Influence on Hospitalization for Respiratory Diseases And Symptoms in the South Israel. Epidemiology. 2006, 17: S423-S424. 10.1097/00001648-200611001-01133.View ArticleGoogle Scholar
- Schwartz J, Samet JM, Patz JA: Hospital admissions for heart disease: the effects of temperature and humidity. Epidemiology. 2004, 15: 755-761. 10.1097/01.ede.0000134875.15919.0f.PubMedView ArticleGoogle Scholar
- Anderson RN, Rosenberg HM: Age standardisation of Death rates: Implementation of the year 2000 standard. Natl Vital Stat Rep. 1998, 47: 1-16.Google Scholar
- Pampalon R, Hamel D, Gamache P, Raymond G: A deprivation index for health planning in Canada. Chronic Dis Can. 2009, 29: 178-191.PubMedGoogle Scholar
- Morin A: Receuil statistique sur la pauvreté et les inégalités socioéconomiques au Québec. 2006, Québec: ISQ, MESS, 134-Google Scholar
- Pampalon R, Raymond G: A deprivation index for health and welfare planning in Quebec. Chronic Dis Can. 2000, 21: 104-113.PubMedGoogle Scholar
- Puderer H: Introducing the dissemination area for 2001 census: an uptade. Edited by: Canada S. 2001, Ottawa: Geography Working Paper Series, 2000: 11-Google Scholar
- Environnement Canada: Canada's National Climate Archive. 2008, Google Scholar
- Hastie TJ, Tibshirani RJ: Generalized Additive Models. 1990, New-York: Chapman & Hall, 1Google Scholar
- Eilers PHC, Marx BD: Flexible Smoothing with B-splines and Penalties. Statist Sci. 1996, 11: 89-120. 10.1214/ss/1038425655.View ArticleGoogle Scholar
- Bélanger D, Gosselin P, Valois P, Abdous B: Vagues de chaleur au Québec méridional: adaptations actuelles et suggestions d'adaptations futures. 2006, Québec: Institut national de santé publique du Québec,Google Scholar
- Bélanger D, Gosselin P, Valois P, Abdous B: Vagues de froid au Québec méridional: adaptations actuelles et suggestions d'adaptations futures. 2006, Québec: Institut national de santé publique du Québec,Google Scholar
- Bonsal BR, Zhang X, Vincent LA, Hogg WD: Characteristics of daily and extreme temperatures over Canada. J Clim. 2001, 14: 1959-1976. 10.1175/1520-0442(2001)014<1959:CODAET>2.0.CO;2.View ArticleGoogle Scholar
- Zhang XB, Vincent LA, Hogg WD, Niitsoo A: Temperature and precipitation trends in Canada during the 20th century. Atmosphere-Ocean. 2000, 38: 395-429.View ArticleGoogle Scholar
- Yagouti A, Boulet G, Vincent L, Vescovi L, Mekis E: Observed changes in daily temperature and precipitation indices for southern Quebec, 1960-2005. Atmosphere-Ocean. 2008, 46: 243-256. 10.3137/ao.460204.View ArticleGoogle Scholar
- Lemmen DS, Warren FJ, Lacroix J, Bush E, Lemmen , Warren : Vivre avec les changements climatiques au Canada: édition 2007. 2008, Ottawa (Ontario): Gouvernement du Canada, 448-View ArticleGoogle Scholar
- Wang TJ, Pencina MJ, Booth SL, Jacques PF, Ingelsson E, Lanier K, Benjamin EJ, D'Agostino RB, Wolf M, Vasan RS: Vitamin D deficiency and risk of cardiovascular disease. Circulation. 2008, 117 (He): 503-11. 10.1161/CIRCULATIONAHA.107.706127.PubMedView ArticleGoogle Scholar
- Zittermann A, Schleithoff SS, Koerfer R: Putting cardiovascular disease and vitamin D insufficiency into perspective. Br J Nutr. 2005, 94 (He): 483-92. 10.1079/BJN20051544.PubMedView ArticleGoogle Scholar
- Grimes DS, Hindle E, Dyer T: Sunlight, cholesterol and coronary heart disease. QJM. 1996, 89 (8): 579-89.PubMedView ArticleGoogle Scholar
- Webb AR, Kline L, Holick MF: Influence of season and latitude on the cutaneous synthesis of vitamin D3: exposure to winter sunlight in Boston and Edmonton will not promote vitamin D3 synthesis in human skin. J Clin Endocrinol Metab. 1988, 67 (2): 373-8. 10.1210/jcem-67-2-373.PubMedView ArticleGoogle Scholar
- Tagaris E, Liao KJ, Delucia AJ, Deck L, Amar P, Russell AG: Potential Impact of Climate Change on Air Pollution-Related Human Health Effects. Environ Sci Technol. 2009, 43: 4979-4988. 10.1021/es803650w.PubMedView ArticleGoogle Scholar
- Zanobetti A, Schwartz J: The Effect of Fine and Coarse Particulate Air Pollution on Mortality: A National Analysis. Environ Health Perspect. 2009, 117: 898-903.PubMedPubMed CentralView ArticleGoogle Scholar
- Cendon S, Pereira LA, Braga AL, Conceicao GM, Cury Junior A, Romaldini H, Lopes AC, Saldiva PH: Air pollution effects on myocardial infarction. Rev Saude Publica. 2006, 40: 414-419. 10.1590/S0034-89102006000300008.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.