A multilevel analysis to explain self-reported adverse health effects and adaptation to urban heat: a cross-sectional survey in the deprived areas of 9 Canadian cities
BMC Public Health volume 16, Article number: 144 (2016)
This study identifies the characteristics and perceptions related to the individual, the dwelling and the neighbourhood of residence associated with the prevalence of self-reported adverse health impacts and an adaptation index when it is very hot and humid in summer in the most disadvantaged sectors of the nine most populous cities of Québec, Canada, in 2011.
The study uses a cross-sectional design and a stratified representative sample; 3485 people (individual-level) were interviewed in their residence. They lived in 1647 buildings (building-level) in 87 most materially and socially disadvantaged census dissemination areas (DA-level). Multilevel analysis was used to perform 3-level models nested one in the other to examine individual impacts as well as the adaptation index.
For the prevalence of impacts, which is 46 %, the logistic model includes 13 individual-level indicators (including air conditioning and the adaptation index) and 1 building-level indicator. For the adaptation index, with values ranging from -3 to +3, the linear model has 10 individual-level indicators, 1 building-level indicator and 2 DA-level indicators. Of all these indicators, 9 were associated to the prevalence of impacts only and 8 to the adaptation index only.
This 3-level analysis shows the differential importance of the characteristics of residents, buildings and their surroundings on self-reported adverse health impacts and on adaptation (other than air conditioning) under hot and humid summer conditions. It also identifies indicators specific to impacts or adaptation. People with negative health impacts from heat rely more on adaptation strategies while low physical activity and good dwelling/building insulation lead to lower adaptation. Better neighbourhood walkability favors adaptations other than air conditioning. Thus, adaptation to heat in these neighbourhoods seems reactive rather than preventive. These first multi-level insights pave the way for the development of a theoretical framework of the process from heat exposure to impacts and adaptation for research, surveillance and public health interventions at all relevant levels.
Health impacts attributed to periods of high summer temperatures are undeniable in all latitudes . To date, several scientific publications have reported on emergency room visits and hospitalizations, mainly by assessing medical consultation records. The risk factors usually refer to individual characteristics, possibly due to the type of data available in health system administrative records. Neighbourhood [2, 3], building or dwelling [4, 5] characteristics as well as coping behaviours [6–10] are considered less frequently, despite their close connection with vulnerability to heat .
Our study exploits the perceptual approach to estimating the adverse health impacts of heat exposure as used or recommended elsewhere [9, 12–15]. This approach has the significant advantage of considering the wide range of factors potentially having a negative effect on health but which are difficult to measure directly [12, 16] and which have to be validated [17–19]. A similar phenomenon has been observed from studies using self-reported health behaviour [20, 21]. The perceptual approach has also been used for the evaluation of environmental exposure . In this regard perception of room temperature is a proximal variable of thermal comfort associated with various health problems [12, 17]. This method has recently been described as an alternative to traditional techniques for measuring ambient temperature .
The perceptual approach is therefore a very promising technique for adding value to health and housing studies . Its use for measuring the health impacts of heat and its determinants is highly relevant from a public health perspective, particularly given that periods of high summer heat will intensify in the near future  and that these events are already affecting the residents of the most disadvantaged neighbourhoods in major cities [2, 3, 25] struggling with the urban heat island effect [26, 27]. These materially and socially disadvantaged census dissemination areas (DAs) in Québec are located within the province’s large urban centers [28, 29] and are characterised by a lack of green space and an abundance of asphalt .
Our study pertains to these geographic pockets of urban poverty localised in the hottest spaces or neighbourhoods. We studied the most disadvantaged DAs in the nine cities of Québec with a population of at least 100 000 in 2010, which constituted 47 % of the total population . In this paper we aim to explain (1) self-reported adverse health impacts and (2) self-reported adaptations other than use of air conditioning at home during very hot and humid summer episodes using a multilevel approach. Dwelling DAs by way of neighbourhoods (DA-level), occupied buildings (building-level), and individuals (individual-level) represent for the analysis three-tiers or levels nested one inside the other. This multilevel approach has been used successfully in the field of health [31, 32] but to our knowledge has not yet been applied to heat-health related problems.
A cross-sectional design with a stratified sample was chosen for this study. A protocol for the study (including a questionnaire) received approval by the research ethics committee of the Centre hospitalier universitaire de Québec prior to commencement of field work.
A two-step selection technique adapted from Vallée et al.  was used to identify representative samples in each of the nine cities. The number of highly disadvantaged DAs per city was identified based on a widely used Canadian deprivation index . Each such DA included minimally one public low-rental housing building. This first step also identified the total number of households to be included in the survey per city. The second step consisted of randomly selecting households to be approached. One adult per household was interviewed.
A total of 3485 individuals living in 1647 buildings and 87 DAs were interviewed. While the overall response rate was 19 %; another statistic is the response rate by question, which was much higher at 95 % or more for more than 90 % of questions asked (for more details, see [35, 36].
Data collection and measures
Questionnaire-based interviews were carried out by 21 interviewers who had each received two days of training. The face-to-face interviews took place in the participant’s home. The development of the draft questionnaire was inspired by relevant population survey tools such as the Statistics Canada’s 2007 questionnaire for the Canadian Community Health Survey and was pre-tested and adjusted prior to data collection. Participants received $10 as payment for their time. For more information see [35, 36].
The first of two dependent variables was based on a proxy variable of participant response to a question about whether “their physical health was negatively affected during very hot and very humid summer conditions”. This variable, referred to as ‘self-reported adverse health impacts’, reflects the self-reported state of overall health in the context of summer heat exposure. A separate and similar question was administered pertaining to their self-reported state of “mental health” in the same context. The preamble and questions were adapted from the Canadian Community Health Surveys going back to 2001 . Both of these key questions were grouped and a binary variable was used to summarize responses. Participants reporting their physical and/or mental health as moderately or greatly adversely affected by very hot and humid weather were considered as the risk group (vs. slightly or not at all).
Second, we constructed an individual adaptation index to very hot and humid summer conditions using Multiple Correspondence Analysis [38, 39]. The heat adaptation index covers 14 easy-to-use and energy-efficient adaptations for reducing the perception of heat (cooling off) or seeking protection from the sun, regardless of location (indoors/outdoors/at home or other). These adaptations are: when the sun is shining, (1) wearing a head covering, (2) closing curtains to keep the home cool; to cool off, (3) taking showers or baths more often than usual, (4) sponging the face and neck, (5) eating iced foods, (6) swimming in a public pool, lake or river, (7) frequenting places other than home for air conditioning, (8) drinking tap water as a main drink, (9) using the balcony and (10) the yard in the evening; to reduce heat sources in the home, (11) switching off the computer when not in use, (12) decreasing dryer and (13) oven use; finally, (14) adopting preventive behaviour according to weather information transmitted through the media or on the Internet (binary variables: 0 = never/rarely, 1 = occasionally/often/always). Seventy-five percent of the total inertia on dimension 1 is explained by these 14 variables (dimension 2: 5 %). Index values range from about –3 to +3 and are almost normally distributed in the case of Dimension 1 more technical details on MCA methods and the proposed index might be found in reference  which explains the development of the index.
The choice of independent variables was based on the scientific literature on human health and oppressive heat [41, 42]. In the current study, they reflect the three levels nested one within the other and considered within the multilevel analysis. The individual-level indicators concern heat exposure in the neighbourhood of residence, the building and the dwelling, air conditioning, index adaptation and city of residence as classified according to the average temperature over the last 30 years  (Additional file 1: Table S1). They also correspond to the states or conditions that increase heat sensitivity and to some indicators related to the access to treatment. For the building-level indicators, the averages of some individual-level variables for each building sampled were used in the absence of census data on buildings. Finally, the DA-level variables were evaluated using averages of some individual-level variables and with data available by DA at the Institut national de santé publique du Québec, including an indicator of intra-urban heat islands  and an index walkability .
F-test, t and Chi-square tests together with classical bivariate analyses were performed. In terms of software, survey procedures such as proc surveylogistic in SAS 9.3 were utilised. To make them more concise, bivariate analyses are presented only for the variables used also in multivariate model.
In the present multivariate analysis, the study of the self-reported adverse health impacts of heat is different from previous papers from the same dataset as we used a multilevel approach rather than generalized estimating equations (GEE) methods [35, 36]. As reported by Chaix and Chauvin and contrary to GEE, multilevel models have the specific objective to explain intergroup variation by higher level variables (here building and DA levels); the contribution of these covariables is thus presented here to explain the self-reported adverse health effects while only individual-level covariables were used in previous GEE analyses. We also added here several individual-level covariables related to heat adaptation (other than air conditioning at home) which were not present in previous analyses [35, 36]. Otherwise, the only previous article on the individual adaptation index to very hot and humid summer conditions  concerned its development and criterion-related validity. The present analysis thus concerns the explanation of the index in a multivariate model within a multilevel approach, including individual-level, building-level and 6 DA-level covariables.
More precisely, for the multivariate analysis, we used MLwiN [46, 47] to fit the three-level model (individual-level nested in building-level nested in DA-level). A logistic regression model served to explain the self-reported adverse health impacts (binary variable), while a linear regression model was utilised to explain the heat adaptation index (continuous variable). For both regressions, an ascending regression procedure enabled the authors to identify the main independent variables associated with the dependent variables. This means that we evaluated first if the building-level, then the DA-level, were statistically significant when added to the null-model. Secondly, we added in an ascending regression procedure the individual-level covariables (when statistically significant in bivariate analyses). Finally, the building-level and DA-level covariables contribution to the model were estimated.
The models were evaluated for the influence of the season in which the interview took place. We make sure in this article, by means of tetrachoric (binary variables) or polychoric (variables with ≥ 3 strata) correlation coefficients, that no independent variables presented (taken two at a time) were strongly correlated (r ≥ 0.6). No statistically significant cross-level interactions were observed (p < 0.05). Consequently no interaction terms were used. Moreover, to obtain more interpretable estimates, some of the individual-level variables were centered when the model included both the variable of individual-level and the average of this variable at building-level or DA-level . Hence, the part attributed to the building mean (or DA) was substracted to estimate only the individuals contribution . These centered variables are noted in Tables 1 and 2. Also, the models that examine a covariable’s influence at two levels are known as contextual or compositional models (see  for more details).
For the linear regression model, the residual variance (variance not explained by the model) in the heat adaptation index has been decomposed into three components: the “within buildings between individuals variance”, the “within DA between buildings variance” and the “between DA variance”. We used the Likelihood ratio test to compare nested models and estimated the percentage of variation found at each level by the variance partition coefficients (VPC) . For the logistic regression model, we did not consider partitioning variances for the perceived adverse health impacts, since a simple VPC is not available for binary variables. Variances and residual variance percentages for health impacts to heat, however, were reported by level to give an idea of the contribution of the addition of a model with regard to the previous model. Moreover, to assess if the response prevalence exhibits more or less variation than would be expected under a binomial distribution, we used the approach of multiplicative over-dispersion by adding a multiplicative scale factor (α) to the variance of the binary response (α = 1, then the prevalence is binomially distributed; α > 1 or α < 1, over-dispersion or under-dispersion respectively) .
Explanation of self-reported adverse health impacts when it is very hot and humid in summer
As already reported [35, 36], the prevalence of self-reported adverse health impacts during very hot and humid summer conditions was 46.0 % (confidence interval, CI: 44.2–47.8) for the 3485 respondents living in the very disadvantaged dissemination areas of the nine most populated cities within Québec.
The random part of the logistic regression model (based on the extra-binomial distribution, since α < 1) shows that the residual variances at the building-level and DA-level decrease significantly as the model evolves (Additional file 2: Table S2). At the building-level, residual variance decreases by 33 % with the addition of the DA-level (M000) to the building-level null model (M00), by 51 % with the addition of individual-level covariables (M000:1) to M000, and by 17 % with the addition of a building-level covariable (M000:12) to M000:1. At the DA-level, residual variance decreases by 48 % from M000 to M000:1 and then stabilizes. Therefore, the use of a 3-level model appears justified.
The 3-level model (M000:12) is adjusted for thirteen individual-level covariables (bivariate analyses presented in the Additional file 3: Table S3) and a building-level covariable. It is not influenced by the season in which the interview took place (p = 0.0580), however, nor by the presence of intra-urban heat islands (66 % of visited buildings were located in a heat island and 32 % within 50 metres of one, p = 0.5380), and none of the DA-level covariables were statistically significant.
Of the fourteen covariables that significantly explain (p < 0.05) the prevalence of impacts (M000:12, Table 1), eleven are risk indicators and three are prevention indicators.
The eleven impact-risk indicators are individual-level. They include: self-reported chronic diseases, especially if there is more than one; long-term medical leave or disability; health problems due to the quality of the air inside the dwelling in the opinion of the respondents; perceived stress that is high or very high in most cases; a neighbourhood perceived as somewhat or very polluted due to urban traffic density; access to air conditioning in the dwelling; higher self-reported heat adaptation according to index (this results in a negative coefficient β in the model); the perception that urban development infrastructure must be added in the dwelling neighbourhood in order to improve adaptation to heat; as well as not practicing physical activities in the past three months, being female and between 45 and 64 years old.
The three impact-prevention indicators include: support in the past year by at least two people living within a radius of 80 km from the place of residence but not in the same neighbourhood; and two covariables for temperature satisfaction inside the dwelling in summer, one being individual-level (correlation with air conditioning at home, r = −.04) and the other being building-level. According to the latter, individuals who are most satisfied with the temperature inside their dwelling in summer (coefficients reduced by −0.475 for each additional satisfaction point; so, the odds of adverse health impacts is reduced by 37.8 % for each additional satisfaction point) and the buildings in which households are on average more satisfied (coefficients reduced by −0.375 for each addition of a satisfaction point; so, the odds of adverse health impacts is reduced by 31.3 % for each additional satisfaction point) are at lower odds of self-reported adverse health impacts to heat.
Explanation of adaptation index when it is very hot and humid in summer
Of the 3485 respondents, 16.6 % (CI 15.3–18.0) were very adaptive, according to the adaptation index (values ≤ −1), 16.7 % (CI 15.3–18.0) were hardly adaptive at all (values ≥ 1), while 66.7 % (CI 65.0–68.4) were between these two poles.
The random part of the linear regression model shows that the 3-level null model is statistically significant (Additional file 4: Table S4). Overall, it also shows that the building-level and DA-level residual variances decrease significantly as the model evolves (see relative changes) and that its adjustment for individual and contextual covariables is significant. Moreover, the 3-step null model explains 83, 14 and 2 % of the total residual variance of the adaptation index to the individual-level, building-level and DA-level, respectively. Residents’ adaptation indices for the same building (and the same DA) are therefore slightly more correlated than if they lived in different buildings.
The 3-level model thus proves to be relevant, but the aggregation rate is low. It is adjusted for ten individual-level covariables (bivariate analyses presented in the Additional file 5: Table S5), one building-level and two DA-level covariables. This model, however, is not influenced by air conditioning at home (p = 0.6748) nor by the season in which the interview took place (p = 0.0797). Of the thirteen covariables that significantly explain (p < 0.05) the adaptation index, four characterize a weaker tendency to adopt behaviours measured by the index (this translates into β positive coefficients in the model), and nine provide incentive to be more adaptive (negative β; Table 2).
The four covariables characterizing a weaker self-reported adaptation according to the index include: not practising physical activity in the past three months, the fact of suffering from at least one functional disability (such as difficulty walking), especially if it is often, and two covariables on satisfaction with the quality of the thermal insulation of the dwelling in summer, one of which is individual-level (correlation with air conditioning at home, r = −.04) and the other building-level. According to the latter, the people most satisfied with the thermal-insulation quality of their dwelling in summer (an increase of 0.79 in the index for each additional point of satisfaction) and the buildings in which households on average are more satisfied (an increase of 0.82 in the index for each additional point of average satisfaction) adopt fewer behaviours measured by the index for adapting to the heat.
Of the nine covariables that encourage more adaptive behaviour according to the index, seven are individual-level and two are DA-level. The first are: being younger than 65 years old, in particular 18 to 44; perceiving the need for more urbanistic infrastructure or services to improve adaptation to heat in the residence’s neighbourhood; feeling the adverse impacts on their health when it is very hot and humid in summer; having met with friends at least once a month and having been supported by at least two persons residing in the same neighbourhood (not the same dwelling) or less than 80 km from the dwelling but not in the same neighbourhood in the last year; and primarily using a car for local travel. The final two reveal that, the lower the residential mobility in a DA (decrease in the index by −0.028 for each additional year to the average duration of residency in the same dwelling) and the better the potential for walking (decrease in the index by −0.013 for each additional walkability point), the more the people who live there adopt behaviours measured by the index of adaptation.
Summary of fixed effects
Among all covariables associated with one dependent variable or another, nine explain only the impact prevalence and eight explain only the adaptation index (Additional file 6: Table S6).
In this study, multilevel analyses made it possible to process building-level and DA-level variance as a source of information , in contrast to other approaches used to analyze the same dataset [35, 36].
The results of random parts of the analyses show that a 3-level model supports the identification of indicators associated with self-reported adverse health impacts and the adaptation index when it is very hot and humid in summer. Therefore, they emphasize the need to develop public health interventions aimed at both residents of very disadvantaged DAs, property managers in these sectors and municipal authorities in large urban centres already struggling with the urban heat island.
The importance of the role of buildings and the environment as risk (or prevention) factors of self-reported health impacts and adaptation to heat is strengthened by the contribution of four building-level and DA-level covariables. However, this number of four variables remains well below the number of individual-level covariables, despite important efforts to measure various aspects of the context (see Additional file 1: Table S1). This observation has been reported in the literature and some explanations have been given on this subject [52–54], including the fact that census data do not characterize neighbourhoods in detail and that a DA is an administrative geographical unit and not a neighbourhood as perceived by respondents.
The effect of buildings and the environment on self-reported adverse health impacts and on the adaptation index could also be covered by some covariables. We cannot make any conclusions regarding this aspect due to the cross-sectional design of our study. Also in the absence of a proven theoretical framework, because although some interesting avenues have been proposed [41, 55, 56], none have been validated. The establishment of such a framework would be very useful for research and public health surveillance, in particular for clarifying the kind of links between covariables; in our study, most are related only to the prevalence of self-reported adverse health impacts or to the adaptation index.
Specifically, in this study, the prevalence of self-reported negative health impacts when it is very hot and humid in summer is 46 %; although high, this could correspond in reality to the highly disadvantaged DAs of the most populous cities in Québec, as discussed elsewhere [35, 36].
As already reported in the literature [41, 57–59], poor physical or mental health is a powerful indicator of risk to heat. Similarly, residential neighbourhoods considered somewhat or very polluted due to high urban traffic density are in all likelihood more paved and therefore hotter. Such highly paved and thus more impermeable neighbourhoods contribute to the negative health impacts of heat [2, 3, 60, 61] and result in the deployment of various adaptations [9, 10].
In addition, hormonal differences may explain the difference in self-reported adverse health impacts by gender [62, 63] for women aged 45 to 64 in particular (menopause period) . Regardless of gender, the 45–64 age group, however, remain at higher odds of impacts according to our results. This is contrary to what is generally reported in this area of research  but plausible in very poor populations such as ours. In fact, the proportion of 45–64 years old having more chronic health problems increases according to a decrease in family income , while seniors having survived the same conditions would in all likelihood have been in better health. Another explanation for the increased odds of heat impacts for the 45–64 group would be greater exposure, mainly because of the requirement to leave home for family or work reasons  compared to older people, who are more confined . Clarifying this issue would be important for public health because those 65 and over are generally the target group for national heatwave plans.
Physical inactivity, a risk factor for several non-communicable diseases , is associated with a higher prevalence of self-reported adverse health impacts. The average level of physical activity declines with age; the result is a downward trend in fitness [64, 65]. It is known that poor physical fitness leads to a low cardiovascular reserve and a low tolerance to humid heat . However, physical inactivity was not associated with self-reported adverse health impacts in multivariate models in our previous articles [35, 36]. It is thus possible that this multilevel analysis highlights a potential contribution of deprived neighbourhoods to physical inactivity [67, 68]. Further studies are needed on this.
In our study, physical inactivity, like functional limitations, is also associated with a lower self-reported adaptation to heat, according to the index. More specifically, these two conditions characterize very poor populations [64, 65], such as our sample drawn from DAs that are farther (in time and distance) from cool greenspaces than more affluent DAs .
Improving the urban development of the poorest sectors of large urban centres would then be more beneficial for countering the health impacts of heat and facilitating adaptation. Over a third of respondents in our study also expressed this need. In addition, this initiative would have many co-benefits [55, 57]. For example, according to our results, DAs that offered better walkability were on average more adaptive according to the index, which could mean a reduction in the use of air conditioning in the home , which contributes to the outside heat and indirectly by exacerbating the air pollution in some cases . Similarly, better equipped DAs could reduce residential mobility and improve levels of well-being in neighbourhoods [71, 72] while encouraging adaptation by means other than air conditioning in the dwelling. These recent hypotheses, however, remain to be confirmed.
As already published [7, 42], but not in our previous multivariate analyses [35, 36], social support is associated with health impacts from heat with this multilevel approach. According to our results in the present article, the network that helps prevent impacts would be less diverse and more likely to be family than the network that promotes adaptation in a way other than air conditioning in the home. Further studies are needed in this respect, as networks with strong links could potentially exacerbate rather than reduce vulnerability to the effects of heat in the elderly, according to some authors .
The relationship between satisfaction with the indoor temperature of the dwelling in summer and the decrease in the prevalence of self-reported adverse health impacts, as well as the relationship between satisfaction with insulation quality and reduced adoption of behaviours measured by index of adaptation, clearly highlights the contribution of the building occupied in reducing exposure to heat, beyond the individual sensitivities. The fact that each of these relationships is observed for both the dwelling (individual-level) and the building (building-level) reinforces the importance of the use of these perceptions for public health surveillance. Especially in highly disadvantaged areas of major urban centres, where there are old buildings whose insulation meets less efficient standards than today’s, unless upgrading has been done since they were built .
Last, in this study the variables were self-reported which is generally considered easy to implement and cost-effective in several situations. As reported elsewhere , the validity of self-reported versus medical-based diagnoses and behaviours has been well established over time, several countries and data collection methods, especially as a tool for predicting future risks and as an epidemiologic survey tool for prevention and public health actions. Perceived health is thus a reliable and valid subjective measurement of the overall state of health and is widely used by statistical and health agencies throughout the world. This self-reporting approach for evaluating environmental exposures is also documented if not as well assessed. Moreover, the objective measures are far from perfect, as already mentioned [12, 17]. From a public health perspective, it remains important to validate these results in other contexts.
Limitations of the study
While the study’s response rate was low, the rate by question was very high. The latter is a another measure of the survey’s response rate . To our knowledge there have been no comparable studies published, particularly with respect to certain key aspects of the sample selection criteria (very disadvantaged DAs in Canada) and therefore evaluating the response rate in light of comparative work is not possible. However, the response rate reported for a qualitative study undertaken in two of the same sample sites as our study  and thus reflecting similar community characteristics such as large urban centres, multiethnic environments, high density, etc., indicates that our response rate could represent what can realistically be obtained in the very disadvantaged DAs of Québec’s largest cities. It is also in line with recent trends in response rates in North America .
In following our ethics protocol, no information was collected from people refusing participation. In response to this limitation, some statistical comparisons with census data available by DA were undertaken at the Institut national de santé publique du Québec. In fact, there are only two minor discrepancies between our sample and census data. Indeed, our sample is older (mean age 53 vs 47) and richer (mean income 30 k vs 24 k). This might be explained by the oversampling of public housing in the study (half of sample). That said, these differences should not influence health impacts and adaptation behaviours because all the proposed models have been adjusted for age, while income was not a significant variable (see Tables 1 and 2). Therefore we can state that despite the low classic response rate, the study samples acceptably represented the populations of the DAs included in the study as well as those of the very disadvantaged DAs in Québec’s nine largest cities, due to the sampling plan that was adopted in the study (for more details, see [35, 36]. It follows that reliable generalizations of this study’s results to deprived neighbourhoods with similar characteristics can be made.
This research, conducted in the most disadvantaged DAs in the most populous cities in Québec in 2011, is highly innovative in many ways. In particular, it shows that the use of a multi-level approach helps better identify the differential contribution of the characteristics of residents, buildings and their surroundings on self-reported adverse health impacts and adaptation (other than air conditioning in the home) when it is very hot and humid during the summer. On the other hand, both at the individual-level and building-level, satisfaction levels with indoor temperature in summer would be good indicators of exposure to heat at home, while the quality of satisfaction levels with their thermal insulation (in summer) are good indicators of the ability to counter exposure to heat. Considered more broadly, it appears that some indicators are more specific to impacts while others primarily characterize adaptation. Finally, people who feel negative health impacts from heat rely more on various mechanical and non-mechanical adaptation strategies as other residents of very disadvantaged DAs. Thus, adaptation to heat might be deployed in response to the impacts of heat, rather than to prevent them.
Further studies are needed to test these hypotheses. In particular, it would prove important to develop a theoretical framework ranging from exposure to self-reported adverse health impacts under very hot and humid summer conditions, including adaptation to heat as well as the characteristics and perceptions related to the life environment. Filling this gap will make a major contribution to research because the establishment of a theoretical framework helps clarify the type of links that exist between the indicators associated with the issues, define how to measure and improve monitoring, and improve surveillance and the public health response. It would also help the eventual development of national databases on housing and the neighbourhood. The establishment of observatories  in public health and climate change would support the achievement of these goals.
Intergovernmental panel on climate change. Summary for policymakers. In: Stocker TF, editor. Climate change 2014: impacts, adaptation and vulnerability. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. New York: Cambridge University Press; 2014. Retrieved from http://www.ipcc.ch/report/ar5/wg2/.
Mitchell R, Popham F. Effect of exposure to natural environment on health inequalities: an observational population study. Lancet. 2008;372(9650):1655–60.
Ngom R, Gosselin P, Blais C, Rochette L. Adaptation to climate change in environmental health through primary prevention: An applied example with green spaces for urbanized regions in the province of Quebec, paper presented at the XVth International Medical Geography Symposium. Michigan: East Lansing; 2013.
McGeehin MA, Mirabelli M. The potential impacts of climate variability and change on temperatuyre-related morbidity and mortality in the United States. Environ Health Perspect. 2001;109(2):185–9.
van den Akker M, Buntinx F, Metsemakers JFM, Roos S, Knottnerus JA. Multimorbidity in general practice: prevalence, incidence, and determinants of Co-occurring chronic and recurrent diseases. J Clin Epidemiol. 1998;51(5):367–75.
Kilbourne EM. Heat-related illness: current status of prevention efforts. Am J Prev Med. 2002;22(4):328–9.
Bouchama A, Dehbi M, Mohamed G, Matthies F, Shoukri M, Menne B. Prognostic factors in heat wave–related deaths: a meta-analysis. Arch Intern Med. 2007;167(20):2170–6.
Vandentorren S, Bretin P, Zeghnoun A, Mandereau-Bruno L, Croisier A, Cochet C, et al. August 2003 heat wave in France: risk factors for death of elderly people living at home. Eur J Public Health. 2006;16(6):583–91.
Alberini A, Gans W, Alhassan M. Individual and public-program adaptation: coping with heat waves in five cities in Canada. Int J Environ Res Public Health. 2011;8(12):4679–701.
White-Newsome JL, Sánchez BN, Parker EA, Dvonch JT, Zhang Z, O’Neill MS. Assessing heat-adaptive behaviors among older, urban-dwelling adults. Maturitas. 2011;70(1):85–91.
Smit B, Wandel J. Adaptation, adaptive capacity and vulnerability. Glob Environ Chang. 2006;16(3):282–92.
Ormandy D, Ezratty V. Health and thermal comfort: from WHO guidance to housing strategies. Energy Policy. 2012;49:116–21.
Anderson M, Carmichael C, Murray V, Dengel A, Swainson M. Defining indoor heat thresholds for health in the UK. Perspect Public Health. 2013;133(3):158–64.
Howden-Chapman P, Matheson A, Crane J, Viggers H, Cunningham M, Blakely T, et al. Effect of insulating existing houses on health inequality: cluster randomised study in the community. BMJ. 2007;334(7591):460.
Berry K, Clarke K-L, Pajot M, Hutton D. Risk perception, health communication, and adaptation to the health impacts of climate change in Canada. In: J Ford, Berrang Ford, Lea, editors. Climate change adaptation in developed nations. Springer; 2011. p. 205-19.
Statistics Canada. Perceived health. 2010. [cited 2013 6 December 2013]; Retrieved from: http://www.statcan.gc.ca/pub/82-229-x/2009001/status/phx-eng.htm.
Ezratty V, Duburcq A, Emery C, Lambrozo J. Liens entre l’efficacité énergétique du logement et la santé des résidents: résultats de l’étude européenne LARES. Environ Risques Santé. 2009;8(6):497–506.
Jamrozik E, Hyde Z, Alfonso H, Flicker L, Almeida O, Yeap B, et al. Validity of self-reported versus hospital-coded diagnosis of stroke: a cross-sectional and longitudinal study. Cerebrovasc Dis. 2014;37(4):256–62.
Fahimi M, Link M, Mokdad A, Schwartz DA, Levy P. Tracking chronic disease and risk behavior prevalence as survey participation declines: statistics from the behavioral risk factor surveillance system and other national surveys. Prev Chronic Dis. 2008;5(3):A80.
Starr GJ, Dal Grande E, Taylor AW, Wilson DH. Reliability of self-reported behavioural health risk factors in a South Australian telephone survey. Aust N Z J Public Health. 1999;23(5):528–30.
Pierannunzi C, Hu SS, Balluz L. A systematic review of publications assessing reliability and validity of the Behavioral Risk Factor Surveillance System (BRFSS), 2004-2011. BMC Med Res Methodol. 2013;13(1):49.
Daniau C, Dor F, Eilstein D, Lefranc A, Empereur-Bissonnet P, Dab W. Étude de la santé déclarée par les personnes riveraines de sources locales de pollution environmentale: une revue. Seconde partie: analyse des résultats et perspectives. Rev Epidemiol Sante Publique. 2013;61(4):388–98.
Anderson BG, Bell ML. Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009;20(2):205–13.
Hansen J, Sato M, Ruedy R. Perception of climate change. Proc Natl Acad Sci. 2012;109(37):E2415–23.
Uejio CK, Wilhelmi OV, Golden JS, Mills DM, Gulino SP, Samenow JP. Intra-urban societal vulnerability to extreme heat: the role of heat exposure and the built environment, socioeconomics, and neighborhood stability. Health Place. 2011;17(2):498–507.
Oke TR. Boundary layer climates. 2nd ed. London: Routledge; 1987. p. 474.
Voogt JA. Urban heat island. In: Douglas I, editor. Encyclopedia of global environmental change - Causes and consequences of global environmental change. New York: Wiley; 2002. p. 660–6.
Canadian Council on Social Development. Poverty by geography urban: urban poverty in Canada, 2000. 2007. Retrieved from http://www.ccsd.ca/images/research/UPP/PDF/UPP-PovertyByGeography.pdf.
Canada S. Dissemination Area (DA). 2012; Retrieved from: https://www12.statcan.gc.ca/census-recensement/2011/ref/dict/geo021-eng.cfm.
Gouvernement du Québec. L’organisation municipale et régionale au Québec en 2013. 2013. Retrieved from http://www.mamrot.gouv.qc.ca/pub/organisation_municipale/organisation_territoriale/organisation_municipale.pdf.
Diez-Roux AV. Multilevel analysis in public health research. Annu Rev Public Health. 2000;21(1):171–92.
Riva M, Gauvin L, Barnett TA. Toward the next generation of research into small area effects on health: a synthesis of multilevel investigations published since July 1998. J Epidemiol Community Health. 2007;61(10):853–61.
Vallée J, Souris M, Fournet F, Bochaton A, Mobillion V, Peyronnie K, et al. Sampling in health geography: reconciling geographical objectives and probabilistic methods. An example of a health survey in Vientiane (Lao PDR). Emerg Themes Epidemiol. 2007;4(1):6.
Pampalon R, Raymond G. A deprivation index for health and welfare planning in Quebec. Chronic Dis Can. 2000;21(3):104–13.
Bélanger D, Gosselin P, Valois P, Abdous B. Perceived adverse health effects of heat and their determinants in deprived neighbourhoods: a cross-sectional survey of nine cities in Canada. Int J Environ Res Public Health. 2014;11(11):11028–53.
Bélanger D, Gosselin P, Valois P, Abdous B. Neighbourhood and dwelling characteristics associated with the self-reported adverse health effects of heat in most deprived urban areas: a cross-sectional study in 9 cities. Health Place. 2015;32:8–18.
Statistics Canada. Canadian community health survey : 2007 Questionnaire. 2007. Retrieved from http://www23.statcan.gc.ca/imdb-bmdi/pub/instrument/3226_Q1_V4-eng.pdf.
Greenacre M. Correspondence analysis in practice. 2nd ed. New York: Chapman & Hall/CRC; 2007.
Asselin L-M. Composite indicator of multidimensional poverty. 2002. Retrieved from http://www.pep-net.org/sites/pep-net.org/files/typo3doc/pdf/CBMS_training/composite_ind.pdf.
Bélanger D, Abdous B, Gosselin P, Valois P. An adaptation index to high summer heat associated with adverse health impacts in deprived neighborhoods. Clim Change. 2015;132(2):279–93.
Kovats RS, Hajat S. Heat stress and public health: a critical review. Annu Rev Public Health. 2008;29:41–55.
Lundgren L, Jonsson A. Assessment of social vulnerability: a literature review of vulnerability related to climate change and natural hazards. 2012. Retrieved from http://www.cspr.se/briefings/1.358869/CSPRBriefing9.pdf.
Chebana F, Martel B, Gosselin P, Giroux J-X, Ouarda TMJ. A general and flexible methodology to define thresholds for heat health watch and warning systems, applied to the province of Québec (Canada). Int J Biometeorol. 2013;57(4):631–44.
Boulfroy E, Khaldoune J, Grenon F, Fournier R, Talbot B. Conservation des îlots de fraîcheur urbains. Description de la méthode suivie pour identifier et localiser les îlots de fraîcheur et de chaleur - Résumé (Rapport 2012-11a). 2012. Retrieved from http://www.cerfo.qc.ca/index.php?id=16&no_cache=1&tx_drblob_pi1[downloadUid]=313.
INSPQ. Indicateur de l’environnement bâti. Retrieved from: https://www.inspq.qc.ca/environnement-bati/indicateurs/indice-de-potentiel-pietonnier.
Rabe-Hesketh S, Skrondal A. Multilevel and longitudinal modeling using stata. 2nd ed. College Station: Stata Press; 2008.
Rasbash J, Charlton C, Browne WJ, Healy M, Cameron B. MLwiN Version 2.1. 2009: Centre for multilevel modelling, University of Bristol.
Steele F. Module 6: regression models for binary responses concepts. [Course custom textbook]. 2009. Retrieved from http://www.bristol.ac.uk/cmm/learning/module-samples/6-concepts-sample.pdf.
Tofighi D. Centering predictor variables in cross-sectional multilevel models: a new look at an old issue. Psychol Methods. 2007;12:121–38.
Diez-Rioux AV. The study of group-level factors in epidemiology: rethinking variables, study designs, and analytical approaches. Epidemiol Rev. 2004;26:104–11.
Chaix B, Chauvin P. L’apport des modèles multiniveau dans l’analyse contextuele en épidémiologie sociale: une revue de la littérature. Rev Epidemiol Sante Publique. 2002;50(5):489–99.
Pickett KE, Pearl M. Multilevel analyses of neighbourhood socioeconomic context and health outcomes: A critical review. J Epidemiol Commun Health (1979-). 2001;55(2):111–22.
Diez-Roux AV. Next steps in understanding the multilevel determinants of health. J Epidemiol Community Health. 2008;62(11):957–9.
Holmes JH, Lehman A, Hade E, Ferketich AK, Gehlert S, Rauscher GH, et al. Challenges for multilevel health disparities research in a transdisciplinary environment. Am J Prev Med. 2008;35(2, Supplement):S182–92.
Rosso AL, Auchincloss AH, Michael YL. The urban built environment and mobility in older adults: a comprehensive review. J Aging Res. 2011;2011:816106.
Wilhelmi O, Wilhelmi M, Hayden MH. Connecting people and place: a new framework for reducing urban vulnerability to extreme heat. Environ Res Lett. 2010;5(1):014021.
Patz JA, Frumkin H, Holloway T, Vimont DJ, Haines A. Climate change: challenges and opportunities for global health. JAMA. 2014;312(15):1565–80.
Luber G, McGeehin M. Climate change and extreme heat events. Am J Prev Med. 2008;35(5):429–35.
Kenny GP, Yardley J, Brown C, Sigal RJ, Jay O. Heat stress in older individuals and patients with common chronic diseases. Can Med Assoc J. 2010;182(10):1053–60.
Smargiassi A, Berrada K, Fortier I, Kosatsky T. Traffic intensity, dwelling value, and hospital admissions for respiratory disease among the elderly in Montreal (Canada): a case-control analysis. J Epidemiol Community Health. 2006;60(6):507–12.
Stone B, Vargo J, Liu P, Habeeb D, DeLucia A, Trail M, et al. Avoided heat-related mortality through climate adaptation strategies in three US cities. PLoS One. 2014;9(6):1–8.
Gagnon D, Kenny GP. Does sex have an independent effect on thermoeffector responses during exercise in the heat? J Physiol. 2012;590(Pt 23):5963–73.
Charkoudian N. Skin blood flow in adult human thermoregulation: how it works, when it does not, and why. Mayo Clin Proc. 2003;78(5):603–12.
National Center for Health Statistics. Health, United States, 2011: with special feature on socioeconomic status and health. 2011. Retrieved from http://www.cdc.gov/nchs/data/hus/hus11.pdf.
Public Health Agency of Canada. Chapter 3: the health and well-being of Canadian seniors. In: The chief public health officer’s report on the state of public health in Canada 2010. 2010. Retrieved from http://www.phac-aspc.gc.ca/cphorsphc-respcacsp/2010/fr-rc/cphorsphc-respcacsp-06-eng.php.
Koppe C, Kovats S, Jendritzky G, Menne B. Heat-waves: risks and responses. 2004. Retrieved from http://www.euro.who.int/__data/assets/pdf_file/0008/96965/E82629.pdf.
Gordon-Larsen P, Nelson MC, Page P, Popkin BM. Inequality in the built environment underlies key health disparities in physical activity and obesity. Pediatrics. 2006;117(2):417–24.
Leal C, Chaix B. The influence of geographic life enviroments on cardiometabolic risk factors: a systematic review, a methodological assessment and a research agenda. Obes Rev. 2011;12(3):217–30.
Ngom R, Gosselin P, Blais C. Reduction of disparities in access to green spaces: their geographic insertion and recreational functions matter. Appl Geogr. 2015. In Press.
Ca VT, Asaeda T, Abu EM. Reductions in air conditioning energy caused by a nearby park. Energy Build. 1998;29(1):83–92.
Oishi S. The psychology of residential mobility: implications for the self, social relationships, and well-being. Perspect Psychol Sci. 2010;5(1):5–21.
Gifford R. Environmental psychology matters. Annu Rev Psychol. 2014;65(1):541–79.
Wolf J, Adger WN, Lorenzoni I, Abrahamson V, Raine R. Social capital, individual responses to heat waves and climate change adaptation: an empirical study of two UK cities. Glob Environ Chang. 2010;20(1):44–52.
Statistics Canada. Canadian community health survey (CCHS) 2011 sampling design and collection. 2011. 16 December 2013]; Retrieved from: http://www12.statcan.gc.ca/nhs-enm/2011/ref/nhs-enm_guide/guide_2-eng.cfm.
Morin P, Leloup X, Baillergeau E, Caillouette J. Habiter en HLM : impacts sur la santé et le bien-être des ménages familiaux. 2010. Retrieved from http://flhlmq.com/files/images-dans-les-textes/Habiter%20en%20HLM%20-%20Impacts%20%20Le%20r%C3%A9sum%C3%A9.pdf.
Choudhury Y, Hussain I, Parsons S, Rahman A, Eldridge S, Underwood M. Methodological challenges and approaches to improving response rates in population surveys in areas of extreme deprivation. Prim Health Care Res Dev. 2012;13(03):211–8.
Hemmings J, Wilkinson J. What Is a public health observatory? J Epidemiol Commun Health (1979-). 2003;57(5):324–6.
We wish to thank the Green Fund in the framework of Action 21 of the Québec government’s 2006–2012 Climate Change Action Plan for funding this study, as well as the Institut national de santé publique du Québec, the Regroupement des offices municipaux du Québec, the municipal housing authorities of the nine largest cities of Québec in 2011, the Société d’habitation du Québec, the Léger Marketing polling and market research firm, and the 3485 participants.
The authors declare that they have no competing interests.
DB, BA, PV and PG contributed to the design of the study; data collection was supervised by DB and PG. Analysis was performed by DB, BA and EALS, with all authors contributing to the writing of the paper. All authors read and approved the final manuscript.
Independent variables by level. (DOCX 19 kb)
3-level multivariate logistic regression model of the self-reported adverse health impacts when it is very hot and humid in the summer: random part. (DOCX 16 kb)
Individual-level covariables associated with the prevalence of self-reported adverse health impacts when it is very hot and humid in summer: bivariate analysis. (DOCX 26 kb)
3-level multivariate linear regression model of the adaptation index when it is very hot and humid in summer: random part. (DOCX 19 kb)
Individual-level covariables associated with the adaptation index when it is very hot and humid in summer. (DOCX 20 kb)
Synthesis of the results regarding the fixed effects of the Impact and Index models. (DOCX 18 kb)
About this article
Cite this article
Bélanger, D., Abdous, B., Valois, P. et al. A multilevel analysis to explain self-reported adverse health effects and adaptation to urban heat: a cross-sectional survey in the deprived areas of 9 Canadian cities. BMC Public Health 16, 144 (2016). https://doi.org/10.1186/s12889-016-2749-y